Papers for Wednesday, Apr 17 2024

ChatGPT has been the most talked-about concept in recent months, captivating both professionals and the general public alike, and has sparked discussions about the changes that artificial intelligence (AI) will bring to the world. As physicists and astrophysicists, we are curious about if scientific data can be correctly analyzed by large language models (LLMs) and yield accurate physics. In this article, we fine-tune the generative pre-trained transformer (GPT) model by the astronomical data from the observations of galaxies, quasars, stars, gamma-ray bursts (GRBs), and the simulations of black holes (BHs), the fine-tuned model demonstrates its capability to classify astrophysical phenomena, distinguish between two types of GRBs, deduce the redshift of quasars, and estimate BH parameters. We regard this as a successful test, marking the LLM's proven efficacy in scientific research. With the ever-growing volume of multidisciplinary data and the advancement of AI technology, we look forward to the emergence of a more fundamental and comprehensive understanding of our universe. This article also shares some interesting thoughts on data collection and AI design. Using the approach of understanding the universe - looking outward at data and inward for fundamental building blocks - as a guideline, we propose a method of series expansion for AI, suggesting ways to train and control AI that is smarter than humans.

We present a catalog of 150 variable stars, including 13 stars with exoplanet candidates. 37 stars were identified as variables for the first time. As a result of a 2.5-year photometric survey of exoplanets, we have obtained and analyzed light curves for almost 50 thousand stars in fields around white dwarfs WD 0009+501 and GRW +708247. Here we describe observations and data processing, the search for variable stars, their cross-identification with other catalogs and classification. The catalog is published in open access and contains the primary classification of variability, light curves and their parameters.

The tidal disruption event (TDE) AT2022cmc represents the fourth known example of a relativistic jet produced by the tidal disruption of a stray star providing a unique probe of the formation and evolution of relativistic jets in otherwise dormant supermassive black holes (SMBHs). Here we present deep, late-time Chandra observations of AT2022cmc extending to $t_{\rm obs} \approx 400$ days after disruption. Our observations reveal a sudden decrease in the X-ray brightness by a factor of $\gtrsim 14$ over a factor of $\approx 2.3$ in time, and a deviation from the earlier power-law decline with a steepening $\alpha \gtrsim 3.2$ ($F_X \propto t^{-\alpha}$), steeper than expected for a jet break, and pointing to the cessation of jet activity at $t_{\rm obs} \approx 215$ days. Such a transition has been observed in two previous TDEs (Swift J1644+57 and Swift J2058+05). From the X-ray luminosity and the timescale of jet shutoff, we parameterize the mass of the SMBH in terms of unknown jet efficiency and accreted mass fraction parameters. Motivated by the disk-jet connection in AGN, we favor black hole masses $\lesssim 10^5 \ \rm M_{\odot}$ (where the jet and disk luminosities are comparable), and disfavor larger black holes (in which extremely powerful jets are required to outshine their accretion disks). We additionally estimate a total accreted mass of $\approx 0.1 \rm \ M_{\odot}$. Applying the same formalism to Swift J1644+57 and Swift J2058+05, we favor comparable black hole masses for these TDEs of $\lesssim$ a few $\times 10^5 \ \rm M_{\odot}$, and suggest that jetted TDEs may preferentially form from lower mass black holes when compared to non-relativistic events, owing to generally lower jet and higher disk efficiencies at higher black hole masses.

We report JWST/NIRSpec observations of a star-forming galaxy at $z=2.76$, MACSJ1149-WR1. We securely detect two temperature-sensitive auroral lines, [SIII]6312 (7.4$\sigma$) and [OII]7320+7331 doublets (10$\sigma$), and tentatively [NII]5755 ($2.3\sigma$) for the first time in an individual galaxy at $z>1$. We perform a detailed analysis of its interstellar media (ISM), and derive electron temperatures, various heavy element abundances (O/H, N/O, S/O, and Ar/O) in the hot ionized region, and the neutral fraction in the warm ionized region. MACSJ1149-WR1 shows a broad feature at the wavelength of HeII 4686, which consists of a broad ($\sim1000$km/s), blue-shifted ($\sim-110$km/s) line component. Taken together with its mildly elevated N/O abundance, we conclude that MACSJ1149-WR1 is experiencing a young starburst ($<10$Myr), likely hosting a large number of Wolf-Rayet stars. None of its spectral features support the presence of AGN, including: $(i)$ the absence of broad components and velocity shifts in Hydrogen recombination lines, $(ii)$ low [FeII]${\rm 1.257 \mu m}$/Pa$\beta$ ratio, and $(iii)$ the absence of high-ionization lines. Our analysis using HeI lines reveals a higher electron temperature and a higher attenuation value, indicating that HeI may probe a smaller spatial scale than HI, presumably the region dominated by the aforementioned Wolf-Rayet stars. The star formation rates derived from various HeI lines broadly agree with those from Hydrogen recombination lines. We thus advocate that HeI can be an excellent, independent probe of multi-phase ISM in the era of JWST.

The increasingly neutral intergalactic gas at $z>6$ impacts the Lyman-$\alpha$ flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Ly$\alpha$ line at $z=6.6$, unseen in any simulation of reionization. Here we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin$^2$ field centered on COLA1. We find 141 galaxies spectroscopically-selected through the [OIII]($\lambda4969,5008$) doublet at $5.35<z<6.95$, with 40 of these sources showing H$\beta$. For COLA1 we additionally detect [OIII]$_{4363}$ and H$\gamma$. We measure a systemic redshift of $z=6.5917$ for COLA1, confirming the double-peak nature of the Ly$\alpha$ profile. This implies that it resides in a highly ionized bubble and that it is leaking ionizing photons with a high escape fraction $f_{\rm esc}{\rm (LyC)}=20$-$50$%, making it a prime laboratory to study Lyman continuum escape in the Epoch of Reionization. COLA1 shows all the signs of a prolific ionizer with a Ly$\alpha$ escape fraction of $81\pm5\%$, Balmer decrement indicating no dust, a steep UV slope ($\beta_{\rm UV}=-3.2\pm 0.4$), and a star-formation surface density $\gtrsim 10\times$ that of typical galaxies at similar redshift. We detect 5 galaxies in COLA1's close environment ($\Delta z<0.02$). Exploiting the high spectroscopic completeness inherent to grism surveys, and using mock simulations that mimic the selection function, we show the that number of detected companions is very typical for a similarly UV-bright ($M_{\rm{UV}}\sim-21.3$) galaxy; that is, the ionized bubble around COLA1 is unlikely due to an excessively large over-density. Instead, the measured ionizing properties suggest that COLA1 by itself might be powering the bubble required to explain its double-peaked Ly$\alpha$ profile ($R_{\rm ion}\approx0.7$ pMpc), with minor contribution from detected neighbours ($-17.5>M_{\rm UV}>-19.5$).

Column density ratios of complex organic molecules are generally constant across protostellar systems with some low-level scatter. However, the scatter in formamide (NH$_2$CHO) to methanol (CH$_3$OH) column density ratio is one of the highest. This larger scatter is sometimes interpreted as evidence of gas-phase formation of NH$_2$CHO. In this work we propose an alternative interpretation in which this scatter is produced by differences in the snowline locations related to differences in binding energies of these species and the small-scale structure of the envelope and the disk system. We also include CH$_3$CN in our work as a control molecule which has a similar binding energy to CH$_3$OH. We use radiative transfer models to calculate the emission from these species in protostellar systems with and without disks. The abundances of these species are parameterized in our models. Then we fit the calculated emission lines to find the column densities as done in real observations. We find a correction factor of ~10 to be multiplied by gas-phase $N_{NH_2CHO}/N_{CH_3OH}$ to give the true abundance ratio of these two species in the ices. We find that models with different physical parameters produce a scatter in $N_{NH_2CHO}/N_{CH_3OH}$, comparable with that of observations. The scatter in $N_{NH_2CHO}/N_{CH_3OH}$ is larger than that of $N_{CH_3CN}/N_{CH_3OH}$ in models consistent with the observations. We show that the scatter in $N_{NH_2CHO}/N_{CH_3OH}$ will be lower if we correct for the difference in sublimation temperatures of these two species in observations of ~40 protostellar systems with ALMA. The scatter in $N_{NH_2CHO}/N_{CH_3OH}$ can be partially explained by the difference in their binding energies. We conclude that gas-phase chemistry routes for NH$_2$CHO are not necessary to explain the larger scatter of $N_{NH_2CHO}/N_{CH_3OH}$ compared with other ratios.

Third-generation (3G) gravitational-wave (GW) detectors like the Einstein Telescope (ET) will observe binary black hole (BBH) mergers at redshifts up to $z\sim 100$. However, unequivocal determination of the origin of high-redshift sources will remain uncertain, due to the low signal-to-noise ratio (SNR) and poor estimate of their luminosity distance. This study proposes a machine learning approach to infer the origins of high-redshift BBHs, specifically differentiating those arising from Population III (Pop. III) stars - likely the first progenitors of stellar-born BBH mergers in the Universe - and those originated from Population I-II (Pop. I-II) stars. We have considered a wide range of state-of-the-art models encompassing current uncertainties on Pop. III BBH mergers. We then estimate parameter errors of detected sources with ET using the Fisher-information-matrix formalism, followed by classification using XGBoost, a machine learning algorithm based on decision trees. For a set of mock observed BBHs, we provide the probability that they belong to the Pop. III class while considering the parameter errors of each source. In our fiducial model, we accurately identify ~10% of detected BBHs originating from Pop. III stars with > 90% precision. Our study demonstrates how machine learning enables to achieve some pivotal aspects of ET science case by exploring the origin of individual high-redshift GW observations. We set the basis for further studies, which will integrate additional simulated populations and account for population modeling uncertainties

