Papers for Wednesday, Oct 23 2024

Modeling strong gravitational lenses is prohibitively expensive for modern and next-generation cosmic survey data. Neural posterior estimation (NPE), a simulation-based inference (SBI) approach, has been studied as an avenue for efficient analysis of strong lensing data. However, NPE has not been demonstrated to perform well on out-of-domain target data -- e.g., when trained on simulated data and then applied to real, observational data. In this work, we perform the first study of the efficacy of NPE in combination with unsupervised domain adaptation (UDA). The source domain is noiseless, and the target domain has noise mimicking modern cosmology surveys. We find that combining UDA and NPE improves the accuracy of the inference by 1-2 orders of magnitude and significantly improves the posterior coverage over an NPE model without UDA. We anticipate that this combination of approaches will help enable future applications of NPE models to real observational data.

We explore the potential of an array of O(100) small fixed telescopes, aligned along a meridian and automated to measure millions of occultations of Gaia stars by minor planets, to constrain gravitational signatures from a "Planet X" mass in the outer solar system. The accuracy of center-of-mass tracking for the occulters is limited by photon noise, uncertainties in asteroid shapes, and Gaia's astrometry of the occulted stars. Using both parametric calculations and survey simulations, we assess the total information obtainable from occultation measurements of main-belt asteroids (MBAs), Jovian Trojans and trans-Neptunian objects (TNOs). We find that MBAs are the optimal target population due to their higher occultation rates and abundance of objects above LSST detection thresholds. A 10-year survey of occultations by MBAs and Trojans using an array of 200 40 cm telescopes at 5 km separation would achieve $5\sigma$ sensitivity to the gravitational tidal field of a $5M_\oplus$ Planet X at 800 AU for $>90\%$ of potential sky locations. This configuration corresponds to an initial cost of $\approx $15 million. While the occultation survey's sensitivity to tidal forces improves rapidly with increasing number of telescopes, sensitivity to a Planet X becomes limited by degeneracy with the uncertain masses of large moonless TNOs. The 200-telescope survey would additionally detect $\approx $1800 TNO occultations, providing detailed shape, size, and albedo information. The survey would also measure the Yarkovski effect on many individual MBAs, measure the masses of many asteroids involved in mutual gravitational deflections, and enable better searches for primordial black holes and departures from General Relativity.

We use the Atacama Large sub/Millimetre Array (ALMA) to efficiently observe spectral lines across Bands 3, 4, 5, 6, 7, and 8 at high-resolution (0.5" - 0.1") for 16 bright southern Herschel sources at $1.5 < z < 4.2$. With only six and a half hours of observations, we reveal 66 spectral lines in 17 galaxies. These observations detect emission from CO (3-2) to CO(18-17), as well as atomic ([CI](1-0), (2-1), [OI] 145 $\mu$m and [NII] 205 $\mu$m) lines. Additional molecular lines are seen in emission (${\rm H_2O}$ and ${\rm H_2O^+}$) and absorption (OH$^+$ and CH$^+$). The morphologies based on dust continuum ranges from extended sources to strong lensed galaxies with magnifications between 2 and 30. CO line transitions indicate a diverse set of excitation conditions with a fraction of the sources ($\sim 35$%) showcasing dense, warm gas. The resolved gas to star-formation surface densities vary strongly per source, and suggest that the observed diversity of dusty star-forming galaxies could be a combination of lensed, compact dusty starbursts and extended, potentially-merging galaxies. The predicted gas depletion timescales are consistent with 100 Myr to 1 Gyr, but require efficient fueling from the extended gas reservoirs onto the more central starbursts, in line with the Doppler-shifted absorption lines that indicate inflowing gas for two out of six sources. This pilot paper explores a successful new method of observing spectral lines in large samples of galaxies, supports future studies of larger samples, and finds that the efficiency of this new observational method will be further improved with the planned ALMA Wideband Sensitivity Upgrade.

The evolution of galaxy sizes in different wavelengths provides unique insights on galaxy build-up across cosmic epochs. Such measurements can now finally be done at $z>3$ thanks to the exquisite spatial resolution and multi-wavelength capability of the JWST. With the public data from the CEERS, PRIMER-UDS, and PRIMER-COSMOS surveys, we measure the sizes of $\sim 3500$ star-forming galaxies at $3 \leqslant z<9$, in 7 NIRCam bands using the multi-wavelength model fitting code GalfitM. The size-mass relation is measured in four redshift bins, across all NIRCam bands. We find that, the slope and intrinsic scatter of the rest-optical size-mass relation are constant across this redshift range and consistent with previous HST-based studies at low-z. When comparing the relations across different wavelengths, the average rest-optical and rest-UV relations are consistent with each other up to $z=6$, but the intrinsic scatter is largest in rest-UV wavelengths compared to rest-optical and redder bands. This behaviour is independent of redshift and we speculate that it is driven by bursty star-formation in $z>4$ galaxies. Additionally, for $3\leqslant z<4$ star-forming galaxies at $\rm M_* > 10^{10} M_{\odot}$, we find smaller rest-$\rm 1\rm\,\mu m$ sizes in comparison to rest-optical (and rest-UV) sizes, suggestive of colour gradients. When comparing to simulations, we find agreement over $\rm M_* \approx 10^{9} - 10^{10} M_{\odot}$ but beyond this mass, the observed size-mass relation is significantly steeper. Our results show the power of JWST/NIRCam to provide new constraints on galaxy formation models.

Parameter estimation from galaxy survey data from the full-shape method depends on scale cuts and priors on EFT parameters. The effects of priors, including the so-called "prior volume" phenomenon have been originally studied in Ivanov et al. (2019) and subsequent works. In this note, we repeat and extend these tests and also apply them to the priors used by D'Amico et al. (2019). We point out that in addition to the "prior volume" effect there is a much more dangerous effect that is largely overlooked: a systematic bias on cosmological parameters due to overoptimistic scale cuts. Unlike the "prior volume" effect, this is a genuine systematic bias due to two-loop corrections that does not vanish with better priors or with larger data volumes. Our study is based on the high fidelity BOSS-like PT Challenge simulation data which offer many advantages over analyses based on synthetic data generated with fitting pipelines. We show that the analysis choices of D'Amico et al.~(2019), especially their scale cuts, significantly bias parameter recovery, overestimating $\sigma_8$ by over $5\%$ (equivalent to $1\sigma$). The bias on measured EFT parameters is even more significant. In contrast, the analysis choices associated with the CLASS-PT code lead to a negligible ($\lesssim 1\%$) bias in the recovery of cosmological parameters based on their best-fit values.

