Papers for Tuesday, May 20 2025

We present the results of the GIGAPYX-4600 image sensor in space environment, more specifically under different types of irradiations (protons and heavy ions). The GIGAPYX-4600 is a state-of-the-art 46M pixel multi-purpose, backside illuminated CMOS image sensor. The sensor features high-speed (200 fps), low-noise (< 2e-rms), rolling shutter readout. It has been fabricated using 65 nm node CMOS technology, making use of capacitive deep trench isolation, thus exhibiting good MTF as well as excellent dark current characteristics. The GIGAPYX image sensor family is meant to be easily scalable thanks to a novel use of stitching technology. The assessed sensor features an impressive 46 Mpixels, but the sensor family is meant to be scaled up to 220 M pixels. The idea of this study was to investigate the radiation hardness of a commercially available off-the-shelf (COTS) image sensor that could be suitable for space-borne applications, such as earth observation or satellite vicinity surveillance. In the near future a radiation hard readout electronic for this sensor family will be made available off-the-shelf by Pyxalis. This presentation will overview the different sensor key performances evolutions after proton irradiation up to a total fluence of 2.3e11 p+/cm${}^2$ (62 MeV): dark current, dark current non-uniformity (DCNU), noise, non-linearity, saturation charge and photo-response non-uniformity (PRNU). As expected, degradations mostly occur on the dark current, DCNU and temporal noise. The orders of magnitude of the degradation are in the range of the already published high performance CIS technologies. Results obtained from heavy ions irradiations will demonstrate that the GIGAPYX is not only latch-up free at least up to 57 Mev.cm2/mg, but resistant to the blooming effects induced by a SET at pixel level thanks to its capacitive deep trench isolations. Various SEE effects have been studied, demonstrating encouraging results for such COTS device in order to fly, in particular in GEO orbit. All these radiations results will be used as inputs in designing space camera based on the GIGAPYX sensors, compatible with multiple-mission types.

Recent observations from the ESA Gaia satellite and with the ESO VLT, have identified the presence of a population of young, 0.5 to 2 Gyr old, stars in the halo and in dwarf spheroidal galaxies surrounding the Milky Way. It suggests that MW dwarf galaxies, currently devoid of gas, had, until recent times, enough gas to sustain a burst of star formation. The recent loss of gas coincides with their arrival in the vicinity of the Milky Way, in agreement with orbital predictions from Gaia that indicate that most dwarf galaxies reached the Milky Way halo less than 3 Gyr years ago. This completely changes the interpretation of their dynamics, mass, and dark matter content.

Neutrinos experience collective flavor conversion in extreme astrophysical environments such as core-collapse supernovae (CCSNe). One manifestation of collective conversion is slow flavor conversion (SFC), which has recently attracted renewed interest owing to its ubiquity across different regions of the supernova environment. In this study, we systematically examine the evolution of kinematic decoherence in a dense neutrino gas undergoing SFC, considering lepton number asymmetries as large as $30\%$. Our findings show that the neutrino gas asymptotically evolves toward a generic state of coarse-grained flavor equilibration which is constrained by approximate lepton number conservation. The equilibration occurs within a few factors of the inverse vacuum oscillation frequency, $\omega^{-1}$, which corresponds to (anti)neutrinos reaching near flavor equipartition after a few kilometers for typical supernova neutrino energies. Notably, the quasi-steady state of the neutrino number densities can be quantitatively described by the neutrino-antineutrino number density ratio $n_{\bar{\nu}_e}/n_{\nu_e}$ alone. Such a simple estimation opens new opportunities for incorporating SFC into CCSN simulations, particularly in regions where SFC develops on scales much shorter than those of collisions.

The existence of axion-like particles (ALPs) can be probed from their signatures in the Cosmic Microwave Background (CMB) due to the photon-ALP resonant conversion over the mass range of ALPs that matches with the effective mass of photons in the plasma in the astrophysical systems. Such a conversion can also occur in the Milky Way halo and disk and can cause a unique spatial and spectral distortion. The signal is highly non-Gaussian and cannot be measured precisely by the usual power-spectrum approach. We devise a new technique to search for this signal from the upcoming full-sky CMB experiment LiteBIRD using its multi-frequency band using a template-based spatial profile of the ALP distortion signal. This technique captures the large-scale non-Gaussian aspects of the ALP distortion signal in terms of a spatial template and makes it possible to search for any non-zero ALP signal. We show that the inference of the ALP coupling using the template-based technique from LiteBIRD can provide constraints on the coupling constant approximately $ g_{a\gamma} < 6.5 \times 10^{-12} \, \mathrm{GeV}^{-1}$ for ALP masses below $10^{-14}$ eV at 95\% confidence interval which is an order of magnitude better than the current bounds from CERN Axion Solar Telescope (CAST) at $g_{a\gamma} < 6.6 \times 10^{-11} \, \mathrm{GeV}^{-1}$, This shows the capability of future multi-band CMB experiment LiteBIRD in opening the discovery space towards physics beyond the standard model.

We are searching for hot, constant-color, offset optical flares in the Zwicky Transient Facility (ZTF) data stream that are ${>}10''$ from any galaxy in public imaging data from the PanSTARRS survey. Here, we present the first discovery from this search: AT 2024puz, a luminous multiwavelength transient offset by $5\,$kpc from a ${\sim}10^8\,M_\odot$ galaxy at $z=0.356$ with a low-moderate star formation rate. It produced luminous $10^{44.79 \pm 0.04}\,{\rm erg\,s}^{-1}$ optical/UV emission that evolved on a ${\sim}20\,$day timescale, as well as $10^{44.12\pm0.03}\,{\rm erg\,s}^{-1}$ X-ray emission with a photon-index $\Gamma=1.7$. No associated radio or millimeter emission was detected. We show that the early-time optical emission is likely powered by reprocessing of high-energy, accretion-powered radiation, with a possible contribution from a shock in a dense circum-transient medium. If the shock is dominant at early-times, the circum-transient medium has a mass ${\sim}0.1-1\,M_\odot$, radius $10^{15}\,$cm, and a density profile shallower than ${\sim}r^{-1}$. A near-infrared excess appears at late-times and is suggestive of reprocessing within a wind or other circum-transient medium. The X-rays are most consistent with a central engine. We suggest that AT 2024puz may be associated with an accretion event onto a $50-10^5\,M_\odot$ BH, where the lower masses are preferred based on the large projected offset from the host galaxy. AT2024puz exhibits properties similar to both luminous fast blue optical transients (LFBOTs) and tidal disruption events (TDEs), but is intermediate between them in its energetics and evolution timescale. This highlights the need for broader exploration of the landscape of hot optical transients to trace their origins.

Pulsar timing arrays (PTAs) are approaching the sensitivity required to resolve for the first time gravitational wave (GW) signals from individual supermassive black hole (SMBH) binaries. However, the large uncertainty in the localisation of the source will make the identification of its host environment challenging. We show how the posterior probability function of binary parameters inferred by the standard GW analysis can be converted into distributions of apparent magnitudes of the host galaxy in the infrared (IR) and optical bands. We do so for two different scenarios: one in which the host is a regular early-type galaxy (ETG), and one in which the binary resides in an active galactic nucleus (AGN). We estimate the binary parameter space PTAs can cover in the near and intermediate future and estimate whether their hosts will be detectable in all-sky electromagnetic surveys. A PTA with a baseline of 20 years and 116 pulsars, resembling the upcoming data release of the International Pulsar Timing Array (IPTA)can detect binaries out to a luminosity distance of 2Gpc (which corresponds to a redshift of $z \approx 0.36$) under the most optimistic scenario for detection, while a PTA with a baseline of 30 years and 200 pulsars can reach out to a maximum distance slightly greater than 3Gpc ($z \approx 0.53$). We find that the host galaxies of all binaries that are detectable by a PTA with a baseline of 20 years are above the threshold for the WISE and SuperCOSMOS, and therefore they are expected to be present in those surveys, if they lie outside the Milky Way plane. 2MASS becomes incomplete for hosts of binary systems more massive than $10^{9.8}$ solar masses located at a luminosity distance from Earth greater than 1 Gpc. The EM surveys become slightly more incomplete when we consider the sensitivity of PTAs with baselines of 25 and 30 years, as PTAs can detect binaries to larger distances.

The halo occupation distribution is a cornerstone in understanding galaxy formation within the large-scale this http URL the other hand, in the context of the large-scale structure of the Universe, voids are regions with singular characteristics, given their extension, the low number of objects that inhabit their interior and their own dynamics. Furthermore, the HOD of these regions exhibits significant deviations from the average, as shown by several studies in semi-analytical models, hydrodynamic simulations, and observations. This paper investigates the temporal evolution of the HOD in cosmic voids from redshift z = 0.5 to the present day. We aim to understand how the void environment shapes galaxy occupation over time and to identify the factors driving the unique HOD observed in these underdense regions. We use the MDPL2-SAG semi-analytic galaxy catalog and identify spherical cosmic voids at different redshifts. We analyze the HOD within voids across ten simulation snapshots, comparing it with the global HOD. Furthermore, we trace the evolution of z = 0 void galaxies and their host halos back in time, examining their formation times and local density evolution. Our results reveal that the HOD within voids is consistently lower than the average HOD at all redshifts considered. Tracking z = 0 void galaxies, we find that their HOD remains lower than average throughout the redshift range, with no significant redshift evolution in this relative difference. Analysis of halo formation histories shows that lower-mass void halos are younger, while massive void halos are older compared to the average. Spatially, lower-mass halos are found in the inner void regions, whereas massive halos are located closer to void boundaries. The local density of massive void halos decreases significantly towards z = 0, contrasting with the relatively stable local density of lower-mass void halos.

Solar-type members of the rich, nearly Solar-age and Solar-metallicity M67 open cluster are systematically investigated for ultraviolet variability. We utilize archival Galaxy Evolution Explorer (GALEX) data which features several imaging observation epochs spanning 5 years. Stars in or suspected of being in binary systems are avoided as well as stars that are blended in GALEX data, leading to a sample of 66 Solar-type stars. We assess variability over a variety of timescales that probe flares and longer-term trends that could be due to rotation and activity cycles. We do not find conclusive evidence for variability and determine that Solar-type members of M67 do not display >30% near-ultraviolet variability over timescales ranging from days to years. Furthermore, within 50-second cadence lightcurves generated for each of the imaging epochs we find no near-ultraviolet flares that are >~2x the quiescent stellar near-ultraviolet emission level; the implied ultraviolet flare rate derived from this study is in mild tension with that derived for stars observed by GALEX in the primary Kepler field. This M67 GALEX study presents one of the most comprehensive ultraviolet datasets currently available for probing continuum emission variability for old Sun-like stars; the planned NASA UVEX mission has the potential to dramatically expand upon this work.

Some carbon-rich Wolf-Rayet stars (WC stars) show an infrared excess from dust emission. Dust forms in the collision of the WC wind with a companion star's wind. As this dust is carried towards the ISM at close to the WCd wind speed and the binary continues through its orbit, a spiral structure forms around the system. The shape depends on the orbital eccentricity and period, as well as stellar parameters like mass-loss rates and terminal wind speeds. Imaging of the WCd binary WR 140 with JWST/MIRI revealed at least 17 concentric dust shells surrounding the binary. We present new JWST imaging of four additional WCd systems (WR 48a, WR 112, WR 125, and WR 137) that were imaged in 2024. In this analysis, we show that the dust is long-lived, surviving for at least 130 years, but more than 300 years in some systems. Longer duration measurements are limited by sensitivity. Regular spacing of dust features confirms the periodic nature of dust formation, consistent with a connection to binary motion. We use these images to estimate the proper motion of the dust, finding the dust to propagate out to the interstellar medium with motion comparable to the wind speed of the WC stars. In addition to these results, we observe unusual structures around WR 48a, which could represent dusty clumps shaped by photoevaporation and wind ablation like young proplyd objects. These results demonstrate that WC dust is indeed long-lived and should be accounted for in galactic dust budgets.

The orbital properties of the (as-yet) small population of hot Jupiters with nearby planetary companions provide valuable constraints on the past migration processes of these systems. In this work, we explore the likelihood that dynamical perturbations could cause nearby inner or outer companions to hot Jupiter to leave the transiting plane, potentially leaving these companions undetected despite their presence at formation. Using a combination of analytical and numerical models, we examine the effects of stellar evolution on hot Jupiter systems with nearby companions and identify several possible outcomes. We find that while inner companions are generally unlikely to leave the transiting plane, outer companions are more prone to decoupling from the hot Jupiter and becoming non-transiting, depending on the system's initial orbital architecture. Additionally, we observe a range of dynamical behaviors, including overall stability, inclination excitation, and, in some cases, instability leading to the ejection or collision of planets. We also show that the effect of stellar obliquity (with respect to the mean planet of the planets) is to amplify these effects and potentially cause outer companions to attain non-mutually-transiting configurations more often. Our results highlight the complex dynamical pathways shaping the architectures of hot Jupiter systems.

