Papers for Tuesday, Feb 25 2025

Rotating black holes are the most powerful source of energy in the known universe, and are the cause of some of the most spectacular and extreme astronomical phenomena. The goal of this article is to analyze in simple terms the physics of energy extraction in rotating black holes. Specifically, the source of said energy, the efficiency of the energy extraction process, and some specific mechanisms that allow said extraction are analyzed. The article is intended primarily for undergraduate students of physics, astronomy and related fields.

Galaxy formation is a complex problem that connects large scale cosmology with small scale astrophysics over cosmic timescales. Hydrodynamical simulations are the most principled approach to model galaxy formation, but have large computational costs. Recently, emulation techniques based on Convolutional Neural Networks (CNNs) have been proposed to predict baryonic properties directly from dark matter simulations. The advantage of these emulators is their ability to capture relevant correlations, but at a fraction of the computational cost compared to simulations. However, training basic CNNs over large redshift ranges is challenging, due to the increasing non-linear interplay between dark matter and baryons paired with the memory inefficiency of CNNs. This work introduces EMBER-2, an improved version of the EMBER (EMulating Baryonic EnRichment) framework, to simultaneously emulate multiple baryon channels including gas density, velocity, temperature and HI density over a large redshift range, from z=6 to z=0. EMBER-2 incorporates a context-based styling network paired with Modulated Convolutions for fast, accurate and memory efficient emulation capable of interpolating the entire redshift range with a single CNN. Although EMBER-2 uses fewer than 1/6 the number of trainable parameters than the previous version, the model improves in every tested summary metric including gas mass conservation and cross-correlation coefficients. The EMBER-2 framework builds the foundation to produce mock catalogues of field level data and derived summary statistics that can directly be incorporated in future analysis pipelines. We release the source code at the official website this https URL.

Ultrahigh-energy cosmic rays (UHECRs) experience deflections as they traverse the Galactic magnetic field (GMF), which must be accounted for when tracing them back to their sources. After briefly summarizing our results on uncertainties in cosmic-ray deflections from the UF23 ensemble of GMF models (Unger & Farrar, 2024), we report a new preliminary fit of the GMF including foreground emission from the Local Bubble. This fit uses the analytic model of Pelgrims et al. (2024) for the magnetic field in the thick shell Galactic bubbles. We also discuss how variations in toroidal halo field modeling account for the key differences between the Jansson & Farrar (2012) GMF model and the UF23 ensemble. We also extend our previous analysis of the origin of the highest-energy "Amaterasu" event observed by the Telescope Array to include the four highest-energy events detected by the Pierre Auger Observatory. Amaterasu and PAO070114 are the UHECR events with the smallest localization uncertainties of 4.7% and 2.4%, respectively. Neither of their back-tracked directions aligns with any compelling candidate for a continuous UHECR accelerator. This strengthens the evidence that at least a fraction of the highest energy events originate from transient sources.

Observed dwarf galaxies tend to have linearly rising rotation curves, which indicate flat density cores in their centers. Furthermore, disk galaxies show a wide range of rotation curves shapes. High resolution simulations of cold collisionless dark matter do not reproduce flat central profiles, or the observed diversity of rotation curve shapes; even hydrodynamic simulations incorporating baryonic feedback cannot do that robustly. However, numerical simulations are not the only way to make predictions about density profiles of equilibrium dark matter halos. A theoretical model based on statistical mechanics shows that maximum entropy solutions for cold collisionless self-gravitating dark matter halos can have a range of inner density profiles, including flat density cores. These theoretical profiles, called DARKexp, have only one shape parameter, and are able to fit the observed rotation curves of galaxies with last measured velocities in the range ~20-200 km/s. Here we present fits to 96 SPARC catalog galaxies, and the Milky Way. DARKexp also provides good fits to the projected stellar density distributions of ultrafaint dwarfs that show cores, suggesting that the dark matter halo hosts could have flat density cores. Thus, DARKexp appears to be able to address the core-cusp problem and the diversity of rotation curves with cold collisionless dark matter alone, without baryonic feedback.

Dark matter (DM) halos simulated via N-body techniques are known to exhibit nearly universal profiles at equilibrium; however, their origin remains uncertain despite thorough investigation. This work aims to probe the origin of simulated DM halo structure by testing DARKexp, a first-principles approach that describes the maximum entropy state of collisionless, self-gravitating systems in equilibrium via their energy distributions, and proposes an entropy functional that is expected to increase during evolution. We fit the DARKexp energy distribution to a set of massive equilibrium halos from the cosmological simulation IllustrisTNG. For the first time, we calculate the entropy of these halos as a function of cosmological time by tracking halos that are in equilibrium at z = 0.0, to z = 3.0 and calculating their entropy at various epochs between. We find that DARKexp provides an excellent fit to the energy distributions of equilibrium DM halos and that such halos exhibit an overall increase in entropy during evolution. Our results indicate that DM halos evolve to become their maximum entropy state at equilibrium and that this state is described by DARKexp.

Recently, it was proposed that an off-center dipole magnetic configuration, together with a non-trivial temperature profile, may be the best model to explain the X-ray light curve of PSR J0030+0451 observed by the Neutron Star Interior Composition Explorer (\emph{NICER}). Using a theoretical model for the electric current density in a force-free pulsar magnetosphere, we compute from first principles the distribution of electric current over the polar cap associated with an off-center magnetic dipole. We then use a simple prescription to compute the resulting temperature distribution, which allows us to derive the observed X-ray light curve. We investigate the role of the volumetric return current region in the polar cap and find that although it does not make a big difference in an aligned dipole case, the difference can be bigger in the case of an off-center dipole. Finally, we apply Markov Chain Monte Carlo (MCMC) fitting to the X-ray light curves of pulsars PSR J0030+0451 and PSR J0437--4715 with and without the volumetric return current, and find that our model can reasonably recover the observed X-ray light curves.

The JCMT Transient Survey recently discovered that the Class 0 protostar HOPS 358 decreased in 350 GHz continuum brightness by $\sim25$% over the course of four years before brightening again for the next four. The JCMT lightcurve can be fit by a long timescale dip lasting roughly eight years. A shorter timescale periodicity is also apparent with a period of 1.75 years and a small 3% amplitude. NEOWise monitoring reveals that the mid-IR wavelength brightness of HOPS 358 follows a similar long-term pattern in time. Here, we present a study of nine epochs of ALMA observations of HOPS 358 taken over the course of the decline and subsequent rise in brightness seen with the JCMT to test whether the variation seen on $\sim15"$ scales, covering both disk and envelope, is also observed on smaller, $<1"$ scales that primarily probe HOPS 358's protostellar disk. We detect both HOPS 358 and its southern companion, HOPS 358B, in our ALMA observations, and find that at least one of the two is varying. Assuming that HOPS 358 is the variable, the light curve has the same shape as that found by the JCMT. Additionally, our high resolution ALMA imaging of HOPS 358 reveals that the disk is warped, with a $16^{\circ}$ warp at a disk radius of 35 au, about halfway through the extent of the disk. The physical origin of the warp along with how it relates to the variability seen towards HOPS 358, however, remain unclear.

The origin of heavy r-process elements in the universe is still a matter of great debate, with a confirmed scenario being neutron star (NS) mergers. Additional relevant sites could be specific classes of events, such as gamma-ray burst (GRB) Supernovae (SNe), where a central engine could push neutron-rich material outwards, contributing to the ejecta of the massive exploding star. Here, we investigate our ability to infer the production of heavy elements in such scenarios, on the basis of the observed nebular emission. We solve the steady-state ionization, level population, and thermal balance, for optically thin ejecta in non-local thermodynamic equilibrium (NLTE), in order to explore the role of heavy elements in cooling the gas, and their imprint in the emergent spectrum a few hundreds days post-explosion. We find that heavy elements would be relevant in the cooling process of the nebula only if they account for at least $\sim1\%$ of the total ejected mass, at the typical kinetic temperatures of a few thousands K. However, even in the absence of such amount, a few $0.1\%$ of the total ejected mass could be instead sufficient to leave a detectable imprint around $\sim1-10~\mathrm{\mu m}$. This wavelength range, which would be relatively clean from features due to light elements, would be instead robustly populated by lines from heavy elements arising from forbidden transitions in their atomic fine structures. Hence, the new generation of telescopes, represented by the James Webb Space Telescope (JWST), will most likely allow for their detection.

The fourth flight of the Focusing Optics X-ray Solar Imager sounding rocket experiment (FOXSI-4) aimed to achieve the first imaging spectroscopic observations of mid-to-large class ( >= GOES C5 class) solar flares, in contrast to the previous three flights that targeted relatively quiet regions of the Sun. To meet the emerging requirements for hard X-ray focal plane detectors for providing simultaneous diagnostics of spectrally (<1 keV FWHM) and spatially (<50 um) separated coronal and chromospheric emissions from solar flares, we developed a new strip-configuration detector called the wide-gap CdTe semiconductor double-sided strip detector (CdTe-DSD). The wide-gap CdTe-DSD employs a unique design principle to enhance position resolution. This enhancement is realized by expanding the gaps between electrodes to induce charge-sharing across adjacent strip electrodes and using this sharing energy information for position reconstruction to a level finer than the strip-pitch. However, the detector response becomes complex and requires consideration of various factors, such as the charge loss due to wider gaps and the dependence on the depth of photon interaction. Thus, we developed an energy reconstruction method that fully leverages the energy information between adjacent strips and from both the cathode and anode sides, achieving an energy resolution of 0.75 keV (FWHM) at 14 keV. Furthermore, we conducted an X-ray scanning experiment using a synchrotron beam at Spring-8 to evaluate the detector response with a fine scale of 10 um. Based on these results, we established a sub-strip position reconstruction method, demonstrating that X-rays interacting at the center of the gap can be determined with an accuracy of 20 um, and even those at the strip center can be determined with an accuracy of 50 um.

We present equivalent widths, improved model atmosphere parameters, and revised abundances for 14 species of 11 elements derived from high resolution optical spectroscopy of 311 metal-poor stars. All of these stars had their parameters previously published by Roederer et al. We use color-Teff relationships calibrated for Gaia and 2MASS photometry to calculate improved effective temperatures. We calculate log of surface gravity values using measurements derived from Gaia parallaxes and other fundamental stellar properties. We perform a standard LTE abundance analysis using MARCS model atmospheres and the MOOG line analysis software to rederive microturbulence velocity parameters, metallicities, and abundances based on O I, Na I, Mg I, Si I, K I, Ca I, Ti I, Ti II, Cr I, Cr II, Fe I, Fe II, Ni I, and Zn I lines using the previously measured equivalent widths. On average, the new Teff values are 310 K warmer, the new log g values are higher by 0.64 dex, and the new [Fe/H] values are higher by 0.26 dex. We apply NLTE corrections to the abundances derived from O I, Na I, Mg I, Si I, K I, Fe I, and Fe II lines. Our sample contains 6 stars with [Fe/H] < -3.5, 28 stars with [Fe/H] < -3.0, and 113 stars with [Fe/H] < -2.5. Our revised abundances for these 311 stars are now in better agreement with those derived by previous studies of smaller samples of metal-poor stars in the Milky Way.

The TRGB-SBF Project team is developing an independent distance ladder using a geometrical calibration of the tip of the red giant branch (TRGB) method in elliptical galaxies that can in turn be used to set the surface brightness fluctuation (SBF) distance scale independent of Cepheid variables and Type Ia supernovae. In this paper, we use JWST TRGB distances calibrated using the megamaser galaxy NGC 4258 to determine a new Cepheid-independent SBF zero point for HST. This new calibration, along with improved optical color measurements from PanSTARRS and DECam, gives an updated value of the Hubble-Lemaitre constant $H_0 = 73.8 \pm 0.7 \pm 2.3$ km/s/Mpc that is virtually identical to the previously published result based on SBF distances calibrated using Cepheids (Blakeslee et al. 2021, arXiv:2101.02221). Future improvements in the geometrical calibration of TRGB and additions to the HST SBF survey will further reduce the random and systematic uncertainties on $H_0$ measured using the extensive HST WFC3/IR SBF dataset. JWST observations of the Coma cluster will soon establish a foundation for the SBF calibration that will extend to much larger distances than is possible with HST.

HeH$^+$ belongs to the class of "reactive" ions that can be destroyed so quickly that chemical formation and destruction rates may compete with inelastic rates and should be considered when solving the statistical equilibrium equations. This so-called chemical "pumping" or "excitation" effect is investigated here for the first time in HeH$^+$. The chemical evolution of HeH$^+$ in NGC 7027 is modeled with the CLOUDY photoionization code using updated reaction rate coefficients. The non-LTE analysis of the three observed HeH$^+$ emission lines is then performed with the CLOUDY and RADEX codes using an extensive set of spectroscopic and inelastic collisional data suitable for the specific high-temperature environment of NGC 7027. In a second approach, chemical formation and destruction rates of HeH$^+$ are implemented in RADEX. This code is combined with MCMC sampling (performed on the RADEX-parameters space) to extract the best-fit HeH$^+$ column density and physical conditions from the observed line fluxes. The CLOUDY and RADEX non-LTE results are found to be in good agreement, and the $\upsilon=1-0 \ P(2)/P(1)$ line ratio is better than 20%. Agreement to better than a factor of 2.3 is obtained when including the reaction between He($2^3S$) and H as an additional source of HeH$^+$. The RADEX/MCMC model with chemical pumping is found to reproduce both the observed line fluxes and the line ratio to 20%. Our results suggest that additional HeH$^+$ lines must be detected in NGC 7027 to better constrain the physical conditions via non-LTE models. Uncertainties in collisional (reactive and inelastic) data of HeH$^+$ have been largely reduced in this work. The three observed lines are not sensitive to chemical pumping while excited "short-lived" levels are significantly overpopulated with respect to a non-LTE model neglecting chemical excitation.

