Papers for Wednesday, Jun 26 2024

We propose a multi-mode bar consisting of mass elements of decreasing size for the direct detection of stimulated absorption of up to kHz gravitons from a neutron star merger and post-merger. We find that the multi-mode detector has normal modes that retain the coupling strength to the gravitational wave of the largest mass-element, while only having an effective mass comparable to the mass of the smallest element. This allows the normal modes to have graviton absorption rates due to the tonne-scale largest mass, while the single graviton absorption process in the normal mode could be resolved through energy measurements of a mass-element in-principle smaller than pico-gram scale. This improves the transduction of the single-graviton signal compared to a single-mode detector, enhancing the feasibility of detecting single gravitons.

The cosmic 21 cm signal serves as a crucial probe for studying the evolutionary history of the Universe. However, detecting the 21 cm signal poses significant challenges due to its extremely faint nature. To mitigate the interference from the Earth's radio frequency interference (RFI), the ground and the ionospheric effects, the Discovering the Sky at the Longest Wavelength (DSL) project will deploy a constellation of satellites in Lunar orbit, with its high-frequency daughter satellite tasked with detecting the global 21 cm signal from cosmic dawn and reionization era (CD/EoR). We intend to employ the Vari-Zeroth-Order Polynomial (VZOP) for foreground fitting and subtracting. We have studied the effect of thermal noise, thermal radiation from the Moon, the Lunar reflection, anisotropic frequency-dependent beam, inaccurate antenna beam pattern, and RFI contamination. We discovered that the RFI contamination can significantly affect the fitting process and thus prevent us from detecting the signal. Therefore, experimenting on the far side of the moon is crucial. We also discovered that using VZOP together with DSL, after 1080 orbits around the Moon, which takes about 103 days, we can successfully detect the CD/EoR 21 cm signal.

Galaxy groups contain the majority of bound mass with a significant portion of baryons due to the combination of halo mass and abundance (Cui 2024). Hence they serve as a crucial missing piece in the puzzle of galaxy formation and the evolution of large-scale structures in the Universe. In observations, mass-complete group catalogues are normally derived from galaxy redshift surveys detected through various three-dimensional group-finding algorithms. Confirming the reality of such groups, particularly in the X-rays, is critical for ensuring robust studies of galaxy evolution in these environments. Recent works have reported numerous optical groups that are X-ray undetected (see, e.g., Popesso et al. 2024), sparking debates regarding the reasons for the unexpectedly low hot gas fraction in galaxy groups. To address this issue, we utilise zoomed-in simulations of galaxy groups from the novel HYENAS project to explore the range of hot gas fractions within galaxy groups and investigate the intrinsic factors behind the observed variability in X-ray emission. We find that the halo formation time can play a critical role -- we see that groups in halos that formed earlier exhibit up to an order of magnitude brighter X-ray luminosities compared to those formed later. This suggests that undetected X-ray groups are preferentially late-formed halos and highlights the connection between gas fraction and halo formation time in galaxy groups. Accounting for these biases in galaxy group identification is essential for advancing our understanding of galaxy formation and achieving precision in cosmological studies.

The statistics of Einstein radii for a sample of strong lenses can provide valuable constraints on the underlying mass distribution. The correct interpretation, however, relies critically on the modelling of the selection of the sample, which has proven to be a limiting factor. This may change thanks to upcoming uniform high-resolution imaging surveys that cover a large fraction of the sky, because they can provide complete lens samples, with well understood selection criteria. To explore how the observed distribution of Einstein radii depends on the galaxy properties, we simulated a realistic complete sample of strong lenses, predicting a number density of lenses of about 2.5 deg$^{-2}$ for a \Euclid-like setup. Such data can break the degeneracy between the stellar initial mass function (IMF) and the inner slope of the density profile of dark matter, without having to rely on additional information from stellar dynamics. We found that a survey covering only 50 deg$^2$ can already provide tight constraints: assuming that the cosmology is known, the dark matter slope is recovered with an uncertainty of $3.5\%$, while the uncertainty in the ratio between the true stellar mass and that inferred from stellar population modelling was found to be $10\%$. These findings highlight the potential of this method when applied to samples of lenses with well-understood selection functions.

Quasi-periodic eruptions (QPEs) are repeated X-ray flares from galactic nuclei. Despite some diversity in the recurrence and amplitude of eruptions, their striking regularity has motivated theorists to associate QPEs with orbital systems. Among the known QPE sources, eRO-QPE2 has shown the most regular flare timing and luminosity since its discovery. We report here on its long-term evolution over $\sim3.3\,$yr from discovery and find that: i) the average QPE recurrence time per epoch has decreased over time, albeit not at a uniform rate; ii) the distinct alternation between consecutive long and short recurrence times found at discovery has not been significant since; iii) the spectral properties, namely flux and temperature of both eruptions and quiescence components, have remained remarkably consistent within uncertainties. We attempted to interpret these results as orbital period and eccentricity decay coupled with orbital and disk precession. However, since gaps between observations are too long, we are not able to distinguish between an evolution dominated by just a decreasing trend, or by large modulations (e.g. due to the precession frequencies at play). In the former case, the observed period decrease is roughly consistent with that of a star losing orbital energy due to hydrodynamic gas drag from disk collisions, although the related eccentricity decay is too fast and additional modulations have to contribute too. In the latter case, no conclusive remarks are possible on the orbital evolution and the nature of the orbiter due to the many effects at play. However, these two cases come with distinctive predictions for future X-ray data: in the former, we expect all future observations to show a shorter recurrence time than the latest epoch, while in the latter we expect some future observations to be found with a larger recurrence, hence an apparent temporary period increase.

We revisit the relation between the variance of three-dimensional (3D) density ($\sigma^{2}_{\rho}$) and that of the projected two-dimensional (2D) column density ($\sigma^{2}_{\Sigma}$) in turbulent media, which is of great importance in obtaining turbulence properties from observations. Earlier studies showed that $\sigma^{2}_{\Sigma / \Sigma_{0}}/\sigma^{2}_{\rho / \rho_{0}} = \mathcal{R}$, where $\Sigma/\Sigma_0$ and $\rho/\rho_0$ are 2D column and 3D volume densities normalized by their mean values, respectively. The factor $\mathcal{R}$ depends only on the density spectrum for isotropic turbulence in a cloud that has similar dimensions along and perpendicular to the line of sight. Our major findings in this paper are as follows. First, we show that the factor $\mathcal{R}$ can be expressed in terms of $N$, the number of independent eddies along the line of sight. To be specific, $\sigma^{2}_{\Sigma / \Sigma_{0}}/\sigma^{2}_{\rho/\rho_{0}}$ is proportional to $\sim 1/N$, due to the averaging effect arising from independent eddies along the line of sight. Second, we show that the factor $\mathcal{R}$ needs to be modified if the dimension of the cloud in the line-of-sight direction is different from that in the perpendicular direction. However, if we express $\sigma^{2}_{\Sigma / \Sigma_{0}}/\sigma^{2}_{\rho / \rho_{0}}$ in terms of $N$, the expression remains same even in the case the cloud has different dimensions along and perpendicular to the line of sight. Third, when we plot $N\sigma^{2}_{\Sigma / \Sigma_{0}}$ against $\sigma^{2}_{\rho / \rho_{0}}$, two quantities roughly lie on a single curve regardless of the sonic Mach number, which implies that we can directly obtain the latter from the former. We discuss observational implications of our findings.

Gravitational lensing by galaxy clusters has emerged as a powerful tool to probe the standard Cold Dark Matter (CDM) paradigm of structure formation in the Universe. Despite the remarkable explanatory power of CDM on large scales, tensions with observations on small scales have been reported. Recent studies find that the observational cross-section for Galaxy-Galaxy Strong Lensing (GGSL) in clusters exceeds the CDM prediction by more than an order of magnitude, and persists even after rigorous examination of various possible systematics. We investigate the impact of modifying the internal structure of cluster dark matter sub-halos on GGSL and report that altering the inner density profile, given by $r^{\gamma}$, to steeper slopes with $\gamma > 2.5$ can alleviate the GGSL discrepancy. Deviating from the $\gamma \sim 1.0$ cusps that CDM predicts, these steeper slopes could arise in models of self-interacting dark matter undergoing core collapse. Our results motivate additional study of sub-halo core collapse in dense cluster environments.

We explore the intriguing phenomenon of time non-locality in the evolution of dark matter and Large Scale Structure (LSS). Recently in\,\cite{Donath:2023sav}, it was shown that time non-locality emerges in bias tracer fluctuations, which are $SO(3)$ scalars in real space, at fifth order in the perturbation expansion in dark matter overdensity. We demonstrate that by breaking the symmetry down to $SO(2)$, which is the case whenever line-of-sight effects become important, such as for flux fluctuations in the Lyman $\alpha$ forest, the temporal non-locality appears at the third order in expansion. Additionally, within the framework of EFTofLSS, we demonstrate that time non-locality manifests in the effective stress tensor of dark matter, which is a second rank tensor under $SO(3)$ transformations, again at the third order in dark matter overdensity. Furthermore, we highlight the effectiveness of the standard $\Pi$ basis\,\cite{Mirbabayi:2014zca} in handling time non-local operators.