The circumgalactic medium (CGM) in $\gtrsim 10^{12}$ $\mathrm{M}_{\odot}$ halos is dominated by a hot phase ($T \gtrsim 10^{6}$ K). While many models exist for the hot gas structure, there is as yet no consensus. We compare cooling flow models, in which the hot CGM flows inward due to radiative cooling, to the CGM of $\sim 10^{12}-10^{13}$ $\mathrm{M}_{\odot}$ halos in galaxy formation simulations from the FIRE project at $z\sim0$. The simulations include realistic cosmological evolution and feedback from stars but neglect AGN feedback. At both mass scales, CGM inflows are typically dominated by the hot phase rather than by the `precipitation' of cold gas. Despite being highly idealized, we find that cooling flows describe $\sim 10^{13}$ $\mathrm{M}_{\odot}$ halos very well, with median agreement in the density and temperature profiles of $\sim 20\%$ and $\sim 10\%$, respectively. This indicates that stellar feedback has little impact on CGM scales in those halos. For $\sim 10^{12}$ $\mathrm{M}_{\odot}$ halos, the thermodynamic profiles are also accurately reproduced in the outer CGM. For some of these lower-mass halos, cooling flows significantly overpredict the hot gas density in the inner CGM. This could be due to multidimensional angular momentum effects not well captured by our 1D cooling flow models and/or to the larger cold gas fractions in these regions. Turbulence, which contributes $\sim 10-40\%$ of the total pressure, must be included to accurately reproduce the temperature profiles. Overall, cooling flows predict entropy profiles in better agreement with the FIRE simulations than other idealized models in the literature.

The AEOS project introduces a series of high-resolution cosmological simulations that model star-by-star chemical enrichment and galaxy formation in the early Universe, achieving 1 pc resolution. These simulations capture the complexities of galaxy evolution within the first ~300 Myr by modeling individual stars and their feedback processes. By incorporating chemical yields from individual stars, AEOS generates galaxies with diverse stellar chemical abundances, linking them to hierarchical galaxy formation and early nucleosynthetic events. These simulations underscore the importance of chemical abundance patterns in ancient stars as vital probes of early nucleosynthesis, star formation histories, and galaxy formation. We examine the metallicity floors of various elements resulting from Pop III enrichment, providing best-fit values for eight different metals (e.g., [O/H] = -4.0) to guide simulations without Pop III models. Additionally, we identify galaxies that begin star formation with Pop II after external enrichment and investigate the frequency of CEMP stars at varying metallicities. The AEOS simulations offer detailed insights into the relationship between star formation, feedback, and chemical enrichment. Future work will extend these simulations to later epochs to interpret the diverse stellar populations of the Milky Way and its satellites.

Gaia characterizes the stellar populations of nearby open clusters with unprecedented precision. We investigate the Böhm-Vitense gap, which has been found as a prominent feature in the stellar sequence of open clusters. Using PARSEC isochrone fitting, we derive astrophysical parameters for more than 1100 stars in Praesepe, identify more than 1100 bona fide single stars in the alpha Persei (Melotte 20) open cluster, and confirm their approximate match in terms of age (approx. 710 Myr and 45 Myr) and metallicity ([M/H] approx. +0.15 dex and +0.13 dex) to the Hyades and Pleiades, respectively. By merging data of the cluster pairs, we improve number statistics. We do not find a clear gap in the combined observational G_abs vs. BP-RP color-magnitude diagram (CMD) in the stellar mass range corresponding to the location of the Böhm-Vitense gap. We reproduce gaps in simulated Hyades-type CMDs randomly drawn from an initial mass function. There is no strong evidence for a discontinuity originating in the transition from radiative to convective energy transport in the stellar photosphere. We conclude that the observed gaps in the stellar sequences of open clusters could be explained by small number statistics and the uneven mass-color relation at the transition from spectral type A to F.

Classifying ionised nebulae in nearby galaxies is crucial to studying stellar feedback mechanisms and understanding the physical conditions of the interstellar medium. This classification task is generally performed by comparing observed line ratios with photoionisation simulations of different types of nebulae (HII regions, planetary nebulae, and supernova remnants). However, due to simplifying assumptions, such simulations are generally unable to fully reproduce the line ratios in observed nebulae. This discrepancy limits the performance of the classical machine-learning approach, where a model is trained on the simulated data and then used to classify real nebulae. In this study, we use a domain-adversarial neural network (DANN) to bridge the gap between photoionisation models (source domain) and observed ionised nebulae from the PHANGS-MUSE survey (target domain). The DANN is an example of a domain-adaptation algorithm, whose goal is to maximise the performance of a model trained on labelled data in the source domain on an unlabelled target domain by extracting domain-invariant features. Our results indicate a significant improvement in classification performance in the target domain when employing the DANN framework compared to a classical neural network (NN) classifier. Additionally, we investigate the impact of adding noise to the source dataset, finding that noise injection acts as a form of regularisation, further enhancing the performances of both the NN and DANN models on the observational data. The combined use of domain adaptation and noise injection improves the classification accuracy in the target domain by 24%. This study highlights the potential of domain adaptation methods in tackling the domain-shift challenge when using theoretical models to train machine-learning pipelines in astronomy.

Observations suggest that filaments in molecular clouds can grow by mass accretion while forming cores via fragmentation. Here we present one of the first large sample studies of filament accretion using velocity gradient measurements of star-forming filaments on the $\sim 0.05$ pc scale with NH$_3$ observations of the Perseus Molecular Cloud, primarily obtained as a part of the GBT Ammonia Survey (GAS). In this study, we find significant correlations between velocity gradient, velocity dispersion, mass per unit length, and the number of cores per unit length of the Perseus filaments. Our results suggest a scenario in which filaments not only grow through mass accretion but also form new cores continuously in the process well into the thermally supercritical regime. Such behavior is contrary to that expected from isolated filament models but consistent with how filaments form within a more realistic cloud environment, suggesting that the cloud environment plays a crucial role in shaping core formation and evolution in filaments. Furthermore, even though velocity gradients within filaments are not oriented randomly, we find no correlation between velocity gradient orientation and the filament properties we analyzed. This result suggests that gravity is unlikely the dominant mechanism imposing order on the $\sim 0.05$ pc scale for dense star-forming gas.

Unravelling galaxy formation theory requires understanding galaxies both at high and low redshifts. A possible way to connect both realms is by studying the oldest stars in the Milky Way (i.e., the proto-Galaxy). We use the $APOGEE$-$Gaia$ surveys to perform a purely chemical dissection of Milky Way (MW) stellar populations, and identify samples of stars likely belonging to proto-Galactic fragments. The metallicity dependence of the distribution of old MW stars in the [Mg/Mn]-[Al/Fe] enables the distinction of at least two populations in terms of their star formation histories: a rapidly evolved population likely associated with the main progenitor system of the proto-MW; and populations characterised by less efficient, slower, star formation. In the Solar neighbourhood these populations are dominated by the $Gaia-Enceladus/Sausage$ accretion debris. In the inner Galaxy, they are largely associated with the $Heracles$ structure. We model the density of chemically defined proto-Galaxy populations, finding that they are well represented by a Plummer model with a scale radius of $a$ ~3.5 kpc, and an oblate ellipsoid with flattening parameters [$p$ ~0.8; $q$ ~0.6]; this finding indicates that the MW plausibly hosts a low-mass, metal-poor, bulge component. We integrate this density for $chemically$ $unevolved$ stars between $-2$ < [Fe/H] < $-0.5$ to obtain a minimum stellar mass for the proto-Galaxy of $M_{*}$ ($r$ < 10 kpc) = $9.1\pm0.2\times10^{8}$ M$_{\odot}$. Our results suggest the proto-Milky Way is at least comprised of two significant fragments: the main $in$ $situ$ progenitor and the $Heracles$ structure.

Gravitational waves (GWs) offer a new observational window into the universe, providing insights into compact objects and cosmic structures. Gravitational lensing, commonly studied in electromagnetic waves, also affects GWs, introducing magnification, time delays, and multiple images. While existing studies focus on static lenses, many astrophysical lenses are dynamic, with time-varying mass distributions such as moving stars or orbiting binaries. We develop a general theoretical framework to describe non-static lenses and demonstrate how they modulate GW signals, inducing unique time-varying amplitude modulations and spectral broadening. By examining uniformly moving and orbiting binary lenses, we show that these modulations provide new observational signatures, enhancing our understanding of lensing objects and the dynamics of the universe. Our findings have important implications for GW astronomy, offering novel ways to probe lens dynamics and improve the interpretations of GW signals.

We propose that a dark matter (DM) spike around the Galactic Center's (GC) supermassive black hole, Sgr A*, could account for most of the bulge's measured 511 keV line intensity while remaining cosmologically compatible. DM annihilation can be the primary source of the 511 keV line emission without violating constraints from disk emission observations and in-flight positron annihilation with the interstellar medium, provided the disk emission is dominated by an astrophysical source of low-energy positrons. We find that a DM mass up to approximately 20 MeV, either with a Gondolo-Silk spike or one softened by stellar heating, could explain the observed 511 keV bulge emission profile. Our proposal can be tested by future observations of the continuum diffuse emission close to the GC.

Most of the dynamical mass loss from star clusters is thought to be caused by the time-variability of the tidal field (``tidal shocks''). Systematic studies of tidal shocks have been hampered by the fact that each tidal history is unique, implying both a reproducibility and a generalisation problem. Here we address these issues by investigating how star cluster evolution depends on the statistical properties of its tidal history. We run a large suite of direct N-body simulations of clusters with tidal histories generated from power spectra of a given slope and with different normalisations, which determine the time-scales and amplitudes of the shocks, respectively. At fixed normalisation (i.e. the same median tidal field strength), the dissolution time-scale is nearly independent of the power spectrum slope. However, the dispersion in dissolution time-scales, obtained by repeating simulations for different realisations of statistically identical tidal histories, increases with the power spectrum slope. This result means that clusters experiencing high-frequency shocks have more similar mass loss histories than clusters experiencing low-frequency shocks. The density-mass relationship of the simulated clusters follows a power-law with slope between 1.08 and 1.45, except for the lowest normalisations (for which clusters effectively evolve in a static tidal field). Our findings suggest that star cluster evolution can be described statistically from a time-series analysis of its tidal history, which is an important simplification for describing the evolution of the star cluster population during galaxy formation and evolution.

We present archival HST and JWST ultraviolet through near infrared time series photometric observations of a massive minimal-contact binary candidate in the metal-poor galaxy WLM ($Z = 0.14 Z_{\odot}$). This discovery marks the lowest metallicity contact binary candidate observed to date. We determine the nature of the two stars in the binary by using the eclipsing binary modeling software (PHysics Of Eclipsing BinariEs; PHOEBE) to train a neural network to fit our observed panchromatic multi-epoch photometry. The best fit model consists of two hot MS stars ($T_1=29800^{+2300}_{-1700}$ K, $M_1=16^{+2}_{-3}~M_{\odot}$, and $T_2=18000^{+5000}_{-5000}$ K, $M_2=7^{+5}_{-3}~M_{\odot}$). We discuss plausible evolutionary paths for the system, and suggest the system is likely to be currently in a contact phase before ultimately ending in a merger. Future spectroscopy will help to further narrow down evolutionary pathways. This work showcases a novel use of data of JWST and HST imaging originally taken to characterize RR Lyrae. We expect time series imaging from LSST, BlackGEM, etc. to uncover similar types of objects in nearby galaxies.

We search for a spatial coincidence between high energy neutrinos detected by the IceCube neutrino detector and ten red dwarfs which have been observed in gamma-rays. For our analysis, we use the unbinned maximum likelihood method to look for a statistically significant excess. We do not find any such spatial association between any of the red dwarfs and Icecube-detected neutrinos. Therefore, we conclude that none of the gamma-ray bright red dwarfs contribute to the diffuse neutrino flux measured by IceCube.

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides deep HST F275W and F435W imaging over four CANDELS fields (GOODS-N, GOODS-S, COSMOS, and EGS). We combine this newly acquired UV imaging with existing HST imaging from CANDELS as well as existing ancillary data to obtain robust photometric redshifts and reliable estimates for galaxy physical properties for over 150,000 galaxies in the $\sim$430 arcmin$^2$ UVCANDELS area. Here, we leverage the power of the new UV photometry to not only improve the photometric redshift measurements in these fields, but also constrain the full redshift probability distribution combining multiple redshift fitting tools. Furthermore, using the full UV-to-IR photometric dataset, we measure the galaxy physical properties by fitting templates from population synthesis models with two different parameterizations (flexible and fixed-form) of the star-formation histories (SFHs). Compared to the flexible SFH parametrization, we find that the fixed-form SFHs systematically underestimate the galaxy stellar masses, both at the low- ($\lesssim10^9 M_\odot$) and high- ($\gtrsim10^{10} M_\odot$) mass end, by as much as $\sim0.5$ dex. This underestimation is primarily due the limited ability of fixed-form SFH parameterization to simultaneously capture the chaotic nature of star-formation in these galaxies.

Massive contact binaries are fundamental test-beds to study properties of close binaries en route to stellar mergers. Progress in their understanding has been largely impeded due to a dearth in their observed counterparts. We accumulate 61 over-contact binary candidates from the OGLE survey of the Magellanic Clouds and form the largest known sample of such systems. We mine the OGLE-IV catalogue for B-type and earlier systems which match our synthetic light curves, and thereafter examine their properties in comparison with our binary populations. Simultaneously, we also investigate the degeneracy in the light curves between over-contact and near-contact binaries and attempt to distinguish these populations. We construct new synthetic populations of such systems from detailed $\texttt{MESA}$ models with initial parameters of M$_\textrm{total,i}=14-80\,\textrm{M}_{\odot}$ and P$_\textrm{i}=0.6-23\,$days, and compute analogous $\texttt{PHOEBE}$ light curves. We then extract OGLE systems that match our synthetic light curves and compare their numbers, and distributions in periods, absolute V-band magnitudes and colours (V-I) with our population predictions. We expect massive over-contact binaries to have smooth equal-depth light curves due to their equal-mass ($q\approx1$) and consequently, equal-temperature components. We mine 47 and 14 predominantly B-type over-contact OGLE binaries from the Large and Small Magellanic Clouds respectively, with periods of $\approx0.6-1\,$d, along with a few wider O-type systems. This sample is consistent with our theoretically expected numbers and colour-magnitude and period distributions. Additionally, we find that binaries nearing contact have periods as large as $4.5\,$d and mass ratios of $q\approx0.4-0.95$, some of which may masquerade as over-contact binaries and thereby lead to potential observational misclassifications.

Observational constraints on classical novae are heavily biased to phases near optical peak and later because of the simple fact that novae are not typically discovered until they become bright. The earliest phases of brightening, coming before discovery, are typically missed, but this is changing with the proliferation of wide-field optical monitoring systems including ZTF, ASAS-SN, and Evryscope. Here, we report on unprecedented observations of the fast nova V1674 Her beginning >10 mag below its optical peak and including high-cadence (2 min.) observations that chart a rise of ~8 mag in just 5 hours. Two clear breaks are identified as the light curve transitions first from rising slowly to rising rapidly, followed by a transition to an even faster, nearly linear rate of increasing flux with time. The depths of the observations allow us to place tight constraints on the size of the photosphere under the assumption of blackbody emission from a white dwarf emitting at its Eddington luminosity. We find that the white dwarf was unlikely to have overflowed its Roche lobe prior to the launch of a fast wind, which poses a challenge for explaining the Fermi $\gamma$-ray detections as the interaction of a fast wind with a slow-torus of gas stripped from the inflated white dwarf envelope by the companion. High-cadence observations of novae from Evryscope and the planned Argus Array can record the diversity of rising nova light curves and help resolve how the interplay between thermonuclear fusion, binary interaction, and shocks power their earliest light.

UGC 2885 (z = 0.01935) is one of the largest and most massive galaxies in the local Universe, yet its undisturbed spiral structure is unexpected for such an object and unpredicted in cosmological simulations. Understanding the detailed properties of extreme systems such as UGC 2885 can provide insight on the limits of scaling relations and physical processes driving galaxy evolution. Our goal is to understand whether UGC 2885 has followed a similar evolutionary path to other high-mass galaxies by examining its place on the fundamental metallicity relation and the star-forming main sequence. We present new observations of UGC 2885 with the CFHT and IRAM 30-m telescopes. These novel data are used to respectively calculate metallicity and molecular hydrogen mass values. We estimate stellar mass (M*) and star formation rate (SFR) based on mid-infrared observations with the Wide-field Infrared Survey Explorer. We find global metallicities Z = 9.28, 9.08 and 8.74 at the 25 kpc ellipsoid from N2O2, R23 and O3N2 indices, respectively. This puts UGC 2885 at the high end of the galaxy metallicity distribution. The molecular hydrogen mass is calculated as M(H2)=(1.89+/-0.24)e11 Msun, the SFR as 1.63+/-0.72 Msun/yr and the stellar mass as (4.83 +/- 1.52)e11 Msun, which gives a star formation efficiency (SFE = SFR/M(H2)) of (8.67+/-4.20)e12/yr. This indicates that UGC 2885 has an extremely high molecular gas content when compared to known samples of star forming galaxies (~100 times more) and a relatively low SFR for its current gas content. We conclude that UGC 2885 has gone through cycles of star formation periods, which increased its stellar mass and metallicity to its current state. The mechanisms that are fueling the current molecular gas reservoir and keeping the galaxy from producing stars remain uncertain. We discuss the possibility that a molecular bar is quenching star forming activity.

The coupling of the angular jitter of the spacecraft and their sub-assemblies with the optical bench and the telescope into the interferometric length readout will be a major noise source in the LISA mission. We refer to this noise as tilt-to-length (TTL) coupling. It will be reduced directly by realignments, and the residual noise will then be subtracted in post-processing. The success of these mitigation strategies depends on an accurate computation of the TTL coupling coefficients. We present here a thorough analysis of the accuracy of the coefficient estimation under different jitter characteristics, angular readout noise levels, and gravitational wave sources. We analyze in which cases the estimates degrade using two estimators, the common least squares estimator and the instrumental variables estimator. Our investigations show that angular readout noise leads to a bias of the least squares estimator, depending on the TTL coupling coefficients, jitter and readout noise level, while the instrumental variable estimator is not biased. We present an equation that predicts the estimation bias of the least squares method due to angular readout noise.

Despite the recent discoveries of planets orbiting stars at all evolutionary stages, the evolution of planetary systems remains poorly understood. Studying planetary systems around red giant branch stars can reveal how main sequence planetary systems can change and evolve into white dwarf systems over time. Decades of radial velocity and transit surveys have yielded the detection of hundreds of planets and planet candidates orbiting evolved stars. These planetary systems have provided important insights into understanding how planetary atmospheres and orbits can be disrupted by stellar evolution, potentially being restructured at late stages, and how planets can be eventually engulfed by their stars, possibly reborn as white dwarf planetary systems. Evolved star targets will reveal planet occurrence at the largest distances and most varied environments across the Galaxy.

Habitable Worlds Observatory (HWO) will search for biosignatures from Earth-size exoplanets in the habitable zones of nearby stars. The wavelength range for biosignatures used by the HabEx and LUVOIR mission concept studies was 200 nm to 2 microns and, as such, this is a candidate wavelength range for HWO. The visible wavelength range (500-1000 nm) provides for detection of water, oxygen, and Raleigh scattering; the near-ultraviolet is valuable for detection of ozone; and the near-infrared enables detection of carbon dioxide and methane for Earth-like atmospheres. Damiano et al. 2023 showed the significant improvement in spectral retrieval reliability when the NUV and NIR are both used with the visible. However, the challenge of the NUV, in addition to the technological and engineering challenges of starlight suppression in the NUV, is the drop in flux of host stars. In the NIR, the challenge is the geometric access to the habitable zone due to the wavelength dependency of the inner working angle limit of coronagraphs. For these reasons, exoplanet yields are lower in the NUV and NIR than in the visible (Morgan et al. 2023, Morgan et al. 2024) and some instrument parameters are more critical for improving NUV and NIR yields than others. In this paper we present a new capability for performing a large number of end-to-end yield modeling simulations to enable large, multivariate parameter sweeps. We utilize this capability to calculate the Visible, NIR, and NUV yield sensitivities to the instrument parameters: aperture diameter, coronagraph core throughput, contrast, and inner working angle (IWA). We find that parameter interactions are important in determining yield, the most important of which is the interaction between contrast and IWA, but that the strength of that interaction is different in each of the three wavebands.

Active galactic nuclei (AGN) are a promising source of the binary black hole (BBH) mergers observed in gravitational waves with LIGO-Virgo-Kagra (LVK). Constraining the AGN channel allows us to limit AGN parameter space (disk density, size, average lifetime) and nuclear star cluster (NSC) parameter space. Constraints on AGN and NSCs have implications for $\Lambda$CDM models of AGN feedback and models of AGN-driven SMBH merger and growth. Here we present several qualitative studies of the AGN channel using new public, open-source, fast, reproducible code \texttt{McFACTS}\footnote{this https URL}:Monte Carlo for AGN channel Testing \& Simulation. We demonstrate several important features for testing the AGN channel, including: i) growth to large mass IMBH is helped by the presence of migration traps or disk boundaries, ii) flat BH initial mass functions highlight hierarchical merger features in the mass spectrum, iii) high rates of dynamical encounters strongly inhibit BBH formation and merger at migration traps, iv) the ($q,\chi_{\rm eff}$) anti-correlation is a strong test of the bias to prograde mergers in the AGN channel, v) spheroid encounters can drive a fraction of mergers with high in-plane spin components ($\chi_{\rm p}$), vi) a high rate of extreme mass ratio inspirals (EMRIs) are driven by an initial population of embedded retrograde BH, vii) Both LVK and LISA are powerful probes of models of AGN disks and their embedded populations.

The most prominent of the IRAS debris disk systems, $\alpha$ Lyrae (Vega), at a distance of 7.7 pc, has been observed by both the NIRCam and MIRI instruments on the James Webb Space Telescope (JWST). This paper describes NIRCam coronagraphic observations which have achieved F444W contrast levels of 3$\times10^{-7}$ at 1\arcsec\ (7.7 au), 1$\times10^{-7}$ at 2\arcsec\ (15 au) and few $\times 10^{-8}$ beyond 5\arcsec\ (38 au), corresponding to masses of $<$ 3, 2 and 0.5 MJup for a system age of 700 Myr. Two F444W objects are identified in the outer MIRI debris disk, around 48 au. One of these is detected by MIRI, appears to be extended and has a spectral energy distribution similar to those of distant extragalactic sources. The second one also appears extended in the NIRCam data suggestive of an extragalactic this http URL NIRCam limits within the inner disk (1\arcsec\ --10\arcsec) correspond to a model-dependent masses of 2$\sim$3 \mj. \citet{Su2024} argue that planets larger even 0.3 MJup would disrupt the smooth disk structure seen at MIRI wavelengths. Eight additional objects are found within 60\arcsec\ of Vega, but none has astrometric properties or colors consistent with planet candidates. These observations reach a level consistent with expected Jeans Mass limits. Deeper observations achieving contrast levels $<10^{-8}$ outside of $\sim$4\arcsec\ and reaching masses below that of Saturn are possible, but may not reveal a large population of new objects.

The magnetic field in the Sun's corona stores energy that can be released to heat the coronal plasma and drive solar eruptions. Measurements of the global coronal magnetic field have been limited to a few snapshots. We present observations using the Upgraded Coronal Multi-channel Polarimeter, which provided 114 magnetograms of the global corona above the solar limb spanning approximately eight months. We determined the magnetic field distributions at different solar radii in the corona, and monitored the evolution at different latitudes over multiple solar rotations. We found varying field strengths from <1 to 20 Gauss within 1.05-1.6 solar radii and a signature of active longitude in the coronal magnetic field. Coronal models are generally consistent with the observational data, with larger discrepancies in high-latitude regions.

We present the results of a search for gravitational-wave transients associated with core-collapse supernova SN 2023ixf, which was observed in the galaxy Messier 101 via optical emission on 2023 May 19th, during the LIGO-Virgo-KAGRA 15th Engineering Run. We define a five-day on-source window during which an accompanying gravitational-wave signal may have occurred. No gravitational waves have been identified in data when at least two gravitational-wave observatories were operating, which covered $\sim 14\%$ of this five-day window. We report the search detection efficiency for various possible gravitational-wave emission models. Considering the distance to M101 (6.7 Mpc), we derive constraints on the gravitational-wave emission mechanism of core-collapse supernovae across a broad frequency spectrum, ranging from 50 Hz to 2 kHz where we assume the GW emission occurred when coincident data are available in the on-source window. Considering an ellipsoid model for a rotating proto-neutron star, our search is sensitive to gravitational-wave energy $1 \times 10^{-5} M_{\odot} c^2$ and luminosity $4 \times 10^{-5} M_{\odot} c^2/\text{s}$ for a source emitting at 50 Hz. These constraints are around an order of magnitude more stringent than those obtained so far with gravitational-wave data. The constraint on the ellipticity of the proto-neutron star that is formed is as low as $1.04$, at frequencies above $1200$ Hz, surpassing results from SN 2019ejj.

The long-term evolution of the outer Solar System is subject to the influence of the giant planets, however, perturbations from other massive bodies located in the region imprint secular signatures, that are discernible in long-term simulations. In this work, we performed an in-depth analysis of the evolution of massive objects Eris, 2015 KH$_{162}$, Pluto, and 2010 EK$_{139}$ (a.k.a. Dziewanna), subject to perturbations from the giant planets and the 34 largest trans-Neptunian objects. We do this by analysing 200, 1 Gyr long simulations with identical initial conditions, but requiring the numerical integrator to take different time steps for each realization. Despite the integrator's robustness, each run's results are surprisingly different, showing the limitations of individual realizations when studying the trans-Neptunian region due to its intrinsic chaotic nature. For each object, we find orbital variables with well-defined oscillations and limits, and others with surprisingly large variances and seemingly erratic behaviors. We found that 2015 KH$_{162}$ is a non-resonant and very stable object that experiences only limited orbital excursions. Pluto is even more stable and we found a new underlying constraining mechanism for its orbit; 2010 EK$_{139}$ is not well trapped in the 7:2 mean motion resonance in the long-term and cannot be trapped simultaneously in von-Zeipel-Lidov-Kozai resonance; and finally, we found that at present Eris's longitude of perihelion is stationary, tightly librating around 190$^\circ$, but unexpectedly loses its confinement, drifting away after 150 Myr, suggesting a missing element in our model.

New line lists for four isotopologues of nickel monohydride, $^{58}$NiH, $^{60}$NiH, $^{62}$NiH, and $^{58}$NiD are presented covering the wavenumber range $<10000$ cm$^{-1}$ ($\lambda > 1$ $\mu$m), $J$ up to 37.5 for transitions within and between the three lowest-lying electronic states, ${X}\,^{2}\Delta$, ${W}\,^{2}\Pi$, and ${V}\,^{2}\Sigma^{+}$. The line lists are applicable for temperatures up to 5000 K. The line lists calculations are based on a recent empirical NiH spectroscopic model [Havalyova et al. J. Quant. Spectrosc. Radiat. Transf., 272, 107800, (2021)] which is adapted for the variational nuclear-motion code Duo. The model consists of potential energy curves, spin-orbit coupling curves, electronic angular momentum curves, spin-rotation coupling curves, $\Lambda$-doubling correction curve for $^2\Pi$ states and Born-Oppenheimer breakdown (BOB) rotational correction curves. New ab initio dipole moment curves, scaled to match the experimental dipole moment of the ground state, are used to compute Einstein A coefficients. The BYOT line lists are included in the ExoMol database at this http URL.

We model the power spectrum of galaxies in the IllustrisTNG simulation when they are weighted by their star formation rate. Such a weighting is relevant in the context of line-intensity mapping (LIM). On intermediate to large scales, the model accounts for non-linear bias of star-forming galaxies and halo exclusion (a 2-halo term). On small scales, it incorporates the weighted distribution of satellite galaxies within haloes (a 1-halo term). The random sampling of satellite galaxies adds a shot noise term to the power spectrum on small scales, and their confinement to haloes introduces a halo shot noise term on large scales. The full model reproduces the measured power spectrum to within a few per cent on all scales, and the fitting parameters have a clear physical meaning. Omitting satellite galaxies from the analysis leads to an underestimation of both the large-scale bias and the mean LIM intensity by approximately 30 per cent each at redshift 1.5. Assigning the LIM intensity of satellites to the centre of their respective haloes affects the power spectrum on scales $k > 0.3$ h Mpc$^{-1}$. We discuss how the LIM power spectrum can be used to constrain cosmology on large scales, and galaxy formation on smaller scales, with our fitting function providing an accurate and well-motivated parametrisation.

We study the sensitivity of presupernova evolution and supernova nucleosynthesis yields of massive stars to variations of the initial composition. We use the solar abundances from Lodders (2009), and compute two different initial stellar compositions: i) scaled solar abundances, and ii) the isotopic galactic chemical history model (GCH) developed by West and Heger (2013b). We run a grid of models using the KEPLER stellar evolution code, with 7 initial stellar masses, 12 initial metallicities, and two for each scaling method to explore the effects on nucleosynthesis over a metallicity range of $-4.0\leq[Z]\leq+0.3$. We find that the compositions from the GCH model better reproduce the weak \emph{s}-process peak than the scaled solar models. The model yields are then used in the OMEGA Galactic Chemical Evolution (GCE) code to assess this result further. We find that initial abundances used in computing stellar structure have more of an impact on GCE results than initial abundances used in the burn network, with the GCH model again being favored when compared to observations. Lastly, a machine learning algorithm was used to verify the free parameter values of the GCH model, which were previously found by West and Heger (2013b) using a stochastic fitting process. The updated model is provided as an accessible tool for further nucleosynthesis studies.

Here we present several use cases for using Generative AI (Gen AI) to improve systems engineering and cognitive knowledge management related to the future of astronomy from a culmination of working meetings and presentations as part of the Gen AI Task Group for the NASA Habitable Worlds Observatory (HWO) Science and Technology Architecture Review Team (START) AI/ML Working Group. Collectively, our group mission statement is "Where is the Human-in-the-loop as Gen AI systems become more powerful and autonomous?" with an emphasis on the ethical applications of Gen AI, guided by using these systems to remove drudgery from human work while simultaneously increasing opportunities for humans to experience more collective creativity and innovation. The HWO mission stands to benefit dramatically from generative models for different data types including text, time series/spectra, and image data. These cover a wide range of applications in science and engineering for HWO, including: mission development acceleration, data analysis and interpretation, enhancing imaging capabilities, anomaly detection, predictive modeling and simulation, data augmentation for machine learning, instrument calibration and optimization, public engagement and education, and assisting in mission planning. As an example, through sensitivity analysis of simulated exoplanet population science data sets of various generative model complexity, we can reverse engineer the measurement uncertainty requirements for HWO instruments to produce data that can constrain population models and thus inform HWO design requirements. This approach to HWO design is one example of a strategy that can ensure that HWO remains AI-ready. Through presenting herein a combination of visionary ideas balanced with grounded validated use case examples, we aim to support the development of a long-term strategy to keep HWO AI-ready as it moves forward.

The spatial-temporal evolution of coronal plasma parameters of the solar outer atmosphere at global scales, derived from solar full-disk imaging spectroscopic observation in the extreme-ultraviolet band, is critical for understanding and forecasting solar eruptions. We propose a multi-slits extreme ultraviolet imaging spectrograph for global coronal diagnostics with high cadence and present the preliminary instrument designs for the wavelength range from 18.3 to 19.8 nm. The instrument takes a comprehensive approach to obtain global coronal spatial and spectral information, improve the detected cadence and avoid overlapping. We first describe the relationship between optical properties and structural parameters, especially the relationship between the overlapping and the number of slits, and give a general multi-slits extreme-ultraviolet imaging spectrograph design process. Themultilayer structure is optimized to enhance the effective areas in the observation band. Five distantly-separated slits are set to divide the entire solar field of view, which increase the cadence for raster scanning the solar disk by 5 times relative to a single slit. The spectral resolving power of the optical system with an aperture diameter of 150 mm are optimized to be greater than 1461. The spatial resolution along the slits direction and the scanning direction are about 4.4''and 6.86'', respectively. The Al/Mo/B4C multilayer structure is optimized and the peak effective area is about 1.60 cm2 at 19.3 nm with a full width at half maximum of about 1.3 nm. The cadence to finish full-disk raster scan is about 5 minutes. Finally, the instrument performance is evaluated by an end-to-end calculation of the system photon budget and a simulation of the observational image and spectra. Our investigation shows that this approach is promising for global coronal plasma diagnostics.

We report the HST WFC3 G141 grism slitless spectroscopy observation of the core region of the Spiderweb protocluster at $z=2.16$. We analyzed the spectra of all objects in a $\sim 2 \times 2 \text{ arcmin}^2$ field of view and identified 40 protocluster members, recovering 19 previously identified H$\alpha$-emitters in addition to revealing 21 new members. The spectra allowed us to identify 11 galaxies with quiescent spectra. Three galaxies with quiescent spectra are possibly still star-forming according to SED fitting, indicating a possible left-over or dust-obscured star formation. We estimate a quiescent fraction of $\sim 50\%$ for $M_\star > 10^{11} M_\odot$. About half of the quiescent galaxies possibly host AGN, hinting at AGN's key role in quenching galaxies in the protocluster environment. These quiescent galaxies have relatively more compact and concentrated light profiles than the star-forming members, but they are not yet as bulge-dominated as local ellipticals. These results are consistent with previous studies that indicate the Spiderweb protocluster is in the maturing stage, with a red sequence that has begun forming.

We analyzed the global and resolved properties of approximately 1,240 nearby star-forming galaxies from the MaNGA survey, comparing compact and extended galaxies -- those with smaller and larger radii ($R_{\rm e}$), respectively -- at a fixed stellar mass ($M_{\ast}$). Compact galaxies typically exhibit lower HI gas fractions, higher dust extinction, higher metallicity, greater mass concentration, and lower angular momentum on a global scale. Radial profiles of stellar mass surface density ($\Sigma_{\ast}$) and star formation rate surface density ($\Sigma_{\rm SFR}$), as functions of the effective radius ($R/$$R_{\rm e}$), reveal that compact galaxies display steeper gradients and higher values, resulting in elevated specific star formation rates (sSFR) in their inner regions compared to their outskirts. At a given $\Sigma_{\ast}$, compact galaxies have higher sSFR than extended galaxies, particularly in low-mass galaxies (log($M_{\ast}$/$M_{\odot}$)$\,\leq\,$10$^{10}$). Additionally, their metallicity profiles differ significantly: extended galaxies have steeper metallicity gradients, while compact galaxies exhibit flatter slopes and higher metallicity at a given $R/$$R_{\rm e}$. After accounting for the dependence of metallicity on $M_{\ast}$ and $\Sigma_{\ast}$, no further correlation with SFR is observed. The combination of higher sSFR and potentially higher star formation efficiency in compact galaxies suggests that their central gas is being rapidly consumed, leading to older stellar populations, as indicated by D$_{n}$(4000) and EW(H$\delta_A$), and resulting in faster central growth. Our results reveal that radial SFR profiles cannot be fully determined by $M_{\ast}$ and $\Sigma_{\ast}$ alone; other factors, such as galaxy size or angular momentum, must be considered to fully understand the observed trends.

We calculate the scalar-induced one-loop correction to the power spectrum of tensor perturbations produced during single-field slow-roll inflation. We find that the correction is given by the square of the product of the slow-roll parameter and the tree-level scalar power spectrum. We also discuss the implications of the logarithmic contribution.

We constrain the fundamental-mode ($f$-mode) oscillation frequencies of nonrotating neutron stars using a phenomenological Gaussian process model for the unknown dense-matter equation of state conditioned on a suite of gravitational-wave, radio and X-ray observations. We infer the quadrupolar $f$-mode frequency preferred by the astronomical data as a function of neutron star mass, with error estimates that quantify the impact of equation of state uncertainty, and compare it to the contact frequency for inspiralling neutron-star binaries, finding that resonance with the orbital frequency can be achieved for the coalescences with the most unequal mass ratio. For an optimally configured binary neutron star merger, we estimate the gravitational waveform's tidal phasing due to $f$-mode dynamical tides as $7^{+2}_{-3}$ rad at merger. We assess prospects for distinguishing $f$-mode dynamical tides with current and future-generation gravitational-wave observatories.

The discrepancies between theoretical and observed spectra, and the systematic differences between various spectroscopic parameter estimates, complicate the determination of atmospheric parameters of M-type stars. In this work, we present an empirical sample of 5105 M-type star spectra with homogeneous atmospheric parameter labels through stellar-label transfer and sample cleaning. We addressed systematic discrepancies in spectroscopic parameter estimates by adopting recent results for Gaia EDR3 stars as a reference standard. Then, we used a density-based spatial clustering of applications with noise to remove unreliable samples in each subgrid of parameters. To confirm the reliability of the stellar labels, a 5-layer neural network was utilized, randomly partitioning the samples into training and testing sets. The standard deviations between the predicted and actual values in the testing set are 14 K for Teff , 0.06 dex for log g, and 0.05 dex for [M/H], respectively. In addition, we conducted an internal cross-validation to enhance validation and obtained precisions of 11 K, 0.05 dex, and 0.05 dex for Teff , log g, and [M/H], respectively. A grid of 1365 high Signal-to-Noise ratio (S/N) spectra and their labels, selected from the empirical sample, was utilized in the stellar parameter pipeline for M-Type stars (LASPM) of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), producing an almost seamless Kiel distribution diagram for LAMOST DR10 and DR11 data. The atmospheric parameters for M-type stars from LAMOST DR11 show improved precision compared to the data from DR9, with improvements (for spectra with S/N higher than 10) from 118 to 67 K in Teff , 0.2 to 0.07 dex in log g, and 0.29 to 0.14 dex in [M/H].

The Ly$\alpha$ forest, a series of HI absorption lines in the quasar spectra, is a powerful tool for probing the large-scale structure of the intergalactic medium. Its three-dimensional (3D) correlation and cross-correlations with quasars allow precise measurements of the baryon acoustic oscillation feature and redshift space distortions at redshifts $z>2$. Understanding small-scale astrophysical phenomena, such as star formation and feedback, is crucial for full-shape analyses. In this study, we measure the 3D auto-power spectrum of the Ly$\alpha$ forest and its cross-power spectrum with halos using hydrodynamic simulations from the GADGET3-OSAKA code, which includes models for star formation and supernova feedback. Across five astrophysical models, we find significant deviations from the Fiducial model, with $5-10\,\%$ differences for wavenumbers $k>2\,h\mathrm{Mpc}^{-1}$ in the Ly$\alpha$ auto-power spectrum. The Ly$\alpha\,\times\,$halo cross-power spectra show even larger deviations, exceeding $10\,\%$ in some cases. Using the fitting models of Arinyo-i-Prats et al. (2015) and Givans et al. (2022), we jointly fit the Ly$\alpha$ auto- and Ly$\alpha$ $\times$ halo cross-power spectra, and assess the accuracy of the estimated $f\sigma_8$ parameter by comparing it with the ground truth from the simulations, while varying the maximum wavenumber $k_\mathrm{max}$ and minimum halo mass $M_h$. Our results demonstrate that the extended model of Givans et al. (2022) is highly effective in reproducing $f\sigma_8$ at $k_\mathrm{max}\leq3.0\,h\mathrm{Mpc}^{-1}$ for $M_h>10^{10.5} M_\odot$, and remains robust against astrophysical uncertainties.

The light curves and spectra of many Type I and Type II supernovae (SNe) are heavily influenced by the interaction of the SN ejecta with circumstellar material (CSM) surrounding the progenitor star. The observed diversity shows that many progenitors have undergone some level of stripping and polluted their CSM shortly before the explosion. The presence of a binary companion and the mass transfer that can ensue offers a mechanism that can give rise to this diversity. We present a set of detailed massive evolutionary models in which the donor star, a Red Supergiant (RSG) is in a wide orbit around a main-sequence companion, and undergoes stable or unstable mass transfer in the later stages of evolution, up to the moment of core collapse. We also discuss some significant physics of these systems that may impact our results, from the presence of pulsations and extended atmospheres in RSGs, to the initial eccentricity of the orbit. The resulting SN types range from Type IIP to H-deficient IIb and H-free Ib. In models undergoing stable mass transfer, the material lost during this process is expected to form a dense CSM surrounding the system by the time of core collapse and give rise to significant interaction effects in the SN light curve and spectra. In the systems with unstable mass transfer mass transfer, the SN may occur during common-envelope evolution. In this case, the progenitor may show significant variability in the last few thousand years before core collapse, and the following SN will likely exhibit strong interaction effects.

We investigate the use of dust-scattering rings from GRB 221009A to measure cosmic curvature. We derive the relationship between scattering angle and time delay in non-flat universes and attempt to constrain cosmological parameters by fitting theoretical predictions to observational data. The results show that this method is limited in its effectiveness due to the small geometric scale of observations. While unable to provide significant constraints, this study introduces a novel approach to cosmological measurements using GRB scattering rings. We discuss the method's limitations and potential future applications with larger-scale observations.

Since the launch of the Solar Dynamics Observatory (SDO) in 2010 and throughout the solar cycle 24, the Sun has produced few tens of xclass flares, which are the most energetic solar events. Those flares are produced in regions where the magnetic flux/energy is large and the magnetic configurations are complex. To provide more insights into the flaring process, we investigate the properties of magnetic null points (MNPs) and their correlation with the energy release sites. During solar cycle 24, we identify 17 xclass flares satisfying selection criteria. From SDO/HMI magnetograms, we perform potential extrapolations around the peak time of the flare to access the 3D coronal magnetic field and thus investigate the existence of coronal MNPs. We then correlate the flaring sites with the existing MNPs using SDO/AIA 171A EUV observations, and deduce their properties (sign, spine, fan). Six active regions out of 10 possess at least one MNP, which is stable and with large magnetic field gradients: this implies that 35% of xclass flares are associated with a MNP; of which 87.5% of MNPs are of positive type. The MNPs associated with the flare sites are predominantly located ata height between 05 and 2 Mm, and with a vertical/radial spine field line. We also find a slight correlation between the MNPs not associated with a flare and negative-type MNPs (55%) within the active region. Regarding the physics of flares, the association between the enhanced intensity at the flaring site and a MNP represents about a third of the possible scenarios for triggering xclass flares.

IRAS 22568+6141 has been classified as a low-ionisation planetary nebula (PN) and presents non-thermal radio continuum emission, which could be a signature of nascent PNe. We present intermediate-resolution long-slit spectra obtained in 2021 and 2023, high-resolution long-slit spectra taken in 2023, and a light curve at the $r$-filter between 1953 and 2019, that reveal changes in IRAS 22568+6141 with timescales of decades and a few years. The object underwent an energetic event around 1990 that suddenly increased its brightness which has been fading since then. A comparison with a published spectrum from 1988 shows an increase of the H$\beta$ flux in 2021 by factor of $\simeq$6 and the [O III] emission lines that were absent in 1988. Between 2021 and 2023 the H$\beta$ flux decreased by a factor of $\simeq$1.7, and the [O III] emission lines almost vanished. These results and the variability observed in other emission lines indicate that IRAS 22568+6141 is recombining and cooling down between 2021 and 2023, and probably since 2005, as suggested by archival radio continuum and mid-IR observations. The intermediate- and high-resolution spectra show that the excitation of the emission lines is dominated by shocks in 2021 and 2023, and, probably, also in 1988, which may be related to the non-thermal radio continuum emission from the object. Although the variability might be due to changes in the physical conditions in the shocks or in a nova-like eruption, it accommodates better to that expected from a late thermal pulse, which is further suggested by a comparison with other similar objects. New observations and monitoring in the coming years are crucial to corroborate the origin of the variability.

We report the detection of an extreme flux decrease accompanied by clear dispersion measure (DM) and rotation measure (RM) variations for pulsar B1929+10 during the 110-minute radio observation with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The radio flux decreases by 2 to 3 orders of magnitude within a rapid time scale of about 20 minutes. Meanwhile, the variations of DM and RM are approximately 0.05 pc cm$^{-3}$ and 0.7 rad m$^{-2}$, respectively. Frequency-dependent analysis of DM indicates an extremely weak chromatic DM feature, which does not notably affect the radiative behavior detected. Moreover, the pulsar timing analysis shows an additional time delay from 100 $\mu$s to 400 $\mu$s in the event. These results are speculated to be due to the eclipse and bend for the radio emission of pulsar B1929+10 by a highly dense outflow from the pulsar. This not only impacts the intrinsic radio emission feature but also affects the pulsar timing behavior. Nevertheless, a plasma lens effect lasting around 20 minutes could also be responsible for the event.

We present a detailed analysis of the thermal X-ray emission from the intracluster medium (ICM) in the cool-core galaxy cluster Abell 2667 ($z=0.23$). Our goal is to detect low-temperature ($<2$ keV) X-ray emitting gas, potentially associated to a cooling flow that connects the hot ICM reservoir to the cold gas phase responsible for star formation and supermassive black hole feeding. We use new deep XMM-Newton EPIC and RGS data, combined with archival Chandra data, to perform a spectral analysis for the core region. We find 1$\sigma$ upper limits to the cooling gas fraction of $\sim$40 $\rm M_{\odot}yr^{-1}$ and $\sim$50-60 $\rm M_{\odot}yr^{-1}$ in the temperature ranges 0.5-1 keV and 1-2 keV, respectively. The lack of OVII, FeXXI-FeXXII, and FeXVII emission lines in the RGS spectra suggest that the fraction of gas cooling below 1 keV is limited to a few tens of $\rm M_{\odot}yr^{-1}$ at most. However, we detect several lines (e.g. SiXIV, MgXII, FeXXIII/FeXXIV, NeX, OVIII$\alpha$) that allow us to estimate a 1$\sigma$ upper limit for turbulent broadening of $\sim$320 km $\rm s^{-1}$, higher that other cool-core clusters such as Abell 1835, implying mechanisms that boost turbulence in Abell 2667's atmosphere. Imaging analysis of Chandra data suggests the presence of a cold front, possibly lined to sloshing or ICM cavities. However, current data do not clearly identify the physical mechanism driving turbulence. These finding indicate that Abell 2667 is similar to other low-redshift cool-core clusters, though the large upper limit on turbulence hints at significant ICM heating, which may suppress cooling for extended periods and contribute to future condensation events.

We simulate the dynamics of about 600 quantum vortices in a spinning-down cylindrical container using a Gross--Pitaevskii model. For the first time, we find convincing spatial-temporal evidence of avalanching behaviour resulting from vortex depinning and collective motion. During a typical avalanche, about 10 to 20 vortices exit the container in a short period, producing a glitch in the superfluid angular momentum and a localised void in the vorticity. After the glitch, vortices continue to depin and circulate around the vorticity void in a similar manner to that seen in previous point-vortex simulations. We present evidence of collective vortex motion throughout this avalanche process. We also show that the effective Magnus force can be used to predict when and where avalanches will occur. Lastly, we comment on the challenge of extrapolating these results to conditions in real neutron stars, which contain many orders of magnitude more vortices.

Matter falling onto black holes, {also called} accretion discs, emit intense high-energy radiation. Accretion discs during {hard to hard intermediate} spectral states also emit bipolar outflows. Radiation drag was supposed to impose the upper limit on the terminal speed. It was later shown that a radiation field around an advective accretion disc imposes no upper limit on speed, about a few hundred of Schwarzschild radius from the disc surface. We {study radiatively driven electron-proton and electron-positron jets, for gemeotrically thick and slim transonic discs} by using numerical simulation. We show that pair-dominated jets can reach ultra-relativistic speeds by radiation driving. We also discuss at what limits radiative acceleration may fail.

In this catalog, we present the results of a systematic study of 199 short gamma-ray bursts (GRBs) detected by Konus-Wind between 2011 January 1 and 2021 August 31. The catalog extends the Second Catalog of short gamma-ray bursts covering the period 1994-2010 by ten years of data. The resulting Konus-Wind short GRB sample includes 494 bursts. From temporal and spectral analyses of the sample, we provide the burst durations, spectral lags, estimates of the minimum variability time scales, rise and decay times, the results of spectral fits with three model functions, the total energy fluences, and the peak energy fluxes of the bursts. We present statistical distributions of these parameters for the complete set of 494 short gamma-ray bursts detected in 1994-2021. We discuss evidence found for an additional spectral component and the presence of extended emission in a fraction of the short GRBs. Finally, we consider the results in the context of the Type I (merger-origin)/Type II (collapsar-origin) classification, and discuss magnetar giant flare sub-sample.

A recent paper (King, 2024) suggested that emission from the central supermassive black holes in high-redshift galaxies must be tightly collimated by the effects of partly expelling a super-Eddington mass supply. I show here that this idea predicts that these galaxies should produce very little detectable rest-frame X-ray emission, appear Compton thick, and show no easily detectable sign of outflows. All of these properties agree with current observations. To produce these effects, the mass supply to the black holes should exceed the Eddington rate by factors 50 - 100, which appears in line with conditions during the early growth of the holes. I note that theoretical derivations of the ratio of black hole mass to host galaxy stellar mass already predict that this should increase significantly at high redshift, in line with recent observations.

Context. Terzan 5 and Liller 1 are the only bulge stellar clusters hosting multi-iron and multi-age stellar populations. They are therefore claimed to constitute a novel class of astrophysical objects: the fossils of massive star-forming clumps that possibly sank to the center of the Milky Way and contributed to the formation of the bulge. This is based on the hypothesis that the ancient clumps were able to retain iron-enriched supernova ejecta, later giving rise to younger and more metal-rich populations. Aims. A way to investigate this scenario is reconstructing their star formation histories (SFHs) and proving a prolonged and multi-episode star formation activity. Methods. Leveraging ground- and space-based high-resolution images, we derived the SFH of Terzan 5 by employing the color-magnitude diagram fitting routine SFERA. Results. The best-fit solution predicts an old, main peak occurred between 12 and 13 Gyr ago that generated 70 % of the current stellar mass, followed by a lower-rate star formation activity with two main additional bursts. Conclusions. These results indicate that Terzan 5, similarly to Liller 1, experienced a prolonged, multiepisode star formation activity, fueled by metal-enriched gas deposited in its central regions, in agreement with the expectations of a self-enrichment scenario in a primordial massive clump.

Context. Transient sky astronomy is entering a new era with the advent of the SVOM mission (Space Variable Objects Monitor), which was successfully launched on the 26th of June, 2024. The primary goal of SVOM is to monitor the hard X-ray sky searching for gamma-ray bursts (GRBs). On top of its on-board follow-up capabilities, SVOM will be backed by its ground segment composed of several facilities, of which the near-infrared imager CAGIRE. Mounted on the robotic telescope COLIBRI, it will be a unique instrument, able to perform fast follow-up of GRB afterglows in J and H bands, an ideal combination to catch high-redshift (z>6) and/or obscured GRBs. Aims. This paper aims at estimating the performances of CAGIRE for GRB near-infrared afterglow detection based on the characteristics of the detector and the specificities of the COLIBRI telescope. Quickly fading GRB afterglows pose challenges that should be addressed by adapting observing strategies to the capabilities of CAGIRE. Methods. We use an end-to-end image simulator to produce realistic CAGIRE images, taking into account results from the characterization of the ALFA detector used by CAGIRE. We implemented a GRB afterglow generator that simulates infrared lightcurves and spectra based on published observation of distant GRBs (z>6). Results. We retrieved the photometry of 9 GRB afterglows in various scenarios covered by CAGIRE. Catching afterglows as early as two minutes after burst allows the identification of a nIR counterpart in the brightest 4 events. When artificially redshifted even further away, these events remain detectable by CAGIRE up to z=9.6 in J band, and z=13.3 in H band, indicating the potential of CAGIRE to be a pioneer in the identification of the most distant GRBs to date.

Rocky exoplanets are faint and difficult to observe due to their small size and low brightness compared to their host star. Despite this, the James Webb Space Telescope (JWST) has allowed us new methods and opportunities to study them. Accurately characterising exoplanet atmospheres could offer insights not only into the planetary demographics of rocky worlds in the universe, but also how our own Earth compares. Previous work simulating the observational spectra of planets with various models of atmospheric composition has constrained conditions under which specific molecules would be observable. We seek a different approach that does not depend on specific molecular features, but rather on correlating "agnostic" spectral characteristics to each other. This study uses 42 synthetic transit transmission spectra in a range matching that of JWST's MIRI instrument. Spectra were de-noised, normalised and split into bands for which the enclosed areas were calculated. Using the HDBSCAN clustering algorithm, separation patterns and clusters were found across these bands, which correlated to various a priori planet characteristics, such as relative gas concentration. Using 0.7 $\mu$m bands between 7 and 12 $\mu$m, we are able to separate reducing vs. oxic atmospheres and constrain CO$_2$ and O$_2$ mixing ratios accurate to within two orders of magnitude in an unknown sample. Through this, we have demonstrated that it is possible to differentiate and analyse exoplanet atmosphere type from agnostic interrogation of their transit spectra with a method we call "Molecule Agnostic Spectral Clustering".

The observed dust rings and gaps in protoplanetary disks could be imprints of forming planets. Even low-mass planets in the one-to-ten Earth-mass regime, that do not yet carve deep gas gaps, can generate such dust rings and gaps by driving a radially-outwards gas flow, as shown in previous work. However, understanding the creation and evolution of these dust structures is challenging due to dust drift and diffusion, requiring an approach beyond previous steady state models. Here we investigate the time evolution of the dust surface density influenced by the planet-induced gas flow, based on post-processing three-dimensional hydrodynamical simulations. We find that planets larger than a dimensionless thermal mass of $m=0.05$, corresponding to 0.3 Earth mass at 1 au or 1.7 Earth masses at 10 au, generate dust rings and gaps, provided that solids have small Stokes numbers} (${\rm St}\lesssim10^{-2}$) and that the disk midplane is weakly turbulent ($\alpha_{\rm diff}\lesssim10^{-4}$). As dust particles pile up outside the orbit of the planet, the interior gap expands with time, when the advective flux dominates over diffusion. Dust gap depths range from a factor a few, to several orders of magnitude, depending on planet mass and the level of midplane particle diffusion. We construct a semi-analytic model describing the width of the dust ring and gap, and then compare it with the observational data. We find that up to 65\% of the observed wide-orbit gaps could be explained as resulting from the presence of a low-mass planet, assuming $\alpha_{\rm diff}=10^{-5}$ and ${\rm St}=10^{-3}$. However, it is more challenging to explain the observed wide rings, which in our model would require the presence of a population of small particles (${\rm St=10^{-4}}$). Further work is needed to explore the role of pebble fragmentation, planet migration, and the effect of multiple planets.

Fast radio bursts (FRBs) are sufficiently energetic to be detectable from luminosity distances up to at least seven billion parsecs (redshift $z > 1$). Probing the maximum energies and luminosities of FRBs constrains their emission mechanism and cosmological population. Here we investigate the maximum energetics of a highly active repeater, FRB 20220912A, using 1,500h of observations. We detect $130$ high-energy bursts and find a break in the burst energy distribution, with a flattening of the power-law slope at higher energy. This is consistent with the behaviour of another highly active repeater, FRB 20201124A. Furthermore, we model the rate of the highest-energy bursts and find a turnover at a characteristic spectral energy density of $E^{\textrm{char}}_{\nu} = 2.09^{+3.78}_{-1.04}\times10^{32}$ erg/Hz. This characteristic maximum energy agrees well with observations of apparently one-off FRBs, suggesting a common physical mechanism for their emission. The extreme burst energies push radiation and source models to their limit.

We present the localization and host galaxy of FRB 20190208A, a repeating source of fast radio bursts (FRBs) discovered using CHIME/FRB. As part of the PRECISE repeater localization program on the EVN, we monitored FRB 20190208A for 65.6 hours at $\sim1.4$ GHz and detected a single burst, which led to its VLBI localization with 260 mas uncertainty (2$\sigma$). Follow-up optical observations with the MMT Observatory ($i\gtrsim 25.7$ mag (AB)) found no visible host at the FRB position. Subsequent deeper observations with the GTC, however, revealed an extremely faint galaxy ($r=27.32 \pm0.16$ mag), very likely ($99.95 \%$) associated with FRB 20190208A. Given the dispersion measure of the FRB ($\sim580$ pc cm$^{-3}$), even the most conservative redshift estimate ($z_{\mathrm{max}}\sim0.83$) implies that this is the lowest-luminosity FRB host to date ($\lesssim10^8L_{\odot}$), even less luminous than the dwarf host of FRB 20121102A. We investigate how localization precision and the depth of optical imaging affect host association, and discuss the implications of such a low-luminosity dwarf galaxy. Unlike the other repeaters with low-luminosity hosts, FRB 20190208A has a modest Faraday rotation measure of a few tens of rad m$^{-2}$, and EVN plus VLA observations reveal no associated compact persistent radio source. We also monitored FRB 20190208A for 40.4 hours over 2 years as part of the ÉCLAT repeating FRB monitoring campaign on the Nançay Radio Telescope, and detected one burst. Our results demonstrate that, in some cases, the robust association of an FRB with a host galaxy will require both high localization precision, as well as deep optical follow-up.

In this work, we analyzed the orbits of more than 35,000 (for 2024) near-Earth objects (NEOs) for the possibility of successive approaches to all pairs of planets: Earth, Venus, and Mars in the time range from 2020 to 2120. We have selected 120 candidates for Earth-Mars, Earth-Venus, Mars-Earth, Mars-Venus, Venus-Earth, and Venus-Mars fast transfers (within 180 days); 2 candidates for double transfers (consecutive approaches with three planets); 10 candidates for multiple transfers, when an asteroid has several consecutive paired approaches to planets in a hundred years.

We reproduce the luminosity functions of the early-time peak optical luminosity, the late-time UV plateau luminosity, and the peak X-ray luminosity of tidal disruption events, using an entirely first-principles theoretical approach. We do this by first fitting three free parameters of the tidal disruption event black hole mass distribution using the observed distribution of late time UV plateau luminosities, using a time-dependent relativistic accretion model. Using this black hole mass distribution we are then, with no further free parameters of the theory, able to reproduce exactly the peak X-ray luminosity distribution of the tidal disruption event population. This proves that the X-ray luminosity of tidal disruption events are sourced from the same accretion flows which produce the late time UV plateau. Using an empirical scaling relationship between peak optical luminosities and black hole masses, itself calibrated using the same relativistic accretion theory, we are able to reproduce the observed peak optical luminosity function, again with no additional free parameters. Implications of these results include that there is no tidal disruption event "missing energy problem", that the optical and X-ray selected tidal disruption event populations are drawn from the same black hole mass distribution, that the early time optical luminosity in tidal disruption events is somewhat simple, at least on the population level, and that future LSST observations will be able to constrain the black hole mass function at low masses.

The effect of baryon physics associated with galaxy formation onto the large-scale matter distribution of the Universe is a key uncertainty in the theoretical modelling required for the interpretation of Stage IV cosmology surveys. We use the FLAMINGO suite of simulations to study the baryon response due to galaxy formation of the total matter power spectrum. We find that it is only well converged for simulation volumes in excess of $(200~Mpc)^3$. We report results for simulations of varying feedback intensity, which either match the X-ray inferred gas fractions in clusters and the z=0 stellar mass function, or shifted versions of the data, as well as for different implementations of AGN feedback. We package our results in the form of a Gaussian process emulator which can rapidly reproduce all the simulations' predictions to better than 1% up to the comoving wavenumber $k = 10~h/Mpc$ and up to z=2 for all the feedback models present in the FLAMINGO suite. We find that the response becomes stronger, the range of scales affected increases, and the position of the minimum of the response moves to smaller scales as the redshift decreases. We find that lower gas fractions in groups and clusters lead to a stronger response and that the use of collimated jets instead of thermally driven winds for AGN feedback enhances the effect. Lowering the stellar masses at fixed cluster gas fractions also increases the magnitude of the response. We find only a small (1% at $k<10~h/Mpc$) dependence of our results on the background cosmology, but a wider range of cosmology variations will be needed to confirm this result. The response we obtain for our strongest feedback models is compatible with some of the recent analyses combining weak lensing with external data. Such a response is, however, in strong tension with the X-ray inferred gas fractions in clusters used to calibrate the FLAMINGO model.

Context: Theoretical models of early accretion during the formation process of massive stars have predicted that HII regions exhibit radio variability on timescales of decades. However, large-scale searches for such temporal variations with sufficient sensitivity have not yet been carried out. Aims: We aim to identify HII regions with variable radio wavelength fluxes and to investigate the properties of the identified objects, especially those with the highest level of variability. Methods: We compared the peak flux densities of 86 ultracompact HII (UCHII) regions measured by the GLOSTAR and CORNISH surveys and identified variables that show flux variations higher than 30% over ~8 yr timespan between these surveys. Results: We found a sample of 38 variable UCHII regions, which is the largest sample identified to date. The overall occurrence of variability is 44$\pm$5%, suggesting that variation in UCHII regions is significantly more common than prediction. The variable UCHII regions are found to be younger than non-variable UCHII regions, all of them meeting the size criterion of hypercompact (HC) HII regions. We studied the 7 UCHII regions (the ``Top7'') that show the highest variability with variations > 100%. The Top7 variable UCHII regions are optically thick at 4--8 GHz and compact, suggesting they are in a very early evolutionary stage of HCHII or UCHII regions. There is a significant correlation between variability and the spectral index of the radio emission. No dependence is observed between the variations and the properties of the sources' natal clumps traced by submillimeter continuum emission from dust, although variable HII regions are found in clumps at an earlier evolutionary stage.

In an era marked by renewed interest in lunar exploration and the prospect of establishing a sustainable human presence on the Moon, innovative approaches supporting mission preparation and astronaut training are imperative. To this end, the advancements in Virtual Reality (VR) technology offer a promising avenue to simulate and optimize future human missions to the Moon. Through VR simulations, tests can be performed quickly, with different environment parameters and a human-centered perspective can be maintained throughout the experiments. This paper presents a comprehensive framework that harnesses VR simulations to replicate the challenges and opportunities of lunar exploration, aiming to enhance astronaut readiness and mission success. Multiple environments with physical and visual characteristics that reflect those found in interesting Moon regions have been modeled and integrated into simulations based on the Unity graphical engine. We exploit VR to allow the user to fully immerse in the simulations and interact with assets in the same way as in real contexts. Different scenarios have been replicated, from upcoming exploration missions where it is possible to deploy scientific payloads, collect samples, and traverse the surrounding environment, to long-term habitation in a futuristic lunar base, performing everyday activities. Moreover, our framework allows us to simulate human-robot collaboration and surveillance directly displaying sensor readings and scheduled tasks of autonomous agents which will be part of future hybrid missions, leveraging the ROS2-Unity bridge. Thus, the entire project can be summarized as a desire to define cornerstones for human-machine design and interaction, astronaut training, and learning of potential weak points in the context of future lunar missions, through targeted operations in a variety of contexts as close to reality as possible.

Radio interferometry is an observational technique used to study astrophysical phenomena. Data gathered by an interferometer requires substantial processing before astronomers can extract the scientific information from it. Data processing consists of a sequence of calibration and analysis procedures where choices must be made about the sequence of procedures as well as the specific configuration of the procedure itself. These choices are typically based on a combination of measurable data characteristics, an understanding of the instrument itself, an appreciation of the trade-offs between compute cost and accuracy, and a learned understanding of what is considered "best practice". A metric of absolute correctness is not always available and validity is often subject to human judgment. The underlying principles and software configurations to discern a reasonable workflow for a given dataset is the subject of training workshops for students and scientists. Our goal is to use objective metrics that quantify best practice, and numerically map out the decision space with respect to our metrics. With these objective metrics we demonstrate an automated, data-driven, decision system that is capable of sequencing the optimal action(s) for processing interferometric data. This paper introduces a simplified description of the principles behind interferometry and the procedures required for data processing. We highlight the issues with current automation approaches and propose our ideas for solving these bottlenecks. A prototype is demonstrated and the results are discussed.

We present coniferest, an open source generic purpose active anomaly detection framework written in Python. The package design and implemented algorithms are described. Currently, static outlier detection analysis is supported via the Isolation forest algorithm. Moreover, Active Anomaly Discovery (AAD) and Pineforest algorithms are available to tackle active anomaly detection problems. The algorithms and package performance are evaluated on a series of synthetic datasets. We also describe a few success cases which resulted from applying the package to real astronomical data in active anomaly detection tasks within the SNAD project.

Recent research provides compelling evidence that the decay of short-lived radioisotopes (SLRs), such as 26Al, provided the bulk of energy for heating and desiccation of volatile-rich planetesimals in the early Solar System. However, it remains unclear whether the early Solar System was highly enriched relative to other planetary systems with similar formation characteristics. While the Solar System possesses an elevated level of SLR enrichment compared to the interstellar medium, determining SLR enrichment of individual protoplanetary disks observationally has not been performed and is markedly more difficult. We use N-body simulations to estimate enrichment of SLRs in star-forming regions through two likely important SLR sources, stellar winds from massive stars and supernovae. We vary the number of stars and the radii of the star-forming regions and implement two models of stellar wind SLR propagation for the radioisotopes 26Al and 60Fe. We find that for 26Al enrichment the Solar System is at the upper end of the expected distribution, while for the more supernovae dependent isotope 60Fe we find that the Solar System is comparatively very highly enriched. Furthermore, combined with our previous research, these results suggest that the statistical role of 26Al-driven desiccation on exoplanet bulk composition may be underestimated in typical interpretations of the low-mass exoplanet census, and that 60Fe is even less influential as a source of heating than previously assumed.

A primary goal of characterizing exoplanet atmospheres is to constrain planetary bulk properties, such as their metallicity, C/O ratio, and intrinsic heat. However, there are significant uncertainties in many aspects of atmospheric physics, such as the strength of vertical mixing. Here we use PICASO and the photochem model to explore how atmospheric chemistry is influenced by planetary properties like metallicity, C/O ratio, $T_{\rm int}$, $T_{\rm eq}$, and $K_{\rm zz}$ in hydrogen-dominated atmospheres. We vary the $T_{\rm eq}$ of the planets between 400 K- 1600 K, across ``cold",``warm," and `hot" objects. We also explore an extensive range of $T_{\rm int}$ values between 30-500 K, representing sub-Neptunes to massive gas giants. We find that gases like CO and CO$_2$ show a drastically different dependence on $K_{\rm zz}$ and C/O for planets with cold interiors (e.g., sub-Neptunes) compared to planets with hotter interiors (e.g., Jupiter mass planets), for the same $T_{\rm eq}$. We also find that gases like CS and CS$_2$ can carry a significant portion of the S- inventory in the upper atmosphere near $T_{\rm eq}$ $\le$ 600 K, below which SO$_2$ ceases to be abundant. For solar C/O, we show that the CO/CH$_4$ ratio in the upper atmospheres of planets can become $\le$1 for planets with low $T_{\rm eq}$, but only if their interiors are cold ($T_{\rm int}$$\le$100 K). We find that photochemical haze precursor molecules in the upper atmosphere show very complex dependence on C/O, $K_{\rm zz}$, $T_{\rm eq}$, and $T_{\rm int}$ for planets with cold interiors (e.g., sub-Neptunes). We also briefly explore fully coupling PICASO and photochem to generate self-consistent radiative-convective-photochemical-equilibrium models.