Precise tracking and measurement of the energy carried by the individual magnetohydrodynamic (MHD) modes has important implications and utility in astrophysical and laboratory plasmas. Previously, this was only achievable in limited linear MHD cases in the $\beta \ll 1$ or $\beta \gg 1$ regimes. In a series of papers, of which this is the third, we introduced the Eigenenergy Decomposition Method (EEDM) and derived exact analytical expressions for the modal energy components--called eigenenergies--of nonlinear 3D disturbances governed by the homogeneous ideal-MHD equations. Here, we extend the method to inhomogeneous ideal-MHD by introducing a source term accounting for gravity, and provide detailed guidelines for applying the decomposition scheme to any general inhomogeneous quasi-linear PDEs that possess a globally conserved quantity, beyond the realm of MHD. Furthermore, we show that the eigenenergies can be used to locate and measure nonlinear mode conversions, which is an additional feature of the method. Finally, we provide well-categorized context for the application of the method to simulations and discuss the possible numerical inaccuracies that may inevitably arise due to discretization. This paper provides a more mature description of the method and its interpretation, and is recommended as the starting point for readers unfamiliar with the method.

Pulsar timing arrays recently found evidence for a gravitational wave background (GWB), likely the stochastic overlap of GWs from many supermassive black hole binaries. Anticipating a continuous gravitational wave (CW) detection from a single binary soon to follow, we examine how well current Bayesian methods can detect CWs and characterize their binary properties by modeling the response of the NANOGrav 15-year pulsar timing array to simulated binary populations. We run Markov Chain Monte Carlo searches for CWs in these datasets and compare them to quicker detection statistics including the optimal signal-to-noise ratio, matched filter detection statistic, and reduced log-likelihood ratio between the signal and noise models calculated at the injected parameters. The latter is the best proxy for Bayesian detection fractions, corresponding to a 50% detection fraction (by Bayes factors >10 favoring a CW detection over noise-only model) at a signal-to-noise ratio of 4.6. Source confusion between the GWB and a CW, or between multiple CWs, can cause false detections and unexpected dismissals. 53% of realistic binary populations consistent with the recently observed GWB have successful CW detections. 82% of these CWs are in the 4th or 5th frequency bin of the 16.03 yr dataset (6.9 nHz and 10.8 nHz), with 95 percentile regions spanning 4nHz-12nHz frequencies, $7-20\times10^9 M_\odot$ chirp masses, 60Mpc-8Gpc luminosity distances, and 18-13,000 sq. deg 68% confidence localization areas. These successful detections often poorly recover the chirp mass, with only 29% identifying the chirp mass accurately to within 1 dex with a 68% posterior width also narrower than 1 dex.

The distribution and physical conditions of molecular gas are closely linked to star formation and the subsequent evolution of galaxies. Emission from carbon monoxide (CO) and its isotopologues traces the bulk of molecular gas and provides constraints on the physical conditions through their line ratios. However, comprehensive understanding on how the particular choice of line modeling approach impacts derived molecular properties remain incomplete. Here, we study the nearby starburst galaxy M82, known for its intense star formation and molecular emission, using the large set of available multi-CO line observations. We present high-resolution (${\sim}85$ pc) emission of seven CO isotopologue lines, including $^{12}$CO, $^{13}$CO, and C$^{18}$O from the $J = 1-0$, $2-1$ and $3-2$ transitions. Using \texttt{RADEX} for radiative transfer modeling, we analyze M82\textsc{\char39}s molecular properties with (i) a one-zone model and (ii) a variable density model, comparing observed and simulated emissions via a minimum $\chi^2$ analysis. We find that inferred gas conditions -- kinetic temperature and density -- are consistent across models, with minimal statistical differences. However, due to their low critical densities (${<}10^{4}$ cm$^{-3}$), low-$J$ CO isotopologue lines do not effectively probe higher density gas prevalent in starburst environments like that of M82. Our results further imply that this limitation extends to high-redshift ($z{\gtrapprox}1$) galaxies with similar conditions, where low-$J$ CO lines are inadequate for density constraints. Future studies of extreme star-forming regions like M82 will require higher-$J$ CO lines or alternative molecular tracers with higher critical densities.

Stellar velocity dispersion ($\sigma$) of massive elliptical galaxies is a key ingredient to breaking the mass-sheet degeneracy and obtaining precise and accurate cosmography from gravitational time delays. The relative uncertainty on the Hubble constant H$_0$ is double the relative error on $\sigma$. Therefore, time-delay cosmography imposes much more demanding requirements on the precision and accuracy of $\sigma$ than galaxy studies. While precision can be achieved with an adequate signal-to-noise ratio (SNR), accuracy depends critically on factors such as elemental abundances and temperature of stellar templates, flux calibration, and wavelength ranges. We carry out a detailed study of the problem using multiple sets of galaxy spectra of massive elliptical galaxies with SNR$\sim$30-160 Å$^{-1}$, and state-of-the-art empirical and semi-empirical stellar libraries and stellar population synthesis templates. We show that the choice of stellar library is generally the dominant source of residual systematic error. We propose a general recipe to mitigate and account for residual uncertainties. We show that sub-percent accuracy can be achieved on individual spectra with our data quality. Covariance between velocity dispersions measured for a sample of spectra can also be reduced to sub-percent levels. We recommend this recipe for all applications that require high precision and accurate stellar kinematics, and make all software publicly available to facilitate its implementation. This recipe will be used in future TDCOSMO collaboration papers.

The Event Horizon Telescope (EHT) observations of Sgr A* resolved the photon ring at the galactic centre, revealing a diameter of $51.8~\mu \text{as}$. The ring-like structure is consistent with that of a Kerr black hole. However, the source of the high bright regions in the image and the time variability remain an open question. Besides the plasma properties and emission models, the spacetime geometry also holds an important role. We present an image depicting the bright hot spots consistent with Sgr A* observations at a wavelength $ \lambda $ = 1.3mm. The image is the result of an eclipsing Schwarzschild black hole situated along the line of sight between the galactic centre and Earth. The separation from both, primary (Sgr A*) and secondary (eclipsing) black holes is $10233\, {\rm AU}$. The central supermassive black hole located at the centre of the galaxy has an observational mass of $4.14\times10^6\,M_{\odot}$ and the secondary eclipsing black hole has an inferred mass of $1035\,M_{\odot}$.

We report the discovery and characterization of ZTF J0112+5827, a new magnetic cataclysmic variable with an orbital period of 80.9 minutes. ROSAT observations revealed X-ray emission with an average flux of $(68.4 \pm 15.7) \times 10^{-14}$ erg s$^{-1}$ cm$^{-2}$ (0.1--2.4 keV). The ZTF light curves show ellipsoidal-like variability in the $g$ band and two prominent humps at phases $\sim$0.0 and $\sim$0.7 in $i$ and $r$ bands. Spectroscopic observations with the Palomar 200-inch telescope revealed cyclotron emission features and strong He II and Balmer emission lines. Doppler tomography shows clear accretion streams with line-of-sight velocities of $\sim$500 km s$^{-1}$, but no accretion disk. Analysis of cyclotron harmonics indicates a magnetic field strength of $38.7^{+1.3}_{-1.1}$ MG, confirming ZTF J0112+5827 as a polar system containing a strongly magnetic white dwarf.

We performed variability analysis of the multiwavelength light curves for the flat-spectrum radio quasar PKS 0727-11. Using the generalized Lomb-Scargle periodogram, we identified a possible quasi-periodic oscillation (QPO) of $\sim$ 168.6 days (persisted for 6 cycles, with a significance of $3.8\sigma$) in the gamma-ray light curve during the flare period (MJD 54687-55738). It is the first time that periodic variations have been detected in this source, and further supported by other methods: weighted wavelet $z$-transform, phase dispersion minimization, REDFIT, autoregressive integrated moving average model, and structure function analysis. Cross-correlation analysis shows that there is a strong correlation between multi-band light variations, indicating that gamma-ray and radio flares may originate from the same disturbance, and the distance between the emission regions of gamma-ray and radio flares is calculated based on the time lag. We demonstrate that QPO arising from the non-ballistic helical jet motion driven by the orbital motion in a supermassive binary black hole is a plausible physical explanation. In this scenario, the estimated mass of the primary black hole is $M\sim3.66\times10^8-5.79\times10^{9}M_\odot$.

In this manuscript, we investigate observational correlations between the properties of gamma-ray bursts (GRBs) across the gamma, X-ray, and optical bands during the prompt and plateau phases of their light curves (LCs). Our analysis includes all GRBs with known redshifts detected by the Neil Gehrels {\it Swift} Observatory ({\it Swift}) and the {\it Fermi} Gamma-ray Space Telescope ({\it Fermi}), as well as ground-based optical telescopes. We identify a tight correlation with the $R^2$ coefficient of $\sim 0.89$ for the three-dimensional Dainotti relation between the luminosity at the end of the plateau, its duration measured by {\it Swift}, and the peak luminosity measured by {\it Fermi} in the 10-1000 keV band. When accounting for redshift evolution, we achieve very small intrinsic scatter $\sigma_{int}=0.25\pm0.04$ ($\sim 43\%$ reduction compared to the previous results). Additionally, we explore correlations involving the optical luminosity at the end of the plateau, yielding promising results. We investigate the clustering of different classes of GRBs in the investigated parameter space and discuss its impact on the correlations. Finally, we discuss the theory supporting the evidence of the plateau emission. We present a new paradigm for the GRB plateau: energy extraction from a quickly rotating black hole via spin-down by a magnetically arrested disk. We compare this model with observations and explain multiple observed features. We predict the plateau luminosity - time anti-correlation and discuss the cosmological evolution within this proposed model. Furthermore, within this new model, we discuss the possible physical origin of the clustering of long and short GRBs in the parameter space of plateau luminosity - time - prompt luminosity.

In the era of astronomical big data, more than one million contact binaries have been discovered. Traditional approaches of light curve analysis are inadequate for investigating such an extensive number of systems. This paper builds on prior research to present an advanced Neural Network model combined with the Markov Chain Monte Carlo algorithm and including spot parameters. This model was applied to 12785 contact binaries selected from All-Sky Automated Survey for Supernovae. By removing those with goodness of fit less than 0.8, we obtained the physical parameters of 12201 contact binaries. Among these binaries, 4332 are A-subtype systems, while 7869 are W-type systems, and 1594 systems have mass ratios larger than 0.72 (H-subtype system). A statistical study of the physical parameters was carried out, and we found that there are two peaks in the mass ratio distribution and that the probability of the presence of spot is about 50%. In addition, the differences in flux between the two light maxima are from $-$0.1 to 0.1. As the orbital period and the temperature of the primary component decrease, the difference between the two light maxima becomes more pronounced. Based on the relationships between transfer parameter and luminosity ratio, as well as between luminosity ratio and mass ratio, we found that A-, W-, and H-type contact binaries are distributed in distinct regions.

We analyze the available observational data on the radial distribution of gas and young stellar populations in the disks of low surface brightness (LSB) galaxies and in the outer regions or the extended disks of normal brightness (HSB) galaxies. These cases involve star formation under special conditions of low volume and surface gas density. There is no well-defined boundary between these subgroups of galaxies that we consider, but in non-dwarf LSB galaxies the rate of current star formation within the wide range of radial distances appears to be higher compared to the outer disks of most of HSB galaxies at similar values of the surface gas density. The factors that could stimulate the compression of the rarefied gas at the periphery of galaxies are briefly discussed. Attention is drawn to the idea that the densities of LSB disks estimated from their brightness may be underestimated.

We investigate a dynamical reconstruction of the dark energy equation of state parameter by assuming that it satisfies a law of motion described by an autonomous second-order differential equation, with the limit of the cosmological constant as an equilibrium point. We determine the asymptotic solutions of this equation and use them to construct two families of parametric dark energy models, employing both linear and logarithmic parametrizations with respect to the scale factor. We perform observational constraints by using the Supernova, the Cosmic Chronometers and the Baryon Acoustic Oscillations. The constraint parameters are directly related with the initial value problem for the law of motion and its algebraic properties. The analysis shows that most of the models fit the observational data well with a preference to the models of the logarithmic parametrization. Furthermore, we introduce a new class of models as generalizations of the CPL model, for which the equilibrium point is a constant value rather than the cosmological constant. These models fit the data in a similar or better way to the CPL and the $\Lambda$CDM cosmological models.

Context. Trojan asteroids of Mars date from an early phase of solar system evolution. Based on a sample of <10 asteroids, MT distribution has been previously shown to be both asymmetric and inhomogeneous with population evolution dominated by thermal radiation forces acting over timescales of Gyr. Remarkably, a single asteroid family associated with (5261) Eureka (H~16) at L5 contains most MTs and all members of the stable population fainter than H=18. Aims. Using the currently available sample of MTs and their orbits, we take a fresh look at this population to re-evaluate these earlier conclusions and to search for additional features diagnostic of the evolutionary history of MTs and the Eureka family in particular. Methods. We perform harmonic analysis on numerical time series of the osculating elements to compile a new proper element catalogue comprising 16 L5 and 1 L4 Mars Trojan asteroids. We then combine sample variance analysis with statistical hypothesis testing to identify clusters in the distribution of orbits and assess their significance. Results. We identify two small clusterings significant at 95% confidence of three H=20-21 asteroids each and investigate their likely origin. One of the clusters is probably the result of rotational breakup of a Eureka family asteroid ~10^8 yr ago. The significantly higher tadpole libration width of asteroids in the other cluster is more consistent with an origin as impact ejecta from Eureka itself and on a timescale comparable to the ~1 Gyr age of its family. We further confirm the previously reported correlations in Eureka family orbital distribution attributed to the long-term action of radiation-driven forces and torques on the asteroids.

