Papers for Friday, Sep 20 2024

The redshift-dependent relation between galaxy stellar mass and star formation rate (the Star-Forming Sequence, or SFS) is a key observational yardstick for galaxy assembly. We use the SAGAbg-A sample of background galaxies from the Satellites Around Galactic Analogs (SAGA) Survey to model the low-redshift evolution of the low-mass SFS. The sample is comprised of 23258 galaxies with H$\alpha$-based star formation rates (SFRs) spanning $6<\log_{10}(\rm M_\star/[M_\odot])<10$ and $z<0.21$ ($t<2.5$ Gyr). Although it is common to bin or stack galaxies at $z \lesssim 0.2$ for galaxy population studies, the difference in lookback time between $z=0$ and $z=0.21$ is comparable to the time between $z=1$ to $z=2$. We develop a model to account for both the physical evolution of low-mass SFS and the selection function of the SAGA survey, allowing us to disentangle redshift evolution from redshift-dependent selection effects across the SAGAbg-A redshift range. Our findings indicate significant evolution in the SFS over the last 2.5 Gyr, with a rising normalization: $\langle {\rm SFR}({\rm M_\star=10^{8.5} M_\odot)}\rangle(z)=1.24^{+0.25}_{-0.23}\ {\rm z} -1.47^{+0.03}_{-0.03}$. We also identify the redshift limit at which a static SFS is ruled out at the 95% confidence level, which is $z=0.05$ based on the precision of the SAGAbg-A sample. Comparison with cosmological hydrodynamic simulations reveals that some contemporary simulations under-predict the recent evolution of the low-mass SFS. This demonstrates that the recent evolution of the low-mass SFS can provide new constraints on the assembly of the low-mass Universe and highlights the need for improved models in this regime.

Numerous studies have reported significant displacements in the coordinates of active galactic nuclei between measurements using radio interferometry techniques and those obtained by the Gaia space observatory. There is a consensus that these discrepancies do indeed manifest astrometrically resolved sub-components of AGN rather than random measurement noise. Among other evidence, it has been reported that AGN with VLBI-to-Gaia displacements (VGD) pointing downstream their parsec-scale radio jets exhibit higher optical polarization compared to sources with the opposite VGD orientation. We aim to verify this previously reported connection between optical polarization and VGD-jet angle using a larger dataset of polarimetric measurements and updated Gaia DR3 positions. We also seek additional evidence supporting the disc-jet dichotomy as an explanation of such a connection by using mm-wave polarization and multi-band optical polarization measurements. We performed optical polarimetric observations of 157 AGN using three telescopes. We cross-match public astrometric data from VLBI and Gaia experiments, obtain corresponding positional displacements, and join this catalog with the polarimetric and jet directions data. AGN with downstream VGD are confirmed to have significantly higher optical fractional polarization than the upstream sample. At the same time, the mm-wavelength polarization of the two samples demonstrates very similar distributions. Our results support the hypothesis that the VLBI-to-Gaia displacements pointing down the radio jet are likely caused by a component in the jet emitting highly polarized synchrotron radiation and dominating in the overall optical emission. The upstream-oriented VGD are likely to be produced by the low-polarized emission of the central engine sub-components dominating in the optical.

Interferometric observations of gravitational microlensing events offer an opportunity for precise, efficient, and direct mass and distance measurements of lensing objects, especially those of isolated neutron stars and black holes. However, such observations were previously possible for only a handful of extremely bright events. The recent development of a dual-field interferometer, GRAVITY Wide, has made it possible to reach out to significantly fainter objects, and increase the pool of microlensing events amenable to interferometric observations by two orders of magnitude. Here, we present the first successful observation of a microlensing event with GRAVITY Wide and the resolution of microlensed images in the event OGLE-2023-BLG-0061/KMT-2023-BLG-0496. We measure the angular Einstein radius of the lens with a sub-percent precision, $\theta_{\rm E} = 1.280 \pm 0.009$ mas. Combined with the microlensing parallax detected from the event light curve, the mass and distance to the lens are found to be $0.472 \pm 0.012 M_{\odot}$ and $1.81 \pm 0.05$ kpc, respectively. We present the procedure for the selection of targets for interferometric observations, and discuss possible systematic effects affecting GRAVITY Wide data. This detection demonstrates the capabilities of the new instrument and it opens up completely new possibilities for the follow-up of microlensing events, and future routine discoveries of isolated neutron stars and black holes.

We present a velocity-resolved reverberation mapping analysis of the hypervariable quasar RM160 (SDSS J141041.25+531849.0) at z = 0.359 with 153 spectroscopic epochs of data representing a ten-year baseline (2013-2023). We split the baseline into two regimes based on the 3x flux increase in the light curve: a 'low state' phase during the years 2013-2019 and a 'high state' phase during the years 2022-2023. The velocity-resolved lag profiles (VRLP) indicate that gas with different kinematics dominates the line emission in different states. The H\b{eta} VRLP begins with a signature of inflow onto the BLR in the 'low state', while in the 'high state' it is flatter with less signature of inflow. The H{\alpha} VRLP begins consistent with a virialized BLR in the 'low state', while in the 'high state' shows a signature of inflow. The differences in the kinematics between the Balmer lines and between the 'low state' and the 'high state' suggests complex BLR dynamics. We find that the BLR radius and velocity (both FWHM and {\sigma}) do not obey a constant virial product throughout the monitoring period. We find that BLR lags and continuum luminosity are correlated, consistent with rapid response of the BLR gas to the illuminating continuum. The BLR kinematic profile changes in unpredictable ways that are not related to continuum changes and reverberation lag. Our observations indicate that non-virial kinematics can significantly contribute to observed line profiles, suggesting caution for black-hole mass estimation in luminous and highly varying quasars like RM160.

Identifying infall motions is crucial for our understanding of accretion processes in regions of star formation. The NH3 (1,1) hyperfine intensity anomaly (HIA) has been proposed to be a readily usable tracer for such infall motions in star-forming regions harboring young stellar objects at very early evolutionary stages. In this paper, we seek to study the HIA toward fifteen infall candidate regions to assess its reliability as an infall tracer. By using deep observations of the NH3 (1,1) transition with the Effelsberg 100 m telescope, HIAs have been identified toward all the targets. Fourteen out of fifteen sources exhibit anomalous intensities either in the inner or outer satellite lines. All the derived HIAs conform to the framework of the existing two models, namely, hyperfine selective trapping (HST) and systematic contraction or expansion motion (CE) models. In our sample of infall candidates, a majority of the HIAs remain consistent with the HST model. Only in three targets, the HIAs are consistent with infall motions under the CE model. Thus HIAs could be used as an infall tracer but seem not highly sensitive to infall motions in our single-dish data. Nevertheless, the emission could be blended with emission from outflow activities. HIAs consistent with the HST model show stronger anomalies with increasing kinetic temperatures (Tk), which is expected by the HST model. On the other hand, HIAs consistent with infall motions show little dependence on Tk. Therefore, HIAs may preferably trace infall of cold gas.

