Papers for Friday, Dec 06 2024

The Roman Space Telescope, equipped with a 2.4 meter primary mirror and optical--NIR wide-field camera, promises to revolutionize our understanding of dark energy, exoplanets, and infrared astrophysics. One of the Roman Core Community Surveys is the High Latitude Time Domain Survey (HLTDS), which will measure more than 10,000 SN Ia light curves but obtain a fraction of this number with spectra. The remaining SNe will have to be photometrically classified to achieve the full potential of the Roman HLTDS. To investigate transient yields and classifications, Rose et al. (in prep.) updated the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) framework (originally developed for the Vera Rubin Observatory) for the Roman HLTDS. This study leverages this Roman Hourglass dataset to train and evaluate the ParSNIP (Parameterized Supernova Identification Pipeline) model. We employ this model to classify various transient types from photometric data, paying particular attention to the types most represented in the dataset: normal SNe Ia, 91bg-like SNe Ia, SNe Iax, and CC SNe. The ParSNIP model performance is assessed through confusion matrices and ROC curves across different redshift ranges. Our analysis reveals that while the model performs robustly at higher redshifts (with the AUC for classification varying between 0.9 and 0.95 in the range 0.5 < z < 2), its accuracy dips at the lowest redshifts of the survey, likely due to limited training data. These findings underscore the importance of ensuring adequate representation of classes in the training set. This work underscores the value of machine learning models for next-generation surveys, paving the way for future studies with the Roman Space Telescope for survey optimization, cosmological forecasts, and synergies with other surveys.

In this work, we leverage CMB data from the Atacama Cosmology Telescope (ACT) and LSS data from the imaging survey conducted by the Dark Energy Spectroscopic Instrument (DESI) to study the distribution of gas around galaxy groups at low redshift, $z \approx 0.3$, via the kinematic Sunyaev-Zel'dovich (kSZ) effect. In particular, we perform velocity-weighted stacking on the photometric Bright Galaxy Sample (BGS) to isolate the monopole and quadrupole of the kSZ signal, orienting the stacked images along 2D filaments identified using the Hessian of the projected gravitational potential. We find a 7.2$\sigma$ detection in the monopole of the signal (i.e., the gas density profile) and a 4$\sigma$ detection in the quadrupole ($m = 2$), constituting the first measurement of the alignment between gas distribution and the cosmic web through the kSZ effect. As it is a linear probe of the local gas density, the kSZ has heightened sensitivity to the warm-hot intergalactic medium (WHIM), which is believed to house the majority of the ``missing baryons.'' Mapping out the gas density at low redshifts, as enabled by our measurements, is crucial for weak lensing surveys, for which the impact of baryons on small scales is a major impediment. We compare the anisotropic signal against two hydrodynamical simulations, TNG300-1 and Illustris, which have very different baryonic feedback prescriptions. We find that the anisotropic signal measured in the data is comparable but slightly larger and more extended compared with the simulations. This suggests that there is excess accretion and feedback taking place through the filaments, hinting at the possible presence of spin-filament alignment of the BGS objects.

Stellar disk truncations, also referred to as galaxy edges, are key indicators of galactic size, determined by the radial location of the gas density threshold for star formation. Accurately measuring galaxy sizes for millions of galaxies is essential for understanding the physical processes driving galaxy evolution over cosmic time. In this study, we aim to explore the potential of the Segment Anything Model (SAM), a foundation model designed for image segmentation, to automatically identify disk truncations in galaxy images. With the Euclid Wide Survey poised to deliver vast datasets, our goal is to assess SAM's capability to measure galaxy sizes in a fully automated manner. SAM was applied to a labeled dataset of 1,047 disk-like galaxies with $M_* > 10^{10} M_{\odot}$ at redshifts up to $z \sim 1$, sourced from the HST CANDELS fields. We 'euclidized' the HST galaxy images by creating composite RGB images, using the F160W (H-band), F125W (J-band), and F814W + F606W (I-band + V-band) HST filters, respectively. Using these processed images as input for SAM, we retrieved various truncation masks for each galaxy image under different configurations of the input data. We find excellent agreement between the galaxy sizes identified by SAM and those measured manually (i.e., by using the radial positions of the stellar disk edges in galaxy light profiles), with an average deviation of approximately $3\%$. This error reduces to about $1\%$ when excluding problematic cases. Our results highlight the strong potential of SAM for detecting disk truncations and measuring galaxy sizes across large datasets in an automated way. SAM performs well without requiring extensive image preprocessing, labeled training datasets for truncations (used only for validation), fine-tuning, or additional domain-specific adaptations such as transfer learning.

The nuclear structure of dusty star-forming galaxies is largely unexplored but harbours critical information about their structural evolution. Here, we present long-baseline Atacama Large (sub-)Millimetre Array (ALMA) continuum observations of a gravitationally lensed dusty star-forming galaxy at $z=2.78$. We use a pixellated lens modelling analysis to reconstruct the rest-frame 230 $\rm\mu$m dust emission with a mean resolution of $\approx55$ pc and demonstrate that the inferred source properties are robust to changes in lens modelling methodology. The central 1 kpc is characterised by an exponential profile, a dual spiral arm morphology and an apparent super-Eddington compact central starburst. We find tentative evidence for a nuclear bar in the central 300 pc. These features may indicate that secular dynamical processes play a role in accumulating a high concentration of cold gas that fuels the rapid formation of a compact stellar spheroid and black hole accretion. We propose that the high spatial resolution provided by long-baseline ALMA observations and strong gravitational lensing will give key insights into the formation mechanisms of massive galaxies.

Galaxy-cluster virial (structure-formation accretion) shock observations are shown to imply $\gtrsim1\%$ magnetization of a layer extending $\gtrsim10^{16}$ Debye lengths downstream, challenging the modelling of high Alfvén-Mach collisionless shocks. Unlike similar shocks in supernova remnants or relativistic shocks in $\gamma$-ray burst afterglows, where macroscopic magnetized layers were detected but purportedly attributed to preexisting or non-resonant cosmic-ray streaming-seeded substructure, the upstream of strong virial shocks is both weakly magnetized and pristine. Hence, some mechanism must generate large-scale and possibly self-similar magnetic sub-structure out of the accreted primordial plasma; such a mechanism may dominate other high-Mach shock systems, too.

We present high angular resolution imaging from the {\sl Hubble Space Telescope} combined with adaptive optics imaging results from the {\sl Keck}-II telescope to determine the mass of the OGLE-2012-BLG-0563L host star and planet to be $M_{\rm host} = 0.801\pm 0.033M_\odot$ and $M_{\rm planet} = 1.116 \pm 0.087 M_{\rm Jupiter}$, respectively, located at a distance of $D_L = 5.46\pm 0.56\,$kpc. There is a close-wide degeneracy in the light curve models that indicates star-planet projected separation of $1.50\pm 0.16\,$AU for the close model and $8.41\pm 0.87\,$AU for the wide model. We used the image-constrained modeling method to analyze the light curve data with constraints from this high angular resolution image analysis. This revealed systematic errors in some of the ground-based light curve photometry that led to an estimate of the angular Einstein Radius, $\theta_E$, that was too large by a factor of $\sim 2$. The host star mass is a factor of 2.4 larger than the value presented in the \citet{fukui15} discovery paper. Although most systematic photometry errors seen in ground-based microlensing light curve photometry will not be repeated in data from the {\sl Roman Space Telescope}'s Galactic Bulge Time Domain Survey, we argue that image constrained modeling will be a valuable method to identify possible systematic errors in {\sl Roman} photometry.

The James Webb Space Telescope (JWST) observations enable the exploration of active galactic nuclei (AGNs) with broad-line emission in the early universe. Despite their clear radiative and morphological signatures of AGNs in rest-frame optical bands, complementary evidence of AGN activity -- such as X-ray emission and UV/optical variability -- remains rarely detected. The weakness of X-rays and variability in these broad-line emitters challenges the conventional AGN paradigm, indicating that the accretion processes or environments around the central black holes (BHs) differ from those of low-redshift counterparts. In this work, we study the radiation spectra of super-Eddington accretion disks enveloped by high-density coronae. Radiation-driven outflows from the disk transport mass to the poles, resulting in moderately optically-thick, warm coronae formed through effective inverse Comptonization. This mechanism leads to softer X-ray spectra and larger bolometric correction factors for X-rays compared to typical AGNs, while being consistent with those of JWST AGNs and low-redshift super-Eddington accreting AGNs. In this scenario, UV/optical variability is suppressed due to photon trapping within super-Eddington disks, while X-ray emissions remain weak yet exhibit significant relative variability. These characteristics are particularly evident in high-redshift AGNs powered by lower-mass BHs with $\lesssim 10^{7-8}~M_\odot$, which undergo rapid mass accretion following overmassive evolutionary tracks relative to the BH-to-stellar mass correlation in the local universe.

The kinematics of the ultra-diffuse galaxy (UDG) NGC1052-DF44 is primarily influenced by the presence of dark matter (DM). In this paper, we conduct a contrasting kinematic study of DF44 within the alternative modified gravity framework. In comparison to NFW DM, we test three alternative gravity models viz Milgromian dynamics (MOND), characterized by a known acceleration scale, a generic $f(R)$ model, assuming an expansion of the Ricci scalar, and a quantum gravity-inspired Renormalization Group correction to General Relativity (RGGR), which involves the running of the gravitational coupling parameter $G$ with the Universe's energy scale. For each gravity model, we evaluate the velocity dispersion (VD) of the galaxy beyond the conventional radial isotropic assumption and extend to two anisotropy scenarios, i.e., constant and Osipkov-Merritt. Our results show that all three gravity models can provide consistent fits to the observed VD of DF44; however, only MOND and RGGR remain competitive with NFW DM. Interestingly, the constant anisotropy scenario in all the models is also found to be competitive with the complete isotropic assumption.

Galaxies are predicted to assemble their stellar haloes through the accretion of stellar material from interactions with their cosmic environment. Observations that trace stellar halo buildup probe the processes that drive galaxy size and stellar mass growth. We investigate stellar halo assembly over $0.2 \leq z \leq 1.1$ in a mass-complete ($M_{\star} \geq 10^{9.5}M_{\odot}$) sample of 242,456 star-forming and 88,421 quiescent galaxies (SFGs and QGs) from the CLAUDS and HSC-SSP surveys. We extract galaxy rest-frame $g$-band surface brightness ($\mu_g$) profiles to study faint, extended emission in galaxy outskirts. We examine trends in galaxy assembly by analyzing the median $\mu_g$ profiles in different SFG and QG \msS ranges with decreasing redshift and connecting evolution in galaxy $\mu_g$ profiles with the underlying stellar mass growth in galaxies. Since $z=1.1$, the majority of evolution in the median $\mu_g$ profiles of galaxies ($\sim$64$\%$ in SFGs and $\sim$71$\%$ in QGs) occurs throughout their stellar halo regions (2-10$R_e$). More massive galaxies assemble stellar halo material more rapidly at $0.2 \leq z \leq 1.1$. Over this period, QGs grow a larger fraction of their stellar haloes than SFGs at fixed $M_{\star}$ (factor of $\sim$1.2). Although star formation can account for the stellar halo growth observed in low-mass SFGs ($10^{9.5}M_\odot \leq M_\star < 10^{10.5}M_\odot$), high-mass SFGs ($M_\star \geq 10^{10.5}M_\odot$) and both low- and high-mass QGs require an additional assembly mechanism. Our results suggest accretion via minor mergers drives additional stellar halo growth in these galaxies. The contribution from accretion is larger in more massive galaxies (over $M_{\star} \geq 10^{9.5}M_{\odot}$), and QGs exhibit larger fractional increases to their ex-situ fractions over $0.2 \leq z \leq 1.1$ than SFGs at fixed $M_{\star}$.