We discuss the possibility that dark energy arises from a strongly-coupled Higgs-Yang-Mills set of interacting fields in the non-perturbative regime. We choose the simplest $SU(2)$ representation, which is compatible with the Cosmological Principle. One of the components of the Higgs doublet act as an effective quintessence scalar field interacting with both gravity and the gauge field. We devise a multiple time scale approach to solve the equations of motion through a hierarchy of the couplings, utilizing exact solutions in terms of Jacobi elliptic functions. We observe that the time scale of variation of the Hubble constant is the slowest one, while, for the scalar field, assuming that its self-coupling is smaller than the coupling of the gauge field, represents an intermediate time scale. From the consistency of the Friedmann equations, we show how the effect of the scalar field is to give origin to dark energy with a proper equation of state corrected by a very small time-dependent term. Finally, agreement with cosmological data for the dark energy density is shown without relying on the fine-tuning of the physical constants of the model.

The Rapid Iterative FiTting (RIFT) parameter inference algorithm provides a simulation-based inference approach to efficient, highly-parallelized parameter inference for GW sources. Previous editions of RIFT have conservatively optimized for robust inference about poorly constrained observations. In this paper, we summarize algorithm enhancements and operating point choices to enable inference for more exceptional compact binaries. Using the previously-reported RIFT/asimov interface to efficiently perform analyses on events with reproducible settings consistent with past work, we demonstrate that the latest version of RIFT can efficiently analyze events with multiple costly models including the effects of precession or eccentricity.

We present SpiderCat, a multi-wavelength catalog of all publicly known compact binary millisecond pulsars (MSPs) in the Galactic field. These systems, colloquially known as "spiders", consist of neutron stars in tight orbits with low-mass companions, which are gradually ablated by the pulsar wind. SpiderCat includes both primary subclasses - redbacks and black widows - distinguished by companion mass, as well as candidates and peculiar systems such as transitional, huntsman and tidarren MSPs. As of this initial release, SpiderCat contains 109 entries: 32 redbacks, 49 black widows, 23 redback candidates, and 5 black widow candidates. In this paper, we compile and summarize key parameters for each system, including spin and orbital properties, and multi-wavelength data from radio, optical, X-ray, and gamma-ray observations. An interactive, publicly accessible web interface (this https URL) enables exploration and visualization of the data. The rapid growth of the number of known spiders, accelerated by the Fermi-LAT survey and its ability to identify MSPs in gamma-rays, has opened the door to population-level studies. Utilizing SpiderCat, we analyze trends in spin period, orbital period, companion mass, emission properties, and spatial distribution. SpiderCat serves as a dynamic, multi-wavelength repository for this unique class of binary pulsars, facilitating new discoveries and constraints on pulsar evolution, particle acceleration, and the neutron star equation of state.

In 2022 and 2023 the hydrogen comae of two long period comets, C/2017 K2 (PanSTARRS) and C/2022 E3 (ZTF), were observed with the Solar Wind ANisotropies (SWAN) all-sky hydrogen Lyman-alpha camera on the SOlar and Heliosphere Observer (SOHO) satellite. SWAN obtains nearly daily full-sky images of the hydrogen Lyman-alpha distribution of the interstellar hydrogen as it passes through the solar system yielding information about the solar wind and solar ultraviolet fluxes that eat away at it by ionization and charge exchange. The hydrogen comae of comets, when of sufficient brightness, are also observed. Water production rates have been calculated over time for each of these comets, covering about 6 months mostly of the post-perihelion period of C/2017 K2 (PanSTARRS) and about 3 months around perihelion of C/2022 E3 (ZTF).

We present the first large-scale, high-resolution simulations of dusty, star formation feedback-driven galactic outflows. Using the Cholla hydrodynamics code, we investigate dust sputtering in these environments for grains ranging in size from $1-0.001~{\mu\mathrm{m}}$. We compare results for two feedback models: one representative of low-redshift nuclear starburst galaxies and one similar to high-redshift main sequence galaxies. In general, our simulations show that multi-phase outflows are capable of safely transporting a vast majority of their dust to large distances ($\sim10~\textrm{kpc}$) from the disk. This work also shows that environmental shielding in cool gas clouds boosts dust survival rates significantly. The evolutionary path of dust depends strongly on grain size. Large grains ($a\geq0.1~{\mu\mathrm{m}}$) can be transported efficiently in all phases. Smaller grains, however, experience significant destruction in the hotter phases. $0.001~{\mu\mathrm{m}}$ grains in particular are quickly sputtered in all but the coolest gas, resulting in these grains strongly tracing the cool phase in outflows. These results may also indicate the importance of in-situ formation mechanisms, such as shattering, for the small dust grains and PAHs observed in emission throughout outflows in nearby galaxies. Surprisingly, we find that the hot phase dominates the transport of dust that survives to populate the CGM.

We present a comprehensive analysis of the digest2 parameters for candidates of the Near-Earth Object Confirmation Page (NEOCP) that were reported between 2019 and 2024. Our study proposes methods for significantly reducing the inclusion of non-NEO objects on the NEOCP. Despite the substantial increase in near-Earth object (NEO) discoveries in recent years, only about half of the NEOCP candidates are ultimately confirmed as NEOs. Therefore, much observing time is spent following up on non-NEOs. Furthermore, approximately 11% of the candidates remain unconfirmed because the follow-up observations are insufficient. These are nearly 600 cases per year. To reduce false positives and minimize wasted resources on non-NEOs, we refine the posting criteria for NEOCP based on a detailed analysis of all digest2 scores. We investigated 30 distinct digest2 parameter categories for candidates that were confirmed as NEOs and non-NEOs. From this analysis, we derived a filtering mechanism based on selected digest2 parameters that were able to exclude 20% of the non-NEOs from the NEOCP while maintaining a minimal loss of true NEOs. We also investigated the application of four machine-learning (ML) techniques, that is, the gradient-boosting machine (GBM), the random forest (RF) classifier, the stochastic gradient descent (SGD) classifier, and neural networks (NN) to classify NEOCP candidates as NEOs or non-NEOs. Based on digest2 parameters as input, our ML models achieved a precision of approximately 95% in distinguishing between NEOs and non-NEOs. Results. Combining the digest2 parameter filter with an ML-based classification model, we demonstrate a significant reduction in non-NEOs on the NEOCP that exceeds 80%, while limiting the loss of NEO discovery tracklets to 5.5%. Importantly, we show that most follow-up tracklets of initially misclassified NEOs are later correctly identified as NEOs.

