Papers for Thursday, Jun 13 2024

The recent IceCube detection of TeV neutrino emission from the nearby active galaxy NGC 1068 suggests that active galactic nuclei (AGN) could make a sizable contribution to the diffuse flux of astrophysical neutrinos. The absence of TeV $\gamma$-rays from NGC 1068 indicates neutrino production in the vicinity of the supermassive black hole, where the high radiation density leads to $\gamma$-ray attenuation. Therefore, any potential neutrino emission from similar sources is not expected to correlate with high-energy $\gamma$-rays. Disk-corona models predict neutrino emission from Seyfert galaxies to correlate with keV X-rays, as they are tracers of coronal activity. Using through-going track events from the Northern Sky recorded by IceCube between 2011 and 2021, we report results from a search for individual and aggregated neutrino signals from 27 additional Seyfert galaxies that are contained in the BAT AGN Spectroscopic Survey (BASS). Besides the generic single power-law, we evaluate the spectra predicted by the disk-corona model. Assuming all sources to be intrinsically similar to NGC 1068, our findings constrain the collective neutrino emission from X-ray bright Seyfert galaxies in the Northern Hemisphere, but, at the same time, show excesses of neutrinos that could be associated with the objects NGC 4151 and CGCG 420-015. These excesses result in a 2.7$\sigma$ significance with respect to background expectations.

Gaia BH3 is the most massive known stellar-origin black hole, with a mass $M_\bullet \approx 33 M_\odot$. Detected from Gaia's astrometry, this black hole resembles those observed via gravitational waves, whose nature is still highly debated. Hosted in a binary system with a companion giant star that is too far away for Roche-lobe mass transfer, this black hole could nonetheless accrete at low levels due to wind-driven mass-loss from its companion star, thus accreting in advection-dominated accretion flow, or ADAF, mode. Using stellar models, we constrain its Eddington ratio in the range $10^{-9} < f_{\rm Edd} < 10^{-7}$, corresponding to radiative efficiencies $5\times10^{-5} < \epsilon < 10^{-3}$, compatible with radiatively inefficient accretion modes. Chandra ACIS-S observed this object and obtained the most sensitive upper bound of its [2-10] keV flux: $F_X < 3.25 \times 10^{-15} \, \rm erg \, s^{-1} \, cm^{-2}$ at $90\%$ confidence level corresponding to $L_{[2-10]}< 2.10 \times 10^{29} \, \rm erg \, s^{-1}$. Using ADAF emission models, we constrained its accretion rate to $f_{\rm Edd}< 4.91 \times 10^{-7}$ at the apastron, in agreement with our theoretical estimate. At the periastron, we expect fluxes $\sim 50$ times larger. Hence, accretion did not significantly contribute to black hole growth over the system's lifetime. Detecting the electromagnetic emission from Gaia BH3 will be fundamental to informing stellar wind and accretion disk models.

In a hierarchical, dark matter-dominated Universe, stellar mass functions (SMFs), galaxy merger rates, star formation histories (SFHs), satellite abundances, and intracluster light, should all be intimately connected observables. However, the systematics affecting observations still prevent universal and uniform measurements of, for example, the SMF and the SFHs, inevitably preventing theoretical models to compare with multiple data sets robustly and simultaneously. We here present our holistic semi-empirical model DECODE (Discrete statistical sEmi-empiriCal mODEl) that converts via abundance matching dark matter merger trees into galaxy assembly histories, using different SMFs in input and predicting all other observables in output in a fully data-driven and self-consistent fashion with minimal assumptions. We find that: 1) weakly evolving or nearly constant SMFs below the knee ($M_\star \lesssim 10^{11} \, M_\odot$) are the best suited to generate star formation histories aligned with those inferred from MaNGA, SDSS, GAMA, and, more recently, JWST; 2) the evolution of satellites after infall only affects the satellite abundances and star formation histories of massive central galaxies but not their merger histories; 3) the resulting SFR-$M_\star$ relation is lower in normalization by a factor of $\sim 2$ with respect to observations, with a flattening at high masses more pronounced in the presence of mergers; 4) the latest data on intracluster light can be reproduced if mass loss from mergers is included in the models. Our findings are pivotal in acting as pathfinder to test the self-consistency of the high-quality data from, e.g., JWST and Euclid.

We present optical and near-infrared (NIR) spectroscopic observations for a sample of $45$ quasars at $6.50 < z \leq 7.64$ with absolute magnitudes at $1450$ Å in the range $-28.82 \leq M_{1450} \leq -24.13$ and their composite spectrum. The median redshift and $M_{1450}$ of the quasars in the sample are $z_{\rm{median}}=6.71$ and $M_{1450,\rm{median}} \simeq -26.1$, respectively. The NIR spectra are taken with echelle spectrographs, complemented with additional data from optical long slit instruments, and then reduced consistently using the open-source Python-based spectroscopic data reduction pipeline PypeIt. The median value of the mean signal-to-noise ratios of the spectra in J, H, and K band (median $\langle \rm{SNR}_{\lambda} \rangle$) is: median $\langle \rm{SNR}_{J} \rangle=9.7$, median $\langle \rm{SNR}_{H} \rangle=10.3$, and median $\langle \rm{SNR}_{K} \rangle=11.7$; demonstrating the good data quality. This work presents the largest medium/moderate-resolution sample of quasars at $z>6.5$ from ground-based instruments. Its homogeneity and reproducibility make it ideally suited for several scientific goals, i.e., the study of the quasar proximity zones and damping wings, the Ly$\alpha$ forest, the intergalactic medium's metal content, as well as other properties such as the distribution of SMBH masses and Eddington ratios. Our composite spectrum is compared to others at both high and low-$z$ from the literature, showing differences in the strengths of many emission lines, probably due to differences in luminosity among the samples, but a consistent continuum slope, which proves that the same spectral features are preserved in quasars at different redshift ranges.

The escape fraction of LyC ionizing radiation $f_{LyC}^{esc}$ is crucial for understanding reionization, yet impossible to measure at z $\gtrsim$ 5.3. Recently, studies have focused on calibrating indirect indicators of $f_{LyC}^{esc}$ at z $\sim$ 0.3, finding that Ly$\alpha$ is closely linked to it. What is still unclear is whether the LyC - Ly$\alpha$ relation evolves with redshift, and if Ly$\alpha$ is truly applicable as an $f_{LyC}^{esc}$ indicator during the reionization epoch. In this study, we investigate seven $-21 \lesssim M_{UV} \lesssim -19$ gravitationally lensed galaxies from the BELLS GALLERY survey at z $\sim$ 2.3. Our targets have rest-frame Ly$\alpha$ equivalent widths ranging from 10 Å to 100 Å and low dust content ($-2.5 \lesssim \beta \lesssim -1.9$), both indicative of high LyC escape. Surprisingly, direct estimates of $f_{LyC}^{esc}$ using Hubble Space Telescope imaging with F275W and F225W reveal that our targets are not LyC emitters, with absolute $f_{LyC}^{esc}$ $\lesssim$ 6.5% at 3$\sigma$ significance (with two sources having absolute $f_{LyC}^{esc}$(3$\sigma$) $\lesssim$ 10% and $\lesssim$ 16%). The low $f_{LyC}^{esc}$, coupled with the high Ly$\alpha$ escape fraction and equivalent width could potentially be attributed to the redshift evolution of the neutral hydrogen column density and dust content. Our results challenge previous studies based on local samples, suggesting that the extrapolation of z ~ 0 Ly$\alpha$-based LyC indirect estimators into the reionization epoch might not be correct.

High resolution (0."037-0."13 [10-35 pc]) e-MERLIN ($\lambda6-18$ cm) and (0."024 [6.5 pc]) ALMA ($\lambda 1.1$ mm) observations have been used to image OH (hydroxyl) and H$_2$CO (formaldehyde) megamaser emission, and HCN 3->2 emission towards the nuclear (<100 pc) region of the luminous infrared galaxy Zw049.057. Zw049.057 hosts a compact obscured nucleus (CON), thus representing a class of galaxies that are often associated with inflow and outflow motions. Formaldehyde megamaser emission is detected towards the nuclear region, <30 pc (<0."1), and traces a structure along the disk major axis. OH megamaser (OHM) emission is detected along the minor axis of the disk, ~30 pc (0."1) from the nucleus, where it exhibits a velocity gradient with extrema of -20 km/s south-east (SE) of the disk and -110 km/s north-west (NW) of the disk. HCN 3->2 emission reveals extended emission, along the disk minor axis out to ~60 pc (0."2). Analysis of the minor axis HCN emission reveals high-velocity features, extending out to 600 km/s, redshifted on the SE side and blueshifted on the NW side. We propose that the high-velocity HCN emission traces a fast >250 km/s and collimated outflow, that is enveloped by a wide-angle and slow ~50 km/s outflow that is traced by the OHM emission. Analysis of the outflow kinematics suggests that the slow wide-angle outflow will not reach escape velocity and instead will fall back to the galaxy disk, evolving as a so-called fountain flow, while the fast collimated outflow traced by HCN emission will likely escape the nuclear region. We suggest that the absence of OHM emission in the nuclear region is due to high densities there. Even though OHMs associated with outflows are an exception to conventional OHM emission, we expect them to be common in CON sources that host both OHM and H$_2$CO megamasers.