The Li-rich stars are a class of rare objects with A(Li) higher than that of other stars in the same evolutionary stage. Their origin is still debated and valuable routes are the Cameron-Fowler mechanism, mass transfer process in a binary system or engulfment of small bodies. Metal-poor ([Fe/H]<-1 dex) stars are only a small fraction of the entire population of Li-rich stars. We observed with MIKE at the Magellan Telescope the metal-poor ([Fe/H]=-3.95+-0.11 dex) giant star HE0057-5959, deriving A(Li)=+2.09+-0.07 dex. Such A(Li) is higher by about 1 dex than that of other stars in the same evolutionary stage. A previous analysis of the same target suggested that its high A(Li) reflects a still ongoing First Dredge-Up process. We revise the nature of HE0057-5959 by comparing its stellar parameters and A(Li) with appropriate stellar evolution models. This comparison rules out that HE0057-5959 is caught during its First Dredge-Up, being this latter already ended according to the parameters of this star. Its A(Li), remarkably higher than the typical lithium plateau drawn by similar giant stars, demonstrates that HE0057-5959 joins the class of the rare metal-poor Li-rich stars. HE0057-5959 is the most metal-poor Li-rich star discovered so far. We consider different scenarios to explain this star. No internal mixing able to activate the Cameron-Fowler mechanism is known for metal-poor stars at this evolutionary stage. Also the engulfment of planets is disfavored because such metal-poor stars should not host planets. Finally, HE0057-5959 is one of the most Na-rich among the Li-rich stars and we found that a strong excess of Na is common to the three Li-rich stars with [Fe/H]<-3 dex. This finding could support the scenario of mass transfer from a massive companion star in a binary system, even if we found no evidence of radial velocity variations.

Quadruply-imaged strongly lensed quasars (quads) are routinely used for measurements of the expansion rate of the Universe with time delays. It has recently been suggested that any quad lens is subject to a Malmquist-like bias that causes the inferred area enclosed within the tangential caustic to be systematically underestimated, and that such a bias might translate into a corresponding bias on the expansion parameter. In this work we extended that analysis by also considering the effect of Eddington bias. We found that the sign and amplitude of the combined bias depend on the functional form of the caustic area distribution of the lens population and on the noise associated with the caustic area estimation process. Based on simulations, we estimated that the corresponding impact on $H_0$ is on the order of a percent or smaller. If the likelihood of the lensing data is known, then the bias can be accounted for when modelling the lens population. However, ignoring the criteria used to select a quad might lead to a bias at the lens modelling stage that causes the inferred caustic area to be overestimated. Such a bias disappears for lens models that are well constrained by the data.

Contribution of resolved and unresolved extragalactic point sources to the low-frequency sky spectrum is a potentially non-negligible part of the astrophysical foregrounds for cosmic dawn 21-cm experiments. The clustering of such point sources on the sky, combined with the frequency-dependence of the antenna beam, can also make this contribution chromatic. By combining low-frequency measurements of the luminosity function and the angular correlation function of extragalactic point sources, we develop a model for the contribution of these sources to the low-frequency sky spectrum. Using this model, we find that the contribution of sources with flux density $>10^{-6}\,$Jy to the sky-averaged spectrum is smooth and of the order of a few kelvins at 50 - $200\,$MHz. We combine this model with measurements of the galactic foreground spectrum and convolve the result with the beam of the conical log-spiral antenna planned as part of the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) project. We find that the contribution of point sources to resultant spectrum is $\sim0.4\%$ of the total foregrounds, but still larger by at least an order of magnitude than the standard predictions for the cosmological 21-cm signal. As a result, not accounting for the point-source contribution leads to a systematic bias in 21-cm signal recovery. We show, however, that in the REACH case, this reconstruction bias can be removed by modelling the point-source contribution as a power law with a running spectral index. We make our code publicly available as a Python package labelled epspy.

Filaments connecting halos are a long-standing prediction of cold dark matter theories. We present a novel detection of the cosmic web emission connecting two massive quasar-host galaxies at cosmic noon in the MUSE Ultra Deep Field (MUDF) using unprecedentedly deep observations that unlock a high-definition view of the filament morphology, a measure of the transition radius between the intergalactic and circumgalactic medium, and the characterization of the surface brightness profiles along the filament and in the transverse direction. Through systematic comparisons with simulations, we validate the filaments' typical density predicted in the current cold dark matter model. Our analysis of the MUDF field, an excellent laboratory for quantitatively studying filaments in emission, opens a new avenue to understanding the cosmic web that, being a fundamental prediction of cosmology, bears key information on the essence of dark matter.

One of the most dominant energy budget of the Universe is Dark Energy, which remains enigmatic since the claim of its existence from the observation of late-time cosmic acceleration. We propose a new way of inferring this by measuring the aging of the Universe using only gravitational wave (GW) signals from coalescing binary compact objects of any masses. We show that the aging of the Universe will lead to a change in the observed chirp mass of the GW sources inferred from different stages of its coalesces by monitoring a coherent source from two far-apart GW frequencies for a few years. We show that coordinated GW detectors which can reach a relative uncertainty on chirp mass measurement of about $10^{-9}$, can measure a tiny departure of about $2\%$ from the dark energy equation of state parameter $w_0=-1$ and its variation with cosmic time by using stellar origin binary black holes up to redshift $z=5$ in 10 years of observation time without using any external calibrator. If the next generation of GW detectors can achieve this precision, then it can open a new window to discover the fundamental nature of dark energy.

Presently, primordial black holes (PBHs) in the asteroid-mass window from $10^{-16}$ M$_\odot$ to $10^{-10}$ M$_\odot$ are a popular dark matter candidate. If they exist, some stars would capture them upon formation, and they would slowly accrete the star over gigayears. Such Hawking stars -- stars with a central PBH -- provide a novel channel for the formation of both sub-Chandrasekhar mass black holes as well as red straggler stars. Here we report on stellar evolution models that extend our previous work to Hawking stars with masses between 0.5 and 1.4 M$_\odot$. We explore three accretion schemes, and find that a wide range of PBHs in the asteroid-mass window can robustly accrete stars as small as 1 M$_{\odot}$ within the age of the Universe. This mechanism of producing sub-solar mass black holes is highly dependent on the assumed accretion physics and stellar metallicity. Lower-metallicity stars are generally accreted more rapidly, suggesting that it may be more likely for sub-Chandrasekhar mass Hawking stars formed in the early universe, such as those in ultra-faint dwarf (UFD) galaxies, to transmute their star into a sub-Chandrasekhar mass black hole within a Hubble time. We present a stellar population synthesis of a Draco II-like UFD galaxy containing Hawking stars and show that the number of red stragglers they produce can qualitatively match the observed population.

Employing the publicly available CosmoLattice code, we conduct numerical simulations of a domain wall network and the resulting gravitational waves (GWs) in a radiation-dominated Universe in the $Z_2$-symmetric scalar field model. In particular, the domain wall evolution is investigated in detail both before and after reaching the scaling regime, using the combination of numerical and theoretical methods. We demonstrate that the total area of closed walls is negligible compared to that of a single long wall stretching throughout the simulation box. Therefore, the closed walls are unlikely to have a significant impact on the overall network evolution. This is in contrast with the case of cosmic strings, where formation of loops is crucial for maintaining the system in the scaling regime. To obtain the GW spectrum, we develop a technique that separates physical effects from numerical artefacts arising due to finite box size and non-zero lattice spacing. Our results on the GW spectrum agree well with Refs. [29, 30], which use different codes. Notably, we observe a peak at the Hubble scale, an exponential falloff at scales shorter than the wall width, and a plateau/bump at intermediate scales. We also study sensitivity of obtained results on the choice of initial conditions. We find that different types of initial conditions lead to qualitatively similar domain wall evolution in the scaling regime, but with important variations translating into different intensities of GWs.

Vision-Language multimodal Models (VLMs) offer the possibility for zero-shot classification in astronomy: i.e. classification via natural language prompts, with no training. We investigate two models, GPT-4o and LLaVA-NeXT, for zero-shot classification of low-surface brightness galaxies and artifacts, as well as morphological classification of galaxies. We show that with natural language prompts these models achieved significant accuracy (above 80 percent typically) without additional training/fine tuning. We discuss areas that require improvement, especially for LLaVA-NeXT, which is an open source model. Our findings aim to motivate the astronomical community to consider VLMs as a powerful tool for both research and pedagogy, with the prospect that future custom-built or fine-tuned models could perform better.

We present the results of the XMM-Newton and NuSTAR observations taken as part of the ongoing, intensive multi-wavelength monitoring program of the Seyfert 1 galaxy Mrk 817 by the AGN Space Telescope and Optical Reverberation Mapping 2 (AGN STORM 2) Project. The campaign revealed an unexpected and transient obscuring outflow, never before seen in this source. Of our four XMM-Newton/NuSTAR epochs, one fortuitously taken during a bright X-ray state has strong narrow absorption lines in the high-resolution grating spectra. From these absorption features, we determine that the obscurer is in fact a multi-phase ionized wind with an outflow velocity of $\sim$5200 km s$^{-1}$, and for the first time find evidence for a lower ionization component with the same velocity observed in absorption features in the contemporaneous HST spectra. This indicates that the UV absorption troughs may be due to dense clumps embedded in diffuse, higher ionization gas responsible for the X-ray absorption lines of the same velocity. We observe variability in the shape of the absorption lines on timescales of hours, placing the variable component at roughly 1000 $R_g$ if attributed to transverse motion along the line of sight. This estimate aligns with independent UV measurements of the distance to the obscurer suggesting an accretion disk wind at the inner broad line region. We estimate that it takes roughly 200 days for the outflow to travel from the disk to our line of sight, consistent with the timescale of the outflow's column density variations throughout the campaign.

The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integral field spectroscopic observations with JWST/MIRI MRS. The H$_2$ line ratios, characteristic of a power-law temperature distribution, indicate that the gas is warmest where it enters the ionisation cone through disk rotation, suggestive of shock excitation. In the nucleus, where the data can be combined with an additional seven ro-vibrational H$_2$ transitions, we find that moderate velocity (30 km s$^{-1}$) shocks in dense ($10^5$ cm$^{-3}$) gas, irradiated by an external UV field ($G_0 = 10^3$), do provide a good match to the full set. The warm molecular gas in the ionisation cone that is traced by the H$_2$ rotational lines has been heated to temperatures $>200$ K. Outside of the ionisation cone the molecular gas kinematics are undisturbed. However, within the ionisation cone, the kinematics are substantially perturbed, indicative of a radial flow, but one that is quantitatively different from the ionised lines. We argue that this outflow is in the plane of the disk, implying a short 50 pc acceleration zone up to speeds of about 400 km s$^{-1}$ followed by an extended deceleration over $\sim$700 pc where it terminates. The deceleration is due to both the radially increasing galaxy mass, and mass-loading as ambient gas in the disk is swept up.