Hydrogen-deficient Carbon (HdC) stars are a class of supergiants with anomalous chemical compositions, suggesting that they are remnants of CO-He white dwarf (WD) mergers. This class comprises two spectroscopically similar subclasses - dusty R Coronae Borealis (RCB) and dustless Hydrogen-deficient Carbon (dLHdC) stars. Both subclasses have a stark overabundance of $^{18}\textrm{O}$ in their atmospheres, but spectroscopic differences between them remain poorly studied. We present high-resolution ($R \approx 75000$) K-band spectra of six RCB and six dLHdC stars, including four newly discovered dLHdC stars, making this the largest sample to date. We develop a semi-automated fitting routine to measure $^{16}\textrm{O}/^{18}\textrm{O}$ ratios for this sample, tripling the number of dLHdC stars with oxygen isotope ratios measured from high resolution spectra. All six dLHdC stars have $^{16}\textrm{O}/^{18}\textrm{O}<1$, while the RCB stars have $^{16}\textrm{O}/^{18}\textrm{O}>4$. Additionally, for the first time, we find a trend of decreasing $^{16}\textrm{O}/^{18}\textrm{O}$ ratios with increasing effective temperature for HdC stars, consistent with predictions of theoretical WD merger models. However, we note that current models overpredict the low $^{16}\textrm{O}/^{18}\textrm{O}$ ratios of dLHdC stars by two orders of magnitude. We also measure abundances of C, N, O, Fe, S, Si, Mg, Na, and Ca for these stars. We observe a correlation between the abundances of $^{14}\textrm{N}$ and $^{18}\textrm{O}$ in our sample, suggesting that a fixed fraction of the $^{14}\textrm{N}$ is converted to $^{18}\textrm{O}$ in these stars via $\alpha$-capture. Our results affirm the emerging picture that the mass ratio/total mass of the WD binary determine whether an RCB or dLHdC is formed post-merger.

Water has proven to be ubiquitously detected in near-infrared (NIR) transmission spectroscopy observations of hot Jupiter atmospheres, including WASP-17b. However, previous analyses of WASP-17b's atmosphere based upon Hubble Space Telescope (HST) and Spitzer data could not constrain the water abundance, finding that sub-solar, super-solar and bimodal posterior distributions were all statistically valid. In this work, we observe one transit of the hot Jupiter WASP-17b using JWST's Near Infrared Imager and Slitless Spectrograph Single Object Slitless Spectroscopy (NIRISS SOSS) mode. We analyze our data using three independent data analysis pipelines, finding excellent agreement between results. Our transmission spectrum shows multiple H$_2$O absorption features and a flatter slope towards the optical than seen in previous HST observations. We analyze our spectrum using both PICASO+Virga forward models and free retrievals. POSEIDON retrievals provide a well-constrained super-solar $\log$(H$_2$O) abundance (-2.96$^{+0.31}_{-0.24}$), breaking the degeneracy from the previous HST/Spitzer analysis. We verify our POSEIDON results with petitRADTRANS retrievals. Additionally, we constrain the abundance of $\log$(H$^-$), -10.19$^{+0.30}_{-0.23}$, finding that our model including H$^-$ is preferred over our model without H$^-$ to 5.1 $\sigma$. Furthermore, we constrain the $\log$(K) abundance (-8.07$^{+0.58}_{-0.52}$) in WASP-17b's atmosphere for the first time using space-based observations. Our abundance constraints demonstrate the power of NIRISS SOSS's increased resolution, precision, and wavelength range to improve upon previous NIR space-based results. This work is part of a series of studies by our JWST Telescope Scientist Team (JWST-TST), in which we use Guaranteed Time Observations to perform Deep Reconnaissance of Exoplanet Atmospheres through Multi-instrument Spectroscopy (DREAMS).

The 2175Å bump is a prominent absorption feature at ultraviolet (UV) wavelengths in dust extinction and attenuation curves. Understanding the relative strength of this feature is important for accurate dust corrections at both low- and high-redshift. This feature is postulated to arise from polycyclic aromatic hydrocarbon (PAH) dust grains; however, the carrier has not been definitively established. We present results on the correlation between the 2175Å feature and PAH abundances in a spatially-resolved manner for 15 local galaxies in the PHANGS-JWST survey that have NUV and mid-IR imaging data from Swift/UVOT and JWST/MIRI, respectively. We find a moderate positive correlation between the 2175Å feature strength and PAH abundance, albeit with large intrinsic scatter. However, most of this trend can be attributed to a stronger negative correlation of both quantities with SFR surface density and specific-SFR (proxies of ionising radiation). The latter trends are consistent with previous findings that both the 2175Å carrier and PAHs are small grains that are easily destroyed by UV photons, although the proxy for PAH abundance could also be influenced by dust heating. When controlling for SFR surface density, we find weaker correlations between the 2175Å feature and PAH abundances, disfavouring a direct link. However, analyses based on spectroscopic measurements of the 2175Å feature and PAH features are required to verify our findings. No significant trends with gas-phase metallicity are found for the 2175Å feature and PAHs, however the metallicity range of our sample is limited. We provide prescriptions for the strength of the 2175Å feature and PAHs in local massive (metal-rich) galaxies with SFR surface density and specific-SFR, however the former should be used with caution since bump strengths measured from Swift/UVOT are expected to be underestimated.

The dominant processes by which galaxies replenish their cold gas reservoirs remain disputed, especially in massive galaxies. Stellar-gas kinematic misalignments offer an opportunity to study these replenishment processes. However, observed distributions of these misalignments conflict with current models of gas replenishment in early-type galaxies (ETGs), with longer relaxation timescales suggested as a possible solution. We use the EAGLE simulation to explore the relaxation of unstable misaligned gas in galaxies with masses of $M_{*}\geqslant \mathrm{10^{9.5}}$ M$_\odot$ between $0<z<1$. We extract misalignments from formation to relaxation providing a sample of $\sim3200$ relaxations. We find relaxation timescales tend to be short-duration, with median lifetimes of $\sim0.5$ Gyr, though with a notable population of unstable misalignments lasting $\gtrsim1$ Gyr. Relaxation time distributions show a log-linear relationship, with $\approx20\%$ of unstable misalignments persisting for $\gtrsim3$ torquing times. Long-lived unstable misalignments are predominantly found in galaxies with higher stellar masses, lower star-forming gas fractions, higher ongoing gas inflow, and which reside in the centres of dense environments. Mergers only cause $\sim17\%$ of unstable misalignments in EAGLE. We conclude that, at least in EAGLE, unstable kinematic misalignments are not predominantly driven by gas-rich minor mergers. Additionally, processes that significantly extend relaxation times are not dominant in the galaxy population. Instead, we see a diverse formation pathway for misalignments such as through hot halo cooling.