We present one of the largest multiwavelength studies of simultaneous optical-to-X-ray spectral energy distributions (SEDs) of unobscured active galactic nuclei (AGN) in the local Universe. Using a representative sample of hard-X-ray-selected AGN from the 70-month Swift/BAT catalog, with optical/UV photometric data from Swift/UVOT and X-ray spectral data from Swift/XRT, we constructed broadband SEDs of 236 nearby AGN (0.001<z<0.3). We employed GALFIT to estimate host galaxy contamination in the optical/UV and determine the intrinsic AGN fluxes. We used an absorbed power law with a reflection component to model the X-ray spectra and a dust-reddened multi-temperature blackbody to fit the optical/UV SED. We calculated total bolometric luminosities ($L_{bol}$), optical-to-X-ray spectral indices ($\alpha_{ox}$), and multiple bolometric corrections (BCs) in the optical, UV, and X-rays. We used black hole masses obtained by reverberation mapping and the virial method to estimate Eddington ratios ($\lambda_{Edd}$) for all our AGN. We confirm the tight correlation between UV and X-ray luminosity for our sample. We observe a significant decrease in $\alpha_{ox}$ with $L_{bol}$ and $\lambda_{Edd}$, suggesting that brighter sources emit more UV photons per X-rays. We report a second-order regression relation between the 2-10 keV BC and $\alpha_{ox}$, which is useful to compute $L_{bol}$ in the absence of multiband SEDs. We also investigate the dependence of optical/UV BCs on the physical properties of AGN and obtain a significant increase in the UV BCs with $L_{bol}$ and $\lambda_{Edd}$, unlike those in the optical, which are constant across five orders of $L_{bol}$ and $\lambda_{Edd}$. We obtain significant dispersions (~0.1-1 dex) in all BCs, and hence recommend using appropriate relations with observed quantities while including the reported scatter, instead of their median values.

We present a series of five new broadband X-ray observations of the ultraluminous X-ray source Holmberg IX X-1, performed by $XMM$-$Newton$ and $NuSTAR$ in coordination. The first three of these show high soft X-ray fluxes but a near total collapse of the high-energy ($\gtrsim$15 keV) emission, previously seen to be surprisingly stable across all prior broadband observations of the source. The latter two show a recovery in hard X-rays, remarkably once again respecting the same stable high-energy flux exhibited by all of the archival observations. We also present a joint analysis of all broadband observations of Holmberg IX X-1 to date (encompassing 11 epochs in total) in order to investigate whether it shows the same luminosity-temperature behaviour as NGC 1313 X-1 (which also shows a stable high-energy flux), whereby the hotter disc component in the spectrum exhibits two distinct, positively-correlated tracks in the luminosity-temperature plane. Holmberg IX X-1 may show similar behaviour, but the results depend on whether the highest energy emission is assumed to be an up-scattering corona or an accretion column. The strongest evidence for this behaviour is found in the former case, while in the latter the new 'soft' epochs appear distinct from the other high-flux epochs. We discuss possible explanations for these new 'soft' spectra in the context of the expected structure of super-Eddington accretion flows around black holes and neutron stars, and highlight a potentially interesting analogy with the recent destruction and re-creation of the corona seen in the AGN 1ES 1927+654.

The early growth of black holes (BHs) in atomic-cooling halos is likely influenced by feedback on the surrounding gas. While the effects of radiative feedback are well-documented, mechanical feedback, particularly from AGN jets, has been comparatively less explored. Building on our previous work that examined the growth of a 100 ${M_\odot}$ BH in a constant density environment regulated by AGN jets, we expand the initial BH mass range from 1 to $10^4$ ${M_\odot}$ and adopt a more realistic density profile for atomic-cooling halos. We reaffirm the validity of our analytic models for jet cocoon propagation and feedback regulation. We identify several critical radii-namely, the terminal radius of jet cocoon propagation, the isotropization radius of the jet cocoon, and the core radius of the atomic-cooling halo-that are crucial in determining BH growth given specific gas properties and jet feedback parameters. In a significant portion of the parameter space, our findings show that jet feedback substantially disrupts the halo's core during the initial feedback episode, preventing BH growth beyond $10^4$ ${M_\odot}$. Conversely, conditions characterized by low jet velocities and high gas densities enable sustained BH growth over extended periods. We provide a prediction for the black hole mass growth as a function of time and feedback parameters. We found that, to form a supermassive BH ($>10^6 {M_\odot}$) within 1 Gyr entirely by accreting gas from an atomic-cooling halo, the jet energy feedback efficiency must be $\lesssim 10^{-4} \dot{M}_{BH} c^2$ even if the seed BH mass is $10^4 {M_\odot}$.

Studying the orbital motion of stars around Sagittarius A* in the Galactic Center provides a unique opportunity to probe the gravitational potential near the supermassive black hole at the heart of our Galaxy. Interferometric data obtained with the GRAVITY instrument at the Very Large Telescope Interferometer (VLTI) since 2016 has allowed us to achieve unprecedented precision in tracking the orbits of these stars. GRAVITY data have been key to detecting the in-plane, prograde Schwarzschild precession of the orbit of the star S2, as predicted by General Relativity. By combining astrometric and spectroscopic data from multiple stars, including S2, S29, S38, and S55 - for which we have data around their time of pericenter passage with GRAVITY - we can now strengthen the significance of this detection to an approximately $10 \sigma$ confidence level. The prograde precession of S2's orbit provides valuable insights into the potential presence of an extended mass distribution surrounding Sagittarius A*, which could consist of a dynamically relaxed stellar cusp comprised of old stars and stellar remnants, along with a possible dark matter spike. Our analysis, based on two plausible density profiles - a power-law and a Plummer profile - constrains the enclosed mass within the orbit of S2 to be consistent with zero, establishing an upper limit of approximately $1200 \, M_\odot$ with a $1 \sigma$ confidence level. This significantly improves our constraints on the mass distribution in the Galactic Center. Our upper limit is very close to the expected value from numerical simulations for a stellar cusp in the Galactic Center, leaving little room for a significant enhancement of dark matter density near Sagittarius A*.

The new tools of gravitational wave and multi-messenger astronomy allow for the study of astrophysical phenomenon in new ways and enables light to be shed on some of the longest-enduring mysteries of high-energy astrophysics. Among the latter stands the Galactic center gamma-ray excess, associated with a source whose nature could be annihilating dark matter or a yet-unresolved population of millisecond pulsars (MSPs). MSPs are most likely asymmetric about their axis of rotation, and are thus thought to also source quasi-monochromatic gravitational waves, that dark matter processes would not emit. Using statistical methods, we simulate realistic MSP population samples with differing morphology and moment of inertia, that could give rise to the gamma-ray excess, and we compute the corresponding gravitational wave signal amplitude and frequency. We find that the gravitational wave signal frequency likely ranges between $\sim$200 and 1400 Hz, and that the collective dimensionless strain from the center of the Galaxy has an amplitude between $10^{-26}$ and $10^{-24}$, thus most likely beyond current and near-term detectors, unless the unresolved MSPs are extraordinarily gamma-ray dim.

Relic galaxies, the oldest ultra-compact massive galaxies (UCMGs), contain almost exclusively "pristine" stars formed during an intense star formation (SF) burst at high redshift. As such, they allow us to study in detail the early mechanism of galaxy assembly in the Universe. Using the largest catalogue of spectroscopically confirmed UCMGs for which a degree of relicness (DoR) had been estimated, the INSPIRE catalogue, we investigate whether or not relics prefer dense environments. The objective of this study is to determine if the DoR, which measures how extreme the SF history was, and the surrounding environment are correlated. In order to achieve this goal, we employ the AMICO galaxy cluster catalogue to compute the probability for a galaxy to be a member of a cluster, and measure the local density around each UCMG using machine learning-based photometric redshifts. We find that UCMGs can reside both in clusters and in the field, but objects with very low DoR (< 0.3, i.e., a relatively extended SF history) prefer under-dense environments. We additionally report a correlation between the DoR and the distance from the cluster centre: more extreme relics, when located in clusters, tend to occupy the more central regions of them. We finally outline potential evolution scenarios for UCMGs at different DoR to reconcile their presence in both clusters and field environments