We investigate the multi-phase structure of gas flows in galaxies. We study 80 galaxies during the epoch of peak star formation ($1.4\leq z\leq2.7$) using data from Keck/LRIS and VLT/KMOS. Our analysis provides a simultaneous probe of outflows using UV emission and absorption features and H$\alpha$ emission. With this unprecedented data set, we examine the properties of gas flows estimated from LRIS and KMOS in relation to other galaxy properties, such as star formation rate (SFR), star formation rate surface density ($\Sigma_{\rm SFR}$), stellar mass (M$_*$), and main sequence offset ($\Delta$MS). We find no strong correlations between outflow velocity measured from rest-UV lines and galaxy properties. However, we find that galaxies with detected outflows show higher averages in SFR, $\Sigma_{\rm SFR}$, and $\Delta$MS than those lacking outflow detections, indicating a connection between outflow and galaxy properties. Furthermore, we find a lower average outflow velocity than previously reported, suggesting greater absorption at the systemic redshift of the galaxy. Finally, we detect outflows in 49% of our LRIS sample and 30% in the KMOS sample, and find no significant correlation between outflow detection and inclination. These results may indicate that outflows are not collimated and that H$\alpha$ outflows have a lower covering fraction than low-ionization interstellar absorption lines. Additionally, these tracers may be sensitive to different physical scales of outflow activity. A larger sample size with a wider dynamic range in galaxy properties is needed to further test this picture.

Recent observations of three nearby black hole low-mass X-ray binaries have suggested possible evidence of dark matter density spikes. It has long been established that cold dark matter should form dense spikes around intermediate-mass and supermassive black holes due to adiabatic compression of the surrounding dark matter halo. Lighter black holes formed from stellar collapse, however, are not expected to develop such spikes, and so it is unclear how the stellar-mass black holes in these binaries could have acquired such features. Given that primordial black holes are expected to form ultra-dense dark matter mini-spikes, in this Letter we explore the possibility that these stellar-mass black holes may actually have a primordial origin.

We study the dynamics of nuclear star clusters, the dense stellar environments surrounding massive black holes in the centers of galaxies. We consider angular momentum diffusion due to two-body scatterings among stellar objects and energy advection due to gravitational wave emission upon interaction with the central massive black hole. Such dynamics is described by a two-dimensional Fokker-Planck equation in energy-angular momentum space. Focusing on the transition between the diffusion-dominated region and the advection-dominated one, we utilize self-similarity to obtain a full solution for the Fokker-Planck equation. This solution provides the density and flux of the stellar objects in nuclear star clusters. This improves the rate estimates for extreme mass-ratio inspirals, and has interesting implications for a new class of galactic center transients called quasi-periodic eruptions.

This paper presents a new framework for understanding the relationship between a galaxy and its circumgalactic medium (CGM). It focuses on how imbalances between heating and cooling cause either expansion or contraction of the CGM. It does this by tracking \textit{all} of the mass and energy associated with a halo's baryons, including their gravitational potential energy, even if feedback has pushed some of those baryons beyond the halo's virial radius. We show how a star-forming galaxy's equilibrium state can be algebraically derived within the context of this framework, and we analyze how the equilibrium star formation rate depends on supernova feedback. We consider the consequences of varying the mass loading parameter etaM = Mdot_wind / Mdot_* relating a galaxy's gas mass outflow rate (Mdot_wind) to its star formation rate (Mdot_*) and obtain results that challenge common assumptions. In particular, we find that equilibrium star formation rates in low-mass galaxies are generally insensitive to mass loading, and when mass loading does matter, increasing it actually results in \textit{more} star formation because more supernova energy is needed to resist atmospheric contraction.

The scaling of galaxy properties with halo mass suggests that feedback loops regulate star formation, but there is no consensus yet about how those feedback loops work. To help clarify discussions of galaxy-scale feedback, Paper I presented a very simple model for supernova feedback that it called the minimalist regulator model. This followup paper interprets that model and discusses its implications. The model itself is an accounting system that tracks all of the mass and energy associated with a halo's circumgalactic baryons--the central galaxy's atmosphere. Algebraic solutions for the equilibrium states of that model reveal that star formation in low-mass halos self-regulates primarily by expanding the atmospheres of those halos, ultimately resulting in stellar masses that are insensitive to the mass-loading properties of galactic winds. What matters most is the proportion of supernova energy that couples with circumgalactic gas. However, supernova feedback alone fails to expand galactic atmospheres in higher-mass halos. According to the minimalist regulator model, an atmospheric contraction crisis ensues, which may be what triggers strong black-hole feedback. The model also predicts that circumgalactic medium properties emerging from cosmological simulations should depend largely on the specific energy of the outflows they produce, and we interpret the qualitative properties of several numerical simulations in light of that prediction.

We address the formation of giant clumps in violently unstable gas-rich disc galaxies at cosmic noon. While these are commonly thought to originate from gravitational Toomre instability, cosmological simulations have indicated that clumps form even in regions where the Toomre $Q$ parameter is well above unity, which should be stable according to linear Toomre theory (Inoue et al., 2016). Examining one of these cosmological simulations, we find that it exhibits an excess in compressive modes of turbulence with converging motions. The energy in converging motions within proto-clump regions is $\sim 70\%$ of the total turbulent energy, compared to $\sim 17\%$ expected in equipartition. When averaged over the whole disc, $\sim 32\%$ of the turbulent energy is in converging motions, with a further $\sim 8\%$ in diverging motions. Thus, a total of $\sim 40\%$ of the turbulent energy is in compressive modes, with the rest in solenoidal modes, compared to the $(1/3)-(2/3)$ division expected in equipartition. By contrast, we find that in an isolated-disc simulation with similar properties, resembling high-$z$ star-forming galaxies, the energy in the different turbulence modes are in equipartition, both in proto-clump regions and over the whole disc. We conclude that the origin of the excessive converging motions in proto-clump regions is external to the disc, and propose several mechanisms that can induce them. This is an additional mechanism for clump formation, complementary to and possibly preceding gravitational instability.

Blazars can be detected from very large distances due to their high luminosity. However, the detection of $\gamma$-ray emission of blazars beyond $z=3$ has only been confirmed for a small number of sources. Such observations probe the growth of supermassive black holes close to the peak of star formation in the history of galaxy evolution. As a result from a continuous monitoring of a sample of 80 $z>3$ blazars with Fermi-LAT, we present the first detection of a $\gamma$-ray flare from the $z=4.31$ blazar TXS 1508+572. This source showed high $\gamma$-ray activity from February to August 2022, reaching a peak luminosity comparable to the most luminous flares ever detected with Fermi -LAT. We conducted a multiwavelength observing campaign involving XMM-Newton, Swift, the Effelsberg 100-m radio telescope and the Very Long Baseline Array. In addition, we make use of the monitoring programs by the Zwicky Transient Facility and NEOWISE at optical and infrared wavelengths, respectively. We find that the source is particularly variable in the infrared band on daily time scales. The spectral energy distribution collected during our campaign is well described by a one-zone leptonic model, with the $\gamma$-ray flare originating from an increase of external Compton emission as a result of a fresh injection of accelerated electrons.

Studying the structures of open clusters is crucial for understanding stellar evolution and galactic dynamics. Based on Gaia DR3 data, we apply the hierarchical clustering algorithm to a young open cluster NGC 6530 and group its members into 5 substructures. By linear tracing with the kinematic information of their members, we find that: Sub 1 is the core of the cluster. It is expanding slowly. Sub 2 consists of less bound members, which began escaping from the core about 0.78 Myr ago. Sub 3 is associated with a young star forming region. It will merge with the core after 0.72 Myr; Sub 4, as an outskirt group, is also moving towards the core, but won't end up falling in. While Sub 5 is composed of less-bound members with field contamination. This work reveals the complex internal structure and evolutionary trends of the cluster NGC 6530. It also shows the potential of the hierarchical clustering algorithm in star cluster structure analysis.

Recent radio surveys have revealed pulsars with dispersion and scattering delays induced by ionized gas that are larger than the rest of the observed pulsar population, in some cases with electron column densities (or dispersion measures; DMs) larger than the maximum predictions of Galactic electron density models. By cross-matching the observed pulsar population against HII region catalogs, we show that the majority of pulsars with $\rm DM > 600$ pc cm$^{-3}$ and scattering delays $\tau(1\ {\rm GHz}) > 10$ ms lie behind HII regions, and that HII region intersections may be relevant to as much as a third of the observed pulsar population. Accounting for HII regions resolves apparent discrepancies where Galactic electron density models place high-DM pulsars beyond the Galactic disk. By comparing emission measures (EMs) inferred from recombination line observations to pulsar DMs, we show that HII regions can contribute tens to hundreds of pc cm$^{-3}$ in electron column density along a pulsar LOS. We find that nearly all pulsars with significant excess (and deficit) scattering from the mean $\tau$-DM relation are spatially coincident with known discrete ionized gas structures, including HII regions. Accounting for HII regions is critical to the interpretation of radio dispersion and scattering measurements as electron density tracers, both in the Milky Way and in other galaxies.

The amount of turbulence in protoplanetary discs around young stars is critical for determining the efficiency, timeline, and outcomes of planet formation. It is also difficult to measure. Observations are still limited, but direct measurements of the non-thermal, turbulent gas motion are possible with the Atacama Large Millimeter/submillimeter Array (ALMA). Using CO(2-1)/$^{13}$CO(2-1)/C$^{18}$O(2-1) ALMA observations of the disc around IM Lup at ~0.4" (~60 au) resolution we find evidence of significant turbulence, at the level of $\delta v_{\rm turb}=(0.18-0.30)$c$_s$. This result is robust against systematic uncertainties (e.g., amplitude flux calibration, midplane gas temperature, disc self-gravity). We find that gravito-turbulence as the source of the gas motion is unlikely based on the lack of an imprint on the rotation curve from a massive disc, while magneto-rotational instabilities and hydrodynamic instabilities are still possible, depending on the unknown magnetic field strength and the cooling timescale in the outer disc.