Estimating the redshifts of distant galaxies is critical for determining their intrinsic properties, as well as for using them as cosmological probes. Measuring redshifts spectroscopically is accurate, but expensive in terms of telescope time, hence it has become common to measure `photometric' redshifts, which are fits to photometry taken in a number of filters using templates of galaxy spectral energy distributions (SEDs). However, most photometric methods rely on optical and near-infrared (NIR) photometry, neglecting longer wavelength data in the far-infrared (FIR) and millimeter. Since the ultimate goal of future surveys is to obtain redshift estimates for all galaxies, it is important to improve photometric redshift algorithms for cases where optical/NIR fits fail to produce reliable results. For specific subsets of galaxies, in particular dusty star-forming galaxies (DSFGs), it can be particularly hard to obtain good optical photometry and thus reliable photometric redshift estimates, while these same galaxies are often bright at longer wavelengths. Here we describe a new method for independently incorporating FIR-to-millimeter photometry to the outputs of standard optical/NIR SED-fitting codes to help improve redshift estimation, in particular of DSFGs. We test our method with the H-ATLAS catalog, which contains FIR photometry from Herschel-SPIRE cross-matched to optical and NIR observations, and show that our approach reduces the number of catastrophic outliers by a factor of three compared to standard optical and NIR SED-fitting routines alone.

We present a detailed reanalysis of the atmospheric properties of WASP-19b, an ultra-hot Jupiter (1.14 M Jup, 1.41 R Jup) orbiting an active Sun-like star every 0.79 day. We reanalyze a transit and secondary eclipse of WASP-19b observed by the Hubble Space Telescope's Wide Field Camera 3 spectrograph (1.1 - 1.7 microns). When combined with Spitzer photometry at longer wavelengths, our analyses indicate the presence of water absorption features in both the planet's transmission and emission spectra, consistent with results from previously published studies. We jointly fit WASP-19b's dayside emission and transmission spectra with a retrieval model in order to constrain its atmospheric composition, and explore the effect of stellar activity on its transmission spectrum in greater depth. We also compare our dayside emission spectrum to predictions from a general circulation model, and conclude that magnetic drag appears to be relatively unimportant in shaping WASP-19b's atmospheric circulation. Lastly, we compare the size of WASP-19b's dayside water absorption feature to the population of hot Jupiters with similar measurements, and show that it is located in the transitional irradiation regime where temperature inversions first begin to emerge. As in previous studies, we find that the current observations provide relatively weak constraints on this planet's atmospheric properties. These constraints could be significantly improved by the addition of spectroscopically resolved observations at longer wavelengths with JWST/NIRSpec PRISM.

Multiphase gas, with hot ($\sim10^6$K) and cold ($\sim10^4$K) gas, is ubiquitous in astrophysical media across a wide range of scales. However, simulating multiphase gas has been a long-standing challenge, due to the large separation between the size of cold gas structures and the scales at which such gas impacts the evolution of associated systems. In this study, we introduce a new subgrid framework for such multiphase gas, MOGLI: Model for Multiphase Gas using Multifluid hydrodynamics, in multifluid AREPO. We develop this approach based on first principles and theoretical results from previous studies with resolved small-scale simulations, leading to a minimal number of free parameters in the formulation. We divide the interactions in the model into three sources: drag, turbulent mixing and cold gas growth. As part of the model, we also include two methods for estimating the local turbulent velocities, one using the Kolmogorov scaling, and the other using the local velocity gradients. We verify the different components of the framework through extensive comparison with benchmark single-fluid simulations across different simulation parameters, such as how resolved the cold gas is initially, the turbulent Mach number, spatial resolution, and random initialisation of turbulence. We test the complete scheme and a reduced version, with and without cold gas growth. We find a very good qualitative and quantitative agreement across the different simulation parameters and diagnostics for both local turbulent velocity estimation methods. We also reproduce behaviour like the cold gas survival criteria as an emergent property. We discuss the applications and possible extensions of MOGLI and demonstrate its capability by running a simulation which would be computationally prohibitive to run as a resolved single-fluid simulation.

Morphological classification of galaxies becomes increasingly challenging with redshift. We apply a hybrid supervised-unsupervised method to classify $\sim 14,000$ galaxies in the CANDELS fields at $0.2 \leq z \leq 2.4$ into spheroid, disk, and irregular systems. Unlike previous works, our method is applied to redshift bins of width 0.2. Comparison between models applied to a wide redshift range versus bin-specific models reveals significant differences in galaxy morphology beyond $z \geq 1$ and a consistent $\sim 25\%$ disagreement. This suggests that using a single model across wide redshift ranges may introduce biases due to the large time intervals involved compared to galaxy evolution timescales. Using the FERENGI code to assess the impact of cosmological effects, we find that flux dimming and smaller angular scales may lead to the misclassification of up to $18\%$ of disk galaxies as spheroids or irregulars. Contrary to previous studies, we find an almost constant fraction of disks ($\sim 60\%$) and spheroids ($\sim 30\%$) across redshifts. We attribute discrepancies with earlier works, which suggest a decreasing fraction of disks beyond $z \sim 1$, to the biases introduced by visual classification. Our claim is further strengthened by the striking agreement to the results reported by Lee et al. (2024) using an objective, unsupervised method applied to James Webb Space Telescope data. Exploring mass dependence, we observe a $\sim 40\%$ increase in the fraction of massive ($M_{\rm stellar} \geq 10^{10.5}{\rm M}_{\odot}$) spheroids with decreasing redshift, well balanced with a decrease of $\sim 20\%$ in the fraction of $M_{\rm stellar} \geq 10^{10.5}{\rm M}_{\odot}$ disks, suggesting that merging massive disk galaxies may form spheroidal systems.

Galaxy formation theory identifies superwinds as a key regulator of star formation rates, galaxy growth, and chemical enrichment. Thermal and radiation pressure are known to drive galactic-scale winds in dusty starbursting galaxies (e.g. M82), but modern numerical simulations have recently highlighted that cosmic-ray (CR) driven winds may be especially important in normal galaxies with modest star formation rate surface densities. However, CR-driven winds have yet to be conclusively observed -- leaving significant uncertainty in their detailed microphysics. We present MeerKAT radio continuum and HI spectral-line observations of one such normal galaxy, NGC 1532; a nearby ($D\sim15\,\mathrm{Mpc}$) and edge-on ($i \gtrsim 80^{\circ}$) spiral galaxy tidally interacting with its smaller elliptical companion, NGC 1531. We find magnetized, highly-ordered radio continuum loops extending $\sim10$ kpc above and below the disk; visibly connecting discrete star-forming regions in the disk with the nucleus. The deep MeerKAT HI observations place an upper limit on the column density of neutral gas coincident with the outflow to $N_\mathrm{HI} \lesssim 3 \times 10^{19}\,\mathrm{cm}^{-2}$. Unlike previously observed outflows -- for which ejected gas and dust can be traced across multiple wavelengths -- the loops in NGC 1532 show no detectable signs of dust or gas coincident with the radio emission far from the disk. We explore multiple possible mechanisms for driving this magnetic wind and favor an explanation where cosmic-ray pressure plays a significant role in launching these outflows.

Diffuse X-ray emission has been detected from a few Galactic globular clusters (GCs), whereas its nature still remains largely unclear. The GC Terzan 5 was previously found to show a significant diffuse thermal X-ray excess from its field, likely contributed by the Galactic background, and a non-thermal component described by a power-law model with photon index $\Gamma \sim 1$. With over 16 times the accumulated Chandra exposure time as in the prior study, we are motivated to reexamine and verify the diffuse X-ray emission from the field of Terzan 5, enabling constraints on its nature. We verify a significant diffuse X-ray excess from the field of Terzan 5 in the band 0.8--3 keV. After constraining the contribution from local X-ray background, we find a diffuse X-ray component that is genuinely associated with Terzan 5, which can be well described by a power-law model. More interestingly, the fitted photon indices show a significant increase from $\Gamma = 1.96 \pm 0.18$ in the inner region to $\Gamma = 3.48 \pm 0.71$ in the outer region. The diffuse X-rays can also be well fitted by a thermal bremsstrahlung model, with plasma temperatures declining from $kT \sim 3$ keV to $kT \sim 1$ keV. We suggest that synchrotron radiation from the combined pulsar winds of Terzan 5's millisecond pulsar population is a possible origin of the observed diffuse X-ray emission, but the the large steepening in the spectra cannot be produced solely by synchrotron cooling. Other radiation processes, like thermal bremsstrahlung, may also contribute to the diffuse X-rays.

We present ALMA observations of the [CI] 492 and 806$\,$GHz fine-structure lines in 25 dusty star-forming galaxies (DSFGs) at $z\,{=}\,4.3$ in the core of the SPT2349$-$56 protocluster. The protocluster galaxies exhibit a median $L^\prime_{[\text{CI}](2-1)}/L^\prime_{[\text{CI}](1-0)}$ ratio of 0.94 with an interquartile range of 0.81-1.24. These ratios are markedly different to those observed in DSFGs in the field (across a comparable redshift and 850$\,\mu$m flux density range), where the median is 0.55 with an interquartile range of 0.50-0.76, and we show that this difference is driven by an excess of [CI](2-1) in the protocluster galaxies for a given 850$\,\mu$m flux density. We estimate gas excitation temperatures of $T_{\rm ex}\,{=}\,59.1^{+8.1}_{-6.8}\,$K for our protocluster sample and $T_{\rm ex}\,{=}\,33.9^{+2.4}_{-2.2}\,$K for the field sample. Our main interpretation of this result is that the protocluster galaxies have had their cold gas driven to their cores via close-by interactions within the dense environment, leading to an overall increase in the average gas density and excitation temperature, and an elevated [CI](2-1) luminosity-to-far-infrared luminosity ratio.

Swirl-shaped flow structures have been observed throughout the solar atmosphere, in both emission and absorption, at different altitudes and locations, and are believed to be associated with magnetic structures. However, the distribution patterns of such swirls, especially their spatial positions, remain unclear. Using the Automated Swirl Detection Algorithm (ASDA), we identified swirls from the high-resolution photospheric observations, centered on Fe I 630.25 nm, of a quiet region near the Sun's central meridian by the Swedish 1-m Solar Telescope. Through a detailed study of the locations of the detected small-scale swirls with an average radius of $\sim$300 km, we found that most of them are located in lanes between mesogranules (which have an average diameter of $\sim$5.4 Mm) instead of the commonly believed intergranular lanes. The squared rotation, expansion/contraction, vector speeds, and proxy kinetic energy are all found to follow Gaussian distributions. Their rotation speed, expansion/contraction speed, and circulation are positively correlated with their radius. These results suggest that photospheric swirls at different scales and locations across the observational 56.5" $\times$ 57.5" field-of-view (FOV) could share the same triggering mechanism at preferred spatial and energy scales. A comparison with previous work suggests that the number of photospheric swirls is positively correlated with the number of local magnetic concentrations, stressing the close relation between swirls and local magnetic concentrations:the number of swirls should positively correlate with the number and strength of local magnetic concentrations.

The O- and B-type (OB-type) pulsating stars are important objects to study the structure and evolution of massive stars through asteroseismology. A large amount of data from various sky surveys provide an unprecedented opportunity to search for and study this kind of variable star. We identify 155 OB-type pulsating stars or candidates, including 38 Oe/Be stars or candidates, from the data observed by TESS, LAMOST, and GAIA, which are almost new. Among the 155 objects, 87 samples are identified as SPB stars including 37 objects with pure low-frequency and 50 objects with both low- and high-frequency pulsation, and 14 samples are identified as BCEP stars with both low- and high-frequency pulsation. The H-R diagram shows that these SPB and BCEP stars are mainly located in their instability regions and in the evolutionary stage of the main-sequence with a mass range of 2.5-20 $M_{\odot}$ and 7-20 $M_{\odot}$. Two special objects show fourier spectra similar to BCEP stars but with different positions in H-R, Period-Temperature (P-T), and Period-Luminosity (P-L) diagrams. Meanwhile, 52 other targets are identified as candidates of OB-type pulsating stars. We also derive the preliminary results of the P-L relation for SPB and BCEP stars, respectively. This work also indicates that in addition to the H-R diagram, P-T and P-L diagrams are also very useful for the classification of SPB and BCEP. Further detailed analysis of these objects can dramatically increase our understanding of theories of evolution and structure for massive OB-type pulsating stars.

This paper reports the discovery of new slowly pulsating B-type stars. Based on the photometric, spectral, and astrometric data of TESS, LAMOST, and Gaia surveys, we have found 286 new slowly pulsating B-type stars (SPB stars) and 21 candidates. Among these, 20 are Be stars or candidates with emission line profiles. It is shown that these SPB stars have luminosities between 40 and 2850 $L_{\odot}$ and effective temperatures ranging from 10000K to 21000K. Their pulsation periods are from 0.14 to 6.5 days with amplitude ranges of 0.2-20 mmag in TESS band. It is indicated that these targets follow the distribution of the SPB stars in the period-luminosity (P-L) and the period-temperature (P-T) diagrams. Their positions on the H-R diagram reveal that most of these pulsators are distributed in the instability region of SPB stars, in the main-sequence evolutionary stage, and with mass ranges of 2.5-7 $M_{\odot}$. However, there are some targets beyond the red edge of the theoretical instability region, which should be caused by the rapid rotation reducing the measured effective temperature. The discovery of these new SPB stars increases the total number by over 60\%, which are significant samples for further investigating the structure and evolution of intermediate-mass and even massive stars by asteroseismology.

The discovery of radio-quiet, X-ray thermally emitting isolated neutron stars (XINSs) in the ROSAT All-Sky Survey revealed a previously overlooked component of the neutron star population. Advancements in X-ray instrumentation and the availability of deep, wide-area optical surveys now enable us to explore XINSs at fainter X-ray fluxes and greater distances. In this study, we investigated candidates selected from the 4XMM-DR9 catalogue using XMM-Newton, focusing on long-term flux stability, spectral characterisation, and astrometry. By leveraging resources from the XMM2ATHENA project -- including updated catalogues, multiwavelength characterisation and machine learning classification -- we refined our understanding of this sample of soft X-ray emitters. Our findings enhance the characterisation of XINS candidates, laying the groundwork for more targeted investigations and future catalogue searches.

We search for transient gamma-ray emission in the energy range from 0.1-300 GeV using data from the Fermi-LAT telescope in coincidence with magnetar flares. For our analysis we look for coincidence with 15 distinct flares from 11 magnetars using two distinct time windows of $\pm$ 1 day and $\pm$ 15 days. For 14 of these flares from 10 magnetars, we do not see any statistically significant gamma-ray emission. However, we see two gamma-ray flares from one magnetar, namely 1E 1048.1-5937, with combined significance of $5\sigma$, observed after about 10 days from the peak of the X-ray flare. However, this magnetar is located close to the galactic plane (with galactic latitude of -0.52\degree) and this signal could be caused by contamination due to diffuse flux from gamma-ray sources in the galactic plane.

$\beta$ Cephei pulsating variable (BCEP) stars are the most massive pulsating variable stars in the main sequence, exhibiting both p- and g-mode pulsations. In this study, we identified 155 BCEP stars or candidates using data from TESS and Gaia, of which 83 were first confirmed as BCEP stars. They have visual magnitudes ranging from 8 to 12 mag and effective temperatures between approximately 20,000 and 30,000 K, while the parallaxes of most targets are between 0.2 and 0.6 mas. The study indicates that these BCEP stars have pulsation periods ranging from 0.06 to 0.31 days, with amplitudes ranging from 0.1 to 55.8 mmag in the TESS band. Additionally, the number of BCEP stars increases as the pulsation amplitude decreases. These targets align with the distribution region of BCEP stars in the luminosity-period (L-P) and temperature-period (T-P) diagrams. We have updated the L-P relation of BCEP stars. The Hertzsprung-Russell (H-R) diagram indicates that these targets are in the main-sequence evolutionary phase, with masses ranging from 7 to 20 $M_{\odot}$ and luminosities between 2800 and 71,000 $L_{\odot}$. They are almost in the theoretical instability region of BCEP stars but as previously reported, this region at the low-mass end (red) is not filled. The distribution of the pulsation constant indicates that the dominant pulsation periods of BCEP stars consist mainly of low-order p-mode pulsations with a high proportion of radial fundamental modes. These BCEP stars are excellent objects for enhancing our understanding of the structure and evolution of massive stars through asteroseismology.

Gas accreting planets embedded in protoplanetary disks are expected to show dust thermal emission from their circumplanetary disks (CPDs). However, a recently reported gas accreting planet candidate, AB Aurigae b, has not been detected in (sub)millimeter continuum observations. We calculate the evolution of dust in the potential CPD of AB Aurigae b and predict its thermal emission at 1.3 mm wavelength as a case study, where the obtained features may also be applied to other gas accreting planets. We find that the expected flux density from the CPD is lower than the 3-sigma level of the previous continuum observation by ALMA with broad ranges of parameters, consistent with the non-detection. However, the expected planet mass and gas accretion rate are higher if the reduction of the observed near-infrared continuum and H-alpha line emission due to the extinction by small grains is considered, resulting in higher flux density of the dust emission from the CPD at (sub)millimeter wavelength. We find that the corrected predictions of the dust emission are stronger than the 3-sigma level of the previous observation with the typical dust-to-gas mass ratio of the inflow to the CPD. This result suggests that the dust supply to the vicinity of AB Aurigae b is small if the planet candidate is not the scattered light of the star but is a planet and has a CPD. Future continuum observations at shorter wavelength are preferable to obtain more robust clues to the question whether the candidate is a planet or not.

Our peculiar velocity imprints a dipole on galaxy density maps derived from redshift surveys. The dipole gives rise to an oscillatory signal in the multipole moments of the observed power spectrum which we indicate as the finger-of-the-observer (FOTO) effect. Using a suite of large mock catalogues mimicking ongoing and future $\textrm{H}\alpha$- and $\textrm{H}\scriptstyle\mathrm{I}$-selected surveys, we demonstrate that the oscillatory features can be measured with a signal-to-noise ratio of up to 7 (depending on the sky area coverage and provided that observational systematics are kept under control on large scales). We also show that the FOTO effect cannot be erased by correcting the individual galaxy redshifts. On the contrary, by leveraging the power of the redshift corrections, we propose a novel method to determine both the magnitude and the direction of our peculiar velocity. After applying this technique to our mock catalogues, we conclude that it can be used to either test the kinematic interpretation of the temperature dipole in the cosmic microwave background or to extract cosmological information such as the matter density parameter and the equation of state of dark energy.

The CO-to-H$_2$ conversion factor ($\alpha_\mathrm{CO}$) is expected to vary with dust abundance and grain size distribution through the efficiency of shielding gas from CO-dissociation radiation. We present a comprehensive analysis of $\alpha_\mathrm{CO}$ and grain size distribution for nearby galaxies, using the PAH fraction ($q_\mathrm{PAH}$) as an observable proxy of grain size distribution. We adopt the resolved observations at 2-kpc resolution in 42 nearby galaxies, where $\alpha_\mathrm{CO}$ is derived from measured metallicity and surface densities of dust and HI assuming a fixed dust-to-metals ratio. We use an analytical model for the evolution of H$_2$ and CO, in which the evolution of grain size distribution is controlled by the dense gas fraction ($\eta$). We find that the observed level of $q_\mathrm{PAH}$ is consistent with the diffuse-gas-dominated model ($\eta=0.2$) where dust shattering is more efficient. Meanwhile, the slight decreasing trend of observed $q_\mathrm{PAH}$ with metallicity is more consistent with high-$\eta$ predictions, likely due to the more efficient loss of PAHs by coagulation. We discuss how grain size distribution (indicated by $q_\mathrm{PAH}$) and metallicity impact $\alpha_\mathrm{CO}$; we however did not obtain conclusive evidence that the grain size distribution affects $\alpha_\mathrm{CO}$. Observations and model predictions show similar anti-correlation between $\alpha_\mathrm{CO}$ and 12+log(O/H). Meanwhile, there is a considerable difference in how resolved $\alpha_\mathrm{CO}$ behaves with $q_\mathrm{PAH}$. The observed $\alpha_\mathrm{CO}$ has a positive correlation with $q_\mathrm{PAH}$, while the model-predicted $\alpha_\mathrm{CO}$ does not have a definite correlation with $q_\mathrm{PAH}$. This difference is likely due to the limitation of one-zone treatment in the model.

We consider the kinematic distances to nearby galaxies obtained by the Numerical Action Method (NAM) based on the Cosmic-flow-3 survey data. NAM distances are compared with 418 high-precision distances measured by the Tip of the Red Giant Branch (TRGB) method using the Hubble Space Telescope. We estimated the average difference <D_NAM - D_TRGB> = -0.30 +- 0.08 Mpc and the standard deviation of 1.57 Mpc. Approximately the same difference in the distance scale is obtained in comparison with less accurate distance estimates through the membership of galaxies in known groups or from the Tully-Fisher relation. We conclude that the NAM method provides distance estimates with an accuracy of 20% within the Local Volume, which is valid for ~90% of the sky, except for the regions of the Virgo cluster and the Coma-I group.

I examine images of 50 planetary nebulae (PNe) with observable post-common envelope evolution (CEE) binary central stars and find that jets are about 40 percent more common than dense equatorial outflows. Because, in some cases, energetic jets can compress an equatorial outflow and because fast jets might disperse early in the PN evolution and avoid detection, the CEE process is likelier to launch jets than to eject a dense equatorial outflow by a larger factor than 1.4. In most cases, the companion, mainly a main sequence star, launches the jets as it accretes mass from the envelope of the giant star. By CEE jets, I also refer to jets launched shortly before the onset of the CEE, likely a grazing envelope evolution phase, and shortly after the CEE. The jets and the accretion of mass by the companion before, during, and after the CEE affect envelope mass removal and the final orbital separation. Most numerical simulations of the CEE ignore jets, and those that include jets omit other processes. Despite the considerable progress in the last decade with tens of hydrodynamical simulations of the CEE, we are still far from correctly simulating the CEE. Including jets in simulations of the CEE requires heavy computer resources, but it must be the next step.

The cosmic ray (CR) flux, as well as the hydrogen flux into the atmosphere of an exoplanet, can change the composition of the atmosphere. Here, we present the CR and hydrogen flux on top of the atmosphere. To do so, we have to study the 3D multifluid MHD structure of astrospheres. We discuss the shock structure of the stellar wind of LHS 1140 using four different models: HD and MHD single-fluid models, as well as multifluid models for both cases, including a neutral hydrogen flow from the interstellar medium. The CR flux in a multifluid model as well as the ionization rate in an exoplanetary atmosphere are also presented. The astrosphere is modeled using the 3D Cronos code, while the CR flux at LHS 1140 b is calculated using both a 1D and a 3D stochastic galactic CR modulation code. Finally, the atmospheric ionization and radiation dose is estimated using the AtRIS code. Results. It is shown that the 3D multifluid positions of the termination shock differ remarkably from those found in the 3D ideal-single fluid hydrodynamic case. CR fluxes computed using a 1D approach are completely different from those calculated using the 3D modulation code and show an essentially unmodulated spectrum at the exoplanet in question. Utilizing these spectra, ionization rates and radiation exposure within the atmosphere of LHS 1140 b are derived. The termination shock, astropause, and bow shock distances must be taken from the 3D multifluid MHD model to determine the CR fluxes correctly. Moreover, because of the tiny astrosphere, the exoplanet is submerged in the neutral hydrogen flow of the interstellar medium, which will influence the exoplanetary atmosphere. A 3D approach to Galactic\0 cosmic ray (GCR) modulation in astrospheres is also necessary to avoid unrealistic estimates of GCR intensities.

In 21cm intensity mapping of the large-scale structure (LSS), regions in Fourier space could be compromised by foreground contamination. In interferometric observations, this contamination, known as the foreground wedge, is exacerbated by the chromatic response of antennas, leading to substantial data loss. Meanwhile, the baryonic acoustic oscillation (BAO) reconstruction, which operates in configuration space to "linearize" the BAO signature, offers improved constraints on the sound horizon scale. However, missing modes within these contaminated regions can negatively impact the BAO reconstruction algorithm. To address this challenge, we employ the deep learning model U-Net to recover the lost modes before applying the BAO reconstruction algorithm. Despite hardware limitations, such as GPU memory, our results demonstrate that the AI-restored 21cm temperature map achieves a high correlation with the original signal, with a correlation ratio of approximately $0.9$ at $k \sim 1 h/Mpc$. Furthermore, subsequent BAO reconstruction indicates that the AI restoration has minimal impact on the performance of the `linearized' BAO signal, proving the effectiveness of the machine learning approach to mitigate the impact of foreground contamination. Interestingly, we demonstrate that the AI model trained on coarser fields can be effectively applied to finer fields, achieving even higher correlation. This success is likely attributable to the scale-invariance properties of non-linear mode coupling in large-scale structure and the hierarchical structure of the U-Net architecture.