The plasma opacity in stars depends mainly on the local state of matter (the density, temperature, and chemical composition at the point of interest), but in supernova ejecta it also depends on the expansion velocity gradient, because the Doppler effect shifts the spectral lines differently in different ejecta layers. This effect is known in the literature as the expansion opacity. The existing approaches to the inclusion of this effect, in some cases, predict different results in identical conditions. In this paper we compare the approaches of Blinnikov (1996) and Friend and Castor (1983) - Eastman and Pinto (1993) to calculating the opacity in supernova ejecta and give examples of the influence of different approximations on the model light curves of supernovae.

Concerning recent published studies exploring the presence or otherwise of a gravitational anomaly at low accelerations in wide binary stars as observed by the {\it Gaia} satellite, Cookson (2024) presents an interesting case. In that study, RMS values of binned relative internal velocities in 1D for wide binaries are compared to Newtonian predictions for that quantity, with the author concluding that the data presented show no indication of any inconsistency with Newtonian expectations. However, the comparison presented is critically flawed, as the Newtonian predictions used refer to wide binaries with mean total masses of 2.0 $M_{\odot}$. This is larger than the 1.56 $M_{\odot}$ value which applies to the data used in said paper. In this short note we correct the error mentioned above and show that the data and error bars as given in Cookson (2024) are in fact inconsistent with Newtonian expectations. Contrary to the assertion by Cookson, the data presented in Cookson (2024) actually show a clear anomaly in the low acceleration gravitational regime, this being consistent with a number of independent studies yielding this same result. Based on the analysis by Cookson (2024), wide binary systems show a clear Milgromian deviation from Newtonian dynamics.