During Big Bang nucleosynthesis (BBN) in the first 15 minutes of the Universe some $^7$Li was created along with isotopes of H and He. The determination of that primordial value of Li can help constrain the conditions at that time. The oldest stars with known ages can be found in globular clusters which have well-determined ages through stellar evolution models. High-resolution spectra have been obtained with the Keck I telescope and HIRES of Li in several unevolved stars in the clusters M13 and M71 with V magnitudes of 17.6 -- 17.9. Abundances of Li have been determined with spectrum synthesis techniques and show a range of a factor of 4. We attribute that spread to differences in initial angular momentum resulting in different amounts of spin-down, related mixing, and destruction of Li. Our results are compared with similar results for main-sequences and turn-off stars in other globular clusters. The range in age of these clusters is 11.2 to 14.2 Gyr for an age span of 3 Gyr. These clusters range in [Fe/H] from $-$0.75 to $-$2.24 corresponding to a factor of 30 in metallicity. The maximum in the Li abundance for these unevolved stars in all eight clusters is the same corresponding to Li/H = 3.16 x 10$^{-10}$ while the predicted Li abundance, based on the deuterium abundance and the BBN predictions, is 5.24 x 10$^{-10}$.

The number of planetary satellites around solid objects in the inner Solar System is small either because they are difficult or unlikely to form, or that they do not survive for astronomical timescales. Here we conduct a pilot study on the possibility of satellite capture from the process of collision-less binary-exchange and show that massive satellites in the range 0.01 - 0.1 Earth masses can be captured by earth-sized terrestrial planets in a way already demonstrated for larger planets both in the Solar System and possibly beyond. In this process, one of the binary objects is ejected, leaving the other object as a satellite in orbit around the planet. We specifically consider satellite capture by an earth in an assortment of hypothetical encounters with large terrestrial binaries at 1 AU around the Sun. In addition, we examine the tidal evolution of captured objects and show that orbit circularization and long-term stability are possible for cases resembling the Earth-Moon system.

Context. Pulsars have a very stable rotation. However, sudden increases in their rotation frequency known as glitches, perturb their evolution. While large glitches are commonly detected, small glitches are harder to detect because of the lack of daily-cadence observations over long periods of time. Aims. We aim to explore the timing behaviour of young pulsars at daily timescales looking for small glitches and other irregularities. This will further our comprehension of the distribution of glitch sizes, which has also consequences for the theoretical modeling of the glitch mechanism. Methods. We observed six pulsars with up to daily cadence during 5 years with the antennas of the Argentine Institute of Radio Astronomy (IAR). We used standard pulsar timing tools to characterise the rotation of pulsars and developed an algorithm to look for small timing events in the data and calculate the changes in $\nu$ and $\dot\nu$ at those epochs. Results. We found that the rotation of pulsars in this dataset is affected by small step changes in $\nu$ and $\dot\nu$. We found three glitches that had not been reported before: two in PSR J1048-5832 with relative sizes $\Delta\nu / \nu= 9.1(4) \times 10^{-10}$ and $\Delta\nu / \nu = 4.5(1) \times 10^{-10}$, and one in the Vela pulsar with a size $\Delta\nu / \nu = 2.0(2) \times 10^{-10}$. We also report new decay terms on the 2021 Vela giant glitch, and on the 2022 giant glitches in PSR J0742-2822 and PSR J1740-3015 respectively. Besides, we found that the red noise contribution significantly diminished in PSR J0742-2822 after its giant glitch in 2022. Conclusions. Our results highlight the importance of high-cadence monitoring with an exhaustive analysis of the residuals to better characterize the distribution of glitch sizes and to deepen our understanding of the mechanisms behind glitches, red noise and timing irregularities.

Quasi-periodic eruptions (QPEs) are high-amplitude, soft X-ray bursts recurring every few hours, associated with supermassive black holes. Many interpretations for QPEs were proposed since their recent discovery in 2019, including extreme mass ratio inspirals and accretion disk instabilities. But, as of today, their nature still remains debated. We perform the first high-resolution X-ray spectral study of a QPE source using the RGS gratings onboard XMM-Newton, leveraging nearly 2 Ms of exposure on GSN 069, the first discovered source of this class. We resolve several absorption and emission lines including a strong line pair near the N VII rest-frame energy, resembling the P-Cygni profile. We apply photoionization spectral models and identify the absorption lines as an outflow blueshifted by $1700-2900$ km/s, with a column density of about $10^{22}$ cm$^{-2}$ and an ionization parameter $\log (\xi$/erg cm s$^{-1})$ of $3.9-4.6$. The emission lines are instead redshifted by up to 2900 km/s, and likely originate from the same outflow that imprints the absorption features, and covers the full $4\pi$ sky from the point of view of GSN 069. The column density and ionization are comparable to the outflows detected in some tidal disruption events, but this outflow is significantly faster and has a strong emission component. The outflow is more highly ionized when the system is in the phase during which QPEs are present, and from the limits we derive on its location, we conclude that the outflow is connected to the recent complex, transient activity of GSN 069 which began around 2010.

The dark matter halo properties, for example, mass, spin and concentration play a significant role in the formation and evolution of bars in disk galaxies. This study highlights the importance of a new parameter: the dark matter halo angular momentum distribution in the central region of disk. We experiment with N-body galaxy models having a disk and dark matter similar to Milky Way-type galaxies. In these models, we vary the discontinuity of the angular momentum distribution of the halo (the total spin is the same for all models). Our N-body experiments suggest that bar forms in all models after a few Gyr of disk evolution. However, in the secular evolution of the bar, as we evolve these models until 9.78 Gyr, the bar gains its strength in the model with the most continuous halo angular momentum distribution, and the bar loses strength for the most discontinuous halo angular momentum distribution. The secular evolution of the bar suggests that box/peanut/x-shaped bulges similar to those found in the Milky Way disk should be more pronounced in halos with continuous halo angular momentum distributions. This study demonstrates the importance of the initial condition setup of galaxy systems, namely the discontinuity in the dark matter halo angular momentum distribution for a given density distribution, on the bar secular evolution in the disk galaxy simulations. Further, this study helps reconcile the conflicting results of bar secular evolution in a high-spinning halo of the recent literature.

The extraplanar diffuse ionised gas is a key component for understanding the feedback processes that connect galactic discs and their halos. In this paper, we present the second study of the BETIS project, which aims to explore the ionisation mechanisms of the eDIG. We use a sample of eight edge-on galaxies observed with MUSE and apply the methodology developed in the first paper of the BETIS project. We found that the vertical and radial profiles of the [NII]/Ha, [SII]/Ha, [OIII]/Hb, and [OI]/Ha ratios depict a complex ionisation structure within galactic halos, influenced by the spatial distribution of HII regions across the galactic plane as observed from our line of sigh, with photon leakage from OB associations constituting the main ionisation source. Our analysis excludes low-mass, hot, and evolved stars as viable candidates for secondary ionisation sources to elucidate the unusual behaviour of the line ratios at greater distances from the galactic midplane. In contrast, we ascertain that shocks induced in the interstellar medium by star formation related feedback mechanisms represent a promising secondary ionisation source of the eDIG. We present a suite of models integrating ionisation mechanisms arising from fast shocks and photoionisation associated with star formation. When applied to the classical BPT diagrams, these models reveal that the ionisation budget of the eDIG ranges from 20% to 50% across our sample, with local variations of up to 20% within individual galaxy halos. This correlates with the presence of filaments and other structural components observed within galaxy halos. The presence of shocks is additionally supported by the observation of high-density, high [OI]/Ha ratios, characteristic of shock-compressed ionised gas, likely induced by feedback from regions of intense SF within the disk. These results are consistent across all galaxies analysed in this sample.

The epoch of reionisation progressed through the emission of ionising photons from galaxies to their local intergalactic medium. In this work, we characterise the dwarf star-forming galaxies as candidates for the source of ionising photons that drove EoR. We investigate the ionising properties and star formation histories of star-forming dwarf galaxies at the last stages of EoR at $4.8<\rm{z}<7$ using observations from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The magnification due to gravitational lensing allows us to probe large dynamic ranges in stellar mass ($2\times 10^{6}\leq\rm{M}_*/\rm{M}_\odot\leq5\times 10^{9}$) and UV magnitudes ($-22.68\leq$M$_{UV}\leq=-15.95$).We find a median UV slope \buv of $-2. 56\pm0.23$ and the production efficiency of ionising photons $\log$ \xiion $=25.39\pm0.6$ for the full sample ($4.8<\rm{z}<7$) with a median stellar mass of $6.3\pm0.5\times10^{7} \rm{M}_\odot$. We find both \buv and \xiion are marginally correlated with the stellar mass of the galaxy, indicating a possible greater contribution of dwarf galaxies to the reionisation of the Universe. We find that on average, galaxies in our sample are experiencing a recent rise/burst of star formation which translates to a higher scatter in \xiion and a large scatter in H$\alpha$ equivalent widths. Finally, we investigate the trends of H$\alpha$ and [OIII]+H$\beta$ EWs with UV magnitude and find M$_{UV}$ is correlated between H$\alpha$ but not with [OIII]+H$\beta$ EWs indicating low metallicities and recent burst in the UV faint galaxies.