During their thermally pulsing phase, Asymptotic Giant Branch (AGB) stars eject material that forms extended dusty envelopes. Visible polarimetric imaging found clumpy dust clouds within two stellar radii of several oxygen-rich stars. Inhomogeneous molecular gas has also been observed in multiple emission lines within several stellar radii of different oxygen rich stars, including W Hya and Mira. At the stellar surface level, infrared images have revealed intricate structures around the carbon semi-regular variable R Scl and in the S-type star $\pi^{\mathrm{1}}$ Gru. Infrared images have also shown clumpy dust structures within a few stellar radii of the prototypical carbon AGB star IRC+10216, and studies of the molecular gas distribution beyond the dust formation zone have also shown complex circumstellar structures. Because of the lack of sufficient spatial resolution, however, the distribution of molecular gas in the stellar atmosphere and the dust formation zone of AGB carbon stars is not known, nor is how it is subsequently expelled. Here we report observations with a resolution of one stellar radius of the recently formed dust and molecular gas in the atmosphere of IRC+10216. Lines of HCN, SiS, and SiC$_2$ appear at different radii and in different clumps, which we interpret as large convective cells in the photosphere, as seen in Betelgeuse. The convective cells coalesce with pulsation causing anisotropies that, together with companions, shape its circumstellar envelope.

The giant molecular cloud complex Sagittarius B2 (Sgr~B2) in the central molecular zone of our Galaxy hosts several high-mass star formation sites, with Sgr~B2(M) and Sgr~B2(N) being the main centers of activity. This analysis aims to comprehensively model each core spectrum, considering molecular lines, dust attenuation, and free-free emission interactions. We describe the molecular content analysis of each hot core and identify the chemical composition of detected sources. Using ALMA's high sensitivity, we aim to characterize the hot core population in Sgr~B2(M) and N, gaining a better understanding of the different evolutionary phases of star formation processes in this complex. We conducted an unbiased ALMA spectral line survey of 47 sources in band 6 (211-275 GHz). Chemical composition and column densities were derived using XCLASS, assuming local thermodynamic equilibrium. Quantitative descriptions for each molecule were determined, considering all emission and absorption features across the spectral range. Temperature and velocity distributions were analyzed, and derived abundances were compared with other spectral line surveys. We identified 65 isotopologs from 48 different molecules, ranging from light molecules to complex organic compounds, originating from various environments. Most sources in the Sgr~B2 complex were assigned different evolutionary phases of high-mass star formation. Sgr~B2(N) hot cores show more complex molecules such as CH$_3$OH, CH$_3$OCHO, and CH$_3$OCH$_3$, while M cores contain lighter molecules such as SO$_2$, SO, and NO. Some sulfur-bearing molecules are more abundant in N than in M. The derived molecular abundances can be used for comparison and to constrain astrochemical models. Inner sources in both regions were generally more developed than outer sources, with some exceptions.

GRB 191019A was a long Gamma-ray burst (GRB) lasting about 65 s and, as such, originally thought to be linked to a core-collapse supernova. However, even though follow-up observations identified the optical counterpart close to the bright nucleus of a nearby ancient galaxy (z=0.248), no associated supernova was found. This led to the suggestion that the burst was caused by the merger of two compact stellar objects, likely in a dense circumnuclear environment. By using a recently developed diagnostic tool based on prompt emission temporal properties, we noticed that GRB 191019A falls among those long GRBs which are associated with compact mergers and with evidence of kilonova light. We thus re-analyzed unpublished GROND multi-color (g'r'i'z'JHK_s) data obtained between 0.4 and 15 days post trigger. Image subtraction confirmed the optical counterpart in all four optical bands, with GROND tracking its fading until 1.5 days post-burst. Incorporating publicly available Swift-XRT data, a joint fit of an afterglow plus a kilonova model revealed a better match than an afterglow-only scenario. The resulting kilonova properties resemble those of AT2017gfo associated with the binary neutron star merger GW170817, with a total ejected mass of about 0.06 solar mass. Contrary to previous findings inferring a high-density circumburst environment (n0=10^7-10^8 cm^-3), our analysis finds standard conditions (n0 = 1 cm^-3), suggesting the long duration of GRB 191019A was intrinsic rather than due to jet interaction with a dense external medium.

The BL Lac object 3C 371 is one of the targets that are regularly monitored by the Whole Earth Blazar Telescope (WEBT) Collaboration to study blazar variability on both short and long timescales. We aim to evaluate the long-term multiwavelength (MWL) behaviour of 3C 371, comparing it with the results derived for its optical emission in our previous study. For this, we make use of the multi-band campaigns organized by the WEBT Collaboration in optical and radio between January 2018 and December 2020, and of public data from Swift and Fermi satellites and the MOJAVE Very Large Interferometry programme. We evaluate the variability shown by the source in each band with the amplitude variability quantification, as well as possible interband correlation using the z-Discrete Correlation Function. We also present a deep analysis of the optical-UV, X-ray and $\gamma$-ray spectral variability. With the MOJAVE data we perform a kinematics analysis, looking for components propagating along the jet, calculating its kinematics parameters. This set of parameters is later used for the interpretation of the source MWL behaviour, modelling the broadband spectral energy distribution (SED) of the source with theoretical blazar emission scenarios.

The properties of the entire neutron star population can be inferred by modeling their evolution, from birth to the present, through pulsar population synthesis. This involves simulating a mock population, applying observational filters, and comparing the resulting sources to the limited subset of detected pulsars. We specifically focus on the magneto-rotational properties of Galactic isolated neutron stars and provide new insights into the intrinsic radio luminosity law by combining pulsar population synthesis with a simulation-based inference (SBI) technique called truncated sequential neural posterior estimation (TSNPE). We employ TSNPE to train a neural density estimator on simulated pulsar populations to approximate the posterior distribution of the underlying parameters. This technique efficiently explores the parameter space by concentrating on regions that are most likely to match the observed data thus allowing a significant reduction in training dataset size. We demonstrate the efficiency of TSNPE over standard neural posterior estimation (NPE), achieving robust inferences of magneto-rotational parameters consistent with previous studies using only around 4% of the simulations required by NPE approaches. Moreover, for the first time, we incorporate data from the Thousand Pulsar Array (TPA) program on MeerKAT, the largest unified sample of neutron stars with consistent fluxes measurement to date, to help constrain the stars' intrinsic radio luminosity. We find that adding flux information as an input to the neural network largely improves the constraints on the pulsars' radio luminosity, as well as improving the estimates on other input parameters.

