Post-starburst galaxies (PSBs) are a population of rapidly quenched galaxies that have experienced a recent starburst. Several aspects of the origin and nature of PSBs are currently debated. For example, some PSBs unexpectedly host substantial molecular gas despite their low star-formation activity. Furthermore, the relative role of galaxy mergers, interaction, and stellar or black hole feedback in the formation of PSBs remains unclear. We investigate the nature of PSBs at z=0.7 and z=1 selected via rest-frame optical photometric properties, using the FIREbox cosmological simulation. The fraction of PSBs in FIREbox is in reasonable agreement with observational estimates. However, only ~10% of PSBs in FIREbox are `true' PSBs with temporarily quenched star-formation activity. Most galaxies selected as PSBs in FIREbox are impostors, having star-formation rates and gas masses that are comparable to those of non-PSB galaxies with similar stellar masses. Furthermore, PSBs rich in molecular gas are almost exclusively impostors. `True' PSBs and impostors differ on average in their mid-to-near infrared ratio but with significant scatter. By combining the impostor frequency among star forming galaxies in FIREbox with observational estimates, we predict that more than 70% of `true' PSBs in the Universe should originate from massive black hole feedback.
Observational and theoretical studies have long held that rapid gas consumption in starbursts is responsible for the formation of quiescent galaxies. However, studies of recently quenched ``post-starburst" galaxies have discovered that a number of them are surprisingly luminous in CO, challenging this assumption. We present deep ALMA CO(2-1) observations of 50 massive ($\log(M_\star/M_\odot)\sim11.2$) post-starburst galaxies from the SQuIGG$\vec{L}$E sample at $z\sim0.7$. We detect a large fraction (27/50) of the galaxies in CO(2-1). Furthermore, we find that the CO luminosity correlates with the age of the recent starburst, suggesting a gas-removal timescale of $\lesssim140$ Myr, an order of magnitude shorter than is implied by their rest optical star formation rates. We perform new spectral energy distribution fits incorporating mid- and far-IR photometry to test whether dust-obscured star formation can explain this trend. We find that while allowing for buried star formation can raise star formation rates by $\sim0.5$ dex, for almost all galaxies it is neither required to fit the observed IR SED, nor is it sufficient to explain the observed depletion trend. Even the combination of significant buried star formation and ULIRG-like $\alpha_{CO}$ is not enough to explain this decay in CO luminosity. Furthermore, there is no strong evidence to support either of those modifications to the depletion time. Therefore, it remains a distinct possibility that the age-CO luminosity trend should not be interpreted as an evolutionary sequence, and that gas-rich SQuIGG$\vec{L}$E galaxies will soon rejuvenate.
Observations from Juno and Cassini suggest that Jupiter and Saturn may possess fuzzy cores -- central regions where the abundance of heavy elements varies smoothly with depth. Such gradients pose a longstanding puzzle for models of planetary evolution and formation, which predict that vigorous convection would homogenize the interior of gas giants within the first $\sim 10^6$--$10^8~\mathrm{years}$ of cooling. Previous 3D simulations and analytic predictions for the propagation of a convection zone into a stable region have demonstrated that the rapid rotation of gas giants can significantly slow convective mixing, but not enough to stop it. Another piece of the puzzle is luminosity. Gas giants cool as they age, and with that comes a declining heat flux over time. Recent ideas suggest that when this declining luminosity is combined with rotational effects, convection may stall. We explore this possibility using 3D hydrodynamic simulations that include both rotation and a surface cooling flux that decreases as $1/t$. Our results demonstrate that, even without rotation, a declining luminosity can suppress mixing sufficiently to preserve an initial compositional gradient in the deep interior of gas giants. If confirmed by more realistic simulations, this may help to explain the long-term survival of fuzzy cores.
The energy range from a few hundred keV to a few MeV includes important probes such as nuclear gamma-rays and the 511 keV annihilation line. However, compared to X-rays and GeV/TeV gamma-rays, this range suffers from lower sensitivity by orders of magnitude. The upcoming NASA SMEX satellite mission Compton Spectrometer and Imager (COSI), scheduled for launch in 2027, is expected to break through this limitation with its Compton telescope utilizing a germanium semiconductor detector, covering the 0.2-5 MeV energy range. In addition to the main instrument, two Background and Transient Observer (BTO) detectors will be installed on COSI. The detectors are NaI(Tl) scintillators coupled with SiPMs, and they are being developed as a student collaboration project. BTO aims to 1) measure background radiation in orbit to maximize COSI's sensitivity and 2) detect GRBs and other gamma-ray transients. For this purpose, it is required to cover the lower-energy range from 30 keV to 2 MeV with < 20% FWHM energy resolution. In addition, large signals and afterglow generated by heavy ions penetrating the NaI(Tl) crystal should be appropriately handled. To address these requirements, we have developed a compact signal processing and data acquisition system comprised of two main components: an analog board and a digital board. The analog board amplifies signals from the SiPM, generates triggers, and performs AD conversion. The digital board features a Microchip SAMV71 microcontroller, and we developed the software to control the analog board, read ADC data via SPI interface, add timestamps, and buffer event data. Through this development, we achieved the required wide bandpass and an energy resolution of 10% FWHM at 662 keV with a processing time of 20 us per event. We also implemented veto signal generation for large signals using a discriminator and an onboard detection algorithm for transient events.
White dwarfs (WDs) represent the final evolutionary stage of most stars, typically originating from progenitor stars with masses below approximately 8 $M_{\odot}$ to 10 $M_{\odot}$. Formation through single-star evolution generally requires at least 25 Myr, with the youngest WDs often near the Chandrasekhar limit of 1.4 $M_{\odot}$. In contrast, WDs formed via binary channels, such as mergers or mass transfer, can develop smaller masses in a shorter timescale and may exhibit unique characteristics, including strong surface magnetic fields and rapid rotation. Accurately determining the ages of these WDs is essential for understanding their formation. A valuable method involves studying WDs in star clusters, where member stars share the same age and chemical composition, allowing for precise constraints on the formation times and metallicities of the WDs' progenitors. Here we report a WD found in the open cluster RSG 5, which is only 35 Myr old. The WD's mass is lower than 1.05 $M_{\odot}$, indicating it may not have formed through single-star evolution. The WD possesses an exceptionally strong surface magnetic field ($\ge 200$ MG), a short rotational period ($\sim 6.5$ min), and, most notably, a co-rotating half-ring of ionized circumstellar debris. This distinctive feature provides evidence for a binary merger origin, a scenario further substantiated by our stellar evolution models.
The measurements of pressure-mode and mixed-mode oscillation frequencies in subgiant stars offer a unique opportunity to probe their internal structures -- from the surface to the deep interior -- and to precisely determine their global properties. We have conducted a detailed asteroseismic analysis of the benchmark subgiant $\mu$ Herculis using eight seasons of radial velocity observations from the SONG telescope, and have determined its mass, radius, age, and surface helium abundance to be $1.105_{-0.024}^{+0.058}$ M$_\odot$, $1.709_{-0.015}^{+0.030}$ R$_\odot$, $8.4_{-0.1}^{+0.4}$ Gyr, and $0.242^{+0.006}_{-0.021}$, respectively. We have demonstrated that simultaneously fitting the helium glitch properties, oscillation frequencies, and spectroscopic observables yields a more accurate inference of the surface helium abundance and hence stellar age. A significant discrepancy between the observed extent of the helium ionization zone and that predicted by stellar models is identified and examined, underscoring potential limitations in the current modelling of stellar interiors. Our analysis confirms that the helium glitch originates from the region between the two stages of helium ionisation, i.e. from the $\Gamma_1$ peak, rather than from the second helium ionisation zone itself. Within the conventional formalism, this implies that the glitch analysis characterises the region located between the two helium ionisation zones.
Star-forming galaxies are in a state of turbulence, with one of the principle components of the turbulence sourced by the constant injection of momentum from supernovae (SNe) explosions. Utilizing high-resolution stratified, gravito-hydrodynamical models of SNe-driven turbulence with interstellar medium (ISM) cooling and heating, we explore how SNe-driven turbulence changes across different galactic conditions, parameterized by the galactic mass and potential, SNe-driving rate, and seeding functions. We show that even though the underlying ISM changes between starburst and Milky Way analogue models, the velocity fluctuations in the turbulence of both models, but not the kinetic energy fluctuations, can be normalized into a universal, single cascade, $u^2(k) \propto k^{-3/2}$, where $u$ is the velocity and $k$ is the wavemode, indicating that the structure of the turbulence is robust to significant changes in the ISM and SNe seeding. Moreover, the cascades connect smoothly from the winds into the galactic disk, pushing the outer-scale of the turbulence, $\ell_{\rm cor}$, to over $\ell_{\rm cor} \approx 6 \ell_0$, where $\ell_0$ is the gaseous scale-height. By providing an analytical model for the sound speed spectrum, $c_s^2(k)$, in the weak-cooling, adiabatic limit, we show that it is the compressible turbulent modes, $u_c$, that control the volume-filling phase structure of the galactic disks in our models, with $c_s^2(k) \propto k^{-2} \propto u_c^2(k)$. This may indicate that galactic turbulence does not only have highly-universal features across different galaxies, but also directly sets the volume-filling hot and warm phase structure of the underlying galactic ISM through turbulent compressible modes.
Spiral arms play a central role in disc galaxies, but their dynamical nature remains a long-standing open question. Azimuthal offsets between molecular gas and star formation are expected if gas crosses spiral arms, as predicted by quasi-stationary density wave theory. In this work, we measure offsets between CO and Halpha peaks in radial bins for 24 galaxies from the PHANGS survey that display a well-delineated spiral structure. The offsets exhibit substantial scatter, implying that star formation is not exclusively initiated at a coherent spiral shock. We define offsets such that positive values mean Halpha peaks lie ahead of CO peaks in the direction of galactic rotation. With this convention, 14 galaxies show mean positive CO-Halpha offsets, typically of a few hundred parsecs. In four of these 14 galaxies (17% of the total), offsets become smaller with increasing radius, as expected for a single quasi-stationary spiral density wave. Ten galaxies (42%) show positive mean offsets but no clear correlation with radius, which is compatible with multiple overlapping modes. In the remaining ten galaxies (42%), we find no significantly positive offsets, which could point to transient dynamical spirals or material arms, where gas and stars co-rotate with the spiral perturbation. Across the full sample, we find mostly positive offsets between CO peaks and the gravitational potential minimum, confirming that gas often crosses the spiral perturbation. For the four galaxies with clear positive offsets and a radial trend, we derived pattern speeds in good agreement with the literature. Overall, our results suggest that even well-delineated spirals in the local Universe can arise from a variety of underlying dynamical mechanisms.
Accurate and parsimonious quantification of baryonic feedback on matter distribution is of crucial importance for understanding both cosmology and galaxy formation from observational data. This is, however, challenging given the large discrepancy among different models of galaxy formation simulations, and their distinct subgrid physics parameterizations. Using 5,072 simulations from 4 different models covering broad ranges in their parameter spaces, we find a unified 2D latent representation. Compared to the simulations and other phenomenological models, our representation is independent of both time and cosmology, much lower-dimensional, and disentangled in its impacts on the matter power spectra. The common latent space facilitates the comparison of parameter spaces of different models and is readily interpretable by correlation with each. The two latent dimensions provide a complementary representation of baryonic effects, linking black hole and supernova feedback to distinct and interpretable impacts on the matter power spectrum. Our approach enables developing robust and economical analytic models for optimal gain of physical information from data, and is generalizable to other fields with significant modeling uncertainty.
We present the first comprehensive catalogs of $^{13}$CO clumps from the Milky Way Imaging Scroll Painting (MWISP) project. By developing an equivalent global detection scheme integrated with the FacetClumps algorithm, we successfully extract 71,661 molecular clumps across a high-resolution $^{13}$CO data cube spanning 2310 deg$^2$ from the MWISP Phase I survey. To determine accurate distances, we design an automatic hierarchical distance decision method using signal regions as fundamental objects, effectively resolving the kinematic distance ambiguity problem and obtaining reliable measurements for 97.94% of the sample. Statistical analysis reveals that 65.3% of clumps are gravitationally bound, accounting for approximately 96.3% of the statistical total mass. Scaling relation analysis across multiple surveys reveals universal power-law behaviors in clump populations. Maser-associated clumps exhibit modified parameter distributions and scaling relations, revealing how active star formation alters clump dynamics and structure. These extensive catalogs establish a foundation for investigating molecular clump properties, star formation processes, and Galactic evolution.
We present cosmological constraints from the Planck PSZ2 cosmological cluster sample, using weak-lensing shear profiles from Dark Energy Survey (DES) data and X-ray observations from the Chandra telescope for the mass calibration. We compute hydrostatic mass estimates for all clusters in the PSZ2 sample with a scaling relation between their Sunyaev-Zeldovich signal and X-ray derived hydrostatic mass, calibrated with the Chandra data. We introduce a method to correct these masses with a hydrostatic mass bias using shear profiles from wide-field galaxy surveys. We simultaneously fit the number counts of the PSZ2 sample and the mass calibration with the DES data, finding $\Omega_\text{m}=0.312^{+0.018}_{-0.024}$, $\sigma_8=0.777\pm 0.024$, $S_8\equiv \sigma_8 \sqrt{\Omega_\text{m} / 0.3}=0.791^{+0.023}_{-0.021}$, and $(1-b)=0.844^{+0.055}_{-0.062}$ for our baseline analysis when combined with BAO data. When considering a hydrostatic mass bias evolving with mass, we find $\Omega_\text{m}=0.353^{+0.025}_{-0.031}$, $\sigma_8=0.751\pm 0.023$, and $S_8=0.814^{+0.019}_{-0.020}$. We verify the robustness of our results by exploring a variety of analysis settings, with a particular focus on the definition of the halo centre used for the extraction of shear profiles. We compare our results with a number of other analyses, in particular two recent analyses of cluster samples obtained from SPT and eROSITA data that share the same mass calibration data set. We find that our results are in overall agreement with most late-time probes, in very mild tension with CMB results (1.6$\sigma$), and in significant tension with results from eROSITA clusters (2.9$\sigma$). We confirm that our mass calibration is consistent with the eROSITA analysis by comparing masses for clusters present in both Planck and eROSITA samples, eliminating it as a potential cause of tension.
this https URL . Published version PDF is attached
this https URL ). Accepted for publication in A&A