Hard X-ray-emitting ($\delta$-type) symbiotic binaries, which exhibit a strong hard X-ray excess, have posed a challenge to our understanding of accretion physics in degenerate dwarfs. RT Cru, which is a member of $\delta$-type symbiotics, shows stochastic X-ray variability. Timing analyses of X-ray observations from XMM-Newton and NuSTAR, which we consider here, indicate hourly fluctuations, in addition to a spectral transition from 2007 to a harder state in 2012 seen with Suzaku observations. To trace the nature of X-ray variability, we analyze the multi-mission X-ray data using principal component analysis (PCA), which determines the spectral components mostly contributing to the flickering behavior and the hardness transition. The Chandra HRC-S/LETG and XMM-Newton EPIC-pn data provide the primary PCA components, which may contain some variable emission features, especially in the soft excess. Additionally, the absorbing column (first order with 50%), along with the source continuum (20%), and a third component (9%) - which likely accounts for thermal emission in the soft band - are the three principal components found in the Suzaku XIS1 observations. The PCA components of the NuSTAR data also correspond to the continuum and possibly emission features. Our findings suggest that the spectral hardness transition between the two Suzaku observations is mainly due to changes in the absorbing material and X-ray continuum, while some changes in the thermal plasma emission may result in flickering-type variations.

Millisecond pulsars, like pulsars, have led to major advances in many areas of astronomy and physics. The discovery in 2023 of a cosmological gravitational wave (GW) stochastic background using millisecond pulsar timing arrays has focused attention on the importance of millisecond pulsars, to both multi-messenger astronomy and cosmology, and for identifying the origin of the GW stochastic background, which is hypothesized to be due to supermassive black hole binaries (SMBHBs). Unlike pulsars, however, for which the details of the discovery are well-known, those of millisecond pulsars are not well known. In particular, the details of the first crucial step in the discovery of millisecond pulsars, namely the discovery of interplanetary scintillation (IPS) in the radio source 4C 21.53, are known only to the author. This article presents a first-hand account of this crucial first step, which resulted ultimately in the discovery of millisecond pulses from this object. A brief description of interplanetary and interstellar scintillation and scattering is given in the Appendix.

RF-ICE is a signal processing platform for the readout of large arrays of superconducting resonators. Designed for flexibility, the system's low digital latency and ability to independently and dynamically set the frequency and amplitude of probe tones in real time has enabled previously-inaccessible views of resonator behaviour, and opened the door to novel resonator control schemes. We introduce a multi-frequency imaging technique, developed with RF-ICE, which allows simultaneous observation of the entire resonance bandwidth. We demonstrate the use of this technique in the examination of the response of superconducting resonators to variations in applied readout current and thermal loading. We observe that, used in conjunction with a conventional frequency sweep at sufficiently large amplitude to induce resonance bifurcation, the multi-frequency imaging technique reveals a resonator response which is not captured by the frequency sweep measurement alone. We demonstrate that equivalent resonant frequency shifts can be achieved using either thermal, optical, or readout loading, and use this equivalence to counteract a change in thermal loading by digitally modulating the readout current through a resonator. We develop and implement a proof-of-concept closed-loop negative electro-quasiparticle feedback algorithm which first sets and then maintains the resonant frequency of a lumped element kinetic inductance detector while the loading on it is varied. Although this simple implementation is not yet suitable to deploy at scale, it demonstrates the utility of this feedback technique to improve linearity while addressing amplifier distortion, resonator response non-uniformity, and crosstalk. It can be applied to kinetic inductors in non-bolometric operation, and sets the stage for future developments.

Recent works have constrained the binary fraction of evolved populations of massive stars in local galaxies such as red supergiants and Wolf-Rayet stars, but the binary fraction of yellow supergiants (YSGs) in the Hertzsprung Gap remains unconstrained. Binary evolution theory predicts that the Hertzsprung Gap is home to multiple populations of binary systems with varied evolutionary histories. In this paper, we develop a method to distinguish single YSGs from YSG plus O- or B-type main sequence binaries using optical and ultraviolet photometry, and then apply this method to identify candidate YSG binaries in the Magellanic Clouds. After constructing a set of combined stellar atmosphere models, we find that optical photometry is, given typical measurement and reddening uncertainties, sufficient to discern single YSGs from YSG+OB binaries if the OB-star is at least $\sim5$M$_{\odot}$ for T$_{\mathrm{eff,YSG}}\sim$ 4000 K, but requires a $\sim$20M$_{\odot}$ OB star for YSGs up to T$_{\mathrm{eff,YSG}}\sim$ 9000 K. For these hotter YSG temperatures, ultraviolet photometry allows binaries with OB companions as small as $\sim$7M$_{\odot}$ to be identified. We use color-color spaces developed from these models to search for evidence of excess blue or ultraviolet light in a set of $\sim$1000 YSG candidates in the Magellanic Clouds. We identify hundreds of candidate YSG binary systems and report a preliminary fraction of YSGs that show a blue/UV color excess of 20-60\%. Spectroscopic follow-up is now required to confirm the true nature of this population.

The Submillimeter Array (SMA) is an array of 8 antennas operating at millimeter and submillimeter wavelengths on Maunakea, Hawaii, operated by the Smithsonian Astrophysical Observatory and Academia Sinica Institute of Astronomy and Astrophysics, Taiwan. Over the past several years, we have been preparing a major upgrade to the SMA that will replace the aging original receiver cryostats and receiver cartridges with all new cryostats and new 230 and 345 GHz receiver designs. This wideband upgrade (wSMA) will also include significantly increased instantaneous bandwidth, improved sensitivity, and greater capabilities for dual frequency observations. In this paper, we will describe the wSMA receiver upgrade and status, as well as the future upgrades that will be enabled by the deployment of the wSMA receivers.

Mature neutron stars are thought to be sufficiently cold that nuclei in the outer layers freeze, solidifying a crust. Crustal elasticity allows the star to support a set of seismic modes, such as torsional oscillations. These axial-parity modes can couple to the polar sector in a number of ways, for example via rotation or a magnetic field. Even in a static, spherically-symmetric star however these modes can couple at non-linear order. In this study, such couplings in the crust are examined for the first time: we derive the axisymmetric perturbation equations for second-order axial eigenfunctions, which are sourced by axial-polar couplings at first-order, and solve the resulting equations in the time domain. Through our studies, we find that the second-order spectrum contains additional oscillation modes, not predicted by the linear analysis with either axial or polar perturbations, which can be excited to relatively large amplitudes.

We study spectra produced by weakly accreting black hole (BH) systems using the semi-analytic advection dominated accretion flow (ADAF) model and the general-relativistic magentohydrodynamic (GRMHD) simulation. We find significant differences between these two approaches related to a wider spread of the flow parameters as well as a much steeper radial distribution of the magnetic field in the latter. We apply these spectral models to the broad-band spectral energy distribution (SED) of the nucleus of M87 galaxy. The standard (in particular, one-dimensional) formulation of the ADAF model does not allow to explain it; previous claims that this model reproduces the observed SED suffer from an inaccurate treatment of the Compton process. The spectra based on GRMHD simulation are in a much better agreement with the observed data. In our GRMHD model, in which we assumed the BH spin $a=0.9$, bulk of radiation observed between the millimeter and the X-ray range is produced in the disk area between 2 and 4 gravitational radii from the BH. In this solution, the synchrotron component easily reproduces the spectral data between the millimeter and the UV range, and the Compton component does not violate the X-rays constraints, for accretion rates not exceeding 0.01 Msun/year and a relatively strong magnetic field, with the plasma $\beta \sim 1$ in the region where radiation is produced. However, the Compton component cannot explain the observed X-ray spectrum. Instead, the X-ray spectrum can be reproduced by a high energy tail of the synchrotron spectrum if electrons have a hybrid energy distribution with a $\sim 5$ per cent nonthermal component.

The discovery that blazars dominate the extra-galactic {\gamma}-ray sky is a triumph in the Fermi era. However, the exact location of {\gamma}-ray emission region still remains in debate. Low-synchrotron-peaked blazars (LSPs) are estimated to produce high-energy radiation through the external Compton process, thus their emission regions are closely related to the external photon fields. We employed the seed factor approach proposed by Georganopoulos et al. It directly matches the observed seed factor of each LSP with the characteristic seed factors of external photon fields to locate the {\gamma}-ray emission region. A sample of 1138 LSPs with peak frequencies and peak luminosities was adopted to plot a histogram distribution of observed seed factors. We also collected some spectral energy distributions (SEDs) of historical flare states to investigate the variation of {\gamma}-ray emission region. Those SEDs were fitted by both quadratic and cubic functions using the Markov-chain Monte Carlo method. Furthermore, we derived some physical parameters of blazars and compared them with the constraint of internal {\gamma}{\gamma}-absorption. We find that dusty torus dominates the soft photon fields of LSPs and most {\gamma}-ray emission regions of LSPs are located at 1-10 pc. The soft photon fields could also transition from dusty torus to broad line region and cosmic microwave background in different flare states. Our results suggest that the cubic function is better than the quadratic function to fit the SEDs.

With the number of supernovae observed expected to drastically increase thanks to large-scale surveys like the Dark Energy Spectroscopic Instrument (DESI), it is necessary that the tools we use to classify these objects keep up with this increase. We previously created Supernova Tagging and Classification (STag) to address this problem by employing machine learning techniques alongside logistic regression in order to assign 'tags' to spectra based on spectral features. STag II is a continuation of this work, which now makes use of model supernova spectra combined with real DESI spectra in order to train STag to better deal with realistic data. We also make use of the rlap score as a trustworthiness cut, making for a more robust and accurate supernova classifier than before.

Both the anomalous magnetic braking of Ap/Bp stars and the surrounding circumbinary disk models can account for the formation of black hole (BH) low-mass X-ray binaries (LMXBs), while the simulated effective temperatures of the donor stars are significantly higher than the observed values. Therefore, the formation of BH LMXBs is not still completely understood. In this work, we diagnose whether the dynamical friction between dark matter and the companion stars can drive BH binaries to evolve toward the observed BH LMXBs and alleviate the effective temperature problem. Assuming that there exists a density spike of dark matter around BH, the dynamical friction can produce an efficient angular momentum loss, driving BH binaries with an intermediate-mass companion star to evolve into BH LMXBs for a spike index higher than $\gamma = 1.58$. Our detailed stellar evolution models show that the calculated effective temperatures can match the observed value of most BH LMXBs for a spike index range of $\gamma = 1.7-2.1$. However, the simulated mass-transfer rates when $\gamma = 2.0$ and $2.1$ are too high to be consistent with the observed properties that BH LMXBs appears as soft X-ray transients. Therefore, the dynamical friction of dark matter can only alleviate the effective temperature problem of those BH LMXBs with a relatively short orbital period.