We analyze JWST NIRSpec+MIRI/MRS observations of the infrared (IR) Polycyclic Aromatic Hydrocarbon (PAH) features in the central regions ($\sim$0.26'' at 6 micron; $\sim$50-440 pc depending on the source) of local luminous IR galaxies. In this work, we examine the effect of nuclear obscuration on the PAH features of deeply obscured nuclei, predominantly found in local luminous IR galaxies, and we compare these nuclei with ``normal'' star-forming regions. We extend previous work to include shorter wavelength PAH ratios now available with the NIRSpec+MIRI/MRS spectral range. We introduce a new diagnostic diagram for selecting deeply obscured nuclei based on the 3.3 and 6.2 micron PAH features and/or mid-IR continuum ratios at $\sim$3 and 5 micron. We find that the PAH equivalent width (EW) ratio of the brightest PAH features at shorter wavelengths (at 3.3 and 6.2 micron) is impacted by nuclear obscuration. Although the sample of luminous IR galaxies used in this analysis is relatively small, we find that sources exhibiting a high silicate absorption feature cluster tightly in a specific region of the diagram, whereas star-forming regions experiencing lower extinction levels occupy a different area in the diagram. This demonstrates the potential of this technique to identify buried nuclei. To leverage the excellent sensitivity of the MIRI imager onboard JWST, we extend our method of identifying deeply obscured nuclei at higher redshifts using a selection of MIRI filters. Specifically, the combination of various MIRI JWST filters enables the identification of buried sources beyond the local Universe and up to z$\sim$3, where other commonly used obscuration tracers such as the 9.7 micron silicate band, are out of the spectral range of MRS. Our results pave the way for identifying distant deeply obscured nuclei with JWST.

During a solar flare, electrons are accelerated to non-thermal energies as a result of magnetic reconnection. These electrons then propagate upwards and downwards from the energy release site along magnetic field lines and produce radio and X-ray emission. On 11 November 2022, an M5.1 solar flare was observed by the Spectrometer/Telescope for Imaging X-rays (STIX) on board Solar Orbiter together with various ground- and space-based radio instruments. The flare was associated with several fine hard X-ray (HXR) structures and a complex set of metric radio bursts (type III, J, and narrowband). By studying the evolution of X-ray, extreme ultraviolet, and radio sources, we aim to study the trajectories of the flare-accelerated electrons in the lower solar atmosphere and low corona. We used observations from the STIX on board Solar Orbiter to study the evolution of X-ray sources. Using radio imaging from the Nançay Radio heliograph (NRH) and the Newkirk density model, we constructed 3D trajectories of 14 radio bursts. Imaging of the HXR fine structures shows several sources at different times. The STIX and NRH imaging shows correlated changes in the location of the HXR and radio source at the highest frequency during the most intense impulsive period. Imaging and 3D trajectories of all the bursts show that electrons are getting accelerated at different locations and along several distinct field lines. The longitude and latitude extent of the trajectories are ~30 arcsec and ~ 152 arcsec. We find that the electrons producing HXR and radio emission have similar acceleration origins. Importantly, our study supports the scenario that the flare acceleration process is temporally and spatially fragmentary, and during each of these small-scale processes, the electron beams are injected into a very fibrous environment and produce complex HXR and radio emission.

In the current galaxy formation paradigm, collisions play a crucial role. A fraction of galaxy collisions results in flyby events, and a galaxy that has passed through another galaxy is called a backsplash galaxy. Such flyby events are of particular interest for explaining the quenching of isolated galaxies. One signature of backsplash galaxies is that they have high velocities relative to their environment, since they do not move in the same flow as surrounding galaxies. This feature can be studied in simulations, but it is also useful to have a theory that can predict the velocities of backsplash galaxies. In this paper, we develop such a theory based on the Zel'dovich approximation and use it to determine the maximal expected velocity of a backsplash galaxy in a given volume.

Protoplanetary disk evolution exhibits trends with stellar mass as well as diversity of structure and lifetime, with implications for planet formation and demographics. We show how varied outcomes can result from evolving structures in the inner disk that attenuate stellar soft X-rays that otherwise drive photoevaporation in the outer disk. The magnetic truncation of the disk around a rapidly rotating T Tauri star is initially exterior to the co-rotation radius and``propeller" accretion is accompanied by an inner magnetized wind, shielding the disk from X-rays. Stellar contraction while rotation is ``locked" to the disk causes the truncation radius to migrate interior to the co-rotation radius, the inner wind to disappear, and photoevaporation to erode a gap in the disk, accelerating its dissipation. This X-ray attenuation scenario explains the trend of longer lifetime, reduced structure and compact size of disks around lower-mass stars. It also explains an observed lower bound and scatter in the distribution of disk accretion rates. Disks that experience early photoevaporation and form gaps can efficiently trap solids at a pressure bump at 1--10 au, triggering giant planet formation, while those with later-forming or no gaps form multiple smaller planets on close-in orbits, a pattern consistent with observed exoplanet demographics.

Context. Substellar objects, including brown dwarfs and free-floating planetary-mass objects, are a significant product of star formation. Their sensitivity to initial conditions and early dynamical evolution makes them especially valuable for studying both planetary and stellar formation processes. Aims. We search for brown dwarfs and isolated planetary mass objects in a young star forming region to better constrain their formation mechanisms. Methods. We take advantage of Euclid unprecedented sensitivity, spatial resolution and wide-field of view to search for brown dwarfs and free-floating planetary mass objects in the LDN 1495 region of the Taurus molecular clouds. We combine the recent Euclid Early Release Observations with older very deep ground-based images obtained over more than 20 yr to derive proper motions and multi-wavelength photometry, and select members based on their morphology and their position in a proper motion diagram and in 9 color-magnitude diagrams. Results. We identify 15 point sources with proper motions, colors, and luminosity consistent with being members of LDN 1495. Six of these objects were already known M9-L1 members. The remaining nine are newly identified sources that could have spectral types ranging from late-M to early-T with masses potentially as low as 1-2 MJup based on their luminosity and according to evolutionary models. However, follow-up observations are needed to confirm their nature, spectral type and membership. If extrapolated to theentire Taurus star forming region, this result suggests the potential presence of several dozen free-floating planetary mass objects.

Context. As the Solar System orbits the Milky Way, it encounters various Galactic environments, including dense regions of the interstellar medium (ISM). These encounters can compress the heliosphere, exposing parts of the Solar System to the ISM, while also increasing the influx of interstellar dust into the Solar System and Earth's atmosphere. The discovery of new Galactic structures, such as the Radcliffe wave, raises the question of whether the Sun has encountered any of them. Aims. The present study investigates the potential passage of the Solar System through the Radcliffe wave gas structure over the past 30 million years (Myr). Methods. We used a sample of 56 high-quality, young ($\leq$ 30 Myr) open clusters associated with a region of interest of the Radcliffe wave to trace its motion back and investigate a potential crossing with the Solar System's past orbit. Results. We find that the Solar System's trajectory intersected the Radcliffe wave in the Orion region. We have constrained the timing of this event to between 18.2 and 11.5 Myr ago, with the closest approach occurring between 14.8 and 12.4 Myr ago. Notably, this period coincides with the Middle Miocene climate transition on Earth, providing an interdisciplinary link with paleoclimatology. The potential impact of the crossing of the Radcliffe wave on the climate on Earth is estimated. This crossing could also lead to anomalies in radionuclide abundances, which is an important research topic in the field of geology and nuclear astrophysics.

Dense cores in massive, parsec-scale molecular clumps are sites that harbor protocluster formation. We present results from observations towards a hub-filament structure of a massive Infrared Dark Cloud (IRDC) G14.225-0.506 using the Atacama Large Millimeter/submillimeter Array (ALMA). The dense cores are revealed by the 1.3 mm dust continuum emission at an angular resolution of $\sim$ 1.5'' and are identified through the hierarchical Dendrogram technique. Combining with the N$_2$D$^+$ 3-2 spectral line emission and gas temperatures derived from a previous NH$_3$ study, we analyze the thermodynamic properties of the dense cores. The results show transonic and supersonic-dominated turbulent motions. There is an inverse correlation between the virial parameter and the column density, which implies that denser regions may undergo stronger gravitational collapse. Molecular outflows are identified in the CO 2-1 and SiO 5-4 emission, indicating active protostellar activities in some cores. Besides these star formation signatures revealed by molecular outflows in the dense cores, previous studies in the infrared, X-ray, and radio wavelengths also found a rich and wide-spread population of young stellar objects (YSOs), showing active star formation both inside and outside of the dense cloud.

Storms are emerging as key drivers in shaping hydrogen-dominated atmospheres. Trace gas condensation can suppress convection and disrupt the distribution of energy and material in hydrogen atmospheres. On Jupiter, the presence of water has been invoked to control the occurrence of large-scale storms; however, the impact of storms on the ammonia and temperature distribution is unknown. We use Juno Microwave Radiometer observations of a large-scale storm in 2017 to study the aftermath of such a storm on the atmosphere. Anomalies in the retrieved ammonia abundance and atmospheric temperature show how storms deplete and heat the upper atmosphere while simultaneously depositing material well below the layers they were triggered at. These observations, aided by simulations, show that the water and ammonia cycles are coupled and that their combined effect plays a key role in explaining the depletion of ammonia in the tropospheres of Jupiter and Saturn.

The chemical compositions of evolved stars in Local Group dwarf spheroidal galaxies (dSphs) provide insight into the galaxy's past star formation and nucleosynthesis. Neutron-capture element abundances are especially interesting. In particular, $s$-process elements can provide a third chemical clock for resolving star formation histories in addition to core collapse and Type Ia supernovae. Likewise, the primary sites of the $r$-process are still areas of extensive research. Until now, the number of stars with neutron-capture element abundances in dSphs has been limited by the need for stars bright enough for high-resolution spectroscopy. We present abundance measurements of the neutron-capture elements Sr, Y, Ba, and Eu with errors $<$ 0.4 dex - as well as new measurements of Mg - in 491 stars in Sculptor, Fornax, Draco, Sextans, and Ursa Minor. The large number of stars in our sample is possible because we used medium-resolution spectra from the DEIMOS spectrograph, assembling the largest homogeneous set of neutron-capture abundances in dwarf spheroidal galaxies to date. By utilizing the abundances of both $s$- and $r$-process elements, we find evidence of an $s$-process contribution at early times in Sculptor from our measurements of [Ba/Fe]. This is a potential signature of $s$-process nucleosynthesis in fast-rotating massive stars. By comparing our measurements of [Eu/Fe] with [Mg/Fe], we show the need for an $r$-process source that has a short delay time to enrich stars in the dSphs. Thus, neutron star mergers are likely not the sole source of $r$-process material in dSphs.

Asteroid shape inversion using photometric data has been a key area of study in planetary science and astronomical this http URL, the current methods for asteroid shape inversion require extensive iterative calculations, making the process time-consuming and prone to becoming stuck in local optima. We directly established a mapping between photometric data and shape distribution through deep neural networks. In addition, we used 3D point clouds to represent asteroid shapes and utilized the deviation between the light curves of non-convex asteroids and their convex hulls to predict the concave areas of non-convex asteroids. We compared the results of different shape models using the Chamfer distance between traditional methods and ours and found that our method performs better, especially when handling special shapes. For the detection of concave areas on the convex hull, the intersection over union (IoU) of our predictions reached 0.89. We further validated this method using observational data from the Lowell Observatory to predict the convex shapes of the asteroids 3337 Milo and 1289 Kuta, and conducted light curve fitting experiments. The experimental results demonstrated the robustness and adaptability of the method

China plans to launch a probe (Tianwen-2) around 2025, mainly for exploring the near-Earth asteroid 2016 HO3 . The mission involves close-range exploration, landing, and mining operations that require three-dimensional modeling of the asteroid, which requires prior knowledge of its material composition and uniformity. This information is crucial in progressive or ground exploration processes. Our research focuses on high-precision intelligent inversion of complex physical properties of asteroids based on spectral data, providing support for further analysis of aster oid materials, density, and structure. We have developed a platform for asteroid spectral classification, albedo estimation, and composition analysis, which includes three types of neural networks based on Transformer attention mechanism: One for spectral classification, achieving a four-class classification accuracy of 97.28% and an eleven-class classification accuracy of 95.69%; second one for albedo estimation, with an average absolute error of 0.0308 in S-type asteroid albedo estimation, and the third one for composition analysis, with a predicted spectral angular distance of only 0.0340 and a root mean square error of 0.1759 for the abundance of end members. These results indicate that our network can provide high-precision asteroid spectral classification, albedo estimation, and composition analysis results. In addition, we utilized the platform to analyze and provide results for six asteroids.

Galaxy clusters act as reservoirs of high-energy cosmic rays (CRs). As CRs propagate through the intracluster medium, they generate diffuse $\gamma$-rays detectable by arrays such as LHAASO. These $\gamma$-rays result from proton-proton ($pp$) collisions of very high-energy cosmic rays (VHECRs) or inverse Compton (IC) scattering of positron-electron pairs created by $p\gamma$ interactions of ultra-high-energy cosmic rays (UHECRs). We analyzed diffuse $\gamma$-ray emission from the Coma, Perseus, and Virgo clusters using LHAASO data. Diffuse emission was modeled as a disk of radius $R_{500}$ for each cluster while accounting for point sources. No significant diffuse emission was detected, yielding 95\% confidence level (C.L.) upper limits on the $\gamma$-ray flux: for WCDA (1-25~TeV) and KM2A ($>25$~TeV), less than $(49.4, 13.7, 54.0)$ and $(1.34, 1.14, 0.40) \times 10^{-14}$~ph~cm$^{-2}$~s$^{-1}$ for Coma, Perseus, and Virgo, respectively. The $\gamma$-ray upper limits can be used to derive model-independent constraints on the integral energy of CRp above 10~TeV (corresponding to the LHAASO observational range $>1$~TeV under the $pp$ scenario) to be less than $(1.96, 0.59, 0.08) \times 10^{61}$~erg. The absence of detectable annuli/ring-like structures, indicative of cluster accretion or merging shocks, imposes further constraints on models in which the UHECRs are accelerated in the merging shocks of galaxy clusters.