In this chapter, we review some of the interesting consequences that tilt between the spin axis of the black hole and angular momentum axis of the accretion disk can have on the dynamics, thermodynamics, and observational appearance of accreting systems, from precessing coronae and jets to standing nozzle shocks and quasi-periodic oscillations. We begin the chapter by examining some of the reasons tilted disks are interesting as well as present arguments for how ubiquitous they may be. We then review the existing simulation results in the literature, broadly dividing them into sections on thick disks, thin disks, and magnetically arrested disks (MADs). We finish by highlighting some of the phenomenology that is unique to tilted disk simulations and discuss how this may apply to observations.

Metallicity offers a unique window into the baryonic history of the cosmos, being instrumental in probing evolutionary processes in galaxies between different cosmic environments. We aim to quantify the contribution of these environments to the scatter in the mass-metallicity relation (MZR) of galaxies. By analysing the galaxy distribution within the cosmic skeleton of the Horizon Run 5 cosmological hydrodynamical simulation at redshift $z = 0.625$, computed using a careful calibration of the T-ReX filament finder, we identify galaxies within three main environments: nodes, filaments and voids. We also classify galaxies based on the dynamical state of the clusters and the length of the filaments in which they reside. We find that the cosmic environment significantly contributes to the scatter in the MZR; in particular, both the gas metallicity and its average relative standard deviation increase when considering denser large-scale environments. The difference in the average metallicity between galaxies within relaxed and unrelaxed clusters is $\approx 0.1 \text{ dex}$, with both populations displaying positive residuals, $\delta Z_{g}$, from the averaged MZR. Moreover, the difference in metallicity between node and void galaxies accounts for $\approx 0.14 \, \text{dex}$ in the scatter of the MZR at stellar mass $M_{\star} \approx 10^{9.35}\,\text{M}_{\odot}$. Finally, both the average [O/Fe] in the gas and the galaxy gas fraction decrease when moving to higher large-scale densities in the simulation, suggesting that the cores of cosmic environments host, on average, older and more massive galaxies, whose enrichment is affected by a larger number of Type Ia Supernova events.

The Gemini High-resolution Optical SpecTrograph (GHOST) is a new echelle spectrograph available on the Gemini-South telescope as of Semester 2024A. We present the first high resolution spectrum of the quasar J1449-1227 (redshift z_em=3.27) using data taken during the commissioning of GHOST. The observed quasar hosts an intervening iron-poor ([Fe/H] = -2.5) damped Lyman alpha (DLA) system at redshift z=2.904. Taking advantage of the high spectral resolving power of GHOST (R~55000), we are able to accurately model the metal absorption lines of the metal-poor DLA and find a supersolar [Si/Fe], suggesting the DLA gas is in an early stage of chemical enrichment. Using simple ionization models, we find that the large range in the C IV/Si IV column density ratio of individual components within the DLA's high ionization absorption profile can be reproduced by several metal-poor Lyman limit systems surrounding the low-ionization gas of the DLA. It is possible that this metal-poor DLA resides within a complex system of metal-poor galaxies or filaments with inflowing gas. The high spectral resolution, wavelength coverage and sensitivity of GHOST makes it an ideal spectrograph for characterizing the chemistry and kinematics of quasar absorption lines.

In the context of the hierarchical model of galaxy evolution, polar ring galaxies (PRGs) are considered the intermediate phase between ongoing mergers and quiescent galaxies. This study explores the globular cluster system (GCS) and its properties in the nearest PRG, NGC4262, serving as a pilot investigation to study GCS in nearby PRGs. We utilize wide and deep field observations of the CFHT as part of the NGVS to investigate the GCS of NGC4262. We presented the first optical image of NGC4262 with an optically faint ring component. The photometric analysis of the GCS displays a distinct color bimodality. We estimate the total number of GCs for NGC4262 to be 266$\pm$16 GCs with a specific frequency of 4.2$\pm$0.8 and a specific mass of 0.23$\pm$0.01, which is relatively high compared to other galaxies of similar mass and environmental conditions. The spatial and azimuthal distributions of subpopulations reveal strong evidence of previous interactions within the host galaxy. The color distribution of GCS in NGC4262 shows a gradient of -0.05$\pm$0.01 within 5.5\arcmin, supporting the notion of past interactions and evolutionary transitions. PRG NGC4262 conforms to the overall trend of the GCS mass with respect to the halo mass. Furthermore, our investigation of the global scaling relations between GCS and host galaxy parameters provides further support for the hypothesis that PRGs are an intermediate phase connecting ongoing mergers and quiescent galaxies.

We present an analysis of the HST COS spectrum of IZw1 aiming to probe the absorbing medium associated with the active galactic nucleus (AGN). We fitted the emission spectrum and performed spectral analysis of the identified absorption features to derive the corresponding ionic column densities and covering fractions of the associated outflows. We employed photoionisation modelling to constrain the total column density and the ionisation parameter of four detected kinematic components. By investigating the implications of the results together with the observed kinematic properties of both emission and absorption features, we derived constraints on the structure and geometry of the absorbing medium in the AGN environment. We find and characterise absorption line systems from outflowing ionised gas in four distinct kinematic components, located at -60, -280, -1950, and -2900 km/s with respect to the source rest frame. While the two slower outflows are consistent with a full covering of the underlying radiation source, the well-constrained doublet line ratios of the faster two, higher column density, outflows suggest partial covering, with a covering fraction of C_f~0.4. The faster outflows show also line-locking in the NV doublet, a signature of acceleration via line absorption. This makes IZw1 possibly the closest object that shows evidence for hosting line-driven winds. The observed -1950 km/s absorption is likely due to the same gas as an X-ray warm absorber. Furthermore, the behaviour in UV and X-ray bands implies that this outflow has a clumpy structure. We find that the highly asymmetric broad emission lines in IZw1, indicative of a collimated, outflowing broad line region, are covered by the absorbing gas. Finally, the strongest UV--X-ray absorber may be connected to some of the blueshifted line emission, indicative of a more spatially extended structure of this ionised medium.

Neutron star mergers (NSMs) produce copious amounts of heavy r-process elements after a time delayed inspiral process. Once NSMs are present in a galaxy, r-process elements, such as Eu, are expected to significantly increase with time. Yet, there has been limited observational data in support of Eu increasing within Local Group galaxies. We have obtained high-resolution Magellan/MIKE observations of 43 metal-poor stars in the Gaia-Sausage/Enceladus tidally disrupted galaxy with $-2.5 < \rm{[Fe/H]} < -1$. For the first time, we find a clear rise in [Eu/Mg] with increasing [Mg/H] within one galaxy. We use a chemical evolution model to study how such a rise can result from the interplay of prompt and delayed r-process enrichment events. Delayed r-process sources are required to explain the rise and subsequent leveling off of [Eu/Mg] in this disrupted galaxy. However, the rise may be explained by delayed r-process sources with either short ($\sim 10$ Myr) or long ($\sim 500$ Myr) minimum delay times. Future studies on the nature of r-process sources and their enrichment processes in the GSE will require additional stars in the GSE at even lower metallicities than the present study.

We analyze a tilt instability of the orbit of an outer planet in a two planet circumbinary system that we recently reported. The binary is on an eccentric orbit and the inner circumbinary planet is on a circular polar orbit that causes the the binary to undergo apsidal precession. The outer circumbinary planet is initially on a circular or eccentric orbit that is coplanar with respect to the binary. We apply a Hamiltonian in quadrupole order of the binary potential to show that the tilt instability is the result of a secular resonance in which the apsidal precession rate of the binary matches the nodal precession rate of the outer planet. Resonance is possible because the polar inner planet causes the apsidal precession of the binary to be retrograde. The outer planet periodically undergoes large tilt oscillations for which we analytically determine the initial evolution and maximum inclination. Following a typically relatively short adjustment phase, the tilt grows exponentially in time at a characteristic rate that is of order the absolute value of the binary apsidal precession rate. The analytic results agree well with numerical simulations. This instability is analogous to the Kozai-Lidov instability, but applied to a circumbinary object. The instability fails to operate if the binary mass ratio is too extreme. The instability occurs even if the outer planet is instead an object of stellar mass and involves tilt oscillations of the inner binary.

With more than 5500 detected exoplanets, the search for life is entering a new era. Using life on Earth as our guide, we look beyond green landscapes to expand our ability to detect signs of surface life on other worlds. While oxygenic photosynthesis gives rise to modern green landscapes, bacteriochlorophyll-based anoxygenic phototrophs can also colour their habitats and could dominate a much wider range of environments on Earth-like exoplanets. Here, we characterize the reflectance spectra of a collection of purple sulfur and purple non-sulfur bacteria from a variety of anoxic and oxic environments. We present models for Earth-like planets where purple bacteria dominate the surface and show the impact of their signatures on the reflectance spectra of terrestrial exoplanets. Our research provides a new resource to guide the detection of purple bacteria and improves our chances of detecting life on exoplanets with upcoming telescopes. Our biological pigment data base for purple bacteria and the high-resolution spectra of Earth-like planets, including ocean worlds, snowball planets, frozen worlds, and Earth analogues, are available online, providing a tool for modellers and observers to train retrieval algorithms, optimize search strategies, and inform models of Earth-like planets, where purple is the new green.

Modern wide field radio surveys typically detect millions of objects. Techniques based on machine learning are proving to be useful for classifying large numbers of objects. The self-organizing map (SOM) is an unsupervised machine learning algorithm that projects a many-dimensional dataset onto a two- or three-dimensional lattice of neurons. This dimensionality reduction allows the user to visualize common features of the data better and develop algorithms for classifying objects that are not otherwise possible with large datasets. To this aim, we use the PINK implementation of a SOM. PINK incorporates rotation and flipping invariance so that the SOM algorithm may be applied to astronomical images. In this cookbook we provide instructions for working with PINK, including preprocessing the input images, training the model, and offering lessons learned through experimentation. The problem of imbalanced classes can be improved by careful selection of the training sample and increasing the number of neurons in the SOM (chosen by the user). Because PINK is not scale-invariant, structure can be smeared in the neurons. This can also be improved by increasing the number of neurons in the SOM. We also introduce pyink, a Python package used to read and write PINK binary files, assist in common preprocessing operations, perform standard analyses, visualize the SOM and preprocessed images, and create image-based annotations using a graphical interface. A tutorial is also provided to guide the user through the entire process. We present an application of PINK to VLA Sky Survey (VLASS) images. We demonstrate that the PINK is generally able to group VLASS sources with similar morphology together. We use the results of PINK to estimate the probability that a given source in the VLASS QuickLook Catalogue is actually due to sidelobe contamination.