For a particular generalised-NFW (gNFW) mass distribution, we derive analytical expressions for its surface mass density, projected mass, and deflection potential, given in terms of three representations of the Fox $H$-function, for which we provide their corresponding power (logarithmic) series expansions. They are handy in computationally intensive tasks. From these results, we obtain closed-form expressions for the super-NFW (sNFW) in terms of complete elliptic functions, which have not yet been reported in the literature. Additionally, we find that, for a fixed $\kappa_0$ (characteristic convergence), when the number of images is maximal (three for the sNFW), the sum of their signed magnification, named $I_{\text{sNFW}}$, varies with the position of the source $y$ inside the radial caustic $y_{\text{c2}}$. The higher variations occur for $\kappa_0\lesssim 1$, while $I_{\text{sNFW}}$ can be considered as constant when $\kappa_0$ is approximately within the range 2-10. This range can be extended depending on the observational uncertainties, since for higher values of $\kappa_0$, the variations in $I_{\text{sNFW}}$ are relatively small. This behaviour is shared with the NFW and Hernquist lenses.

We present the idea that replacing the cosmological constant $\Lambda$ in the $\Lambda$CDM model by a distribution of walls, with very low tension compared to what one would expect from new physics, could help explaining the tension in the Hubble constant fits in the Standard Cosmological Model. Using parameters from our model for dark matter as macroscopic pearls, we can get a promising order of magnitude for the correction to the Hubble constant estimated from observations of the cosmic microwave background. Our model is on the borderline of failing by predicting too large extra fluctuations as a function of direction in the cosmological microwave background radiation. However, imagining the bubbles in the voids to have come from more somewhat smaller big bubbles also occurring outside the big voids may help. We estimate that, in order to have big volumes of the new vacuum in intergalactic space, a very high temperature is needed and that such regions would be likely to get cooled, freeze and shrink down to the degenerate form of dark matter if hitting some ordinary matter, as is likely in the denser parts of the Universe. We also review our model for dark matter, and develop the understanding of the stopping of the dark matter particles in the shielding of say the DAMA-LIBRA underground experiment and the counting rate this experiment observes. We manage to obtain a consistent fit with a mass $M = 2 \cdot 10^{-18}$kg = $10^9$GeV and radius $R = 10^{-10}$m for the dark matter pearls, corresponding to a tension in the domain wall of $S= (8 {\rm MeV})^3$.

Survey strategies for upcoming exoplanet direct imaging missions have considered varying assumptions of prior knowledge. Precursor radial velocity surveys could have detected nearby exo-Earths and provided prior orbit and mass constraints. Alternatively, a direct imaging mission performing astrometry could yield constraints on orbit and phase angle of target planets. Understanding the impact of prior mass and orbit information on planetary characterization is crucial for efficiently recognizing habitable exoplanets. To address this question, we use a reflected-light retrieval tool to infer the atmospheric and bulk properties of directly imaged Earth-analogs while considering varying levels of prior information and signal-to-noise ratio (SNR). Because of the strong correlation between the orbit-related parameters and the planetary radius, prior information on the orbital distance and planetary phase yield tight constraints on the planetary radius: from $R_{\rm{p}}=2.95^{+2.69}_{-1.95}~R_{\oplus}$ without prior knowledge, to $R_{\rm{p}}=1.01^{+0.33}_{-0.19}~R_{\oplus}$ with prior determination of the orbit for $\rm{SNR}=20$ in the visible/near-infrared spectral range, thus allowing size determination from reflected light observations. However, additional knowledge of planet mass does not notably enhance radius ($R_{\rm{p}}=0.98^{+0.17}_{-0.14}~R_{\oplus}$) or atmospheric characterization. Also, prior knowledge of the mass alone does not yield a tight radius constraint ($R_{\rm{p}}=1.64^{+1.29}_{-0.80}~R_{\oplus}$) nor improves atmospheric composition inference. By contrast, because of its sensitivity to gas column abundance, detecting a Rayleigh scattering slope or bounding Rayleigh opacity helps to refine gas mixing ratio inferences without requiring prior mass knowledge. Overall, apart from radius determination, increasing the SNR is more beneficial than additional prior observations.

Rocky planets in compact configurations are the most common ones around M dwarfs. Many disks around very low mass stars (between 0.1 and 0.5 M$_\odot$) are rather compact and small (without observable substructures and radius less than 20 au), which favours the idea of an efficient radial drift that could enhance planet formation in compact orbits. We aim to investigate the potential formation paths of the observed close-in rocky exoplanet population around M dwarfs, assuming that planet formation could take place in compact disks with an efficient dust radial drift. We developed N-body simulations that include a sample of embryos growing by pebble accretion exposed to planet-disk interactions, star-planet tidal interactions and general relativistic corrections. For a star of 0.1 M$_\odot$ we considered different gas disk viscosity and initial embryo distributions. We also explore planet formation by pebble accretion around stars of 0.3 and 0.5 M$_\odot$. Lastly, for each stellar mass, we run simulations that include a sample of embryos growing by planetesimal accretion. Our main result is that the sample of simulated planets that grow by pebble accretion in a gas disk with low viscosity ($\alpha=10^{-4}$) can reproduce the close-in low-mass exoplanet population around M dwarfs in terms of multiplicity, masses and semi-major axis. Furthermore, we found that a gas disk with high viscosity ($\alpha=10^{-3}$) can not reproduce the observed planet masses. Also, we show that planetesimal accretion favours the formation of smaller planets than the ones formed by pebble accretion. Rocky planet formation around M dwarfs can take place in compact and small dust disks driven by an efficient radial drift in a gas disk with low viscosity. This result points towards a new approach in the direction of the disk conditions needed for rocky planet formation around very low mass stars.

Exoplanet direct imaging using adaptive optics (AO) is often limited by non-common path aberrations (NCPAs) and aberrations that are invisible to traditional pupil-plane wavefront sensors (WFSs). This can be remedied by focal-plane (FP) WFSs that characterize aberrations directly from a final science image. Photonic lanterns (PLs) can act as low-order FPWFSs with the ability to direct some light to downstream science instruments. Using a PL on the SEAL (Santa Cruz Extreme AO Laboratory) high-contrast imaging testbed, we demonstrate (1) linear ranges and (2) closed-loop control. Additionally, we simulate the use of the PL in a multi-wavefront sensor AO system, in which multiple WFSs feed back to the same common-path deformable mirror. Building on previous multi-WFS AO demonstrations on SEAL, we simulate a modulated pyramid WFS to sense aberrations of high spatial order and large amplitude, and the PL to sense low order aberrations including NCPAs. We assess adaptive optics performance in this setting using three different PL wavefront reconstruction algorithms. We also provide a new method to experimentally identify the propagation matrix of a PL, making advanced model-based algorithms practical. This work demonstrates the role of photonic technologies and multi-stage wavefront sensing in the context of extreme AO and high contrast imaging.

We present novel insights into the interplay between tidal forces and star formation in interacting galaxies before their first pericentre passage. We investigate seven close pair galaxies devoid of visible tidal disturbances, such as tails, bridges, and shells. Using integral field spectroscopy (IFS) data of extended Calar Alto Legacy Integral Field Area (eCALIFA), we unveil a previously unreported phenomenon: H alhpa emission, a proxy for recent star formation, exhibits a significant enhancement in regions facing the companion galaxy, reaching up to 1.9 times higher flux compared to opposite directions. Notably, fainter companions within pairs display a more pronounced one-sided H alpha excess, exceeding the typical range observed in isolated galaxies with 2 sigma confidence level. Furthermore, the observed H alpha excess in fainter companion galaxies exhibits a heightened prominence at the outer galactic regions. These findings suggest that tidal forces generated before the first pericentre passage exert a stronger influence on fainter galaxies due to their shallower potential wells by their brighter companions. This unveils a more intricate interplay between gravitational interactions and star formation history within interacting galaxies than previously understood, highlighting the need further to explore the early stages of interaction in galaxy evolution.

Detecting large-scale flux ropes (FRs) embedded in interplanetary coronal mass ejections (ICMEs) and assessing their geoeffectiveness are essential since they can drive severe space weather. At 1 au, these FRs have an average duration of 1 day. Their most common magnetic features are large, smoothly rotating magnetic fields. Their manual detection has become a relatively common practice over decades, although visual detection can be time-consuming and subject to observer bias. Our study proposes a pipeline that utilizes two supervised binary-classification machine learning (ML) models trained with solar wind magnetic properties to automatically detect large-scale FRs and additionally determine their geoeffectiveness. The first model is used to generate a list of auto-detected FRs. Using the properties of southward magnetic field the second model determines the geoeffectiveness of FRs. Our method identifies 88.6\% and 80\% large-scale ICMEs (duration $\ge 1$ day) observed at 1 au by Wind and Sun Earth Connection Coronal and Heliospheric Investigation (STEREO) mission, respectively. While testing with a continuous solar wind data obtained from Wind, our pipeline detected 56 of the 64 large-scale ICMEs during 2008 - 2014 period (recall= 0.875) but many false positives (precision= 0.56) as we do not take into account any additional solar wind properties than the magnetic properties. We found an accuracy of 0.88 when estimating the geoeffectiveness of the auto-detected FRs using our method. Thus, in space weather now-casting and forecasting at L1 or any planetary missions, our pipeline can be utilized to offer a first-order detection of large-scale FRs and geoeffectiveness.