The calculation of the emerging radiation from a model atmosphere requires knowledge of the emissivity and absorption coefficients, which are proportional to the atomic level population densities of the levels involved in each transition. Due to the intricate interdependency of the radiation field and the physical state of the atoms, iterative methods are required in order to calculate the atomic level population densities. A variety of different methods have been proposed to solve this problem, which is known as the Non-Local Thermodynamical Equilibrium (NLTE) problem. In this study we have developed a Jacobian-Free Newton-Krylov method (JFNK) to solve multi-level NLTE radiative transfer problems. Using the Rybicki & Hummer (1992) method as a reference (Rybicki, G. B. & Hummer, D. G. 1992, A&A, 262, 209), our results show that our JFNK solver can achieve up to a factor two speed up when using local approximate operators / preconditioners, while also achieving a lower residual error in the statistical equilibrium equations. Another advantage of this method is that the addition of charge conservation and partial redistribution effects should be straight forward. Our method can help accelerating the calculation of the emerging spectra from numerical models and also the reconstruction of chromospheric datasets through NLTE inversions.

We present a versatile transit simulator aimed at generating light curves for arbitrarily shaped objects transiting stars. Utilizing a Monte Carlo algorithm, it accurately models the stellar flux blocked by these objects, producing precise light curves. The simulator adeptly handles realistic background stars, integrating effects such as tidal distortions and limb darkening, alongside the rotational dynamics of transiting objects of arbitrary geometries. We showcase its wide-ranging utility through successful simulations of light curves for single and multi-planet systems, tidally distorted planets, eclipsing binaries and exocomets. Additionally, our simulator can simulate light curves for hypothetical alien megastructures of any conceivable shape, providing avenues to identify interesting candidates for follow-up studies. We demonstrate applications of the simulator in modeling a Dyson Swarm in construction, Dyson rings and Dyson disks, discussing how tidally locked Dyson disks can be distinguished from planetary light curves.

BL and UV Ceti are a nearby (2.7 pc) binary system with similar masses, spectral types, and rapid rotation rates, but very different magnetic activity. UV Ceti's much stronger large-scale magnetic field may cause this difference, highlighting key unanswered questions about dynamo processes in fully convective objects. Here we present multi-epoch characterization of the radio spectrum of UV Ceti spanning 1-105 GHz, exhibiting flared emission similar to coronal activity, auroral-like emission analogous to planetary magnetospheres, and slowly-varying persistent emission. Radio observations are a powerful means to probe the role that the large-scale magnetic field of UV Ceti has in non-thermal particle acceleration, because radio-frequency phenomena result from both the activity of small-scale field features as well as large-scale auroral current systems. We find temporal variability at all bands observed, and a hint of rotational modulation in the degree of circular polarization up to 40 GHz. The persistent component of the emission is fairly constant from 1-105 GHz, making optically thick emission or optically thin gyrosynchrotron from electrons with an isotropic pitch angle distribution unlikely. We discuss the possibility of emission mechanisms analogous to Jupiter's radiation belts.

This review summarizes popular unsupervised learning methods, and gives an overview of their past, current, and future uses in astronomy. Unsupervised learning aims to organise the information content of a dataset, in such a way that knowledge can be extracted. Traditionally this has been achieved through dimensionality reduction techniques that aid the ranking of a dataset, for example through principal component analysis or by using auto-encoders, or simpler visualisation of a high dimensional space, for example through the use of a self organising map. Other desirable properties of unsupervised learning include the identification of clusters, i.e. groups of similar objects, which has traditionally been achieved by the k-means algorithm and more recently through density-based clustering such as HDBSCAN. More recently, complex frameworks have emerged, that chain together dimensionality reduction and clustering methods. However, no dataset is fully unknown. Thus, nowadays a lot of research has been directed towards self-supervised and semi-supervised methods that stand to gain from both supervised and unsupervised learning.

We report two binary systems, LAMOST J035540+381550 (hereafter J035540) and LAMOST J035916+400732 (hereafter J035916), identified through the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution survey (LRS). Each of these two systems contains an M-type star orbiting with a invisible compact object candidate. Follow-up spectroscopic observations of Palomar 200-inch telescope (P200) enhance radial velocity measurements. We use radial velocities from LAMOST and P200, as well as light curves from Zwicky Transient Facility (ZTF) to constrain orbital parameters. The masses of the visible M-type stars are estimated by fitting the MIST isochrones and SEDs. The mass functions for the unseen companions are: $0.22\pm 0.01 M_{\odot}$ for J035540 and $0.16\pm 0.01 M_{\odot}$ for J035916. With the orbital and stellar parameters derived above and assuming the orbital inclination is 90 degree (edge-on), we find that the minimum masses of the invisible companions exceeds that of the visible stars. The single-lined feature and the dynamical evidence suggest the presence of compact objects. J035540's ZTF light curve, modeled with PHOEBE, yields a compact object mass of $0.70^{+0.12}_{-0.05} M_{\odot}$. For J035916, ellipsoidal modulation analysis constrains the light curve amplitude, yielding a compact object mass range of $0.57-0.90 M_{\odot}$. The mass estimates indicate that both are likely white dwarfs. These findings underscore the efficiency of optical time-domain surveys and dynamical methods in identifying faint, massive white dwarfs, along with other compact objects in binaries.

We present proper motion (PM) measurements for a sample of 23 massive star clusters within the Large Magellanic Cloud using multi-epoch data from the Hubble Space Telescope (HST). We combined archival data from the ACS/WFC and WFC3/UVIS instruments with observations from a dedicated HST programme, resulting in time baselines between 4.7 and 18.2 yr available for PM determinations. For bright well-measured stars, we achieved nominal PM precisions of 55 $\mu$as/yr down to 11 $\mu$as/yr . To demonstrate the potential and limitations of our PM data set, we analysed the cluster NGC 1850 and showcase a selection of different science applications. The precision of the PM measurements allows us to disentangle the kinematics of the various stellar populations that are present in the HST field. The cluster has a centre-of-mass motion that is different from the surrounding old field stars and also differs from the mean motion of a close-by group of very young stars. We determined the velocity dispersion of field stars to be 0.128 +/- 0.003 mas/yr (corresponding to 30.3 +/- 0.7 km/s). The velocity dispersion of the cluster inferred from the PM data set most probably overestimates the true value, suggesting that the precision of the measurements at this stage is not sufficient for a reliable analysis of the internal kinematics of extra-galactic star clusters. Finally, we exploit the PM-cleaned catalogue of likely cluster members to determine any radial segregation between fast and slowly-rotating stars, finding that the former are more centrally concentrated. With this paper, we also release the astro-photometric catalogues for each cluster.

We investigate the chemical evolution of complex organic molecules (COMs) in turbulent disks using gas-ice chemical reaction network simulations. We trace trajectories of dust particles considering advection, turbulent diffusion, gas drag, and vertical settling, for 10$^6$ yrs in a protoplanetary disk. Then, we solve a gas-ice chemical reaction network along the trajectories and obtain the temporal evolution of molecular abundances. We find that the COM abundances in particles can differ by more than two orders of magnitude even when the UV fluence (i.e., the time integral of UV flux) received by the particles are similar, suggesting that not only the UV fluence but also the time variation of the UV flux does matter for the evolution of COMs in disks. The impact of UV fluence on molecular abundances differs between oxygen-bearing and nitrogen-bearing COMs. While higher UV fluence results in oxygen being locked into CO$_2$, leading to reduced abundances of oxygen-bearing COMs such as CH$_3$OCH$_3$, mild UV exposure can promote their formation by supplying the precursor radicals. On the other hand, nitrogen is not locked up into specific molecules, allowing the formation of nitrogen-bearing COMs, particularly CH$_3$NH$_2$, even for the particle that receives the higher UV fluence. We also find that the final COM abundances are mostly determined by the inherited abundances from the protostellar core when the UV fluence received by dust particles is less than a critical value, while they are set by both the inherited abundances and the chemistry inside the disk at higher UV fluence.

The gravitational wave (GW) signal from binary black hole (BBH) mergers is a promising probe of Population III (Pop III) stars. To fully unleash the power of the GW probe, one important step is to understand the relative importance and features of different BBH evolution channels. We model two channels, isolated binary stellar evolution (IBSE) and nuclear star cluster-dynamical hardening (NSC-DH), in one theoretical framework based on the semi-analytical code A-SLOTH, under various assumptions on Pop III initial mass function (IMF), initial binary statistics and high-$z$ nuclear star clusters (NSCs). The NSC-DH channel contributes $\sim 8-95\%$ of Pop III BBH mergers across cosmic history, with higher contributions achieved by initially wider binary stars, more top-heavy IMFs, and more abundant high-$z$ NSCs. The dimensionless stochastic GW background (SGWB) produced by Pop III BBH mergers has peak values $\Omega^{\rm peak}_{\rm GW}\sim 10^{-11}-8\times 10^{-11}$ around observer-frame frequencies $\nu\sim 10-100\ \rm Hz$. The Pop III contribution can be a non-negligible ($\sim 2-32\%$) component in the total SGWB at $\nu\lesssim 10\ \rm Hz$. The estimated detection rates of Pop III BBH mergers by the Einstein Telescope are $\sim 6-230\ \rm yr^{-1}$ and $\sim 30-1230\ \rm yr^{-1}$ for the NSC-DH and IBSE channels, respectively. Pop III BBH mergers in NSCs are more massive than those from IBSE, so they dominate the Pop III SGWB below $20$ Hz in most cases. Besides, the detection rate of Pop III BBH mergers involving at least one intermediate-mass BH above $100\ \rm M_\odot$ by the Einstein Telescope is $\sim 0.5-200\ \rm yr^{-1}$ in NSCs but remains below $0.1\ \rm yr^{-1}$ for IBSE.