We present a catalog of unobscured QSO candidates in the southern hemisphere from the early interim data of the KMTNet Synoptic Survey of Southern Sky (KS4). The KS4 data covers $\sim2500\,{\rm deg}^{2}$ sky area, reaching 5$\sigma$ detection limits of $\sim$22.1-22.7 AB mag in the $BVRI$ bands. Combining this with available infrared photometric data from the surveys covering the southern sky, we select the unobscured QSO candidates based on their colors and spectral energy distributions (SEDs) fitting results. The final catalog contains 72,964 unobscured QSO candidates, of which only 0.4% are previously identified as QSOs based on spectroscopic observations. Our selection method achieves an 87% recovery rate for spectroscopically confirmed bright QSOs at $z<2$ within the KS4 survey area. In addition, the number count of our candidates is comparable to that of spectroscopically confirmed QSOs from the Sloan Digital Sky Survey in the northern sky. These demonstrate that our approach is effective in searching for unobscured QSOs in the southern sky. Future spectro-photometric surveys covering the southern sky will enable us to discern their true nature and enhance our understanding of QSO populations in the southern hemisphere.

In this study, we present the pulse profile of the unique and the second brightest gamma-ray burst GRB 230307A, and analyze its temporal behavior using a joint GECAM--Fermi/GBM time-resolved spectral analysis. The utilization of GECAM data is advantageous as it successfully captured significant data during the pile-up period of the Fermi/GBM. We investigate the evolution of its flux, photon fluence, photon flux, peak energy, and the corresponding hardness-intensity and hardness-flux correlations. The findings within the first 27 seconds exhibit consistent patterns reported previously, providing valuable insights for comparing observations with predictions from the synchrotron radiation model invoking an expanding shell. Beyond the initial 27 seconds, we observe a notable transition in the emitted radiation, attributed to high latitude emission (HLE), influenced by the geometric properties of the shells and the relativistic Doppler effects. By modeling the data within the framework of the large-radius internal shock model, we discuss the required parameters as well as the limitations of the model. We conclude that a more complicated synchrotron emission model is needed to fully describe the observational data of GRB 230307A.