We present early-time, hour-to-day cadence spectroscopy of the nearby type II supernova (SN II) 2024ggi, which was discovered at a phase when the SN shock just emerged from the red-supergiant (RSG) progenitor star. Over the first few days after the first light, SN\,2024ggi exhibited prominent narrow emission lines formed through intense and persistent photoionization of the nearby circumstellar material (CSM). In the first 63 hours, spectral lines of He, C, N, and O revealed a rapid rise in ionization, as a result of the progressive sweeping-up of the CSM by the shock. The duration of the IIn-like spectra indicates a dense and relatively confined CSM distribution extending up to $\sim\,4\,\times\,10^{14}$\,cm. Spectral modeling reveals a CSM mass loss rate at this region exceeding $5 \times\, 10^{-3}$\,\Msun\,yr$^{-1}$\,is required to reproduce low-ionization emissions, which dramatically exceeds that of an RSG. Analyzing H$\alpha$ emission shift implies the velocity of the unshocked outer CSM to be between 20 and 40 \kms, matching the typical wind velocity of an RSG. The differences between the inner and outer layers of the CSM and an RSG progenitor highlight a complex mass loss history before the explosion of SN 2024ggi.

Line intensity mapping (LIM) is a novel observational technique in astrophysics that utilizes the integrated emission from multiple atomic and molecular transition lines from galaxies to probe the complex physics of galaxy formation and evolution, as well as the large-scale structure of the universe. Modeling multiple line luminosities of galaxies with varying masses or their host halo masses poses significant uncertainty due to the lack of observational data across a wide redshift range and the intricate nature of astrophysical processes, making them challenging to model analytically or in simulations. While future experiments aim to measure multiple line intensities up to $z\sim 8$ across a wide volume using tomographic methods, we leverage publicly available datasets from the CMB experiment Planck and the galaxy survey eBOSS to constrain the CO(3-2) emission from galaxies. We correlate galaxies from eBOSS data onto the full-sky CO(2-1) map produced by Planck and report the first measurement of the average CO(3-2) intensity, $I_{CO} = 45.7 \pm 14.2\, \mathrm{Jy/sr}$ at $z\sim 0.5$ with $3.2\sigma$ confidence. Our findings demonstrate that stacking methods are already viable with existing observations from CMB experiments and galaxy surveys, and are complementary to traditional LIM experiments.

We investigate the optical variability of low-redshift ($0.15< z\leq0.4$) active galactic nuclei using the multi-epoch data from the Zwicky Transient Facility. We find that a damped random walk model well describes the ensemble structure function in the $g$ band. Consistent with previous studies, more luminous active galactic nuclei tend to have a steeper structure function at a timescale less than the break timescale and smaller variability amplitude. By comparing the structure functions in the optical with the mid-infrared obtained from the Wide-field Infrared Survey Explorer, we derive the size of the dusty torus using a toy model for the geometry of the torus. The size of the torus positively correlates with the luminosity of the active nucleus, following a relation that agrees well with previous studies based on reverberation mapping. This result demonstrates that the structure function method can be used as a powerful and highly efficient tool to examine the size of the torus.

Sulfur chemistry in the formation process of low-mass stars and planets remains poorly understood. The protoplanetary disks (PPDs) are the birthplace of planets and its distinctive environment provides an intriguing platform for investigating models of sulfur chemistry. We analyzed the ALMA observations of CS 7-6 transitions in the HD 163296 disk and perform astrochemical modeling to explore its sulfur chemistry. We simulated the distribution of sulfur-containing molecules and compared it with observationally deduced fractional column densities. We have found that the simulated column density of CS is consistent with the observationally deduced fractional column densities, while the simulated column density of C$_2$S is lower than the observationally deduced upper limits on column densities. This results indicate that we have a good understanding of the chemical properties of CS and C$_2$S in the disk. We also investigated the influence of the C/O ratio on sulfur-containing molecules and found that the column densities of SO, SO$_2$, and H$_2$S near the central star are dependent on the C/O ratio. Additionally, we found that the $N$[CS]/$N$[SO] ratio can serve as a promising indicator of the disk's C/O ratio in the HD 163296. Overall, the disk of HD 163296 provides a favorable environment for the detection of sulfur-containing molecules.

We investigate the star formation process within the central 3.3 kpc region of the nearby luminous infrared Seyfert NGC 7469, probing scales ranging from 88 to 330 pc. We combine JWST/MIRI imaging with the F770W filter, with CO(2-1) and the underlying 1.3 mm dust continuum data from ALMA, along with VLA radio continuum observations at 22 GHz. NGC 7469 hosts a starburst ring which dominates the overall star formation activity. We estimate a global star formation rate SFR $\sim 11.5$ $\rm M_{\odot}~yr^{-1}$ from the radio at 22 GHz, and a cold molecular gas mass M(H2) $\sim$ 6.4 $\times$ $\rm 10^9 M_{\odot}$ from the CO(2-1) emission. We find that the 1.3 mm map shows a morphology remarkably similar to those traced by the 22 GHz and the 7.7 $\rm \mu m$ polycyclic aromatic hydrocarbon (PAH) emission observed with JWST. The three tracers reproduce the morphology of the starburst ring with good agreement. We further investigate the correlations between the PAHs, the star formation rate and the cold molecular gas. We find a stronger correlation of the PAHs with the star formation than with the CO, with steeper correlations within the starburst ring ($n > 2$) than in the outer region ($n < 1$). We derive the correlation between the star formation rate and the cold molecular gas mass surface densities, the Kennicutt-Schmidt star formation law. Comparisons with other galaxy populations, including starburst galaxies and active galactic nuclei, highlighted that NGC 7469 exhibits an intermediate behavior to the Kennicutt-Schmidt relations found for these galaxy populations.

Exoplanets are celestial bodies orbiting stars beyond our Solar System. Although historically they posed detection challenges, Kepler's data has revolutionized our understanding. By analyzing flux values from the Kepler Mission, we investigate the intricate patterns in starlight that may indicate the presence of exoplanets. This study investigates a novel approach for exoplanet classification using Spiking Neural Networks (SNNs) applied to data obtained from the NASA Kepler mission. SNNs offer a unique advantage by mimicking the spiking behavior of neurons in the brain, allowing for more nuanced and biologically inspired processing of temporal data. Experimental results demonstrate the efficacy of the proposed SNN architecture, excelling in various performance metrics such as accuracy, F1 score, precision, and recall.

Spectroscopic observations of the far-infrared [O III] and [C II] lines present a pathway to explore the mechanisms of the emergence of massive galaxies in the epoch of reionization and beyond, which is one of the most fundamental questions in astronomy. To address this question, the Far-Infrared Nebular Emission Receiver (FINER) project is developing two wideband dual-polarization sideband-separating heterodyne receivers at 120--210 GHz and 210--360 GHz for the Large Millimeter Telescope (LMT) in Mexico. Compared with Atacama Large Millimeter/submillimeter Array (ALMA), LMT provides 40% of ALMA's light-collecting area and a similar atmospheric transmittance, but FINER plans to have an instantaneous intermediate frequency (IF) of 3--21 GHz per sideband per polarization which is five times wider than current ALMA's bandwidth. Therefore, FINER is going to offer cutting-edge spectral scanning capability in the next several years. The project is currently in an active development phase. In this proceeding, the latest development status for FINER, including the optics, wideband waveguide components as well as low-noise superconductor-insulator-superconductor (SIS) mixers is reported.

Unveiling the emergence and prevalence of massive/bright galaxies during the epoch of reionization and beyond, within the first 600 million years of the Universe, stands as a pivotal pursuit in astronomy. Remarkable progress has been made by JWST in identifying an immense population of bright galaxies, which hints at exceptionally efficient galaxy assembly processes. However, the underlying physical mechanisms propelling their rapid growth remain unclear. With this in mind, millimeter and submillimeter-wave spectroscopic observations of redshifted far-infrared spectral lines, particularly the [O III] 88 micron and [C II] 158 micron lines, offers a crucial pathway to address this fundamental query. To this end, we develop a dual-polarization sideband-separating superconductor-insulator-superconductor (SIS) mixer receiver, FINER, for the Large Millimeter Telescope (LMT) situated in Mexico. Harnessing advancements from ALMA's wideband sensitivity upgrade (WSU) technology, FINER covers radio frequencies spanning 120-360 GHz, delivering an instantaneous intermediate frequency (IF) of 3-21 GHz per sideband per polarization, which is followed by a set of 10.24 GHz-wide digital spectrometers. At 40% of ALMA's light-collecting area, the LMT's similar atmospheric transmittance and FINER's 5 times wider bandwidth compared to ALMA culminate in an unparalleled spectral scanning capability in the northern hemisphere, paving the way for finer spectral-resolution detection of distant galaxies.

For efficient spectroscopic redshift identification of early galaxies in the northern hemisphere, we aim to combine the Large Millimeter Telescope (LMT) with a wide-band heterodyne receiver, FINER, which will cover radio frequencies of 120--360 GHz and offer a 3--21 GHz intermediate frequency (IF) per sideband and polarization. To take full advantage of such wide IFs, we present a novel 10.24-GHz-wide digital spectrometer, DRS4 (Elecs Industry Co., Ltd.). It incorporates 20.48 Gsps samplers with an FPGA-based digital signal processing module. To mitigate the noise contamination from the image sideband, it is equipped with a digital sideband separation function to improve the sideband rejection up to 25 dB. Laboratory performance evaluations show that it exhibits an Allan time of at least ~100 s and a total power dynamic range of at least 7 dB. These results demonstrate its capability of instantaneously wide-band spectroscopy toward high-redshift galaxies with position-switching observations.