We investigate the relation between galaxy structure and star formation rate (SFR) in a sample of $\sim2.9\times10^{4}$ central galaxies with $z<0.0674$ and axial ratios $b/a>0.5$. The star-forming main sequence (SFMS) shows a bend around the stellar mass of $M_\ast\leq{}M_c=2\times10^{10}{}M_{\odot}$. At $M_\ast\leq{}M_c$ the SFMS follows a power-law $\text{SFR}\propto{}M_\ast^{0.85}$, while at higher masses it flattens. $M_c$ corresponds to a dark matter halo mass of $M_\text{vir}\sim{}10^{11.8}M_{\odot}$ where virial shocks occurs. Some galaxy structure (e.g., half-light radius, $R_e$) exhibits a non-monotonic dependence across the SFMS at a fixed $M_\ast$. We find $\text{SFR}\propto{R_e^{-0.28}}$ at fixed $M_\ast$, consistent with the global Kennicutt-Schmidt (KS) law. This finding suggests that galaxy sizes contribute to the scatter of the SFMS. However, at $M_\ast>M_c$ the relationship between SFR and $R_e$ diminishes. Low-mass galaxies above the mean of the SFMS have smaller radii, exhibit compact and centrally concentrated profiles resembling green valley (GV) and quiescent galaxies at the same mass, and have higher $M_{\text{H}_2}/M_\text{HI}$. Conversely, those below the SFMS exhibit larger radii, lower densities, have no GV or quiescent counterparts at their mass and have lower $M_{\text{H}_2}/M_\text{HI}$. The above data suggest two pathways for quenching low-mass galaxies, $M_\ast\leq{}M_c$: a fast one that changes the morphology on the SFMS and a slow one that does not. Above $M_c$, galaxies below the SFMS resemble GV and quiescent galaxies structurally, implying that they undergo a structural transformation already within the SFMS. For these massive galaxies, CG are strongly bimodal, with SFMS galaxies exhibiting negative color gradients, suggesting most star formation occurs in their outskirts, maintaining them within the SFMS.

We use the magnetic field components measured by Zeeman Doppler imaging (ZDI) to calculate the stellar surface force and torque due to magnetic stresses for the fast rotators $\sigma$ Ori E, 36 Lyn and CU Vir, and the slow rotator $\tau$ Sco. If we assume the stars have spherical photospheres, the estimated torques give spin down time scales no larger than $7 \times 10^5$ yr. For $\sigma$ Ori E, the predicted spin down time scale, $\simeq 6000$ yr, is much less than the observationally measured time scale of $\simeq 10^6$ yr. However, for CU Vir, we find that the spin down time scale from its ZDI map is $7 \times 10^5$ yr in good agreement with its average rate of spin down from 1960 to 2010. With the exception of $\tau$ Sco, the net force due to magnetic stresses at the stellar surface are large compared to the surface-integrated pressure. We discuss possible reasons for the large values of the forces (and torques), and suggest that the likely explanation is that rotation and the magnetic stresses create significant departures from spherical symmetry.

Massive stars (those larger than 8 solar masses at formation) have radiative envelopes that cannot sustain a dynamo, the mechanism that produces magnetic fields in lower-mass stars. Despite this, approximately 7\% of massive stars have observed magnetic fields, the origin of which is debated. We used multi-epoch interferometric and spectroscopic observations to characterize HD 148937, a binary system of two massive stars. We found that only one star is magnetic and that it appears younger than its companion. The system properties and a surrounding bipolar nebula can be reproduced with a model in which two stars merged (in a previous triple system) to produce the magnetic massive star. Our results provide observational evidence that magnetic fields form in at least some massive stars through stellar mergers.

Sgr A lobes are a pair of 15-pc-sized bipolar bubbles with co-aligned major axes perpendicular to the Galactic plane found in X-ray and radio observations of the Galactic center (GC). Their elusive origin is a vital ingredient in understanding the ongoing high energy processes at the GC. Here we perform a suite of hydrodynamic simulations to explore the tidal disruption event (TDE) outflow model for the origin of the Sgr A lobes. By following the outflow evolution in the circumnuclear medium, we find that TDE outflows naturally produce bipolar lobes delimited by forward shocks, and the dense postshock shell contributes significantly to the lobe's X-ray emission. Our fiducial run reproduces the morphology, density, temperature, and X-ray surface brightness distribution of the observed Sgr A lobes reasonably well. The current lobe age is ~3300 yr. Our model further predicts that the uplifted wake flow breaks through the ejecta-induced shock front, producing a shock-enclosed head region which is relatively dim in X-rays compared to other lobe regions. Both the dim head region and the predicted limb-brightening feature of the lobes are slightly inconsistent with current X-ray observations, and may be used to further test our model with more sensitive future X-ray observations. While light narrow jets and massive wide winds from TDE events usually do not reproduce the observed oval-shaped morphology of the lobes, the TDE outflow in our fiducial run is massive and yet narrow. Whether it is a jet or wind remains uncertain and future simulations covering both the outflow acceleration region and its pc-scale evolution would be very helpful in determining whether the Sgr A lobes indeed originate from a pair of TDE jets or winds.

Millimeter and sub-millimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 $\mu$m and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can exhibit different correlations between observing wavelength and polarization morphology and fraction. We explore possible origins for these differences in polarization, including differences in dust populations and protostar properties. For RY Tau and MWC 480, which are consistent with scattering at both wavelengths, we present models of the scattering polarization from several dust grain size distributions. These models aim to reproduce two features of the observational results for these disks: (1) both disks have an observable degree of polarization at both wavelengths and (2) the polarization fraction is higher at 3 mm than at 870 $\mu$m in the centers of the disks. For both disks, these features can be reproduced by a power-law distribution of spherical dust grains with a maximum radius of 200 $\mu$m and high optical depth. In MWC 480, we can also reproduce features (1) and (2) with a model containing large grains ($a_{max}$ = 490 $\mu$m ) near the disk midplane and small grains ($a_{max}$ = 140 $\mu$m) above and below the midplane.

Active galactic nuclei (AGNs) are known to exhibit optical/UV variability and most of them can be well modeled by the damped random walks. Physical processes that are not related to the accretion disk, such as tidal disruption events (TDE) or moving foreground dusty clouds, can cause flare-like and eclipse-like features in the optical light curve. Both long-term and high-cadence monitoring are needed to identify such features. By combining the Sloan Digital Sky Survey (SDSS), Panoramic Survey Telescope, and Rapid Response System (Pan-STARRS) with the Zwicky Transient Facility (ZTF) survey, we are able to identify a rare sample (11) out of the SDSS quasar catalog (around 83, 000). These quasars exhibit more or less constant brightness but show rapid optical variation in the ZTF DR2 epochs. To investigate the possible origins of these flare/eclipse-like variabilities, we propose the second epoch spectroscopic observations with the Gran Telescopio CANARIAS (GTC). We find that the change in accretion rate plays a significant role in these quasar variabilities. Among them, we identify two Changing-Look Active Galactic Nuclei (CL-AGN) candidates: SDSS J1427+2930 and SDSS J1420+3757. The luminosity change of the former may be caused by the enhanced SMBH accretion or the tidal disruption event, while the latter is more related to the change in the accretion rate.

We report measurements characterizing the performance of a kinetic inductance detector array designed for a wavelength of 25 microns and very low optical background level suitable for applications such as a far-infrared instrument on a cryogenically cooled space telescope. In a pulse counting mode of operation at low optical flux, the detectors can resolve individual 25-micron photons. In an integrating mode, the detectors remain photon noise limited over more than six orders of magnitude in absorbed power from 70 zW to 200 fW, with a limiting NEP of 4.6 x 10^-20 W/rtHz at 1 Hz. In addition, the detectors are highly stable with flat power spectra under optical load down to 1 mHz. Operational parameters of the detector are determined including the efficiency of conversion of the incident optical power into quasiparticles in the aluminum absorbing element and the quasiparticle self-recombination constant.

The S-shaped swing of the linear polarization position angle (PPA) observed in many pulsars can be interpreted by the rotating vector model (RVM). However, efforts to fit the RVM for a large sample of pulsars observed with the MeerKAT telescope as a part of the Thousand-Pulsar-Array (TPA) programme, only succeeded for about half the cases. High time-resolution studies suggest that the failed cases arise due to the presence of orthogonal polarization modes, or highly disordered distribution of PPA points. One such example is PSR~J1645-0317. Recently it has been shown that the RVM can be recovered in this pulsar by using only time samples which are greater than 80% linearly polarized. In this work we test this novel approach on the brightest 249 pulsars from the TPA sample, of which 177 yield sufficient highly polarized samples to be amenable to our method. Remarkably, only 9 of these pulsars (5%) now fail to fit the RVM as opposed to 59% from the original analysis. This result favours the paradigm that the underlying mechanism is coherent curvature radiation.

This paper introduces a novel dynamical model, building upon the existing dynamical model for Deimos in the current numerical ephemerides, which only encompasses the simple libration effects of Deimos. The study comprehensively incorporates the rotational dynamics of Deimos influenced by the torque exerted by the major celestial bodies (Mars, the Sun) in the solar system within the inertial space. Consequently, a full dynamical model is formulated to account for the complete coupling between the rotation and orbit of Deimos. Simultaneously, employing precision orbit determination methods used for artificial satellites, we develop an adjustment model for fitting data to the complete model. The 12-order Adams--Bashforth--Moulton (ABM) integration algorithm is employed to synchronously integrate the 12 state variables of the full model to obtain the orbit of Deimos. Numerical simulation results indicate that the full dynamical model and adjustment model are stable and reliable. Compared to the simple model, the polar axis of Deimos in the inertial space exhibits a more complex oscillation in the full model. This work further advances the current dynamical model for Deimos and establishes the foundational model for the generation of a new set of precise numerical ephemerides for Deimos.