The Earth sits inside a 300pc-wide void that was carved by a series of supernova explosions that went off tens of millions of years ago, pushing away interstellar gas and creating a bubble-like structure. The 60Fe peak deposits found in the deep-sea crust have been interpreted by the imprints left by the ejecta of supernova explosions occurring about 2-3 and 5-6 Myr ago. It is likely that the 60Fe peak at about 2-3 Myr originated from a supernova occurring in the Upper Centaurus Lupus association in Scorpius Centaurus (140 pc) or the Tucana Horologium association (70 pc). Whereas, the 5-6 Myr peak is likely attributed to the solar system's entrance into the bubble. In this Letter, we show that the supernova source responsible for synthesizing the 60Fe peak deposits 2-3 Myr ago was also likely a Galactic PeVatron source. We demonstrate that this supernova can consistently explain the "knee" in the cosmic-ray spectrum and the large-scale anisotropy between 100 TeV and 100 PeV. Matching the intensity and shape of the cosmic-ray spectrum allows us to place stringent constraints on the cosmic-ray energy content from the supernova as well as on the cosmic-ray diffusion coefficient. Making use of such constraints we provide a robust estimate of the temporal variation of terrestrial ionizing cosmic radiation levels and discuss their implications in the development of early life on Earth by plausibly influencing the mutation rate and, as such, conceivably assisting in the evolution of complex organisms.

We analyze the variability of the parsec-scale jet directions in active galactic nuclei (AGNs). Our analysis involves 317 AGNs at frequencies ranging from 2 to 43 GHz, and is made possible by developing an automatic inner jet direction measurement procedure. We find strong significant variations in a one quarter of these AGNs; the effect is likely ubiquitous, and not detected in the rest due to a limited sensitivity and observations epoch coverage. Average apparent jet rotation speeds range from 0.21 deg/yr at 2 GHz to 1.04 deg/yr at 43 GHz. This strong frequency dependence indicates that the variability cannot be explained by jet components propagating ballistically without acceleration: more complex jet shapes or patterns are required. Still, we demonstrate that the apparent direction changes are predominantly caused by the jet nozzle rotations, and not by individual components propagating transversely to the jet. In this work, we focus on variability scales much longer than the times of observations, that is > 50 years. Using our measurements, we bound potential periods to less than 1000 years in the source rest frame for 90% AGNs in the sample. This allows us to constrain mechanisms causing these variations if they are periodic, such as instabilities, disk-driven precession, or binary black hole effects.

With close to three years of observations in hand, the NEID Earth Twin Survey (NETS) is starting to unearth new astrophysical signals for a curated sample of bright, radial velocity (RV)-quiet stars. We present the discovery of the first NETS exoplanet, HD 86728 b, a $m_p\sin i = 9.16^{+0.55}_{-0.56}\ \rm{M}_\oplus$ planet on a circular, $P=31.1503^{+0.0062}_{-0.0066}$ d orbit, thereby confirming a candidate signal identified by Hirsch et al. (2021). We confirm the planetary origin of the detected signal, which has a semi-amplitude of just $K=1.91^{+0.11}_{-0.12}$ m s$^{-1}$, via careful analysis of the NEID RVs and spectral activity indicators, and we constrain the mass and orbit via fits to NEID and archival RV measurements. The host star is intrinsically quiet at the $\sim1$ m s$^{-1}$ level, with the majority of this variability likely stemming from short-timescale granulation. HD 86728 b is among the small fraction of exoplanets with similar masses and periods that have no known planetary siblings.

Spitzer 'hidden' observations of the background are used to construct a catalog of 4,090 spectra and examine the signature of polycyclic aromatic hydrocarbon (PAH) molecules and their connection to extinction by dust. A strong positive correlation is recovered between WISE12, E(B-V), and the 11.2 $\mu$m PAH band. For 0.06 $\leq$ E(B-V) $\leq$ 5.0, correlations of the 6.2, 11.2, and 12.7 $\mu$m PAH band are positive with E(B-V). Three dust temperature regimes are revealed. Correlations with WISE12 are well-constrained and that with 12.7/11.2 is flat. Decomposition with the NASA Ames PAH IR Spectroscopic Database reveals a tentative positive correlation between the 6.2/11.2 and the PAH ionization fraction, while that with 12.7/11.2 is slightly negative, suggesting PAH structural changes. The relation with PAH size and 6.2/11.2 is negative, while that with 12.7/11.2 is positive. Averaging spectra into five E(B-V) and three T$_{\rm dust}$ bins shows an evolution in PAH emission and variations in 12.7/11.2. Database-fits show an increase in $f_{\rm i}$ and the PAH ionization parameter $\gamma$, but a more stable large PAH fraction. While the largest $\gamma$s are associated with the highest T$_{\rm dust}$, there is no one-to-one correlation. The analysis is hampered by low-quality data at short wavelengths. There are indications that PAHs in the more-diffuse backgrounds behave differently from those in the general interstellar medium. However, they are often still associated with larger scale filamentary cloud-like structures. The spectra and auxiliary data have been made available through the Ames Background Interstellar Medium Spectral Catalog and may guide JWST programs.

The Compton Spectrometer and Imager (COSI) is an upcoming NASA Small Explorer satellite mission scheduled for launch in 2027 and designed to conduct an all-sky survey in the energy range of 0.2-5 MeV. Its instrument consists of an array of germanium detectors surrounded on four sides and underneath by active shields that work as anticoincidence system (ACS) to reduce the contribution of background events in the detectors. These shields are composed of bismuth germanium oxide (BGO), a scintillator material, coupled with Silicon photomultipliers, aimed to collect optical photons produced from interaction of ionizing particles in the BGO and convert them into an electric signal. The reference simulation framework for COSI is MEGAlib, a set of software tools based on the Geant4 toolkit. The interaction point of the incoming radiation, the design of the ACS modules and the BGO surface treatment change the light collection and the overall shielding accuracy. The use of the Geant4 optical physics library, with the simulation of the scintillation process, is mandatory for a more realistic evaluation of the ACS performances. However, including the optical processes in MEGAlib would dramatically increase the computing time of the COSI simulations. We propose the use of a response function encoding the energy resolution and 3D light yield correction based on a separate Geant4 simulation of the ACS that includes the full optical interaction. We present the verification of the Geant4 optical physics library against analytical computations and available laboratory measurements obtained using PMTs as readout device, as a preparatory phase for the simulation of the COSI ACS response.

The dynamics of a black hole traveling through a plasma -- a general relativistic extension of the classic Bondi-Hoyle-Lyttleton (BHL) accretion problem -- can be related to a variety of astrophysical contexts, including the aftermath of binary black hole mergers in gaseous environments. We perform three-dimensional general relativistic magnetohydrodynamics simulations of BHL accretion onto a rotating black hole for an incoming flow with a toroidal (inclined) magnetization with respect to the spin axis of the black hole. Irrespective of inclination but dependent on the wind speed, we find that the accretion flow onto the black hole can become magnetically arrested, launching an intermittent (and sometimes striped) jet. The upstream ram pressure of the wind bends the jet, and confines the angular extent into which the magnetic flux tubes ejected from quasi-periodic eruptions are released. Recoil from magnetic flux eruptions drives strong oscillations in the accretion plane, resulting in jet nutation at the outer radii and occasionally ripping off the inner part of the accretion disk. In addition to dynamical friction, the black hole experiences a perpendicular drag force analogous to the Magnus effect. Qualitative effects of the incoming magnetic field orientation, the strength of the magnetization, and the incoming wind speed are investigated as well.