Recent studies suggest that filamentary structures are representative of the initial conditions of star formation in molecular clouds and support a filament paradigm for star formation, potentially accounting for the origin of the stellar initial mass function (IMF). Using Herschel imaging observations of the California giant molecular cloud, we aim to further investigate the filament paradigm for low- to intermediate-mass star formation and to better understand the exact role of filaments in the origin of stellar masses. Using the multiscale, multiwavelength extraction method getsf, we identify starless cores, protostars, and filaments in the Herschel data set and separate these components from the background cloud contribution to determine accurate core and filament properties. Both the prestellar core mass function (CMF) and the distribution of filament masses per unit length or filament line mass function (FLMF) are consistent with power-law distributions at the high-mass end, $\Delta N/\Delta {\rm log}M\propto M^{-1.4 \pm 0.2}$ at $M > 1\,M_\odot$ for the CMF and $\Delta N/\Delta {\rm log} {M}_{\rm line} \propto {M}_{\rm line}^{-1.5\pm0.2}$ for the FLMF at $M_{\rm line} > 10\,M_\odot {\rm pc^{-1}}$, which are both consistent with the Salpeter power-law IMF. Based on these results, we propose a revised model for the origin of the CMF in filaments, whereby the global prestellar CMF in a molecular cloud arises from the integration of the CMFs generated by individual thermally supercritical filaments within the cloud. Our findings support the existence a tight connection between the FLMF and the CMF/IMF and suggests that filamentary structures represent a critical evolutionary step in establishing a Salpeter-like mass function.

Close encounter between a star and a supermassive black hole (SMBH) results in the tidal disruption of the star, known as a tidal disruption event (TDE). Recently, a few TDEs, e.g., ASASSN-15oi and AT2018hyz, have shown late-time (hundreds of days after their UV/optical peaks) radio flares with radio luminosities of $10^{38\sim39}$ erg/s. The super-Eddington fallback or accretion in a TDE may generate a mass outflow. Here we investigate a scenario that the late-time radio flares come from the interaction of the outflow with the circum-nuclear gaseous clouds, in addition to the slow-evolving emission component due to the outflow-diffuse medium interaction. We calculate the associated radio temporal and spectral signatures and find that they reproduce well the observations. The outflows have the inferred velocity of 0.2$\sim0.8$ c, the total mass of $10^{-3}\sim10^{-1}$ $\mathrm{M_{\odot}}$ and the ejection duration of a month to a year. The distances of the clouds to the SMBH are $0.1\sim1$ pc. This scenario has advantages in explaining the long delay, sharpness of the rise and the multiplicity of the late radio flares. Future observations may build up a much larger sample of late-time radio flares and enable their use as a probe of the TDE physics and the host circumnuclear environment.

In this work we compute the rates and numbers of different types of stars and phenomena (SNe, novae, white dwarfs, merging neutron stars, black holes) that contributed to the chemical composition of the Solar System. Stars die and restore the newly formed elements into the interstellar gas. This process is called "chemical evolution". In particular, we analyse the death rates of stars of all masses, dying either quiescently or explosively. These rates and total star numbers are computed in the context of a revised version of the two-infall model for the chemical evolution of the Milky Way, which reproduces fairly well the observed abundance patterns of several chemical species, as well as the global solar metallicity. We compute also the total number of stars ever born and still alive as well as the number of stars born up to the formation of the Solar System and with a mass and metallicity like the Sun. This latter number will account for all the possible existing Solar Systems which can host life in the solar vicinity. Among all the stars (from 0.8 to 100 M$_{\odot}$) born and died from the beginning up to the Solar System formation epoch, which contributed to its chemical composition, 93.00$\%$ are represented by stars dying as single white dwarfs (without interacting significantly with a companion star) and originating in the mass range 0.8-8 M$_{\odot}$, while 5.24$\%$ are neutron stars and 0.73$\%$ are black holes, both originating from SNe core-collapse (M>8 M$_{\odot}$); 0.64$\%$ are Type Ia SNe and 0.40$\%$ are nova systems, both originating from the same mass range as the white dwarfs. The number of stars similar to the Sun born from the beginning up to the Solar System formation, with metallicity in the range 12+log(Fe/H)= 7.50 $\pm$ 0.04 dex is 3.1732$\cdot$ 10$^{7}$, and in particular our Sun is the 2.6092$\cdot$ 10$^7$-th star of this kind, born in the solar vicinity.

We explore the evolution of various properties of multiple-population globular clusters (GCs) for a broad range of initial conditions. We simulated over 200 GC models using the MOCCA Monte Carlo code and find that present-day properties (core and half-light radii, ratio of the number of second-generation (SG) stars to the total number of stars, NSG/NTOT) of these models cover the observed values of these quantities for Milky Way GCs. Starting with a relatively small value of the SG fraction (NSG/NTOT ~ 0.25) and a SG system concentrated in the inner regions of the cluster, we find, in agreement with previous studies, that systems in which the first-generation (FG) is initially tidally filling or slightly tidally underfilling best reproduce the observed ratios of NSG/NTOT and have values of the core and half-light radii typical of those of many Galactic globular clusters. Models in which the FG is initially tidally underfilling retain values of NSG/NTOT close to their initial values. These simulations expand previous investigations and serve to further constrain the viable range of initial parameters and better understand their influence on present-day GC properties. The results of this investigation also provide the basis for our future survey aimed at building specific models to reproduce the observed trends (or lack thereof) between the properties of multiple stellar populations and other clusters properties.

Jupiters south polar region was observed by JWST Mid Infrared Instrument in December 2022. We used the Medium Resolution Spectrometer mode to provide new information about Jupiters South Polar stratosphere. The southern auroral region was visible and influenced the atmosphere in several ways. 1: In the interior of the southern auroral oval, we retrieved peak temperatures at two distinct pressure levels near 0.01 and 1 mbar, with warmer temperatures with respect to non auroral regions of 12 pm 2 K and 37 pm 4 K respectively. A cold polar vortex is centered at 65S at 10 mbar. 2: We found that the homopause is elevated to 590+25-118 km above the 1-bar pressure level inside the auroral oval compared to 460+60-50 km at neighboring latitudes and with an upper altitude of 350 km in regions not affected by auroral precipitation. 3: The retrieved abundance of C2H2 shows an increase within the auroral oval, and it exhibits high abundances throughout the polar region. The retrieved abundance of C2H6 increases towards the pole, without being localized in the auroral oval, in contrast with previous analysis. We determined that the warming at 0.01 mbar and the elevated homopause might be caused by the flux of charged particles depositing their energy in the South Polar Region. The 1 mbar hotspot may arise from adiabatic heating resulting from auroral driven downwelling. The cold region at 10 mbar may be caused by radiative cooling by stratospheric aerosols. The differences in spatial distribution seem to indicate that the hydrocarbons analyzed are affected differently by auroral precipitation.

The multi-scaled solar magnetic field consists of two major components: active regions (ARs) and magnetic network. Unraveling the cycle-dependent properties and interrelations of these components is crucial for understanding the evolution of the solar magnetic field. In this study, we investigate these components using magnetic power spectra derived from high-resolution and continuous synoptic magnetograms since cycle 23 onwards. Our results show that the size of the magnetic network ranges from 26 Mm to 41 Mm without dependence on the solar cycle. The power of the network field ($P_{NW}$) accounts for approximately 20\% of the total power during any phase of solar cycles. In contrast to the AR power ($P_{AR}$), $P_{NW}$ displays a weaker cycle dependence, as described by the relationship $P_{NW}$ $\approx$ 0.6* $P_{AR}$ + 40. The power-law index between AR sizes and magnetic network sizes presents a strong anti-correlation with the activity level. Additionally, our study indicates that in the absence of sunspots on the solar disc, the magnetic power spectra remain time-independent, consistently exhibiting similarity in both shape and power. This study introduces a new method to investigate the properties of the magnetic network and provides magnetic power spectra for high-resolution simulations of the solar magnetic field at the surface at various phases of solar cycles.

We are probing the gravitational force perpendicular to the Galactic plane at the position of the Sun with a sample of red giants whose measurements are taken from the DR3-Gaia catalogue. Measurements far out of the Galactic plane up to 3.5 kpc allow us to determine directly the total mass density where dark matter is dominant, while stellar and gas densities are very low. In a complementary way, we are also using a new determination of the local baryonic mass density to help us determine the density of dark matter in the Galactic plane at the solar position. We obtain for the local mass density of dark matter $\rho_\mathrm{dm}$=0.0128$\pm $0.0008 M$_\odot$ pc$^{-3}$ = 0.486 $\pm$0.030 Gev cm$^{-3}$, for the flattening of the gravitational potential of the dark halo $q_\mathrm{\phi,h}$=0.843$\pm0.035$, and for its density $q_\mathrm{\rho,h}$=0.781$\pm$0.055.

Despite recent progress in the spectroscopic characterization of individual exoplanets, the atmospheres of key ultra-hot Jupiters (UHJs) still lack comprehensive investigations. These include WASP-178b, one of the most irradiated UHJs known to date. We observed the dayside emission signal of this planet with CRIRES$^+$ in the spectral K-band. By applying the cross-correlation technique and a Bayesian retrieval framework to the high-resolution spectra, we identified the emission signature of $^{12}$CO (S/N = 8.9) and H$_2$O (S/N = 4.9), and a strong atmospheric thermal inversion. A joint retrieval with space-based secondary eclipse measurements from TESS and CHEOPS allows us to refine our results on the thermal profile and thus to constrain the atmospheric chemistry, yielding a solar to super-solar metallicity (1.4$\pm$1.6 dex) and a solar C/O ratio (0.6$\pm$0.2). We infer a significant excess of spectral line broadening and identify a slight Doppler-shift between the $^{12}$CO and H$_2$O signals. These findings provide strong evidence for a super-rotating atmospheric flow pattern and suggest the possible existence of chemical inhomogeneities across the planetary dayside hemisphere. In addition, the inclusion of photometric data in our retrieval allows us to account for stellar light reflected by the planetary atmosphere, resulting in an upper limit on the geometric albedo (0.23). The successful characterization of WASP-178b's atmosphere through a joint analysis of CRIRES$^+$, TESS, and CHEOPS observations highlights the potential of combined studies with space- and ground-based instruments and represents a promising avenue for advancing our understanding of exoplanet atmospheres.