In a recent study, Banerjee et al. (2021) produced an atlas of all major emission lines found in a large sample of 115 Galactic field Be stars using the 2-m Himalayan Chandra Telescope (HCT) facility located at Ladakh, India. This paper presents our further exploration of these stars to estimate the electron density in their discs. Our study using Balmer decrement values indicate that their discs are generally optically thick in nature with electron density (n_e) in their circumstellar envelopes (CEs) being in excess of 10^13 cm^-3 for around 65% of the stars. For another 19% stars, the average n_e in their discs probably range between 10^12 cm^-3 and 10^13 cm^-3. We noticed that the nature of the H{\alpha} and H\b{eta} line profiles might not influence the observed Balmer decrement values (i.e. D_34 and D_54) of the sample of stars. Interestingly, we also found that around 50% of the Be stars displaying D_34 greater than 2.7 are of earlier spectral types, i.e. within B0 -B3.

Context. Dust polarization observations of the massive protocluster G31.41+0.31 carried out at ~1'' (~3750 au) resolution with the SMA at 870 microm have revealed one of the clearest examples to date of an hourglass-shaped magnetic field morphology in the high-massregime. Additionally, ~0.24'' (~900 au) resolution observations with ALMA at 1.3 mm have confirmed these results. The next step is to investigate whether the magnetic field maintains its hourglass-shaped morphology down to circumstellar scales. Aims. To study the magnetic field morphology toward the four (proto)stars A, B, C, and D contained in G31.41+0.31 and examine whether the self-similarity observed at core scales (1'' and 0.24'' resolution) still holds at circumstellar scales, we carried out ALMA observations of the polarized dust continuum emission at 1.3 mm and 3.1 mm at an angular resolution of ~0.068'' (~250 au), sufficient to resolve the envelope emission of the embedded protostars. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6) and 97.5 GHz (Band 3) with a synthesized beam of 0.072'' x 0.064''. We carried out polarization observations at two different wavelengths to confirm that the polarization traces magnetically aligned dust grains and is not due to dust self-scattering. Results. The polarized emission and the direction of the magnetic field obtained at the two wavelengths are basically the same, except for an area between the embedded sources C and B. In such an area, the emission at 1.3 mm could be optically thick and affected by dichroic extinction. ...

We estimate the coronal density of Capella using the O VII and Fe XVII line systems in the soft X-ray regime that have been observed over the course of the Chandra mission. Our analysis combines measures of error due to uncertainty in the underlying atomic data with statistical errors in the Chandra data to derive meaningful overall uncertainties on the plasma density of the coronae of Capella. We consider two Bayesian frameworks. First, the so-called pragmatic-Bayesian approach considers the atomic data and their uncertainties as fully specified and uncorrectable. The fully-Bayesian approach, on the other hand, allows the observed spectral data to update the atomic data and their uncertainties, thereby reducing the overall errors on the inferred parameters. To incorporate atomic data uncertainties, we obtain a set of atomic data replicates, the distribution of which captures their uncertainty. A principal component analysis of these replicates allows us to represent the atomic uncertainty with a lower-dimensional multivariate Gaussian distribution. A $t$-distribution approximation of the uncertainties of a subset of plasma parameters including a priori temperature information, obtained from the temperature-sensitive-only Fe XVII spectral line analysis, is carried forward into the density- and temperature-sensitive O VII spectral line analysis. Markov Chain Monte Carlo based model fitting is implemented including Multi-step Monte Carlo Gibbs Sampler and Hamiltonian Monte Carlo. Our analysis recovers an isothermally approximated coronal plasma temperature of $\approx$5 MK and a coronal plasma density of $\approx$10$^{10}$ cm$^{-3}$, with uncertainties of 0.1 and 0.2 dex respectively.

From the IllustrisTNG-50 simulation, a sample of 836 central disk galaxies with tiny stellar halos is chosen to study the inherent evolution of galaxies driven by nature. These galaxies are classified as compact, normal, or extended by referencing their locations on the mass-size ($M_\star-R_{\rm 1/2}$) diagram. This research demonstrates the distinctive evolutionary pathways of galaxies with different sizes in IllustrisTNG simulations, primarily driven by nature. It is confirmed that disk galaxies inherit the angular momentum of their parent dark matter halos. More compact galaxies form earlier within halos possessing lower specific angular momentum through heightened star formation during the early phase at redshifts above 2. During the later phase, the size of extended galaxies experiences more pronounced growth by accreting gas with high angular momentum. Additionally, we reveal that many key characteristics of galaxies are linked to their mass and size: (1) compact galaxies tend to exhibit higher metal content, proportional to the potential well $\frac{M_\star}{R_{\rm 1/2}}$, (2) compact galaxies host more massive bulges and black holes, and higher central concentration. Furthermore, our analysis indicates that galaxies of all types continue to actively engage in star formation, with no evident signs of quenching attributed to their varying sizes and angular momenta.

Direct atmospheric retrievals of exoplanets at high-resolution have recently allowed for a more detailed characterisation of their chemistry and dynamics from the ground. By monitoring the longitudinal distribution of species, as well as the varying vertical temperature structure and dynamics between the limbs of WASP-76b, across multiple transits, we aim to enhance our understanding of the 3D nature and chemical/dynamical evolution of such objects over the timescales of months/years. We present retrievals of three VLT/ESPRESSO observations of the ultra-hot Jupiter WASP-76b, including one not yet reported in the literature, from which we constrain the atmospheric abundances, vertical temperature structure, and atmospheric dynamics for both the leading and trailing limbs of the atmosphere separately, via novel rotational broadening kernels. We confirm the presence of VO recently reported in the atmosphere of WASP-76b. We find a uniform longitudinal distribution of Fe and Mg across the limbs of the atmosphere, for each of our transits, consistent with previous works as well as with stellar values. We constrain substellar Na/Fe and Cr/Fe ratios across each of our transits, consistent with previous studies of WASP-76b. Where constrained, V/Fe and VO/Fe ratios were also found to be broadly consistent between the limbs of the atmosphere for each of the transits, as well as with previous studies. However, for two of the transits, both V and VO were unconstrained in the leading limb, suggesting possible depletion due to recombination and condensation. The consistency of our constraints across multiple high-resolution observations, as well as with previous studies using varying modelling/retrieval frameworks and/or instruments, affirms the efficacy of high-resolution ground-based retrievals of exoplanetary atmospheres.

In this paper, we apply the machine learning clustering algorithm Density Based Spatial Clustering of Applications with Noise (DBSCAN) to study the membership of stars in twelve open clusters (NGC~2264, NGC~2682, NGC~2244, NGC~3293, NGC~6913, NGC~7142, IC~1805, NGC~6231, NGC~2243, NGC 6451, NGC 6005 and NGC 6583) based on Gaia DR3 Data. This sample of clusters spans a variety of parameters like age, metallicity, distance, extinction and a wide parameter space in proper motions and parallaxes. We obtain reliable cluster members using DBSCAN as faint as $G \sim 20$ mag and also in the outer regions of clusters. With our revised membership list, we plot color-magnitude diagrams and we obtain cluster parameters for our sample using ASteCA and compare it with the catalog values. We also validate our membership sample by spectroscopic data from APOGEE and GALAH for the available data. This paper demonstrates the effectiveness of DBSCAN in membership determination of clusters.

Due to the high quality constraints available for the Sun, we can carry out combined analyses using neutrino, spectroscopic and helioseismic observations. Such studies lay the ground for future improvements of key physical components of solar and stellar models, such as the equation of state, radiative opacities or prescriptions for macroscopic transport processes of chemicals which are then used to study other stars in the Universe. We study the existing degeneracies in solar models using the recent high-metallicity spectroscopic abundances by comparing them to helioseismic and neutrino data and discuss how their properties are impacted by changes in various physical ingredients. We carry out a detailed study of solar models computed with a high-metallicity composition from the literature based on averaged-3D models that was claimed to solve the solar problem. The properties of the solar models are significantly affected by using the recent OPLIB opacities and the inclusion of macroscopic transport. The properties of the standard solar models computed using the OPAL opacities are similar to those using the OP opacities. We show that a modifying the temperature gradient just below the base of the convective zone is required to erase the discrepancies in solar models, particularly in the presence of macroscopic mixing. This can be simulated by a local increase of opacity of a few percent. We conclude that the existing degeneracies and issues in solar modelling are not erased by an increase in the solar metallicity in contradiction to was suggested in recent papers. Therefore, standard solar models cannot be used as an argument for a high metallicity composition. While further work is required to improve solar models, we note that direct helioseismic inversions indicate a low metallicity in the convective envelope, in agreement with spectroscopic analyses based on full 3D models.

Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way.

While SNRs have been considered the most relevant Galactic CR accelerators for decades, CCSNe could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV gamma-rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. The multi-wavelength observation of SN 2023ixf, a Type II SN in the nearby galaxy M101, opens the possibility to constrain CR acceleration within a few days after the collapse of the RSG stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV. We performed a maximum-likelihood analysis of gamma-ray data from the Fermi Large Area Telescope up to one month after the SN explosion. We searched for high-energy emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration. We do not find significant gamma-ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN. Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. A more sophisticated, inhomogeneous characterisation of the shock and the progenitor's environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times.

Much has been written about the representation of noisy linear 2-ports. Here we present a theory of noisy $N$-ports. We show how in the general case there are $(2N)!/(N!)^2$ equivalent representations and give the transformations relating them. We also discuss singular cases in which some of the transformations are not possible as well as how to measure the noise properties of an $N$-port. This work is motivated by the REACH experiment to observe the global 21 cm signal for which modelling noise with exquisite precision is essential for a reliable calibration.