This research paper has light curve modelling, spectroscopy, and detailed asteroseismic studies for five detached eclipsing binaries with a $\gamma$ Dor component that have been found so far by the sky surveys. The objective is to study the pulsational characteristics of the oscillating stars of the systems, as well as to estimate their absolute parameters and their relation to the pulsational frequencies and enrich the sample size of this kind of system, which has been relatively poor up to today. The physical properties of these systems are compared with other similar cases and the locations of their components are plotted in the Hertzsprung-Russell (HR) diagrams. In order to discover the pulsation frequencies, the photometric data are analysed using eclipsing binary modelling methods and their residuals are subjected to Fourier analysis. Finally, we evaluated that there is a possible relationship between the pulsation periods of the primaries and the orbital period ($P_{\rm orb}$), the force exerted per gram on the surface of the pulsating star ($f$/M$_{\rm puls}$), and the fractional radius of the pulsating star ($r_{\rm puls}$) of the systems, using a sample of 39 eclipsing binary systems with $\gamma$ Dor type primaries. We utilised them to construct a unique observational $\gamma$ Dor instability strip, delineated by a lower limit of 6,850 K at the red edge and an upper limit of 7,360 K at the blue edge on the zero-age main sequence. The majority of 39 pure $\gamma$ Dor stars are located within the region that is covered by this observational strip.

The apparent slipping motion of flare loops is regarded as a key feature of the 3D magnetic reconnection in the solar flares. The slippage with a super-Alfvénic speed could be defined as slipping-running reconnection while the slippage with a sub-Alfvénic speed is called slipping reconnection. Due to the limitation of the observational instrument temporal resolution, the apparent slippage of the flare loop footpoints along the flare ribbons with super-Alfvénic speed is quite rare to our knowledge. In this paper, we report a unique event that exhibits not only the sub-Alfvénic slippage, but also the quasiperiodic super-Alfvénic slippage of ribbon substructures during a C3.4-class flare (SOL2023-01-18-T15:23), using the high temporal resolution observations of the Interface Region Imaging Spectrograph ($\sim$2 s). The super-Alfvénic slippage with a speed of up to $\sim$ 1688 km s$^{-1}$ is directly observed in this study. The calculated period of the apparent super-Alfvénic slippage in both ribbons is between 8.4 and 11.9 seconds. This work provides the first observational evidence of the periodicity for the slipping-running magnetic reconnection.

The non-detection of periodicity related to rotation challenges the magnetar model for fast radio bursts (FRBs). Moreover, a bimodal distribution of the burst waiting times is widely observed in hyper-active FRBs, a significant deviation from the exponential distribution expected from stationary Poisson processes. By combining the epidemic-type aftershock sequence (ETAS) earthquake model and the rotating vector model (RVM) involving the rotation of the magnetar and orientations of the spin and magnetic axes, we find that starquake events modulated by the rotation of FRB-emitting magnetar can explain the bimodal distribution of FRB waiting times, as well as the non-detection of periodicity in active repeating FRBs. We analyze data from multiple FRB sources, demonstrating that differences in waiting time distributions and observed energies can be explained by varying parameters related to magnetar properties and starquake dynamics. Our results suggest that rotation-modulated starquakes on magnetars can possibly be a unified source for FRBs. Notably, we find that active repeaters tend to have small magnetic inclination angles in order to hide their periodicity. We also show that our model can reproduce the waiting time distribution of a pulsar phase of the galactic magnetar SGR J1935+2154 with a larger inclination angle than the active repeaters, which could explain the detection of spin period and the relatively low observed energy for FRBs from the magnetar. The spin periods of active repeaters are not well constrained, but most likely fall in the valley region between the two peaks of the waiting time distributions.

This thesis investigates the impact of dark matter on neutron star properties, focusing on mass, radius, and tidal deformability. Using two-fluid and single-fluid models, dark matter is incorporated into the equation of state (EOS) via a Relativistic Mean Field (RMF) approach. The study finds that increasing dark matter content reduces the maximum mass, radius, and tidal deformability. Bayesian inference, supported by LIGO-Virgo gravitational wave data and NICER mass-radius measurements, refines these models. Despite dark matter's influence, the semi-universal C-Love relation remains valid. Machine learning techniques effectively classify dark matter-admixed neutron stars. The thesis also explores a sigma-cut potential in the EOS, which stiffens the EOS at high densities, favoring larger radii and lower f-mode frequencies. The study of non-radial oscillations, particularly f- and p-modes, highlights their sensitivity to neutron star composition and EOS. These findings enhance our understanding of neutron star interiors and dark matter's role, emphasizing the need for further observational and theoretical advancements.

We present a study of the dynamics of multi-component models of spiral galaxies at various stages of grand merging. Numerical models include a self-consistent account of the dynamics of collisionless stellar subsystems and N-body dark matter, as well as gaseous components. Calculation of gas heating and cooling processes allows us to consider a wide temperature range from 80 to 100 thousand degrees. Using the Smoothed-particle hydrodynamics method to solve hydrodynamic equations makes it possible to track the evolution of the gas of each galaxy, calculating the content of gas components of each object in the process of complex interchange of matter. The gravitational interaction is determined in a direct way by summing the contributions according to Newton's law, which minimizes the modeling error. This requires significant computing resources using graphics accelerators on hybrid computing platform CPU + multi-GPUs. The study aims to reconcile theoretical models with the morphology and kinematics of a number of observed systems. In particular, Taffy-type objects are considered, where two galaxies are connected by a gas bridge with a characteristic small-scale gas structure after the disks pass through each other in an approximately flat orientation. Examples of such systems are the observed pairs of galaxies UGC 12914/UGC 12915, NGC 4490/NGC 4485, UGC813/UGC816 etc.