Coronal mass ejections (CMEs) on stars can change the stars' magnetic field configurations and mass loss rates during the eruption and propagation and therefore, may affect the stars' rotation properties on long time-scales. The dynamics of stellar CMEs and their influence on the stellar angular momentum loss rate are not yet well understood. In order to start investigating these CME-related aspects on other stars, we conducted a series of magnetohydrodynamic simulations of CMEs on a solar-type star of moderate activity levels. The propagation and evolution of the CMEs were traced in the three-dimensional outputs and the temporal evolution of their dynamic properties (such as masses, velocities, and kinetic energies) were determined. The simulated stellar CMEs are more massive and energetic than their solar analog, which is a result of the stronger magnetic field on the surface of the simulated star than that of the Sun. The simulated CMEs display masses ranging from ~10^16 g to ~10^18 g and kinetic energies from ~10^31 erg to ~10^33 erg. We also investigated the instantaneous influence of the CMEs to the star's angular momentum loss rate. Our results suggest that angular momentum can either be added to or be removed from the star during the evolution of CME events. We found a positive correlation between the amplitude of the angular momentum loss rate variation and the CME's kinetic energy as well as mass, suggesting that more energetic/massive CMEs have higher possibility to add angular momentum to the star.

The prevalence of light element enhancement in the most metal-poor stars is potentially an indication that the Milky Way has a metallicity floor for star formation around $\sim$10$^{-3.5}$ Z$_{\odot}$. We propose that this metallicity floor has its origins in metal-enriched star formation in the minihalos present during the Galaxy's initial formation. To arrive at this conclusion, we analyze a cosmological radiation hydrodynamics simulation that follows the concurrent evolution of multiple Population III star-forming minihalos. The main driver for the central gas within minihalos is the steady increase in hydrostatic pressure as the halos grow. We incorporate this insight into a hybrid one-zone model that switches between pressure-confined and modified free-fall modes to evolve the gas density with time according to the ratio of the free-fall and sound-crossing timescales. This model is able to accurately reproduce the density and chemo-thermal evolution of the gas in each of the simulated minihalos up to the point of runaway collapse. We then use this model to investigate how the gas responds to the absence of H$_{2}$. Without metals, the central gas becomes increasingly stable against collapse as it grows to the atomic cooling limit. When metals are present in the halo at a level of $\sim$10$^{-3.7}$ Z$_{\odot}$, however, the gas is able to achieve gravitational instability while still in the minihalo regime. Thus, we conclude that the Galaxy's metallicity floor is set by the balance within minihalos of gas-phase metal cooling and the radiation background associated with its early formation environment.

Galaxy interactions in groups can lead to intense starbursts and the activation of active galactic nuclei (AGNs). The stripped gas from the outer disk can lead to star-forming clumps along the tidal tails or sometimes tidal dwarf galaxies. We investigate the impact of interaction on various galaxy properties, including morphology, star formation rates, and chemical composition in the galaxy group AM\,1054-325 using Multi Unit Spectroscopic Explorer (MUSE) data. We conduct a comprehensive spatially and spectrally resolved investigation of the star formation rate, star formation histories, metallicity, and AGN activity. The galaxy subgroup AM\,1054-325A shows multiple star-forming clumps in H$\alpha$ emission along the western tidal tail, which are formed due to tidal stripping. These clumps also have higher metallicities. AM\,1054-325B is quenched and shows disturbed gas kinematics and the signature of gas accretion in the H$\alpha$ map. The specific star formation along the tidal tail is higher, contributing to the galaxy's overall stellar mass growth.

We have started a measurement campaign of numerous methanol isotopologs in low-lying torsional states in order to provide extensive line lists for radio astronomical observations from an adequate spectroscopic model and to investigate how the intricate vibration-torsion-rotation interactions manifest themselves in the spectra of different isotopic species. After CD$_3$OH and CD$_3$OD, we turn our focus to CH$_3$OD, which is an important species for studying deuteration in prestellar cores and envelopes that enshroud protostars. Notably, deuteration is frequently viewed as a diagnostic tool for star formation. The measurements used in this study were obtained in two spectroscopic laboratories and cover large fractions of the 34 GHz--1.35 THz range. As done in previous studies, we employed a torsion-rotation Hamiltonian model for our analysis that is based on the rho-axis method. The resulting model describes the ground and first excited torsional states of CH$_3$OD well up to quantum numbers $J \leqslant 51$ and $K_a \leqslant 18$. We derived a line list for radio astronomical observations from this model that is accurate up to at least 1.35~THz and should be sufficient for all types of radio astronomical searches for this methanol isotopolog in these two lowest torsional states. This line list was applied to a reinvestigation of CH$_3$OD in data from the Protostellar Interferometric Line Survey of IRAS 16293--2422 obtained with the Atacama Large Millimeter/submillimeter Array. The new accurately determined value for the column density of CH$_3$OD implies that the deuteration in methanol differs in its two functional groups by a factor of $\sim$7.5.

We present a complete numerical model of the afterglow of a laterally-structured relativistic ejecta from radio to very high energy (VHE). This includes a self-consistent calculation of the synchrotron radiation, with its maximum frequency, and of Synchrotron Self-Compton (SSC) scatterings taking into account the Klein-Nishina regime. The attenuation due to pair production is also included. This model is computationally-efficient, allowing for multi-wavelength data fitting. As a validation test, the radiative model is used to fit the broad-band spectrum of GRB 190114C at 90 s up to the TeV range. The full model is then used to fit the afterglow of GW 170817 and predict its VHE emission. We find that the SSC flux at the peak was much dimmer than the upper limit from H.E.S.S. observations. However, we show that either a smaller viewing angle or a larger external density would make similar off-axis events detectable in the future at VHE, even above 100 Mpc with the sensitivity of the CTA. Large external densities are expected in the case of fast mergers, but the existence of a formation channel for such binary neutron stars is still uncertain. We highlight that VHE afterglow detections would help probing efficiently such systems.

We present a multi-wavelength analysis of the young star cluster Berkeley 59 (Be 59) based on the $Gaia$ data and deep infrared (IR) observations with the 3.58-m Telescopio Nazionale Galileo and $Spitzer$ space telescope. The mean proper motion of the cluster is found to be $\mu$$_\alpha$cos$\delta$ $\sim$ -0.63 mas yr$^{-1}$ and $\mu$$_\delta$ $\sim$ -1.83 mas yr$^{-1}$ and the kinematic distance of the cluster, $\sim$ 1 kpc, is in agreement with previous photometric studies. Present data is the deepest available near-IR observations for the cluster so far and reached below 0.03 M$_\odot$. The mass function of the cluster region is calculated using the statistically cleaned color-magnitude diagram and is similar to the Salpeter value for the member stars above 0.4 M$_\odot$. In contrast, the slope becomes shallower ($\Gamma$ $\sim$ 0.01 $\pm$ 0.18) in the mass range 0.04 - 0.4 M$_\odot$, comparable to other nearby clusters. The spatial distribution of young brown dwarfs (BDs) and stellar candidates shows a non-homogeneous distribution. This suggests that the radiation feedback from massive stars may be a prominent factor contributing to the BD population in the cluster Be 59. We also estimated the star-to-BD ratio for the cluster, which is found to be $\sim$ 3.6. The Kolomogorov-Smirnov test shows that stellar and BD populations significantly differ, and stellar candidates are near the cluster center compared to the BDs, suggesting mass segregation in the cluster toward the substellar mass regime.

We developed a modular silicon photomultiplier camera to detect Earth-skimming PeV to EeV tau neutrinos with the imaging atmospheric Cherenkov technique. We built two cameras, a 256-pixel camera with S14161-6050HS SiPMs for the Trinity Demonstrator located on Frisco Peak, Utah, and a 512-pixel camera with S14521-6050AN SiPMs for the EUSO-SPB2 Cherenkov Telescope. The front-end electronics are based on the eMUSIC ASIC, and the camera signals are sampled and digitized with the 100MS/s and 12-bit AGET system. Both cameras are liquid-cooled. We detail the camera concept and the results from characterizing the SiPMs, bench testing, and calibrating the two cameras.

The yellow hypergiant star V509 Cas is currently undergoing an extreme phase of evolution. Having experienced eruptive mass-loss outbursts in the 20th century, the star's effective temperature reached record high values in the early 2000s. However, since then, the star's behaviour has displayed an unprecedented level of stability. In spite of that, the star could be traversing through the 'yellow void' instability region. To describe the current evolutionary state of V509 Cas, we analysed its variability using photometric and spectroscopic data collected over recent years. By comparing our findings with historical records, we aim to determine whether the star's surface shows signs of stabilisation. Additionally, we investigate the variability of emission components in the wings of certain spectral lines to highlight the contribution of the circumstellar gaseous disc to this phenomenon. Our spectroscopic monitoring observations were carried out at Tartu Observatory over the course of seven years, supplemented by echelle spectra obtained at the Nordic Optical Telescope, as well as publicly available photometric data from Gaia, AAVSO, and AAVSO's Bright Star Monitor programme. We estimated the variability of effective temperature and radial velocity from the spectral time series and correlated it with the brightness variability of V509 Cas. The results indicate that the star's average brightness level has remained stable throughout the observed period, with an amplitude of variability ~0.1 mag. While the amplitude of short-term temperature fluctuations has decreased compared to the early 2000s, the variability of the radial velocity remains similar to historical values from the early 20th century. Moreover, we show how the variable radial velocity affects the emission components in some absorption lines (e.g. Sc II) and how that follows the hypothesis of a disc surrounding the star.