The explosion mechanism of a core-collapse supernova is a complex interplay between neutrino heating and cooling (including the effects of neutrino-driven convection), the gravitational potential, and the ram pressure of the infalling material. To analyze the post-bounce phase of a supernova, one can use the generalized Force Explosion Condition (FEC+), which succinctly formalizes the interplay among these four phenomena in an analytical condition, consistent with realistic simulations. In this paper, we use the FEC+ to study the post-bounce phase of 341 spherically symmetric simulations, where convection is included through a time-dependent mixing length approach. We find that the accretion of the Si/O interface through the expanding shock can significantly change the outcome of the supernova by driving the FEC+ above the explosion threshold. We systematically explore this by (i) artificially smoothing the pre-supernova density profile, and (ii) artificially varying the mixing length. In both cases, we find that large-enough density contrasts at the Si/O interface lead to successful shock revival only if the FEC+ is already close to the explosion threshold. Furthermore, we find that the accretion of the Si/O interface has a substantial effect on the critical condition for supernova explosions, contributing between 5\% and 15\%, depending on how pronounced the density contrast at the interface is. Earlier studies showed that convection affects the critical condition by 25--30\%, which demonstrates that the accretion of the Si/O interface through the shock can play a nearly comparable role in influencing shock dynamics.

We investigate the 1D plane-parallel front connecting the warm ($10^4$ K) and hot ($10^6$ K) phases of the circumgalactic medium (CGM), focusing on the influence of cosmic rays (CRs) in shaping these transition layers. We find that cosmic rays dictate the thermal balance while other fluxes (thermal conduction, radiative cooling, and gas flow) adjust to compensate. We compute column density ratios for selected transition temperature ions and compare them with observational data. While most of our models fail to reproduce the observations, a few are successful, although we make no claims for their uniqueness. Some of the discrepancies may indicate challenges in capturing the profiles in cooler, photoionized regions, as has been suggested for by previous efforts to model thermal transition layers.