Accurately determining neutrino masses is a main objective of contemporary cosmology. Since massive neutrinos affect structure formation and evolution, probes of large scale structure are sensitive to the sum of their masses. In this work, we explore future constraints on $\sum m_\nu$ utilizing line-intensity mapping (LIM) as a promising emerging probe of the density of our Universe, focusing on the fine-structure [CII] line as an example, and compare these constraints with those derived from traditional galaxy surveys. Additionally, we perform a multi-tracer analysis using velocity tomography via the kinetic Sunyaev-Zeldovich and moving lens effects to reconstruct the three-dimensional velocity field. Our forecasts indicate that the next-generation AtLAST detector by itself can achieve $\sigma_{\Sigma m_\nu} \sim 50$ meV sensitivity. Velocity tomography will further improve these constraints by 4%. Incorporating forecasts for CMB-S4 and DESI-BAO in a comprehensive multi-tracer analysis, while setting a prior on the optical depth to reionization $\tau$ derived using 21-cm forecasted observations, to break degeneracies, we find that a $\gtrsim5\sigma$ detection of $\sum m_\nu\!\sim\! 60$ meV, under the normal hierarchy, is within reach with LIM. Even without a $\tau$ prior, our combined forecast reaches $\sigma_{\Sigma m_\nu} \!\sim\! 18$ meV.

Irradiation by energetic ions, electrons, and UV photons induces sputtering and chemical processes (radiolysis) in the surfaces of icy moons, comets, and icy grains. Laboratory experiments, both of ideal surfaces and of more complex and realistic analog samples, are crucial to understand the interaction of surfaces of icy moons and comets with their space environment. This study shows the first results of mass spectrometry measurements from porous water ice regolith samples irradiated with electrons as a representative analogy to water-ice rich surfaces in the solar system. Previous studies have shown that most electron-induced H2O radiolysis products leave the ice as H2 and O2 and that O2 can be trapped under certain conditions in the irradiated ice. Our new laboratory experiments confirm these findings. Moreover, they quantify residence times and saturation levels of O2 in originally pure water ice. H2O may also be released from the water ice by irradiation, but the quantification of the released H2O is more difficult and the total amount is sensitive to the electron flux and energy.

In this work, optical observations of the nova V5584 Sgr are presented. These observations cover different phases including pre-maximum, early decline, and nebular. The spectra are dominated by hydrogen Balmer, Fe II, and O I lines with P-Cygni profiles in the early phase, which are subsequently observed in complete emission. The presence of numerous Fe II lines and low ejecta velocity aligns with the Fe II type nova classification. From optical and NIR colors it is clear that this nova manifests dust formation in the ejecta. The dust temperature and mass were estimated from a spectral energy distribution (SED) fit to the JHK band magnitudes and the WISE data. Light curve analysis shows t$_2$ and t$_3$ values of $\sim$ 26 and $\sim$ 48 days, classifying the nova as moderately fast. The physical and chemical properties during early decline and later phases were evaluated using the photoionization code CLOUDY. The best-fit model parameters from two epochs of multiwavelength spectra are compatible with a hot white dwarf source with a roughly constant luminosity of $\sim$ (2.08 $\pm$ 0.10) $\times$ 10$^{36}$ erg s$^{-1}$. We find an ejected mass of $\sim$ (1.59 $\pm$ 0.04) $\times$ 10$^{-4}$M$_{\odot}$. Abundance analysis indicates that the ejecta is significantly enriched relative to solar values, with O/H = 30.2, C/H = 10.8, He/H = 1.8, Mg/H = 1.68, Na/H = 1.55, and N/H = 45.5 in the early decline phase, and O/H = 4.5, Ne/H = 1.5, and N/H = 24.5 in the nebular phase.

The three-dimensional velocity structure of the shock-heated Si-reach and S-reach ejecta were reconstructed in Tycho supernova remnant from Doppler-shifted lines. The vector components along the line of sight were restored from the spatially resolved spectral analysis of the Doppler shifts of Si XIII and S XV lines. The components in the plane of the sky were derived from analysis of the proper motion of the remnant's edge at different azimuths. This has been done by using the data of X-ray observations from Chandra observatory as well as the radio data from the Very Large Array. Differences in Doppler velocities over the Tycho's SNR are of the order of thousands of km/s. The speed of the ejecta on the opposite sides of the remnant as a three-dimensional object differs on 20-30%. There are asymmetries and differences in the spatial distributions between the Si-reach and S-reach ejecta components. Namely, the level of isotropy is higher in Si while the vector components directed outward of the observer are larger in S. This puts limitations on the level of deviation of the internal structure of the progenitor star from the ideal layered structure as well as on the level of asymmetries in supernova explosion.

Since the first detection by the DASI experiment in 2002, measurements of the polarization of the cosmic microwave background (CMB) have grown into an important role in testing our understanding of conditions in the early universe and cosmology. The field has seen rapid experimental progress, driven in large part by the desire to make increasingly precise measurements of CMB polarization. Precise measurements of the CMB polarization anisotropies contain as much information as the CMB temperature anisotropy, and promise to unlock new tests of physics and the standard cosmological model. In this chapter, we first discuss how polarization is produced through Thomson scattering, after which types of polarization patterns are connected to the cosmological sources. Finally, we briefly discuss the experimental hardware that enables these measurements.