Accurate flux density calibration is essential for precise analysis and interpretation of observations across different observation modes and instruments. In this research, we firstly introduce the flux calibration model incorporated in HIFAST pipeline, designed for processing HI 21-cm spectra. Furthermore, we investigate different calibration techniques and assess the dependence of the gain parameter on the time and environmental factors. A comparison is carried out in various observation modes (e.g. tracking and scanning modes) to determine the flux density gain ($G$), revealing insignificant discrepancies in $G$ among different methods. Long-term monitoring data shows a linear correlation between $G$ and atmospheric temperature. After subtracting the $G$--Temperature dependence, the dispersion of $G$ is reduced to $<$3% over a one-year time scale. The stability of the receiver response of FAST is considered sufficient to facilitate HI observations that can accommodate a moderate error in flux calibration (e.g., $>\sim5\%$) when utilizing a constant $G$ for calibration purposes. Our study will serve as a useful addition to the results provided by Jiang et al. (2020). Detailed measurement of $G$ for the 19 beams of FAST, covering the frequency range 1000 MHz -- 1500 MHz can be found on the HIFAST homepage: this https URL.

This paper delves into the late-time accelerated expansion of the universe and the evolution of cosmic structures within the context of a specific \( f(R, L_m) \) gravity model, formulated as \( f(R, L_m) = \lambda R + \beta L_m^\alpha + \eta \). To study the cosmological viability of the model, we employed the latest cosmic measurement datasets: i) 57 observational Hubble parameter data points (\texttt{OHD}); ii) 1048 distance moduli data points (\texttt{SNIa}); iii) a combined dataset (\texttt{OHD+SNIa}); and large scale structure datasets, including iv) 14 growth rate data points (\texttt{f}); and v) 30 redshift space distortion data points (\texttt{f}$\sigma_8$). These datasets facilitated the constraint of the \( f(R, L_m) \)-gravity model via MCMC simulations, followed by a comparative analysis with the \(\Lambda\)CDM model. A comprehensive statistical analysis has been conducted to evaluate the \( f(R, L_m) \)-gravity model's efficacy in explaining both the accelerated expansion of the universe and the growth of cosmic structures.

The Near-InfraRed Planet Searcher or NIRPS is a precision radial velocity spectrograph developed through collaborative efforts among laboratories in Switzerland, Canada, Brazil, France, Portugal and Spain. NIRPS extends to the 0.98-1.8 $\mu$m domain of the pioneering HARPS instrument at the La Silla 3.6-m telescope in Chile and it has achieved unparalleled precision, measuring stellar radial velocities in the infrared with accuracy better than 1 m/s. NIRPS can be used either stand-alone or simultaneously with HARPS. Commissioned in late 2022 and early 2023, NIRPS embarked on a 5-year Guaranteed Time Observation (GTO) program in April 2023, spanning 720 observing nights. This program focuses on planetary systems around M dwarfs, encompassing both the immediate solar vicinity and transit follow-ups, alongside transit and emission spectroscopy observations. We highlight NIRPS's current performances and the insights gained during its deployment at the telescope. The lessons learned and successes achieved contribute to the ongoing advancement of precision radial velocity measurements and high spectral fidelity, further solidifying NIRPS' role in the forefront of the field of exoplanets.

The quasar 3C 286, a well-known calibrator source in radio astronomy, was found to exhibit exceptional multiwavelength properties. Its rich and complex optical emission-line spectrum revealed its narrow-line Seyfert 1 (NLS1) nature. Given its strong radio emission, this makes 3C 286 one of the radio-loudest NLS1 galaxies known to date. 3C 286 is also one of very few known compact steep-spectrum (CSS) sources detected in the gamma-ray regime. Observations in the X-ray regime, rarely carried out so far, revealed evidence for variability, raising the question if driven by the accretion disk or jet. 3C 286 is also well known for its damped Lyman alpha system from an intervening absorber at z = 0.692, triggering a search for the corresponding X-ray absorption along the line-of-sight. Here, we present new observations in the radio, X-ray, optical and UV band. The nature of the X-ray variability is addressed. Spectral evidence suggests that it is primarily driven by the accretion disk (not the jet), and the X-ray spectrum is well fit by a powerlaw plus soft excess model. The radio flux density and polarization remain constant at the Effelsberg telescope resolution, reconfirming the use of 3C 286 as radio calibrator. The amount of reddening/absorption along the line-of-sight {\em{intrinsic}} to 3C 286 is rigorously assessed. None is found, validating the derivation of a high Eddington ratio (L/L-Edd ~ 1) and of the very high radio-loudness index of 3C 286. Based on the first deep Chandra image of 3C 286, tentative evidence for hard X-ray emission from the SW radio lobe is reported. A large variety of models for the gamma-ray emission of 3C 286 is briefly discussed.

Blueberry galaxies (BBs) are fainter, less massive, and lower redshift counterparts of the Green pea galaxies. They are thought to be the nearest analogues of the high redshift Lyman Alpha (Ly$\alpha$) emitters. We report the interferometric imaging of HI 21 cm emission from a Blueberry galaxy, J1509+3731, at redshift, z = 0.03259, using the Giant Metrewave Radio Telescope (GMRT). We find that this Blueberry galaxy has an HI mass of $M_{\text{HI}} \approx 3\times 10^8 \, M_{\odot}$ and an HI-to-stellar mass ratio $M_{\text{HI}}/M_* \approx$ 2.4. Using SFR estimates from the H$\beta$ emission line, we find that it has a short HI depletion time scale of $\approx 0.2$ Gyr, which indicates a significantly higher star-formation efficiency compared to typical star-forming galaxies at the present epoch. Interestingly, we find an offset of $\approx 2$ kpc between the peak of the HI 21 cm emission and the optical centre which suggests a merger event in the past. Our study highlights the important role of mergers in triggering the starburst in BBs and their role in the possible leakage of Lyman-$\alpha$ and Lyman-continuum photons which is consistent with the previous studies on BB galaxies.

The aim of this study is to probe the sub-mJy polarized source population with LOFAR. We present the method used to stack LOFAR polarization datasets, the resulting catalog of polarized sources, and the derived polarized source counts. The ELAIS-N1 field was selected for a polarimetric study at 114.9-177.4 MHz. A total area of 25 deg2 was imaged at 6"- resolution in the Stokes Q and U parameters. Alignment of polarization angles was done both in frequency and in Faraday space before stacking datasets from 19 eight-hour-long epochs. A search for polarized sources was carried out in the final, stacked dataset, and the properties of the detected sources were examined. The depolarization level of sources known to be polarized at 1.4 GHz was quantified. A one-sigma noise level of 19 {\mu}Jy/beam was reached in the central part of the field after stacking. Twenty-five polarized sources were detected above 8\sigma, five of which had not been detected in polarization at any other radio frequencies before. Seven additional polarized components were found by lowering the threshold to 6\sigma at positions corresponding to sources known to be polarized at 1.4 GHz. In two radio galaxies, polarization was detected from both radio lobes, so the final number of associated radio continuum sources is 31. The detected sources are weakly polarized, with a median degree of polarization of 1.75% for the sample of sources detected in polarized emission. The sources previously detected in polarization at 1.4 GHz are significantly depolarized at 150 MHz. The catalog is used to derive the polarized source counts at 150 MHz. This is the deepest and highest-resolution polarization study at 150 MHz to date.

We search for atmospheric constituents for the UHJ WASP-178 b with two ESPRESSO transits using the narrow-band and cross-correlation techniques, focusing on the detections of NaI, H$\alpha$, H$\beta$, H$\gamma$, MgI, FeI and FeII. Additionally, we show parallel photometry used to obtain updated and precise stellar, planetary and orbital parameters. We report the resolved line detections of NaI (5.5 and 5.4 $\sigma$), H$\alpha$ (13 $\sigma$), H$\beta$ (7.1 $\sigma$), and tentatively MgI (4.6 $\sigma$). In cross-correlation, we confirm the MgI detection (7.8 and 5.8 $\sigma$) and additionally report the detections of FeI (12 and 10 $\sigma$) and FeII (11 and 8.4 $\sigma$), on both nights separately. The detection of MgI remains tentative, however, due to the differing results between both nights, as well as compared with the narrow-band derived properties. None of our resolved spectral lines probing the mid- to upper atmosphere show significant shifts relative to the planetary rest frame, however H$\alpha$ and H$\beta$ exhibit line broadenings of 39.6 $\pm$ 2.1 km/s and 27.6 $\pm$ 4.6 km/s, respectively, indicating the onset of possible escape. WASP-178 b differs from similar UHJ with its lack of strong atmospheric dynamics in the upper atmosphere, however the broadening seen for FeI (15.66 $\pm$ 0.58 km/s) and FeII (11.32 $\pm$ 0.52 km/s) could indicate the presence of winds in the mid-atmosphere. Future studies on the impact of the flux variability caused by the host star activity might shed more light on the subject. Previous work indicated the presence of SiO cloud-precursors in the atmosphere of WASP-178 b and a lack of MgI and FeII. However, our results suggest that a scenario where the planetary atmosphere is dominated by MgI and FeII is more likely. In light of our results, we encourage future observations to further elucidate these atmospheric properties.