At optical/ultraviolet energies, blazars display an underlying thermal (unpolarized) contribution from the accretion disc, torus and line emitting regions, diluting the polarized emission from the jet-component. Optical polarimetry can be used to disentangle the thermal and non-thermal components, and place constraints on the particle populations and acceleration mechanisms responsible for the non-thermal emission. We present the results of a linear optical spectropolarimetric observing campaign of 18 blazars (6 BLLs and 12 FSRQs) undertaken with the Southern African Large Telescope between 2016 and 2022. This was done to observe these systems during flaring states, as well as long term monitoring of PKS1510-089, AP Lib and PKS 1034-293. The observations traced the frequency dependence of the degree and angle of polarization, as well as changes in the spectral line strengths. We investigated possible correlations between the polarization and other observed characteristics for the sources. While an indication of correlation was found between the frequency dependence and the average level of polarization for some sources, a correlation was not found for the population as a whole. These results highlight that continuous observations and in-depth modelling of polarization and its frequency dependence is required to obtain a more holistic view of TeV blazars.

The blue compact dwarf galaxy I Zw 18 is one of the most metal-poor ($Z \sim 3% Z_{\sun}$) star-forming galaxies in the local Universe. Its evolutionary status has sparked debate within the astronomical community. We aim to investigate the stellar populations of I Zw 18 in the near-IR using JWST/NIRCam's high spatial resolution and sensitivity. Additionally, we aim to derive the galaxy's spatially resolved star formation history (SFH) over the last 1 Gyr and provide constraints for older epochs. We used DOLPHOT to measure positions and fluxes of point sources in the F115W and F200W filters' images of I Zw 18. To derive I Zw 18's SFH, we applied the color-magnitude diagram (CMD) fitting technique SFERA 2.0, using two independent sets of stellar models. Our analysis reveals three main stellar populations: one younger than $\sim30$ Myr, mainly in the northwest star-forming (SF) region; an intermediate-age population ($\sim 100 - 800$ Myr) in the southeast SF region; and a red and faint population linked to the underlying halo, older than 1 Gyr and possibly as old as 13.8 Gyr. The main body of the galaxy shows a very low star formation rate (SFR) of $\sim 10^{-4} M_{\odot} \text{yr}^{-1}$ between 1 and 13.8 Gyr ago. In the last billion years, I Zw 18 shows increasing SF, with strong bursts around $\sim10$ and $\sim100$ Myr ago. Component C mirrors the main body's evolution but with lower SFRs. Our findings confirm that I Zw 18 contains stars of all ages, indicating it is not a young galaxy but has an old stellar halo, similar to other BCDs. The low SF activity over the past billion years supports the "slow cooking" dwarf scenario, explaining its low metal content. Currently, the galaxy is undergoing its strongest SF episode ($\sim 0.6 M_{\odot} \text{yr}^{-1}$) mainly in the northwest region, likely due to a recent gravitational interaction with Component C.

Supernova remnant (SNR) detection along the Galactic plane poses a number of challenges. The SNR G309.8+00.0 lies exactly on the Galactic plane, with its center coinciding with galactic latitude (b)=0 deg. In this paper we report the first detection of the SNR G309.8+00.0 in X-rays and $\gamma$ rays, using stacked data from the first four consecutive extended ROentgen Survey Imaging Telescope Array (eROSITA) -- on board the Russian-German Spektrum Roentgen Gamma (SRG) -- all-sky surveys (eRASS:4) and $\sim15.5$ yr of Pass 8 data recorded from Fermi-LAT, respectively. The SNR appears to have an elliptical shape of 0.43 x 0.32 deg in size in both radio synchrotron and X-ray data. The SNR's emission exhibits a shell-like morphology and good spatial correlation in both energy bands. The X-ray emission was solely detected in the 1-2 keV energy band (subject to strong absorption at soft X-rays) and the spectral analysis results of eRASS:4 data present a purely thermal SNR with a high absorption column density $3.1_{-0.5}^{+0.7}\cdot10^{22}~\mathrm{cm^{-2}}$ and a temperature of $0.34\pm0.1$ keV. In combination with optical extinction data, the absorption column density values derived from the remnant's spectral analysis support a remnant's distance greater than 6 kpc, rather than a 3.12 kpc distance as reported in the literature, and yield an age of $1-3.5\cdot10^5$ yr. Employing $\sim15.5$ yr of Fermi-LAT $\gamma$-ray data at and around the remnant's vicinity, we confirm the detection of the to-date unidentified 4FGL J1349.5-6206c source that can either be modeled as a single source or a conglomerate of multiple distinct source components and we argue that the SNR G309.8+00.0 likely represents at least a significant portion (if not all) of the emission from the 4FGL J1349.5-6206c $\gamma$-ray source.

The merger of compact binary stars produces short gamma-ray bursts (sGRBs), involving channels such as neutron star - neutron star (BNS) and neutron star - black hole (NS-BH). The association between sGRB 170817A and gravitational wave GW 170817 provides reliable evidence for the BNS channel. The spatial distribution and merger rate differ between BNS mergers and NS-BH mergers. Some speculations suggest that sGRBs with extended emission (EE) may represent another distinct population. We compared the offset distributions of these two types of samples and found that they follow the same distribution. Utilizing non-parametric methods, we investigated the origin of these burst types in terms of their formation rate. We examined the luminosity function and formation rate of sGRBs without assuming any specific model.The luminosity function can be described as $\psi(L_{0}) \propto L_{0}^{-0.09 \pm 0.01}$ for $L_{0} < L_0^b$ for standard sGRBs and $\psi(L_{0}) \propto L_{0}^{-0.09 \pm 0.01}$ for $L_{0} < L_0^b$ for sGRBs with EE. The formation rate is characterized as $\rho(z) \propto (1 + z)^{-4.21 \pm 0.22}$ for $z < 0.8$ and $\rho(z) \propto (1 + z)^{-0.22 \pm 0.47}$ for $0.8 < z < 3$ for standard sGRBs, while for sGRBs with EE, it is $\rho(z) \propto (1 + z)^{-4.30 \pm 0.13}$ for $z < 0.8$ and $\rho(z) \propto (1 + z)^{-0.33 \pm 0.66}$ for $0.8 < z < 3$. Based on these findings, we suggest that there is no significant difference in the progenitor stars of sGRBs with and without EE, considering the spatial offset and formation rate perspectives.

The imaging atmospheric Cherenkov technique provides potentially the highest angular resolution achievable in astronomy at energies above the X-ray waveband. High-resolution measurements provide the key to progress on many of the major questions in high-energy astrophysics, including the sites of particle acceleration to PeV energies. The potential of the next-generation CTA observatory in this regard can be realised with the help of improved algorithms for the reconstruction of the air-shower direction and energy. Hybrid methods combining likelihood-fitting techniques with neural networks represent a particularly promising approach and have recently been applied to the reconstruction of astrophysical neutrinos. Here, we present the FreePACT algorithm, a hybrid reconstruction method for IACTs. In this, making use of the neural ratio estimation technique from the field of likelihood-free inference, the analytical likelihood used in traditional image likelihood fitting is replaced by a neural network that approximates the charge probability density function for each pixel in the camera. The performance of this algorithm is demonstrated using simulations of the planned CTA southern array. For this setup, FreePACT provides significant performance improvements over analytical likelihood techniques, with improvements in angular and energy resolution of 25% or more over a wide energy range and an angular resolution as low as 40 arcseconds at energies above 50 TeV for observations at 20 degrees zenith angle. It also yields more accurate estimations of the uncertainties on the reconstructed parameters and speeds up the reconstruction compared to analytical likelihood techniques while showing the same stability with respect to changes in the observation conditions. Therefore, the FreePACT method is a promising upgrade over the current state-of-the-art likelihood event reconstruction techniques.

Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ component from the accreted stellar halos. The low surface brightness of stellar halos makes it inaccessible with traditional integrated light spectroscopy. In this work, using a novel technique, we study the kinematics of the stellar halo of the edge-on galaxy NGC 4945. We couple new deep Multi Unit Spectroscopic Explorer spectroscopic observations with existing Hubble Space Telescope imaging data to spectroscopically measure the line-of-sight (LOS) heliocentric velocity and velocity dispersion in two fields at a galactocentric distance of 12.2 kpc (outer disk field) and 34.6 kpc (stellar halo field) along NGC 4945 major axis, by stacking individual spectra of red giant branch and asymptotic giant branch stars. We obtain a LOS velocity and dispersion of 673+/-11 km/s and 73+/-14 km/s, respectively, for the outer disk field. This is consistent with the mean HI velocity of the disk at that distance. For the halo field we obtain a LOS velocity and dispersion of 519+/-12 km/s and 42+/-22 km/s. The halo fields' velocity measurement is within ~40 km/s from the systemic LOS velocity of NGC 4945, which is 563 km/s, suggesting that its stellar halo at 34.6 kpc along the major axis is counter-rotating and is of likely accretion origin. This provides the first ever kinematic measurement of the stellar halo of a Milky Way-mass galaxy outside the Local Group from its resolved stellar population, and establishes a powerful technique for measuring the velocity field of the stellar halos of nearby galaxies.

We implement a local model for a spherical collapsing/expanding gas cloud into the Athena++ magnetohydrodynamic code. This local model consists of a Cartesian periodic box with time-dependent geometry. We present a series of benchmark test problems, including non-linear solutions and linear perturbations of the local model, confirming the code's desired performance. During a spherical collapse, a horizontal shear flow is amplified, corresponding to angular momentum conservation of zonal flows in the global problem; wave speed and amplitude of sound waves increase in the local frame, due to the reduction in the characteristic length scale of the box, which can lead to an anisotropic effective sound speed in the local box. Our code conserves both mass and momentum to machine precision. This numerical implementation of the local model has potential applications to the study of local physics and hydrodynamic instabilities during protostellar collapse, providing a powerful framework for better understanding the earliest stages of star and planet formation.