We analyze the magnetic evolution of solar active region (AR) NOAA 11476 that, between 9 and 10 May 2012, produced a series of surge-type eruptions accompanied by GOES X-ray class M flares. Using force-free models of the AR coronal structure and observations in several wavelengths, in previous works we studied the detailed evolution of those eruptions, relating them to the characteristic magnetic topology of the AR and reconstructing the involved reconnection scheme. We found that the eruptions were due to the ejection of minifilaments, which were recurrently ejected and reformed at the polarity inversion line of a bipole that emerged in the middle of the positive main AR magnetic polarity. The bipole was observed to rotate for several tens of hours before the events. In this article we analyze, for the full AR and the rotating bipole, the evolution of a series of magnetic parameters computed using the Helioseismic and Magnetic Imager (HMI) vector magnetograms. We combine this analysis with estimations of the injection of magnetic energy and helicity obtained using the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM) method that determines, from vector magnetograms, the affine velocity field constrained by the induction equation. From our results, we conclude that the bipole rotation was the main driver that provided the magnetic energy and helicity involved in the minifilament destabilizations and ejections. The results also suggest that the observed rotation is probably due to the emergence of a kinked magnetic flux rope with negative writhe helicity.

Understanding how giant and terrestrial planets form and evolve, what is their internal structure and that of their atmosphere, represents one of the major challenges of modern astronomy, which is directly connected to the ultimate search for life at the horizon 2040. However, several astrophysical, biological and technological obstacles must be overcome. From the astrophysical point of view, it is indeed crucial to understand the mechanisms of formation and evolution of giant planets, including planet and disk interactions, which will completely sculpt the planetary architectures and thus dominate the formation of terrestrial planets, especially in regions around the host star capable of supporting life. It is also important to develop dedicated instrumentation and techniques to study in their totality the population of giant and terrestrial planets, but also to reveal in the near future the first biological clues of life in the atmospheres of terrestrial planets. In that perspective, direct imaging from ground-based observatories or in space is playing a central role in concert with other observing techniques. This review introduces the genesis of this observing technique, the main instrumental innovation and challenges, stellar targets and surveys, to then present the main results obtained so far about the physics and the mechanisms of formation and evolution of young giant planets and planetary system architectures. I will then present the exciting perspectives offered by the upcoming generation of planet imagers about to come online, particularly on the future extremely large telescopes. On the timescale of a human Life, we may well be witnessing the first discovery of an exoplanet and the first detection of indices of life in the atmosphere of a nearby exo-Earth!

We propose a new model of scalar field dark matter interacting with dark energy. Adopting a fluid description of the dark matter field in the regime of rapid oscillations, we find that the equation of state for dark matter is non-zero and even becomes increasingly negative at late times during dark energy domination. Furthermore, the speed of sound of dark matter is non-vanishing at all length scales, and a non-adiabatic pressure contribution arises. The results indicate that there are still unexplored possible interactions within the dark sector that lead to novel background effects and can impact structure formation processes.

Decay and dispersal of the tilted Bipolar Magnetic Regions (BMRs) on the solar surface are observed to produce large-scale poloidal field, which acts as the seed for the toroidal field and, thus, the next sunspot cycle. However, various properties of BMR, namely, the tilt, time delay between successive emergences, location, and flux, all have irregular variations. Previous studies show that these variations can lead to changes in the polar field. In this study, we first demonstrate that our 3D kinematic dynamo model, STABLE, reproduces the robust feature of the surface flux transport (SFT) model, namely the variation of the generated dipole moment with the latitude of the BMR position. Using STABLE in both SFT and dynamo modes, we perform simulations by varying the individual properties of BMR and keeping their distributions the same in all the cycles as inspired by the observations. We find that randomness due to the distribution in either the time delay or the BMR latitude produces negligible variation in the polar field and the solar cycle. However, randomness due to BMR flux distribution produces substantial effects, while the scatter in the tilt around Joy law produces the largest variation. Our comparative analyses suggest that the scatter of BMR tilt around Joy law is the major cause of variation in the solar cycle. Furthermore, our simulations also show that the magnetic field-dependent time delay of BMR emergence produces more realistic features of the magnetic cycle, consistent with observation.

We quantify the stellar rotation of galaxies by computing the $\lambda_{R}$ parameter, a proxy for the stellar angular momentum in a sample of 106 galaxies with redshift 0.1 $<$ z $<$ 0.8 and stellar masses from $\sim$10$^{7.5}$ to 10$^{11.8}$ M$_{\odot}$. The sample is located in the CANDELS/GOODS-S and COSMOS fields, and it was observed by various MUSE surveys. We create stellar velocity and velocity dispersion maps using a full-spectrum fitting technique, covering spatially $\sim$2$R_{e}$ for the galaxies. We study the impact of the atmospheric seeing on the spin parameter and apply corrections when pertinent. Through the analysis of the $\lambda_{R}-\epsilon$ diagram, we notice that the fraction of round and massive galaxies increases with redshift. We lack galaxies with $\lambda_{R}$ < 0.1 in the sample and we find only one potential, but uncertain, low-mass slow rotator at z $\sim0.3$. Moreover, we do not see an evident evolution or trend in the stellar angular momentum with redshift. We characterize the sample environment using two indicators: a local estimator based on the Voronoi tesselation method, and a global estimator derived by the use of the Friends-of-Friends algorithm. We find no correlation between the environment and $\lambda_{R}$ given that we are not probing dense regions or massive galaxy structures. We also analyze the kinematic maps of the sample finding that about 40$\%$ of galaxies are consistent with being regular rotators, having rotating stellar discs with flat velocity dispersion maps, while $\sim20\%$ of galaxies have complex velocity maps and can be identified as non-regular rotators in spite of their $\lambda_{R}$ values. For the remaining galaxies the classification is uncertain. As we lack galaxies with $\lambda_{R}$< 0.1, we are not able to identify when galaxies become slow rotators within the surveyed environments, area and redshift range.

We study metal enrichment originating from stellar wind and supernovae in low-metallicity clouds by performing three-dimensional radiation hydrodynamics simulations. We find that metals ejected from stellar wind are accumulated, leading to subsequent star formation in the nitrogen-enriched gas. During this early phase, the N/O ratios are similar to observed nitrogen-enriched galaxies (${\rm [N/O]}\gtrsim0.5$). Then, once supernovae occur, the N/O ratios decrease significantly. If the duration of star formation is comparable to the timescale of SNe, the mass fraction of nitrogen-enriched stars reaches half the mass of star clusters. We suggest that the mass of the star cluster needs to exceed $\sim 10^6~M_{\odot}$ to have multiple populations due to stellar wind, considering the condition for massive star cluster formation and the timescales of stellar evolution.

This study presents a comprehensive analysis of the spectral properties of 886 pulsars across a wide frequency range from 20MHz to 343.5GHz, including a total of 86 millisecond pulsars. The majority of the pulsars exhibit power-law behavior in their spectra, although some exceptions are observed. Five different spectral models, namely simple power-law, broken power-law, low-frequency turn-over, high-frequency cut-off, and double turn-over, were employed to explore the spectral behaviors. The average spectral index for pulsars modeled with a simple power-law is found to be -1.64 +/-0.80, consistent with previous studies. Additionally, significant correlations between the spectral index and characteristic parameters are observed particularly in millisecond pulsars, while no strong correlation is observed in normal pulsars. Different models show variations in the most influential characteristic parameters associated with the spectral index, indicating diverse dominant radiation mechanisms in millisecond pulsars.Finally, this study identifies 22 pulsars of the Gigahertz-peaked Spectra (GPS) type for the first time based on the Akaike information criterion.

Recent discoveries of copious amounts of dust in quiescent galaxies (QGs) at high redshifts ($z\gtrsim 1-2$) challenge the conventional view that these objects have poor interstellar medium (ISM) in proportion to their stellar mass. We use the SIMBA cosmological simulation to explore the evolution of dust and cold gas content in QGs in relation to the quenching processes affecting them. We track the changes in the ISM dust abundance across the evolutionary history of QGs identified at $0 \lesssim z \lesssim2$ in the field and cluster environments. The QGs quench via diverse pathways, both rapid and slow, and exhibit a wide range of times elapsed between the quenching event and cold gas removal (from $\sim650$ Myr to $\sim8$ Gyr). We find that quenching modes attributed to the feedback from active galactic nuclei (AGN) do not affect dust and cold gas within the same timescales. Remarkably, QGs may replenish their dust content in the quenched phase primarily due to internal processes and marginally by external factors such as minor mergers. The key mechanism for re-formation of dust is prolonged grain growth on gas-phase metals, it is effective within $\sim100$ Myr after the quenching event, and rapidly increases the dust-to-gas mass ratio in QGs above the standard values ($\delta_{\rm DGR}\gtrsim1/100$). As a result, despite heavily depleted cold gas reservoirs, roughly half of QGs maintain little evolution in their ISM dust with stellar age within the first 2 Gyr following the quenching. Overall, we predict that relatively dusty QGs ($M_{\rm dust}/M_{\star}\gtrsim10^{-3}-10^{-4}$) arise from both fast and slow quenchers, and are prevalent in systems of intermediate and low stellar masses ($9<\log(M_{\star}/M_{\odot})<10.5$). This prediction poses an immediate quest for observational synergy between e.g., James Webb Space Telescope (JWST) and the Atacama Large Millimeter Array (ALMA).

