Papers for Thursday, Jun 26 2025

This thesis constitutes an observational effort to advance our understanding of planet formation and evolution across different regimes. To do so, we exploited Kepler/K2 and TESS data in combination with ground-based observations, serving as a preparation for the PLATO mission. This endeavour led to the discovery and statistical validation of four planets, the confirmation and characterisation of two other planets, the detection of 14 new planet candidates, and the revision of the properties of 25 previously reported planets. Among them, we find TOI-244 b to be a super-Earth with a density below what would be expected for an Earth-like composition, and we thus refer to TOI-244 b as a low-density super-Earth (LDSE). We also confirm the transiting super-Neptune TOI-5005 b. In this system, we find photometric variability that matches the planetary orbital period, suggesting the existence of magnetic star-planet interactions (MSPIs). We also detect similar synchronised signals in the eccentric HD 118203 system, which appear and disappear across different orbits, in agreement with the well-known 'on/off' nature of MSPIs. At a population level, we find that LDSEs tend to be hosted by stars with subsolar metallicities and that the lowest dense LDSEs tend to receive low insolation fluxes. We also identify an over-density of Neptunian planets in orbits corresponding to periods between ~3.2 and 5.7 days, which we refer to as the Neptunian 'ridge'. The 'ridge', coinciding with the orbital range of the hot Jupiter pileup, suggests a common path in the evolution of the closest giant planets. The low-insolation trend of LDSEs suggests that the largest fraction of volatile elements in these planets is directly exposed to the received irradiation from their host stars. Solving these two problems will greatly contribute to completing the full picture of the evolution of planetary systems.

We present a model that explains some extreme variability phenomena observed in active galactic nuclei (AGN). In this model, an orbiting companion interacts with the accretion disk surrounding the central supermassive black hole (SMBH). This interaction excites spiral density waves, leading to bursts of mass inflow lasting a few orbital timescales, whose mass content depends sensitively on the sphere of influence of the orbiting companion. To explain changing-look (CL) AGN, we find that lighter SMBH companions are necessary, while we generally exclude stellar companions. The moderately massive black hole perturber must be on a highly eccentric orbit in order to account for the non-repeating nature of most CL AGN. When applied to quasi-periodic eruptions (QPEs), we find that stars with highly eccentric orbits and close pericenter passages can produce accretion flares with QPE-like characteristics. For many QPEs, low-mass main-sequence stars or stripped-envelope stars are required. Although moderately massive black hole perturbers could also match the observed properties of QPEs, the gravitational-wave merger timescales of such binary systems are prohibitively short.

Clusters of galaxies host spectacular diffuse radio sources, extending over scales from 100 kpc to several Mpcs. These sources, with extremely faint surface brightness ($\mu$Jy/arcsec$^2$ level), are not tied to individual galaxies but trace synchrotron emission from large-scale magnetic fields and relativistic particles within the intracluster environment. Here, we report the discovery of a candidate radio mini-halo in SpARCS104922.6+564032.5, the most distant cool-core galaxy cluster identified to date at $z=1.709$, using deep LOFAR 120-168 MHz observations. We show that this emission originates from diffuse cluster-associated processes rather than unresolved AGN or star-forming galaxies. The diffuse radio emission coincides spatially with the X-ray emission of the hot intracluster medium and has a radio power of $P_{\rm 150~MHz}=49.8^{+14.7}_{-11.7} \times10^{24}$ W Hz$^{-1}$, exhibiting striking similarities to low-redshift radio mini-halos. This discovery doubles the redshift of previously known mini-halos, challenging models of inverse Compton losses and indicating the presence of strong magnetic fields, enhanced turbulence in high-redshift clusters, or active hadronic processes that require a cosmic ray-to-thermal energy ratio of 0.07 within 200 kpc, assuming a clumped distribution with spatial correlations among the gas, cosmic rays, and magnetic field that partially compensate for cosmological redshift dimming. It further implies that magnetic fields are efficiently amplified to $\sim$$10~\mu$G levels within a Mpc$^3$ volume during the epoch of cluster formation before $z\sim2$. These findings provide critical insights into high-redshift cluster physics and emphasize the transformative potential of next-generation radio surveys, such as those with the SKA and ngVLA, in exploring the early evolution of galaxy clusters.

Supermassive black holes (SMBHs) grow primarily through gas accretion, observed as active galactic nuclei (AGNs). While mergers can drive luminous AGN episodes, secular processes may fuel a substantial portion of cosmic BH growth. Whether these mechanisms can sustain high BH accretion rates remains uncertain. This study aims to identify the secular mechanism driving high BH accretion rates, by targeting a galaxy with moderately massive SMBH, high central gas densities, accretion rates of a few percent of the Eddington limit, and gas kinematics resolved close to black hole-dominated scales. A blind search led to the identification of NGC 4593, which is representative of the AGN population driving BH mass density growth since $z=1$. A prominent single-arm ("$m=1$") molecular gas spiral with ${\rm log}\,M_{\rm mol}/{\rm M}_\odot=8.1\pm0.3$ extends from 1.3$\,$kpc down to the SMBH's sphere of influence ($1.7^{+0.5}_{-0.2}\,{\rm pc}$). Star formation in the spiral is inefficient (${\rm SFR} = 4.9 \times 10^{-2} \,{\rm M}_\odot {\rm /yr}$, $\langle t_{\rm dep}\rangle=3.9 \pm 0.6\,\rm{Gyr}$), whereas inflow rate exceeds the SFR by two orders of magnitude and is sufficient to sustain the current BH accretion rate for $\geq$35$\,$Myr, enabling $\sim$10% SMBH growth. A young, lopsided stellar component (${\rm log}\,M_\star/{\rm M}_\odot=7.5-9.3$) exerts torques on the molecular gas, likely driving the gas inflow. This young stellar component may serve as both a cause and a product of sustained gas funnelling towards the SMBH. These findings provide evidence for the sustained secular $m=1$ feeding mode at high BH accretion rates, linking kpc-scale gas dynamics to the black hole's sphere of influence. This mechanism, consistent with simulation predictions, may represent a key contributor to SMBH growth in luminous AGNs since cosmic noon.

Modelling the gravitational interaction between an eccentric perturber and a differentially shearing gas disc is a longstanding problem with various astrophysical applications, ranging from the evolution of planetary systems to the migration of black holes embedded in AGN discs. Recent work has promoted the use of a global, linear, modal approach for calculating the excited wake and the resulting feedback on the perturber's orbital evolution. In this work we perform a complementary suite of targeted hydrodynamic simulations to test this linear framework across a range of disc temperature and density profiles. In particular, we push from circular orbits to highly eccentric trajectories for which the perturber moves supersonically with respect to the background gas. We find remarkable agreement between our simulations and the linear methodology across a range of diagnostics -- lending support to the predicted wake morphologies, complex radial torque density profiles, and torque reversal phenomena, which occur when the eccentricity exceeds the local scale height. In contrast, comparison with previous fitting functions exposes noticeable discrepancies, cautioning against their indiscriminate use in studies which explore a wide range of perturber eccentricities, in varied disc structures. Our simulations also probe the fundamentally nonlinear effects of shock induced angular momentum deposition and coorbital horseshoe drags, which exhibit clear variations with eccentricity. Finally, this careful comparison between linear theory and numerics provides a detailed benchmark for planet-disc interaction problems and therefore we have provided a repository of our linear calculations for use as a rigorous test of future numerical investigations.

Models of Stellar Population Synthesis (SPS) provide a predictive framework for the spectral energy distribution (SED) of a galaxy. SPS predictions can be computationally intensive, creating a bottleneck for attempts to infer the physical properties of large populations of individual galaxies from their SEDs and photometry; these computational challenges are especially daunting for near-future cosmology surveys that will measure the photometry of billions of galaxies. In this paper, we explore a range of computational techniques aimed at accelerating SPS predictions of galaxy photometry using the JAX library to target GPUs. We study a particularly advantageous approximation to the calculation of galaxy photometry that speeds up the computation by a factor of 50 relative to the exact calculation. We introduce a novel technique for incorporating burstiness into models of galaxy star formation history that captures very short-timescale fluctuations with negligible increase in computation time. We study the performance of Hamiltonian Monte Carlo (HMC) algorithms in which individual chains are parallelized across independent GPU threads, finding that our pipeline can carry out Bayesian inference at a rate of approximately $1,000$ galaxy posteriors per GPU-minute. Our results provide an update to standard benchmarks in the literature on the computational demands of SPS inference; our publicly available code enables previously-impractical Bayesian analyses of large galaxy samples, and includes several standalone modules that could be adopted to speedup existing SPS pipelines.

Solid state CO$_2$ has been detected throughout the outer solar system, even at temperatures where crystalline CO$_2$ is unstable, requiring that the CO$_2$ be trapped in a separate host material. The Saturnian satellites provide an ideal laboratory for the study of this trapped CO$_2$, allowing us to examine objects with identical insolation, but with a range of environments, ice exposure, organic abundance, and formation locations. Here, we present JWST spectra of 8 mid-sized satellites of Saturn, including Mimas, Enceladus, Tethys, Dione, and Rhea interior to Titan, and Hyperion, Iapetus, and Phoebe exterior. The $\sim$4.26 $\mu$m CO$_2$ $\nu_3$ band is detected on each satellite, and the $\sim$2.7 $\mu$m $\nu_1+\nu_3$ band is detected on all but Phoebe and the leading hemisphere of Iapetus. Based on the wavelength shifts of these bands, we find four separate types of trapped CO$_2$ on the satellites. On the inner satellites, CO$_2$ appears trapped in amorphous ice sourced from Saturn's E-ring, and a second component of CO$_2$ is associated with the dark material most prominent on the trailing hemispheres of Dione and Rhea. On the outer satellites, CO$_2$ appears to be produced by irradiation of organics on Phoebe, which are then transported to the dark leading hemisphere of Iapetus and onto the dark regions of Hyperion. CO$_2$ is also trapped by water ice on the trailing hemisphere of Iapetus and on Hyperion. These observations point to the continued need for laboratory studies to better understand the sources and trapping mechanisms of CO$_2$ throughout the outer solar system.

We present 0.8 to 5 $\mu$m JWST spectra of four $\sim$20 km diameter Jupiter Trojans known to have albedos elevated above the values typical in the remaining Trojan population. The spectra of these four high albedo Jupiter Trojans are all similar, with red slopes in the optical-IR transition region, a break to lower slopes at 1.3 $\mu$m, and a broad absorptions from 2.8 to 4 $\mu$m. The 0.8 to 2.5 $\mu$m spectra of these objects match the spectra of neither the well-known "red" and "less-red" Jupiter Trojans nor of any known asteroid taxonomic class. The reflecticity of these objects does not rise redward of 4 $\mu$m, a property that is seen in the previous JWST observations of Jupiter Trojans only in Polymele. Indeed, the high albedo Jupiter Trojan spectra are a good match to that of Polymele, and Polymele is both the smallest Jupiter Trojan in the previous JWST sample and has the highest albedo of the objects in that sample. We conclude that Polymele and the other high albedo Jupiter Trojans represent a third class of Jupiter Trojans not represented in the more heavily-studied larger objects and are perhaps the products of recent disruptions. The Lucy flyby of Polymele in September 2027 will give a direct view of one of this new class of Jupiter Trojans.