We present an analysis of interplanetary scintillation (IPS) observations conducted with the Arecibo 305-m radio telescope during the minimum phase at the end of solar cycle 24 and the onset of solar cycle 25. These observations span a broad frequency range of ~300 to 3100 MHz, encompassing the P-, L-, and S-bands, and covered heliocentric distances from ~5 to 200 solar radii. The dynamic spectrum of the scintillations obtained at L-band shows a systematic decrease in the scintillation index from the lowest to the highest frequency, offering valuable insight into the influence of the solar wind density microstructures responsible for scintillation. Analyses of the scintillation index ($m$) for multiple sources at L-band, along with near-simultaneous observations of selected sources covering the P-, L-, and S-bands, clearly demonstrate a wavelength dependence of $m \propto \lambda^\omega$, which inherently leads to a dependence of $m$ on the Fresnel scale, when considering the effective distance to the scattering screen, $z$. The index $\omega$ ranges between $\sim$1 and 1.8. The average $\omega$ value of a source, determined from observations made on different days, exhibits variability across sources. The results on the radial dependence of scintillation agree with earlier IPS measurements. The temporal power spectra obtained over the wide frequency range exhibit a power-level evolution in accordance with the wavelength dependence, and a broadening with increasing observation frequency. Furthermore, the increased temporal-frequency rounding of the `Fresnel knee' in the spectrum with the observing frequency suggests a novel phenomenon: an increase in anisotropy as the scale size of the density-turbulence structure decreases.

In astrochemical exploration, infrared (IR) spectroscopy is vital for understanding the composition and structure of ice in various space environments. This article explores the impact of incident angles on IR spectroscopy, focusing on molecular components present in interstellar and circumstellar ice mantles such as CO, CO$_2$, H$_2$O, CH$_3$OH, NH$_3$, CH$_4$, H$_2$S. The experiment involves changing the angle at which the infrared beam hits the surface used for ice deposition. It is important to measure the density of the ice layer accurately, especially for experiments that involve using different angles in infrared spectroscopy. Furthermore, the experimental methodology allowed us to derive the {\it effective} refraction index values in the infrared range for each ice component. Existing corrections typically consider geometric configurations but overlook the refractive index of the ice ($n$), a factor dependent on ice composition. The study reveals that the incident angle and the refractive index, determine the pathlength of the IR beam across the ice sample. This insight challenges conventional corrections, impacting the integrated absorption values of the IR bands and column densities. In addition, for certain ice components, variations in the incidence angle affect the longitudinal (LO) and transverse (TO) optical modes of the ice, leading to observable changes in the IR band profiles that provide information on the amorphous or crystalline structure of the ice. The practical implications of this work apply to experimental setups where normal IR measurements are unfeasible. Researchers using, for example, the standard 45$^{\circ}$ angle for IR spectroscopy, will benefit from a more accurate estimation of ice column density.

With a new probabilistic technique for sampling interstellar object (ISO) orbits with high efficiency, we assess the observability of ISOs under a realistic cadence for the upcoming Vera Rubin Observatory's Legacy Survey of Space and Time (LSST). Using the Ōtautahi-Oxford population model, we show that there will be complex on-sky structure in the pattern of direction and velocity revealed by the detected ISO population, with the expected enhanced northern flux complicating efforts to derive population parameters from the LSST's predominately southern footprint. For luminosity functions with slopes of $2.5\leq q_s\leq 4.0$, the most discoverable ISOs have $H_r\simeq 14.6-20.7$; for previously estimated spatial densities, between 6 and 51 total ISOs are expected. The slope of the luminosity function of ISOs will be relatively quickly constrained. Discoveries are evenly split around their perihelion passage and are biased to lower velocities. After their discovery by LSST, it will be rare for ISOs to be visible for less than a month; most will have $m_r \leq 23$ for months, and the window for spectroscopic characterization could be as long as two years. These probabilistic assessments are robust against model or spatial density refinements that change the absolute numbers of ISO discoveries.

We use FUV spectra of 36 T Tauri stars, predominately from $\textit{Hubble Space Telescope}$'s ULLYSES program, to examine the kinematic properties of fluorescent H$_2$ emission lines for evidence of disk outflows. Leveraging improvements to the $\textit{HST}$-COS wavelength solution, we co-add isolated lines within four fluorescent progressions ([$\textit{v'}$,$\textit{J'}$] = [1,4], [1,7], [0,2], and [3,16]) to improve signal-to-noise, and we fit each co-added line profile with one or two Gaussian components. Of the high S/N line profiles (S/N $\geq$ 12 at the peak of the profile), over half are best fit with a combination of a broad and narrow Gaussian component. For profiles of the [1,4] and [1,7] progressions, we find a systematic blue-shift of a few km s$^{-1}$ between the broad and narrow centroid velocities and stellar radial velocities. For the [0,2] progression, we find centroid velocities consistently blueshifted with respect to stellar radial velocities on the order of -5 km s$^{-1}$ for the single and narrow components, and -10 km s$^{-1}$ for the broad components. Overall, the blueshifts observed in our sample suggest that the molecular gas traces an outflow from a disk wind in some sources, and not solely disk gas in Keplerian rotation. The low-velocity systematic blue-shifts, and emitting radii as inferred from line FWHMs, observed in our sample are similar to those observed with optical [O I] surveys of T Tauri stars. We estimate H$_2$ mass-loss rates of 10$^{-9}$ to 10$^{-11}$ $M_{\odot}$ yr$^{-1}$, but incomplete knowledge of wind parameters limits comparisons to global models.

Photointeractions of ultra-high energy cosmic rays (UHECRs) in astro-physical scenarios are in general of stochastic nature and are often modeled with Monte Carlo methods to obtain the form of the distributions resulting from a sequence of interactions. These distributions are non trivial because the products resulting from each interaction as well as the number and distances covered by the secondary nuclear species are all random. In this work, a stochastic approach based on the theory of matrix exponential distributions is employed to describe the cascade distributions analytically and illustrate their potential for tracing the individual history of UHECRs, including inside the source. This analytic description has the advantage of better precision and considerably reduced computational cost in contrast to Monte Carlo codes, while requiring the same inputs: the interaction rates, the multiplicity, and the energy distributions of secondaries from a single interaction. The description of the composition evolution from in-source to extragalactic propagation (currently performed in separate simulations in the literature) is achieved here as a continuous distribution, using a gamma-ray burst scenario inspired by the event GRB170817A. Finally, the potential for locating a source based on the reconstructed UHECR origin employing this description is discussed under simplified general assumptions.

We study the potential bias in the identification of fast radio burst (FRB) host galaxies due to radio localisation uncertainty. Using a sample of FRBs localised to typically 0.5'' by the Australian Square Kilometre Array Pathfinder (ASKAP), we artificially increase the localisation uncertainty up to 10'', and re-run the Probabalistic Association of Transients to their Hosts (PATH) algorithm to determine the most likely host galaxy. We do not find evidence of a significant change in identified hosts until the localisation precision is worsened to 2'' or greater.

Image denoising based on deep learning has witnessed significant advancements in recent years. However, existing deep learning methods lack quantitative control of the deviation or error on denoised images. The neural networks Self2Self is designed for denoising single-image, training on it and denoising itself, during which training is costly. In this work we explore training Self2Self on an astronomical image and denoising other images of the same kind, which is suitable for quickly denoising massive images in astronomy. To address the deviation issue, the abnormal pixels whose deviation exceeds a predefined threshold are restored to their initial values. The noise reduction includes training, denoising, restoring and named TDR-method, by which the noise level of the solar magnetograms is improved from about 8 G to 2 G. Furthermore, the TDR-method is applied to galaxy images from the Hubble Space Telescope and makes weak galaxy structures become much clearer. This capability of enhancing weak signals makes the TDR-method applicable in various disciplines.

We examine the effects of thermal conduction on relativistic, magnetized, viscous, advective accretion flows around rotating black holes considering bremsstrahlung and synchrotron cooling processes. Assuming the toroidal component of magnetic fields as the dominant one, we self-consistently solve the steady-state fluid equations to derive the global transonic accretion solutions for a black hole of spin $a_{\rm k}$. Depending on the model parameters, the magnetized accretion flow undergoes shock transitions and shock-induced global accretion solutions persist over a wide range of model parameters including the conduction parameter ($\Upsilon_{\rm s}$), plasma-$\beta$, and viscosity parameter ($\alpha_{\rm B}$). We find that the shock properties -- such as shock radius ($r_{\rm s}$), compression ratio ($R$), and shock strength ($S$) -- are regulated by $\Upsilon_{\rm s}$, plasma $\beta$, and $\alpha_{\rm B}$. Furthermore, we compute the critical conduction parameter ($\Upsilon_{\rm s}^{\rm cri}$), a threshold beyond which shock formation ceases to exist, and investigate its dependence on plasma-$\beta$ and $\alpha_{\rm B}$ for both weakly rotating ($a_{\rm k} \rightarrow 0$) and rapidly rotating ($a_{\rm k} \rightarrow 1$) black holes. Finally, we examine the spectral energy distribution (SED) of the accretion disc and observe that increased thermal conduction and magnetic field strength lead to more luminous emission spectra from black hole sources.

Jets emanating from active galactic nuclei (AGN) represent some of the most formidable particle accelerators in the universe, thereby emerging as viable candidates responsible for the detection of ultra-high-energy cosmic rays (UHECRs). If AGN jets indeed serve as origins of UHECRs, then the diffuse flux of these cosmic rays would be dependent on the power and duty cycle of these jets, which are inherently connected to the nature of black hole accretion flows. In this article, we present our cosmological semi-analytic framework, JET(Jets from Early Times), designed to trace the evolution of jetted AGN populations. This framework serves as a valuable tool for predictive analyses of cosmic ray energy density and, potentially, neutrino energy density. By using JET, we model the formation and evolution of galaxies and supermassive black holes (SMBHs) from z = 20 to z = 1, incorporating jet formation and feedback mechanisms and distinguishing between various accretion states determined by the SMBH Eddington ratios. The implications of different SMBH growth models on predicting cosmic ray flux are investigated. We provide illustrative examples demonstrating how the associated diffuse UHECR fluxes at the source may vary in relation to the jet production efficiencies and the selected SMBH growth model, linking cosmological models of SMBH growth with astroparticle backgrounds.

We present TESSELLATE, a dedicated pipeline for performing an untargeted search documenting all variable phenomena captured by the TESS space telescope. Building on the TESSreduce difference imaging pipeline, TESSELLATE extracts calibrated and reduced photometric data for every full frame image in the TESS archive. Using this data, we systematically identify transient, variable and non-sidereal signals across timescales ranging from minutes to weeks. The high cadence and wide field of view of TESS enables us to conduct a comprehensive search of the entire sky to a depth of ~17 $m_i$. Based on the volumetric rates for known fast transients, we expect there to be numerous Fast Blue Optical Transients and Gamma Ray Burst afterglows present in the existing TESS dataset. Beyond transients, TESSELLATE can also identify new variable stars and exoplanet candidates, and recover known asteroids. We classify events using machine learning techniques and the work of citizen scientists via the Zooniverse Cosmic Cataclysms project. Finally, we introduce the TESSELLATE Sky Survey: a complete, open catalog of the variable sky observed by TESS.

Sunspot engravings and measurements in 1660-1676 are analyzed to retrieve sunspot area and heliocoordinates. Based on these data, we revise the Hoyt and Schatten (The role of the sun in climate change, 1997) hypothesis of long-lived sunspots during the Maunder minimum as a sign of weakened convection. Historical reports also clarify what each observer defined as a sunspot and the purpose of the observations. The reconstructed longitudes of sunspots allow us to evaluate the rotation rate, revealing that the historical rotation profile resembles that of long-lived sunspot groups in the modern era.

We present a statistical study of the neutral atomic hydrogen (HI) gas extending into the circumgalactic medium perpendicular to the disk for 7 edge-on galaxies with inclinations above $85^{\circ}$ from the FEASTS program with a $3\sigma$ ($20\,\text{km}\,\text{s}^{-1}$) column density ($N_{\text{HI}}$) depth of $5\times10^{17} \text{cm}^{-2}$. We develop two photometric methods to separate the extraplanar HI from the disk component, based on existing interferometric data and parametric modeling of the disk flux distribution respectively. With both methods, the FEASTS data exhibit clear extended wings beyond the disk along the minor axis. The extraplanar HI accounts for 5% to 20% of the total HI mass and extends to $20\text{-}50$ kpc at $N_{\text{HI}}=10^{18} \text{cm}^{-2}$. We find a tight positive correlation between vertical extensions of the extraplanar HI and total HI mass $M_\text{HI}$. The iso-density shape of HI at $N_{\text{HI}}=10^{18} \text{cm}^{-2}$ has an average axis ratio of $0.56\pm0.11$. The off-disk $N_{\text{HI}}$ profiles of these edge-on galaxies well represent the lower envelop of previous Lyman-$\alpha$ absorption measurements at low-redshift. Our results suggest that at $N_{\text{HI}}=5\times10^{17} \text{cm}^{-2}$, the HI extends considerably further than the known thin and thick disks in the vertical direction, but still remains much flattener than a spherical distribution, consistent with theoretical expectations that outflow, circulation, and accretion should have different impacts in these two directions. We show the tension of our results with Illustris and TNG predictions, highlighting the constraining power of our results for future simulations.

Collisionless shocks are one of the most powerful particle accelerators in the Universe. In the heliosphere, type II solar radio bursts are signatures of electrons accelerated by collisionless shocks launched at the Sun. Spectral observations of these bursts show a variety of fine structures often composing multiple type II lanes. The origin of these lanes and structures is not well understood and has been attributed to the inhomogeneous environment around the propagating shock. Here, we aim to determine the large-scale local structures near a coronal shock wave using high-resolution radio imaging observations of a complex type II radio burst observed on 3 October 2023. By using inteferometric imaging from the Low Frequency Array (LOFAR), combined with extreme ultraviolet observations, we investigate the origin of multiple type II lanes at low frequencies (30--80~MHz) relative to the propagating shock wave. We identify at least three radio sources at metric wavelengths corresponding to a multi-lane type II burst. The type II burst sources propagate outwards with a shock driven by a coronal mass ejection. We find a double radio source that exhibits increasing separation over time, consistent with the expansion rate of the global coronal shock. This suggests that the overall shock expansion is nearly self-similar, with acceleration hotspots forming at various times and splitting at a rate proportional to the shock's expansion. Our results show the importance of increased spatial resolution in determining either the small-scale spatial properties in coronal shocks or the structuring of the ambient medium. Possible shock corrugations or structuring of the upstream plasma at the scale of 10$^5$~km can act as hotspots for the acceleration of suprathermal electrons.

According to the standard inside-out galaxy formation scenario, galaxies first form a dense core and then gradually assemble their outskirts. This implies that galaxies with similar central stellar mass densities might have evolutionary links. We use the UVJ color-color diagram to select quiescent galaxies in the redshift interval from 0.5 to 2.5 and classify them into different subsamples based on their central stellar mass densities, stellar mass, morphological type and redshift. We then infer the intrinsic axis ratios $\mu_{B/A}$ and $\mu_{C/A}$ of different subsamples based on the apparent axis ratio $q$ distributions, where A, B, and C refers to, respectively, the major, intermediate and minor axis of a triaxial ellipsoidal model. We find that 1) massive quiescent galaxies have typical intrinsic shapes similarly close to thick oblate structures, with $\mu_{B/A} \gtrsim 0.9$, regardless of stellar mass, redshift, or central stellar mass densities, and 2) galaxies at higher redshift are systematically thinner than their lower-redshift counterparts, and 3) when splitting the sample into early type and late type with Sersic indices, ETGs at higher redshift are slightly more prolate (smaller average $\mu_{B/A}$) than those at lower redshift. Minor mergers of galaxies may have played important roles in the structural evolution of quiescent galaxies found in this work.

The ultra-relativistic jet released in gamma-ray bursts (GRBs) is expected to produce ultra-high-energy cosmic rays (UHECRs), prompt gamma-ray emission and hence prompt high-energy neutrinos by photopion interactions. In this work, we calculate the time-integrated neutrino spectrum during the expansion of jets by taking into account the time evolution of cosmic ray and secondary spectra and neutrino production. We numerically solve the continuity equations for nucleons, pions, and muons for their spectral evolution. Since pion and muon damping weakens as the jet expands, the neutrino production at large radii at high energies may dominate that around the jet energy dissipation radius. Compared with the usually adopted approaches that only consider neutrino production around the energy dissipation radius, the overall UHE neutrino fluence integrated over time can be significantly larger, and the flavor fraction of electron neutrinos as function of neutrino energy is different at UHE, due to neutrino production at radii much larger than the energy dissipation radius. The faster magnetic field decay leads to larger UHE neutrino fluence, and the UHE neutrino spectra is weakly dependent on the energy dissipation radius and the jet Lorentz factor. Observations of prompt EeV neutrinos from GRBs by the next-generation neutrino telescopes, e.g., GRAND and IceCube-Gen2, will test the hypothesis of GRBs as UHECR sources and probe the physics of GRB jets.

GRS 1758-258 and 1E 1740.7-2942 are two long-known persistent black hole binaries in the Galactic Center region. Using INTEGRAL's extensive monitoring of the Galactic Center and Bulge, we studied their temporal and spectral evolutions in the 30-610 keV energy range from March 2003 through April 2022 with the IBIS/ISGRI gamma-ray telescope. Our analyses found that the sources typically had Comptonized spectra, though not always with the same parameters. The spectral states with more than 8 Ms of observation time show deviations from a Comptonized spectrum above ~200 keV or a "hard tail" that extends up to at least 600 keV. The origin of this component remains debated with the most popular scenarios being synchrotron emission from the jet or Comptonization in a hybrid thermal/non-thermal plasma. Anyway, the GRS 1758-258 and 1E 1740.7-2942 spectra are acceptably described by CompTT+po (jet) and Eqpair (hybrid Comptonization) scenarios. To differentiate between the two scenarios, we calculated the Spearman correlation coefficient comparing 30-50 keV count rates with those in higher energy bands (50-100, 100-300, and 300-600 keV). The count rates below 300 keV are strongly correlated, indicating those photons arise from the same physical process. Above 300 keV the count rates are either anti-correlated or not correlated with the 30-50 keV count rates for GRS 1758-258, which suggests that the photons originate from a different physical process. For 1E 1740.7-2942, the level of correlation is unclear due to scatter in the data points. However, the 300-600 keV count rates are consistent with a constant value. This disfavors the hybrid Comptonization scenario for both sources.

IceCube measures a diffuse neutrino flux comparable to the Waxman-Bahcall bound, which suggests the possibility that the ultra-high energy cosmic rays (UHECRs) have a common origin with diffuse high energy neutrinos. We propose high energy gamma-ray and/or neutrino observations toward the arrival directions of UHECRs to search for the sources and test this possibility. We calculate the detection probability of gamma-ray/neutrino sources, and find that the average probability per UHECR of >10 EeV is $\sim$10% if the sensitivity of the gamma-ray or neutrino telescope is $\sim$10$^{-12}$ erg cm$^{-2}$s$^{-1}$ and the source number density is $\sim$10$^{-5}$ Mpc$^{-3}$. Future gamma-ray and neutrino observations toward UHECRs, e.g., by LHAASO-WCDA, CTA, IceCube/Gen2, are encouraged to constrain the density of UHECR sources or even identify the sources of UHECRs.

We report the discovery of a kpc-scale lacy filamentary structure in low surface-brightness H$\alpha$ and [OIII]$\lambda5007$ emission around the low-redshift, extremely metal-poor, and compact reionisation-era analogue SDSS J1044+0535. We identify seven elliptical arcs in H$\alpha$ emission at $\mathrm{SB}_\mathrm{H\alpha} \sim 1 - 2 \times 10^{-18}$erg\,s$^{-1}$cm$^{-2}$arcsec$^{-2}$. These arcs extent 6-9 kpc outwards from the star-forming complexes. We interpret these features as limb-brightened giant-shells that bound egg-shaped kpc-scale bubbles. These shells are five to eight times larger than the known giant-shells around nearby star-forming dwarfs. Kinematic maps reveal a gradient perpendicular to the major-axis and line broadening in the outskirts. The latter, when interpreted due to line-splitting of the expanding shells, provides an estimate of the kinetic energy of the bubbles that can be well reconciled with them being blown by the mass dominating 15-20 Myr old stellar population.

Next-generation galaxy surveys promise unprecedented precision in testing gravity at cosmological scales. However, realising this potential requires accurately modelling the non-linear cosmic web. We address this challenge by exploring conditional generative modelling to create 3D dark matter density fields via score-based (diffusion) and flow-based methods. Our results demonstrate the power of diffusion models to accurately reproduce the matter power spectra and bispectra, even for unseen configurations. They also offer a significant speed-up with slightly reduced accuracy, when flow-based reconstructing the probability distribution function, but they struggle with higher-order statistics. To improve conditional generation, we introduce a novel multi-output model to develop feature representations of the cosmological parameters. Our findings offer a powerful tool for exploring deviations from standard gravity, combining high precision with reduced computational cost, thus paving the way for more comprehensive and efficient cosmological analyses

Using 6 Milky Way analogues with two different numerical resolutions from the Auriga simulation, we investigate the total mass, spatial distribution and kinematics of the earliest stars relics in the Milky Way at $z=0$. These relics (second generation stars) formed over a wide redshift range, from about $z=22$ to $z=4$, with an average formation redshift of $z \sim 10.0$, and comprise about $2\times10^{-5}$ of the entire galactic stellar population. The disk and bulge components host only a small fraction of these relics, contributing less than $12$ percent in total. The stellar halo, in particular the outer stellar halo of which galactic radius $r>30$ kpc, hosts the largest fraction (about 46 percent on average), with an average of one relic star for per $4,000$ to $10,000$ stars, making it a promising region for observational searches. Additionally, around $18$ percent of the earliest stars relics are found in satellite galaxies, with smaller and older satellite galaxies tending to contain a higher proportion of these stars. Thus, low-mass and early-formed satellite galaxies are also ideal targets for finding such relics, although some satellite galaxies may lack them entirely. The spatial distribution and kinematics of these stars show good numerical convergence across different simulation resolutions. Our results provide valuable guidance for searches of the earliest stars relics and offer insights for interpreting findings from ongoing and future stellar archaeology surveys.

We investigate the occurrence of accretion bursts, dust accumulation, and the prospects for planetesimal formation in a gravitationally unstable magnetized protoplanetary disk with globally suppressed but episodically triggered magnetorotational instability (MRI), particularly in young intermediate-mass stars (YIMSs) but with a brief comparison to low-mass counterparts. We use numerical magnetohydrodynamics simulations in the thin-disk limit (FEOSAD code) to model the formation and long-term evolution of a gravitationally unstable magnetized protoplanetary disk, including dust dynamics and growth, since the collapse of a massive slowly-rotating prestellar cloud core. Massive gas concentrations and dust rings form within the inner disk region owing to the radially varying efficiency of mass transport by gravitational instability (GI). These rings are initially susceptible to streaming instability (SI). However, gradual warming of the dust rings, thanks to high opacity and GI-induced influx of matter increases the gas temperature above a threshold for the MRI to develop via thermal ionization of alkaline metals. The ensuing MRI bursts destroy the dust rings, making planetesimal formation via SI problematic. In the later evolution phase, when the burst activity starts to diminish, SI becomes inefficient because of growing dust drift velocity and more extended inner dead zone, both acting to reduce the dust concentration below the threshold for SI to develop. Low-mass objects appear to be less affected by these adverse effects. Our results suggest that disks around young intermediate-mass stars may be challenging environments for planetesimal formation via SI. This may explain the dearth of planets around stars with $M_\ast > 3.0 \,$$M_\odot$.

The shapes of dark matter haloes are sensitive to both cosmology and baryon physics, but are difficult to measure observationally. A promising way to constrain them is to use the positions of satellite galaxies as tracers of the underlying dark matter, but there are typically too few galaxies per halo for reliable shape estimates, resulting in biased shapes. We present a method to model sampling noise to correct for the shape bias. We compare our predicted median shape bias with that obtained from the FLAMINGO suite of simulations and find reasonable agreement. We check that our results are robust to resolution effects and baryonic feedback. We also explore the validity of our bias correction at various redshifts and we discuss how our method might be applied to observations in the future. We show that median projected halo axis ratios are on average biased low by 0.31 when they are traced by only 5 satellites. Using the satellite galaxies, the projected host halo axis ratio can be corrected with a residual bias of ~ 0.1, by accounting for sampling bias. Hence, about two-thirds of the projected axis ratio bias can be explained by sampling noise. This enables the statistical measurement of halo shapes at lower masses than previously possible. Our method will also allow improved estimates of halo shapes in cosmological simulations using fewer particles than currently required.

High-velocity clouds (HVCs) in the Galactic center have garnered significant attention due to their mysterious formation, potentially linked to starburst events or supermassive black hole activity in the region. However, it remains challenging to explain the observed column density and velocity distribution of HVCs. The discovery of high-velocity molecular clouds (HVMCs), which are denser and more massive, adds to this complexity. To address this, we conduct three-dimensional numerical simulations to explore the origin and magneto-hydrodynamic evolution of HVCs in the context of a starburst in the Galactic center. By incorporating magnetic fields and an initial tangential velocity for the clouds, our simulation results align with the observed properties of HVCs, supporting the notion that these clouds can originate from a starburst process. In addition, ~5% of the total mass of initial clouds can survive after 3.5 Myr, as a result, the following star formation will be more efficient than a feedback process that destroys all cool clouds.

As the closest known active galactic nucleus, Centaurus A (Cen A) provides a rich environment for astrophysical exploration. It has been observed across wavelengths from radio to gamma rays, and indications of ongoing particle acceleration have been found on different scales. Recent measurements of very-high-energy (VHE) gamma-rays ($>240$ GeV) by the HESS observatory have inferred the presence of ultra-relativistic electrons along Cen A's jet, yet the underlying acceleration mechanism remains uncertain. Various authors have proposed that jet substructures, known as knots, may serve as efficient particle accelerators. In this study, we investigate the hypothesis that knots are the particle acceleration sites along Cen A's jets. We focus on stationary knots, and assume that they result from interactions between the jet and the stellar winds of powerful stars. By combining relativistic hydrodynamic simulations and shock acceleration theory with the radio and X-ray data, we compare theoretical predictions with morphological and spectral data from different knots. We estimate the maximum electron energy and the resulting VHE gamma-ray emission. Our findings suggest that electrons accelerated at the knots are responsible for the gamma-ray spectrum detected in the VHE band.

Recently, Leitinger et al. (2023) and Cadelano et al (2024}) have shown that, in some Milky Way globular clusters (GCs), pristine red giant branch (RGB) stars are more centrally concentrated than enriched RGB stars. This contradicts most multiple-population formation scenarios, which predict that the enriched population 2p should initially be more concentrated than the pristine population 1P. We analyze a MOCCA GC model that exhibits a higher spatial concentration of 1P RGB stars than 2P RGB stars at 13 Gyr. The MOCCA models assume the asymptotic giant branch (AGB) pollution scenario, where 2P stars are initially more concentrated than 1P stars. Our results indicate that the observed spatial distributions of multiple populations, and possibly their kinematics, are significantly shaped by dynamical interactions. These interactions preferentially eject 2P RGB progenitors from the central regions, leading to a transient over-concentration of 1P RGB stars at late times. This effect is particularly relevant for GCs with present-day of a few $10^5 M_{\odot}$, which have retained only about 10 - 20 percent of their initial mass. Such clusters may appear dynamically young due to heating from a black hole subsystem, even if they have undergone significant mass loss and dynamical evolution. Additionally, the relatively small number of RGB stars in these clusters suggests that interpreting the spatial distributions of multiple populations solely from RGB stars may lead to biased conclusions about the overall distribution of 2P and 1P. The apparent over-concentration of the 1P relative to the 2P is likely a transient effect driven by the preferential removal of 2P RGB progenitors via strong dynamical encounters. MOCCA models of multiple stellar populations based on the AGB scenario may explain anomalous features observed in some Galactic GCs, such as NGC 3201 and NGC 6101.

The origin of the cluster of S-stars located in the Galactic Centre is tied to the supermassive black hole Sagittarius A*, but exactly how is still debated. In this paper, we investigate whether the Hills mechanism can simultaneously reproduce both the S-star cluster's properties and the observed number of hypervelocity stars. To do so, we forward-model the capture and disruption of binary stars originating from the nuclear star cluster (NSC) and the clockwise disk (CWD). We find that the ratio of evolved to main-sequence S-stars is highly sensitive to the origin of the binaries, and that neither the injection of binaries from the CWD nor from the NSC exclusively can reproduce all observations. However, when considering the injection of binaries from both locations, we are able to reproduce all the observations simultaneously, including the number of observed hypervelocity stars, the evolutionary stage of the S-stars, their luminosity function and the distribution of their semi-major axes. The implications of our findings include that ~90% of hypervelocity stars ejected over the past ~10 Myr should originate from the CWD, that the main-sequence S-stars originated in the CWD, and that the evolved S-stars originated in an old stellar population such as the NSC.

With the advent of surveys like $Euclid$ and $Vera$ $C.$ $Rubin$, astrophysicists will have access to both deep, high-resolution images, and multi-band images. However, these two conditions are not simultaneously available in any single dataset. It is therefore vital to devise image deconvolution algorithms that exploit the best of the two worlds and that can jointly analyze datasets spanning a range of resolutions and wavelengths. In this work, we introduce a novel multi-band deconvolution technique aimed at improving the resolution of ground-based astronomical images by leveraging higher-resolution space-based observations. The method capitalizes on the fortunate fact that the $Vera$ $C.$ $Rubin$ $r$-, $i$-, and $z$-bands lie within the $Euclid$ $VIS$ band. The algorithm jointly deconvolves all the data to turn the $r$-, $i$-, and $z$-band $Vera$ $C.$ $Rubin$ images to the resolution of $Euclid$ by enabling us to leverage the correlations between the different bands. We also investigate the performance of deep learning-based denoising with DRUNet to further improve the results. We illustrate the effectiveness of our method in terms of resolution and morphology recovery, flux preservation, and generalization to different noise levels. This approach extends beyond the specific $Euclid$-$Rubin$ combination, offering a versatile solution to improve the resolution of ground-based images in multiple photometric bands by jointly using any space-based images with overlapping filters.

This work proposes a saliency-based attribution framework to evaluate and compare 10 state-of-the-art explainability methods for deep learning models in astronomy, focusing on the classification of radio galaxy images. While previous work has primarily emphasized classification accuracy, we prioritize model interpretability. Qualitative assessments reveal that Score-CAM, Grad-CAM, and Grad-CAM++ consistently produce meaningful attribution maps, highlighting the brightest regions of FRI and FRII galaxies in alignment with known astrophysical features. In contrast, other methods often emphasize irrelevant or noisy areas, reducing their effectiveness.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT telescope. The survey provides deep narrow-band images for 385 galaxies hosting star forming HII regions. We identify individual HII regions and measure their main physical properties such as Halpha luminosity, equivalent diameter, and electron density with the purpose of deriving standard relations as reference for future local and high-z studies of HII regions in star forming systems in different environments. For this purpose we use a complete sample of ~ 13.000 HII regions of luminosity L(Halpha)>= 10^37 erg s^-1 to derive the main statistical properties of HII regions in unperturbed systems, identified as those galaxies with a normal HI gas content (64 objects). These are the composite Halpha luminosity function, equivalent diameter and electron density distribution, and luminosity-size relation. We also derive the main scaling relations between several parameters representative of the HII regions properties (total number, luminosity of the first ranked regions, fraction of the diffuse component, best fit parameters of the Schechter luminosity function measured for individual galaxies) and those characterising the properties of the host galaxies (stellar mass, star formation rate and specific star formation rate, stellar mass and star formation rate surface density, metallicity, molecular-to-atomic gas ratio, total gas-to-dust mass ratio). We briefly discuss the results of this analysis and their implications in the study of the star formation process in galaxy discs.

Providing a detailed picture of the Sagittarius (Sgr) stream offers important constraints on the build-up of the Galactic halo as well as its gravitational potential at large radii. While several attempts have been made to model the structure of the Sgr stream, no model has yet been able to match all the features observed for the stream. Moreover, for several of these features, observational characterisation of their properties is rather limited, particularly at large distances. The aim of this work is to investigate the kinematics of the Sgr stream outermost spur feature using blue horizontal branch (BHB) stars. Candidate BHB stars were selected by combining two approaches; one capitalising on Pan-STARRS1 3$\Pi$ griz and u photometry taken as part of UNIONS, the other using Pristine Survey CaHK and SDSS ugr photometry. Follow-up optical spectra are obtained using ESO/VLT/FORS2 to confirm their BHB nature and obtain line-of-sight (LOS) velocities. Of our 25 candidates, 20 stars can be confirmed as bona fide BHB stars. Their LOS velocities, together with the 3D positions of these stars qualitatively match well with Sgr model predictions and trace the outer apocentre of the trailing arm and its spur feature very nicely. The quantitative offsets that are found between our data and the different models can be used to provide information about the Galactic gravitational potential at large distances. We present a first, tentative, analysis in this direction, showing that the model of Vasiliev et al. (2021) would provide better agreement with our observations if the enclosed mass of the Milky Way within 100 kpc were lowered to $(5.3\!\pm\!0.4)\!\times\!10^{11}$ M$_\odot$ (versus $(5.6\!\pm\!0.4)\!\times\!10^{11}$ M$_\odot$). Our selection of BHB stars provides a new view on the outermost structure in 3D positions and LOS velocities of the Sgr debris.

Standard cosmological weak lensing analyses using cosmic shear are inevitably sensitive to small-scale, non-linear clustering from low-redshift structures. The need to adequately model the clustering of matter on this non-linear regime, accounting for both gravitational and baryonic effects, adds significant uncertainty to weak lensing studies, particularly in the context of near-future Stage-IV datasets. In this paper, inspired by previous work on so-called ``nulling'' techniques, we present a general method that selects the linear combinations of a given tomographic cosmic shear dataset that are least sensitive to small-scale non-linearities, by essentially suppressing the contribution from low-redshift structures. We apply this method to the latest public cosmic shear data from the Dark Energy Survey, DES-Y3, that corresponds to 3 years of observation, and show: a) that a large fraction of the signal is dominated by the single mode that is most affected by non-linear scales, and b) that removing this mode leads to a $\sim1\sigma$ upwards shift in the preferred value of $S_8\equiv\sigma_8\sqrt{\Omega_M/0.3}$, alleviating the tension with current CMB data. However, the removal of the most contaminated mode also results in a significant increase in the statistical uncertainties. Taking this into account, we find this shift to be compatible with a random fluctuation caused by removing this most-contaminated mode at the $\sim1.4\sigma$ level. We also show that this technique may be used by future Stage-IV surveys to mitigate the sensitivity of the final constraints to baryonic effects, trading precision for robustness.

The local stellar halo of the Milky Way is known to contain the debris from accreted dwarf galaxies and globular clusters, in the form of stellar streams and over-densities in the space of orbital properties (e.g. integrals of motion). While several over-densities have been uncovered and characterised dynamically using Gaia data, their nature is not always clear. Especially for a full understanding of the smaller halo substructures, the kinematic information from Gaia needs to be coupled with chemical information. In this work, we combine Gaia data with targeted high-resolution UVES spectroscopy of five small substructures that were recently discovered in the local halo, namely the ED-2, 3, 4, 5 and 6; the ED streams. We present the chemical abundances measured from our newly obtained UVES spectra (20 stars) and from archival UVES spectra (9 stars). We compare these with homogeneously derived abundances from archive spectra of 12 Gaia Enceladus (GE) stars. The chemical abundances of all 5 substructures suggest that they are of accreted origin, except for two stars that present a high [{\alpha}/Fe] at high [Fe/H] more in line with an in-situ origin. All but ED-2 present a significant spread in [Fe/H] suggestive of a dwarf galaxy origin. ED-3 and ED-4 tend to exhibit lower [{\alpha}/Fe] compared to GE stars. ED-5 and ED-6 are consistent with the GE chemical track and could be high-energy tails of GE that were lost earlier in the accretion process. We present new elemental abundances for 5 ED-2 stars, including more elements for the Gaia BH3 companion star. Our findings are in line with the picture that ED-2 is a disrupted ancient star cluster.

Previous studies suggested that the back-reaction of super-Hubble cosmological fluctuations could lead to a dynamical relaxation of the cosmological constant. Moreover, this mechanism appears to be self-regulatory, potentially leading to an oscillatory behavior in the effective dark energy. Such an effect would occur in any cosmological model with super-Hubble matter fluctuations, including the standard $\Lambda\text{CDM}$ model, without the need for additional mechanisms. The recent DESI data, which supports the possibility of dynamical dark energy, has renewed interest in exploring scenarios leading to such oscillatory behavior. In this work we propose a parametrization to account for the impact of super-Hubble fluctuations on the background energy density of the Universe. We model the total effective cosmological constant as the sum of a constant piece and an oscillating contribution. We perform a preliminary comparison of the background dynamics of this model with the recent radial BAO data from DESI. We also discuss the status of the $H_0$ tension problem in this model. While this analysis serves as a first step in testing the scenario, we emphasize the importance of considering this effect, as it is naturally expected to be present and could offer valuable insights in light of the recent observational data.

GJ 436b is the archetype warm Neptune exoplanet. The planet's thermal emission spectrum was previously observed via intensive secondary eclipse campaigns with Spitzer. The atmosphere has long been interpreted to be extremely metal-rich, out of chemical equilibrium, and potentially tidally heated. We present the first panchromatic emission spectrum of GJ 436b observed with JWST's NIRCAM (F322W2 and F444W) and MIRI (LRS) instruments between 2.4 and 11.9 $\mu$m. Surprisingly, the JWST spectrum appears significantly fainter around 3.6 $\mu$m than that implied by Spitzer photometry. The molecular absorption features in the spectrum are relatively weak, and we only find tentative evidence of CO$_2$ absorption at 2$\sigma$ significance. Under the assumption of a day-side blackbody, we find $T_{\rm day}$=662.8$\pm$5.0 K, which is similar to the zero Bond albedo equilibrium temperature. We use it to obtain a 3$\sigma$ upper limit on the Bond albedo of $A_B{\le}$0.66. To understand the spectrum we employ 1D radiative-convective models but find that atmospheric constraints depend strongly on model assumptions. If thermochemical equilibrium is assumed, we find a cloudy metal-enriched atmosphere (metallicity $\ge$ 300$\times$solar). We employ 1D photochemical modeling to show that the observed spectrum is also consistent with a cloud-free, relatively lower-metallicity atmosphere (metallicity $\ge$ 80$\times$solar) with a cold internal temperature ($T_{\rm int}$$\sim$60 K). These are much lower metallicities and internal temperatures than inferences from Spitzer photometry. The low $T_{\rm day}$ and non-detection of transmission features at high spectral resolution does suggest a role for cloud opacity, but this is not definitive.

Henrietta Swan Leavitt's discovery of the relationship between the period and luminosity (hereafter the Leavitt Law) of 25 variable stars in the Small Magellanic Cloud, published in 1912, revolutionized cosmology. These variables, eventually identified as Cepheids, became the first known "standard candles" for measuring extragalactic distances and remain the gold standard for this task today. Leavitt measured light curves, periods, and minimum and maximum magnitudes from painstaking visual inspection of photographic plates. Her work paved the way for the first precise series of distance measurements that helped set the scale of the Universe, and later the discovery of its expansion by Edwin Hubble in 1929. Here, we re-analyze Leavitt's first Period-Luminosity relation using observations of the same set of stars but with modern data and methods of Cepheid analysis. Using only data from Leavitt's notebooks, we assess the quality of her light curves, measured periods, and the slope and scatter of her Period-Luminosity relations. We show that modern data and methods, for the same objects, reduce the scatter of the Period-Luminosity relation by a factor of two. We also find a bias brightward at the short period end, due to the non-linearity of the plates and environmental crowding. Overall, Leavitt's results are in excellent agreement with contemporary measurements, reinforcing the value of Cepheids in cosmology today, a testament to the enduring quality of her work.

This work studies interference patterns created by simple lens models (point mass, Chang-Refsdal, and binary lens) in the wave optics regime, primarily in the context of lensing of gravitational waves (GWs) in the LIGO band at frequencies around 100 Hz. We study how the interference patterns behave close to the caustic curves which mark the high magnification regions in conventional geometric optics. In addition, we also look at the formation of highly de-amplified regions in the amplification maps close to caustics and how they differ under wave and geometric optics. We see that for a source close to caustics, the oscillations in the amplification factor (their amplitude and location of crests and troughs) can differ significantly in wave optics compared to geometric optics. As we move away from caustics, the wave optics amplification factor starts to converge towards geometric optics one, especially the frequencies at which crests and through occur in the amplification factor, although the amplitude of these oscillations can still be considerably different. For Chang-Refsdal and binary lens with ${\sim}100\:{\rm M_\odot}-200\:{\rm M_\odot}$ can introduce significant de-amplification at frequencies ${\sim}100$ Hz when the source is close to caustics which may help us distinguish such lenses from the point mass lens.

TianQin is a dedicated geocentric mission for space-based gravitational wave (GW) detection. Among its core technologies, the drag-free and pointing control subsystem (DFPCS) - consisting of suspension, drag-free and pointing controls - keeps the two test masses (TMs) centered and aligned within their housings while maintaining drag-free conditions and precise telescope pointing along the laser-arm directions. This results in orbit-attitude coupled dynamics for the constellation. The coupling is made more prominent due to satellite self-gravity, which requires compensation from DFPCS and generally makes the satellites deviate from pure free-fall orbits. Previous studies assumed that the orbit and attitude dynamics could be decoupled in numerical simulation, neglecting the back-action from the closed-loop control to orbit propagation. To address this, we develop a comprehensive model that can propagate the full 9-body (6 TMs + 3 satellites, orbits and attitudes) dynamics inter-dependently under the inter-satellite pointing and drag-free conditions. This paper is threefold. First, we reassess the applicability of the two TMs and telescope pointing scheme to TianQin using the new model, and confirm the previous conclusion. Second, to meet the constellation stability requirements, it is found that the DC common self-gravity in the flight direction should be minimized, or kept close for the three satellites. Finally, we simulate the long-range light path between two TMs/satellites with a precision of sub-pm/Hz$^{1/2}$, and the results support the decoupling of the closed-loop dynamics and high-precision orbit for computational efficiency. The method is instrumental to other future missions where the orbit-attitude coupling needs careful consideration.

We introduce a benchmark to evaluate the capability of AI to solve problems in theoretical physics, focusing on high-energy theory and cosmology. The first iteration of our benchmark consists of 57 problems of varying difficulty, from undergraduate to research level. These problems are novel in the sense that they do not come from public problem collections. We evaluate our data set on various open and closed language models, including o3-mini, o1, DeepSeek-R1, GPT-4o and versions of Llama and Qwen. While we find impressive progress in model performance with the most recent models, our research-level difficulty problems are mostly unsolved. We address challenges of auto-verifiability and grading, and discuss common failure modes. While currently state-of-the art models are still of limited use for researchers, our results show that AI assisted theoretical physics research may become possible in the near future. We discuss the main obstacles towards this goal and possible strategies to overcome them. The public problems and solutions, results for various models, and updates to the data set and score distribution, are available on the website of the dataset this http URL.

Fluctuations in ultralight dark matter produce significant metric perturbations, which may be detected by monitoring the arrival times of light from millisecond pulsars. While searches using this technique are already underway, they do not consistently account for the statistical properties of the dark matter field. The statistics of this field depend on the velocity dispersion of dark matter and, consequently, its coherence length. In the mass range relevant for pulsar timing arrays, the coherence length is comparable to separations between pulsars, making it crucial to incorporate its effects into the analysis. This work presents a consistent statistical method for gravitational direct detection of ultralight dark matter. Our key result is the derivation of the two-point function of the metric fluctuations, which we apply to pulsar timing and discuss its implementation in future searches.

In this paper, we extend the functional approach for calculating the EFT likelihood by applying the saddle-point expansion. We demonstrate that, after suitable reformulation, the likelihood expression is consistent with the path integral required to be computed in the theory of false vacuum decay. In contrast to the saddle-point approximation, the application of the saddle-point expansion necessitates more nuanced considerations, particularly concerning the treatment of the negative eigenvalues of the second derivative of the action at the saddle point. We illustrate that a similar issue arises in the likelihood calculation, which requires approximating the original integral contour through the combination of the steepest descent contours in the field space. Gaussian distribution as a working example, we then concentrate on the calculation of the EFT likelihood and propose a general procedure for computing the likelihood via saddle-point expansion method for arbitrary partition functions. Precise computation of the likelihood will benefit Bayesian forward modeling, thereby enabling more reliable theoretical predictions.

The Hamilton-Jacobi method offers a natural and concise framework for describing inflation, with implications that extend to the reheating phase. Additionally, reheating plays a crucial role in constraining the observationally viable parameter space of inflationary models. In this study, we employ the Hamilton-Jacobi approach to investigate reheating predictions within non-minimally coupled inflation models, comparing the metric and Palatini formulations. Our results show that the coupling effect suppresses the tensor-to-scalar ratio, aligning predictions with the Planck CMB and BICEP/Keck data in both formulations. Additionally, reheating predictions in the Palatini formulation are more sensitive to coupling strength variations, leading to a stronger suppression of the tensor-to-scalar ratio. This highlights a key difference in reheating dynamics between the metric and Palatini formulations.

This work studies the influence of scalar dark matter on the structural properties of strange quark stars (SQS) within a one-fluid framework, considering Yukawa interactions between dark matter and quark matter. Contributions from perturbative QCD, Yukawa interaction between scalar dark matter and quarks, and Bose-Einstein condensation of dark matter are included in the model. We first determine the allowable range of Yukawa interaction coupling by imposing the stability condition for strange quark matter (SQM). Using this range, we derive the equation of state (EOS) for different fractions of dark matter within the total pressure of SQS. These fractions are constrained by the tidal deformability limit from GW170817. The presence of dark matter alters the EOS, leading to changes in the mass-radius relationship, tidal deformability, and stability of SQS. We demonstrate that increasing the mass of dark matter softens the EOS, whereas higher fractions of dark matter lead to stiffer EOSs. We also explore the reasons behind this behavior. Our EOSs not only describe massive objects, such as PSR J0952-0607 and PSR J2215+5135, but also satisfy the tidal deformability constraint from GW170817. These results reveal that incorporating dark matter modifies the EOS, enabling the support of higher stellar masses while maintaining consistency with observational data.

Assuming that (1) the universe underwent a post-inflationary accelerated expansion phase driven by a fluid with equation of state $P=w\rho$ and $-1<w<-1/3$, that (2) the cosmic horizon in an accelerating, quasi-de Sitter universe has a temperature inversely proportional to the proper size of the horizon, and that (3) we are static observers, we calculate the frozen-in density of a stable particle of mass $m$ produced by the cosmic horizon that does not undergo any number-changing processes in the late universe. We find that, as a function of the equation of state and the temperature when radiation domination starts and the quasi-de Sitter phase ends, the mass of the dark matter producing the observed cosmological abundance via this mechanism ranges from 10 keV up to close to the Planck scale.

We explore a novel mechanism for dark matter production through the formation of light black holes from the collapse of dark baryons in confining SU(N) gauge theories in the large $N$ limit. While glueballs and mesons cannot form black holes under physically reasonable conditions, we prove that for appropriate ranges of the confinement scale, quark masses, number of colors $N$, and dark sector temperature, dark baryons can produce Planck-scale black hole relics in the early universe. Assuming the relics are stable, the abundance of both the dark baryon black hole population directly arising at confinement and that frozen in from dark glueball and meson pair annihilation are exponentially suppressed in $N$, leading to an upper limit $N\lesssim 100$ and of a few hundreds Planck units in mass for models where the black hole relics are the entirety of the dark matter. We present a detailed numerical study of the parameter space where this scenario is realized.

We investigate the influence of an early matter-dominated era in cosmic history on the dynamics of cosmic strings and the resulting stochastic gravitational waves. Specifically, we examine the case where this era originates from the dark matter dilution mechanism within the framework of the minimal left-right symmetric model. By numerically solving the Boltzmann equations governing the energy densities of the relevant components, we meticulously analyze the modifications to the cosmological scale factor, the number density of cosmic string loops, and the gravitational wave spectrum. Our results reveal that the early matter-dominated era causes a characteristic suppression in the high-frequency regime of the gravitational wave spectrum, providing distinct and testable signatures for future ground-based interferometer experiments.

In a recent work [Phys. Rev. D 110, 023520 (2024)] a baryogenesis scenario was proposed where our visible Universe is a 3-brane coevolving with a hidden 3-brane in a multidimensional bulk. This model introduced a new pseudo-scalar boson. In the present paper, it is shown that this boson can exist today as a relic from the Big Bang and constitute a minor dark matter component. Additionally, one identifies a one-loop process that allows for a non-zero, albeit very small, photon-boson scattering amplitude. One explores this phenomenon and suggests that observations of a subtle light scattering around blue giant stars could constrain the properties of the scalar boson. These constraints would, in turn, provide valuable insights into the parameters of the proposed baryogenesis model.

Some astrobiological models suggest that molecular clouds may serve as habitats for extraterrestrial life. This study reviews recent theoretical work addressing the physical and biochemical prerequisites for life in such environments, with particular focus on three subjects: (1) bioenergetic pathways under extreme low-temperature conditions; (2) the emergence and preservation of biomolecular chirality; and (3) detection methodologies for potential biosignatures. In this paper, we formally introduce the molecular cloud biology concept, which integrates all physicochemical and metabolic processes hypothesized to sustain life within molecular clouds. As a potential branch of astrobiology, molecular cloud biology warrants interdisciplinary collaborative research to validate its foundational assumptions and explore its scientific implications.

In this paper, we conduct an in-depth investigation into the optical images of boson stars with the solitonic potential. In the context of a celestial source and a thin accretion disk, the optical characteristics of the soliton boson star have been derived. Considering the influence of the initial scalar field $\psi_0$ and a larger coupling parameter $\alpha$ (the weak coupling case), the optical images of boson stars primarily exhibit direct and lensed images. The results demonstrate that variations in $\psi_0$ and $\alpha$ influence the image size, whereas the observer's inclination angle $\theta$ has a substantial impact on the image shape. In contrast, when the coupling parameter $\alpha$ is small (the strong coupling case), a sub-annular structure emerges within the Einstein ring for a spherically symmetric light source. In the presence of a thin accretion disk, higher-order gravitational lensing images emerge, indicating that photons are capable of orbiting the equatorial plane of the boson star multiple times. We also analyze how the effective potential and redshift factor depend on the parameters $\psi_0$, $\alpha$, and $\theta$. The results indicate that at smaller values of $\theta$, gravitational redshift is the dominant effect, resulting in an optical image featuring a bright ring surrounding a comparatively dim central region. At larger values of $\theta$, the Doppler effect becomes more pronounced, resulting in a substantial brightness disparity between the left and right sides of the optical image. These findings offer robust theoretical underpinnings for differentiating solitonic boson stars from black holes via high-resolution astronomical observations.

With the detection of gravitational waves (GW) in recent years, many papers that use simulated GW datasets to constrain dark energy models have been published. However, these works generally consider GW generated by massive astrophysical objects, such as mergers of black holes or neutron stars. In this paper, we analyse the evolution and spectrum of GW generated in the inflationary epoch, assuming a standard slow-roll single-field inflationary scenario. Three models for the background are considered: a field theory model of interacting dark energy - the Interacting Holographic Tachyonic Model, the Holographic Dark Energy Model, and the {\Lambda}CDM. The results show significant dependence on the cosmological model, especially for the spectrum, and in particular show that an interaction between dark energy and dark matter can leave a significant imprint. Therefore, future primordial gravitational waves (PGW) datasets could be very useful for constraining dark energy models, including to probe an interaction in the dark sector of the universe.

Experiments at the US National Ignition Facility (NIF) [Döppner et al., Nature {\bf 618}, 270-275 (2023)] have created highly compressed hot hydrogen-like Be plasmas. Published analyses of the the NIF experiment have used finite-$T$ multi-atom density-functional theory (DFT) with Molecular dynamics (MD), and Path-Integral Monte Carlo (PIMC) simulations. These methods are very expensive to implement and often lack physical transparency. Here we (i) relate their results to simpler first-principles average-atom results, (ii) establish the feasibility of rapid data analysis, with good accuracy and gain in physical transparency, and (iii) show that the NIF experiment reveals high-$T$ spin-singlet pairing of hydrogen-like Be ions with near neighbours. Our analysis predicts such stabilization over a wide range of compressed densities for temperatures close to two million Kelvin. Calculations of structure factors $S(k)$ for electrons or ions, the Raleigh weight and other quantities of interest to X-ray Thomson scattering are presented. We find that the NIF data at the scattering wavevector $k_{sc}$ of 7.89 Å$^{-1}$ are more consistent with a density of $20\pm2$ g/cm$^3$, mean ionization $\bar{Z}=$3.25, at a temperature of $\simeq$ 1,800,000 K than the 34 g/cm$^3, \bar{Z}=3.4$ proposed by the NIF team. The relevance of ion-electron coupled-modes in studying small $k_{sc}$ data is indicated.

We summarize our recent achievements in the description of neutrino oscillations in various gravitational fields. After the short review of the previous studies of neutrinos in gravitational fields, we consider the neutrinos propagation and oscillations in two gravitational backgrounds. First, we discuss neutrino spin oscillations in their gravitational scattering off a supermassive black hole surrounded by a thick magnetized accretion disk. Second, we study neutrino flavor oscillations in stochastic gravitational waves. We also consider applications of the obtained results for oscillations of astrophysical neutrinos.

In a cosmological setting, particle production is ubiquitous. It may occur as a consequence of the expansion of the background or because a field couples to other degrees of freedom that evolve with time. The process is well understood in the context of quantum field theory, and calculable as long as the produced quanta are weakly interacting. For extended objects like strings and membranes a second quantized formulation is much less developed than for particles. In this light, we revisit particle production from a first quantized perspective. We show how to obtain occupation numbers, both from the vacuum persistence amplitude and from the Green's function. We also derive the much less studied but phenomenologically interesting two-particle wavefunction of the produced quanta.

We adopt an effective action inspired by asymptotically safe gravity, in which the effective gravitational constant is parameterized as $G(\epsilon) = G_{N} [1 + \tilde{\omega} (G_{N}^{2} \epsilon)^{\alpha}]^{-1}$, where $G_{N}$ and $\epsilon$ denote Newton's gravitational constant and the energy density of the matter field, respectively, with two dimensionless model parameters, $\tilde{\omega}$ and $\alpha$. Within this framework, we investigate the complete gravitational collapse of a homogeneous ball of perfect fluid and find that the singularity is completely resolved for $\alpha > 1$ but not for $1/2 \le \alpha \le 1$. The case $0 < \alpha < 1/2$ is inconsistent with asymptotic safety. Moreover, we note that although the singularity cannot be fully resolved for $\alpha = 1$, it is significantly weakened by quantum gravity effects. Furthermore, we successfully construct a static exterior metric which, together with the interior solution, describes the dynamical formation of regular black holes in an asymptotically flat spacetime for the perfectly resolved case $\alpha > 1$. The resulting regular black hole, obtained as the final static state, contains a de Sitter core and admits a static metric fully expressible in terms of the Lerch transcendent for general cases and in elementary functions for certain values of $\alpha$, including \alpha = 2$. We also discuss the formation of gravastars and the late-time evaporation process of the regular black holes.

Using direct numerical simulations of forced rotating turbulence, we study the effect of rotation on the growth rate and the saturation level of the small-scale dynamo. For slow rotation rates, increasing the rotation rate reduces both the growth rate and the saturation level. Once the rotation rate crosses a threshold, large-scale vortices are formed which enhance the growth rate and the saturation level. Below this threshold, the suppression of the small-scale dynamo with increasing rotation is explained by the fact that at scales close to, but smaller than, the forcing scale, rotating turbulence is one-dimensionalized, with the velocity component along the rotation axis being larger than the other two components. This is due to the rotational destabilization of vortices produced by the forcing function. While the rotational effect on the growth rate becomes small at high Re, the ratio of the steady-state magnetic to kinetic energies remains suppressed by up to 35% as compared to the non-rotating case.

How much energy is required to unbind baryons from the cosmological structures that originally bind them? This tutorial article explains why trying to answer this question using just a halo model can be misleading. Instead, it recommends parsing the universe into ``bound domains,'' which are the gravitationally bound structures that ultimately become widely separated islands as the universe evolves. It explains why a bound domain's potential well was about as deep ~1 Gyr after the Big Bang as it is now, and it outlines how future research might take advantage of a bound-domain approach to make progress on some open questions about the baryon distributions in and around galaxy groups and clusters.

The calculation of gray-body factors is essential for understanding Hawking radiation and black hole thermodynamics. While the formalism developed by Chandrasekhar is effective for static black holes, it faces significant challenges in Kerr spacetimes, particularly in the superradiant regime, where a specific choice of coordinates introduces numerical inaccuracies. To address these limitations, an alternative method based on re-scaling radial coordinates and employing Frobenius-like expansions has been investigated. We compare the gray-body factors obtained for a near-maximally rotating black hole using both methods and find that the Chandrasekhar formalism systematically overestimates the values in the superradiant regime compared to well-established analytical results. Specifically, for a spin parameter of $a_* = 0.999$, the Chandrasekhar method yields values approximately twice as large as the correct result. Since this approach has been implemented in \texttt{BlackHawk}, we assess the impact of these discrepancies on constraints derived from gamma-ray observations of highly spinning primordial black holes.

We examine the sterile neutrino dark matter production in the primordial plasma with lepton asymmetry unequally distributed over different neutrino flavors. We argue that with the specific flavor fractions, one can mitigate limits from the Big Bang Nucleosynthesis on the sterile-active neutrino mixing angle and sterile neutrino mass. It happens due to cancellation of the neutrino flavor asymmetries in active neutrino oscillations, which is more efficient in the case of inverse hierarchy of active neutrino masses and does not depend on the value of CP-phase. This finding opens a window of lower sterile-active mixing angles. Likewise, we show that, with lepton asymmetry disappearing from the plasma at certain intermediate stages of the sterile neutrino production, the spectrum of produced neutrinos becomes much colder, which weakens the limits on the model parameter space from observations of cosmic small-scale structures (Ly-$\alpha$ forest, galaxy counts, etc.). This finding reopens the region of lighter sterile neutrinos. The new region may be explored with the next generation of X-ray telescopes searching for the inherent peak signature provided by the dark matter sterile neutrino radiative decays in the Galaxy.

The Parameter Estimation (PE) for Gravitational Waves (GW) merger events relies on a waveform model calibrated using numerical simulations. Within the Bayesian framework, this waveform model represents the GW signal produced during the merger and is crucial for estimating the likelihood function. However, these waveform models may possess systematic errors that can differ across the parameter space. Addressing these errors in the current data analysis pipeline is an active area of research. This work presents a framework for accounting for uncertainties in waveform modeling. We introduce two parametrizations, relative and absolute errors in the phase of the waveform, to modify the base waveform model, which can account for uncertainties. When the waveform errors are known, those error budgets can be used as a prior distribution in the Bayesian framework. We also show that conservative priors can be used to quantify uncertainties in waveform modeling without any knowledge of waveform error budgets. By conducting zero-noise injections and recoveries, we demonstrate through PE results that even 1-2% of errors in relative phase to the actual waveform model can introduce biases in the recovered parameters. These biases can be corrected when we account for waveform uncertainties within the PE framework. By injecting a series of precessing waveform models and using the nonspinning model for recovery, we show that our method can account for the missing physics by making the posterior samples broad enough to account for bias. We also present a Python package that is easily integrated with the publicly available GW analysis tool PyCBC and can be used to do PE with the parametrization presented in this paper.

The astrophysical gravitational wave background in the nanohertz (nHz) band is expected to be primarily composed of the superposition of signals from binaries of supermassive black holes. The spatial discreteness of these sources introduces shot noise, which, in certain regimes, would overwhelm efforts to measure the anisotropy of the gravitational wave background. In this work, we explicitly demonstrate, starting from first principles, that cross-correlating a gravitational wave background map with a sufficiently dense galaxy survey can mitigate this issue. This approach could potentially reveal underlying properties of the gravitational wave background that would otherwise remain obscured. We quantify the shot noise level and show that cross-correlating the gravitational wave background with a galaxy catalog improves by more than two orders of magnitude the prospects for a first detection of the background anisotropy by a gravitational wave observatory operating in the nHz frequency range, provided it has sufficient sensitivity.

The spins of black holes in binaries measured with gravitational waves provide insights about the formation, evolution, and dynamics of these systems. The imprint of spin in the inspiral, where the black holes are well-separated, is understood through analytic equations for the binary dynamics. During the merger phase, the binary dynamics can only be studied with numerical relativity simulations. Though such simulations provide an exact solution (to within numerical error), the imprint of the full six spin degrees of freedom on the signal is not transparent. In the absence of analytic expressions for the merger, here we propose a waveform-based approach. Leveraging a neural network to efficiently calculate mismatches between waveforms, we identify regions in the parameter space of spins and mass ratio that result in low mismatches and thus similar waveforms. We map these regions with a Gaussian fit, thus identifying correlations between the mass ratio and spins and quantifying their strength. For low-mass, inspiral-dominated systems, we recover the known physical imprint: larger aligned spins are correlated with more equal masses as they have opposite effects on the inspiral length. For high-mass, merger-dominated signals, a qualitatively similar correlation is present, though its shape is altered and strength decreases with increasing total mass. Correlations between in-plane spins and mass ratio follow a similar trend, with their shape and strength altered as the mass increases. Waveform-based correlation mapping can motivate effective spin parameters and reveal the imprint of spins on signals for which no simple analytic descriptions exist.

Due to the sheer complexity of the Laser Interferometer Space Antenna (LISA) space mission, data gaps arising from instrumental irregularities and/or scheduled maintenance are unavoidable. Focusing on merger-dominated massive black hole binary signals, we test the appropriateness of the Whittle-likelihood on gapped data in a variety of cases. From first principles, we derive the likelihood valid for gapped data in both the time and frequency domains. Cheap-to-evaluate proxies to p-p plots are derived based on a Fisher-based formalism, and verified through Bayesian techniques. Our tools allow to predict the altered variance in the parameter estimates that arises from noise mismodeling, as well as the information loss represented by the broadening of the posteriors. The result of noise mismodeling with gaps is sensitive to the characteristics of the noise model, with strong low-frequency (red) noise and strong high-frequency (blue) noise giving statistically significant fluctuations in recovered parameters. We demonstrate that the introduction of a tapering window reduces statistical inconsistency errors, at the cost of less precise parameter estimates. We also show that the assumption of independence between inter-gap segments appears to be a fair approximation even if the data set is inherently coherent. However, if one instead assumes fictitious correlations in the data stream, when the data segments are actually independent, then the resultant parameter recoveries could be inconsistent with the true parameters. The theoretical and numerical practices that are presented in this work could readily be incorporated into global-fit pipelines operating on gapped data.