Clusters of galaxies can be identified from peaks in weak lensing aperture mass maps constructed from weak lensing shear catalogs. Such purely gravitational cluster selection considerably differs from traditional cluster selections based on baryonic properties of clusters. In this review, we present the basics and applications of weak lensing shear-selected cluster samples. Detailed studies of baryonic properties of shear-selected clusters shed new light on cluster astrophysics. The purely gravitational selection suggests that the selection function can be quantified more easily and robustly, which is crucial for deriving accurate cosmological constraints from the abundance of shear-selected clusters. The recent advance of shear-selected cluster studies is driven by the Subaru Hyper Suprime-Cam survey, in which more than 300 shear-selected clusters with the signal-to-noise ration greater than 5 are identified. It is argued that various systematic effects in the cosmological analysis can be mitigated by carefully choosing the set-up of the analysis, including the choice of the kernel functions and the source galaxy sample.

The abundance and distribution of metal in asteroid surfaces can be constrained from thermal emission measurements at radio wavelengths, informing our understanding of planetesimal differentiation processes. We observed the M-type asteroid (22) Kalliope and its moon Linus in thermal emission at 1.3, 9, and 20 mm with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA) over most of Kalliope's rotation period. The 1.3 mm data provide ~30 km resolution on the surface of Kalliope, while both the 1.3 and 9 mm data resolve Linus from Kalliope. We find a thermal inertia for Kalliope of 116$^{+326}_{-91}$ J m$^{-2}$ s$^{-0.5}$ K$^{-1}$ and emissivities of 0.65$\pm$0.02 at 1.3 mm, 0.56$\pm$0.03 at 9 mm, and 0.77$\pm$0.02 at 20 mm. Kalliope's millimeter wavelength emission is suppressed compared to its centimeter wavelength emission, and is also depolarized. We measure emissivities for Linus of 0.73$\pm$0.04 and 0.85$\pm$0.17 at 1.3 and 9 mm respectively, indicating a less metal-rich surface composition for Linus. Spatial variability in Kalliope's emissivity reveals a region in the northern hemisphere with a high dielectric constant, suggestive of enhanced metal content. These results are together consistent with a scenario in which Linus formed from reaggregated ejecta from an impact onto a differentiated Kalliope, leaving Kalliope with a higher surface metal content than Linus, which is distributed heterogeneously across its surface. The low emissivity and lack of polarization suggest a reduced regolith composition where iron is in the form of metallic grains and constitutes ~25% of the surface composition.