Torsional Alfvén waves do not produce any intensity variation and are, therefore, challenging to observe with imaging instruments. Previously, Alfvén wave observations were reported throughout all the layers of the solar atmosphere using spectral imaging. We present an observation of a torsional Alfvén wave detected in an inverted y-shape structure observed with the HRIEUV telescope of the EUI instrument onboard Solar Orbiter in its 174 Å channel. The feature consists of two footpoints connected through short loops and a spine with a length of 30 Mm originating from one of the footpoints. In the current work, we also make use of the simultaneous observations from two other instruments onboard Solar Orbiter. The first one is PHI that is used to derive the magnetic configuration of the observed feature. The second one is SPICE that provided observations of intensity maps in different lines including Ne VIII and C III lines. We also address the issues of the SPICE point spread function and its influence on the Doppler maps via performed forward modeling analysis. The difference movie shows clear signatures of propagating rotational motions in the spine. Doppler maps obtained with SPICE show strong signal in the spine region. Under the assumption that the recovered point spread function is mostly correct, synthesized raster images confirm that this signal is predominantly physical. We conclude that the presented observations are compatible with an interpretation of either propagating torsional Alfvén waves in a low coronal structure or untwisting of a flux rope. This is the first time we see signatures of propagating torsional motion in corona as observed by the three instruments onboard Solar Orbiter.