The 21cm line from the spin-flip transition of neutral hydrogen (HI) provides a unique window into the Epoch of Reionization (EoR), the final phase transition of our Universe. The Square Kilometre Array (SKA) enables precise measurements of 21cm fluctuations that trace ionization, temperature, and density fluctuations of the intergalactic medium (IGM). Nevertheless, a direct reconstruction of the timeline of the EoR in terms of the progress of ionization remains an ongoing challenge due to the highly non-Gaussian nature and thus intractable likelihood of the 21cm signal. Here, we present EoRFlow, a simulation-based inference (SBI) framework for reconstructing the global neutral hydrogen fraction $x_{\mathrm{HI}}(z)$ directly from 2D cylindrically averaged power spectra (2DPS) of the 21cm signal. We validate our method on realistic mock datasets for SKA-Low. Bypassing the need for explicit likelihood formulations, our approach enables fast, unbiased posterior estimation of the $x_{\mathrm{HI}}$ evolution in narrow redshift slices, allowing for piecewise reconstruction of the global reionization history. By directly inferring the reionization history from 21cm power spectra, our framework provides a scalable and robust path forward for 21cm cosmology in the SKA era.

We present a near-infrared spectroscopic analysis (0.9-2.4 micron) of gravity indices for 57 ultracool dwarfs (spectral types M5.5 to L0), including exoplanet hosts TRAPPIST-1, SPECULOOS-2, SPECULOOS-3, and LHS 3154. Our dataset includes 61 spectra from the SpeX and FIRE spectrographs. Using gravity-sensitive indices such as FeH absorption (at 0.99, 1.20, and 1.55 microns), the VO band at 1.06 microns, the H-band continuum, and alkali lines like K I (at 1.17 and 1.25 microns), we investigate correlations between surface gravity, stellar metallicity, and the presence of close-in transiting planets. All four planet-hosting stars show intermediate-gravity spectral signatures despite indicators of field age. However, a volume-corrected logistic regression reveals no significant association between gravity class and planet occurrence. Among individual indices, FeH_z is the most promising tracer of planet-hosting status. We tentatively identify a correlation between FeH_z (0.99 micron) and planet presence at the 2-sigma level, though this may reflect observational biases including transit probability, small-number statistics, and detection sensitivity. More robustly, we find a significant anti-correlation between FeH_z and metallicity ([Fe/H]) at 3.3 sigma. A Kruskal-Wallis test shows no significant metallicity difference across gravity classes, suggesting the observed FeH_z-metallicity trend is not driven by bulk metallicity differences. We propose this anti-correlation reflects interplay between age, gravity, and composition: higher-metallicity objects may be systematically younger with lower gravities, suppressing FeH absorption. While our results only hint at a link between gravity-related characteristics and planet occurrence among late-M dwarfs, they underscore the need for caution when using spectral diagnostics to infer properties of planet-hosting ultracool dwarfs.

The study of planetesimal debris accreted by white dwarfs offers unique insights into the composition of exoplanets. Using far-ultraviolet and optical spectroscopy, we have analysed the composition of planetesimals accreted by three metal enriched H-dominated white dwarfs with effective temperatures of T_eff = 20 000 K. WD 0059+257 is accreting an object composed of 71.8 +/- 7.9 per cent Fe and Ni by mass, indicating a large core mass fraction of 69 per cent, similar to that of Mercury. We model this planetesimal as having a differentiated Earth-like composition with 65 per cent of its mantle stripped, and we find this mass loss can be caused by vaporisation of the planetesimal's mantle during post-main sequence evolution. The tentative S detection in WD 0059+257 is a possible clue to the nature of the light element in planetary cores, including that of the Earth. The volatile-rich composition of WD 1943+163 is consistent with accretion of a carbonaceous chondrite-like object, but with an extreme Si depletion. WD 1953-715 accretes a planetesimal which contains 64 +/- 21 per cent of O in the form of ices, likely H2O. This body therefore requires an initial orbit at formation beyond a radial distance of > 100 au for ice survival into the white dwarf phase. These three planetary enriched white dwarfs provide evidence of differing core fractions, volatile budgets, and initial orbital separations of the accreted planetesimals, all of which help us understand their formation and evolutionary history.

The young M-star TWA 7 hosts a bright and near face-on debris disk, which has been imaged from the optical to the submillimeter. The disk displays multiple complex substructures such as three disk components, a large dust clump, and spiral arms, suggesting the presence of planets to actively sculpt these features. The evidence for planets in this disk was further strengthened with the recent detection of a point-source compatible with a Saturn-mass planet companion using JWST/MIRI at 11 $\mu$m, at the location a planet was predicted to reside based on the disk morphology. In this paper, we present new observations of the TWA 7 system with JWST/NIRCam in the F200W and F444W filters. The disk is detected at both wavelengths and presents many of the same substructures as previously imaged, although we do not robustly detect the southern spiral arm. Furthermore, we detect two faint potential companions in the F444W filter at the 2-3$\sigma$ level. While one of these companions needs further followup to determine its nature, the other one coincides with the location of the planet candidate imaged with MIRI, providing further evidence that this source is a sub-Jupiter mass planet companion rather than a background galaxy. Such discoveries make TWA 7 only the second system, after $\beta$ Pictoris, in which a planet predicted by the debris disk morphology has been detected.

We present the first high-precision proper motion catalog, tied to the International Celestial Reference System (ICRS), of infrared astrometric reference stars within R $\leq$ 25" (1 pc) of the central supermassive black hole at the Galactic center (GC). This catalog contains $\sim$2,900 sources in a highly extinguished region that is inaccessible via Gaia. New astrometric measurements are extracted from HST observations (14 epochs, 2010 - 2023) and transformed into the ICRS using 40 stars in common with Gaia-DR3. We implement a new method for modeling proper motions via Gaussian Processes that accounts for systematic errors, greatly improving measurement accuracy. Proper motion and position measurements reach precisions of $\sim$0.03 mas/yr and $\sim$0.11 mas, respectively, representing a factor of $\sim$20x improvement over previous ICRS proper motion catalogs in the region. These measurements define a novel HST-Gaia reference frame that is consistent with Gaia-CRF3 to within 0.025 mas/yr in proper motion and 0.044 mas in position, making it the first ICRS-based reference frame precise enough to probe the distribution of extended mass within the orbits of stars near SgrA*. In addition, HST-Gaia provides an independent test of the radio measurements of stellar masers that form the basis of current GC reference frames. We find that the HST-Gaia and radio measurements are consistent to within 0.041 mas/yr in proper motion and 0.54 mas in position at 99.7% confidence. Gaia-DR4 is expected to reduce the HST-Gaia reference frame uncertainties by another factor of $\sim$2x, further improving the reference frame for dynamical studies.

Context. Open clusters are groups of stars formed from the same cloud of gas and cosmic dust. They play an important role in studying stars' formation and evolution and understanding galaxies' structure and dynamics. Aims. The main objective of this work is to identify stars that belong to open clusters using astrometric data from Gaia EDR3 and spectroscopic data from APOGEE DR17. Furthermore, we investigate the metallicity gradients and orbital properties of the open clusters in our sample. Methods. By applying the HDBSCAN clustering algorithm to this data, we identified observed stars in our galaxy with similar dynamics, chemical compositions, and ages. The orbits of the open clusters were also calculated using the GravPot16 code. Results. We found 1987 stars belonging to 49 open clusters, and 941 of these stars have probabilities above 80 % of belonging to open clusters. Our metallicity gradient presents a two-slope shape for two measures of different Galactic center distances, the projected Galactocentric distance, and the guiding center radius to the Galactic Center, as already reported in previous work. However, when we separate the open clusters by age, we observe no significant difference in the metallicity gradient slope beyond a certain distance from the Galactic center. Our results show a shallower gradient for clusters younger than 2 Gyr than those older than 2 Gyr. All our OCs dynamically assemble the Disk-like population very well, and they are in prograde orbits, which is typical for disk-like populations. Some OCs resonate with the Galactic bar at the Lagrange points L4 and L5.

We present a machine learning (ML) framework for the detection of wide binary star systems and the identification of deviations from standard gravity using Gaia DR3 data. By training supervised ML models on established wide binary catalogues, we efficiently classify wide binaries and employ clustering and nearest neighbour search to pair candidate systems. We recast the problem of identifying systematic deviations from standard gravity as a supervised anomaly detection. Our approach incorporates data preprocessing techniques such as SMOTE, correlation analysis, and PCA, and achieves high accuracy and recall in both binary classification and anomaly detection tasks. The resulting publicly available code enables rapid, scalable, and customizable analysis of wide binaries and gravitational anomalies, complementing conventional analyses and providing a valuable resource for future astrophysical studies.

Rotation plays an important role in stellar evolution. However, the mechanisms behind the transport of angular momentum in stars at various stages of their evolution are not well understood. To improve our understanding of these processes, it is necessary to measure and validate the internal rotation profiles of stars across different stages of evolution and mass regimes. Our aim is to constrain the internal rotation profile of the 12Msun $\beta$ Cep pulsator HD 192575 from the observed pulsational multiplets and the asymmetries of their component frequencies. We update the forward asteroseismic modelling of HD 192575 based on new TESS observations. We invert for the rotation profile from the symmetric part of the splittings, and compute the multiplet asymmetries due to the Coriolis force and stellar deformation treated perturbatively. We compare the computed asymmetries with the observed asymmetries. Our new forward asteroseismic modelling is in agreement with previous results, but with increased uncertainties, partially due to increased frequency precision, requiring us to relax certain constraints. Ambiguity in the mode identification is the main source of the uncertainty, which also affects the inferred rotation profiles. Almost all acceptable rotation profiles occur in the regime below 0.4/d and favour weak radial differential rotation, with a ratio of core and envelope rotation below two. We find that the quality of the match between the observed and theoretically predicted mode asymmetries is strongly dependent on the mode identification and the internal structure of the star. Our results offer the first detailed rotation inversion for a $\beta$ Cep pulsator. They show that the rotation profile and the mode asymmetries provide a valuable tool to further constrain the evolutionary properties of HD 192575, and in particular the details of angular momentum transport in massive stars.

We use legacy data from the Hubble Space Telescope (HST) and the International Ultraviolet Explorer (IUE) to search for a spectral line from the spontaneous decay of axion-like particle (ALP) dark matter. The HST data consists of blank sky observations taken with the Faint Object Spectrograph in the 165-240~nm wavelength range, while the IUE data consists of observations of the Virgo Cluster obtained with the long- and short-wavelength spectrographs, covering 195-325~nm and 123-200~nm, respectively. We set a 95\% C.L. upper limit on the ALP-photon coupling $g_{a\gamma} < 10^{-11}~\mathrm{GeV}^{-1}$ across the full ALP mass range. Notably, we rule out values of $g_{a\gamma}$ above $1.6 \times 10^{-12}~\mathrm{GeV}^{-1}$ for ALP masses between 12.4 and 14.5\,eV, improving upon previous limits by an order of magnitude.

In exoplanet research, the focus is increasingly on identifying Earth analogs, planets similar in density and habitability potential. As the number of rocky exoplanets grows, parallel discussions have emerged on system architectures and Galactic environments that may support life, drawing comparisons to our own Earth. This has brought renewed attention to the concept of the Galactic Habitable Zone (GHZ) as a broader context for interpreting the diversity of planetary environments. This study is the first to use detailed chemical evolution models to investigate the impact of stellar migration, modeled through a parametric approach, on the GHZ. Our findings reveal that stellar migration significantly enhances the number of stars capable of hosting habitable planets in the outer Galactic regions, with an increase of up to a factor of five at 18 kpc relative to a baseline value of unity at 6 kpc. Furthermore, we explore a novel scenario where the presence of gas giant planets increases the probability for the formation of terrestrial ones. We find that this increased probability is higher in the inner Galactic disc, but is also mitigated by stellar migration. In particular, at the present time, the number of FGK stars hosting terrestrial planets with minimum habitability conditions in the ring centered at 4 kpc is approximately 1.4 times higher than in scenarios where gas giants are assumed to hinder the formation and evolution of Earth-like planets. Without stellar migration, this factor increases to 1.5. Even larger ratios are predicted for terrestrial planets orbiting retired A stars, reaching 2.8 in models with stellar migration and 3.3 in models without it.