Polarization of electromagnetic waves carries a large amount of information about their astrophysical emitters and the media they passed through, and hence is crucial in various aspects of astronomy. Here we demonstrate an important but long-overlooked depolarization mechanism in astrophysics: when the polarization vector of light travels along a non-planar curve, it experiences an additional rotation, in particular for radio waves. The process leads to depolarization, which we call `geometric' depolarization (GDP). We give a concise theoretical analysis of the GDP effect on the transport of radio waves in a randomly inhomogeneous plasma under the geometrical optics approximation. In the case of isotropic scattering in the coronal plasma, we show that the GDP of the angle-of-arrival of the linearly polarized radio waves propagating through the turbulent plasma cannot be ignored. The GDP effect of linearly polarized radio waves can be generalized to astrophysical phenomena, such as fast radio bursts and stellar radio bursts, etc. Our findings may have a profound impact on the analysis of astrophysical depolarization phenomena.

Most stars are born in the crowded environments of gradually forming star clusters. Dynamical interactions between close-passing stars and the evolving UV radiation fields from proximate massive stars are expected to sculpt the protoplanetary discs in these clusters, potentially contributing to the diversity of planetary systems that we observe. Here, we investigate the impact of cluster environment on disc demographics by implementing simple protoplanetary disc evolution models within $N$-body simulations of gradual star cluster formation. We consider a range of star formation efficiency per free-fall time, $\epsilon_{\rm ff}$, and mass surface density of the natal cloud environment, $\Sigma_{\rm cl}$, both of which affect the overall duration of cluster formation. We track the interaction history of all stars to estimate the dynamical truncation of the discs around stars involved in close encounters. We also track external photoevaporation of the discs due to the ionizing radiation field of the nearby high- and intermediate-mass ($> 5 M_\odot$) stars. We find that $\epsilon_{\rm ff}$, $\Sigma_{\rm cl}$, and the degree of primordial binarity have major influences on the masses and radii of the disc population. In particular, external photo-evaporation has a greater impact than dynamical interactions in determining the fate of discs in our clusters.

A fraction of Galactic stars have compact companions which could be white dwarfs (WDs), neutron stars (NSs) or stellar-mass black holes (SBHs). In a detached and edge-on binary system including a main-sequence star and a compact object (denoted by WDMS, NSMS, and BHMS systems), the stellar brightness can change periodically due to self-lensing or eclipsing features. The shape of a self-lensing signals is a degenerate function of stellar radius and compact object's mass because the self-lensing peak strongly depends on the projected source radius normalized to Einstein radius. Increasing the inclination angle $i$ changes the self-lensing shape from a strict top-hat model to one with slow-increasing edges. We simulate stellar light curves due to these binary systems which are observed by NASA's Transiting Exoplanet Survey Satellite (TESS) telescope and evaluate the efficiencies to detect their periodic signatures using two sets of criteria (i)SNR$>3$ and $N_{\rm{tran}}>1$ (Low-Confidence, LC), and (ii) SNR$>5$ and $N_{\rm{tran}}>2$ (High-Confidence, HC). The HC efficiencies for detecting WDMS, NSMS, and BHMS systems with the inclination angle $i<20^{\circ}$ during different time spans are $5$-$7\%$, $4.5$-$6\%$, and $4$-$5\%$, respectively. Detecting lensing-induced features is possible in only $\lesssim3\%$ and $\lesssim33\%$ of detectable WDMS and NSMS events. The detection efficiencies for closer source stars with higher priorities are higher and drop to zero for $b\gtrsim R_{\star}$, where $b\simeq \tan(i) a$ is the impact parameter($a$ is the semi-major axis). We predict the numbers of WDs, NSs, and SBHs that are discovered from the TESS Candidate Target List stars are $15$-$18$, $6$-$7$, and $<1$.

The MMT Adaptive optics exoPlanet characterization System (MAPS) is currently in its engineering phase, operating on-sky at the MMT Telescope on Mt. Hopkins in southern Arizona. The MAPS Adaptive Secondary Mirror's actuators are controlled by a closed loop modified PID control law and an open loop feed-forward law, which in combination allows for faster actuator response time. An essential element of achieving the secondary's performance goals involves the process of PID gain tuning. To start, we briefly discuss the design of the MAPS ASM and its actuators. We then describe the actuator positional control system and control law. Next, we discuss a few of the issues that make ASM tuning difficult. We then outline our initial attempts at tuning the actuator controllers and discuss the use of actuator positional power spectra for both tuning and determining the health and failure states of individual actuators. We conclude by presenting the results of our latest round of tuning configuration trials, which have been successful at decreasing mirror latency, increasing operational mirror modes and improving image PSF.

We present the Deciphering the Interplay between the Interstellar medium, Stars, and the Circumgalactic medium (DIISC) Survey. This survey is designed to investigate the correlations in properties between the circumgalactic medium (CGM), the interstellar medium (ISM), stellar distributions, and young star-forming regions. The galaxies were chosen to have a QSO sightline within 3.5 times the HI radii probing the disk-CGM interface. The sample contains 34 low-redshift galaxies with a median stellar mass of 10$^{10.45}~\rm M_{\odot}$ probed at a median impact parameter of $\rho=55~kpc$. The survey combines ultraviolet spectroscopic data from the Cosmic Origins Spectrograph aboard the Hubble Space Telescope with HI 21 cm hyperfine transition imaging with the Very Large Array (VLA), ultraviolet imaging from Galaxy Evolution Explorer (GALEX), and optical imaging and spectroscopy with the MMT and Vatican Advanced Technology Telescope. We describe the specific goals of the survey, data reduction, high-level data products, and some early results. We present the discovery of a strong inverse correlation, at a confidence level of 99.99%, between Lyman $\alpha$ equivalent width, $\rm W_{Ly\alpha}$, and impact parameter normalized by the HI radius ($\rho/R_{HI}$). We find $\rho/R_{HI}$ to be a better empirical predictor of Lyman $\alpha$ equivalent width than virial radius normalized impact parameter ($\rho/R_{vir}$) or parameterizations combining $\rho,~R_{vir}$, stellar mass, and star formation rate. We conclude that the strong anticorrelation between the Lyman $\alpha$ equivalent width and $\rho/R_{HI}$ indicates that the neutral gas distribution of the CGM is more closely connected to the galaxy's gas disk rather than its stellar and dark matter content.

The X-ray Timing Instrument as part of the Neutron Star Interior Composition Explorer has the potential to examine the time-domain properties of compact objects in regimes not explored by previous timing instruments, due to its combination of high effective area and timing resolution. We consider the effects of instrumental deadtime at a range of effective countrates in a series of observations of the X-ray binary GX 339-4 to determine what effect deadtime has on photometric and Fourier frequency-domain products. We find that there are no significant inconsistencies across the functional detectors in the instrument, and that in the regimes where instrumental deadtime is a limiting factor on observations that previous approaches to dealing with deadtime, as applied to RXTE and other detectors, are still appropriate, and that performing deadtime corrections to lightcurves before creating Fourier products are not necessary at the count rates considered in our analysis.

Context. Fast Radio Burst 20180916B is a repeating FRB whose activity window has a 16.34 day periodicity that also shifts and varies in duration with the observing frequency. Recently, arXiv:2205.09221 reported the FRB has started to show secular Rotation Measure (RM) increasing trend after only showing stochastic variability around a constant value of $-114.6$ rad m$^{-2}$ since its discovery. Aims. We aim to further study the RM variability of FRB 20180916B. The data comes from the ongoing campaigns of FRB 20180916B using the upgraded Giant Metrewave Radio Telescope (uGMRT). The majority of the observations are in Band 4, which is centered at 650 MHz with 200 MHz bandwidth. Methods. We apply a standard single pulse search pipeline to search for bursts. In total, we detect 116 bursts with $\sim$36 hours of on-source time spanning 1200 days, with two bursts detected during simultaneous frequency coverage observations. We develop and apply a polarization calibration strategy suited for our dataset. On the calibrated bursts, we use QU-fitting to measure RM. Lastly, we also measure various other properties such as rate, linear polarization fraction and fluence distribution. Results. Of the 116 detected bursts, we could calibrate 79 of them. From which, we observed in our early observations the RM continued to follow linear trend as modeled by arXiv:2205.09221. However, our later observations suggest the source switch from the linear trend to stochastic variations around a constant value of $-58.75$ rad m$^{-2}$. We also study cumulative rate against fluence and note that rate at higher fluences (> 1.2 Jy ms) scales as $\gamma = -1.09(7)$ whereas that at lower fluences (between 0.2 and 1.2 Jy ms) only scales as $\gamma = -0.51(1)$, meaning rate at higher fluence regime is steeper than at lower fluence regime.

Recent observations of high-redshift galaxies have revealed starburst galaxies with excessive amounts of nitrogen, well above that expected in standard evolutionary models. The Sunburst Arc galaxy, particularly its young and massive star cluster, represents the closest ($z=2.4$) and brightest of these as a strongly lensed object. In this work, we study the chemical history of this star cluster to determine the origin of the elevated gas-phase nitrogen using a chemical evolution model. Our model includes the enrichment of OB stars through stellar winds and core-collapse supernovae assuming that massive stars ($M>25$ $M_\odot$) collapse directly into black holes at the end of their lives. We fit the model parameters to the observed chemical abundances of the Sunburst Arc cluster: O/H, C/O, and N/O. We find that the observed chemical abundances can be explained by models featuring intense star formation events, characterized by rapid gas accretion and high star formation efficiencies. Additionally, the stellar population contributing to the gas enrichment must exclude Wolf-Rayet stars. These conditions might be present in other nitrogen-rich objects as their similar chemical abundances suggest a common history. As previous studies have proposed the presence of Wolf-Rayet stars in the new nitrogen-rich objects, further research using chemodynamic modeling is necessary to ascertain the true nature of these objects.

GRB 230812B, detected by the Gamma-Ray Integrated Detectors (GRID) constellation mission, is an exceptionally bright gamma-ray burst (GRB) with a duration of only 3 seconds. Sitting near the traditional boundary ($\sim$ 2 s) between long and short GRBs, GRB 230812B is notably associated with a supernova (SN), indicating a massive star progenitor. This makes it a rare example of a short-duration GRB resulting from stellar collapse. Our analysis, using a time-evolving synchrotron model, suggests that the burst has an emission radius of approximately $10^{14.5}$~cm. We propose that the short duration of GRB 230812B is due to the combined effects of the central engine's activity time and the time required for the jet to break through the stellar envelope. Our findings provide another case that challenges the conventional view that short-duration GRBs originate exclusively from compact object mergers, demonstrating that a broader range of durations exists for GRBs arising from the collapse of massive stars.

We present new results of a numerical study of the pumping of 4.8 GHz and 14.5 GHz maser of o-Formaldehyde in the presence of a free-free radiation field. It is shown that in the presence of a free-free radiation field inversion of not only the 4.8 GHz transition, but also the 14.5 GHz transition and other doublet state transitions occur. Further results are presented to illustrate how, as a consequence of the pumping scheme, the inversion of the 4.8 GHz and 14.5 GHz transitions respond to the free-free radiation fields associated with HII regions with different emission measures and levels of geometric dilution with respect to the masing region. We also discuss the criticism raised in the past by various authors against the pumping of the 4.8 GHz Formaldehyde masers by a free-free radiation field. It is argued that the rarity of the Formaldehyde masers is not to be ascribed to the pumping scheme but to other factors such as, e.g., the evolution of the associated HII region or the chemical evolution of the star forming region which determines the Formaldehyde abundance or a combination of both.

Solar irradiance variations across various timescales, from minutes to centuries, represents a potential natural driver of past regional and global climate cold phases. To accurately assess the Sun's effect on climate, particularly during periods of exceptionally low solar activity known as grand minima, an accurate reconstruction of solar forcing is essential. While direct measurements of Total Solar Irradiance (TSI) only began in the late 1970s with the advent of space radiometers, indirect evidence from various historical proxies suggests that the Sun's magnetic activity has undergone possible significant fluctuations over much longer timescales. Employing diverse and independent methods for TSI reconstruction is essential to gaining a comprehensive understanding of this issue. This study employs a semi-empirical model to reconstruct TSI over the past millennium. Our approach uses an estimated open solar magnetic field ($F_{o}$), derived from cosmogenic isotope data, as a proxy for solar activity. We reconstruct the cyclic variations of TSI, due to the solar surface magnetic features, by correlating $F_{o}$ with the parameter of active region functional form. Instead, we obtain the long-term TSI trend by applying the Empirical Mode Decomposition (EMD) algorithm to the reconstructed $F_{o}$ to filter out the 11-year and 22-year solar variability. We prepare a reconstructed TSI record, spanning 971 to 2020 CE. The estimated departure from modern TSI values occurred during the Spörer Minimum (around 1400 CE), with a decrease of approximately 2.3 $W m^{-2}$. A slightly smaller decline of 2.2 $W m^{-2}$ is reported during the Maunder Minimum, between 1645 and 1715 CE.

The anions C$_7$N$^-$ and C$_{10}$H$^-$ are the two longest of the linear (C,N)-bearing and (C,H)-bearing chains which have so far been detected in the Interstellar Medium. In order to glean information on their collision-induced rotational state-changing processes, we analyse the general features of new ab initio potentials describing the interaction of both linear anions with H$_2$, one of the most abundant partners in their ISM environment. We employ an artificial neural network fit of the reduced-dimensionality potential energy surface for C$_7$N$^-$...H$_2$ interaction and discuss in detail the spatial features in terms of multipolar radial coefficients. For the C$_{10}$H$^-$...H$_2$ interaction we use the initial grid of two dimensional raw points to generate by quadrature the Legendre expansion directly, further including the long-range terms as discussed in the main text. Quantum scattering calculations are employed to obtain rotationally inelastic cross sections, for collision energies in the range of 10$^{-4}$ to 400 cm$^{-1}$. From them we generate the corresponding inelastic rate coefficients as a function of temperature covering the range from 10 to 50 K. The results for the rate coefficients for the longest cyanopolyyne are compared with the earlier results obtained for the smaller terms of the same series, also in collision with H$_2$. We obtain that the inelastic rate coefficients for the long linear anions are all fairly large compared with the earlier systems. The consequences of such findings on their non-equilibrium rotational populations in interstellar environments are illustrated in our conclusions.

The Kepler space mission provided high-quality light curves for more than 16 000 red giants. The global stellar oscillations extracted from these light curves carry information about the interior of the stars. Several hundred red giants were found to have low amplitudes in their dipole modes (i.e. they are suppressed dipole-mode stars). A number of hypotheses (involving e.g. a magnetic field, binarity, or resonant mode coupling) have been proposed to explain the suppression of the modes, yet none has been confirmed. We aim to gain insight into the mechanism at play in suppressed dipole-mode stars by investigating the mode properties (linewidths, heights, and amplitudes) of the radial oscillation modes of red giants with suppressed dipole modes. We selected from the literature suppressed dipole-mode stars and compared the radial-mode properties of these stars to the radial-mode properties of stars in two control samples of stars with typical (i.e. non-suppressed) dipole modes. We find that the radial-mode properties of the suppressed dipole-mode stars are consistent with the ones in our control samples, and hence not affected by the suppression mechanism. From this we conclude that (1) the balance between the excitation and damping in radial modes is unaffected by the suppression, and by extrapolation the excitation of the non-radial modes is not affected either; and (2) the damping of the radial modes induced by the suppression mechanism is significantly less than the damping from turbulent convective motion, suggesting that the additional damping originates from the more central non-convective regions of the star, to which the radial modes are least sensitive.

Binary systems comprise approximately 15 per cent of the near-Earth asteroid population, yet thermal-infrared data are often interpreted for these bodies as if they are single objects. Thermal-IR light curves of binary asteroids (3905) Doppler and (175706) 1996 FG3 are analysed using an adaptation of the Advanced Thermophysical Model, deriving new constraints on their thermal inertias as $\Gamma = 114 \pm 31 \mathrm{J} \mathrm{m}^{-2} \mathrm{K}^{-1} \mathrm{s}^{-1/2}$ and $\Gamma = 142 \pm 6 \mathrm{J} \mathrm{m}^{-2} \mathrm{K}^{-1} \mathrm{s}^{-1/2}$, respectively. We determine that this adapted model is suitable for binary systems where their primary rotation to secondary orbit period ratios can be approximately characterised by small integers. Objects with more complex orbital states require a model with alternative temperature convergence methodologies. Thermal inertia is shown to have a strong effect on binary thermophysical light curve morphology, introducing significant modulations both inside and outside of mutual event times. The depth of eclipse events are shown to be suppressed at longer wavelengths due to the sensitivity to cooler parts of the surface, meanwhile surface roughness is shown to have little effect on the thermal light curve morphology. A proof of concept model for the (65803) Didymos system is demonstrated, showing how such a binary model could be used to study the system during the European Space Agency's Hera mission, and the applicability of this adapted model to NASA's Lucy mission is also briefly discussed.

During the core collapse of a massive star just before its supernova explosion, the amplification of asymmetric motions by the standing accretion shock instability (SASI) imprints on the neutrino flux and the gravitational waves a frequency signature carrying direct information on the explosion process. The physical interpretation of this multi-messenger signature requires a detailed understanding of the instability mechanism. A perturbative analysis is used to characterize the properties of SASI, and assess the effect of the region of neutronization above the surface of the proto-neutron star. The eigenfrequencies of the most unstable modes are compared to those obtained in an adiabatic approximation where neutrino interactions are neglected above the neutrinosphere. The differential system is solved analytically using a Wronskian method and approximated asymptotically for a large shock radius. The oscillation period of SASI is well fitted with a simple analytic function of the shock radius, the radius of maximum deceleration and the mass of the proto-neutron star. The oscillation period is weakly dependent on the parametrized cooling function which however affects the SASI growth rate. The general properties of SASI eigenmodes are described using an adiabatic model. In this approximation the eigenvalue problem is formulated as a self-forced oscillator. The forcing agent is the radial advection of baroclinic vorticity perturbations and entropy perturbations produced by the shock oscillation. The differential system defining the eigenfrequencies is reduced to a single integral equation. Its analytical approximation sheds light on the radially extended character of the region of advective-acoustic coupling. The simplicity of this adiabatic formalism opens new perspectives to investigate the effect of stellar rotation and non-adiabatic processes on SASI.

During the late stages of giant planet formation, protoplanets are surrounded by a circumplanetary disk and an infalling envelope of gas and dust. For systems with sufficient cooling, material entering the sphere of influence of the planet falls inward and approaches ballistic conditions. Due to conservation of angular momentum, most of the incoming material falls onto the disk rather than directly onto the planet. This paper determines the spectral energy distributions of forming planets in this stage of evolution. Generalizing previous work, we consider a range of possible geometries for the boundary conditions of the infall and determine the two-dimensional structure of the envelope, as well as the surface density of the disk. After specifying the luminosity sources for the planet and disk, we calculate the corresponding radiative signatures for the system, including the emergent spectral energy distributions and emission maps. These results show how the observational appearance of forming planets depend on the input parameters, including the instantaneous mass, mass accretion rate, semimajor axis of the orbit, and the planetary magnetic field strength (which sets the inner boundary condition for the disk). We also consider different choices for the form of the opacity law and attenuation due to the background circumstellar disk. Although observing forming planets will be challenging, these results show how the observational signatures depend on the underlying properties of the planet/disk/envelope system.

Studying high-energy cosmic-ray air showers through the radio emission produced by their secondary particles is a well-established technique. However, due to the increasing size and density of the radio arrays, analyses are running into computational limits, as these rely on Monte Carlo simulations to model the emission. To address this, we have been developing template synthesis. With this method, we use semi-analytical expressions to describe how the radio emission from an air shower depends on the shower age and the position of the antenna with respect to the shower. These expressions are extracted from a set of microscopic simulations, thus benefiting from their accuracy. Once obtained, we can use these relations to synthesise the emission from an air shower with any longitudinal profile, by using a single Monte Carlo simulation as an input. Previously we have demonstrated that this hybrid approach can synthesise the radio emission from air showers and agrees with results from microscopic simulations within 10%. The method was however limited to a specific geometry. Here we present our first step towards generalising template synthesis across geometries. We found a set of scaling relations which correct for the shower geometry as well as the viewing angle under which the radiation is observed. This allows us to reformulate the semi-analytical relations in a way that does not longer depend on the geometry, significantly reducing the number of parameters that need to be fitted. We apply these scaling relations to a simulation library of CORSIKA showers with a zenith angle of 50 degrees. We then extract the semi-analytical expressions required for template synthesis, and use them to synthesise the emission from air showers with lower zenith angles. We investigate the accuracy by comparing both to microscopic simulations as well as the single geometry version of template synthesis.

Over the last few decades, radio detection has become one of the standard techniques to study high-energy cosmic-ray air showers. For the purpose of analysing the data, we heavily rely on Monte Carlo simulations. Upcoming dense radio array experiments such as LOFAR2.0 and SKA will, however, reach the limit of what is computationally feasible with these. Other techniques are available, based on macroscopic quantities, but their accuracy has thus far not been adequate to use them in precision analyses. In this contribution we present the latest update on the template synthesis approach, a hybrid model using both micro- and macroscopic inputs to synthesise the radio emission for an air shower with an arbitrary longitudinal profile. The method starts from the emission of a given shower and employs semi-analytical relations which only depend on the atmospheric depth at shower maximum and antenna position in order to transform it. Core to the template synthesis approach is the slicing of the atmosphere. By considering the radio emission from each slice separately, we only need to explicitly account for shower age effects. In previous work it was shown this could be done over a wide range of primary energy and across primary types for vertical air showers, with an accuracy of 10%. Here, we generalise the method to other zenith angles. We investigate the potential to synthesise between different geometries using a data set consisting of several hundreds of CORSIKA showers with primary energies ranging from $10^{17}$ eV to $10^{19}$ eV.

We present a novel way to synthesise the radio emission from extensive air showers. It is a hybrid approach which uses a single microscopic Monte-Carlo simulation to generate the radio emission from a shower with a different longitudinal evolution, primary particle type and energy. The method employs semi-analytical relations which only depend on the shower parameters to transform the radio signal in the simulated antennas. We apply this method to vertical air showers with energies ranging from $10^{17}$ eV to $10^{19}$ eV and compare the results with CoREAS using two different metrics. In order to gauge the performance over our simulation set, we subsequently use every shower in the set as a template to synthesise the emission from the other showers. Depending on the scoring metric, template synthesis reconstructs the radio emission with an accuracy of 5 to 10%.

The magnetic field strength and its topology play an important role in understanding the formation, evolution, and dynamics of the solar corona. Also, it plays a significant role in addressing long-standing mysteries such as coronal heating problem, origin and propagation of coronal mass ejections, drivers of space weather, origin and acceleration of solar wind, and so on. Despite having photospheric magnetograms for decades, we do not have reliable observations of coronal magnetic field strengths today. To measure the coronal magnetic field precisely, the spectropolarimetry channel of the Visible Emission Line Coronagraph (VELC) on board the Aditya-L1 mission is designed. Using the observations of coronal emission line Fe XIII [10747{Å~}], it is possible to generate full Stokes maps (I, Q, U, and V) that help in estimating the Line-of-Sight (LOS) magnetic field strength and to derive the magnetic field topology maps of solar corona in the Field of View (FOV) (1.05 -- 1.5~R$_{\odot}$). In this article, we summarize the instrumental details of the spectropolarimetry channel and detailed calibration procedures adopted to derive the modulation and demodulation matrices. Furthermore, we have applied the derived demodulation matrices to the observed data in the laboratory and studied their performance.

We present a comprehensive analysis of the broadband spectral energy distribution (SED) of the extreme high-energy peaked BL Lac (EHBL) source, 1ES 0229$+$200. Our study utilizes near-simultaneous data collected at various epochs between September 2017 and August 2021 (MJD: 58119$-$59365) from different instruments, including {\em AstroSat}$-$UVIT, SXT, LAXPC, {\em Swift}$-$UVOT, {\em Fermi}-LAT, and MAGIC. We investigate the one-zone synchrotron and synchrotron self-Compton (SSC) model, employing diverse particle distributions such as the log parabola, broken power law, power law with a maximum electron energy $\gamma$, energy-dependent diffusion (EDD), and energy-dependent acceleration (EDA) models to fit the broadband SED of the source. Our findings indicate that both peaks in the SED are well described by the one-zone SSC model across all particle distribution models. We estimate the jet power for different particle distributions. The estimated jet power for broken power law particle distributions is found to be on the order of $10^{47}$ ($10^{44}$) erg s$^{-1}$ for a minimum electron energy $\gamma_{min}$ $\sim$10 (10$^4$). However, for intrinsically curved particle energy distributions (e.g., log parabola, EDD, and EDA models), the estimated jet power is $\sim$10$^{44}$ erg s$^{-1}$. The SED fitting at five epochs enables us to explore the correlation between the derived spectral parameters of various particle distribution models. Notably, the observed correlations are inconsistent with the predictions in the power-law with a maximum $\gamma$ model, although the EDD and EDA models yield the correlations as expected. Moreover, the estimated physical parameter values are consistent with the model assumptions.

In this conference summary I first provide some historical notes on my own collaboration with You-Hua Chu and on the discovery of X-rays from superbubbles by Margarita Rosado S. I then consider the central subject of the conference, stellar feedback, and how interactions between stellar winds and the interstellar medium can limit or enhance the effects of feedback compared to models including only supernova explosions. Finally, I review results in other areas covered by the conference, including planet and star formation, nova and supernova remnants, different topics in stellar evolution and interaction with the interstellar medium, star clusters, observational surveys, and observational and numerical techniques.

Transitional millisecond pulsars (tMSPs) represent a dynamic category of celestial sources that establish a crucial connection between low-mass X-ray binaries and millisecond radio pulsars. These systems exhibit transitions from rotation-powered states to accretion-powered ones and vice versa, highlighting the tight evolutionary link expected by the so-called recycling scenario. In their active phase, these sources manifest two distinct emission modes named high and low, occasionally punctuated by sporadic flares. Here, we present high-time-resolution spectroscopic observations of the binary tMSP J1023+0038, in the sub-luminous disc state. This is the first short-timescale (~ 1 min) optical spectroscopic campaign ever conducted on a tMSP. The campaign was carried out over the night of June 10, 2021 using the Gran Telescopio Canarias. The optical continuum shows erratic variability, without clear evidence of high and low modes or of orbital modulation. Besides, the analysis of these high-temporal-cadence spectroscopic observations reveals, for the first time, evidence for a significant (up to a factor of ~ 2) variability in the emission line properties (equivalent width and full width half maximum) over a timescale of minutes. Intriguingly, the variability episodes observed in the optical continuum and in the emission line properties seem uncorrelated, making their origin unclear.

Phosphorus (P) is an important element for the chemical evolution of galaxies and many kinds of biochemical reactions. Phosphorus is one of the crucial chemical compounds in the formation of life on our planet. In an interstellar medium, phosphine (PH$_{3}$) is a crucial biomolecule that plays a major role in understanding the chemistry of phosphorus-bearing molecules, particularly phosphorus nitride (PN) and phosphorus monoxide (PO), in the gas phase or interstellar grains. We present the first confirmed detection of phosphine (PH$_{3}$) in the asymptotic giant branch (AGB) carbon-rich star IRC+10216 using the Atacama Large Millimeter/Submillimeter Array (ALMA) band 6. We detect the $J$ = 1$_{0}$$-$0$_{0}$ rotational transition line of PH$_{3}$ with a signal-to-noise ratio (SNR) of $\geq$3.5$\sigma$. This is the first confirmed detection of phosphine (PH$_{3}$) in the ISM. Based on LTE spectral modelling, the column density of PH$_{3}$ is (3.15$\pm$0.20)$\times$10$^{15}$ cm$^{-2}$ at an excitation temperature of 52$\pm$5 K. The fractional abundance of PH$_{3}$ with respect to H$_{2}$ is (8.29$\pm$1.37)$\times$10$^{-8}$. We also discuss the possible formation pathways of PH$_{3}$ and we claim that PH$_{3}$ may be created via the hydrogenation of PH$_{2}$ on the grain surface of IRC+10216.

We present the final results of the Apodized Pupil Lyot Coronagraph (APLC) on the High-contrast imager for Complex Aperture Telescopes (HiCAT) testbed, under NASA's Strategic Astrophysics Technology program. The HiCAT testbed was developed over the past decade to enable a system-level demonstration of coronagraphy for exoplanet direct imaging with the future Habitable Wolds Observatory. HiCAT includes an active, segmented telescope simulator, a coronagraph, and metrology systems (Low-order and Mid-Order Zernike Wavefront Sensors, and Phase Retrieval camera). These results correspond to an off-axis (un-obscured) configuration, as was envisioned in the 2020 Decadal Survey Recommendations. Narrowband and broadband dark holes are generated using two continuous deformable mirrors (DM) to control high order wavefront aberrations, and low-order drifts can be further stabilized using the LOWFS loop. The APLC apodizers, manufactured using carbon nanotubes, were optimized for broadband performance and include the calibrated geometric aperture. HiCAT is, to this date, the only testbed facility able to demonstrate high-contrast coronagraphy with a truly segmented aperture, as is required for the Habitable World Observatory, albeit limited to ambient conditions. Results presented here include $6\times 10^{-8}$ (90% CI) contrast in 9% bandpass in a 360 deg dark hole with inner and outer working angles of $4.4 \lambda/D_{pupil}$ and $11 \lambda/D_{pupil}$ . Narrowband contrast (3% bandpass) reaches $2.4\times 10^{-8}$ (90% confidence interval).

To improve the performance of full-shape analyses of large-scale structure, we consider using a halo occupation distribution (HOD)-informed prior for the effective field theory (EFT) nuisance parameters. We generate 320 000 mock galaxy catalogs using 10 000 sets of HOD parameters across 32 simulation boxes with different cosmologies. We measure and fit the redshift-space power spectra using a fast emulator of the EFT model, and the resulting best-fit EFT parameter distributions are used to create the prior. This prior effectively constrains the EFT nuisance parameter space, limiting it to the space of HOD-mocks that can be well fit by a EFT model. We have tested the stability of the prior under different configurations, including the effect of varying the HOD sample distribution and the inclusion of the hexadecapole moment. We find that our HOD-informed prior and the cosmological parameter constraints derived using it are robust. While cosmological fits using the standard EFT prior suffer from prior effects, sometimes failing to recover the true cosmology within Bayesian credible intervals, the HOD-informed prior mitigates these issues and significantly improves cosmological parameter recovery for $\Lambda$CDM and beyond. This work lays the foundation for better full-shape large-scale structure analyses in current and upcoming galaxy surveys, making it a valuable tool for addressing key questions in cosmology.

The cosmic microwave background (CMB) spectrum is an extraordinary tool for exploring physics beyond the Standard Model. The exquisite precision of the measurement makes it particularly sensitive to small effects caused by hidden sector interactions. In particular, CMB spectral distortions can unveil the existence of dark photons which are kinetically coupled to the standard photon. In this work, we use the COBE-FIRAS dataset to derive accurate and robust limits on photon-to-dark-photon oscillations for a large range of dark photon masses, from $10^{-10}$ to $10^{-4}$ eV. We consider in detail the redshift dependence of the bounds, computing CMB distortions due to photon injection/removal using a Green's function method. Our treatment improves on previous results, which had set limits studying energy injection/removal into baryons rather than photon injection/removal, or ignored the redshift evolution of distortions. The difference between our treatment and previous ones is particularly noticeable in the predicted spectral shape of the distortions, a smoking gun signature for photon-to-dark-photon oscillations. The characterization of the spectral shape is crucial for future CMB missions, which could improve the present sensitivity by orders of magnitude, exploring regions of the dark photon parameter space that are otherwise difficult to access.

We implement a simple, main beam correction in the maximum-likelihood, parametric component separation approach, which allows on accounting for different beamwidths of input maps at different frequencies without any preprocessing. We validate the approach on full-sky and cut-sky simulations and discuss the importance and impact of the assumptions and simplifications. We find that, in the cases when the underlying sky model is indeed parametric, the method successfully recovers component spectral parameters and component maps at the pre-defined resolution. The improvement on the precision of the estimated spectral parameters is found to be minor due to the redness of the foreground angular spectra, however the method is potentially more accurate, in particular if the foreground properties display strong, spatial variability, as it does not assume commutation of the beam smoothing and mixing matrix operators. The method permits a reconstruction of the CMB map with a resolution significantly superior to that of the lowest resolution map used in the analysis and with the nearly optimal noise level, facilitating exploitation of the cosmological information contained on angular scales, which would be otherwise inaccessible. The method preserves all the advantages of a pixel-domain implementation of the parametric approach, and, as it deals with the beams in the harmonic domain, it can also straightforwardly account for spatially stationary map-domain noise correlations.

We discuss a mechanism of primordial black hole (PBH) formation that does not require specific features in the inflationary potential, revisiting previous literature. In this mechanism, a light spectator field evolves stochastically during inflation and remains subdominant during the post-inflationary era. Even though the curvature power spectrum stays small at all scales, rare perturbations of the field probe a local maximum in its potential, leading to non-Gaussian tails in the distribution of curvature fluctuations, and to copious PBH production. For a concrete axion-like particle (ALP) scenario we analytically determine the distribution of the compaction function for perturbations, showing that it is characterized by a heavy tail, which produces an extended PBH mass distribution. We find the ALP mass and decay constant to be correlated with the PBH mass, for instance, an ALP with a mass $m_a = 5.4 \times 10^{14}$ eV and a decay constant $f_a = 4.6 \times 10^{-5} Mpl$ can lead to PBHs of mass $M_{\rm PBH} = 10^{21}$ g as the entire dark matter (DM) of the universe, and is testable in future PBH observations via lensing in the NGRST and mergers detectable in the LISA and ET Gravitational Waves (GW) detectors. We then extend our analysis to mixed ALP and PBH dark matter and Higgs-like spectator fields. We find that PBHs cluster strongly over all cosmological scales, clashing with CMB isocurvature bounds. We argue that this problem is shared by all PBH production from inflationary models that depend solely on large non-Gaussianity without a peak in the curvature power spectrum and discuss possible remedies.