We present a comprehensive analysis of globular cluster (GC) formation and evolution across the $34^3$ Mpc$^3$ volume of the E-MOSAICS galaxy formation simulations. Defining GCs as surviving, high-mass ($>10^5$ M$_\odot$) clusters, we analyse their formation histories as a function of their metallicity and host galaxy mass, also distinguishing between central and satellite galaxies. The redshift of peak GC formation rate increases weakly with galaxy mass, decreases with metallicity, and does not differ between centrals and satellites. The epoch of peak GC formation precedes that of the stars by a factor of $1.1{-}1.6$, primarily due to `downsizing', i.e. low-mass galaxies form their stars later. Consequently, this offset decreases with galaxy mass, leading to nearly coeval stellar and GC populations in massive galaxies ($>10^{11}$ M$_\odot$). GCs themselves do not exhibit strong downsizing, because they predominantly formed at early cosmic epochs conducive to the formation (through high gas pressures) and survival (through high galaxy merger and GC migration rates) of massive, compact stellar systems. The total GC formation rate in the volume peaks at $z\approx 2.5$, shortly before star formation peaks at $z\approx 2$, but well after the general cluster formation rate at $z\approx 4$, reflecting a survivor bias where surviving GCs formed more recently. We find that GC formation commenced early, at $z>10$, such that the results of this work may provide a framework for interpreting direct observations of proto-GC formation with the JWST, especially as these observations accumulate to enable statistical studies.

The Main Belt asteroid (203) Pompeja shows evidence of extreme variability in visible and near-infrared spectral slope with time. The observed spectral variability has been hypothesized to be attributed to spatial variations across Pompeja's surface. In this scenario, the observed spectrum of Pompeja is dependent on the geometry of the Sun and the observer relative to the asteroid's spin pole and surface features. Knowledge of the rotational spin pole and shape can be gleaned from light curves and photometric measurements. However, dense light curves of Pompeja are only available from two apparitions. Further, previous estimates of Pompeja's sidereal period are close to being Earth-commensurate, making ground-based light curves difficult to obtain. To overcome these difficulties, we implement a pipeline to extract a dense light curve of Pompeja from cutouts of TESS Full Frame Images. We succeeded in obtaining a dense light curve of Pompeja covering $\sim$22 complete rotations. We measure a synodic period of $P_{syn} =24.092 \pm 0.005$ hours and amplitude of 0.073 $\pm$ 0.002 magnitudes during Pompeja's 2021 apparition in the TESS field of view. We use this light curve to refine models of Pompeja's shape and spin pole orientation, yielding two spin pole solutions with sidereal periods and spin pole ecliptic coordinates of $P_{\mathrm{sid}, 1} = 24.0485 \pm 0.0001$ hours, $\lambda_1 = 132^{\circ}$, $\beta_1 = +41^{\circ}$ and $P_{\mathrm{sid}, 2} = 24.0484 \pm 0.0001$ hours, $\lambda_2 =307^{\circ}$, $\beta_2 =+34^{\circ}$. Finally, we discuss the implications of the derived shape and spin models for spectral variability on Pompeja.

The measurements of the CMB have determined the cosmological parameters with high accuracy, and the observation of the flatness of space have contributed to the status of the concordance $\Lambda$CDM model. However, the cosmological constant $\Lambda$, necessary to close the model to critical density, remains an open conundrum. We explore the observed late-time accelerated expansion of the Universe, where we consider that the Friedmann equation describes the expansion history of FLRW universes in the local reference frame of freely falling comoving observers, which perceive flat, homogeneous and isotropic space in their local inertial system, where, as a consequence of the equivalence principle, special relativity applies. We use this fact to propose an extension to $\Lambda$CDM, incorporating the initial conditions of the background universe, comprising the initial energy densities as well as the initial post big bang expansion rate. The observed late-time accelerated expansion is then attributed to a kinematic effect akin to a dark energy component. Choosing the same $\Omega_{m,0} \simeq 0.3$ as $\Lambda$CDM, its equation of state $w_{de} \simeq -0.8$. Furthermore, we include the impact on the expansion history caused by the cosmic web of the late Universe, once voids dominate its volume, and find that the initially constant $w_{de}$ becomes time-dependent, evolving to a value of $w_{de} \simeq -0.9$ at the present. While this impact by voids is minor, it is sufficient to provide a solution to the Hubble tension problem. We use CLASS to calculate the expansion history and power spectra of our extension and compare our results to concordance $\Lambda$CDM and to observations. We find that our model agrees well with current data, in particular with the final data release PR4 of the Planck mission, where it explains the reported spatial curvature of $\Omega_{k,0} = - 0.012 \pm 0.010$.

The combination of detection techniques enhances our ability to identify companions orbiting nearby stars. We employed high-contrast imaging to constrain mass and separation of possible companions responsible for the significant proper motion anomalies of the nearby stars HIP 11696, HIP 47110 and HIP 36277. These targets were observed using the LBT's high-contrast camera, SHARK-NIR, in H-band using a Gaussian coronagraph, and with the LMIRCam instrument in the L'-band and using a vAPP coronagraph. Both observations were conducted simultaneously. Additionally, constraints at short separations from the host star are derived analyzing the renormalized unit weight error (RUWE) values from the Gaia catalogue. We find that the companion responsible for the anomaly signal of HIP 11696 is likely positioned at a distance from 2.5 to 28 astronomical units from its host. Its mass is estimated to be between 4 and 16 Jupiter masses, with the greater mass possible only at the upper end of the separation range. Similar limits were obtained for HIP 47110 where the companion should reside between 3 and 30 au with a mass between 3 and 10 MJup. For HIP 36277, we identified a faint stellar companion at large separation, though it might be substellar depending on the assumed age for the star. Considering the older age, this object accounts for the absolute value of the PMa vector but not for its direction. Additionally, we found a substellar candidate companion at a closer separation that could explain the PMa signal, considering a younger age for the system.

We introduce the Bias-free Extragalactic Analysis for Cosmic Origins with NIRCam (BEACON) survey, a JWST Cycle2 program allocated up to 600 pure-parallel hours of observations. BEACON explores high-latitude areas of the sky with JWST/NIRCam over $\sim100$ independent sightlines, totaling $\sim0.3$deg$^2$, reaching a median F444W depth of $\approx28.2$AB mag (5$\sigma$). Based on existing JWST observations in legacy fields, we estimate that BEACON will photometrically identify 25--150 galaxies at $z>10$ and 500--1000 at $z\sim7$--10 uniquely enabled by an efficient multiple filter configuration spanning $0.9$--5.0$\mu$m. The expected sample size of $z>10$ galaxies will allow us to obtain robust number density estimates and to discriminate between different models of early star formation. In this paper, we present an overview of the survey design and initial results using the first 19 fields. We present 129 galaxy candidates at $z>7$ identified in those fields, including 11 galaxies at $z>10$ and several UV-luminous ($M_{\rm UV}<-21$mag) galaxies at $z\sim8$. The number densities of $z<13$ galaxies inferred from the initial fields are overall consistent with those in the literature. Despite reaching a considerably large volume ($\sim10^5$Mpc$^3$), however, we find no galaxy candidates at $z>13$, providing us with a complimentary insight into early galaxy evolution with minimal cosmic variance. We publish imaging and catalog data products for these initial fields. Upon survey completion, all BEACON data will be coherently processed and distributed to the community along with catalogs for redshift and other physical quantities.

We have obtained near-infrared ($0.80-2.45\mu$m) spectra of the recurrent nova LMCN 1968-12a on two occasions during its 2024 August eruption. This is the first near-infrared spectroscopy of an extragalactic nova. The initial spectrum, on day 8.48, caught the nova in the coronal phase, with the [SiX] $1.43\mu$m line being extremely strong. This line had a luminosity of $\sim95$L$_\odot$, and is clearly a very powerful coolant. Its presence, together with the absence of [SiIX] 1.56$\mu$m, implies a coronal temperature $\gtrsim3\times10^6$K, possibly amongst the highest recorded coronal temperature in a nova eruption. With the exception of the [SiX] line, the near-infrared spectra are remarkable for being devoid of metal lines. We suggest that this is due, in part, to the exceptionally high temperature of the coronal gas, causing ions, whose emission lines would normally appear in the near-infrared spectrum, to be collisionally ionised to higher stages.

The dynamics of giant planet magnetospheres is controlled by a complex interplay between their fast rotation, their interaction with the solar wind, and their diverse internal plasma and momentum sources. In the ionosphere, the Hall and Pedersen conductances are two key parameters that regulate the intensity of currents coupling the magnetosphere and the ionosphere, and the rate of angular momentum transfer and power carried by these currents. We perform a comparative study of Hall and Pedersen conductivities and conductances in the four giant planets of our Solar System - Jupiter, Saturn, Uranus and Neptune. We use a generic ionospheric model (restraining the studied ions to H3+, CH5+, and meteoric ions) to study the dependence of conductances on the structure and composition of these planets' upper atmospheres and on the main ionization sources (photoionization, ionization by precipitating electrons, and meteoroid ablation). After checking that our model reproduces the conclusions of Nakamura et al. (2022, this https URL) at Jupiter, i.e. the contribution of meteoric ions to the height-integrated conductances is non-negligible, we show that this contribution could also be non-negligible at Saturn, Uranus and Neptune, compared with ionization processes caused by precipitating electrons of energies lower than a few keV (typical energies on these planets). However, because of their weaker magnetic field, the conductive layer of these planets is higher than the layer where meteoric ions are mainly produced, limiting their role in magnetosphere-ionosphere coupling.

We explore the radio emission of spectroscopically confirmed, X-ray weak, Broad Line AGN (BLAGN, or type 1) selected with JWST in the GOODS-N field, one of the fields with the best combination of deep radio observations and statistics of JWST-selected BLAGN. We use deep radio data at different frequencies (144\,MHz, 1.5\,GHz, 3\,GHz, 5.5\,GHz, 10\,GHz), and we find that none of the 22 sources investigated is detected at any of the aforementioned frequencies. Similarly, the radio stacking analysis does not reveal any detection down to an rms of $\sim 0.2\mu$Jy beam$^{-1}$, corresponding to a $3\sigma$ upper limit at rest frame 5 GHz of $L_{5GHz}=2\times10^{39}$ erg s$^{-1}$ at the mean redshift of the sample $z\sim 5.2$. We compared this and individual sources upper limits with expected radio luminosities estimated assuming different AGN scaling relations. For most of the sources the radio luminosity upper limits are still compatible with expectations for radio-quiet (RQ) AGN; nevertheless, the more stringent stacking upper limits and the fact that no detection is found would suggest that JWST-selected BLAGN are weaker than standard AGN even at radio frequencies. We discuss some scenarios that could explain the possible radio weakness, such as free-free absorption from a dense medium, or the lack of either magnetic field or a corona, possibly as a consequence of super-Eddington accretion. These scenarios would also explain the observed X-ray weakness. We also conclude that $\sim$1 dex more sensitive radio observations are needed to better constrain the level of radio emission (or lack thereof) for the bulk of these sources. The Square Kilometer Array Observatory (SKAO) will likely play a crucial role in assessing the properties of this AGN population.

Recent simulations of Be stars in misaligned binary systems have revealed that misalignment between the disc and binary orbit can cause the disc to undergo Kozai-Lidov (KL) oscillations or disc-tearing. We build on our previous suite of three-dimensional smoothed particle hydrodynamics simulations of equal-mass systems by simulating eight new misaligned Be star binary systems, with mass-ratios of 0.1 and 0.5, or equal-mass systems with varying viscosity. We find the same phenomena occur as previously for mass ratios of 0.5, while the mass ratio of 0.1 does not cause KL oscillations or disc-tearing for the parameters examined. With increased viscosity in our equal-mass simulations, we show that these phenomena and other oscillations are damped out and do not occur. We also briefly compare two viscosity prescriptions and find they can produce the same qualitative disc evolution. Next, we use the radiative transfer code HDUST to predict observable trends of a KL oscillation, and show how the observables oscillate in sync with disc inclination and cause large changes in the polarization position angle. Our models generate highly complex line profiles, including triple-peak profiles that are known to occur in Be stars. The mapping between the SPH simulations and these triple-peak features gives us hints as to where they originate. Finally, we construct interferometric predictions of how a gap in the disc, produced by KL oscillations or disc-tearing, perturbs the visibility versus baseline curve at multiple wavelengths, and can cause large changes to the differential phase profile across an emission line.

The production mechanism of fast radio bursts (FRBs) remains elusive, and potential correlations between burst occurrence times and various burst properties may offer important clues. Among them, the spectral peak frequency is particularly important because it may encode direct information about the physical conditions and environment at the emission site. Analyzing over 4,000 bursts from the three most active sources -- FRB 20121102A, FRB 20201124A, and FRB 20220912A -- we measure the two-point correlation function $\xi(\Delta t, \Delta\nu_{\:\rm peak}\:)$ in the two-dimensional space of time separation $\Delta t$ and peak frequency shift $\Delta\nu_{\:\rm peak}\:$ between burst pairs. We find a universal trend of asymmetry about $\Delta\nu_{\:\rm peak}\:$ at high statistical significance; $\xi(\Delta\nu_{\:\rm peak}\:)$ decreases as $\Delta\nu_{\:\rm peak}\:$ increases from negative to positive values in the region of short time separation ($\Delta t < 0.3$ s), where physically correlated aftershock events produce a strong time correlation signal. This indicates that aftershocks tend to exhibit systematically lower peak frequencies than mainshocks, with this tendency becoming stronger at shorter $\Delta t$. We argue that the "sad trombone effect" -- the downward frequency drift observed among sub-pulses within a single event -- is not confined within a single event but manifests as a statistical nature that extends continuously to independent yet physically correlated aftershocks with time separations up to $\Delta t \sim 0.3$ s. This discovery provides new insights into underlying physical processes of repeater FRBs.

A recent paper (ref. 1) used infrared images of Io acquired by the Juno/JIRAM instrument to derive a latitudinal dependence of the spectral radiance and conclude that such latitudinal dependence is consistent with a magma ocean model. We challenge their conclusions, and we draw attention to some potential issues with their analysis. In this letter, we will use three arguments to show that: (1) the (ref. 1) paper uses saturated data; (2) the M-filter of the JIRAM imager is only a weak and incomplete proxy for the total power output; and finally (3) even assuming that the radiance was correctly estimated, the latitudinal dependence of the 4.8-$\mu$m spectral radiance is not statistically significant. These facts, taken together, demonstrate that the results presented in (ref. 1) are not sufficient to confirm consistency with a magma ocean model on Io.

The neutral sodium resonance doublet (Na i D) has been detected in the upper atmosphere of several close-in gas giants, through high-resolution transmission spectroscopy. We aim to investigate whether its variability is linked to the planets' properties, the data quality, or the accuracy of the system parameters used. Using the public code SLOPpy, we extracted the transmission spectrum in the Na i D region of ten gas giants for which a large number of HARPS-N observations are available. We modelled the absorption signals found, performing an MCMC analysis, and converted the measured absorption depth to the corresponding atmospheric height over which most sodium absorption occurs. While two targets (GJ 436 b and KELT-7 b) show no Na i D feature, we found variability in the transmission spectrum of the other targets. Three of them (HD 209458 b, WASP-80 b, and WASP-127 b) present absorption on only some nights, while in the other five targets (HD 189733 b, KELT-9 b, KELT-20 b, WASP-69 b, and WASP-76 b), a significant absorption signal is present on most of the nights analysed. Except for WASP-69 b, the measured absorption depths lead to a ratio of the two Na I D depths that is compatible with or slightly larger than one. As was expected from literature, the relative atmospheric height follows an empirical exponential trend as a function of a scaled product of the planet's equilibrium temperature and surface gravity. We confirm the sodium detection on HD 189733 b, KELT-9 b, KELT-20 b, WASP-69 b, and WASP-76 b. The signal detected in WASP-127 b requires further observations for definitive confirmation. We exclude a planetary origin for the signals found on HD 209458 b and WASP-80 b. The sodium absorption variability does not appear to be related to planetary properties, but rather to data quality, sub-optimal data treatment, or stellar activity.

The accurate forecasting of solar flares is considered a key goal within the solar physics and space weather communities. There is significant potential for flare prediction to be improved by incorporating topological fluxes of magnetogram datasets, without the need to invoke three-dimensional magnetic field extrapolations. Topological quantities such as magnetic helicity and magnetic winding have shown significant potential towards this aim, and provide spatio-temporal information about the complexity of active region magnetic fields. This study develops time-series that are derived from the spatial fluxes of helicity and winding that show significant potential for solar flare prediction. It is demonstrated that time-series signals, which correlate with flare onset times, also exhibit clear spatial correlations with eruptive activity; establishing a potential causal relationship. A significant database of helicity and winding fluxes and associated time series across 144 active regions is generated using SHARP data processed with the ARTop code that forms the basis of the time-series and spatial investigations conducted here. We find that a number of time-series in this dataset often exhibit extremal signals that occur 1-8 hours before a flare. This, publicly available, living dataset will allow users to incorporate these data into their own flare prediction algorithms.

Aims: We investigate the impact of galaxy mergers on the Lyman Continuum (LyC) radiation escape, fesc, from high-redshift galaxies. Methods: We post-process ~ 6e5 galaxies (redshift 5.2 < z < 10) extracted from the TNG50 cosmological simulation using a physically motivated analytic model for LyC escape. Results: Galaxies that have not experienced a merger for the last ~ 700 Myr have an average fesc ~ 3%, which increases to up to 14% immediately following a merger. The strongest effect can be observed in galaxies with stellar masses of ~ 1e7 Msun. We attribute the increase in the escape fraction to two main factors: (i) accretion of metal-poor gas onto the central region of a galaxy, which feeds star formation and LyC emission; and (ii) displacement of neutral gas relative to star-forming regions, which reduces the optical depth to LyC photons. We additionally examine how proximity to other galaxies influences LyC escape, finding that galaxies with more neighbors tend to have more frequent mergers, and thus a higher LyC leakage. However, galaxies in overdense regions tend to have a larger LyC escape fraction independently from mergers, because of their higher gas inflow, and consequent increase in the star formation rate. The increase in both mergers and gas inflow could contribute to low-mass galaxies ionizing proximity zones of high-z Ly-alpha leakers recently observed with JWST.

We examine images of the protoplanetary disk 114--426 with JWST/NIRCam in 12 bands. This large disk is oriented edge-on with a dark midplane flanked by lobes of scattered light. The outer edges of the midplane are seen in silhouette against the Orion Nebula, providing a unique opportunity to study planet-forming material in absorption. We discover a dip in the scattered light of the disk at 3\,$\micron$ -- compelling evidence for the presence of water ice. The 3\,$\micron$ dip is also seen in the silhouette of the disk, where we quantify the ice abundance with models of pure absorption and avoid the complications of disk scattering effects. We find grain ice-to-refractory mass ratios of up to $\sim$0.2, maximum grain sizes of 0.25 to 5\,$\micron$, and a total dust plus ice mass of 0.46\,$M_\oplus$ in the silhouette region. We also discover excess absorption in the NIRCam bands that include the Paschen $\alpha$ line, suggesting there may be excited atomic hydrogen in the disk. Examining the morphology of the scattered light lobes reveals that they are laterally offset from each other and exhibit a brightness asymmetry that flips with wavelength -- both evidence for a tilted inner disk in this system.

Exoplanet atmospheric modeling is advancing from chemically diverse one-dimensional (1D) models to three-dimensional (3D) global circulation models (GCMs), which are crucial for interpreting observations from facilities like the James Webb Space Telescope (JWST) and Extremely Large Telescope (ELT). However, maintaining chemical diversity in models, especially in GCMs, is computationally expensive, limiting their complexity. Optimizing the number of reactions and species can address this tradeoff, but transparent and efficient methods for such optimization are lacking in current exoplanet literature. We aim to develop a systematic approach for reducing chemical networks in exoplanetary atmospheres while balancing accuracy and computational efficiency. Our data-driven method selects optimal reduced chemical networks based on accuracy and computational efficiency metrics. This approach can optimize networks for similar planets simultaneously, assign weights to prioritize accuracy or efficiency, and is applicable when including photochemistry. We base our method on sensitivity analysis of a typical 1D chemical kinetics model, applying principal component analysis to the sensitivities. To achieve fast and reliable network reduction, we utilize a genetic algorithm, a machine-learning optimization method that mimics natural selection. We present three schemes tailored for different priorities (accuracy, computational efficiency, and adaptability to photochemistry) that demonstrate improved performance and reduced computational costs. Our genetic algorithm-based method, the first to reduce a chemical network including photochemistry in exoplanet research, offers a versatile and efficient approach to enhance both accuracy and computational efficiency.

Recent studies based on the relativistic mean field (RMF) model found certain nuclear empirical parameters, in particular the nucleon effective mass, to be strongly correlated with observable properties of Neutron Stars (NSs), such as the frequencies of $f-$mode oscillations. This shows the potential to constrain the values of effective mass from future observations of $f-$modes. One of our primary goals of this work is to investigate whether such correlations are physical or an artifact of the underlying nuclear model. To test this, we perform a comparative study of the correlations between NS astrophysical observables and nuclear physics parameters using two different equation of state models based on RMF theory and non-relativistic Meta-Modelling (MM) scheme. The nuclear meta-model does not assume any underlying nuclear model and therefore allows us to test the model dependence of the results. The calculations of the $f-$mode characteristics are performed within the relativistic Cowling approximation. We use state-of-the-art nuclear microscopic calculations at low density and multi-messenger astrophysical data at high-density within a Bayesian-inspired scheme to constrain the parameter space of the nuclear models. From the posterior distribution, we probe the underlying correlations among nuclear parameters and with NS observables. We find that the correlation between the symmetry energy and its slope is physical, while that of the nucleon effective mass with NS observables is model-dependent. The study shows that the effective mass governs the high density behaviour in RMF models, while in the MM it is controlled by the higher order saturation parameters, and hence the effective mass in the MM is not correlated with NS observables. The findings of this investigation are interesting both for astrophysics as well as nuclear physics communities.

The thermal Sunyaev-Zel'dovich (tSZ) effect arises from inverse Compton scattering of low energy photons onto thermal electrons, proportional to the integrated electron pressure, and is usually observed from galaxy clusters. However, we can expect that the Epoch of Reionization (EoR) also contributes to this signal, but that contribution has not been previously evaluated. In this work we analyse a suite of fully-coupled radiation-hydrodynamics simulations based on RAMSES-CUDATON to calculate and study the tSZ signal from the Reionization Epoch. We construct lightcones of the electron pressure in the intergalactic medium for $6<z$ to calculate the resulting Compton y-parameters. We vary the box sizes, resolutions and star formation parameters to investigate how these factors affect the tSZ effect. We produce plots of maps and distributions of y, as well as angular temperature power spectra of the tSZ signal obtained from integrating the lightcones constructed for each simulation. We find that the tSZ signal from reionization is generally sub-dominant to the post-reionization one at larger scales ($\ell< 10^4$), but can contribute non-trivially and potentially contaminate the measured signals. At scales probed by current experiments like SPT ($\ell\sim10^3-10^4$), we find that the tSZ signal power spectrum from reionization contributes at roughly a percent level compared to the current templates, with the quadratic Doppler effect contributing an additional $\sim10\%$ to the tSZ signal. At smaller scales the tSZ from reionization peaks and can potentially dominate the total signal and is thus a potentially much more important contribution to take into account in any future, more sensitive experiments.

The red supergiant (RSG) problem, which describes the apparent lack of high-luminosity progenitors detected in Type II supernova (SN) pre-images, has been a contentious topic for two decades. We re-assess this problem using a new RSG population of the Milky Way supplemented with RSGs from other galaxies in the Local Group. In particular, we quantify the uncertainties inherent to assumptions made regarding the star's temperature or spectral type and the corresponding bolometric correction. We find that only M3 or later RSGs reproduce the steepness seen from the SN II pre-imaged sample. To assess the significance of the RSG problem, we build a metallicity-weighted cumulative luminosity distribution of M3 or later RSGs and directly compare it to the luminosity distribution of SN II pre-imaged progenitors. We find no evidence of missing high-luminosity pre-imaged progenitors since the uncertainties on the pre-imaged SN progenitors and single-band derived luminosity are too large to meaningfully infer population differences.

The coldest metal-poor population made of T and Y dwarfs are archaeological tracers of our Galaxy because they are very old and have kept the pristine material. The optical properties of these objects are important to characterise their atmospheric properties. We aim at characterising further the optical properties of ultracool metal-poor population with deep far-red optical images and parallax determinations. We solve trigonometric parallaxes of five metal-poor T dwarf candidates using 2-year monitoring with Calar-Alto 3.5-m telescope. We obtain $z'$-band photometry for the other 12 metal-poor T dwarf candidates using the 10.4-m GTC, the 8.2-m VLT, and the DES, increasing the sample of T subdwarfs with optical photometry from 12 to 24. We report a 3-$\sigma$ limit for the Accident in five optical bands using the 10.4-m GTC. We confirm four T subdwarfs and the Accident as a Y subdwarf, and propose two more Y subdwarf candidates. We emphasise that the $z_{PS1}-W1$ colour combining with the $W1-W2$ colour could break the metallicity-temperature degeneracy for T and possibly for Y dwarfs. The $z_{PS1}-W1$ colour shifts redward when metallicity decreases for a certain temperature, which is not predicted by state-of-the-art ultracool models. The Accident has the reddest $z_{PS1}-W1$ colour among our sample. The $z_{PS1}-W1$ colour will be useful to search for other examples of this cold and old population in upcoming and existing deep optical and infrared large-area surveys.

We conducted a review of the fundamental aspects of describing and detecting the Baryon Acoustic Oscillation (BAO) feature in galaxy surveys, emphasizing the optimal tools for constraining this probe based on the type of observation. Additionally, we included new results with two spectroscopic datasets to determine the best-fit model for the power spectrum, $P(k)$. Using the framework described in a previous analysis, we applied this to a different sub-sample of the BOSS survey, specifically galaxies with redshifts $0.3<z<0.65$. We also examined the eBOSS dataset with redshifts $0.6<z<1.0$, adjusting the number of parameters in the traditional polynomial fit to account for the higher redshift range. Our results showed that the dilation scale parameter $\alpha$ derived from the BOSS dataset had smaller error bars compared to the eBOSS dataset, attributable to the larger number of luminous red galaxies (LRGs) in the BOSS sample. We also compared our findings with other surveys such as WiggleZ, DES Y6, and DESI III, noting that photometric surveys typically yield larger error bars due to their lower precision. The DESI III results were in good agreement with ours within $1\sigma$, with most bins close to unity. The variation of $\alpha$ with respect to the redshift is an unresolved issue in the field, appearing in both three-dimensional and angular tomographic analyses.

Three-dimensional hydrodynamical simulations of common envelope evolution are often terminated soon after the initial dynamical plunge of the companion transitions into a long-lasting post-dynamical inspiral with slowly varying semi-major axis, $a_\text{b}$. This premature termination is often due to insufficient numerical resolution and challenges associated with the softening of the gravitational potential of the two cores. In this work, we use statically-refined 3D hydrodynamical simulations to study binaries orbiting inside a common envelope, exploring the effects of varying numerical resolution, $\delta$, gravitational potential softening prescriptions, and the associated softening lengthscale, $\epsilon$. We find that quantities such as the binary inspiral timescale or the volume-averaged shearing rate typically converge to asymptotic values only for $\epsilon \le 0.1 a_\text{b}$ and $\delta \le 6 \times 10^{-3}a_\text{b}$ with smaller $\epsilon$ requiring correspondingly smaller $\delta$. After a few tens of binary orbits, the two cores become surrounded by a corotating, nearly hydrostatic gas structure, resembling the shared envelope of a contact binary. We propose that this structure is responsible for the slowing down of the dynamical inspiral, leading to an asymptotic inspiral timescale of approximately $10^5$ orbital periods for a binary mass ratio $q=1/3$, and approximately $10^6$ orbital periods for a binary mass ratio $q=1$. By investigating kinetic helicity, we argue that the magnetic field is unlikely to organize into large-scale structures via the usual $\alpha$--effect during the post-dynamical phase. Even in the absence of magnetic fields, we observe intermittent polar outflows collimated by partially centrifugally evacuated polar funnels. (abridged)

Recent results from Type Ia Supernovae (SNe), baryon acoustic oscillations (BAO), and the cosmic microwave background (CMB) indicate 1) potentially discrepant measurements of the matter density $\Omega_m$ and Hubble constant $ H_0 $ in a $\Lambda$CDM model when data are analyzed individually, and 2) hints of dynamical dark energy in a $w_0w_a$CDM model when data are combined in a joint analysis. We examine whether underlying dynamical dark energy cosmologies favored by real data would result in biases in $\Omega_m$ and $ H_0 $ for each probe when analyzed individually in a $\Lambda$CDM framework. We generate mock datasets in $w_0w_a$CDM cosmologies, fit the individual probes under the $\Lambda$CDM model, and find that expected biases in $\Omega_m$ are $\sim 0.03$. Notably, the $\Omega_m$ differences between probes are consistent with the values observed in the real datasets. We also observe that mock DESI BAO datasets generated in the $ w_0w_a $CDM cosmologies will lead to a biased measurement of $ H_0 $ higher by ($\sim1.2$km/s/Mpc) when fitted under $\Lambda$CDM, appearing to mildly improve the Hubble tension, but as the true underlying $H_0$ is lower, the tension is in fact worsened. We find that the $\Omega_m$ discrepancies, the high BAO $ H_0 $ relative to CMB, and the joint dynamical dark energy signal itself are all related effects that could be explained simultaneously with either new physics or new systematics. While we find it is possible to unite many of the discrepancies seen in recent analyses along a single axis, our results underscore the importance of understanding systematic differences in datasets, as they have unique impacts in different cosmological parameter spaces.

Hot dust-obscured galaxies (Hot DOGs), are a family of hyper-luminous, heavily obscured quasars. A number of studies have shown that these objects reside in significantly overdense regions of the Universe based on the identification of companions at optical through far-IR wavelengths. Here we present further characterization of their environments by studying the surface density of Lyman break galaxy (LBG) candidates in the vicinity of three Hot DOGs. For two of them, WISE J041010.60-091305.2 at z=3.631 and WISE J083153.25+014010.8 at z=3.912, we identify the candidate LBG companions using deep observations obtained with Baade/IMACS. For the third, WISE J224607.56-052634.9 at z=4.601, we re-analyse previously published data obtained with Gemini-S/GMOS-S. We optimise the LBG photometric selection criteria at the redshift of each target using the COSMOS2020 catalog. When comparing the density of LBG candidates found in the vicinity of these Hot DOGs with that in the COSMOS2020 catalog, we find overdensities of $\delta=1.83\pm 0.08$ ($\delta' = 7.49\pm 0.68$), $\delta=4.67\pm 0.21$ ($\delta' = 29.17\pm 2.21$), and $\delta = 2.36\pm 0.25$ ($\delta' = 11.60\pm 1.96$) around W0410-0913, W0831+0140, and W2246-0526, respectively, without (with) contamination correction. Additionally, we find that the overdensities are centrally concentrated around each Hot DOG. Our analysis also reveals that the overdensity of the fields surrounding W0410-0913 and W0831+0140 declines steeply beyond physical scales of $\sim$2 Mpc. If these overdensities evolve to clusters by z=0, these results suggest that the Hot DOG may correspond to the early formation stages of the brightest cluster galaxy. We were unable to determine if this is also the case for W2246-0526 due to the smaller field of view of the GMOS-S observations. Our results imply that Hot DOGs may be excellent tracers of protoclusters.

We present a pattern emerging from stellar obliquity measurements in single-star systems: planets with high planet-to-star mass ratios ($M_{\rm p}/M{_*}$$>$ $2\times10^{-3}$) -- such as super-Jupiters, brown dwarf companions, and M-dwarfs hosting Jupiter-like planets -- tend to be aligned, even around hot stars. This alignment represents a 3.7$\sigma$ deviation from the obliquity distribution observed in systems with lower mass ratios ($M_{\rm p}/M{_*}$$<$ $2\times10^{-3}$), which predominantly include Jupiters and sub-Saturns. The only known outlier system, XO-3, exhibits misalignment confirmed via our newly collected Rossiter-McLaughlin effect measurement ($\lambda=41.8^{+2.1}_{-2.0}$ degrees). However, the relatively large $\textit{Gaia}$ Renormalized Unit Weight Error (RUWE) of XO-3 suggests that it may harbor an undetected binary companion, potentially contributing to its misalignment. Given that tidal realignment mechanisms are weak for hot stars, the observed alignment in high $M_{\rm p}/M{_*}$ systems is likely $\textit{primordial}$ rather than resulting from tidal interactions. One possible explanation is that only dynamically isolated planets can continue accreting gas and evolve into super-Jupiters while maintaining their primordial alignment. Conversely, planets formed in compact configurations may be unable to grow beyond the gap-opening mass, for which our work suggests an empirical boundary $M_{\rm p}/M{_*}$$=$ $2\times10^{-3}$, identified between aligned high $M_{\rm p}/M{_*}$ systems and misaligned low $M_{\rm p}/M{_*}$ systems, with dynamical instabilities contributing to the diverse spin-orbit misalignments observed in the latter.

The next generation of galaxy surveys has the potential to substantially deepen our understanding of the Universe. This potential hinges on our ability to rigorously address systematic uncertainties. Until now, diagnosing systematic effects prior to inferring cosmological parameters has been out of reach in field-based implicit likelihood cosmological inference frameworks. As a solution, we aim to diagnose a variety of systematic effects in galaxy surveys prior to inferring cosmological parameters, using the inferred initial matter power spectrum. Our approach is built upon a two-step framework. First, we employ the Simulator Expansion for Likelihood-Free Inference (SELFI) algorithm to infer the initial matter power spectrum, which we utilise to thoroughly investigate the impact of systematic effects. This investigation relies on a single set of N-body simulations. Second, we obtain a posterior on cosmological parameters via implicit likelihood inference, recycling the simulations from the first step for data compression. For demonstration, we rely on a model of large-scale spectroscopic galaxy surveys that incorporates fully non-linear gravitational evolution and simulates multiple systematic effects encountered in real surveys. We provide a practical guide on how the SELFI posterior can be used to assess the impact of misspecified galaxy bias parameters, selection functions, survey masks, inaccurate redshifts, and approximate gravity models on the inferred initial matter power spectrum. We show that a subtly misspecified model can lead to a bias exceeding $2\sigma$ in the $(\Omega_\mathrm{m},\sigma_8)$ plane, which we are able to detect and avoid prior to inferring the cosmological parameters. This framework has the potential to significantly enhance the robustness of physical information extraction from full-forward models of large-scale galaxy surveys such as DESI, Euclid, and LSST.

Millicharged particles are generic in theories of dark sectors. A cosmic or local abundance of them may be produced by the early universe, stellar environments, or the decay or annihilation of dark matter/dark energy. Furthermore, if such particles are light, these production channels result in a background of millicharged radiation. We show that light-shining-through-wall experiments employing superconducting RF cavities can also be used as ``direct deflection" experiments to search for this relativistic background. The millicharged plasma is first subjected to an oscillating electromagnetic field of a driven cavity, which causes charge separation in the form of charge and current perturbations. In turn, these perturbations can propagate outwards and resonantly excite electromagnetic fields in a well-shielded cavity placed nearby, enabling detection. We estimate that future versions of the existing Dark SRF experiment can probe orders of magnitude of currently unexplored parameter space, including millicharges produced from the Sun, the cosmic neutrino background, or other mechanisms that generate a thermal abundance with energy density as small as $\sim 10^{-4}$ that of the cosmic microwave background.

The blue loop stage of intermediate mass stars has been called a "magnifying glass", where even seemingly small effects in prior stages of evolution, as well as assumptions about stellar composition, rotation, and convection, produce discernible changes. As such, blue loops, and especially the existence and properties of Cepheids, can serve as a laboratory where feebly connected Beyond Standard Model particles such as axions can be gainfully studied. We undertake a careful study of the effects of these putative particles on the blue loop, paying close attention to the evolution of the core potential and the hydrogen profile. Our simulations, performed with MESA, place bounds on the axion-photon coupling using the galactic Cepheid S Mus, with dynamically-determined mass of $6 M_\odot$, as a benchmark. The effects of varying convective overshoot on the core potential and hydrogen profile, and the ensuing changes in the axion constraints, are carefully studied. Along the way, we explore the "mirror principle" induced by the hydrogen burning shell and contrast our results with those existing in the literature. Less conservative (but more stringent) bounds on the axion-photon coupling are given for a $9 M_\odot$ model, which is the heaviest that can be simulated if overshoot is incorporated, and tentative projections are given for a $12 M_\odot$ model, which is approximately the heaviest tail of the mass distribution of galactic Cepheids determined by pulsation models using Gaia DR2. Our main message is that the reliable simulation and observation (ideally, through dynamical mass determination) of massive Cepheids constitutes an important frontier in axion searches, challenges in modeling uncertainties in the microphysics of the blue loop stage notwithstanding.

Light feebly-coupled bosonic particles can efficiently extract the rotational energy of rapidly spinning black holes on sub-astrophysical timescales via a phenomenon known as black hole superradiance. In the case of light axions, the feeble self-interactions of these particles can lead to a non-linear coupled evolution of many superradiant quasi-bound states, dramatically altering the rate at which the black hole is spun down. In this work, we extend the study of axion superradiance to higher order states, solving for the first time the coupled evolution of all states with $n \leq 5$ in the fully relativistic limit (with $n$ being the principle quantum number). Using a Bayesian framework, we re-derive constraints on axions using the inferred spins of solar mass black holes, demonstrating that previously adopted limit-setting procedures have underestimated current sensitivity to the axion decay constant $f_a$ by around one order of magnitude, and that the inclusion to higher order states allows one to reasonably capture the evolution of typical high-spin black holes across a much wider range of parameter space, thereby allowing constraints to be extended to more massive axions. We conclude with an extensive discussion on the systematics associated with spin inference from x-ray observations.

Axions and other putative feebly interacting particles (FIPs) with a mass of tens to several hundreds of keVs can be produced in stellar cores with a Lorentz boost factor $E_a/m_a\lesssim 10$. Thus, starburst galaxies such as M82 are efficient factories of slow axions. Their decay $a\rightarrow\gamma\gamma$ would produce a large flux of X-ray photons, peaking around $100$ keV and spread around the galaxy by an angle that can be relatively large. We use observations of the Nuclear Spectroscopic Telescope Array (NuSTAR) mission to show that the absence of these features can constrain $30-500$ keV axion masses into uncharted regions for axion-photon coupling of $g_{a\gamma}\sim 10^{-10}-10^{-12}\,\rm GeV^{-1}$. Our argument can be applied to other heavy FIPs and astrophysical sources that are hot enough to produce them, yet cold enough to avoid large boost factors which slow down the decay.

We propose a first-order theory of dissipative fluids in the trace-fixed particle frame, which is similar to Eckart's frame except that the temperature is determined by fixing the trace of the stress-energy tensor. Our theory is hyperbolic and causal provided a single inequality holds. For low wave numbers, the expected damped modes in the shear, acoustic, and heat diffusion channels are recovered. Stability of global equilibria with respect to all wave numbers is also analyzed. The conditions for hyperbolicity, causality and stability are satisfied for a simple gas of hard spheres or disks.

In this paper, we propose a novel electron-assisted Baryogenesis scenario that does not require explicit B-L violation, which is essential for the traditional Leptogenesis mechanism. This scenario is based on the assumption of high-scale electroweak symmetry restoration, which implies that the electron Yukawa interaction, crucial for the mechanism, does not reach thermal equilibrium before the electroweak sphaleron process is quenched in the early universe. Primordial charge asymmetries for chiral electrons, which can be generated through various mechanisms such as axion inflation, the evaporation of primordial black holes, or the CP-asymmetric decays of a heavy Higgs doublet, serve as the initial condition for the amplification of the baryon asymmetry through transport equations. Right-handed electron asymmetry is almost irrelevant to the baryon asymmetry due to high-scale electroweak symmetry restoration, leading to both a non-zero baryon asymmetry and the electron asymmetry. We dub this mechanism as the Eogenesis.

Let two test particles A and B revolving about a spinning primary along ideally identical orbits in opposite directions be considered. From the general expressions of the precessions of the orbital inclination induced by the post-Newtonian gravitomagnetic and Newtonian quadrupolar fields of the central object, it turns out that the Lense-Thirring inclination rates of A and B are equal and opposite, while the Newtonian ones due to the primary's oblateness are identical. Thus, the difference of the inclination shifts of the two orbiters would allow, in principle, to cancel out the classical effects by enhancing the general relativistic ones. The conditions affecting the orbital configurations that must be satisfied for this to occur and possible observable consequences in the field of Earth are investigated. In particular, a scenario involving two spacecraft in polar orbits, branded POLAr RElativity Satellites (POLARES) and reminiscent of an earlier proposal by Van Patten and Everitt in the mid-1970s, is considered. A comparison with the ongoing experiment with the LAser GEOdynamics Satellite (LAGEOS) and LAser RElativity Satellite (LARES) 2 is made.

Interaction rates of neutrinos and antineutrinos within a QED plasma determine the dynamics of their decoupling in the early universe. We show how to define the relevant double-differential production, annihilation, and scattering rates at NLO. Integrating over these rates with specific weights, other quantities from the literature can be obtained, such as energy transfer rates, or a neutrino interaction rate. We show that NLO corrections to the energy transfer rates are as small as those that enter the previously determined neutrino interaction rate, and are therefore expected to have only a small influence on the neutrino decoupling parameter, $N_{\rm eff}\,$. We also provide a tabulation and fast interpolation routine for all double-differential rates, in order to allow for their use in non-approximate kinetic equations, which may further reduce the systematic uncertainties of the Standard Model prediction for $N_{\rm eff}\,$

It is notoriously difficult to construct a stable non-singular bouncing cosmology that avoids all possible instabilities throughout the entire evolution of the universe. In this work, we explore whether a non-singular bounce driven by a specific class of modifications of General Relativity, the vector-tensor generalized Proca theories, can be constructed without encountering any pathologies in linear perturbation theory. We find that such models unavoidably lead to either an instability in the matter sector or strong coupling in the scalar one. As our analysis is performed in a gauge-independent way, this result can be cast in the form of a no-go theorem for non-singular bounces with generalized Proca. In contrast to the no-go theorem found for Horndeski theories, however, it cannot be evaded by considering beyond generalized Proca theory. At the core of our result lies the non-dynamical nature of the temporal component of the vector field, which renders it an ill-suited mediator for a bouncing solution.

This study investigates the radial oscillations of hybrid neutron stars, characterized by a composition of hadronic external layers and a quark matter core. Utilizing a density-dependent relativistic mean-field model that incorporates hyperons and baryons for describing hadronic matter, and a density-dependent quark model for quark matter, we analyze the ten lowest eigenfrequencies and their corresponding oscillation functions. Our focus lies on neutron stars with equations-of-state involving N, N + $\Delta$, N + H, and N + H + $\Delta$, featuring a phase transition to quark matter. Emphasizing the effects of a slow phase transition at the hadron-quark interface, we observe that the maximum mass is attained before the fundamental mode's frequency decreases for slow phase transitions. This observation implies the stability of stellar configurations with higher central densities than the maximum mass, called Slow Stable Hybrid Stars (SSHSs), even under small radial perturbations. The length of these SSHS branch depends upon the energy density jump between two phases and the stiffness of the quark EoS.

Axions have long been considered plausible candidates for dark matter. The axion dark matter emitted from cosmic strings after the Peccei-Quinn (PQ) symmetry breaking in the early Universe was extensively simulated. In this work, we study dark matter and gravitational waves through the lattice simulation of the Axion-Higgs string. We gave the dark matter overproduction and the Big Bang nucleosynthesis bounds on the axion decay constant $f_a$ and the axion mass $m_a$ for axion-like particles, and found that the predicted gravitational wave spectra cannot be probed by the dataset of the current pulsar timing array experiments.