In this {\em Letter}, we explore the formation of the mass-gap black hole-neutron star (mgBHNS) merger detected in gravitational wave (GW) event, i.e., GW230529, from the isolated binary evolution channel, and study potential signatures of its electromagnetic signals. By adopting the `delayed' supernova prescription and reasonable model realizations, our population synthesis simulation results can simultaneously match the inferred event rate densities of GW230529-like mgBHNS and total BHNS mergers, as well as the population distribution of the BH mass in BHNS mergers reported by the LIGO-Virgo-KAGRA Collaboration. Thus, we conclude that the recently-discovered mgBHNS merger, GW230529, can be explained through the isolated binary evolution channel. Considering the equation of states of AP4 and DD2, the probabilities that GW230529 can make tidal disruption are $12.8\%$ and $63.2\%$, respectively. If GW230529 is a disrupted event, the associated kilonova is predicted to have an apparent magnitude of $\sim23-24\,{\rm{mag}}$, and hence, can be detected by the present survey projects and LSST. Since GW230529 could be an off-axis event inferred from the GW observation, its associated gamma-ray burst (GRB) might be too dim to be observed by $\gamma$-ray detectors, interpreting the lack of GRB observations. The detection of GW230529 confirms the existence of mgBHNS mergers formed through the isolated binary evolution channel, suggesting that BHNS mergers are still likely to be multimessenger sources that emit GWs, GRBs, and kilonovae. Although mgBHNS mergers account for $\sim60\%$ cosmological BHNS population, we find that $\gtrsim90\%$ disrupted BHNS mergers are expected to originate from mgBHNS mergers.

A minimalist approach to the linear stability problem in fluid dynamics is developed that ensures efficiency by utilizing only the essential elements required to find the eigenvalues for given boundary conditions. It is shown that the problem is equivalent to a single first-order ordinary differential equation, and that studying the argument of the unknown complex function in the eigenvalue space is sufficient to find the dispersion relation. The method is applied to a model for relativistic magnetized astrophysical jets.

The capability of Cosmic Inflation to explain the latest Cosmic Microwave Background and Baryonic Acoustic Oscillation data is assessed by performing Bayesian model comparison within the landscape of nearly three-hundred models of single-field slow-roll inflation. We present the first Bayesian data analysis based on the third-order slow-roll primordial power spectra. In particular, the fourth Hubble-flow function $\epsilon_4$ remains unbounded while the third function verifies, at two-sigma, $\epsilon_{3}\in[-0.4,0.5]$, which is perfectly compatible with the slow-roll predictions for the running of the spectral index. We also observe some residual excess of $B$-modes within the BICEP/Keck data favoring, at a non-statistically significant level, non-vanishing primordial tensor modes: $\log(\epsilon_{1}) > -3.9$, at $68\%$ confidence level. Then, for 283 models of single-field inflation, we compute the Bayesian evidence, the Bayesian dimensionality and the marginalized posteriors of all the models' parameters, including the ones associated with the reheating era. The average information gain on the reheating parameter $R_\mathrm{reh}$ reaches $1.3 \pm 0.18$ bits, which is more than a factor two improvement compared to the first Planck data release. As such, inflationary model predictions cannot meet data accuracy without specifying, or marginalizing over, the reheating kinematics. We also find that more than $40\%$ of the scenarios are now strongly disfavored, which shows that the constraining power of cosmological data is winning against the increase of the number of proposed models. In addition, about $20\%$ of all models have evidences within the most probable region and are all favored according to the Jeffreys' scale of Bayesian evidences.

Fast X-ray Transients (FXTs) are extragalactic bursts of soft X-rays first identified >10 years ago. Since then, nearly 40 events have been discovered, although almost all of these have been recovered from archival Chandra and XMM-Newton data. To date, optical sky surveys and follow-up searches have not revealed any multi-wavelength counterparts. The Einstein Probe, launched in January 2024, has started surveying the sky in the soft X-ray regime (0.5-4 keV) and will rapidly increase the sample of FXTs discovered in real time. Here, we report the first discovery of both an optical and radio counterpart to an FXT, the fourth source publicly released by the Einstein Probe. We discovered a fast-fading optical transient within the 3 arcmin localisation radius of EP240315a with the all-sky optical survey ATLAS, and our follow-up Gemini spectrum provides a redshift, z=4.859+/-0.002. Furthermore, we uncovered a radio counterpart in the S-band (3.0 GHz) with the MeerKAT radio interferometer. The optical (rest-frame UV) and radio luminosities indicate the FXT most likely originates from either a long gamma-ray burst or a relativistic tidal disruption event. This may be a fortuitous early mission detection by the Einstein Probe or may signpost a mode of discovery for high-redshift, high-energy transients through soft X-ray surveys, combined with locating multi-wavelength counterparts.

One of the leading mechanisms invoked to explain the existence of the radius valley is atmospheric mass loss driven by X-ray and extreme-ultraviolet irradiation, with this process stripping the primordial envelopes of young, small planets to produce the observed bimodal distribution. We present an investigation into the TOI-431 and $\nu^2$ Lupi planetary systems, both of which host planets either side of the radius valley, to determine if their architectures are consistent with evolution by the XUV mechanism. With $\textit{XMM-Newton}$, we measure the current X-ray flux of each star, and see evidence for a stellar flare in the TOI-431 observations. We then simulate the evolution of all of the transiting planets across the two systems in response to the high-energy irradiation over their lifetimes. We use the measured X-ray fluxes as an anchor point for the XUV time evolution in our simulations, and employ several different models of estimating mass loss rates. While the simulations for TOI-431b encountered a problem with the initial calculated radii, we estimate a likely short ($\sim$ Myr) timespan for primordial envelope removal using reasonable assumptions for the initial planet. $\nu^2$ Lupi b is likely harder to strip, but is achieved in a moderate fraction of our simulations. None of our simulations stripped any of the lower density planets of their envelope, in line with prediction. We conclude that both systems are consistent with expectations for generation of the radius valley through XUV photoevaporation.

Microquasars are Black Hole X-ray binaries (BHXB) which can eject material in the form of a bipolar jet, similarly to quasars, but at much smaller scales. Their high-energy emission comes from an accretion disk (~ 1 keV) and from a hot "corona" near the black hole that up-scatters photons from the disk in the hard X-ray domain (1--100 keV). A high-energy component above 150 keV has been detected in bright sources and its precise origin is still unknown: it could come either from Compton scattering of disk photons on coronal relativistic non-thermal electrons (a.k.a hybrid Comptonization), or from the synchrotron emission from the very base of the compact jet. The measurement of polarization above 150 keV can provide valuable insights into the processes at play as we expect higher polarization fraction due to synchrotron emission from the jets (up to 70 % with a very ordered magnetic field). We use the INTEGRAL/IBIS telescope to measure the soft gamma-ray polarization of the Crab Nebula and the BHXB Swift J1727.8-1613.

$ $The first public 0.9-4.4 $\mu$m NIRCam images of the North Ecliptic Pole (NEP) Time Domain Field (TDF) uncovered many galaxies that display point-source features in their cores as seen in the longer wavelength filters. We visually identified a sample of 66 galaxies ($\sim$1 galaxy per arcmin$^2$) with point-like cores and fit their spectral energy distributions (SED)s using EAZY and CIGALE to characterize the sample's active galactic nucleus (AGN) and host galaxy parameters. Single-template fitting best fits $70\%$ of the sample with a Seyfert-blended SED. With CIGALE we compute the median fractional AGN contribution to the 0.1-30.0 $\mu$m flux to be $0.30\pm0.06$, and that $56\%$ of the 66 galaxies have star-formation rates in the starburst range whereas the remainder are near the star-formation main sequence. There are Very Large Array (VLA) 3 GHz detections for 24/66 galaxies, implying some combination of AGN emission and vigorous star formation. We present a novel sample selection procedure in tandem to the morphological sample selection based on objects' light profiles at 4.4 $\mu$m. This procedure identifies a parameter space which probes point-source features and automatically recovers our visual sample with minimal contamination from both stars and brighter galaxies. This procedure may be used in other extant and future NIRCam images to streamline the search for galaxies with point-like cores. The morphological approach to recognizing AGN is being resurrected by the James Webb Space Telescope (JWST) by virtue of its superb angular resolution at infrared wavelengths.

V583 Lyr is an extremely low mass ratio Algol-type binary with an orbital period of 11.2580 days. We determined an effective temperature of T_{eff1} = 9000 \pm 350 K from newly observed spectra, which might be an underestimate due to binary mass transfer. The binary mass ratio q = 0.1 \pm 0.004 and the orbital inclination i = 85.5° are determined based on the assumption that the secondary fills its Roche lobe and rotates synchronously. The radial velocity curve is obtained from time series spectra, allowing for improved estimation of stellar masses and radii: M1 = 3.56 \pm 0.5 Msun, R1 = 2.4 \pm 0.2 Rsun; and M2 = 0.36 \pm 0.02 Msun, R2= 6.9 \pm 0.4 Rsun. The variations in the double-peaked H_{\alpha} emission indicate the formation of a stable disk during mass transfer. V583 Lyr appears to be a post-mass-reversal system, according to the estimated mass transfer using O-C period analysis. Its orbital period is slowly increasing, from which the rate of mass accretion by the primary star is estimated to be dM1/dt = 3.384 \times10^{-8} Msun/yr. The pulsation analysis was conducted on the residuals of the light curve. The primary component was found to be a g-mode pulsating star with 26 frequencies extracted lower than 9 d^{-1}. The frequency groups and rotational splitting properties of the g-mode were studied in detail. This study provides compelling evidence for an accretion disk surrounding the g-mode pulsating primary.

Solar spicules are plasma jets observed in the interface region between the visible solar surface and the corona. At any given time, there are millions of spicules present all over the Sun. While various models attempt to elucidate their origin and characteristics, here, we consider the one driven by the magneto-convection undulations. The radiative magneto-hydrodynamical (rMHD) equations are solved using PENCIL CODE with a spatial resolution of 16 km using various magnetic field strengths. The obtained rMHD simulation data are investigated to unveil the various trends in spicular properties as function of the applied magnetic fields. The important outcome of this study is the finding of a consistent reduction in both the number density and the maximum height reached by spicules as magnetic field strength increases. We also use parabolic fitting on time-distance curves of spicules that are taller than $75^\mathrm{th}$ percentile in the distribution, in order to find a relation between deceleration of the spicule tip and the magnetic field strength. Our results offer insights into the response of solar spicules to magnetic field strength.

This study explores the metal enrichment signatures attributed to the first generation of stars (PopIII) in the Universe, focusing on the E-XQR-30 sample. We aim to identify traces of Pop III metal enrichment by analyzing neutral gas in the interstellar medium of primordial galaxies and their satellite clumps, detected in absorption. To chase the chemical signature of PopIII stars, we studied metal absorption systems in the E-XQR-30 sample, selected through the detection of the OI absorption line at 1302A. The OI line is a reliable tracer of HI and allowed us to overcome the challenges posed by the Lyman-$\alpha$ forest's increasing saturation at redshifts above $\sim5$ to identify Damped Lyman-$\alpha$ systems (DLA). We detected and analyzed 29 OI systems at $z\geq5.4$, differentiating between proximate DLAs (PDLA) and intervening DLAs. Voigt function fits were applied to obtain ionic column densities, and relative chemical abundances were determined for 28 systems. These were then compared with the predictions of theoretical models. Our findings expand the study of OI systems at $z\geq5.4$ fourfold. No systematic differences were observed in the average chemical abundances between PDLAs and intervening DLAs. The chemical abundances in our sample align with literature systems at $z>4.5$, suggesting a similar enrichment pattern for this class of absorption systems. A comparison between these DLA-analogues at $4.5<z<6.5$ with a sample of very metal-poor DLAs at $2<z<4.5$ shows in general similar average values for the relative abundances, with the exception of [C/O], [Si/Fe] and [Si/O] which are significantly larger for the high-$z$ sample. Furthermore, the dispersion of the measurements significantly increases in the high-redshift bin. This increase is predicted by the theoretical models and indicates a potential retention of PopIII signatures in the probed gas. (Abridged)

We use 3D computer modelling to investigate the timescales and radiative output from maser flares generated by the impact of shock-waves on astronomical unit-scale clouds in interstellar and star-forming regions, and in circumstellar regions in some circumstances. Physical conditions are derived from simple models of isothermal hydrodynamic (single-fluid) and C-type (ionic and neutral fluid) shock-waves, and based on the ortho-H$_2$O 22-GHz transition. Maser saturation is comprehensively included, and we find that the most saturated maser inversions are found predominantly in the shocked material. We study the effect on the intensity, flux density and duration of flares of the following parameters: the pre-shock level of saturation, the observer's viewpoint, and the shock speed. Our models are able to reproduce observed flare rise times of a few times 10 days, specific intensities of up to 10$^5$ times the saturation intensity and flux densities of order $100(R/d)^2$Jy from a source of radius $R$ astronomical units at a distance of $d$ kiloparsec. We found that flares from C-type shocks are approximately 5 times more likely to be seen by a randomly placed observer than flares from hydrodynamically shocked clouds of similar dimensions. We computed intrinsic beaming patterns of the maser emission, finding substantial extension of the pattern parallel to the shock front in the hydrodynamic models. Beaming solid angles for hydrodynamic models can be as small as $1.3\times 10^{-5}$sr, but are an order of magnitude larger for C-type models.

We present an analysis of rest-frame UV continuum slopes, $\beta$, using a sample of 1011 galaxies at $6.5<z<13$ from the EPOCHS photometric sample collated from the GTO PEARLS and public ERS/GTO/GO (JADES, CEERS, NGDEEP, GLASS) JWST NIRCam imaging across $178.9~\mathrm{arcmin}^2$ of unmasked blank sky. We correct our UV slopes for the photometric error coupling bias using $200,000$ power law SEDs for each $\beta=\{-1,-1.5,-2,-2.5,-3\}$ in each field, finding biases as large as $\Delta\beta\simeq-0.55$ for the lowest SNR galaxies in our sample. Additionally, we simulate the impact of rest-UV line emission (including Ly$\alpha$) and damped Ly$\alpha$ systems on our measured $\beta$, finding biases as large as $0.5-0.6$ for the most extreme systems. We find a decreasing trend with redshift of $\beta=-1.51\pm0.08-(0.097\pm0.010)\times z$, with potential evidence for Pop.~III stars or top-heavy initial mass functions (IMFs) in a subsample of 68 $\beta+\sigma_{\beta}<-2.8$ galaxies. At $z\simeq11.5$, we measure an extremely blue $\beta(M_{\mathrm{UV}}=-19)=-2.73\pm0.06$, deviating from simulations, indicative of low-metallicity galaxies with non-zero Lyman continuum escape fractions $f_{\mathrm{esc, LyC}}\gtrsim0$ and minimal dust content. The observed steepening of $\mathrm{d}\beta/\mathrm{d}\log_{10}(M_{\star}/\mathrm{M}_{\odot})$ from $0.22\pm0.02$ at $z=7$ to $0.81\pm0.13$ at $z=11.5$ implies that dust produced in core-collapse supernovae (SNe) at early times may be ejected via outflows from low mass galaxies. We also observe a flatter $\mathrm{d}\beta/\mathrm{d}M_{\mathrm{UV}}=0.03\pm0.02$ at $z=7$ and a shallower $\mathrm{d}\beta/\mathrm{d}\log_{10}(M_{\star} / \mathrm{M}_{\odot})$ at $z<11$ than seen by HST, unveiling a new population of low mass, faint, galaxies reddened by dust produced in the stellar winds of asymptotic giant branch (AGB) stars or carbon-rich Wolf-Rayet binaries.

We suggest that "Godzilla", an intriguing source in the lensed Sunburst galaxy at $z=2.37$, is a young super star cluster powering a compact nebula within gravitationally trapped stellar ejecta. Employing HST photometry and spectroscopy from MUSE and X-Shooter at VLT, we infer physical and chemical properties of the cluster and nebula, finding Godzilla is young (4-6 Myr), massive ($\sim 10^{6-7}M_\odot$), a stellar metallicity $Z \simeq 0.25Z_\odot$, and has the FUV component more compact than a few pc. The nebula gas is significantly enriched with N and He, indicating stellar wind material, and has highly elevated O relative to the sub-solar stellar metallicity, which indicates entrainment of CCSNe ejecta. The high gas density $n_{\rm e} \simeq 10^{7-8}{\rm cm}^{-3}$ implies a highly pressurized intracluster environment. We propose the high pressure is due to CCSN-driven supersonic turbulence in warm, self-shielding gas, which has accumulated in the cluster center after runaway radiative cooling and is dense enough to resist removal by CCSNe. The nebula gas shows sub-solar C/O, Ne/O and Si/O values, which may reflect the CCSN element yields for initial stellar masses $>40M_\odot$. A comparison to element yield synthesis models for young star clusters shows that the gas abundance pattern is consistent with complete retention and mixture of stellar winds and CCSNe ejecta until the inferred cluster age. The O and He enhancement we find may have implications for the formation of multiple stellar populations in globular clusters, as Godzilla likely has already formed second-generation stars prior to the onset of CCSNe and evolved star winds, in order not to contradict the non-observation of O and large He enhancement in second-generation stars.

Light curves serve as a valuable source of information on stellar formation and evolution. With the rapid advancement of machine learning techniques, it can be effectively processed to extract astronomical patterns and information. In this study, we present a comprehensive evaluation of deep-learning and large language model (LLM) based models for the automatic classification of variable star light curves, based on large datasets from the Kepler and K2 missions. Special emphasis is placed on Cepheids, RR Lyrae, and eclipsing binaries, examining the influence of observational cadence and phase distribution on classification precision. Employing AutoDL optimization, we achieve striking performance with the 1D-Convolution+BiLSTM architecture and the Swin Transformer, hitting accuracies of 94\% and 99\% correspondingly, with the latter demonstrating a notable 83\% accuracy in discerning the elusive Type II Cepheids-comprising merely 0.02\% of the total dataset.We unveil StarWhisper LightCurve (LC), an innovative Series comprising three LLM-based models: LLM, multimodal large language model (MLLM), and Large Audio Language Model (LALM). Each model is fine-tuned with strategic prompt engineering and customized training methods to explore the emergent abilities of these models for astronomical data. Remarkably, StarWhisper LC Series exhibit high accuracies around 90\%, significantly reducing the need for explicit feature engineering, thereby paving the way for streamlined parallel data processing and the progression of multifaceted multimodal models in astronomical applications. The study furnishes two detailed catalogs illustrating the impacts of phase and sampling intervals on deep learning classification accuracy, showing that a substantial decrease of up to 14\% in observation duration and 21\% in sampling points can be realized without compromising accuracy by more than 10\%.

We present a statistical analysis of the local, approximately 50-100 pc scale, H-alpha emission at the locations of recent (less than 125 years) supernovae (SNe) in nearby star-forming galaxies. Our sample consists of 32 SNe in 10 galaxies that are targets of the PHANGS-MUSE survey. We find that 41% (13/32) of these SNe occur coincident with a previously identified HII region. For comparison, HII regions cover 32% of the area within 1 kpc of any recent SN. Contrasting this local covering fraction with the fraction of SNe coincident with HII regions, we find a statistical excess of 7.6% +/- 8.7% of all SNe to be associated with HII regions. This increases to an excess of 19.2% +/- 10.4% when considering only core-collapse SNe. These estimates appear to be in good agreement with qualitative results from new, higher resolution HST H-alpha imaging, which also suggest many CCSNe detonate near but not in HII regions. Our results appear consistent with the expectation that only a modest fraction of stars explode during the first 5 Myr of the life of a stellar population, when H-alpha emission is expected to be bright. Of the HII region associated SNe, 8% (11/13) also have associated detected CO(2-1) emission, indicating the presence of molecular gas. The HII region associated SNe have typical Av extinctions approximately equal to 1 mag, consistent with a significant amount of pre-clearing of gas from the region before the SNe explode.

We present recent JWST NIRCam imaging observations of SPT0615-JD (also known as the Cosmic Gems Arc), lensed by the galaxy cluster SPT-CL J0615-5746. The 5-arcsec-long arc is the most highly magnified $z>10$ galaxy known, straddling the lensing critical curve and revealing five star clusters with radii $\sim 1$ pc or less. We measure the full arc to have F200W 24.5 AB mag, consisting of two mirror images, each 25.3 AB mag with a magnification $\mu \sim 60$ (delensed 29.7 AB mag, $M_{UV} = -17.8$). The galaxy has an extremely strong Lyman break F115W$-$F200W $>3.2$ mag ($2\sigma$ lower limit), is undetected in all bluer filters ($< 2\sigma$), and has a very blue continuum slope redward of the break ($\beta = -2.7 \pm 0.1$), resulting in a photometric redshift $z_{phot} = 10.2 \pm 0.2$ (95% confidence) with no significant likelihood below $z < 9.8$. Based on SED fitting to the total photometry, we estimate an intrinsic stellar mass of $M_{*} \sim 2.4 - 5.6 \times 10^{7} M_{\odot}$, young mass-weighted age of $\sim 21 - 79$ Myr, low dust content ($A_V < 0.15$), and a low metallicity of $\lesssim 1\%~Z_{\odot}$. We identify a fainter third counterimage candidate within 2.2 arcsec of the predicted position, lensed to AB mag 28.4 and magnified by $\mu \sim 2$, suggesting the fold arc may only show $\sim60$% of the galaxy. SPT0615-JD is a unique laboratory to study star clusters observed within a galaxy just 460 Myr after the Big Bang.

Almost 50% of galaxies in the local Universe are in clusters or groups coexisting with both hot and cold gas components. In the present study, we observationally probed the cold-gas content of X-ray-selected massive galaxy clusters with spectroscopic redshift measured from the SDSS/SPIDERS survey. This paper focuses on the most massive structures: galaxy clusters with a mean mass of M$_{500c}$ = 2.7$\times 10^{14}$ M$_{\odot}$. We used a large number of background quasar optical spectra from SDSS DR16 to probe the diffuse T$=$10$^4$K gas in their intracluster medium. We first analysed a sample of spectra with known MgII absorbers, and then blindly stacked about 16,000 archival spectra at the redshifts of the foreground galaxy clusters. We tentatively ($3.7 \sigma$ significance) detect MgII in the clusters with an equivalent width EW(MgII $\lambda$2796) of 0.056$\pm$0.015 Å, corresponding to a column density of log [N(MgII)/cm$^{-2}$]=12.12$\pm0.1$. We tested our methodology by generating 22,000 mock SDSS spectra with MgII absorbers from Illustris-TNG50 cosmological magnetohydrodynamical simulations, combining photo-ionisation modelling and ray tracing. We also performed bootstrapping stacking at different cluster redshifts and stacked quasar spectra with no intervening clusters in the line of sight to measure the significance of our detection. These results are in line with the findings of recent, similar observational studies but challenge predictions from Illustris-TNG simulations. Together, our findings indicate that large amounts of cold gas may be found in the most massive structures of the Universe.

We present a new inelastic dark matter search: neutron stars in dark matter-rich environments capture inelastic dark matter which, for interstate mass splittings between about $45 - 285 \ \rm MeV$, will annihilate away before becoming fully trapped inside the object. This means a sizable fraction of the dark matter particles can annihilate while being outside the neutron star, producing neutron star-focused gamma-rays and neutrinos. We analyze this effect for the first time and target the neutron star population in the Galactic Center, where the large dark matter and neutron star content makes this signal most significant. Depending on the assumed neutron star and dark matter distributions, we set constraints on the dark matter-nucleon inelastic cross-section using existing H.E.S.S. observations. We also forecast the sensitivity of upcoming gamma-ray and neutrino telescopes to this signal, which can reach inelastic cross-sections as low as $\sim 2 \times 10^{-47} \ \rm cm^2$.

We consider the scattering of low-mass halo dark-matter particles in the hot plasma of the Sun, focusing on dark matter that interact with ordinary matter through a dark-photon mediator. The resulting ``solar-reflected'' dark matter (SRDM) component contains high-velocity particles, which significantly extend the sensitivity of terrestrial direct-detection experiments to sub-MeV dark-matter masses. We use a detailed Monte-Carlo simulation to model the propagation and scattering of dark-matter particles in the Sun, including thermal effects, with special emphasis on ultralight dark-photon mediators. We study the properties of the SRDM flux, obtain exclusion limits from various direct-detection experiments, and provide projections for future experiments, focusing especially on those with silicon and xenon targets. We find that proposed future experiments with xenon and silicon targets can probe the entire ``freeze-in benchmark,'' in which dark matter is coupled to an ultralight dark photon, including dark-matter masses as low as $\mathcal{O}$(keV). Our simulations and SRDM fluxes are publicly available.

We perform in this work an analysis of the background dynamics for $\alpha$-attractor models in the context of loop quantum cosmology. Particular attention is given to the determination of the duration of the inflationary phase that is preceded by the quantum bounce in these models. From an analysis of the general predictions for these models, it is shown that we can be able to put constraints in the parameter $\alpha$ of the potentials and also on the quantum model itself, in special the Barbero-Immirzi parameter.

We propose a novel leptogenesis scenario in the gauged $U(1)_{L_{\mu}-L_{\tau}}$ model. Achieving successful leptogenesis in the $U(1)_{L_{\mu}-L_{\tau}}$ symmetric phase is challenging due to the absence of a CP phase, caused by restriction from the gauge symmetry. To overcome this issue, we introduce an additional global symmetry, $U(1)_{B-L}$, and a scalar field $\Phi$ responsible for breaking this symmetry. Through the kinetic misalignment mechanism, the majoron field associated with $U(1)_{B-L}$ symmetry breaking has a kinetic motion in the early universe. Subsequently, time-dependent majoron field background induces the background CP phase dynamically, leading to successful leptogenesis in the $U(1)_{L_{\mu}-L_{\tau}}$ symmetric phase. Furthermore, majoron itself serves as a dark matter candidate in this scenario. As one of the phenomenological applications, we consider the model that can also explain the muon $g-2$ anomaly.

This study undertakes a reconsideration of the potential for a first-order electroweak phase transition, focusing on the next-to-minimal two Higgs doublet model (N2HDM). Our exploration spans diverse parameter spaces associated with the phase transition, with a particular emphasis on examining the generation of stochastic Gravitational Waves (GW) resulting from this transition. The obtained results are meticulously compared against data from prominent gravitational wave observatories, and the possibility of their detection in the future GW observations have been established. In passing by we analyse the strength of the phase transition through the production of entropy during the electroweak phase transition.

We present novel findings concerning the parameter space of axion stars, extended object forming in dense dark matter environments through gravitational condensation. We emphasize their formation within the dense minihalos that potentially surround primordial black holes and in axion miniclusters. Our study investigates the relation between the radius and mass of an axion star in these dense surroundings, revealing distinct morphological characteristics compared to isolated scenarios. We explore the implications of these results when applied to gravitational microlensing from extended objects and gravitational wave detection from compact binaries, leading to insights on the observational constraints from such ``halo'' axion stars. We provide a constraint on the fraction of the galactic population of axion stars from their contribution to the microlensing events from the EROS-2 survey, using the numerical resolution of the Schrödinger-Poisson equation.

We study the dynamical evolution of superradiant instabilities of rotating black holes for multiple axion fields with comparable masses, motivated by string theory constructions, which typically exhibit a large number of light axions, with a broad range of masses. We show, in particular, that even though superradiant clouds for the heavier axion species grow faster, they are eventually reabsorbed by the black hole as the latter amplifies the lighter axion field(s), analogously to the dynamics of different species competing for the same resources in an ecosystem. We also incorporate in our study the effects of accretion and gravitational wave emission. We further demonstrate that the existence of multiple axion species with comparable masses may have a substantial impact on the stochastic gravitational wave background produced by axion clouds around black hole binary merger remnants, which could be probed with planned detectors.

Ultralight bosons can condense to form the so-called bosonic clouds around spinning black holes by superradiance instability. When quantum effects are taken into account, the classical black holes were replaced by exotic compact objects including area quantized black holes. In this work, we consider the superradiant instabilities of massive scalar fields around area quantized Kerr black hole. We introduce the reflectivity of area quantized black hole possesses the distinct discrete feature, and the scalar fields have the superradiant modes solution only within the specific mass range. In addition, the area quantization may terminate the superradiance when the black hole spins down, or even suppress the formation of the bosionic cloud.