{We investigate the multi-phase gas surrounding QSOs traced by 33 absorption lines (e.g., Ly$\alpha$, C\,\textsc{iv}, Fe\,\textsc{ii}, Mg\,\textsc{ii}, etc.) in the stacked spectra of background sources, using the early data release from the Dark Energy Spectroscopic Instrument. Our analysis reveals that the equivalent width (\( W \)) of metal absorption lines decreases with increasing redshift, following an overall trend described by $W \propto (1+z)^{-4.0\pm 2.2}$. This indicates cosmic enrichment of the gas associated with QSOs, with a more pronounced enrichment observed at larger scales (a few Mpc), reflecting a rapid metal enrichment of cosmic web. Additionally, the \( W \) of these absorption lines decreases with distance ($D$) from QSOs, which can be effectively characterized by a two-halo model. Compared to the projected two point correlation function of galaxies at similar redshifts, low-ionization ions exhibit similar clustering scales, while high-ionization ions show a significantly more extended spatial distribution. We also find that $W_{\text{FeII}}/W_{\text{MgII}}$ increases towards lower redshifts and larger scales, which can be attributed to evolving star formation histories and/or changes in initial mass function for galaxies. By leveraging multiple absorption tracers, we conduct the first comprehensive investigation of diffuse, multiphase gas from the circumgalactic medium to cosmological scales, offering new insights into baryon cycles and the transport of metals throughout cosmic time.

Optical interferometric image reconstruction is a challenging, ill-posed optimization problem which usually relies on heavy regularization for convergence. Conventional algorithms regularize in the pixel domain, without cognizance of spatial relationships or physical realism, with limited utility when this information is needed to reconstruct images. Here we present PIRATES (Polarimetric Image Reconstruction AI for Tracing Evolved Structures), the first image reconstruction algorithm for optical polarimetric interferometry. PIRATES has a dual structure optimized for parsimonious reconstruction of high fidelity polarized images and accurate reproduction of interferometric observables. The first stage, a convolutional neural network (CNN), learns a physically meaningful prior of self-consistent polarized scattering relationships from radiative transfer images. The second stage, an iterative fitting mechanism, uses the CNN as a prior for subsequent refinement of the images with respect to their polarized interferometric observables. Unlike the pixel-wise adjustments of traditional image reconstruction codes, PIRATES reconstructs images in a latent feature space, imparting a structurally derived implicit regularization. We demonstrate that PIRATES can reconstruct high fidelity polarized images of a broad range of complex circumstellar environments, in a physically meaningful and internally consistent manner, and that latent space regularization can effectively regularize reconstructed images in the presence of realistic noise.

Planets residing within the hot-Neptune Desert are rare, and studying their atmospheres can provide valuable insights into their formation and evolutionary processes. We present the atmospheric characterization of the first known ultra-hot Neptune, LTT-9779 b, using transmission spectroscopic observations obtained with the HST/WFC3 G141 and G102 grisms. Using the Iraclis pipeline and TauREx3 retrieval code, we find that LTT-9779 b likely possesses a H/He-dominated primary atmosphere with an opaque aerosol layer and the pure cloudy, flat-line model is rejected with approximately 2.7-$\sigma$ confidence. Although we do not find conclusive evidence supporting the presence of any molecular species, we place 95% confidence level upper limits on the volume mixing ratios (VMRs) of hydroxyl radical (OH) and iron hydride (FeH) at $7.18\times10^{-2}$ and $1.52\times10^{-8}$, respectively. Notably, the retrieval results are inconsistent with predictions from equilibrium chemistry models, which favor higher $\rm H_2O$ abundances over OH. This discrepancy suggests that disequilibrium processes, such as photochemistry or vertical mixing, may have altered the atmospheric composition. Comparisons between HST, Spitzer and JWST data reveal no evidence of temporal variations in the atmospheric composition of the terminator region. Our results highlight the need for higher-resolution spectroscopy and secondary eclipse observations to resolve LTT-9779 b's temperature-pressure (T-P) profile and chemical inventory definitively.

Coronal simulations of the solar maximum struggle with poor numerical stability and low computational efficiency since the magnetic field is more complex and stronger and coronal structures evolve more rapidly. This paper aims to enhance the numerical stability of the time-evolving COCONUT coronal model to mitigate these issues, to evaluate differences between the time-evolving and quasi-steady-state coronal simulation results, and to assess the impact of spatial resolution on global MHD coronal modelling of solar this http URL enhancing the positivity-preserving property of the time-evolving COCONUT, we employ it to simulate the evolution of coronal structures from the solar surface to 0.1 AU over two CRs around the May 2024 solar storm event. These simulations are performed on unstructured meshes containing 6.06, 1.52, and 0.38 M cells to assess the impact of grid resolution. We also conduct a quasi-steady-state coronal simulation, treating the solar surface as a rigidly rotating spherical shell, to demonstrate the impact of magnetic flux emergence and cancellation in global coronal simulations. Comparison with observations further validates the reliability of this this http URL paper demonstrates that incorporating magnetic field evolution in inner-boundary conditions can significantly improve the fidelity of global MHD coronal simulations around solar maximum. The simulated magnetic field strength using a refined mesh with 6.06 M cells can be more than 40% stronger than that in the coarser mesh with 0.38 M cells. A time step of 5 minutes and the mesh containing 1.5 M cells can effectively capture the evolution of large-scale coronal structures and small-sized dipoles. Thus, this model shows promise for accurately conducting real-time global coronal simulations of solar maximum, making it suitable for practical applications.

Asteroseismology, the study of stellar vibrations, is a method which can probe the structure deformation and internal rotation of stars. Salient among the seismic inferences of rotation from TESS observations are TIC 408165734, whose equatorial rotation rate is 10\% faster than the pole, and TIC 307930890, which has significant radial shear and shows a decreasing spin rate outward through its envelope. We also measure structural deformation in fifteen stars, nine of which are oblate, a finding consistent with expectations for relatively fast-rotating, non-magnetic stars. The difference between polar and equatorial radii in TIC 47639058 is 130 times larger than that for the Sun. The remaining six stars display splittings consistent with a prolate shape (surprisingly), possibly indicating the presence of equatorial toroidal magnetic fields. These inferences provide constraints for numerical simulations and new insights to guide theories of $\delta$ Scuti structure and rotation.

WISE 0855-0714 is the coldest known brown dwarf, located 2.28 pc from the solar system. Discovered by the Wide-Field Infrared Survey Explorer (WISE) in 2014 (Luhman 2014), the object is of interest to scientists because of its low temperature ($\approx270$ K), proximity to the solar system, small mass ($\sim3-10\: M_{J}$), and high proper motion. The first observations of W0855 by WISE in 2010 are heavily contaminated by a background source. With 10.5 years of observations following the NEOWISE reactivation in 2013 (Mainzer et al., 2014), we present a robust analysis of W0855's flux and color unobstructed by this background source. We obtain W1 = 19.3 and W1-W2 = 5.4 magnitudes with an error of 0.37 magnitudes.

We investigate a sample of 14 galaxy clusters from the GOGREEN and GCLASS (GG) spectroscopic datasets within the redshift range $(0.87 \leq z \leq 1.37)$ and cluster masses $\mathrm{M}_{200} \gtrsim 2\times 10^{14}$ \hm. Using the highly effective GalWeight technique for cluster membership assignment developed by our own team, we derive the dynamical parameters of these clusters through the virial mass estimator. We examine the velocity dispersion-cluster mass relation $(\sigma \mathrm{MR})$ for the GG cluster sample. We find, $\log{\sigma_{200}} = (2.94\pm0.02) + (0.37\pm0.07)\log{\mathrm{M}_{200}}$ with an intrinsic scatter of $(\sigma_\mathrm{int} = 0.02 \pm 0.02)$. Our results demonstrate that the $(\sigma \mathrm{MR})$ relation is consistent with predictions from cosmological simulations, highlighting the reliability of the GalWeight technique for cluster membership assignment. Furthermore, the $(\sigma \mathrm{MR})$ validates the robustness of the virial mass estimator in accurately recovering cluster masses and associated parameters. Importantly, our findings confirm that velocity dispersion can be used directly to estimate cluster mass without relying on dynamical mass estimators.

In hot and ultra-hot Jupiters, stellar irradiation is a primary driver of atmospheric circulation and the wave structures that sustain it. We aim to investigate how variations in radiative and dynamical timescales influence global flow regimes, atmospheric circulation efficiency, and the interplay of wave structures across a sample of hot Jupiters. In particular, we explore a previously predicted transition in the global flow regime, where enhanced stellar irradiation suppresses the smaller-scale wave and eddy features that feed into superrotating jets and ultimately leads to simpler, day-to-night dominated flows. We simulate a suite of eight well-studied hot Jupiters with the THOR general circulation model, spanning equilibrium temperatures from about $1100$ K to $2400$ K. We develop a wavelet-based analysis method to decompose simulated wind fields into their underlying wave modes, which we validate on analytical examples. As a preliminary exploration of the flow regime of ultra-hot Jupiters, we perform an additional simulation for WASP-121b, where the mean molecular weight is set to represent an atmosphere dominated by atomic hydrogen. Our results confirm that increasing stellar irradiation diminishes atmospheric heat redistribution efficiency and weakens contributions from smaller-scale modes that are critical to sustain superrotation. As equilibrium temperatures rise, large-scale modes dominate the atmospheric circulation, driving a transition from jet-dominated flows toward day-to-night circulation. Additionally, by artificially lowering the mean molecular weight, we partially restore circulation efficiency and reintroduce a more complex, multi-scale flow pattern. These findings refine our understanding of how atmospheric circulation evolves with increasing irradiation and composition changes, offering a more nuanced framework for interpreting hot and ultra-hot Jupiter atmospheres.

We consider a basic planetary system composed by a Sun like star, a Jupiter-like planet an a Neptune-like planet in a wide range of orbital configurations not limited to the hierarchical case. We present atlases of resonances showing the domains of approx. 1300 mutual mean-motion resonances (MMRs) and their link to chaotic and regular dynamics. Following a semi-analytical method for the study of the secular dynamics we found two regimes for equilibrium configurations: one for low mutual inclinations were equilibrium is related to oscillations of the difference between the pericenter longitudes around 0 or 180 degrees, and another for high mutual inclinations where the equilibrium is given by defined values of the argument of the pericenters equal to integer multiples of 90 degrees. By numerical integration of the full equations of motion we calculate the fundamental frequencies of the systems in their diverse configurations and study their dependence with the orbital elements. According to the analysis of the fundamental frequencies we found two dynamical regimes depending on the initial mutual inclination and the limit between the two regimes occurs at some critical inclination 30<ic<40 defined by the occurrence of the secular resonance g1=g2. For i<ic the dynamics is analogue a the classic secular model for low (e,i) with well defined three fundamental frequencies and free and forced modes, conserving quasi constant the mutual inclination. For i>ic the dynamics is completely different with increasing changes in mutual inclination and emerging combinations of the fundamental frequencies and, depending on the case, dominated by the secular resonance or the vZLK mechanism.

We have measured the dust emissivity index beta for 21 infrared-bright sources (including several gravitationally lensed galaxies) at 1.5 < z < 4.2 using Atacama Large Millimeter/submillimeter Array (ALMA) 101-199 GHz data sampling the Rayleigh-Jeans side of the SED. These data are largely insensitive to temperature variations and therefore should provide robust measurements of beta. We obtain a mean beta of 2.2 with a standard deviation of 0.6 that is at the high end of the range of values that had previously been measured in many galactic and extragalactic sources. We find no systematic variation in beta versus redshift. We also demonstrate with a subset of our sources that these higher beta values have significant implications for modelling dust emission and in particular for calculating dust masses or the wavelength at which dust becomes optically thick.

The He I line at 1.08 $\mu$m is a valuable tracer of atmospheric escape in exoplanet atmospheres. We expand past networks used to predict the absorbing He(2$^3S$) by including, firstly, processes that involve H$_2$ and some molecular ions and, secondly, the interaction of photoelectrons with the atmosphere. We survey the literature on the chemical-collisional-radiative processes that govern the production-loss of He(2$^3S$). We simulate the atmospheric outflow from the Neptune-sized GJ 436 b by coupling a hydrodynamic model that solves the bulk properties of the gas and a Monte Carlo model that tracks the energy degradation of the photoelectrons. We identify Penning ionization of H as a key He(2$^3S$) loss process at GJ 436 b and update its rate coefficient to a value consistent with the most recent available cross sections. The update affects notably the predicted strength of the He I line. For GJ 436 b, photoelectron-driven processes (mainly ionization and excitation) modify the He(2$^3S$) population in layers too deep to affect the in-transit spectrum. The situation might be different for other atmospheres though. The spectral energy distribution of GJ 436 has a strong effect on the predicted in-transit signal. The published non-detections of the He I line for GJ 436 b are reasonably consistent with our model predictions for a solar-metallicity atmosphere when the model adopts a recently proposed spectral energy distribution. The interpretation of the He I line at 1.08 $\mu$m is model-dependent. Our revised network provides a general framework to extract more robust conclusions from measurements of this line, especially in atmospheres where H$_2$ remains abundant to high altitudes. We will explore additional, previously-ignored processes in future work.

The detection of SO2 in both WASP-39 b and WASP-107 b recently brought more attention to the modeling of photochemistry in exoplanets. However, sulfur kinetics data is lacking in the literature for the full C/H/O/N/S system. The networks used to model this chemistry neglect the coupling between sulfur and other C/H/O/N species. We aimed to integrate sulfur kinetics to our previously developed C_0-C_2/H/O/N chemical network, with the inclusion of its coupling to carbon and nitrogen chemistry, for conditions between 500 - 2500 K and 100 - 10^-6 bar, and any atomic composition. To achieve this reliability, we used a dual approach, deriving the network from other available combustion networks and from original ab initio calculations where data was lacking. This was performed together with an extensive validation of the network on 1606 experimental measurements from the literature on the combustion and pyrolysis of multiple sulfur compounds such as H2S, CH3SH, CS2 and OCS. To examine the consequences of this new chemical network on exoplanets atmospheric studies, we generated abundance profiles for GJ 436 b, GJ 1214 b, HD 189733 b, HD 209458 b, WASP-39 b and WASP-107 b using the 1D kinetic model FRECKLL, calculated the corresponding transmission spectra using TauREx 3.1 and compared these results with other chemical networks used in exoplanet modeling with sulfur. The coupling between carbon and sulfur chemistry is found to be impactful on both abundance profiles and observables, with CH2S being its key species. CS2 abundance is found to be probably much higher than anticipated in current kinetic networks for exoplanets. The detection of CS2 in TOI-270 d highlights the importance of using validated chemical networks to improve the reliability of our models, particularly in the JWST era. Combustion and pyrolysis data are largely available tools that reveal to be very useful for this task.

Recent multi-wavelength observations by JWST and ALMA are unveiling both ionized and neutral ISM components in high-redshift ($z>6$) galaxies. In this work, we investigate the origin of rest-frame far-infrared [OIII]88 $\mu$m and [CII]158 $\mu$m emission by performing zoom-in cosmological simulations of dwarf-galaxy progenitors at $z=9-13$. Our simulations incorporate on-the-fly radiative transfer at sub-pc ($\sim$ 0.1 pc) resolution, allowing us to resolve the multi-phase ISM. We compute emission lines on a cell-by-cell basis, taking into account local temperature, density, metallicity, radiation field strength, column density, and spectral hardness of radiation bins. We find that [OIII] predominantly arises from centrally located ionizing bubbles with temperatures of $\sim (1-5)\times 10^4\,\mathrm{K}$ and high ionization parameters of $\log U_{\mathrm{ion}} \simeq -1.5$. In contrast, [CII] is produced in the surrounding dense neutral regions at $\sim 5\times 10^3\,\mathrm{K}$, which are heated by strong FUV radiation from the central stellar clusters. This spatial arrangement leads to large local variations in [OIII]/[CII], ranging from $\sim$ 100 to 0.01. Our galaxy reproduces the global ratio [OIII]/[CII]$\sim5-30$, consistent with recent ALMA detections at $z>6$ without invoking enhanced O/C abundance ratios. We further derive that [OIII]/[CII] linearly scales with the mass and density ratios of ionized to neutral gas, $M_{\rm HII}/M_{\rm HI}$ and $n_{\rm HII}/n_{\rm HI}$ and show that the [OIII]/[CII] ratio typically changes from 5.7 to 0.3 from high-z to low-z. For future synergies of JWST and ALMA, we derived $M_{\rm HII}/M_{\rm HI}$ for observed $z >6$ galaxies using ${\rm H}\beta$ and [CII] and show the validity of our scaling relations.

In this work we use particle-in-cell (PIC) numerical simulations to study interaction of a spatially uniform electron beam with a rotational magnetic hole in a form of a Harris current sheet. We vary width of the Harris current sheet to investigate how this affects the quasi-linear relaxation, i.e. plateau formation of the bump-on-tail unstable electron beam. We find that when width of the Harris current sheet approaches and becomes smaller than the electron gyro-radius, quasi-linear relaxation becomes hampered and a positive slope in the electron velocity distribution function (VDF) persists. We explain this by the effects of non-conservation of electron magnetic moment, which, as recent works suggest, can maintain the positive slope of the VDF. In part, this can explain why some electron beams in the solar wind travel much longer distances than predicted by the quasi-linear theory, at least in those cases when the electron beams slide along the current sheets that are abundant when the different-speed solar wind streams interact with each other.

Observations of fast X-ray transients (FXRTs) with the Einstein Probe have successfully led to the discovery of some unusual extragalactic optical transients. EP241021a is a newly discovered FXRT that was featured by a significant bump around ten days in both optical and X-ray bands. This timescale and the exceptionally high peak bolometric luminosity up to $\sim \rm 10^{44}erg~s^{-1}$ of the optical bump make it somewhat similar to fast blue optical transients, but still distinctive from them by its relatively red color. We then suggest that the multi-wavelength bump of EP241021a could represent an explosion-type transient, while the underlying power-law decaying component of the optical and X-ray emission as well as the total radio emission are produced by a moderately relativistic jet. By fitting the observed multi-wavelength light curves, it is found that the explosion ejecta that produce the thermal optical emission can have a mass of $\sim0.03~M_{\odot}$, an expanding velocity of $\sim0.25~c$, and an optical opacity of $\sim12~\rm cm^2g^{-1}$, which was continuously powered by a rapidly rotating and highly magnetized neutron star (NS; i.e., a magnetar). In addition to heating the explosion ejecta, the magnetar also provided the dominant contribution to the observed X-ray rebrightening through the non-thermal emission of its wind. These properties suggest that the explosion may result from a catastrophic collapse/merger of a compact star system, which led to the formation of a millisecond magnetar, and the possible progenitor could be an accreting white dwarf (WD) or a binary consisting of double WDs, double NSs, or a WD and an NS.

Recent JWST observations have revealed a growing population of galaxies at $z>4$ with elevated nitrogen-to-oxygen ratios. These "N/O-enhanced" galaxies (NOEGs) exhibit near- to super-solar N/O at sub-solar O/H, clearly deviating from the well-established scaling relation between N/O and O/H observed in local galaxies. The origin of this abundance anomaly is unclear. Interestingly, local globular clusters also exhibit anomalous light-element abundances, whose origin remains debated. In this work, we compare the chemical abundance patterns of 22 known NOEGs at $0\lesssim z\lesssim 12$ -- primarily discovered with JWST -- to those observed in local globular clusters. We find striking similarities in the abundances of C, N, O, Fe, and He between the two populations. The similar abundance patterns support the scenario in which globular cluster stars formed within proto-cluster environments -- similar to those traced by NOEGs -- that were self-enriched. Indeed, the enhancement in N/O in early galaxies appears to be only found in dense stellar environments with $\Sigma _{\star}\gtrsim 10^{2.5}~M_\odot~{\rm pc^{-2}}$, as expected for the progenitors of globular clusters in the Milky Way, and similar to those of star clusters identified in strongly lensed high-redshift galaxies. Furthermore, we find a tentative positive correlation between N/O ratios and stellar mass among NOEGs. The apparent high occurrence rate of NOEGs at high redshift is consistent with the picture of cluster-dominated star formation during the early stages of galaxy evolution. Measuring chemical abundances across diverse stellar environments in high-redshift galaxies will be crucial for elucidating the connection between NOEGs and globular clusters.

The upcoming space-based gravitational wave (GW) observatory, LISA, is expected to detect GW signals from supermassive black hole (SMBH) mergers occurring at high redshifts. However, understanding the origin and growth of SMBHs in the early Universe remains an open problem in astrophysics. In this work, we utilize the First Light And Reionization Epoch Simulations (FLARES), a suite of cosmological hydrodynamical zoom-in simulations, to study SMBH mergers at $5 \lesssim z \lesssim 10$ across a wide range of environments. Most mergers in FLARES involve secondary SMBHs near the seed mass ($m_{seed} \approx 1.5 \times 10^{5} M_{\odot}$) while primary SMBHs span up to $10^{9} M_{\odot}$, resulting in mass ratios from $q \sim 10^{-4}$ to $1$, with a peak at $q \sim 1$. The number of mergers increases rapidly towards lower redshifts, and the comoving total number density scales with overdensity as $n_{merger} = 10^{-3.80} (1 + \delta)^{4.56}$. Denser regions host more massive mergers, with higher merger redshifts and lower mass ratios. Within the FLARES redshift range, LISA is expected to detect mergers with $10^{5} \lesssim M_{tot} / M_{\odot} \lesssim 10^{8}$ and $q \gtrsim 10^{-2}$, corresponding to a detection rate of 0.030 $yr^{-1}$ for events with signal-to-noise ratio $SNR \geq 10$. Our study demonstrates the sensitivity of GW predictions at high redshifts to SMBH seed models and merger time delays, highlighting the need for improved modeling in future cosmological simulations to maximize LISA's scientific return.

The formation and early evolution of Jupiter played a pivotal role in sculpting the large-scale architecture of the solar system, intertwining the narrative of Jovian early years with the broader story of the solar system's origins. The details and chronology of Jupiter's formation, however, remain elusive, primarily due to the inherent uncertainties of accretionary models, highlighting the need for independent constraints. Here we show that by analyzing the dynamics of Jupiter's satellites concurrently with its angular momentum budget, we can infer Jupiter's radius and interior state at the time of proto-solar nebula's dissipation. In particular, our calculations reveal that Jupiter was $2$ to $2.5$ times as large as it is today, 3.8 million years after the formation of the first solids in the solar system. Our model further indicates that young Jupiter possessed a magnetic field of approximately $B_{\rm{J}}^{\dagger} \approx 21$ mT (a factor of $\sim50$ higher than its present-day value) and was accreting material through a circum-Jovian disk at a rate of $\dot{M} = 1.2-2.4$ Jupiter masses per million years. Our findings are fully consistent with the core-accretion theory of giant planet formation and provide an evolutionary snapshot that pins down properties of the Jovian system at the end of the protosolar nebula's lifetime.

In the bulge of M31, the Chandra observations discovered a possible black hole (BH) ultracompact X-ray binary (UCXB) Seq.1 with an orbital period of 7.7 minutes and a maximum X-ray luminosity $L_{\rm X}=1.09^{+0.02}_{-0.01}\times10^{38}~ \rm erg\,s^{-1}$ in the $0.5-8$ keV band. The minimum orbital period of the BH UCXBs predicted by the standard magnetic braking (MB) model is longer than 8.3 minutes. In this work, we investigate whether the convection- and rotation-boosted (CARB) MB prescription can account for the formation of a BH UCXB like Seq.1. Our detailed stellar evolution models indicate that the CARB MB law can drive isolated BH-main sequence (MS) binaries to evolve toward BH UCXBs with an orbital period of $7.7~ \rm minutes$, in which a low-mass white dwarf transfers the material onto a BH in a short-term mass transfer episode, producing an X-ray luminosity of $10^{38}~\rm erg\,s^{-1}$. We also obtain an initial parameter space of BH-MS binaries as the progenitors of Seq.1 in the donor-star masses and orbital periods plane, which can be applied to future population synthesis simulations. If Seq.1 is indeed a BH UCXB, future spaceborne gravitational wave (GW) detectors can detect the low-frequency GW signals from this source, and a tidal disruption event will be expected after 0.12 Myr.

One of the most mysterious results from observations of the James Webb Space Telescope (JWST) is the detection of numerous, high-redshift, very red, extremely compact, broad-line sources termed ``little red dots'' (LRDs). It is unclear whether the LRDs belong to an active galactic nucleus (AGN) or simply a collection of very compact star clusters. We build spectral energy distributions (SEDs) for 29 LRDs at $z \approx 3-8.5$ based on JWST photometric and spectroscopic observations. We find that the V-shaped SEDs of these LRDs exhibit a roughly similar break frequency at $\nu_{\rm b}=10^{14.96\pm0.06}$ Hz, which corresponds to $\lambda_{\rm b}=3287_{-424}^{+487} \textÅ$ in the rest frame. We propose that this unique SED can be explained by the combination of an inner standard disk and an outer gravitationally unstable accretion disk with Toomre parameter $Q\sim1$. The outer disk has a temperature of $\sim2000-4000$ K for typical AGN parameters, which can well reproduce the near-infrared to optical bump as observed in LRDs. This model can naturally explain the strong infrared to optical emission and the V-shaped SED with a similar break frequency $\simeq 10^{15}$ Hz for LRDs without invoking strong dust extinction or unusual stellar contribution from a host galaxy. Most LRDs stay in sub-Eddington state based on the SED modeling, which are intrinsically weak in X-rays.

We apply two machine learning methods, a CNN deep-leaning model and a gradient-boosting decision tree, to estimate individual cluster optical depths from observed properties derived from multiple complementary datasets. The models are trained and tested with simulated N-body derived halo catalogs and synthetic full-sky CMB maps designed to mirror data from the DESI and Simons Observatory experiments. Specifically, the thermal Sunyaev-Zel'dovich (tSZ) and CMB lensing convergence, along with cluster virial mass estimates are used as features to train the machine learning models. The predicted optical depths are combined with kinematic Sunyaev-Zel'dovich (kSZ) measurements to estimate individual cluster radial peculiar velocities. The method is shown to recover an unbiased estimate of the pairwise velocity statistics of the simulated cluster sample. The model's efficacy is demonstrated for halos with mass range $10^{13} M_{\odot} < M_{200} < 10^{15} M_{\odot}$ over a redshift range $0<z<1$, and validated in the presence of primary CMB, instrument noise, lensing convergence noise, and potential uncertainties in halo virial mass estimates. We apply the method to ACT CMB data, using ACT DR4 component-separated maps for tSZ and CMB lensing and ACT DR5 maps for kSZ, in conjunction with galaxy clusters observed in the SDSS DR15 spectroscopic survey. We demonstrate that the machine learning approach is an effective one to analyze data from current and upcoming CMB experiments such as Simons Observatory and CCAT, and galaxy surveys, such as DESI and Roman, for which the pairwise velocity statistics can provide valuable insights into the properties of neutrinos and gravity on cosmic scales.

High-velocity clouds (HVCs), which are gas clouds moving at high velocity relative to the galactic disk, may play a critical role in galaxy evolution, potentially supplying gas to the disk and triggering star formation. In this study, we focus on the nearby face-on barred spiral galaxy M83, where high spatial resolution, high-sensitivity CO (1-0) data are available. We identified molecular clouds and searched for clouds with velocities deviating by more than 50km/s from the disk velocity field as HVCs. A total of 10 HVCs were detected -nine redshifted and one blueshifted -- clearly highlighting an asymmetry in their velocity distribution. These HVCs have radii of 30-80 pc, masses on the order of $10^5 M_\odot$, and velocity dispersions of 3-20 km/s, displaying a tendency toward higher velocity dispersion compared to disk molecular clouds in M83. Most of the HVCs do not overlap with the candidates of supernova remnants, and the energy needed to drive HVCs at such high velocities exceeds single supernova energy. Together with the asymmetry in their velocity distribution, we thus conclude that most of the HVCs found in this study are inflow from outside the M83's disk.

Mass loss on the red giant branch (RGB) influences stellar evolution, properties of stellar populations, and Galactic chemical enrichment, yet remains poorly constrained observationally. Current models provide limited insight into how stellar properties, particularly how metallicity and mass, affects RGB mass loss. Here, I introduce a new observational approach that uses the age-velocity-dispersion relation and the lower-mass boundary of red giants as precise evolutionary markers. These markers, informed by Galactic evolution, allow us to construct observational isochrones for field stars. By comparing masses of RGB stars and red clump (RC) stars at the same age in the Kepler sample, I derive empirical measurements of integrated RGB mass loss at several points in age and metallicity. Combining these new observational measurements with open-cluster studies, I showed that the integrated mass loss on the RGB decreases with metallicity, and may also decrease with stellar mass. The average mass loss rate, which accounts for RGB lifetimes and the initial mass differences between RGB and RC stars at the same age, also show a similar trend. These findings challenge current mass-loss prescriptions widely adopted in stellar evolutionary models, since none of them is able to produce the observed mass-loss trend without widely adjusting free parameters. This highlights an urgent need to revise mechanisms that govern RGB mass loss.

To fully understand the star formation process, we are compelled to study it in a variety of environments. Of particular interest are how star formation and the resulting initial mass function (IMF) vary as a function of metallicity. We have observed an embedded young cluster in Sh2-284 (hereafter S284), the HII region associated with the open cluster Dolidze 25 using JWST/NIRCam with the aim to study star formation in a metal-poor, i.e., about 1/3 of solar, environment. In particular, we aim to measure the peak of the IMF. Using JWST NIRCam photometry, we identified the embedded cluster S284-EC1 and resolved its low-mass content. By comparison with pre-main sequence evolutionary tracks we determine the mass and extinction for the individual cluster members. Extinction limited samples are created based on the distribution of extinction and the completeness of the data. For the region with a completeness of 50% or higher, we have fitted a log-normal distribution to the IMF. Adopting a fiducial age of 1 Myr of the members, the peak of the IMF is found to be at mc = 0.16+-0.02Msun, which is significantly smaller than the peak mass measured in local young clusters, such as mc = 0.26+0.11-0.07 Msun in the Orion Nebula Cluster (Gennaro & Robberto 2020), or the local Galactic disk value of mc = 0.25 Msun (Chabrier 2005). We have found evidence for IMF variation as a function of metallicity, i.e., the peak of the IMF shifts to lower masses as one goes from solar to 1/3 solar metallicity. However, we caution that the result is sensitive to the assumed age of the stellar population, i.e., with peak mass rising if an age older than 1 Myr is adopted. This study further motivates the need for expanded samples of low-metallicity regions and their content to enable more comprehensive measures of the IMF in such environments.

The recent measurements from the Atacama Cosmology Telescope (ACT) favor a higher value of the scalar spectral index $n_s$ compared to the Planck data, challenging many well-established inflationary models. In this work, we investigate the viability of polynomial potential inflation in light of the latest ACT data, focusing on the quadratic ($n=2$) and cubic ($n=3$) models. By exploring the parameter space and deriving constraints on the model coefficients, we find that the cubic model can provide a good fit to the data, while the quadratic model struggles to simultaneously accommodate the ACT data and the requirement of sufficient inflation. Our findings show that higher-order polynomial potentials remain viable for describing cosmic inflation, demonstrating how recent precision measurements continue to refine our understanding of the early Universe.

The system TIC 470710327 is comprised of three main-sequence OB stars, with an inner compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. With the tertiary mass of 14.5-16 $M_{\odot}$ and both components in the inner eclipsing binary with individual masses of 6--7 and 5.5-6.3 $M_{\odot}$, it is currently the most massive compact system known. The formation scenario of such a compact triple is uncertain. It has been suggested that `2 + 2' quadruple dynamics can lead to a stellar merger in the initially more massive binary and finally result in a highly magnetised tertiary. Our study confirms the presence of a kG-order magnetic field in the tertiary and the slow rotation typical for massive magnetic stars. We conclude that finding massive merger candidates by studies of dynamics in compact, multiple-star systems is an efficient way to understand the evolution of massive stellar multiplicity and the generation of magnetic fields.

The early growth of high-redshift quasars and their host galaxies raises critical questions about their cosmic evolution. We exploit the angular resolution and sensitivity of NIRCam to investigate the host galaxies of 31 quasars at $4\lesssim z\lesssim7$ drawn from multiple JWST surveys. Using a new multi-band forward-modeling code (\textsc{GalfitS}) that incorporates physically motivated priors, we securely detect and quantify the host emission in 30 objects, while simultaneously characterizing the nuclear spectral energy distribution. The host galaxies of high-redshift quasars are $\sim 0.3$~dex more compact than star-forming galaxies of comparable mass. A striking dichotomy emerges: luminous ``blue'' quasars ($L_{5100}\gtrsim10^{45}\,{\rm erg\,s^{-1}}$) reside in bulge-dominated galaxies ($n \approx 5$) and exhibit a narrow range of ultraviolet nuclear slopes (median $\beta_{\rm UV} \approx -1.4$), while fainter ``red'' quasars inhabit disk-like hosts ($n\approx 1$) and display a broad range of slopes ($\beta_{\rm UV}\approx-2$ to 4). These two populations differ markedly in their black hole-to-stellar mass ratios, with high-luminosity quasars showing $M_{\mathrm{BH}}/M_\ast = 1.2\%$ compared to $4.7\%$ for lower luminosity sources, placing them collectively $\sim$0.6~dex above the local $M_{\mathrm{BH}}-M_\ast$ relation. This offset likely reflects rapid black hole growth in early gas-rich environments, where feedback from the active galactic nucleus becomes effective only after substantial gas depletion. Our findings suggest that the observed dichotomy, whether due to intrinsic spectral differences or dust extinction, fundamentally shapes the coevolution of supermassive black holes and their host galaxies in the early Universe.

A significant source of noise for Imaging Atmospheric Cherenkov Telescopes (IACTs), which are designed to measure air showers caused by astrophysical gamma rays, is optical light emitted from the night sky. This Night Sky Background (NSB) influences IACT operating times and their sensitivity. Thus, for scheduling observations and instrument simulation, an accurate estimate of the NSB is important. A physics-driven approach to simulating wavelength-dependent, per-photomultiplier-pixel NSB was developed. It includes contributions from scattered moonlight, starlight, diffuse galactic light, zodiacal light, and airglow emission. It also accounts for the absorption and scattering of optical light in the atmosphere and telescope-specific factors such as mirror reflectivity, photon detection efficiency, and focal length. The simulated results are corrected for pointing inaccuracies and individual pixel sensitivities and compared to data from the High Energy Stereoscopic System (H.E.S.S.) IACT array. The software package developed for this analysis will be made publicly available. Validation against H.E.S.S.\ data shows small deviations from the prediction, attributable to airglow and atmospheric variability. Per-Pixel predictions provide a good match to the data, with the relative 90\% error range being [-21\%, 19\%]. Compared to the existing standard modelling approach of assuming a constant background, where the relative 90\% error range was [-64\%, 48\%], this is a significant improvement.

These lecture notes provide an overview of high-energy astrophysical processes involving axions, axion-like particles (ALPs), and other weakly interacting slim particles (WISPs) focusing on their potential observational signatures in astrophysical environments. After introducing key concepts in high-energy astrophysics, we present the fundamental properties of WISPs, emphasizing their phenomenological implications. Particular attention is given to ALP-photon conversion in strong magnetic fields and the possible decay signatures of ALPs in sources such as active galactic nuclei, galaxy clusters, and cosmic-ray accelerators. These effects can lead to distinctive modifications in astrophysical spectra, spatial distributions, and polarization patterns, providing unique probes of physics beyond the Standard Model. We discuss their role in dark matter scenarios and their potential impact on high-energy observations. The lecture series is supplemented by hands-on tutorials, including exercises on axion electrodynamics and an analysis of gamma-ray data from NGC 1275 to search for ALP-photon conversion signatures.

Shocks in the solar corona can accelerate electrons that in turn generate radio emission known as type II radio bursts. The characteristics and morphology of these radio bursts in the dynamic spectrum reflect the evolution of the shock itself, together with the properties of the local corona where it propagates. In this work, we study the evolution of a complex type II radio burst showing a multi-lane structure, to find the locations where the radio emission is produced and relate them to the properties of the local environment. Using radio imaging, we track separately each lane composing the type II burst and relate the position of the emission to the properties of the ambient medium such as density, Alfven speed, and magnetic field. We show that the radio burst morphology in the dynamic spectrum changes with time and it is related to the complexity of the local environment. The initial stage of the radio emission show a single lane in the spectrum, while the latter stages of the radio signature evolve in a multi-lane scenario. The radio imaging reveals how the initial stage of the radio emission separates with time into different locations along the shock front as the density and orientation of the magnetic field change along the shock propagation. At the time where the spectrum shows a multi-lane shape, we found a clear separation of the imaged radio sources. By combining radio imaging with the properties of the local corona, we described the evolution of a type II radio burst and, for the first time, identified three distinct radio emission regions above the CME front. Two regions were located at the flanks, producing earlier radio emission than the central position, in accordance with the complexity of density and Alfven speed values in the regions where radio emission is generated. This unprecedented observation provides new insights into the nature of multi-lane type II bursts.

As the fields of stellar and exoplanetary study grow and revolutionary new detection instruments are created, it is imperative that a homogeneous, precise source of stellar parameters is available. This first work of the gr8stars collaboration presents the all-sky magnitude limited sample of 5645 bright FGKM dwarfs, along with homogeneously derived spectroscopic parameters of a subset of 1716 targets visible from the Northern hemisphere. We have collected high-resolution archival and new spectra from several instruments. Spectrosocpic parameters are determined using the PAWS pipeline, employing both the curve-of-growth equivalent width method, and the spectral synthesis method. We achieve median uncertainties of 106K in stellar effective temperature, 0.08 dex in surface gravity, and 0.03 dex in metallicity. This paper also presents photometric stellar parameters for these dwarfs, determined using SED fitting. The full gr8stars sample selection, including derived spectroscopic and photometric parameters, is made available through an interactive online database. We also perform a kinematic analysis to classify these stars according to their Galactic component.

While marginal in mass terms, dust grains play an outsized role in both the physics and observation of the interstellar medium (ISM). However, explicit modelling of this ISM constituent remains uncommon in large cosmological simulations. In this work, we present a model for the life-cycle of dust in the ISM that couples to the forthcoming COLIBRE galaxy formation model, which explicitly simulates the cold ISM. We follow 6 distinct grain types: 3 chemical species, including carbon and two silicate grains, with 2 size bins each. Our dust model accounts for seeding of grains from stellar ejecta, self-consistent element-by-element metal yields and depletion, grain size transfer and destruction of dust in the ISM. We detail the calibration of this model, particularly the use of a clumping factor, to account for unresolved gas clouds in which dust readily evolves. We present a fiducial run in a 25$^3$~cMpc$^3$ cosmological volume that displays good agreement with observations of the cosmic evolution of dust density, as well as the $z=0$ galaxy dust mass function and dust scaling relations. We highlight known tensions between observational datasets of the dust-to-gas ratio as a function of metallicity depending on which metallicity calibrator is used; our model favours higher-normalisation metallicity calibrators, which agree with the observations within 0.1~dex for stellar masses $>10^9 \; {\rm M_\odot}$. We compare the grain size distribution to observations of local galaxies, and find that our simulation suggests a higher concentration of small grains, associated with more diffuse ISM and the warm-neutral medium (WNM), which both play a key role in boosting H$_2$ content. Putting these results and modelling approaches in context, we set the stage for upcoming insights into the dusty ISM of galaxies using the COLIBRE simulations.

Terrestrial planets within the Venus zone surrounding M dwarf stars can retain surface ice caps on the perpetual dark side if atmospheric heat transport is inefficient, {as suggested by previous global climate simulations \citep[e.g.,][]{leconte2013}.} This condition is {proposed} to play a role in the potential regional habitability of these planets. However, the amount of surface ice may be limited by considering the water condensed from the steam atmosphere in a runaway greenhouse state, and the physical mechanism for triggering the condensation process is not clear. Here, we use a two-column moist radiative-convective-subsiding model to investigate the water condensation process on tidally locked planets from the runaway greenhouse state. We find that the water condensation process is characterized by two distinct equilibrium states under the same {incoming stellar flux}. The initiation of condensation corresponds to a warm, unstable state exhibiting positive Planck feedback, whereas the termination phase corresponds to a cold, stable state exhibiting negative Planck feedback. We further show that the surface water mass in the collapsed state {decreases with the} incoming stellar flux, background surface pressure, and optical thickness of non-condensible greenhouse gases, with a global equivalent depth of less than $\sim 20$ cm. Our two-column approach provides a straightforward way to understand the water evolution on Venus zone planets around M dwarfs.

Seismic isolation is crucial for gravitational wave detectors as it minimizes ground vibrations, enabling the detection of faint gravitational wave signals. An active seismic isolation platform for precision measurement experiments is described. The table features piezo actuation along five degrees of freedom: three translational actuations and two tip-tilt degrees of freedom along the horizontal axes. It is stiff in rotation about the vertical axes. A seismometer is used to sense table motion. Piezo actuators are used to suppress seismic noise with feedback control bandwidth of 0.3 to 3 Hz. Suppression levels ranging from 21 to 36 dB of seismic noise within the frequency range of 0.5 to 1.3 Hz are demonstrated, as measured by a witness seismometer on the table, with the suppression direction along the axis of the longitudinal translation of the suspended mirror on the table. The suppression results in 1 $\mathrm{\mathrm{nm/\sqrt{Hz}}}$ residual horizontal motion at 1 Hz. Limitations such as tilt-to-translation coupling that prevent actuation over the desired range of 0.03 to 3 Hz are discussed.

Emission from within the plunging region of black hole accretion flows has recently been detected in two X-ray binary systems. There is, furthermore, a possible discrepancy between the inferred spins of gravitational wave and electromagnetically detected black holes. Motivated by these two results we demonstrate, using theoretical calculations, numerical simulations and observational data, that the inclusion of emission from within the innermost stable circular orbit (ISCO) results in a black hole with a low spin producing a thermal continuum X-ray spectrum that mimics that produced by a much more rapidly rotating black hole surrounded by a disk with no emission from within the ISCO. We demonstrate this explicitly using the observed X-ray spectrum of a canonical soft-state high mass X-ray binary system M33 X-7. A vanishing ISCO temperature model requires a high spin $a_\bullet = 0.84\pm0.05$, as has been found previously in the literature. However, a disk around a Schwarzschild black hole can equally well (in fact slightly better) describe the data, provided that photons emitted from within the plunging region are included, and the ISCO stress is in line with that seen in numerical simulations of the accretion process. We then present an analysis of two further soft-state X-ray binaries (MAXI J1820+070 and MAXI J0637$-$430) which require the presence of intra-ISCO emission at high statistical significance. These two sources sit on the low-spin moderate-stress part of the degeneracy exhibited by M33 X-7, suggesting that when high quality data are available the high-spin low-stress region of parameter space is ruled out. We discuss how future advances in numerical simulations and data modelling will be essential to determining the spin of X-ray binary black holes which may well be systematically lower than current continuum fitting methods suggest.

The Atacama Large Millimeter Array (ALMA) has revolutionised the field of solar millimetre astronomy with its high angular resolution and cadence. However, with a limited field of view (FOV), targeted observations of highly dynamic phenomena such of flares are challenging. A large aperture single-dish telescope with a large FOV, such as the future Atacama Large Aperture Submillimeter Telescope (AtLAST), would prove useful in observing such phenomena, as one could scan the full solar disk on shorter timescales. We aimed to explore what FOVs, detector counts, and scan strategies are suitable for AtLAST to push the required full-disk scan times below 1 minute, enabling regular observations of dynamic solar phenomena. Utilising the maria code, we were able to simulate solar observations with AtLAST, and thoroughly explored how instrumental properties and scanning strategies affect the full-disk observations in the planned frequency bands. We find the double-circle scan pattern, currently employed at ALMA for full-disk mapping to also be an acceptable way of scanning the Sun with AtLAST. Using small to intermediately sized instruments (1000 - 50,000 detector elements), the estimated observational cadence would be less than 1 minute across AtLAST's frequency range with a reasonable pixel spacing. Using instruments with larger FOVs ($\gtrapprox 0.25^\circ$, equivalent to $\gtrapprox$ 1 R$_\odot$), we find a simple circular scan to be more efficient, achieving cadences on second time scales, but requiring more detector elements ($\gtrapprox$ 100,000). We find that a large FOV single-dish telescope such as AtLAST could provide the solar millimetre community with hitherto unachievable observations, namely full-disk observations at high cadence and adequate resolution. With cadences potentially down to seconds, such an instrument would be ideal in the study of quickly evolving solar phenomena.

We derive the star formation history (SFH) and chemical evolution history (CEH) of the Ursa Minor (UMi) dwarf spheroidal galaxy (dSph). We detect two distinct stellar populations that exist over 6 times half-light radius from its center. The results are obtained by applying a newly developed algorithm to the deep and wide-field photometric dataset taken with Hyper Suprime-Cam on the Subaru Telescope. The algorithm employs the genetic algorithm and the simulated annealing to minimize a $\chi^{2}$ value between the observed color-magnitude diagram (CMD) and synthetic CMD generated from the stellar isochrones. The age and metallicity resolutions are set to $0.5$ Gyr and $0.1$ dex, respectively. The accuracy assessment with mock galaxies shows that it returns the peaks of metallicity distributions and star formation period within 1 $\sigma$ of input value in the case of a single population. In tests with two populations, two distinct metallicity peaks are identified without an offset from the input values, indicating the robustness of this algorithm. The detected two populations in the UMi dSph have the metallicity peaks of [Fe/H] = $-2.2$ and $-2.5$; the metal-rich population started its star formation about 1 Gyr later than the metal-poor one. The SFH of both metal-rich and metal-poor populations varies with distance from the center of the UMi dSph, without any age-gradients. These results suggest that the UMi dSph underwent a complex formation process, contrary to the simple formation history of dwarf galaxies previously thought.

Spatially resolved maps of stellar populations and nebular emission are key tools for understanding the physical properties and evolutionary stages of galaxies. We aim to characterize the spatially resolved stellar population and emission line properties of galaxies in the M101 group using Javalambre Photometric Local Universe Survey (J-PLUS) data. The datacubes first go through pre-processing steps, which include masking, noise suppression, PSF homogenization, and spatial binning. The improved data are then analyzed with the spectral synthesis code \alstar, which has been previously shown to produce excellent results with the unique 12 bands filter system of J-PLUS and S-PLUS. We produce maps of stellar mass surface density ($\Sigma_\star$), mean stellar age and metallicity, star formation rate surface density ($\Sigma_{\rm SFR}$), dust attenuation, and emission line properties such as fluxes and equivalent widths of the main optical lines. Relations among these properties are explored. All galaxies exhibit a well-defined age-$\Sigma_\star$ relation, except for the dwarfs. Similarly, all of the galaxies follow local $\Sigma_\star$-$\Sigma_{\rm SFR}$ star-forming MS relations, with specific star formation rates that grow for less massive systems. A stellar $\Sigma_\star$-metallicity relation is clearly present in M101, while other galaxies have either flatter or undefined relations. Nebular metallicities correlate with $\Sigma_\star$ for all galaxies. This study demonstrates the ability of J-PLUS to perform IFS-like analysis of galaxies, offering robust spatially resolved measurements of stellar populations and emission lines over large fields of view. The M101 group analysis showcases the potential for expanding such studies to other groups and clusters, contributing to the understanding of galaxy evolution across different environments.

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pre-tabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes, these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus' instability criteria for TI can be used to define a critical cooling rate, $\Lambda_c$, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore extend Balbus' criteria to account for thermal conduction. Upon applying a $\Lambda_c$-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

The nature of the Diffuse Interstellar Band (DIB) carriers is perhaps the most studied, longest standing, unresolved problem in astronomy. While four bands have been associated with the fullerene cation (C^+_60) the vast majority (> 550) remain unidentified. This works is an attempt to provide a conceptual framework for the typical energy transitions that are central to explaining the origin of DIBs, however, it does not make an association between these transitions and any particular DIBs. The effect of quantum confinement on excitons is used, including charge transfer excitons, to construct a generic basis for the electronic transitions that could, in principle, be coherent with the energies associated with DIBs. In this model the carriers are carbon nanodots (CNDs) modelled as nanodiamonds and a-C(:H) nanoparticles. These preliminary results seem to show that particle size dependent effects in nanodiamond and a-C(:H) CNDs could be consistent with the positions of, and intervals between, some of the DIBs. One particular strength of the model is that predicts single bands from the majority of single-size particles and, at most, two bands from some of these same carriers. In the latter case the two bands come from different transitions and may or may not inter-correlate, depending upon the local environment. This generic framework indicates that the size dependent fundamental transitions in CNDs could provide a viable scenario for the origin of some DIB-type bands. While this work does not identify a single DIB, it furnishes a conceptual view for the DIB origin, and suggests that a more refined exploration of quantum confinement size effects, and exciton physics within the astronomical domain might prove fruitful. This work also hints at the requirement for stable configurations for particular size domains in order to explain DIB wavelength stability.

The formation of the giant binary-like system composed by the Corona Borealis and Abell 2142 superclusters is an intriguing conundrum of the formation of large scale structures since, from the observational point of view, it represents a rare peculiarity in the distribution of massive galaxy superclusters. Having a configuration similar to a giant binary system interconnected by a huge filament, it is likely one if not the unique case to date observed in the Local Universe. So a question arise: how and when did it form? Here, the CrB and A2142 system has been hypothesized to be a descendent of a primordial binary system composed of two close galaxy proto-clusters, the CrB and A2142, originated from two independent collapsing process within a dense cloud of galaxies. Then, assuming that at a certain point the decoupling occurred due to the interplay of gravity-antigravity of DE, they began to move away in radial motion from each other following the accelerate expansion of the Universe. In the context of the LCDM model a Newtonian approximation of the two body motion in presence of DE has been applied using current physical parameters of the CrB and A2142 system with the aim to calculate the look-back time at decoupling. Its compatibility with the era of formation of the primordial galaxy proto-clusters has been tested. The event of the binary decoupling happened in a time larger than the age of the Universe! Of course, after such a result the advanced hypothesis on the origins of the CrB-A2142 system can be rejected, but not the test of compatibility since the look-back time at the error lower limit is 12.5 Gys ( z = 4.5) largely consistent with observations.

The origin of compact millimeter (mm) continuum emission from radio-quiet AGNs (RQAGNs) is still not fully understood. Changing-state AGNs (CSAGNs) display rapid and strong variability, which can allow us to investigate the origin of the mm emission. We present here the results of the first study of the mm continuum variability of a CSAGN using archival ALMA band 6 ($\sim 230$ GHz) observations of NGC 1566 obtained in 2014-2023. We find a positive correlation between the mm and X-ray flux with an intrinsic scatter of 0.05 dex ($1\sigma$), suggesting a common origin. The mm spectral index ($\alpha_{\rm mm}$) is found in the range of $0.13\pm0.38$ to $-0.26\pm0.53$, consistent with a compact optically thick synchrotron source. No significant correlation was found between the $\alpha_{\rm mm}$ and the mm flux. The mm/X-ray ratio also shows no clear link to the Eddington ratio but is higher in the low-accretion state. We discuss several scenarios about the origin of the mm emission in NGC 1566. We find that synchrotron emission in the magnetized X-ray corona appears to be the most probable origin of mm emission, confirming that mm emission can be used as a tracer of AGN activity in RQAGNs.

This study investigates the radio spectral properties of \textit{K}$_{S}$-selected star-forming galaxies (SFGs) in the XMM-LSS field using extensive multiwavelength data. By employing various diagnostics, SFGs are distinguished from quiescent galaxies and AGN across seven redshift bins ($\rm{0.1\leq\,\textit{z}\,\leq\,3.0}$). The broadband radio frequency spectral energy distribution is analysed at observer-frame frequencies from 144 to 1500 MHz using median stacking techniques correcting for median flux boosting. We investigate the relationship between the radio spectral index, $\alpha$ (where $S\propto\nu^{\alpha}$) and redshift ($z$). Our analysis reveals no significant inverse correlation between $\alpha$ and $z$, indicating that the radio spectrum remains independent with varying redshift. We fit the stacked median radio SEDs with a power law (\textit{PL}), curved power law (\textit{CPL}) and double power law (\textit{DPL}) models. For the \textit{DPL} and \textit{CPL} models, we observe a consistent steepening of the low-frequency spectral index across all redshift bins. For the \textit{CPL} model, the curvature term $q$ is greater than zero in all redshift bins. Model comparisons indicate that spectra are generally well fitted by all the models considered. At 1500 MHz, SFGs display both a steep synchrotron component and a flat free-free emission component, with a thermal fraction consistently around 11$\%$ to 18$\%$. Further deep radio observations, with higher resolution to better deal with source blending and confusion noise and wider frequency coverage to better separate non-thermal and thermal radio emission, are required to reveal the detailed physical processes, thus clarifying the nature of radio sources.

Measurement of the growth rate of structures ($\fsig$) with Type Ia supernovae (\sns) will improve our understanding of the nature of dark energy and enable tests of general relativity. In this paper, we generate simulations of the 10 year \sn\ dataset of the Rubin-LSST survey, including a correlated velocity field from a N-body simulation and realistic models of \sns\ properties and their correlations with host-galaxy properties. We find, similar to SN~Ia analyses that constrain the dark energy equation-of-state parameters $w_0w_a$, that constraints on $\fsig$ can be biased depending on the intrinsic scatter of \sns. While for the majority of intrinsic scatter models we recover $\fsig$ with a precision of $\sim13 - 14\%$, for the most realistic dust-based model, we find that the presence of non-Gaussianities in Hubble diagram residuals leads to a bias on $\fsig$ of about $\sim-20\%$. When trying to correct for the dust-based intrinsic scatter, we find that the propagation of the uncertainty on the model parameters does not significantly increase the error on $\fsig$. We also find that while the main component of the error budget of $\fsig$ is the statistical uncertainty ($>75\%$ of the total error budget), the systematic error budget is dominated by the uncertainty on the damping parameter, $\sigma_u$, that gives an empirical description of the effect of redshift space distortions on the velocity power spectrum. Our results motivate a search for new methods to correct for the non-Gaussian distribution of the Hubble diagram residuals, as well as an improved modeling of the damping parameter.

Directly imaged exoplanets in wide orbits challenge current gas giant formation theories. They need to form quickly and acquire enough material before the disk dissipates, which cannot be accommodated by in-situ formation by core accretion. We search for wide separation ($>$ 100 au) planetary-mass companions with the Young Suns Exoplanet Survey (YSES). Here, we present a planetary-mass candidate companion discovered in the survey. We conducted follow-up observations of the candidate system after the first epoch observations and obtained six epochs of observations for this system between 2018 and 2024, and integral field spectroscopy of the stellar component. We report the detection of a candidate companion with H=22.04 $\pm$ 0.13 mag at a projected separation of 730 $\pm$ 10 au away from the primary star. High angular resolution imaging observations of the central star show it is a visual binary. Acceleration data, orbital fitting, spectral energy distribution fitting and radial velocity differences all suggest that there is at least one more unresolved low-mass stellar companion in this system. The planetary-mass candidate shows a significant proper motion comparable to that of the primary star. We estimate an age of 19-28 Myr for the primary star. We cannot confirm the companionship of the candidate due to the unknown barycentre of the stars. Long-term imaging and radial velocity monitoring of the central stars, along with spectroscopy of the candidate companion, are key to resolving the nature of this system. If confirmed, the candidate companion would have a mass of 3-5 Mj estimated with the ATMO evolutionary model. It would be another cold low-mass planet imaged similar to 51 Eri b and AF Lep b. Its extremely wide separation from the host star would challenge the formation theory of gas giant exoplanets.

The knowledge of the point spread function (PSF) of the Mikhail Pavlinsky Astronomical Roentgen Telescope - X-ray Concentrator (ART-XC) telescope aboard the Spectrum-Roentgen-Gamma (SRG) observatory plays an especially crucial role in the detection of point X-ray sources in the all-sky survey and the studies of extended X-ray objects with low surface brightness. In this work, we calibrate the far off-axis shape of the ART-XC PSF using in-flight data of Sco X-1 and the Crab Nebula, in all-sky survey or scan mode, respectively. We demonstrate that the so-called "slewing" ART-XC PSF (in contrast to the on-axis PSF), in convolution with the detector pixels, is consistent with ground calibration performed at the Marshall Space Flight Center, and can be used to model the PSF up to large off-axis distances in all-sky survey or scan modes. The radial profile of the Crab Nebula in the 4-12 keV band shows an extended structure out to ~150" and is consistent with Sco X-1 at larger off-axis angles. Finally, we performed an analytic parametrization of the slewing ART-XC PSF as a function of energy.

Among the thousands of observed fast radio bursts (FRBs), a few sources exhibit exceptionally high burst activity observable by many telescopes across a broad range of radio frequencies. Almost all of these highly active repeaters have been discovered by CHIME/FRB, due to its daily observations of the entire Northern sky as a transit radio telescope. FRB 20240114A is a source discovered and reported by CHIME/FRB to the community in January 2024; given its low declination, even the detection of a few bursts hints at a high burst rate. Following the community announcement of this source as a potentially active repeater, it was extensively followed up by other observatories and has emerged as one of the most prolific FRB repeaters ever observed. This paper presents the five bursts CHIME/FRB observed from FRB 20240114A, with channelized raw voltage data saved for two bursts. We do not observe changes in the DM of the source greater than ~1.3 pc cm$^{-3}$ in our observations over nearly a year baseline. We find an RM of ~ +320 rad m$^{-2}$. We do not find evidence for scattering at the level of < 0.3 ms in the bursts, and we find no evidence for astrophysical scintillation. In our observations of FRB 20240114A, we see a burst rate ~49x higher than the median burst rate of apparent non-repeaters also discovered by CHIME/FRB. Each discovery of highly active FRBs provides a valuable opportunity to investigate whether there is a fundamental difference between repeating and apparently non-repeating sources.

The $\Gamma$ growth model is an effective parameterization employed across various scientific disciplines and scales to depict growth. It has been demonstrated that the cosmic star formation rate density (CSFRD) can also be described broadly by this pattern, i.e. $\frac{dM(T)}{dT} = M_{z,0}\, \times \frac{\beta^{\alpha}}{\Gamma(\alpha)} \, T^{\alpha-1} e^{-\beta \, T }$ M$_{\odot}$ Gyr$^{-1}$, where $M_{z,0}$ is the stellar mass at $z$ = 0, $\alpha = 3.0$, $\beta = 0.5 $ Gyr$^{-1}$ and $T$ describes time. We use the identical $\Gamma$ growth pattern given by the CSFRD to extend the present day (z = 0) stellar mass bins $M_{\ast}(T)$ of the Galaxy Stellar Mass Function (GSMF) and investigate if we are able to reproduce observations for the high redshift GSMFs. Surprisingly, our scheme describes successfully the evolution of the GSMF over 13.5 Gyrs, especially for objects with intermediate and low masses. We observe some deviations that manifest {\it solely} at very high redshifts ($z > 1.5$, i.e. more than 9.5 Gyr ago) and {\it specifically} for very small and exceedingly massive objects. We discuss the possible solutions (e.g. impacts of mergers) for these offsets. Our formalism suggests that the evolution of the GSMF is set by simple (few parameters) and physically motivated arguments. The parameters $\beta$ and $\alpha$ are theoretically consistent within a multi-scale context and are determined from the dynamical time scale ($\beta$) and the radial distribution of the accreting matter ($\alpha$). We demonstrate that both our formalism and state-of-the-art simulations are consistent with recent GSMFs derived from JWST data at high redshifts.

Understanding the origin and evolution of magnetic fields on cosmological scales opens up a window into the physics of the early Universe. Numerical simulations of such fields require a careful treatment to faithfully solve the equations of magnetohydrodynamics (MHD) without introducing numerical artefacts. In this paper, we study the growth of the magnetic fields in controlled kinematic dynamo setups using both smoothed particle hydrodynamics implementations in the SWIFT code. We assess the quality of the reconstructed solution in the Roberts flow case against the reference implementation in the Pencil code and find generally a good agreement. Similarly, we reproduce the known features of the more complex ABC flow. Using a simple induction-diffusion balance model to analyse the results, we construct an "overwinding" trigger metric to locally detect regions where the magnetic diffusion cannot counteract the expected induction because of limitations in the method's ability to resolve magnetic field gradients. This metric is then used to identify the necessary resolution and resistivity levels to counteract the overwinding problem. We finally apply this metric to adiabatic cosmological simulations and discuss the resolution requirements needed to resolve the growth of the primordial fields without artefacts.

Although most white dwarf stars have hydrogen-dominated atmospheres, a significant fraction have atmospheres in which hydrogen is spectroscopically absent, with the fraction of hydrogen-free atmospheres varying with effective temperature. Estimates of the total mass of hydrogen, $M_H$, in the stellar envelope from either asteroseismology or spectral evolution are at odds with predicted values from theoretical stellar evolution modeling. Recent work has found that models in the early post AGB phase of evolution can exhibit thermally and dynamical unstable behavior. Here we investigate whether this Early Post AGB Instability (EPAGBI) can help resolve the conflict in $M_H$ values determined from white dwarf spectral evolution, analysis of DAV pulsations and canonical stellar evolution modeling, by evolving solar composition models of mass $1$ and $2 M_\odot$ through the AGB phases and to the white dwarf cooling track. The $M_H$ values at the end of the calculations are in the range consistent with asteroseismological determinations. However, we caution that, because hydrodynamic behavior is not included in our modeling, it is possible that all hydrogen would be removed. Thus, it is unclear whether the occurrence of the EPAGBI resolves the discrepancy between predictions of stellar evolution modeling and the asteroseismological hydrogen envelope mass determinations. The major impact of EPAGBIs is that they cause loops in the HRD. For models of AGB-departure mass 0.567 and 0.642 $M_\odot$, it takes approximately 100 and 10 yr for a single HRD loop, respectively. Such loops might be detectable in a long-term monitoring program, or perhaps by their imprint on planetary nebula morphology imparted by the cyclically varying mass loss rate.

During the early stages of star formation, accretion processes such as infall from the envelope and molecular streamers, and ejection of matter through winds and jets take place simultaneously. The Class 0/I binary [BHB2007] 11 shows evidence for accretion and ejection at the scales of the circumbinary disk and the inner close binary. Recent H$_2$CO observations showed two elongated structures with hints of outflowing motion almost perpendicular to the main CO outflow, which is launched from the circumbinary disk. With the aim of verifying the nature of these elongated structures, we analyze the line emission of H$^{13}$CO$^+$, CCH, c-C$_3$H$_2$ and SiO observed with ALMA within the Large Program FAUST. These molecules trace material moving at velocities close to the ambient cloud velocity. The images of H$^{13}$CO$^+$, CCH, c-C$_3$H$_2$ show the elongated structures, whose gas kinematics are consistent with outflowing motions and with rotation in the opposite sense to the main CO outflow. The derived mass loss rate from these large-scale structures is $(1.8\pm0.5)\times10^{-6}M_{\odot}\textrm{ yr}^{-1}$, in agreement with those measured in outflows driven by Class 0/I protostars. The SiO image reveals compact emission close to the binary system, with a slight elongation aligned with the larger-scale structures. This suggests that SiO is released from the sputtering of dust grains in the shocked material at the base of the potential new outflow, with a relative abundance of $\geq(0.11-2.0)\times10^{-9}$. However, higher angular and spectral resolution observations are needed to accurately estimate the outflow launching radius and its powering source. Given the location and the abundance of the SiO emission, we propose that the second outflow may be launched from inside the circumbinary disk, likely by the less massive companion, which is actively accreting material from its surroundings.

Polarization behaviour shows a transition in the pulsar population, where energetic sources with higher spin-down energy loss, $\dot{E} > 10^{34}$ erg~s$^{-1}$, often have fractional linear polarisation ($L/I$) close to 100\%, while below this range $L/I$ is usually lower than 50\%. The polarisation behaviour has been primarily studied at higher frequencies above 1 GHz, and in this work we explore the single pulse polarisation behaviour in pulsars with $\dot{E} > 5\times10^{33}$ erg~s$^{-1}$ at a lower frequency range of 300-750 MHz. The polarisation behaviour can be divided into two categories, the first with $L/I>$ 70\% where the polarisation position angle (PPA) follows a single track, and a second group with $L/I <$ 70\% and scattered PPA behaviour with or without orthogonal modes. However, there are some single pulses in the first category that also have lower $L/I$ and exhibit the presence of two polarisation modes along orthogonal tracks. The radio emission in pulsars arises due to coherent curvature radiation (CCR) from charge bunches, which develops due to non-linear instabilities in the pulsar plasma forming charge separated envelope solitons. The CCR excites orthogonally polarised X and O modes oriented perpendicular and parallel to the magnetic field line planes, that detach in the plasma and propagate independently. The O-mode is seven times stronger than the X-mode but gets damped in the medium. We show that incoherent mixing of the X and O modes with different levels of damping can reproduce the observed polarisation features in the energetic pulsar population.

As solar coronal mass ejections (CMEs) propagate through the heliosphere, they expend energy in heating protons to compensate for the cooling that occurs due to expansion. CME propagation models usually treat energy dissipation implicitly via a polytropic index ($\delta$). Here we calculate the power dissipation implied by a given $\delta$ and compare it with the power available in the turbulent velocity fluctuations. We make this comparison using near-Earth {\em in-situ} observations of 27 of the most geoeffective CMEs ($D_{\rm st} < -75$ nT) in solar cycle 24. For $\delta = 5/3$, the power in the turbulent velocity fluctuations is $\approx 54$\% smaller than what would be required to maintain the proton temperature at the observed values. If the power in the turbulent cascade is assumed to be fully expended in local proton heating, the most probable value for $\delta$ is 1.35. Our results contribute to a better understanding of CME energetics, and thereby to improved CME propagation models and estimates of Earth arrival times.

We present a new method for estimating galaxy cluster masses using weak-lensing magnification. The effect of weak-lensing magnification introduces a correlation between the position of foreground galaxy clusters and the density of background sources. Therefore, cluster masses can be inferred through observations of these correlations. In this work, we introduce a method that allows us to considerably reduce noise correlations between different radial bins of the cluster magnification signal via a Wiener filtering of our observed magnification field on large scales. This method can reduce the uncertainty on the estimated galaxy cluster mass and it can also be applied to cluster mass estimation for weak-lensing shear. The method was applied to Hyper-Suprime Cam galaxies and CAMIRA clusters detected within the Hyper-Suprime Cam survey (HSC). With HSC data, we find that our filtering method significantly reduces the correlation of noise between radial magnification bins. The estimated cluster mass is consistent between the filtered and unfiltered methods, with similar errors between the two methods as our current measurement errors contain significant contributions from the irreducible shot-noise. For deeper surveys, the effects of shot noise will be less important and this method will lead to greater improvements on the estimated cluster mass.

Recent JWST observations of the temperate sub-Neptune K2-18 b have been interpreted as suggestive of a liquid water ocean with possible biological activity. Signatures of DMS and DMDS have been claimed in the near-infrared (using the NIRISS and NIRSpec instruments) and mid-infrared (using MIRI). However, the statistical significance of the atmospheric imprints of these potential biomarkers has yet to be quantified from a joint analysis of the entire planet spectrum. We test the robustness of the proposed DMS/DMDS detections by simultaneously modeling the NIRISS and NIRSpec observations jointly with the MIRI spectrum, considering different data reductions and modeling choices. We use three well-tested pipelines to re-reduce the JWST observations, and two retrieval codes to analyze the resulting transmission spectra as well as previously published data. Our joint analysis of the panchromatic (0.6 - 12 um) spectrum of K2-18 b finds insufficient evidence for the presence of DMS and/or DMDS in the atmosphere of the planet. Furthermore, other molecules containing methyl functional groups (e.g., ethane) with absorption bands similar to DMS/DMDS provide an equally good fit to the data. We find that any marginal preferences are the result of limiting the number of molecules considered in the model and oversensitivity to small changes between data reductions. Our results confirm that there is no statistical significance for DMS or DMDS in K2-18 b's atmosphere. While previous works have demonstrated this on MIRI or NIRISS/NIRSpec observations alone, our analysis of the full transmission spectrum does not support claims of potential biomarkers. Using the best-fitting model including DMS/DMDS on the published data, we estimate that ~25 more MIRI transits would be needed for a 3-sigma rejection of a flat line relative to DMS/DMDS features in the planet's mid-infrared transmission spectrum.

Primordial scalar perturbations that reenter the horizon after inflation may induce a second-order Gravitational Wave spectrum with information about the primordial Universe on scales inaccessible to Cosmic Microwave Background experiments. In this work, we develop a general framework for the study of Scalar-Induced Gravitational Waves in Palatini $f(R)$ gravity that was proven to successfully realise inflation and quintessence, and consider the case of the Starobinsky-like model as an example. A regime of radiation domination with a subdominant matter component is assumed, allowing for a well-motivated perturbative approach to the modified gravity corrections. We calculate the kernel function and the density spectrum numerically and find accurate analytical expressions. The spectral density, which may be tested across a wide range of frequencies by upcoming Gravitational Wave experiments, is shown to differ from the General Relativity and metric $f(R)$ gravity predictions under certain conditions. We comment on previous results in the literature regarding the metric formulation and make special emphasis on the potential of these distinctive features of the spectrum to probe the two formalisms of gravity.

Oscillations in the primordial bispectrum are sourced by a range of inflationary phenomena, including features in the inflaton potential and interactions with massive fields through the Cosmological Collider scenario. These signatures offer a powerful window into early-universe physics. In this work, we study how oscillations of the form $\lim_{q\ll k}B(q,k)\propto \cos(\mu \ln(q/k))$ impact the non-linear squeezed matter bispectrum. Using a suite of $N$-body simulations with non-Gaussian initial conditions, we show that non-linear evolution significantly damps these oscillations, effectively erasing the signal on scales $k \gtrsim 0.3~h/{\rm Mpc}$ at redshift $z=0$. This damping is well-described by the Zel'dovich approximation and can be modeled deep into the non-linear regime using non-perturbative separate universe simulations. Promisingly, we show that reconstruction techniques developed for baryon acoustic oscillation (BAO) analyses can largely undo this damping, improving constraints on the amplitude (phase) of oscillations in the primordial squeezed bispectrum by up to a factor of five (four) at $z=0$. We also discuss several challenges with modeling the non-linear evolution of the squeezed bispectrum in the Cosmological Collider scenario, where the bispectrum is suppressed by a factor of $(q/k)^{3/2}$ relative to the template studied here. Our findings pave the way for future searches for oscillatory bispectra using large-scale structure data.

We investigate the impact of curvature corrections to Starobinsky inflation in light of the latest observational results from the Atacama Cosmology Telescope (ACT). While the pure Starobinsky model remains a compelling candidate for cosmic inflation, we explore how the higher-order curvature terms $R^3$, $R^{4}$ and $R^{3/2}$ modify the inflationary predictions. Using the scalar-tensor formulation of $f(R)$ gravity, we derive the effective scalar potentials and compute the resulting scalar tilt $n_{s}$ and tensor-to-scalar ratio $r$. We show that those curvature corrections can shift the predictions to align better with the ACT data, thus providing a possible resolution to a minor discrepancy between the standard Starobinsky model and ACT observations. Our findings suggest that the modified Starobinsky models with the higher-order curvature terms offer a viable pathway to reconciling inflationary predictions with precision cosmological measurements. At the same time, measuring or constraining primordial tensor modes can help to discriminate between these corrections.

We report the existence of horizonless compact object solutions supported by dust in the Minimal Exponential Measure (MEMe) model, a theory which modifies the couplings between gravity and matter without introducing dynamical degrees of freedom. For a perfect fluid source, the field equations for the MEMe model can be rewritten as the Einstein field equations sourced by a perfect fluid with a transformed equation of state, which can endow a sufficiently dense cloud of dust with an effective pressure. The resulting dust-supported horizonless compact objects can have masses below $\sim 10^{-11}~M_\odot$, making them suitable as MACHOs comprising a significant mass fraction for dark matter. A necessary condition for the existence of these compact object solutions is that the single free parameter in the MEMe model is positive-valued. Additionally, we find that this positive sign for the parameter can provide a mechanism for suppressing the formation of (primordial) black holes from the gravitational collapse of matter below a certain mass scale.

Correlated noise sources, particularly magnetic noise, form a risk to future gravitational-wave searches aimed at detecting the gravitational-wave background. Potential noise contamination is investigated by making noise projections which typically rely on an accurate measurement of the coupling strength of the noise to the detector. To make these projections, we inject, for the first time, broadband, coherent magnetic noise between two gravitational-wave detectors, LIGO Hanford and LIGO Livingston, separated by several thousands of kilometers. We describe the noise injection as well as its impact on the analysis pipelines and investigate the accuracy of noise projection techniques used in the past decade. Finally, we present a proof-of-concept demonstration of noise subtraction using Wiener filtering, while also highlighting potential risks associated with this method. This unique data set with correlated noise caused by magnetic field fluctuations in two gravitational-wave detectors, as well as in an array of witness sensors, provides an excellent testing ground for additional future studies. Ultimately, this study demonstrates that Wiener filtering is effective and can be applied in the eventual detection of the gravitational-wave background by the LIGO-Virgo-KAGRA Collaboration.

The symmetry energy expansion is a useful way to parameterize the properties of dense matter near nuclear saturation density, and much work has been done to connect physical quantities like the neutron star radius and the core-crust transition density to the symmetry energy parameters. In this work, I connect the weak-interaction-driven bulk viscosity in neutron-proton-electron ($npe$) matter to the symmetry parameters by calculating the susceptibilities of dense matter in terms of the symmetry energy. I use this result to calculate the resonant-peak value of the bulk viscosity as a function of density, finding that it strongly depends on $L$, as does the minimum bulk-viscous dissipation timescale. Also resulting from this calculation is a formula for finding the conformal points of the zero-temperature equation of state. Finally, I determine for which values of the symmetry parameters the maximum r-mode stable rotation frequency of an $npe$-matter neutron star is smaller than the Kepler frequency, in the high-temperature conditions where bulk viscosity is the dominant dissipation mechanism.

A plasma void forms downstream of the Moon when the solar wind impacts the lunar surface. This void gradually refills as the solar wind passes by, forming the lunar wake. We investigate this refilling process using a fully kinetic particle-in-cell (PIC) simulation. The early stage of refilling follows plasma-vacuum interaction theory, characterized by exponential decay of plasma density into the wake, along with ion acceleration and cooling in the expansion direction. Our PIC simulation confirms these theoretical predictions. In the next stage of the refilling process, the counter-streaming supersonic ion beams collide, generating Debye-scale electrostatic shocks at the wake's center. These shocks decelerate and thermalize the ion beams while heating electrons into flat-top velocity distributions along magnetic field lines. Additionally, fast magnetosonic waves undergo convective growth via anomalous cyclotron resonance as they co-propagate with temperature-anisotropic ion beams toward the wake's center. Electromagnetic ion cyclotron waves may also be excited through normal cyclotron resonance, counter-propagating with these anisotropic ion beams. Our findings provide new insights into the kinetic aspects of lunar wake refilling and may enhance interpretation of spacecraft observations.

On November 6, 2024, the Parker Solar Probe flew past Venus to make the first accurate electric field measurement in the nightside Venusian magnetosphere. To achieve this result, the electric field antennas were current biased in a way never before experienced by an electric field detector. This biasing requirement, that the positive bias current in the Venus shadow be about equal to the electron thermal current, is discussed and illustrated. About one minute of useful electric f ield data in the eight-minute nightside magnetosphere crossing was obtained, during which the only feature observed was a few Hz signal. This result, along with the magnetic field measurements, showed that there were few if any electromagnetic waves, such as low-frequency electromagnetic turbulence or whistlers, in the nightside crossing. Instead, a few Hertz, purely electrostatic signal was found. This suggests that the interaction of the solar wind with an unmagnetized body having an ionosphere may be different from that of previously studied magnetized bodies. In the sunlit flanks, many electromagnetic wave modes were observed. These results describe the first step in the proper technique for future measurements of electric fields in shadow.

Hierarchical granular piles composed of aggregates are key structural features in both geoscience and planetary science, from fault gouge in seismic zones to the internal structures of comets. Although experimental studies have suggested a multi-step evolution in their packing structure, this hypothesis has lacked numerical validation. In this study, we performed large-scale numerical simulations using the discrete element method to investigate the compressive behavior of hierarchical granular piles. We successfully reproduced and confirmed a three-stage evolution process: (i) rearrangement of the aggregate packing structure, (ii) plastic deformation of small aggregates, and (iii) elastic deformation of constituent particles. Additionally, we developed a semi-analytical model for the compression curve, offering insights into the compressive stages and structural dynamics. Our findings have applications in modeling the internal density profiles of comets and in understanding the early thermal evolution of small icy bodies.

The composite geometry and spectral anisotropy of the solar wind turbulence are very important topics in the investigations of solar wind. In this work, we use the magnetic field and plasma data from Wind spacecraft measured during 1995 January to 2023 December, which covers more than two solar cycles, to systematically investigate these subjects in the context of solar-cycle variability. The so-called spectrum ratio test and spectrum anisotropy test are employed to determine the three-dimensional (3D) geometry of the solar wind turbulence. Both the tests reveal that the solar wind turbulence is dominated by the two-dimensional (2D) component (~80% by turbulence energy). More interestingly, we find that the fraction of slab turbulence increases with the rising sunspot number, and the correlation coefficient between the slab fraction and the sunspot number is 0.61 (ratio test result) or 0.65 (anisotropy test result). This phenomenon suggests that the increasing solar activity (signified by sunspot number) causes increasing slab component in the solar wind turbulence. The relationship between spectral anisotropy and solar activity is discussed and explained. The enhancement of slab fraction is associated with the intensified interplanetary magnetic field magnitude and the increased Alfven speed during the rise phases of the solar cycles. Our findings will be very helpful for achieving a better understanding of the 3D composite geometry and spectral anisotropy of the solar wind turbulence, and especially of their solar-cycle variability.

Gravitational wave (GW) observations offer a promising probe of new physics associated with a first-order electroweak phase transition. Precision studies of the Higgs potential, including Fisher matrix analyses, have been extensively conducted in this context. However, significant theoretical uncertainties in the GW spectrum, particularly those due to renormalization scale dependence in the conventional daisy-resummed approach, have cast doubt on the reliability of such precision measurements. These uncertainties have been highlighted using the Standard Model Effective Field Theory (SMEFT) as a benchmark. To address these issues, we revisit Fisher matrix analyses based on the daisy-resummed approach, explicitly incorporating renormalization scale uncertainties. We then reassess the prospects for precise new physics measurements using GW observations. Adopting the SMEFT as a benchmark, we study the effects of one-loop RGE running of dimension-six operators on the Higgs effective potential via the Higgs self-couplings, top Yukawa coupling, and gauge couplings, in addition to the SMEFT tree-level effects. We find that future GW observations can remain sensitive to various dimension-six SMEFT effects, even in the presence of renormalization scale uncertainties, provided that the SMEFT $(H^{\dagger}H)^3$ operator is precisely measured, e.g., by future collider experiments.

Neutron stars (NSs) are excellent laboratories for testing gravity theories as they are strongly self-gravitating bodies and have rich observational phenomena. However, strong-gravity effects in NS could be degenerate with their equation of state (EOS) which is largely unknown. Fortunately, there exist the so-called universal relations among the NS macroscopic quantities that are found to be insensitive to the underlying EOS. Studying the origin of these relations can lead to a better understanding of NSs and the gravitational interaction. We develop a new perspective of view to analyze the I-C and I-Love universal relations for NSs. At the linear order, we separate the deviation of the universal relations into two parts, where one is the EOS perturbation while the other only depends on the background star structures. The smallness of the second part fully determines the universality of the relation irrespective of the form of the first part. We discuss the validity of our linear approximation when considering the difference among realistic EOSs. Our study can be regarded as a new frame for quantitative representation of the universality and may provide new insights to the universal relations of NSs.

The abundance ratios of radioactive elements U/Th and stable elements Pb/Os from the $r$-process are found to have a strong correlation. This correlation is quite robust with respect to astrophysical conditions. The U/Th-Pb/Os correlation is then applied to provide customized initial abundance ratios U/Th from the observed abundance ratios Pb/Os for six $r$-process enhanced metal-poor stars respectively. Ages of these six metal-poor stars are predicted by the U/Th chronometer, which are approximately between $11$ and $15$ Gyr. Their ages are compatible with the cosmic age of 13.8 billion years predicted from the cosmic microwave background radiation.

We explore the phenomenological implications of the latest Atacama Cosmology Telescope (ACT) DR6 observations, in combination with Planck 2018, BICEP/Keck 2018, and DESI, on the physics of inflation and post-inflationary reheating. We focus on the $\alpha$-attractor class of inflationary models (both E- and T-models) and consider two reheating scenarios: perturbative inflaton ($\phi$) decay ($\phi \rightarrow bb$) and inflaton annihilation ($\phi \phi \rightarrow bb$) into Standard Model (SM) bosonic particles ($b$). By solving the Boltzmann equations, we derive bounds on key reheating parameters, including the reheating temperature, the inflaton equation of state (EoS), and the inflaton-SM coupling, in light of ACT data. To accurately constrain the coupling, we incorporate the Bose enhancement effect in the decay width. To ensure the validity of our perturbative approach, we also identify the regime where nonperturbative effects, such as parametric resonance, become significant. Additionally, we include indirect constraints from primordial gravitational waves (PGWs), which can impact the effective number of relativistic species, $\Delta N_{\rm eff}$. These constraints further bound the reheating temperature, particularly in scenarios with a stiff EoS. Finally, we analyze dark matter (DM) production through purely gravitational interactions during reheating and determine the allowed mass ranges consistent with the constrained reheating parameter space and recent ACT data. Our results highlight the constraining power of next-generation CMB datasets on inflationary dynamics, reheating, and the DM parameter space.