Astronomical observatories have been identified as substantial contributors to the carbon footprint of astrophysical research. Being part of the collaboration that currently develops the Medium-Sized Telescope (MST) of the Cherenkov Telescope Array, a ground-based observatory for very-high-energy gamma rays that will comprise 64 telescopes deployed on two sites, we assessed the environmental impacts of one MST on the Northern site by means of a Life Cycle Assessment. We identified resource use and climate change as the most significant impacts, being driven by telescope manufacturing and energy consumption during operations. We estimate life cycle greenhouse gas emissions of 2,660 +/- 274 tCO2 equivalent for the telescope, 44% of which arise from construction, 1% from on-site assembly and commissioning, and 55% from operations over 30 years. Environmental impacts can be reduced by using renewable energies during construction and operations, use of less electronic components and metal casting, and use of recycled materials. We propose complementing project requirements with environmental budgets as an effective measure for impact management and reductions.

CO$_2$ is a major component of the icy mantles surrounding dust grains in planet and star formation regions. Understanding its photodesorption is crucial for explaining gas phase abundances in the coldest environments of the interstellar medium irradiated by vacuum-UV (VUV) photons. Photodesorption yields determined experimentally from CO$_2$ samples grown at low temperatures (T=15~K) have been found to be very sensitive to experimental methods and conditions. Several mechanisms have been suggested for explaining the desorption of CO$_2$, O$_2$ and CO from CO$_2$ ices. In the present study, the cross sections characterizing the dynamics of photodesorption as a function of photon fluence (determined from released molecules in the gas phase) and of ice composition modification (determined in situ in the solid phase) are compared for the first time for different photon flux conditions (from 7.3$\times 10^{12}$~photon/s/cm$^2$ to 2.2$\times 10^{14}$~photon/s/cm$^2$) using monochromatic synchrotron radiation in the VUV range (on the DESIRS beamline at SOLEIL). This approach reveals that CO and O$_2$ desorption are decorrelated from that of CO$_2$. CO and O$_2$ photodesorption yields depend on photon flux conditions and can be linked to surface chemistry. By contrast, the phodesorption yield of CO$_2$ is independent of the photon flux conditions and can be linked to bulk ice chemical modification, consistently with an indirect desorption induced by electronic transition (DIET) process.

We introduce an improved method for decomposing the emission of active galactic nuclei (AGN) and their host galaxies using templates from principal component analysis (PCA). This approach integrates prior information from PCA with a penalized pixel fitting mechanism which improves the precision and effectiveness of the decomposition process. Specifically, we have reduced the degeneracy and over-fitting in AGN-host decomposition, particularly for those with low signal-to-noise ratios (SNR), where traditional methods tend to fail. By applying our method to 76,565 SDSS Data Release 16 quasars with $z<0.8$, we achieve a success rate of $\approx$ 94%, thus establishing the largest host-decomposed spectral catalog of quasars to date. Our fitting results consider the impact of the host galaxy on the overestimation of the AGN luminosity and black hole mass ($M_{\rm BH}$). Furthermore, we obtained stellar velocity dispersion ($\sigma_*$) measurements for 4,137 quasars. The slope of the $M_{\rm BH}-\sigma_*$ relation in this subsample is generally consistent with previous quasar studies beyond the local universe. Our method provides a robust and efficient approach to disentangle the AGN and host galaxy components across a wide range of SNRs and redshifts.

For the first time, in March 2024, the transient Galactic black hole candidate Swift J151857.0-572147 experienced an outburst. Using publicly available archived {\it Insight}-HXMT data, we analyze the timing and spectral features of this source. Through model fitting of the power density spectrum, we were able to extract the properties of quasi-periodic oscillations, and based on those properties, we have determined that the QPOs are of type C. We also conclude that the shock instabilities in the transonic advective accretion processes surrounding black holes may be the source of the QPOs. This shock instability could produce variabilities of flux up to 48 keV, as we checked from the QPO energy dependence. High-frequency QPO is not observed during this period. In the broad energy band of $2-100$ keV, simultaneous data from the three on-board instruments of \textit{Insight}-HXMT were used to perform the spectral analysis. A combination of models, including broken power-law, multi-color disk-blackbody continuum, interstellar absorption, and reflection in both neutral and ionized medium were needed for spectral fitting to obtain the best fit. We discovered that at the beginning of the analysis period, the source was in an intermediate state and was transitioning toward the softer states based on the spectral features. It has a hydrogen column density of $(4.3-6.9) \times 10^{22}$ cm$^{-2}$.

SPHERE, operating at the VLT since 2014, is currently one of the high-contrast instruments with a higher performance. Its adaptive optics system, known as SAXO, will be upgraded to SAXO+, which features the addition of a second stage of adaptive optics. This stage will use a near-infrared pyramid wavefront sensor to record images of fainter exoplanets around redder stars. In this work, we compare the performance of SAXO and SAXO+. We look for the optimal values of the key system parameters of SAXO+ for various science cases and turbulence conditions. We performed numerical simulations using COMPASS, an end-to-end adaptive optics simulation tool. We simulated perfect coronagraph images of an on-axis point source, and we minimized the residual starlight intensity between 3 and 5 ${\lambda}/D$ as a performance criterion. The explored parameter space includes science cases, turbulence conditions, and key system parameters. In every science case and turbulence condition, SAXO+ reduces the residual starlight intensity inside the correction zone of the second stage by a factor of ten compared to SAXO. The optimal first stage gain is lower for SAXO+ than for SAXO alone. We quantified the gain in performance of SAXO+ when changing the second stage frequency from 2 kHz to 3 kHz, and we conclude that 2 kHz may be sufficient for most realistic conditions. We give the optimal first stage gain as well as the first and second stage frequencies for every seeing, coherence time, and science case. Finally, we find that a 2 ${\lambda_{\mathrm{WFS}}}/D$ pyramid modulation radius is a good trade-off between performance and robustness against varying turbulence conditions. This study shows that the future SAXO+ system will outperform the current SAXO system in all studied cases.

Context. Both the stellar activity and the accretion processes of young stellar objects can induce variations in their radial velocity (RV). This variation is often modulated on the stellar rotation period and may hide a RV signal from a planetary or even a stellar companion. Aims. The aim of this study is to detect the companion of HQ Tau, the existence of which is suspected based on our previous study of this object. We also aim to derive the orbital elements of the system. Methods. We used multi-variate Gaussian process regression on the RV and the bisector inverse slope of a six-month high-resolution spectroscopic follow-up observation of the system to model the stellar activity. This allowed us to extract the Keplerian RV modulation induced by the suspected companion. Results. Our analysis yields the detection of a $\sim$50 M$_{\rm jup}$ brown dwarf companion orbiting HQ Tau with a $\sim$126 day orbital period. Although this is consistent with the modulation seen on this dataset, it does not fit the measurements from our previous work three years earlier. In order to include these measurements in our analysis, we hypothesise the presence of a third component with orbital elements that are consistent with those of the secondary according to our previous analysis (M$_{\rm B}$ $\sim$48 M$_{\rm jup}$, P$_{\rm orb,B}$ $\sim$126 days), and a $\sim$465 M$_{\rm jup}$ tertiary with a $\sim$767 day orbital period. However, the hypothesis of a single companion with M$_{\rm B}$ $\sim$188 M$_{\rm jup}$ and P$_{\rm orb}$ $\sim$247 days can fit both datasets and cannot be completely excluded at this stage of the analysis. Conclusions. At minima, HQ Tau is a single-lined spectroscopic binary, and several factors indicate that the companion is a brown dwarf and that a third component is responsible for larger RV variation on a longer timescale.

On the route towards merging neutron stars and stripped-envelope supernovae, binary population synthesis predicts a large number of post-interaction systems with massive stars that have stripped off their outer layers. Yet, observations of such stars in the intermediate-mass regime below the Wolf-Rayet masses are rare. Using X-Shooting ULLYSES (XShootU) data, we discovered three partially stripped star + Be/Oe binaries in the Magellanic Clouds. We analyzed the UV and optical spectra using the PoWR model atmosphere code by superimposing model spectra corresponding to each component. The estimated current masses of the partially stripped stars fall within the intermediate mass range of 4-8 $M_{\odot}$. These objects are overluminous for their stellar masses, matching core He-burning luminosities. Their Be/Oe secondaries have much higher masses than their stripped primaries (mass ratio > 2). All three partially stripped stars show significant nitrogen enrichment and carbon and oxygen depletion on their surfaces. Additionally, one of our sample stars exhibits significant helium enrichment. Our study provides the first comprehensive determination of the wind parameters of partially stripped stars in the intermediate mass range. The wind mass-loss rates of these stars are found to be on the order of $10^{-7} M_\odot$ yr$^{-1}$, which is over ten times higher than that of OB stars of the same luminosity. Current evolutionary models characterizing this phase typically employ OB or WR mass-loss rates, which underestimate or overestimate stripped stars' mass-loss rates by an order of magnitude. Binary evolution models indicate that the observed primaries had initial masses of 12-17 $M_{\odot}$, making them potential candidates for stripped-envelope supernovae that form neutron stars. If they survive the explosion, these systems may become Be X-ray binaries and later double neutron stars.

In our previous study, we introduced a machine-learning technique, namely CMBFSCNN, for the removal of foreground contamination in cosmic microwave background (CMB) polarization data. This method was successfully employed on actual observational data from the Planck mission. In this study, we extend our investigation by considering the CMB lensing effect in simulated data and utilizing the CMBFSCNN approach to recover the CMB lensing B-mode power spectrum from multi-frequency observational maps. Our method is first applied to simulated data with the performance of CMB-S4 experiment. We achieve reliable recovery of the noisy CMB Q (or U) maps with a mean absolute difference of $0.016\pm0.008\ \mu$K (or $0.021\pm0.002\ \mu$K) for the CMB-S4 experiment. To address the residual instrumental noise in the foreground-cleaned map, we employ a "half-split maps" approach, where the entire dataset is divided into two segments sharing the same sky signal but having uncorrelated noise. Using cross-correlation techniques between two recovered half-split maps, we effectively reduce instrumental noise effects at the power spectrum level. As a result, we achieve precise recovery of the CMB EE and lensing B-mode power spectra. Furthermore, we also extend our pipeline to full-sky simulated data with the performance of LiteBIRD experiment. As expected, various foregrounds are cleanly removed from the foreground contamination observational maps, and recovered EE and lensing B-mode power spectra exhibit excellent agreement with the true results. Finally, we discuss the dependency of our method on the foreground models.