The radius of a planet is a fundamental parameter that probes its composition and habitability. Precise radius measurements are typically derived from the fraction of starlight blocked when a planet transits its host star. The wide-field Transiting Exoplanet Survey Satellite (TESS) has discovered hundreds of new exoplanets, but its low angular resolution means that the light from a star hosting a transiting exoplanet can be blended with the light from background stars. If not fully corrected, this extra light can dilute the transit signal and result in a smaller measured planet radius. In a study of hundreds of TESS planet discoveries using deblended light curves from our validated methodology, we show that systematically incorrect planet radii are common in the literature: studies using various public TESS photometry pipelines have underestimated the planet radius by a weighted median of $6.1\% \pm 0.3\%$, leading to a $\sim20\%$ overestimation of planet density. The widespread presence of these biases in the literature has profoundly shaped-and potentially misrepresented-our understanding of the exoplanet population. Addressing these biases will refine the exoplanet mass-radius relation, reshape our understanding of exoplanet atmospheric and bulk composition, and potentially inform prevailing planet formation theories.

Interplanetary (IP) shock waves are greatly interesting, as they represent significant phenomena in near-Earth space and are direct drivers of geomagnetic and radiation storms. Moreover, various data and parameters are being explored for the identification and characterization of these waves. The spatial dimensions of shock waves vary significantly with the conditions and parameters of the propagating medium. For example, the radii of curvature of the shock wave fronts can vary by several hundred Earth radii or more in the inner heliosphere. In this study, we improved the semi-automated identification of shock waves by analyzing the solar wind and IP magnetic field (MF) parameters. More precisely, we analyzed the data recorded by the Proton Alpha Sensor of the Solar Wind Analyzer (SWA-PAS) and magnetometer (MAG) onboard the Solar Orbiter (SolO) mission. These data were collected and analyzed during the SolO journey around the Sun in 2023 at distances of 0.29...0.95 AU. Employing the developed algorithm, we identified over 40 IP-shock waves that occurred in 2023 using SWA-PAS and MAG. Additionally, we determined and presented a list of shock types and their basic parameters, kinetic and magnetohydrodynamic. The compression ratios, plasma beta ${\beta}_{us}$, angle between the shock normal and upstream MF ${\theta}{_B}{_n}$, Mach number, and others were among these parameters. Furthermore, we investigated the statistical distributions of the ${\theta}{_B}{_n}$ and ${\beta}_{us}$ parameters in the upstream region. Finally, the dependence of number of the identified shock waves as a function of the distance away from the Sun was explored; the number of shocks increased gradually with the increasing heliocentric distance.

Aims. We investigate the gas reservoirs, star formation (SF) properties, and environment of the ultra-diffuse galaxy GAMA526784 to understand its formation history, the efficiency of molecular gas conversion into stars, and the possible role of an interacting companion in shaping its morphology. Methods. We analyse low and high-resolution CO observations to constrain the molecular gas content, compare with HI data, and examine the SF efficiency of GAMA526784. The potential influence of a newly identified nearby dwarf galaxy is assessed using photometric and spatial information. Results. GAMA526784 exhibits a regular HI reservoir (M_HI/M* = 2.88) but only upper limits on its molecular gas mass (M_H2(5sigma)/M* < 0.23). The HI reservoir and CO non-detection can be explained by several mechanisms: (1) predominance of CO-dark H2, invisible to CO observations but contributing to SF; (2) a time delay in HI-to-H2 conversion following a recent interaction; or (3) elevated turbulence inhibiting gas collapse. An identified companion, found at a projected distance of 48 kpc, shows similar colours and lies in the direction of young star clusters in GAMA526784, indicating a possible association. We hypothesise this companion may have triggered the formation of star clusters in GAMA526784 through a high-velocity encounter. Conclusions. Our findings suggest GAMA526784 may have undergone an interaction that influenced its gas reservoirs and SF activity. The presence of a nearby companion agrees with predictions of a high-speed encounter, potentially offering a rare example of such an interaction in progress. We suggest this encounter may have shaped the system's recent evolution. Future observations, particularly targeting molecular gas tracers beyond CO and resolved HI maps, will be crucial in determining the extent of GAMA526784's cold gas reservoir and the nature of its recent star formation.

We present a preliminary study exploring whether the stellar orbital circularity of simulated galaxies, available from precomputed catalogs in the IllustrisTNG project, can be used as a proxy for broad morphological classification. We focus on the publicly available "Stellar Circularities, Angular Momenta, Axis Ratios" catalog, which enables a simple kinematic decomposition of the stellar component into disk and spheroid subsystems. By validating this approach against the more detailed five-component kinematic decomposition in TNG50, we confirm that the circularity-based disk fraction correlates most strongly with the thin disk, while the bulge fraction broadly represents the combined contribution of classical bulges and stellar halos. We then apply this decomposition to galaxies in the TNG100 simulation at redshift $z=0$ and identify a data-motivated threshold of $\mathrm{F_{disk}} = 0.25$ to distinguish early- and late-type galaxies. This threshold, lower than the commonly adopted value of $0.4$, better captures the diversity of disk-dominated systems and avoids excluding galaxies with moderately prominent disks. Additionally, we identify irregular or morphologically complex systems based on galaxies with low total disk and spheroid mass fractions. Using this classification, we recover a morphology-density relation that is broadly consistent with observations: late-type galaxies dominate in the field, while early-type galaxies are the most prevalent morphological type in clusters. Our results demonstrate that stellar circularity alone can serve as an accessible and computationally efficient morphological proxy. We also discuss the potential for this classification to support machine learning efforts as a baseline or training set for future morphological studies.