We propose a novel approach using neural networks (NNs) to differentiate between cosmological models, especially in the case where they are nested and the additional model parameters are close to zero, making it difficult to discriminate them with traditional approaches. Our method complements Bayesian analyses for cosmological model selection, which heavily depend on the chosen priors and average the unnormalized posterior over potentially large prior volumes. By analyzing simulated realistic data sets of the growth rate of the large scale structure (LSS) of the Universe, based on current galaxy-clustering survey specifications, for the cosmological constant and cold dark matter ($\Lambda$CDM) model and the Hu-Sawicki $f(R)$ model, we demonstrate the potential of NNs to enhance the extraction of meaningful information from cosmological LSS data. We find that the NN can successfully distinguish between $\Lambda$CDM and the $f(R)$ models, by predicting the correct model with approximately $97\%$ overall accuracy, thus demonstrating that NNs can maximise the potential of current and next generation surveys to probe for deviations from general relativity.

We present a systematic search for tidal disruption events (TDEs) using radio data from the Variables and Slow Transients (VAST) Pilot Survey conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). Historically, TDEs have been identified using observations at X-ray, optical, and ultraviolet wavelengths. After discovery, a few dozen TDEs have been shown to have radio counterparts through follow-up observations. With systematic time-domain radio surveys becoming available, we can now identify new TDEs in the radio regime. A population of radio-discovered TDEs has the potential to provide several key insights including an independent constraint on their volumetric rate. We conducted a search to select variable radio sources with a single prominent radio flare and a position consistent within 2$\sigma$ of the nucleus of a known galaxy. While TDEs were the primary target of our search, sources identified in this search may also be consistent with active galactic nuclei exhibiting unusual flux density changes at the timescales probed, uncharacteristically bright supernovae, or a population of gamma-ray bursts. We identify a sample of 12 radio-bright candidate TDEs. The timescales and luminosities range from ~6 to 230 days and ~10$^{38}$ to 10$^{41}$ erg s$^{-1}$, consistent with models of radio emission from TDEs that launch relativistic jets. After calculating the detection efficiency of our search using a Monte Carlo simulation of TDEs, and assuming all 12 sources are jetted TDEs, we derive a volumetric rate for jetted TDEs of 0.80$^{+0.31}_{-0.23}$ Gpc$^{-3}$ yr$^{-1}$, consistent with previous empirically estimated rates.

Early-generation in-situ dust detectors in near-Earth space have reported the occurrence of clusters of sub-micron dust particles that seemed unrelated to human spaceflight activities. In particular, data from the impact ionization detector onboard the HEOS-2 satellite indicate that such swarms of particles occur throughout the Earth's magnetosphere up to altitudes of 60,000 km -- far beyond regions typically used by spacecraft. Further account of high-altitude clusters has since been given by the GEO-deployed GORID detector, however, explanations for the latter have so far only been sought in GEO spaceflight activity. This perspective piece reviews dust cluster detections in near-Earth space, emphasizing the natural swarm creation mechanism conjectured to explain the HEOS-2 data -- that is, the electrostatic disruption of meteoroids. Highlighting this mechanism offers a novel viewpoint on more recent near-Earth dust measurements. We further show that the impact clusters observed by both HEOS-2 and GORID are correlated with increased geomagnetic activity. This consistent correlation supports the notion that both sets of observations stem from the same underlying phenomenon and aligns with the hypothesis of the electrostatic breakup origin. We conclude that the nature of these peculiar swarms remains highly uncertain, advocating for their concerted investigation by forthcoming dust science endeavours, such as the JAXA/DLR DESTINY+ mission.

We report the discovery of a compact star-forming galaxy at $z=9.380$ in the GOODS-North field (named GN-z9p4) which shows numerous strong UV-optical emission lines and a single UV line, NIV] 1486. This makes GN-z9p4 the third-highest redshift N-emitter known to date. We determine the nebular abundances of H, C, N, O and Ne, size, and other physical properties of this object, and compare them to those of the other N-emitters known so far and to other star-forming galaxies. Using the direct method we find a metallicity 12+log(O/H)$=7.37 \pm 0.15$, one of the lowest among the N-emitters. The N/O abundance ratio is highly super-solar, and C/O and Ne/O normal compared to other galaxies at low metallicity. We show that the compactness of GN-z9p4 (with effective radius $118\pm16$ pc at 2 micron) and other N-emitters translates into very high stellar mass and SFR surface densities, which could be a criterium to identify other N-emitters. Future studies and larger samples are needed to understand these recently discovered, rare, and enigmatic objects.

Gravitational lensing data is frequently collected at low resolution due to instrumental limitations and observing conditions. Machine learning-based super-resolution techniques offer a method to enhance the resolution of these images, enabling more precise measurements of lensing effects and a better understanding of the matter distribution in the lensing system. This enhancement can significantly improve our knowledge of the distribution of mass within the lensing galaxy and its environment, as well as the properties of the background source being lensed. Traditional super-resolution techniques typically learn a mapping function from lower-resolution to higher-resolution samples. However, these methods are often constrained by their dependence on optimizing a fixed distance function, which can result in the loss of intricate details crucial for astrophysical analysis. In this work, we introduce $\texttt{DiffLense}$, a novel super-resolution pipeline based on a conditional diffusion model specifically designed to enhance the resolution of gravitational lensing images obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Our approach adopts a generative model, leveraging the detailed structural information present in Hubble Space Telescope (HST) counterparts. The diffusion model, trained to generate HST data, is conditioned on HSC data pre-processed with denoising techniques and thresholding to significantly reduce noise and background interference. This process leads to a more distinct and less overlapping conditional distribution during the model's training phase. We demonstrate that $\texttt{DiffLense}$ outperforms existing state-of-the-art single-image super-resolution techniques, particularly in retaining the fine details necessary for astrophysical analyses.

In our previous work, we introduced a method that combines two unsupervised algorithms: DBSCAN and GMM. We applied this method to 12 open clusters based on Gaia EDR3 data, demonstrating its effectiveness in identifying reliable cluster members within the tidal radius. However, for studying cluster morphology, we need a method capable of detecting members both inside and outside the tidal radius. By incorporating a supervised algorithm into our approach, we successfully identified members beyond the tidal radius. In our current work, we initially applied DBSCAN and GMM to identify reliable members of cluster stars. Subsequently, we trained the Random Forest algorithm using DBSCAN and GMM-selected data. Leveraging the random forest, we can identify cluster members outside the tidal radius and observe cluster morphology across a wide field of view. Our method was then applied to 15 open clusters based on Gaia DR3, which exhibit a wide range of metallicity, distances, members, and ages. Additionally, we calculated the tidal radius for each of the 15 clusters using the King profile and detected stars both inside and outside this radius. Finally, we investigated mass segregation and luminosity distribution within the clusters. Overall, our approach significantly improved the estimation of the tidal radius and detection of mass segregation compared to previous work. We found that in Collinder 463, low-mass stars do not segregate in comparison to high-mass and middle-mass stars. Additionally, we detected a peak of luminosity in the clusters, some of which were located far from the center, beyond the tidal radius.

We search for physically consistent realizations of evolving dark energy suggested by the cosmological fit of DESI, Planck and Supernovae data. First we note that any lagrangian description of the standard Chevallier-Linder-Polarski (CLP) parametrization for the dark energy equation of state $w_{\rm eff}$, allows for the addition of a cosmological constant. We perform the cosmological fit finding new regions of parameter space that however continue to favour dark energy with $w_{\rm eff}<-1$ at early times that are challenging to realize in theoretically consistent theories. Next, in the spirit of effective field theories, we consider the effect of higher order terms in the Taylor expansion of the equation of state around the present epoch. We find that non-linear corrections of the equation of state are weakly constrained, thus opening the way to scenarios that differ from CLP at early times, possibly with $w_{\rm eff}>-1$ at all times. We present indeed scenarios where evolving dark energy can be realized through quintessence models. We introduce in particular the ramp model where dark energy coincides with CLP at late times and approximates to a cosmological constant at early times. The latter model provides a much better fit than $\Lambda$CDM, and only slightly worse than $w_0w_a$CDM, but with the big advantage of being described by a simple and theoretically consistent lagrangian of a canonical quintessence model.

The classical Cepheids CE Cas A, CE Cas B, CF Cas, and CG Cas are likely members of the binary open cluster comprising NGC 7790 and Berkeley 58. The clusters are of comparable age and in close proximity, as deduced from differentially dereddened $UuB_PBVGR_P$ photometry, and Cepheid period-age relations. Gaia DR3 astrometric and spectroscopic solutions for the clusters are likewise consistent. Conversely, the seemingly adjacent open cluster NGC 7788 is substantially younger and less distant.

In the past decade, the diversity of strong lens modeling methods has exploded, from being purely analytical to pixelated and non-parametric, or based on deep learning. We embrace this diversity by selecting different software packages and use them to blindly model independently simulated Hubble Space Telescope imaging data. To overcome the difficulties arising from using different codes and conventions, we use the COde-independent Organized LEns STandard (COOLEST) to store, compare and release all models in a self-consistent and human-readable manner. From an ensemble of six modeling methods, we study the recovery of the lens potential parameters and properties of the reconstructed source. In particular, we simulate and infer parameters of an elliptical power-law mass distribution with external shear for the lens while each modeling method reconstructs the source differently. We find that overall, both lens and source properties are recovered reasonably well, but systematic biases arise in all methods. Interestingly, we do not observe that a single method is significantly more accurate than others, and the amount of bias largely depends on the specific lens or source property of interest. By combining posterior distributions from individual methods using equal weights, the maximal systematic biases on lens model parameters inferred from individual models are reduced by a factor of 5.4 on average. We investigate a selection of modeling effects that partly explain the observed biases, such as the cuspy nature of the background source and the accuracy of the point spread function. This work introduces, for the first time, a generic framework to compare and ease the combination of models obtained from different codes and methods, which will be key to retain accuracy in future strong lensing analyses.