The North-South asymmetry in the number density and bulk velocity of stars in the solar neighborhood provides valuable insights into the formation and evolution of the Milky Way disk. Our objective is to investigate the wave-like disk oscillations of mono-age stellar populations in the Solar neighbourhood using data from Gaia Data Release 3. We have selected a comprehensive sample of main sequence turn off stars. The ages of these stars can be accurately determined using isochrone fitting methods. Our findings indicate that the north-south density and mean vertical velocity asymmetries remain consistent across all age groups.The uniformity of perturbations across all subsamples suggests that all populations are responding to the same external influence, which likely affects them irrespective of their age. Moreover, the fact that these perturbations appear consistently implies they could be either ongoing or recent. Regarding vertical velocity dispersions, we observe that older stars exhibit larger dispersions.

Spectropolarimetry is a powerful tool to investigate the central regions of active galactic nuclei (AGNs) as polarization signatures are key to probing magnetic field structure, evolution, and the physics of particle acceleration in jets. Optical linear polarization of blazars is typically greater than a few percent, indicating the emission is dominated by nonthermal synchrotron radiation, while polarization less than a few percent is common for other type 1 AGNs. We present a spectropolarimetric study of PKS 0637-75 and PKS 1510-089 to determine how the head-on orientation of a jet and dominant emission processes influence polarimetric variations in the broad lines and continuum. Observations were obtained biweekly from the Robert Stobie Spectrograph on the Southern African Large Telescope. Variability in the continuum polarization is detected for both PKS 0637-75 and PKS 1510-089, with a total average level of 2.5% +/- 0.1% and 7.5% +/- 0.1%, respectively. There is no clear polarization in the broad Balmer emission lines and weak polarization in Mg II as the average level across all observations is 0.2% +/- 0.1% for Hbeta, 0.2% +/- 0.3% for Hgamma, and 0.6% +/- 0.2% for Mg II. We find that polarization measurements confirm the conclusions drawn from spectral energy distribution modeling of the disk-jet contributions to the emission as optical polarization and time variability for PKS 0637-75 are shown to be dominated by accretion disk emission while those of PKS 1510-089 are due to both disk and jet emission, with greater jet contribution during flaring states.

We present the design and validation of a variable temperature cryogenic blackbody source, hereinafter called a cold load, that will be used to characterize detectors to be deployed by CMB-S4, the next-generation ground-based cosmic microwave background (CMB) experiment. Although cold loads have been used for detector characterization by previous CMB experiments, this cold load has three novel design features: (1) the ability to operate from the 1 K stage of a dilution refrigerator (DR), (2) a 3He gas-gap heat switch to reduce cooling time, and (3) the ability to couple small external optical signals to measure detector optical time constants under low optical loading. The efficacy of this design was validated using a 150 GHz detector array previously deployed by the Spider experiment. Thermal tests showed that the cold load can be heated to temperatures required for characterizing CMB-S4's detectors without significantly impacting the temperatures of other cryogenic stages when mounted to the DR's 1 K stage. Additionally, optical tests demonstrated that external signals can be coupled to a detector array through the cold load without imparting a significant optical load on the detectors, which will enable measurements of the CMB-S4 detectors' optical time constants.

The neutrino signal from the next galactic core-collapse supernova will provide an invaluable early warning of the explosion. By combining the burst trigger from several neutrino detectors, the location of the explosion can be triangulated minutes to hours before the optical emission becomes visible, while also reducing the rate of false-positive triggers. To enable multi-messenger follow-up of nearby supernovae, the SuperNova Early Warning System 2.0 (SNEWS 2.0) will produce a combined alert using a global network of neutrino detectors. This paper describes the trigger publishing and alert formation framework of the SNEWS 2.0 network. The framework is built on the HOPSKOTCH publish-subscribe system to easily incorporate new detectors into the network, and it implements a coincidence system to form alerts and estimate a false-positive rate for the combined triggers. The paper outlines the structure of the SNEWS 2.0 software and the initial testing of coincident signals.

I examine recent theoretical studies and observations of the recent core-collapse supernova (CCSN) SN 20224ggi and find that the likely explanation for its dense, compact circumstellar material is an effervescent model, where parcels or streams of gas are uplifted by stellar convection and pulsation and fall back. The effervescent zone exists alongside the regular wind from the red supergiant (RSG) progenitor of SN 2024ggi. I find that an extended wind-acceleration zone encounters some difficulties in accounting for the required CSM mass. Recent modelling finds the explosion energy of SN 2024ggi to be >1e51 erg and up to 2e51 erg. I examine this explosion energy against a recent study of the delayed neutrino explosion mechanism and find that this mechanism has difficulties in accounting for the required energy. This suggests that the explosion was caused by the jittering jets explosion mechanism (JJEM), adding to other recent pieces of evidence supporting the JJEM, particularly point-symmetric CCSN remnants.

The Event Horizon Telescope has produced resolved images of the supermassive black holes Sgr A* and M87*, which present the largest shadows on the sky. In the next decade, technological improvements and extensions to the array will enable access to a greater number of sources, unlocking studies of a larger population of supermassive black holes through direct imaging. In this paper, we identify 12 of the most promising sources beyond Sgr A* and M87* based on their angular size and millimeter flux density. For each of these sources, we make theoretical predictions for their observable properties by ray tracing general relativistic magnetohydrodynamic models appropriately scaled to each target's mass, distance, and flux density. We predict that these sources would have somewhat higher Eddington ratios than M87*, which may result in larger optical and Faraday depths than previous EHT targets. Despite this, we find that visibility amplitude size constraints can plausibly recover masses within a factor of 2, although the unknown jet contribution remains a significant uncertainty. We find that the linearly polarized structure evolves substantially with Eddington ratio, with greater evolution at larger inclinations, complicating potential spin inferences for inclined sources. We discuss the importance of 345 GHz observations, milli-Jansky baseline sensitivity, and independent inclination constraints for future observations with upgrades to the Event Horizon Telescope (EHT) through ground updates with the next-generation EHT (ngEHT) program and extensions to space through the Black Hole Explorer (BHEX).

We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~\mu{\text{m}}$) and its evolution with redshift up to $z=2.5$. Sérsic profiles are measured for $\sim$ $26\,000$ galaxies with stellar masses $M_\star > 10^9~{\text{M}}_\odot$ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys, using coordinates, redshifts, colors and stellar mass estimates from the COSMOS2020 catalog. The new rest-frame near-IR effective radii are generally smaller than previously measured rest-frame optical sizes, on average by 0.14~dex, with no significant dependence on redshift. For quiescent galaxies this size offset does not depend on stellar mass, but for star-forming galaxies the offset increases from -0.1~dex at $M_\star = 10^{9.5}~{\text{M}}_\odot$ to -0.25~dex at $M_\star > 10^{11}~{\text{M}}_\odot$. That is, we find that the near-IR stellar mass - size relation for star-forming galaxies is flatter in the rest-frame near-IR than in the rest-frame optical at all redshifts $0.5<z<2.5$. The general pace of size evolution is the same in the near-IR as previously demonstrated in the optical, with slower evolution ($R_{\text{e}} \propto (1+z)^{-0.7}$) for $L^*$~star-forming galaxies and faster evolution ($R_{\text{e}} \propto (1+z)^{-1.3}$) for $L^*$~quiescent galaxies. Massive ($M_\star>10^{11}~{\text{M}}_\odot$) star-forming galaxies evolve in size almost as fast as quiescent galaxies. Low-mass ($M_\star<10^{10}~{\text{M}}_\odot$)~quiescent galaxies evolve as slow as star-forming galaxies. Our main conclusion is that the size evolution narrative as it has emerged over the past two decades does not radically change when accessing with JWST the rest-frame near-IR, a better proxy of the underlying stellar mass distribution.

In this study, spectral, age, kinematic, and orbital dynamical analyses were conducted on metal-poor and high proper-motion (HPM) stars, HD 8724 and HD 195633, selected from the Solar neighborhood. This analysis combines detailed abundance measurements, kinematics, and orbital dynamics to determine their origin. Standard 1D local thermodynamic equilibrium analysis provides a fresh determination of the atmospheric parameters: $T_{\rm eff}=$4700$\pm$115 K, $\log g=$ 1.65$\pm$0.32 cgs, [Fe/H]=-1.59$\pm$0.04 dex, and a microturbulent velocity $\xi=$ 1.58$\pm$0.50 km s$^{\rm -1}$ for HD 8724 and $T_{\rm eff}=$6100$\pm$205 K, $\log g=$3.95$\pm$0.35 cgs, [Fe/H]=-0.52$\pm$0.05 dex, and $\xi=$1.26$\pm$0.50 km s$^{\rm -1}$ for HD 195633. The ages were estimated using a Bayesian approach (12.25 Gyr for HD 8724 and 8.15 Gyr for HD 195633). The escape scenarios of these stars from 170 candidate globular clusters (GCs) in the Galaxy were also investigated because of their chemical and physical differences (HPM and metal-poor nature). Accordingly, the calculated probability of encounter ($59\%$) for HD 8724 at a distance of five tidal radius suggests that star HD 8724 may have escaped from NGC 5139 ($\omega$ Cen), supported by its highly flattened orbit and may belong to a sub-population of this GC. Conversely, HD 195633's kinematics, age, and metal abundances point towards an escape from the bulge GC NGC 6356.