The dusty interstellar medium (ISM) of the Milky Way is distributed in a complex, cloudy structure. It is fundamental to the radiation balance within the Milky Way, provides a reaction surface to form complex molecules, and is the feedstock for future generations of stars and planets. The life cycle of interstellar dust is not completely understood, and neither are its structure nor composition. The abundance, composition, and structure of dust in the diffuse ISM can be determined by combining infrared, optical and ultraviolet spectroscopy. JWST enables measurement of the faint absorption of ISM dust grains against bright stars at kiloparsec distances across the infrared spectrum. Here we present an overview of the project `Webb Investigation of Silicates, Carbons, and Ices' (WISCI) along with interpretation of two targets, GSC 08152-02121 and CPD-59 5831. Observations of 12 WISCI target stars were taken by JWST, the Hubble Space Telescope, Himalayan Chandra Telescope, and the Very Large Telescope. We use these to characterize the targets' spectral types and calculate their line-of-sight extinction parameters, $A_{\rm V}$ and $R_{\rm V}$. We find absorption in the JWST spectra of GSC 08152-02121, and CPD-59 5831 associated with carbonaceous dust around 3.4 and 6.2 micron and amorphous silicates at 9.7 micron. In GSC 08152-02121 we also find indications of absorption by trapped water around 3 micron. This first look from WISCI demonstrates the line-of-sight variability within the sample, and the program's potential to identify and correlate features across ultraviolet to mid-infrared wavelengths.

The reliable detection of the global 21-cm signal, a key tracer of Cosmic Dawn and the Epoch of Reionization, requires meticulous data modelling and robust statistical frameworks for model validation and comparison. In Paper I of this series, we presented the Beam-Factor-based Chromaticity Correction (BFCC) model for spectrometer data processed using BFCC to suppress instrumentally induced spectral structure. We demonstrated that the BFCC model, with complexity calibrated by Bayes factor-based model comparison (BFBMC), enables unbiased recovery of a 21-cm signal consistent with the one reported by EDGES from simulated data. Here, we extend the evaluation of the BFCC model to lower amplitude 21-cm signal scenarios where deriving reliable conclusions about a model's capacity to recover unbiased 21-cm signal estimates using BFBMC is more challenging. Using realistic simulations of chromaticity-corrected EDGES-low spectrometer data, we evaluate three signal amplitude regimes -- null, moderate, and high. We then conduct a Bayesian comparison between the BFCC model and three alternative models previously applied to 21-cm signal estimation from EDGES data. To mitigate biases introduced by systematics in the 21-cm signal model fit, we incorporate the Bayesian Null-Test-Evidence-Ratio (BaNTER) validation framework and implement a Bayesian inference workflow based on posterior odds of the validated models. We demonstrate that, unlike BFBMC alone, this approach consistently recovers 21-cm signal estimates that align with the true signal across all amplitude regimes, advancing robust global 21-cm signal detection methodologies.

Follow-up observations of neutrino events have been a promising method for identifying sources of very-high-energy cosmic rays. Neutrinos are unambiguous tracers of hadronic interactions and cosmic rays. On June 15, 2020, IceCube detected a neutrino event with an 82.8% probability of being astrophysical in origin. To identify the astrophysical source of the neutrino, we used X-ray tiling observations to identify potential counterpart sources. We performed additional multiwavelength follow-up with NuSTAR and the VLA in order to construct a broadband spectral energy distribution (SED) of the most likely counterpart. From the SED, we calculate an estimate for the neutrinos we expect to detect from the source. While the source does not have a high predicted neutrino flux, it is still a plausible neutrino emitter. It is important to note that the other bright X-ray candidate sources consistent with the neutrino event are also radio-quiet AGN. A statistical analysis shows that 1RXS J093117.6+033146 is the most likely counterpart (87.5%) if the neutrino is cosmic in origin and if it is among X-ray detectable sources. This results adds to previous results suggesting a connection between radio-quiet AGN and IceCube neutrino events.

The Zwicky Transient Facility Census of the Local Universe survey yielded a sample of 330 Type IIP supernovae (SNe) with well-constrained peak luminosities. In paper I (arXiv:2502.19493), we measured their luminosity function and volumetric rate. Here (paper II), we present the largest systematic study of lightcurve properties for Type IIP SNe from a volume-limited survey, analyzing a selected subset of 129 events, including 16 low-luminosity Type IIP (LLIIP) SNe with M${r,peak} \geq -16$ mag. We find that plateau slope correlates with peak brightness, with many LLIIP SNe showing positive slopes--suggesting smaller progenitor radii and distinct density profiles compared to brighter Type IIP SNe. The plateau duration shows only a weak dependence on peak brightness, likely suggesting binary interaction. One SN exhibits a plateau-to-tail drop of >3.5 mag, consistent with an electron-capture or failed SN with very low or zero nickel mass. We derive explosion and progenitor parameters of the entire Type IIP SN sample using semi-analytical and radiation-hydrodynamical models. Based on radiation-hydrodynamical model fitting, LLIIP SNe are characterized by low nickel masses (0.001-0.025 $\mathrm{M_\odot}$), low explosion energies (0.1-0.28 $\times 10^{51}$ erg), low ejecta masses ($8.1^{+0.8}_{-1.7}$ $\mathrm{M\odot}$), and ZAMS masses below 11 $\mathrm{M_\odot}$. In comparison, the full Type IIP SN sample spans a wider range with nickel masses (0.001-0.222 $\mathrm{M_\odot}$), explosion energies (0.10-4.43 $\times 10^{51}$ erg), ejecta masses (5.4-24.8 $\mathrm{M_\odot}$), and ZAMS masses (9.3-16.7 $\mathrm{M_\odot}$). We find strong correlations between peak brightness, explosion energy, and nickel mass that extend to the low-luminosity end. We conclude that LLIIP SNe represent the faint, low-energy end of the Type IIP population and originate from the lowest-mass core-collapse progenitors.

Red clump (RC) stars are reliable standard candles for studying the structure and evolution of the Milky Way. In this study, we present empirical calibrations of RC absolute magnitudes in the Mephisto (v, g, r, i) and CSST (g, r, i) photometric systems using a high-purity sample of 25,059 RC stars cross-matched between APOGEE and Gaia DR3 XP spectra. Through synthetic photometry and polynomial fitting, we find that RC absolute magnitudes exhibit strong dependencies on effective temperature and metallicity, with the strongest variations observed in bluer bands and progressively decreasing towards redder wavelengths. In particular, the Mephisto v band exhibits the highest sensitivity, with variations reaching up to 2.0 mag across the metallicity range (-1.0 dex <[Fe/H] <0.5 dex) and the temperature range (4500-5200 K). The calibrations achieve high precision for all bands, enabling accurate determination of RC absolute magnitudes and distances. Furthermore, we evaluate the metallicity estimation capabilities of both systems using a Random Forest-based method, achieving a precision of 0.12 dex for Mephisto and 0.14 dex for CSST under typical photometric uncertainties (< 0.01 mag). These results provide robust tools for distance and metallicity determinations, supporting future Galactic structure studies with Mephisto and CSST data.

The Galactic center region was observed with the XRISM X-ray observatory during the performance verification phase in 2024 and a point-like X-ray source was detected with the X-ray imager Xtend at a position of (RA, Dec)=(17h46m10.8s, -29°00'21''), which is thus named XRISM J174610.8-290021. This source was bright in February to March and showed time variations in count rate by more than one order of magnitude in one week. The 2-10 keV X-ray luminosity was ~$10^{35}$ erg/s for the assumed distance of 8 kpc. However, after six months, it was below the detection limit. We found a hint of periodicity of 1537 s from timing analysis. The XRISM/Xtend spectrum has emission lines from helium-like iron (Fe He$\alpha$) at 6.7 keV and hydrogen-like iron (Fe Ly$\alpha$) at 6.97 keV; their intensity ratio is unusual with the latter being four times stronger than the former. If the emission is of thermal origin, the ionization temperature estimated from the iron-line intensity ratio is ~30 keV, which is inconsistent with the electron temperature estimated from the thermal bremsstrahlung, ~7 keV. Spectral models of magnetic cataclysmic variables, which are often seen in the Galactic center in this luminosity range, are found to fail to reproduce the obtained spectrum. By contrast, we found that the spectrum is well reproduced with the models of low-mass X-ray binaries containing a neutron star plus two narrow Gaussian lines. We consider that the source is intrinsically bright reaching $10^{37}$ erg/s, but is blocked from direct view due to a high inclination and only the scattered emission is visible. The photoionized plasma above the accretion disk with an ionization parameter of ~$10^{5}$ may explain the unusual iron line ratio. We further discuss the potential contribution of point sources of the type of XRISM J174610.8-290021 to the diffuse Galactic center X-ray emission.

The origin of Galactic cosmic rays (CRs), particularly around the knee region ($\sim$3 PeV), remains a major unsolved question. Recent observations by LHAAASO suggest that the knee is shaped mainly by protons, with a transition to heavier elements at higher energies. Microquasars -- compact jet-emitting sources -- have emerged as possible PeV CR accelerators, especially after detections of ultrahigh-energy gamma rays from these systems. We propose that the observed proton spectrum (hard below a few PeV, steep beyond) arises from the reacceleration of sub-TeV Galactic CRs via shear acceleration in large-scale microquasar jet-cocoon structures. Our model also naturally explains the observed spectrum of energies around a few tens of PeV by summing up heavier nuclei contributions. Additionally, similar reacceleration processes in radio galaxies can contribute to ultrahigh-energy CRs, bridging Galactic and extragalactic origins. Combined with low-energy CRs from supernova remnants and galaxy clusters around the second knee region, this scenario could provide a unified explanation for CRs across the entire energy spectrum.

Hybrid morphology radio sources (HyMoRSs) are a rare subclass of radio galaxies that display a Fanaroff-Riley type I (FR I) morphology on one side of the central supermassive black hole (SMBH) and a type II (FR II) morphology on the other. These unique sources provide valuable insights into the physical mechanisms behind the FR dichotomy. In this study, we report the discovery of thirty-six new HyMoRSs, marking the largest collection of such sources in the southern sky to date, using data obtained from the MeerKAT absorption line survey (MALS). The identified HyMoRSs exhibit moderate radio luminosities, with a median value of 4.4 $\times 10^{24}$ W Hz$^{-1}$ at 1.4 GHz, and are located within the redshift range 0.04 < z < 1.34. To gain further insights into their origin, we also investigate the mid-infrared properties and environments of their host galaxies. Notably, nine out of the thirty-six sources are situated near the centers of galaxy clusters, including one with giant radio jets that extend over 811 kpc. Our analysis reveals that the majority of HyMoRSs are high-excitation radio galaxies (HERGs), hosted by actively star-forming galaxies that exhibit elevated star formation rates. Furthermore, our findings suggest that hybrid morphology may arise from FR II jets being deflected by a dense, cluster-like environment, along with orientation effects that make one jet appear FR I-like. As our candidates are selected through visual inspection of MALS radio maps, higher-resolution follow-up observations are still necessary to confirm the nature of their morphologies.

We present the most extensive high-time resolution multi-band campaign to date on the candidate transitional millisecond pulsar (tMSP) 3FGL J1544.6-1125 in the sub-luminous disk state, with coordinated observations from the radio to the X-ray band. While XMM-Newton and NuSTAR X-ray light curves exhibit the characteristic high- and low-mode bimodality, the source faintness prevents firm evidence for similar bimodality in the ultraviolet and near-infrared light curves, presented here for the first time. A re-analysis of archival XMM-Newton/OM data reveals an optical flare without an X-ray counterpart, likely originating from the outer accretion disk or the companion star. During our observations, no radio emission was detected, with a 3$\sigma$ flux density upper limit of 8 $\mu$Jy at 6 GHz. While past works have already reported radio variability in the source, this limit is a factor of 3.5 below the average value measured in 2019 in similar conditions, underscoring significant radio variability despite the relatively stable X-ray flux. Simultaneous optical light curves in five filters with GTC/HiPERCAM revealed flickering and dipping activities that resemble the observed X-ray variability, along with a reddening trend at lower fluxes. The latter is consistent with discrete mass ejections that disrupt the inner flow and reduce both X-ray and optical fluxes, thereby driving the high-to-low-mode switches. This suggests a common origin for most optical and X-ray emission at the boundary region between the pulsar wind and the inner disk, as also supported by our modelling of the spectral energy distribution in the high mode. Overall, our findings reinforce the mini-pulsar nebula picture for tMSPs in the sub-luminous state and demonstrate how coordinated, high-time resolution, multi-wavelength campaigns are essential to probe the processes governing rapid mode switches in these systems.

Molecular gas plays a critical role in regulating star formation and nuclear activity in galaxies. Sulphur bearing molecules, such as H2S, are sensitive to the physical and chemical environments in which they reside and are potential tracers of shocked, dense gas in galactic outflows and active galactic nuclei (AGN). We aim to investigate the origin of H2S emission and its relation to dense gas and outflow activity in the central regions of nearby infrared luminous galaxies. We present ALMA Band 5 observations of the ortho H2S 1(1,0) 1(0,1) transition in three nearby galaxies: NGC 1377, NGC 4418, and NGC 1266. We perform radiative transfer modelling using RADEX to constrain the physical conditions of the H2S emitting gas and compare the results to ancillary CO and continuum data. We detect compact H2S emission in all three galaxies, arising from regions smaller than approximately 150 parsecs. The H2S spectral profiles exhibit broad line wings, suggesting an association with outflowing or shocked gas. In NGC 4418, H2S also appears to be tracing gas that is counterrotating. A peculiar red shifted emission feature may correspond to inflowing gas, or possibly a slanted outflow. RADEX modelling indicates that the H2S emitting gas has high densities (molecular hydrogen density greater than 10^7 cm^-3) and moderately warm temperatures (between 40 and 200 Kelvin). The derived densities exceed those inferred from CO observations, implying that H2S traces denser regions of the interstellar medium.

We have surveyed all conventional methods proposed or conceivable for obtaining resolved images of an Earth-like exoplanet. Generating a 10x10 pixel map of a 1 R_Earth world at 10 pc demands ~0.85 uas angular resolution and photon collection sufficient for SNR >= 5 per "micro-pixel". We derived diffraction-limit and photon-budget requirements for: (1) large single-aperture space telescopes with internal coronagraphs; (2) external starshades; (3) space-based interferometry (nulling and non-nulling); (4) ground-based extremely large telescopes with extreme adaptive optics; (5) pupil-densified "hypertelescopes"; (6) indirect reconstructions (rotational light-curve inversion, eclipse mapping, intensity interferometry); and (7) diffraction occultation by Solar System bodies. Even though these approaches serve their primary goals - exoplanet discovery and initial coarse characterization - each remains orders of magnitude away from delivering a spatially resolved image. In every case, technology readiness falls short, and fundamental barriers leave them 2-5 orders of magnitude below the astrometric resolution and photon-budget thresholds needed to map an Earth analog even on decadal timescales. Ultimately, an in situ platform delivered to <= 0.1 AU of the target could, in principle, overcome both diffraction and photon-starvation limits - but such a mission far exceeds current propulsion, autonomy, and communications capabilities. By contrast, the Solar Gravitational Lens - providing on-axis gain of ~1e10 and inherent uas-scale focusing once a spacecraft reaches >=550 AU - remains the only near-term, scientifically and technologically viable means to acquire true, resolved surface images and spatially resolved spectroscopy of Earth-like exoplanets in our stellar neighborhood.

Using data from the MIRIS Pa{\alpha} Galactic Plane Survey (MIPAPS), we present a Pa{\alpha} 1.87 um line image of the entire Galactic plane within the latitude range of -3° < b < +3°, revealing numerous Pa{\alpha} features. Based on the MIPAPS Pa{\alpha} image and the WISE H II region catalog, we compile a catalog of 1489 Pa{\alpha} emission-line sources in the Galactic plane within 90° < l < 330°. By comparing our Pa{\alpha} images with H{\alpha} images constructed from the IPHAS and VPHAS+ survey data, we demonstrate the advantages of Pa{\alpha} line observations. We identify 902 Pa{\alpha} sources associated with H II regions, and newly confirm 619 H II region candidates as definitive H II regions through Pa{\alpha} or H{\alpha} detections. We also identify 261 extended and 326 point-like Pa{\alpha} sources not included in the WISE catalog, most of which have H{\alpha} counterparts in the IPHAS or VPHAS+ images. A search of the SIMBAD database indicates that these sources originate from diverse object types. By measuring Pa{\alpha} and H{\alpha} fluxes, we estimate the E(B-V) color excesses derived from extended emissions for 138 Pa{\alpha} sources, showing good agreement with values obtained from spectrophotometry of ionizing stars in previous studies. Futhermore, we calculate total Lyman continuum luminosities for 42 Pa{\alpha} sources, providing constraints on the distances to H II regions and the spectral types of their ionizing stars. These results highlight the scientific potential of Pa{\alpha} line observations and the benefits of combining multiple hydrogen recombination lines in exploring ionized regions.

The evolution of planetary systems around white dwarfs is crucial to understanding the presence of planetary material in the atmospheres of white dwarfs. These systems uniquely probe exoplanetary compositions. Periodic signals in the photometry of a handful of white dwarfs suggest material blocking the star, potentially from disintegrating planetesimals. Tidal evolution followed by scattering can bring planetesimals onto close-in orbits that would have been within the envelope of the white dwarf progenitor. The orbital period distribution of planetesimals undergoing tidal evolution will peak at short-period (nearly) circularized orbits (~ 10 hour-1 day), with a rising tail towards long-period highly eccentric orbits (~ 100 day). This prediction is generally consistent with the observed white dwarf transiting systems. In order for the planetesimal on the 4.5 hour period around WD 1145+017 to be explained by the tidal evolution of a planetesimal, that planetesimal must have an ultimate tensile strength comparable to that of iron meteorites.

Planetary material accreted by white dwarfs provides unique insights regarding exoplanetary composition. The evolutionary pathways of planetary bodies around white dwarfs are crucial to understanding the presence of close-in planetary material, observed in the form of pollutants in the atmospheres of white dwarfs and planetary material transiting white dwarfs. Periodic transits around white dwarfs potentially reveal the existence of close-in planetary bodies undergoing dust production. Tidal interactions can bring planetesimals that have been gravitationally perturbed onto long-period highly eccentric orbits around white dwarfs towards shorter orbital periods and smaller eccentricities. Tidal interactions may also induce melting and volcanism in these planetesimals, potentially being a mechanism for dust and debris production, the result of which may be seen in transit. Tidally induced volcanism may be triggered in a wide parameter space: for a 100km-sized rocky planetesimals perturbed to a pericentre distance $\lesssim$ 0.01AU ($\gtrsim$ twice its Roche limit), both on long-period (~ 100day) highly eccentric orbits and short-period (~ 10hr) near circular orbits. We comment on the potential link between the resultant volcanic ejecta and observed optical transits.

Between 1987 and 2023 we carried out single-dish monitoring campaigns of water-maser emission at 22 GHz in the circumstellar envelopes (CSEs) of four stars: RCas, oCet, RLeo and ChiCyg. The variability in integrated flux in RCas and oCet followed that in the optical with the same period, but with a lag of about one-third in phase. RLeo was too often below our sensitivity threshold for us to determine a radio period. No maser was detected in ChiCyg. The variability in the masers in RCas has a distinctive pattern that confirms the existence of a zone in the CSE with favourable conditions for maser excitation. The total flux, modulated by the pulsations of the star, gradually increases to a maximum, which is followed by a similar decrease. This takes about 20 years. The pattern is repeated after an interval of quiescence of several years. During its decline, the variation in the flux resembles a damped harmonic oscillator. There are two dominant emission components that move almost tangentially on either side of the star with respect to the observer. The redshifted one is consistent with an origin in a single cloud that lived for at least about eight years, falling back towards the star. The blue component, moving in the CSE hemisphere nearest to us, has no drift in the line-of-sight velocity and appears to originate in a time series of short-living clouds with line-of-sight velocities within ca. 1 km/s of each other. No systematic drifts are found in RLeo and oCet. A few bursts of emission were detected in RCas and RLeo that lasted about a year and caused an increase in the flux density by 1-2 orders of magnitude. The velocity range of the maser emission is < 10 km/s and is clearly correlated with the bolometric luminosity of the stars.

The study of the evolution of Martian atmosphere and its response to EUV irradiation is an extremely important topic in planetary science. One of the dominant effects of atmospheric losses is the photochemical escape of atomic oxygen from Mars. Increasing the magnitude of the irradiation changes the response of the atmosphere. The purpose of the current paper is to analyze the effects of enhanced EUV irradiation on the escape rates of oxygen atoms. We have used the solar flare of 2017 September 10 as the baseline flare intensity and varied the intensity of the flare from a factor of 3 up to 10 times the baseline flare. We see an increase in the escape flux by 40% for flares up to 5x the intensity of the baseline flare. However, beyond this point, the increase in escape flux tapers off, reaching only about 17% above the baseline. At 10x the baseline flare intensity, the escape flux decreases by nearly 25% compared to the escape rate of the original flare. We also found that the total escape amount of hot O peaks at 7x the original flare intensity. Additionally, we have studied the effects of the time scales over which the flare energy is delivered. We find that energy dissipative processes like radiative cooling and thermal collisions do not come into play instantaneously. The escape flux from higher intensity flares dominate initially, but as time progresses, energy dissipative processes have a significant effect on the escape rate.

The aim of this work is to obtain precise atmospheric parameters and chemical abundances automatically for solar twins and analogs to find signatures of exoplanets, as well as to assess how peculiar the Sun is compared to these stars and to analyze any possible fine structures in the Galactic thin disk. We developed a neural network (NN) algorithm using Python to obtain these parameters for a sample of 99 solar twins and solar analogs previously studied in the literature from normalized high-quality spectra from HARPS, with a resolving power of R $\sim$ 115000 and a signal-to-noise ratio S/N > 400. We obtained precise atmospheric parameters and abundance ratios [X/Fe] of 20 chemical elements (Li, C, O, Na, Mg, Al, Si, S, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, and Ba). The results are in line with the literature, with average differences and standard deviations of $(2 \pm 27)$ K for T$_{\rm eff}$, $(0.00 \pm 0.06)$ dex for log g, $(0.00 \pm 0.02)$ dex for [Fe/H], $(-0.01 \pm 0.05)$ km s$^{-1}$ for microturbulence velocity, $(0.02 \pm 0.08)$ km s$^{-1}$ for the macro turbulence velocity, and $(-0.12 \pm 0.26)$ km s$^{-1}$ for the projected rotational velocity (vsin$i$). Regarding the chemical abundances, most of the elements agree with the literature within 0.01 - 0.02 dex. The abundances were corrected from the effects of the Galactic chemical evolution and analyzed with the condensation temperature (T$_{\rm cond}$) to verify whether the stars presented depletion of refractories compared to volatiles. We found that the Sun is more depleted in refractory elements compared to volatiles than 89% of the studied solar analogs, with a significance of 9.5$\sigma$ when compared to the stars without detected exoplanets. We also found the possible presence of three subpopulations in the solar analogs: one Cu-rich, one Cu-poor, and the last one slightly older and poor in Na.

Context: Intermediate-mass asymptotic giant branch (AGB) stars influence Galactic lithium evolution by ejecting surface material (including Li) via stellar winds. Internal gravity waves (IGW), excited by convective motions, drive turbulent mixing in non-convective regions, altering stellar surface and wind chemistry. Aims: Investigate IGW-induced extra-mixing in the radiative zone between thermal pulses and convective envelopes of AGB stars and its impact on Li production. Derive the total Li contribution from intermediate-mass AGB stars using stellar models and initial mass functions. Methods: Construct stellar models (from zero-age main sequence to AGB end) with MESA, incorporating IGW-induced mixing and element diffusion. Calculate Li yields for stars of varying masses and metallicities using grids and population synthesis. Results: IGW triggers extra-mixing during He-shell flashes, transporting Be-7 from the radiative zone to the convective envelope, where it decays into Li-7. The positive effect of IGW on Li yield decreases with initial stellar mass but increases with metallicity. Most AGB stars (3.5-7.5 solar masses) produce positive Li yields. The total Li yield with IGW mixing (approximately 15 solar masses) is twice that without, contributing about 10 percent to Galactic Li. Conclusions: Through this extra-mixing mechanism induced by IGW, AGB stars can achieve a maximum A(Li) exceeding 5 and intermediate-mass AGB stars significantly contribute to Li in the Galactic ISM. These findings underscore the crucial role of IGW in stellar evolution, particularly in enhancing Li production.

The E$_{\rm p,i}$--L$_{\rm iso}$ correlation of long gamma-ray bursts (LGRBs) is regarded as a fundamental correlation for standardizing LGRBs to probe cosmology and constrain LGRB physics. However, this correlation may be affected by potential selection effects which are likely overlooked in the current small LGRB redshift sample. In this work, we simulate a large LGRB sample that reflects the observed situation, aiming to study the impact of peak flux $P$ on the observed LGRB E$_{\rm p,i}$--L$_{\rm iso}$ correlation. We find that the overall (E$_{\rm p,i}$, L$_{\rm iso}$) distribution, which will directly affect the best-fit result of the correlation, is significantly dependent on the value of $P$. This indicates that the impact of peak flux selection should be carefully considered in the studies and applications of the E$_{\rm p,i}$--L$_{\rm iso}$ correlation. Notably, we show that our simulated data can reproduce the observed $P$ distribution only if some dependence of (E$_{\rm p,i}$, L$_{\rm iso}$) is included in the simulation. This is an indication that the (E$_{\rm p,i}$, L$_{\rm iso}$) connection is a crucial property of LGRBs. We also find that GRBs with high peak flux in the low-E$_{\rm p,i}$ & L$_{\rm iso}$ region are not the straightforward extrapolation of the GRB population in the higher-E$_{\rm p,i}$ & L$_{\rm iso}$ region. Selecting four bursts with $ L_{\rm iso}\le10^{50}$ erg s$^{-1}$, $E_{\rm p, i}\le10^{2.5}$ keV, and $P\ge10^{0.5}$ ph cm$^{-2}$ s$^{-1}$, we find two bursts, GRB 060614 and GRB 191019A, which may not be associated with the theoretical massive-star origin of LGRBs. This suggests that combining $P$ with the position in the E$_{\rm p,i}$--L$_{\rm iso}$ diagram may be used to indicate alternative LGRB origins.

PBjam is an open-source software package for measuring mode frequencies of solar-like oscillators. These frequencies help constrain stellar evolution models to precisely estimate masses, radii, and ages of stars. The overall aim of PBjam is to simplify this process to the point where it may be done by non-experts or performed on thousands of stars with minimal interaction. The initial release of PBjam was restricted to only identifying modes of $\ell=0$ and $\ell=2$, since these are the simplest to treat consistently across different stellar evolutionary stages. Here we introduce a new set of three separate models which lets PBjam automatically identify $\ell=1$ modes in stars that experience varying degrees of coupling between p- and g-modes. These include a simple asymptotic relation for p-modes which can be applied to main-sequence stars, a matrix formalism aimed at treating frequency dependent coupling in sub-giants, and a uniform coupling model which is suitable for red giants. These models follow the Bayesian methodology established in the first release of PBjam, where a large set of previous observations is used to construct a nonparametric prior probability density for the new set of model parameters. This extension allows PBjam to build a more complete description of the power due to oscillations across a wider range of evolutionary stages.

A binary system (G3425) consisting of a massive unseen component and a red giant star on a nearly circular orbit was recently discovered. The formation of such a system is puzzling because orbital stability generally breaks down due to the large mass loss from the system caused by the SN explosion while forming the unseen component. Analytical solutions of the variable-mass two-body problem suggest that the explosion should have occurred when the component was close to its apocenter to explain the near-circular remnant system. This provides a strong constraint on the total mass and orbital configuration of the progenitor system. The nearly circular orbit of G3425 rules out type II SN scenarios and allows only for a fine-tuned SN~Ib/c explosion to occur when the secondary was close to its apocenter. Such a scenario, although possible, is highly unlikely. However, the most likely scenario is a failed SN that produced a black hole, for which no additional constraints on the position of the secondary are needed. We propose that the unseen component of G3425 is a mass-gap black hole with a mass constrained between the theoretical minimum for failed supernova progenitors (4 MSun) and the observed upper limit (4.4 MSun). Our analysis can be applied to any wide binary system containing an unseen component on a nearly circular orbit.

The lightest Kaluza-Klein (KK) dark matter particles and primordial black holes (PBHs) emerge as plausible candidates for dark matter. In scenarios where dark matter is a mix of KK particles and PBHs, PBHs can attract surrounding KK dark matter particles post-formation, leading to the creation of ultracompact dark matter halos (UCMHs). The distribution of KK dark matter particles within UCMHs tends to be steeper than in classical dark matter halo structures, such as the Navarro-Frenk-White model. Consequently, the annihilation rate of KK dark matter particles in UCMHs is significant. The high-energy photons resulting from the annihilation of KK particles in UCMHs contribute to the extragalactic gamma-ray background (EGB). Leveraging data from the $\mathtt{Fermi\text{-}LAT}$ experiment, we have derived, for the first time, upper limits on the cosmological abundance of PBHs in the context of KK dark matter annihilation. For a KK dark matter mass range of $500\le m_{\rm B^{(1)}}\le 1500$ GeV, which aligns with the observed present abundance of dark matter, the conservative limits on the fraction of dark matter in PBHs, for massive PBHs with $M_{\rm PBH}\gtrsim 10^{-11}M_{\odot}$, are $f_{\rm PBH} \lesssim 2\times 10^{-5}$.

Lyman-$ \alpha $ emission, which owing to its resonant nature strongly couples the emergent line profile to gas kinematics, has become a key observable for probing outflows from star-forming galaxies in the early universe. Typically, Lyman-$ \alpha $ profiles are explained in terms of idealised expanding shell models, while empirically-motivated interpretations of Lyman-$ \alpha $ profiles that make use of comparisons with kinematics revealed by other spectral features (e.g. metal emission and absorption lines) are limited. Here, we analyse the kinematics of emission/absorption lines in the spectra of 338 gravitationally-lensed Lyman-$ \alpha $ sources revealed via blind spectroscopy of galaxy clusters with the Multi-Unit Spectroscopic Explorer (MUSE). Using metal emission lines to measure systemic redshifts, we confirm that the Lyman-$ \alpha $ profiles are consistent with outflowing gas. In cases where interstellar metal absorption lines are detected, we find their kinematics and line ratios to be indicative of clumpy outflows. We investigate the differences between sources with single- and double-peaked Lyman-$ \alpha $ profiles, finding the latter to drive weaker outflows yet show spectral signatures suggestive of younger stellar ages ($ < 4\,$Myr compared with $ > 20\,$Myr for sources with single-peaked Lyman-$ \alpha $ profiles). We argue that double-peaked Lyman-$ \alpha $ profiles may reflect weaker feedback in extremely young starbursts due to the delayed onset of supernova explosions. We also investigate whether outflow properties can be reliably inferred by fitting simple expanding shell models. We find that these models yield unphysical parameters in many cases even when they can reproduce the Lyman-$ \alpha $ profiles to high precision, and urge caution when inferring outflow properties from such models.

Ultrahot Jupiters are predicted to show inverted temperature-pressure (T-P) profiles in the presence of optical absorbers such as TiO and VO. An inverted T-P profile of KELT-7b was recently detected, in line with these predictions, but such diagnoses are known to be model-dependent. We used CHEOPS, TESS, and literature data to characterize the atmosphere of KELT-7b, reassess its T-P profile, measure its albedo, and search for distortions in its CHEOPS transit light curve due to stellar rotation. We jointly fitted CHEOPS and TESS data to measure the occultation depths and modeled CHEOPS transits including gravity darkening. Emission and transmission retrievals were performed, and the albedo was calculated in the CHEOPS and TESS passbands. Thermochemical-equilibrium retrievals yield a non-inverted T-P profile, while free-chemistry retrievals yield an inverted profile with likely unphysical TiO/VO abundances. A 3D GCM supports a TiO-driven inversion. We report a low geometric albedo of $A_\mathrm{g} = 0.05 \pm 0.06$, consistent with inefficient heat redistribution and supported by a GCM with magnetic drag. CHEOPS data provide no constraint on the sky-projected orbital obliquity. Retrieval results strongly depend on the chemical framework. Free-chemistry fits are better but risk unphysical solutions for ultrahot Jupiters. We applied a coherent stellar variability correction to CHEOPS and TESS data; future observations would benefit from similar treatment.

Stars are fossils that retain the history of their host galaxies. Carbon and heavier elements are created inside stars and are ejected when they die. From the spatial distribution of elements in galaxies, it is therefore possible to constrain the physical processes during galaxy formation and evolution. This approach, Galactic archaeology, has been popularly used for our Milky Way Galaxy thanks to a vast amount of data from the Gaia satellite and multi-object spectrographs, and now can also be applied to very distant galaxies with the James Webb Space Telescope (JWST) - extra-galactic archaeology. In these studies the most important factor is the input stellar physics, namely nucleosynthesis yields and binary physics, which predominantly determine the model predictions. In this review I give a summary of stellar nucleosynthesis, and how they are tested with the observations in the Milky Way. Then I show how chemical enrichment of galaxies can be calculated, and show some results with the latest nucleosynthesis yields.

We present ALMA CO(1-0) and CO(3-2) observations of a powerful young radio galaxy, PKS 0023-26, hosted by a far-infrared bright galaxy. The galaxy has a luminous optical AGN and a very extended distribution of molecular gas. We used these observations (together with available CO(2-1) data) to trace the impact of the AGN across the extent of the radio emission and beyond on scales of a few kpc. Despite the strength of the optical AGN, the kinematics of the cold molecular gas is strongly affected only in the central kpc, and is more weakly affected around the northern lobe. We found other signatures of the substantial impact of the radio AGN, however. Most notably, extreme line ratios of the CO transitions in a region aligned with the radio axis indicate conditions very different from those observed in the undisturbed gas at large radii. The non-detection of CO(1-0) at the location of the core of the radio source implies extreme conditions at this location. Furthermore, on the scale of a few kpc, the cold molecular gas appears to be wrapped around the northern radio lobe. This suggests that a strong jet-cloud interaction has depleted the northern lobe of molecular gas, perhaps as a result of the hot wind behind the jet-induced shock that shreds the clouds via hydrodynamic instabilities. The higher gas velocity dispersion and molecular excitation that we observed close to this location may then be the result of a milder interaction in which the expanding jet cocoon induces turbulence in the surrounding interstellar medium. These results highlight that the impact of an AGN can manifest itself not only in the kinematics of the gas, but also in molecular line ratios and in the distribution of the gas. Although the radio plasma and the cold molecular gas are clearly coupled, the kinetic energy that is transferred to the ISM is only a small fraction of the energy available from the AGN.

Radio pulsating white dwarf (WD) systems, known as WD pulsars, are non-accreting binary systems where the rapidly spinning WD interacts with a low-mass companion producing pulsed non-thermal emission that can be observed across the entire electromagnetic spectrum. Only two such systems are known: AR Sco and eRASSU J191213.9$-$441044. Here we present the discovery of a third WD pulsar, SDSS J230641.47$+$244055.8. The optical spectrum is dominated by molecular bands from an M-dwarf companion, with additional narrow emission lines from the Balmer series and He I. The long-term optical light-curve folded on its orbital period ($P_\mathrm{orb} = 3.49$ h) exhibits large scatter (roughly 10 per cent). High-cadence photometry reveals a short period signal, which we interpret to be the spin period of the WD primary ($P_\mathrm{spin} \simeq 92$ s). The WD spin period is slightly shorter than that of AR Sco ($\rm \sim 117$ s), the WD pulsar prototype. Time-resolved spectroscopy reveals emission from the irradiated companion and Na I absorption lines approximately tracing its centre of mass, which yields a binary mass function of $f(M) \simeq 0.2 {\rm M_\odot}$. The H$\alpha$ emission includes a low-amplitude broad component, resembling the energetic emission line flashes seen in AR Sco. Using spectral templates, we classify the companion to be most likely a $\rm M4.0\pm 0.5$ star with $T_\mathrm{\rm eff} \approx 3300$ K. Modelling the stellar contribution constrains the secondary mass ($0.19\,{\rm M_\odot}\lesssim M_2\lesssim 0.28\,{\rm M_\odot}$), system distance ($\simeq1.25\,{\rm kpc}$), and inclination ($i \simeq 45-50^\circ$). We discuss the proposed evolutionary scenarios and summarize the observational properties of all three known WD pulsars, establishing a benchmark for identifying and classifying future members of this emerging class.

This paper provides an overview of the spectroscopic data obtained by the JWST Guaranteed Time Observations (GTO) program 1192, "The Physics and Chemistry of PDR Fronts", including an inventory of the spatially resolved dust, gas, and molecular content in the Horsehead nebula and the NW filament of NGC~7023. We demonstrate the unique capability of this high spatial resolution data set to elucidate the evolution of gas and dust at the interface between stars and their natal clouds at the scale at which the physics and chemistry occur. The Disassociation Regions (PDRs) in the Horsehead nebula and the North West (NW) filament of NGC 7023 were mapped with a spectral resolution 1000-3000 and a spatial resolution of ~2e-4 pc between 0.97-28um. Spectra extracted from template regions yield a large number of atomic, ionized, and molecular lines. Full line lists and extracted spectra for all 10 regions are provided through CDS. Absorption from H$_2$O, CO$_2$, and CO ices are identified in 3 regions in NGC 7023. In this overview, we have focused on the spectral region between 3 and 5 um which is dominated by emission from aromatic and aliphatic carbon bonds to illustrate the power of the data set. We confirm the entrainment of aromatic carbonaceous species in the photo-evaporative flow from the PDR surface into the H\,{\sc{ii}} region in the Horsehead. No aliphatic emission is present in the outflow, indicating the complete removal of aliphatic bonds when exposed to strong UV fields. There is a clear detection of deuterium substitution in the carbon bonds. Aliphatic D-substitution is more efficient relative to aromatic D-substitution, ranging from N$_{D}$/N$_{H}$ ~ 0.1-0.3 for aliphatics compared to ~0.03 for the aromatics.

Symbiotic stars, binary pairs with a cool giant fueling accretion onto a hot compact companion, offer unique insights to our understanding of stellar evolution. Yet, only a few hundred symbiotic stars are confirmed. Here, we report on a new search for symbiotic star candidates in Gaia Data Release 3 (GDR3), based entirely on the archive's astrometric, photometric, and spectroscopic information. To begin our search, we identified known symbiotic stars in GDR3 and assessed their absolute magnitude and colors, which are dominated by the cool giant. We also considered measures of astrometric quality that might be affected by binary motion in these systems. Finally, from those sources with Gaia spectroscopic data, we built a low-resolution spectral template that characterizes the unique features of these systems, including H$\alpha$ emission from interaction with the giant's wind and radiation from the hot star. We then queried the full GDR3 archive for sources with spectroscopic data that are bright (< 17 mag in G-band), have modest relative parallax uncertainties (< 20%), and fall within a region of color-magnitude space characteristic of red giants, keeping only sources with spectra that quantitatively match our template. A machine-learning algorithm, trained on known symbiotic stars, produced a new catalog of 1,674 sources. From cross-matches with infrared and X-ray surveys, we present 25 of these sources as particularly compelling candidates for new symbiotic stars.

We consider dynamics and multi-frequency emission patterns of relativistic van Allen belts - particles trapped in the magnetosphere of neutron stars and white dwarths. We account for synchrotron radiative losses and effects of relativistic beaming of radiation. The system is non-Hamiltonian (non-energy conserving): this results in a wide non-scalable variety of spectral and temporal behaviors. There are three types of trapped particles' trajectories: (i) oscillating (particles experience multiple bounces between magnetic bottles); (ii) precipitating (particles fall onto the star with finite transverse momentum); (iii) freezing (particles lose their transverse motion before falling onto the star). The separation between regimes (i) and (ii) depends both on the ratio of the bounce time to cooling time at magnetic equator $\tau_{ 0} $, $\eta _0 = R_0/( c \tau_0) \leq 1$, as well as the initial pitch angle $\alpha_0$; regimes (i) and (ii) are separated at $ \alpha_{0, crit} \sim \eta_0^{3/10}$. Resulting emission patterns show large variety: single or double peaked, and/or flat hat with sharp walls. Multi-frequency profiles - in optical and X-ray bands - can be used to get information about physical (magnetic field strength, injection point) and geometrical (dipolar angle and the line of sight) properties.

Observations of some quiescent black hole X-ray binaries have revealed an excess of mid-infrared (MIR) emission above that expected from their donor stars. In one system, V404 Cygni, this excess has been variously suggested to arise from the accretion disk, circumbinary material, or a compact relativistic jet. Here we present simultaneous James Webb Space Telescope (JWST), Atacama Large Millimeter/submillimeter Array (ALMA), and complementary multi-wavelength observations undertaken to resolve this uncertainty. We observed large-amplitude 21 $\mu$m variability on short timescales with JWST, particularly a dramatic flare which swiftly rose to $\approx 2.4$ mJy, over 10 times the lowest observed MIR flux density. Similar variability was simultaneously observed from radio to X-ray wavelengths with other facilities throughout the campaign. This variability and the flat radio/mm/MIR spectral index ($\alpha = 0.04 \pm 0.01$) suggest that the MIR excess in V404 Cyg does not arise from the accretion disk or circumbinary material but is instead dominated by synchrotron radiation from a jet which persists into quiescence. This result reinforces the ubiquity of the disk-jet connection in accreting black holes across a range of masses and accretion rates.

We present the first lattice simulations of the nonlinear evolution after axion inflation by self-consistently incorporating currents arising from Schwinger pair production. The tachyonically amplified gauge fields trigger the growth of Schwinger currents, leading to universal values for the conductivity and magnetic field at the onset of strong backreaction and subsequent quenching of gauge field production. We show that the Schwinger effect (prematurely) saturates gauge field production, thereby diminishing the prospects of axion inflation magnetogenesis as a viable solution for blazar observations.

The delayed conversion of a neutron star (NS) into a highly magnetized hybrid star (HS) ejects ~ 0.01Msun of outer NS layers with kinetic energy ~ 10^50 erg (Keranen et al. 2005); the Quark-Nova (QN) ejecta. If this occurs in isolation (or after the SN ejecta had dissipated), the rotational energy of the HS powers a transient with rise and decay times of a few days, consistent with Luminous Fast Blue Optical Transients (LFBOTs) (Ouyed 2025). Breakup of the QN ejecta creates optically thin fragments that release trapped radiation over tens of minutes timescales, producing optical flares with luminosities comparable to the LFBOT peak value. X-ray emission arises from HS spin-down energy leaking through the gradually transparent ejecta, while radio emission is due to synchrotron radiation from the ejecta's interaction with the ambient medium. Our model reproduces key features of AT2022tsd and other LFBOTs, including AT2020xnd, AT2020mrf, and AT2018cow. Flares are absent or dwarfed when fragmentation occurs before or during the LFBOT peak emission phase, as likely in AT2018cow. Our model provides a novel explanation for aged NSs with long spin periods that exhibit youthful magnetar-like activity, such as ASKAP J1832-0911 (Wang et al. 2024). It also carries implications for heavy-element enrichment in the universe, given the formation rate of LFBOTs and the neutron-rich composition of the QN ejecta (Jaikumar et al. 2007, Ouyed et al. 2009).

While JWST has provided us with the opportunity to probe the atmospheres of potentially-habitable planets, observations of the TRAPPIST-1 system have shown us that active stars severely complicate efforts at studying their planets. GJ 12b is a newly-discovered temperate (Teq ~ 300 K), Earth-sized (Rp = 0.96 +/- 0.05 Earth radii) planet orbiting an inactive M dwarf that might be a good alternate to the TRAPPIST-1 planets for atmospheric characterization. In this paper, we use MAROON-X radial velocities to measure a mass of 0.71 +/- 0.12 Earth masses for GJ 12b. We also find moderate evidence that the planet has an eccentric (e ~ 0.16) orbit. GJ 12b's mass results in a planetary density comparable to or less dense than Earth, possibly indicating the presence of water or a low bulk iron mass fraction. With its low mass, GJ 12b is likely within reach of JWST transmission spectroscopy observations, making it an excellent target for determining the location of the cosmic shoreline. Its low mass may mean that the planet could have trouble retaining its primary atmosphere during the star's active pre-main-sequence phase. However, if it has a heightened eccentricity, it may be able to sustain a secondary atmosphere through tidally-induced volcanism.

Small temperate planets are prime targets for exoplanet studies due to their possible similarities with the rocky planets in the Solar System. M dwarfs are promising hosts since the planetary signals are within our current detection capabilities. Gliese 12 b is a Venus-sized temperate planet orbiting a quiet M dwarf. We present here the first precise mass measurement of this small exoplanet. We performed a detailed analysis using HARPS-N, ESPRESSO, and CARMENES radial velocities, along with new and archival TESS, CHEOPS, and MuSCAT2/3 photometry data. From fitting the available data, we find that the planet has a radius of $R_\mathrm{p} = 0.904^{+0.037}_{-0.034} \,\mathrm{R_\oplus}$ and a mass of $M_\mathrm{p} = 0.95^{+0.26}_{-0.27} \,\mathrm{M_\oplus}$ (a $3.6\sigma$ measurement of the semi-amplitude $K=0.70^{+0.19}_{-0.20}\,\mathrm{m\,s^{-1}}$), and is on an orbit with a period of $12.761421 \pm 0.000047\,\mathrm{d}$. A variety of techniques were utilised to attenuate stellar activity signals. Gliese 12 b has an equilibrium temperature of $T_\mathrm{eq}=315 \pm 7\,\mathrm{K}$, assuming an albedo of zero, and a density consistent with that of Earth and Venus ($\rho_\mathrm{p}=7.0^{+2.3}_{-2.1}\,\mathrm{g\,cm^{-3}}$). We find that Gliese 12 b has a predominantly rocky interior and simulations indicate that it is unlikely to have retained any of its primordial gaseous envelope. The bulk properties of Gliese 12 b place it in an extremely sparsely-populated region of both mass--radius and density--$T_\mathrm{eq}$ parameter space, making it a prime target for follow-up observations, including Lyman-$\alpha$ studies.

Analyzing filaments and cores in molecular clouds is key to understanding galactic star formation and its environmental dependence. This paper studies the properties and distribution of dense cores and filaments in the Ophiuchus molecular cloud, with a focus on the L1688 hub-filament system (HFS) and its star formation potential. We extracted sources and filaments from Herschel images and a 13.5 arcsec resolution surface density map using the getsf method, identified prestellar cores among the extracted sources, evaluated core mass segregation, and constructed the core mass function (CMF). We derived properties of the filaments from their radial surface density profiles, constructed the filament linear density function (FLDF), and assessed the mass distribution in the L1688 HFS to estimate the core and filament formation efficiencies (CFE, FFE). The results suggest that the filament-dominated core formation is a key mechanism in star formation within the system.

We present a 6GHz radio continuum image of the galaxy cluster Abell 2744 ($z = 0.307$) obtained with the Karl G. Jansky Very Large Array (VLA) as part of the VLA Frontier Fields program, whose goal is to explore the radio continuum emission from high-redshift galaxies that are magnified by foreground, massive galaxy clusters. With an rms noise of $\approx 1 \mu$Jy beam$^{-1}$, at the phase center, and sub-arcsec angular resolution ($\theta_{1/2}=0.82$ arcsec), this is the deepest and most detailed radio image of Abell 2744 ever obtained. A total of 93 sources are detected with a peak signal-to-noise ratio $\geq5$, of which 46 have optical/near-infrared (IR) counterparts with available redshift, magnification ($\mu$), and stellar mass (${M}_*$) estimates. The radio sources are distributed over a redshift from 0.15 to 3.55, with a median redshift value of $z = 0.93$, and a median stellar mass of ${M}_{*}= 2.3 \times 10^{10}\,\rm{M}_{\odot}$. A comparison between the radio-based star formation rates (SFRs) and those derived from ultraviolet-to-near IR data reveals that the radio SFRs are, on average, an order of magnitude higher than the ultraviolet-to-near IR SFRs. We look for radio counterparts of the so-called "Little Red Dots (LRDs)'' galaxies at $z\approx6$ in Abell 2744, but find no significant detections. After stacking, we derive a 3$\sigma$ upper limit to the 6GHz radio luminosity of LRDs of $4.1\times 10^{39}\,\rm erg\,s^{-1}$. Finally, we present a sample of 22 moderately/strongly lensed galaxies ($\mu \gtrsim 2$) in the VLA Frontier Fields survey, which is adequate to zoom into star formation processes of main sequence galaxies at $z\approx 1-2$.

The interpretation of the ultra-high-energy cosmic ray spectrum (UHECRs) and composition suggests a suppression of the flux below $\sim$1 EeV, as observed by the Pierre Auger Observatory and Telescope Array. A natural explanation for this phenomenon involves magnetic confinement effects. We investigate the possibility that UHECRs self-generate the magnetic turbulence necessary for such confinement via current-driven plasma instabilities. Specifically, we show that the electric current produced by escaping UHECRs can excite a non-resonant streaming instability in the surrounding plasma. This instability reduces the diffusion coefficient in the source environment, effectively trapping particles with energies $E \lesssim 0.6$ EeV $\mathcal{L}_{45}^{1/2} R_{\text{Mpc}}^{-1} \lambda_{10}^{2}$ for times exceeding the age of the Universe. Here, $\mathcal{L}_{45}$ is the source luminosity in units of $10^{45}$ erg/s, $R_{\text{Mpc}}$ is the radial size in Mpc, and $\lambda_{10}$ is the intergalactic magnetic field coherence length in units of 10 Mpc. We discuss in detail the conditions, in terms of source luminosity, initial magnetic field, and the environment in which this phenomenon occurs, that need to be fulfilled in order for self-confinement to take place near a source of UHECRs. By modeling a population of UHECR sources with a luminosity function typical of extragalactic gamma-ray sources, we connect the spectrum of escaping particles to the luminosity distribution. Furthermore, we calculate the contribution of these confined particles to cosmogenic neutrino production, finding consistency with current observational constraints. Our results suggest that self-induced turbulence may play an important role in shaping the UHECR spectrum and, in particular, may account for the flux suppression near their sources, offering a promising framework for interpreting current observations.

We investigate the impact of dark matter (DM) annihilation on the global 21-cm signal during the dark ages and cosmic dawn eras. The 21-cm line provides a complementary probe for studying the nature of dark matter beyond standard cosmological observables. In the standard $\Lambda$CDM framework, the expected absorption amplitude of the dark ages global 21-cm signal is approximately $-42\, \rm mK$. However, energy injection from DM annihilation can significantly heat and ionize the intergalactic medium, potentially altering or even erasing this absorption feature. We evaluate the thermal and ionization history of the gas to derive an upper bound on $f_\chi^2 \langle \sigma v \rangle / M_\chi$ using the dark ages signal, which is free from astrophysical uncertainties. After incorporating observational and theoretical uncertainties arising from future lunar-based experiments and variations in cosmological parameters, respectively -- we obtain a conservative upper limit of $f_\chi^2\langle\sigma v\rangle/M_\chi \lesssim 10^{-28}~\rm cm^3\,s^{-1}\,\rm GeV^{-1}$. This constraint is nearly an order of magnitude stronger than the bounds derived from Planck 2018 data.

With the astrophysics community working towards the first observations and characterizations of Earth-like exoplanets, interest in space-based nulling interferometry has been renewed. This technique promises unique scientific and technical advantages by enabling direct mid-infrared observations. However, concept studies of nulling interferometers often overlook the impact of systematic noise caused by instrument perturbations. Earlier research introduced analytical and numerical models to address instrumental noise and, building on these results, we reproduce key simulations and report that the noise in the differential output of nulling interferometers follows a non-Gaussian distribution. The presence of non-Gaussian noise challenges the validity of classical hypothesis tests in detection performance estimates, as their reliance on Gaussian assumptions leads to overconfidence in detection thresholds. For the first time, we derive the true noise distribution of the differential output of a dual Bracewell nulling interferometer, demonstrating that it follows iterative convolutions of Bessel functions. Understanding this noise distribution enables a refined formulation of hypothesis testing in nulling interferometry, leading to a semi-analytical prediction of detection performance. This computationally efficient instrument model, implemented in a publicly available codebase, is designed for integration into science yield predictions for nulling interferometry mission concepts. It will play a key role in refining key mission parameters for the Large Interferometer For Exoplanets (LIFE).

Exoplanet surveys like Kepler, TESS, and K2 have shown that planetary systems are common in our galaxy. These surveys, along with several others, have identified thousands of planetary candidates, with more than five thousand having already been confirmed. Many of these planetary systems host multiple planets. As we discover more multiplanet systems, notable trends begin to appear in the data. We use kernel density estimation (KDE) to analyze the period ratios of adjacent planet pairs in multiplanet systems in the most recent Kepler, TESS and K2 data, paying particular attention to pairs in first order mean motion resonance (MMR). We compare a recent Kepler catalog with the DR25 data release. We also compare TESS and K2 against the latest Kepler catalog. To verify the significance of our findings against selection bias, we perform Monte Carlo simulations of multiplanet systems, finding an excess of planet pairs near the 2.2, 2 (2:1), and 1.5 (3:2) period ratios in the Kepler, K2, and TESS catalogs, all exceeding the 99% confidence interval when constrained to orbital periods less than 25 days. We identify two planet pairs orbiting M dwarf stars in a very tight ratio, as well as two likely misidentified planet pairs.

Noise maps from CMB experiments are generally statistically anisotropic, due to scanning strategies, atmospheric conditions, or instrumental effects. Any mis-modeling of this complex noise can bias the reconstruction of the lensing potential and the measurement of the lensing power spectrum from the observed CMB maps. We introduce a new CMB lensing estimator based on the maximum a posteriori (MAP) reconstruction that is minimally sensitive to these instrumental noise biases. By modifying the likelihood to rely exclusively on correlations between CMB map splits with independent noise realizations, we minimize auto-correlations that contribute to biases. In the regime of many independent splits, this maximum closely approximates the optimal MAP reconstruction of the lensing potential. In simulations, we demonstrate that this method is able to determine lensing observables that are immune to any noise mis-modeling with a negligible cost in signal-to-noise ratio. Our estimator enables unbiased and nearly optimal lensing reconstruction for next-generation CMB surveys.

Atmospheric layer structure is a primary factor affecting the precision of single-point satellite positioning. The assumption of electromagnetic wave rectilinear propagation hinders the accurate implementation of ionospheric and tropospheric corrections, whereas curvilinear positioning methods fully account for ray deflection. This study aims to derive partial derivative formulas for theoretical travel time with respect to latitude, longitude, elevation, and velocity models by formulating electromagnetic wave travel time equations under a coordinate-based one-dimensional layered velocity model. Subsequently, a linearized LSQR method is employed to invert station coordinates, receiver clock biases, and electromagnetic wave velocities at the bottom of the ionosphere and troposphere using over six sets of observations. This replaces conventional ionospheric/tropospheric pseudorange corrections in single-point positioning, establishing a novel spherical coordinate refraction travel time positioning method. The classical straight-line pseudorange positioning is reformulated into a time-of-flight positioning approach, and the positioning accuracy differences between straight-line and spherical coordinate refraction travel time methods are compared. By integrating classical ionospheric and tropospheric models to construct corresponding refractive index models and combining them with curvilinear ray tracing methods, the inherent theoretical limitations of positioning can be effectively mitigated.

The equivalence principle and the inverse-square law of gravity could be violated at short distances ($10^{-6}$ to $10^{-12}$ meters) by scalars sporting a coupling $g_N$ to nucleons and mass $\mathrm{eV}\lesssim m_\phi\lesssim\rm MeV$. We show for the first time that stringent bounds on the existence of these scalars can be derived from the observed cooling of nearby isolated neutron stars (NSs). Although NSs can only be used to set limits comparable to the classic SN 1987A cooling bound in the case of pseudoscalars such as the QCD axion, the shallow temperature dependence of the scalar emissivity results in a huge enhancement in the effect of $\phi$ on the cooling of cold NSs. As we do not find evidence of exotic energy losses, we can exclude couplings down to $g_N\lesssim 5 \times 10^{-14}$. Our new bound supersedes all existing limits on scalars across six orders of magnitude in $m_\phi$. These conclusions also extend to Higgs-portal models, for which the bound on the scalar-Higgs mixing angle is $\sin\theta\lesssim 6\times 10^{-11}$.

The QCD axion is a well-motivated hypothetical particle beyond the Standard Model (SM) and a compelling dark matter candidate. Its relic abundance is highly sensitive to the thermal history of the universe when the temperature is around the QCD confinement scale. Meanwhile, the NANOGrav Collaboration has reported evidence for a stochastic gravitational wave background, which could originate from a supercooled first-order phase transition (FOPT) with a nucleation temperature around the O(MeV-GeV) scale. We explore how such an FOPT might alter the evolution of the QCD axion. Our findings suggest that it could induce the axion to go through a short stage of mini kinetic misalignment. Moreover, in some parameter regime, the formation of QCD axion domain walls becomes generically expected. This has intriguing implications for both the existence of the QCD axion and the FOPT interpretation of the NANOGrav signal.

Magnetic reconnection and turbulence in magnetically-dominated environments have been proposed as important nonthermal particle acceleration mechanisms that generate high energy particles and associated emissions. While the acceleration to high energy that produces the power-law energy distribution has drawn strong interest, recent studies actively discuss pre-acceleration, or injection, to a sufficient energy for a sustained and prolonged Fermi-like acceleration. The injection process is important for determining the fraction of nonthermal particles and energy partition between thermal and nonthermal particles. We review recent advances in understanding the injection mechanisms responsible for populating these nonthermal power-law spectra, and conclude with an outlook for studies and applications of injection models.

We present Bayesian inference results from an extensive injection-recovery campaign to test the validity of three state of the art quasicircular gravitational waveform models: \textsc{SEOBNRv5PHM}, \textsc{IMRPhenomTPHM}, \textsc{IMRPhenomXPHM}, the latter with the \textsc{SpinTaylorT4} implementation for its precession dynamics. We analyze 35 strongly precessing binary black hole numerical relativity simulations with all available harmonic content. Ten simulations have a mass ratio of $4:1$ and five, mass ratio of $8:1$. Overall, we find that \textsc{SEOBNRv5PHM} is the most consistent model to numerical relativity, with the majority of true source properties lying within the inferred 90\% credible interval. However, we find that none of the models can reliably infer the true source properties for binaries with mass ratio $8:1$ systems. We additionally conduct inspiral-merger-ringdown (IMR) consistency tests to determine if our chosen state of the art waveform models infer consistent properties when analysing only the inspiral (low frequency) and ringdown (high frequency) portions of the signal. For the simulations considered in this work, we find that the IMR consistency test depends on the frequency that separates the inspiral and ringdown regimes. For two sensible choices of the cutoff frequency, we report that \textsc{IMRPhenomXPHM} can produce false GR deviations. Meanwhile, we find that \textsc{IMRPhenomTPHM} is the most reliable model under the IMR consistency test. Finally, we re-analyze the same 35 simulations, but this time we incorporate model accuracy into our Bayesian inference. Consistent with the work in Hoy et al. 2024 [arXiv: 2409.19404], we find this approach generally yields more accurate inferred properties for binary black holes with less biases compared to methods that combine model-dependent posterior distributions based on their evidence, or with equal weight.

The undergraduate-led Polarization-modUlated Laser Satellite Experiment (PULSE-A) at the University of Chicago seeks to demonstrate the feasibility of circular polarization shift keyed satellite-to-ground laser communication. PULSE-A's low-cost open-source bus serves as the backbone of the mission and has been designed in tandem with the Payload, with design driven by strict requirements for pointing accuracy, component alignment, power demand, and thermal stability. This work presents the design and testing of the PULSE-A bus. The spacecraft bus was designed to fill two major needs: (1) to meet the requirements of the PULSE-A mission, and (2) to be easily configurable for future missions that desire enhanced capabilities over other low-cost open-source designs. At its core, the bus features dual BeagleBone Black Industrial compute units, selected for their flight heritage, integrated via a PC/104 header standard. PULSE-A implements Goddard Space Flight Center's core Flight System (cFS), which takes a modular software architecture approach and is built in C. The use of C as the primary language aligns with the expertise of the University of Chicago's Computer Science department, allowing for ease of development by PULSE-A's undergraduate flight software team. The CubeSat structure utilizes Gran Systems' 3U frame, modified to accommodate openings for various ports and deployable components. Inside, the avionics stack uses the PC/104 standard quad rails, which terminate in PULSE-A's custom-designed Payload Box that houses all of the Payload components and optical fiber runs. This work also covers the techniques and iterative engineering processes used to develop the thermal control and dissipation mechanisms for the specific requirements, under volume, mass, and temperature-range constraints.

The lifetime of primordial black holes (PBHs), which were formed in the early universe, can be extended by the memory burden effect. Light PBHs may be in existence today and become a candidate for dark matter (DM). We assume that DM is made of thermally produced weakly interacting massive particles (WIMPs), WIMPs produced via Hawking radiation of PBHs, and survived PBHs via the memory burden effect. In addition, feebly interacting massive particles (FIMPs) are considered as an alternative to WIMPs. Various similar studies have been conducted either with or without considering the memory burden effect. In this study, we simultaneously account for thermally produced WIMPs or FIMPs and the memory burden effect of PBHs to explain the relic abundance of DM. We show that this abundance is highly sensitive to the memory burden effect in PBHs.

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In $\beta$-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously-varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-$\beta$ equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in $\beta$-equilibrium, the pressure in the mixed phase is not constant even for the case of fully-local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill-defined when applied to the mixed phase.