Papers for Wednesday, Dec 04 2024

We use the star formation rate-stellar mass relation ($SFR - {M_ * }$) of galaxies to measure their luminosity distances as well as estimate cosmological parameters. We obtain a sample of 341 high-redshift galaxies at $5 < z < 14$, and 51 others at $1.8 < z < 3.5$ from the JWST observations, which can be used to investigate the correlation between the star formation rate and stellar mass, and determine their cosmological distance. The relation $SFR \propto {\kern 1pt} {({M_ * }/{M_ \odot })^\gamma }$ can be applied to check the nonlinear correlation between the star formation rate and stellar mass ${M_ * }$ of galaxies and give their distance. On the other hand, it also can be used to test if the $SFR - {M_ * }$ relation depends on the absolute UV magnitude ${M_{UV}}$, and the data suggest that the $SFR $ is positively correlated with stellar mass, the $SFR - {M_ * }$ relation of the high-redshift galaxies has an absolute magnitude dependence. Finally, we apply a combination of galaxies and Type Ia supernova (SNIa) Pantheon to test the property of dark energy concerning whether or not its density deviates from the constant, and give the statistical results.

Determining the internal structure of Uranus is a key objective for planetary science. Knowledge of Uranus's bulk composition and the distribution of elements is crucial to understanding its origin and evolutionary path. In addition, Uranus represents a poorly understood class of intermediate-mass planets (intermediate in size between the relatively well studied terrestrial and gas giant planets), which appear to be very common in the Galaxy. As a result, a better characterization of Uranus will also help us to better understand exoplanets in this mass and size regime. Recognizing the importance of Uranus, a Keck Institute for Space Studies (KISS) workshop was held in September 2023 to investigate how we can improve our knowledge of Uranus's internal structure in the context of a future Uranus mission that includes an orbiter and a probe. The scientific goals and objectives of the recently released Planetary Science and Astrobiology Decadal Survey were taken as our starting point. We reviewed our current knowledge of Uranus's interior and identified measurement and other mission requirements for a future Uranus spacecraft, providing more detail than was possible in the Decadal Survey's mission study and including new insights into the measurements to be made. We also identified important knowledge gaps to be closed with Earth-based efforts in the near term that will help guide the design of the mission and interpret the data returned.

We have conducted a sensitivity analysis on the mid-infrared spectral decomposition of galaxies and the modeling of the PAH emission spectrum with the NASA Ames PAH Infrared Spectroscopic Database (PAHdb) to assess the variance on the average galaxy PAH population properties under a grid of different modeling parameters. We find that the SL and SL+LL Spitzer-IRS decomposition with PAHFIT provides consistent modeling and recovery of the 5-15 $\mu$m PAH emission spectrum. For PAHdb modeling, application of a redshift to the calculated spectra to account for anharmonic effects introduces a $15\%$-$20\%$ variance on the derived parameters, while its absence improves the fits by $\sim13\%$. The 4.00-$\alpha$ release of PAHdb achieves the complete modeling of the 6-15 $\mu$m PAH spectrum, including the full 6.2 $\mu$m band, improving the average fitting uncertainty by a factor of 2. The optimal PAHdb modeling configuration requires selection of pure PAHs without applying a redshift to the bands. Although quantitatively the PAHdb-derived parameters change under different modeling configurations or database versions, their variation follows a linear scaling, with previously reported trends remaining qualitatively valid. PAHdb modeling of JWST observations, and JWST observations smoothed and resampled to the Spitzer-IRS resolution and dispersion have consistent PAHdb derived parameters. Decomposition with different codes, such as PAHFIT and CAFE, produce PAH emission spectra with noticeable variation in the 11-15~$\mu$m region, driving a $\sim7\%$ difference in the neutral PAH fraction under PAHdb modeling. A new library of galaxy PAH emission templates is delivered to be utilized in galaxy SED modeling.

The bulk abundances of exoplanetesimals can be measured when they are accreted by white dwarfs. Recently, lithium from the accretion of exoplanetesimals was detected in relatively high levels in multiple white dwarfs. There are presently three proposed hypotheses to explain the detection of excess lithium in white dwarf photospheres: Big Bang and Galactic nucleosynthesis, continental crust, and an exomoon formed from spalled ring material. We present new observations of three previously known lithium-polluted white dwarfs (WD J1824+1213, WD J2317+1830, and LHS 2534), and one with metal pollution without lithium (SDSS J1636+1619). We also present atmospheric model fits to these white dwarfs. We then evaluate the abundances of these white dwarfs and two additional lithium-polluted white dwarfs that were previously fit using the same atmospheric models (WD J1644$-$0449 and SDSS J1330+6435) in the context of the three extant hypotheses for explaining lithium excesses in polluted white dwarfs. We find Big Bang and Galactic nucleosynthesis to be the most plausible explanation of the abundances in WD J1644$-$0449, WD J1824+1213, and WD J2317+1830. SDSS J1330+6435 will require stricter abundances to determine its planetesimal's origins, and LHS 2534, as presently modeled, defies all three hypotheses. We find the accretion of an exomoon formed from spalled ring material to be highly unlikely to be the explanation of the lithium excess in any of these cases.

Venus's climatic history provides powerful constraint on the location of the inner-edge of the liquid-water habitable zone. However, two very different histories of water on Venus have been proposed: one where Venus had a temperate climate for billions of years, with surface liquid water, and the other where a hot early Venus was never able to condense surface liquid water. Here we offer a novel constraint on Venus's climate history by inferring the water content of its interior. By calculating the present rate of atmospheric destruction of H$_2$O, CO$_2$ and OCS, which must be restored by volcanism to maintain atmospheric stability, we show Venus's interior is dry. Venusian volcanic gases have at most a 6% water mole fraction, substantially drier than terrestrial magmas degassed at similar conditions. The dry interior is consistent with Venus ending its magma ocean epoch desiccated and thereafter having had a long-lived dry surface. Volcanic resupply to Venus's atmosphere therefore indicates that the planet has never been `liquid-water' habitable.

UBV CCD observations of standard stars selected from Landolt (2009, 2013) were performed using the 1-meter telescope (T100) of the TÜBİTAK National Observatory equipped with a back-illuminated and UV enhanced CCD camera and Bessell UBV filters. Observations span a long time from the years 2012 to 2024, 50 photometric nights in total. Photometric measurements were used to find the standard transformation relations of the T100 photometric system. The atmospheric extinction coefficients, zero points and transformation coefficients of each night were determined. It could not be found time dependence of the secondary extinction coefficients. However, it was determined that the primary extinction coefficients decreased until the year 2019 and increased after that year. It could not be found a strong seasonal variation of the extinction coefficients. Small differences in seasonal median values of them were used to attempt to find the atmospheric extinction sources. We found calculated minus catalogue values for each standard star, $\Delta(U-B)$, $\Delta(B-V)$ and $\Delta V$. Means and standard deviations of $\Delta(U-B)$, $\Delta(B-V)$ and $\Delta V$ were estimated to be 1.4$\pm$76, 1.9$\pm$18 and 0.0$\pm$36 mmag, respectively. We found that our data well matched Landolt's standards for $V$ and $B-V$, i.e. there are no systematic differences. However, there are systematic differences for $U-B$ between the two photometric systems, which is probably originated from the quantum efficiency differences of the detectors used in the photometric systems, although the median differences are relatively small ($|\Delta(U-B)|$< 50 mmag) for stars with $-0.5<U-B~{(\rm mag)} <1.6$ and $0.2<B-V~{(\rm mag)} <1.8$. As an overall result, we conclude that the transformation relations found in this study can be used for standardized photometry with the T100 photometric system.

JWST has revealed a population of compact `Little Red Dots' (LRDs) at $z\gtrsim4$, with red rest-frame optical and blue UV colors. These objects are likely compact dusty starbursts or heavily reddened AGNs, playing a pivotal role in early black hole growth, dust production, and stellar assembly. We introduce a new photometric selection to identify LRDs over a broad range in redshifts and rest-frame UV-to-NIR colors enabling a more complete census of the population. This method identifies 248 LRDs with F444W$<27$ mag over 263 arcmin$^2$ in the JADES, PRIMER-COSMOS, and UDS fields with MIRI coverage, increasing the number density by $\times$1.7 compared to previous samples, suggesting that previous census were underestimated. Most LRDs are detected in MIRI/F770W but only 7% (17) are detected in F1800W. We use MIRI-based rest-frame [1$-$3 $\mu$m] colors to trace dust emission. F1800W-detected LRDs have a median [1$-$3 $\mu$m]$=1.5$ mag, with a broad scatter indicative of diverse dust emission properties. About 20% exhibit [1$-$3 $\mu$m]$<1$ mag colors consistent with negligible dust emission, but the majority show significant dust emission at 3 $\mu$m (f$^{\rm dust}_{3\mu m}\lesssim0.8$) from the galaxy ISM or a hot-dust-deficient AGN torus. A correlation between bluer UV-to-NIR colors and stronger IR emission suggests that the bluest LRDs may resemble unobscured QSOs. We report a LRD at $z_{\rm spec}=3.1386$, detected in MIRI, Spitzer/MIPS, and Herschel/PACS. Its IR SED rises steeply at $\lambda_{\rm rest}>6~\mu$m and peaks near $\sim40~\mu$m, providing the first direct evidence of warm dust emission (T$=50-100$ K) in a LRD.

Cosmological analyses of redshift space clustering data are primarily based on using luminous ``red'' galaxies (LRGs) and ``blue'' emission line galaxies (ELGs) to trace underlying dark matter. Using the large high-fidelity high-resolution MillenniumTNG (MTNG) and Astrid simulations, we study these galaxies with the effective field theory (EFT)-based field level forward model. We confirm that both red and blue galaxies can be accurately modeled with EFT at the field level and their parameters match those of the phenomenological halo-based models. Specifically, we consider the state of the art Halo Occupation Distribution (HOD) and High Mass Quenched (HMQ) models for the red and blue galaxies, respectively. Our results explicitly confirm the validity of the halo-based models on large scales beyond the two-point statistics. In addition, we validate the field-level HOD/HMQ-based priors for EFT full-shape analysis. We find that the local bias parameters of the ELGs are in tension with the predictions of the LRG-like HOD models and present a simple analytic argument explaining this phenomenology. We also confirm that ELGs exhibit weaker non-linear redshift-space distortions (``fingers-of-God''), suggesting that a significant fraction of their data should be perturbative. We find that the response of EFT parameters to galaxy selection is sensitive to assumptions about baryonic feedback, suggesting that a detailed understanding of feedback processes is necessary for robust predictions of EFT parameters. Finally, using neural density estimation based on paired HOD-EFT parameter samples, we obtain optimal HOD models that reproduce the clustering of Astrid and MTNG galaxies.

If primordial black holes (PBH) are present in the early universe, their contribution to the energy budget grows relative to that of radiation and generically becomes dominant unless the initial abundance is exponentially small. This black hole domination scenario is largely unconstrained for PBHs with masses $\lesssim 10^9\,\mathrm{g}$, which evaporate prior to Big Bang nucleosynthesis. However, if the era of PBH domination is sufficiently long, the PBHs form clusters and can merge appreciably within these objects. We calculate the population statistics of these clusters within the Press-Schechter formalism and find that, for a wide range of PBH masses and Hubble rates at the onset of PBH domination, the mergers within PBH clusters can exhibit runaway behavior, where the majority of the cluster will eventually form a single black hole with a mass much greater than the original PBH mass. These mergers can dramatically alter the PBH mass distribution and leave behind merged relic black holes that evaporate after Big Bang nucleosynthesis and yield various observational signatures, excluding parameter choices previously thought to be viable

JWST has allowed the discovery of a significant population of galaxies at z > 10. Our understanding of the astrophysical properties of these ultra-high redshift galaxies relies on fitting templates, developed using astrophysical models representing our current understanding of high-redshift galaxies. In this work, the highest confidence recent JWST spectroscopic observations are used to evaluate the performance of several high-redshift templates based on two tests: (1) comparing photometric redshifts against spectroscopic redshifts; and (2) comparing the reconstructed spectral energy distributions against observed SEDs. Strict sample selection and error-propagation by bootstrapping is employed to make results robust towards future JWST systematics mitigation. It is shown that some templates perform adequately at high redshift prediction, given a sample selection by observational filters and depth. Other templates work better for SED fitting, but a few objects remain unrepresented in their spectra. We conclude that although templates are usable, models are not yet able to reliably extract astrophysical properties.

MIR spectra imply considerable chemical diversity in the inner regions of protoplanetary discs: some are H2O-dominated, others by CO2. Sublimating ices from radially drifting dust grains are often invoked to explain some of this diversity, particularly the H2O-rich discs. We use a 1D protoplanetary disc evolution code to model how radially drifting dust grains that transport ices inwards to snowlines impact the chemistry of the inner regions of protoplanetary discs. We explore differences between smooth discs and those where radial drift is impeded by dust trapping outside gas gaps and quantify the effects of gap location and formation time. Discs evolve through an initial H2O-rich phase due to sublimating ices, followed by a CO2-rich phase as H2O vapour advects onto the star and CO2 advects into the inner disc from its snowline. The inclusion of traps hastens the transition between the phases, raising the CO2/H2O ratio; gaps opened early or close-in produce lower increases by blocking more CO2 ice from reaching the inner disc. This leads to a potential correlation between CO2/H2O and gap location that occurs on Myr timescales for fiducial parameters. We produce synthetic spectra from the models which we analyse with 0D LTE slab models to understand how this evolution may be expressed observationally. Whether the evolution can be retrieved depends on the contribution of dust grains to the optical depth: dust that couples to the gas after crossing the H2O snowline can add to the continuum optical depth and obscure the delivered H2O, largely hiding the evolution in its visible column density. However, the CO2/H2O visible column density ratio is only weakly sensitive to dust continuum obscuration. This suggests it may be a clearer tracer of the impact of transport on chemistry than individual column densities for spectra that show weak features probing deep enough in the disc. (Abridged)

This study compiles stellar populations and internal properties of ultra-diffuse galaxies (UDGs) to highlight correlations with their local environment, globular cluster (GC) richness, and star formation histories. Complementing our sample of 88 UDGs, we include 36 low-surface brightness dwarf galaxies with UDG-like properties, referred to as NUDGes (nearly-UDGs). All galaxies were studied using the same spectral energy distribution fitting methodology to explore what sets UDGs apart from other galaxies. We show that NUDGes are similar to UDGs in all properties except for being, by definition, smaller and having higher surface brightness. We find that UDGs and NUDGes show similar behaviours in their GC populations, with the most metal-poor galaxies hosting consistently more GCs on average. This suggests that GC content may provide an effective way to distinguish extreme galaxies within the low surface brightness regime alongside traditional parameters like size and surface brightness. We confirm previous results using clustering algorithms that UDGs split into two main classes, which might be associated with the formation pathways of a puffy dwarf and a failed galaxy. The clustering applied to the UDGs+NUDGes dataset yields an equivalent result. The difference in mass contained in the GC system suggests that galaxies in different environments have not simply evolved from one another but may have formed through distinct processes.

Horizon-scale imaging with the Event Horizon Telescope (EHT) has provided transformative insights into supermassive black holes but its resolution and scope are limited by ground-based constraints such as the size of the Earth, its relatively slow rotation, and atmospheric delays. Space-based very long baseline interferometry (VLBI) offers the capability for studying a larger and more diverse sample of black holes. We identify a number of nearby supermassive black holes as prime candidates for horizon-scale imaging at millimeter wavelengths, and use source characteristics such as angular size, sky distribution, and variability timescales to shape the design of a space-based array. We identify specific metrics that serve as key predictors of image fidelity and scientific potential, providing a quantitative basis for optimizing mission design parameters. Our analysis demonstrates that the optimal configuration requires two space-based elements in high Earth orbits (HEO) that are not coplanar and are apparently counter-rotating. Our results delineate the key requirements for a space-based VLBI mission, enabling detailed studies of black hole shadows, plasma dynamics, and jet formation, advancing black hole astrophysics beyond the current capabilities of the EHT.

Recent JWST observations have allowed for the first time to obtain comprehensive measurements of the ionizing photon production efficiency $\xi_\text{ion} $ for a wide range of reionization-epoch galaxies. We explore implications for the inferred UV luminosity functions and escape fractions of ionizing sources in our suite of simulations. These are run with the GPU-based radiative transfer code ATON-HE and are calibrated to the XQR-30 Lyman-alpha forest data at $5<z<6.2$. For our fiducial source model, the inferred ionizing escape fractions increase from (6.1, 5.4, 4.9)% at $z=6$ to (14.4, 23.8, 29.4)% at $z=10$ for our (Fiducial, Early, Extremely Early) models in good agreement with extrapolations of lower redshift escape fraction measurements. Extrapolating observed luminosity functions beyond the resolution limit of the simulations to faint sources with $M_\text{UV}=-11$ increases the inferred escape fractions by a factor $\sim 1.5$ at $z=10$. For our oligarchic source model, where no ionizing photons are emitted in faint sources, the inferred escape fractions increase from 10% at $z=6$ to uncomfortably large values $>50$% at $z> 10$, disfavouring the oligarchic source model at very high redshift. The inferred effective clumping factors in our simulations are in the range of $3-6$, suggesting consistency between the observed ionizing properties of reionization-epoch galaxies and the ionizing photon budget in our simulations.

Increasing the angular resolution of an interferometric array requires placing its elements at large separations. This often leads to sparse coverage and introduces challenges to reconstructing images from interferometric data. We introduce a new interferometric imaging algorithm, KRISP, that is based on kernel methods, is statistically robust, and is agnostic to the underlying image. The algorithm reconstructs the complete Fourier map up to the maximum observed baseline length based entirely on the data without tuning by a user or training on prior images and reproduces images with high fidelity. KRISP works efficiently for many sparse array configurations even in the presence of significant image structure as long as the typical baseline separation is comparable to or less than the correlation length of the Fourier map, which is inversely proportional to the size of the target image.

Significant mass loss in the red supergiant (RSG) phase has great influence on the evolution of massive stars and their final fate as supernovae. We present near-infrared interferometric imaging of the circumstellar environment of the dust-enshrouded RSG WOH G64 in the Large Magellanic Cloud. WOH G64 was observed with the GRAVITY instrument at ESO's Very Large Telescope Interferometer (VLTI) at 2.0--2.45 micron. We succeeded in imaging the innermost circumstellar environment of WOH G64 -- the first interferometric imaging of an RSG outside the Milky Way. The reconstructed image reveals elongated compact emission with a semimajor and semiminor axis of ~2 and ~1.5 mas (~13 and 9 stellar radii), respectively. The GRAVITY data show that the stellar flux contribution at 2.2 micron at the time of our observations in 2020 is much lower than predicted by the optically and geometrically thick dust torus model based on the VLTI/MIDI data taken in 2005 and 2007. We found a significant change in the near-infrared spectrum of WOH G64: while the (spectro)photometric data taken at 1--2.5 micron before 2003 show the spectrum of the central RSG with H2O absorption, the spectra and JHK' photometric data taken after 2016 are characterized by a monotonically rising continuum with very weak signatures of H2O. This spectral change likely took place between December 2009 and 2016. On the other hand, the mid-infrared spectrum obtained in 2022 with VLT/VISIR agrees well with the spectra obtained before 2007. The compact emission imaged with GRAVITY and the near-infrared spectral change suggest the formation of hot new dust close to the star, which gives rise to the monotonically rising near-infrared continuum and the high obscuration of the central star. The elongation of the emission may be due to the presence of a bipolar outflow or effects of an unseen companion.

Observations of tidal disruption events (TDEs) on a timescale of years after the main flare show evidence of continued activity in the form of optical/UV emission, quasi-periodic eruptions, and delayed radio flares. Motivated by this, we explore the time evolution of these disks using semi-analytic models to follow the changing disk properties and feeding rate to the central black hole (BH). We find that thermal instabilities typically begin $\sim150-250\,{\rm days}$ after the TDE, causing the disk to cycle between high and low accretion states for up to $\sim10-20\,{\rm yrs}$. The high state is super-Eddington, which may be associated with outflows that eject $\sim10^{-3}-10^{-1}\,M_\odot$ with a range of velocities of $\sim0.03-0.3c$ over a span of a couple of days and produce radio flares. In the low state, the accretion rate slowly grows over many months to years as continued fallback accretion builds the mass of the disk. In this phase, the disk may reach luminosities of $\sim10^{41}-10^{42}\,{\rm erg\,s^{-1}}$ in the UV as seen in some late-time observations. We highlight the importance of the iron-opacity "bump" at $\approx2\times10^5\,{\rm K}$ in generating sufficiently high luminosities. This work suggests that joint optical/UV observations with radio monitoring could be key for following the disk state as the radio flares are produced.

Binary black hole (BBH) evolution in the discs of active galactic nuclei (AGN) is a promising channel for gravitational wave (GW)-driven mergers. It is however unclear whether binaries interacting with the surrounding disc undergo orbital contraction or expansion. We develop a simple analytic model of accreting BBHs in AGN discs to follow the orbital evolution from the disc-dominated regime at large separations into the GW-driven regime at small separations (the coupled `disc+GW'-driven evolution). We obtain that accreting binaries expand in thick discs with aspect ratio greater than a critical value ($> h_\mathrm{crit}$); whereas accreting binaries contract and eventually merge in thin discs ($< h_\mathrm{crit}$). Interestingly, accreting BBHs can experience faster mergers compared to non-accreting counterparts, with a non-monotonic dependence on the disc aspect ratio. The orbital contraction is usually coupled with eccentricity growth in the disc-dominated regime, which lead to accelerated inspirals in the GW-driven regime. We quantify the resulting BBH merger timescales in AGN discs ($\tau_\mathrm{merger} \sim 10^5 - 10^7$ yr) and estimate the associated GW merger rates ($\mathcal{R} \sim (0.2 - 5) \, \text{Gpc}^{-3} \text{yr}^{-1}$). Overall, accreting binaries may efficiently contract and merge in thin discs, hence this particular BBH-in-AGN channel may provide a non-negligible contribution to the observed GW merger event rate.

Data derived from general relativistic MHD simulations of accretion onto black holes can be used as input to a post-processing scheme that predicts the radiated spectrum. Combining a relativistic Compton-scattering radiation transfer solution in the corona with detailed local atmosphere solutions incorporating local ionization and thermal balance within the disk photosphere, it is possible to study both spectral formation and intrinsic spectral variability in the radiation from relativistic accretion disks. With this method, we find that radiatively-efficient systems with black holes of $10M_\odot$ accreting at $\approx 0.01$ in Eddington units produce spectra very similar to those observed in the hard states of X-ray binaries. The spectral shape above 10keV is well-described by a power-law with an exponential cutoff. Intrinsic turbulent variations lead to order-unity fluctuations in bolometric luminosity, logarithmic spectral slope fluctuations over a range $\approx 0.15$, and factor of 2 changes in the cut-off energy on timescales $\sim 50\, (M_{\rm BH}/10 M_\odot)$ms. Within the corona, the range of gas temperature spans more than an order of magnitude. The wide distribution of temperatures is central to defining the spectrum: the logarithmic spectral slope is harder by $\sim 0.3$ and the cut-off energy larger by a factor $\sim 10 - 30$ than if the coronal temperature everywhere were its mass-weighted mean.

As revealed by Hubble in 1928, our Universe is expanding. This discovery was fundamental to widening our horizons and our conception of space, and since then determining the rate at which our Universe is expanding has become one of the crucial measurements in cosmology. At the beginning of this century, these measurements revealed the unexpected behavior that this expansion is accelerating and allowed us to have a first glimpse of the dark components that constitute $\sim$95\% of our Universe. Cosmic chronometers represent a novel technique to obtain a cosmology-independent determination of the expansion of the Universe, based on the differential age dating of a population of very massive and passively evolving galaxies. Currently, with this new cosmological probe it is possible to constrain the Hubble parameter with an accuracy of around 5\% at $z\sim0.5$ up to 10-20\% at $z\sim2$. In this Chapter, the cosmic chronometers approach is presented, describing the method and how an optimal sample can be selected; it is then discussed how the most recent measurements of the expansion history of the Universe have been obtained with this approach, as well as the cosmological constraints that can be derived. Particular attention will be given to the systematics involved in this approach and the treatment to properly take them into account. We conclude by presenting forecasts that show how future spectroscopic surveys will significantly boost the accuracy of this method and open the possibility to a percent determination of the Hubble constant, making cosmic chronometers a powerful independent tool to derive information on the expansion history of the Universe.

The James Webb Space Telescope (JWST) has revealed an unexpectedly high abundance of UV luminous galaxies at redshifts $z\gtrsim 10$, challenging `standard' galaxy formation models. This study investigates the role of rapidly rotating (massive) stars undergoing chemically homogeneous evolution (CHE) in reconciling this potential tension. These stars are more compact, hotter, and exhibit enhanced UV emission. We find that the rest-frame UV luminosity of star-forming galaxies can be significantly enhanced by a factor of $\sim 3-6$ when CHE stars above a minimum initial mass of $m_{\star,\min}^{\rm CHE}\sim 2-10\ \rm M_\odot$ account for more than half of the total stellar mass following a Salpeter initial mass function. As a result, the UV luminosity functions observed at $z\sim 12-16$ can be reproduced with less extreme values of star formation efficiency and UV luminosity stochastic variability. Our results highlight the potential of CHE in explaining the UV-bright galaxy populations detected by JWST and call for future work to explore the broader astrophysical implications of CHE and its associated phenomena in the early universe, such as gamma-ray bursts, compact object binaries, and metal enrichment.

ALMA-IMF imaged 15 massive protoclusters down to a resolution of of 2 kau scales, identifying about 1000 star-forming cores. The mass and luminosity of these cores, which are fundamental physical characteristics, are difficult to determine, a problem greatly exacerbated at the distances >2 kpc of ALMA-IMF protoclusters. We combined new datasets and radiative transfer modeling to characterize these cores. We estimated their mass-averaged temperature and the masses these estimates imply. For 1/6 of the sample, we measured the bolometric luminosities, implementing deblending corrections when necessary. We used spectral energy distribution (SED) analysis obtained with the PPMAP Bayesian procedure, which aims to preserve the best angular resolution of the input data. We extrapolated the luminosity and dust temperature images provided by PPMAP at 2.5" resolution to estimate those of individual cores, which were identified at higher angular resolution. To do this, we applied approximate radiative transfer relationships between the luminosity of a protostar and the temperature of its surrounding envelope and between the external heating of prestellar cores and their temperatures. For the first time, we provide data-informed estimates of dust temperatures for 883 cores identified with ALMA-IMF: 17-31 K and 28-79 K (5th and 95th percentiles, up to 127 K) for the 617 prestellar and 266 protostellar cores, respectively. We also measured protostellar luminosities spanning 20-80 000 Lsun. For hot cores, we estimated systematically lower temperatures than studies based on complex organic molecules. We established a mass-luminosity evolutionary diagram, for the first time at the core spatial resolution and for a large sample of high-mass protostellar cores. The ALMA-IMF data favor a scenario in which protostars accrete their mass from a larger mass reservoir than their host cores.

Solar filaments are spectacular objects in the solar atmosphere, consisting of accumulations of cool, dense, and partially ionized plasma suspended in the hot solar corona against gravity. The magnetic structures that support the filament material remain elusive, partly due to the lack of high resolution magnetic field measurements in the chromosphere and corona. In this study, we reconstruct the magnetic structures of a solar intermediate filament using EUV observations and two different methods, to follow the injection of hot material from a B-class solar flare. Our analysis reveals the fine-scale magnetic structures of the filament, including a compact set of mutually wrapped magnetic fields encasing the cool filament material, two groups of helical magnetic structures intertwining with the main filament, and a series of arched magnetic loops positioned along the filament. Additionally, we also find that the northern footpoints of the helical structures are rooted in the same location, while their southern footpoints are rooted in different areas. The results obtained in this study offer new insights into the formation and eruption mechanisms of solar filaments.

We report the spectroscopic analysis of ten satellite galaxy candidates in the sphere of influence of the Sombrero galaxy (M104, NGC4594), based on data obtained with IFUM (Integral Field Units for Magellan). Based on their newly-observed recessional velocities, we confirm that nine of these candidates are satellite galaxies of M104, with one being a background dwarf galaxy. All ten dwarfs have stellar masses $2\times10^{7}\,M_{\odot}$ to $1\times10^{9}\,M_{\odot}$ and mean weighted metallicities $-1.7<\langle{[\mathrm{M/H}]}\rangle<-0.3$. Although these dwarfs are predominantly old, with stellar populations $\sim5-11\,$Gyr. However, this sample contains a local example of a low-mass "Green Pea" candidate, it exhibits extreme optical emission features and broad emission line features ($\sigma\sim250\,\mathrm{km\,s^{-1}}$) reminiscent of high-redshift Ly$\alpha$/LyC photon leaking galaxies. Using the newly-acquired recessional velocities of the nine satellites of M104, we find no evidence of coherent satellite motions unlike other nearby $L_*$ galaxy environments. Given the small sample, this results does not statistically rule out such coherent motions. There remain 60 satellite candidates of M104 for which future spectroscopy can more reliably test for such motion. Using the observed dwarf galaxies as tracers of the gravitational potential of M104, we estimate the dynamical mass of M104, $M_{dyn}=(12.4\pm6.5)\times10^{12}M_{\odot}$, and find that, making a reasonable estimate of M104's gas mass, $>90\%$ of its baryons are missing. These results agree with previous measurements of M104's dynamical mass.

Stellar-mass black holes descend from high-mass stars, most of which had stellar binary companions. However, the number of those binary systems that survive the binary evolution and black hole formation is uncertain by multiple orders of magnitude. The survival rate is particularly uncertain for massive stars with low-mass companions, which are thought to be the progenitors of most black hole X-ray binaries. We present a search for close black hole companions (separations less than 20 solar radii) to AFGK-type stars in TESS, i.e. the non-accreting counterparts to and progenitors of low-mass X-ray binaries. Such black holes can be detected by the tidally induced ellipsoidal deformation of the visible star, and the ensuing photometric light-curve variations. From an initial sample of 4.7 million TESS stars, we have selected 457 candidates for such variations. However, spectroscopic followup of 251 of them shows that none are consistent with a close black hole companion. On the basis of this non-detection, we determine (2 $\sigma$ confidence) that fewer than one in $10^5$ Solar-type stars in the Solar neighbourhood host a short-period black hole companion. This upper limit is in tension with a number of ``optimistic'' population models in the literature that predict short-period black hole companions around one in $10^{4-5}$ stars. Our limits are still consistent with other models that predict only a few in $10^{7-8}$.

Studying the early stages of transient events provides crucial information about the fundamental physical processes in cataclysmic variables (CVs). However, determining an appropriate observation mode immediately after the discovery of a new transient presents challenges due to significant uncertainties regarding its nature. We developed a framework designed for autonomous decision making in prompt follow-up observations of CVs using the Kanata 1.5-m telescope. The system, named Smart Kanata, first estimates the class probabilities of variable star types using a generative model. It then selects the optimal observation mode from three possible options based on the mutual information calculated from the class probabilities. We have operated the system for ~300 days and obtained 21 samples, among which automated observations were successfully performed for a nova and a microlensing event. In the time-series spectra of the nova V4370 Oph, we detected a rapid deepening of the absorption component of the H_alpha line. These initial results demonstrate the capability of Smart Kanata in facilitating rapid observations and improving our understanding of outbursts and eruptions of CVs and other galactic transients.

Atmospheric studies are essential for elucidating the formation history, evolutionary processes, and atmospheric dynamics of exoplanets. High-resolution transmission spectroscopy offers the advantage of detecting subtle variations in stellar spectral profiles, thereby enabling the identification of the sources of observed signals. In this study, we present the transmission spectra of the exoplanet WASP-77Ab, a hot Jupiter with a 1.36-day orbital period around a G8 host star with $V=11.29$ mag. These observations were conducted using the high-resolution spectrograph ESPRESSO at the Very Large Telescope over three transit events. We analyze the Rossiter-McLaughlin effect for WASP-77A and determine a projected spin-orbit angle of ${\lambda = 16.131^{\circ}}^{+2.106}_{-2.324}$, indicating that the planet's orbit is nearly aligned. Following the generation of transmission spectra for the three nights, we model and correct for center-to-limb variation and the Rossiter-McLaughlin effects. In the residual transmission spectra, we detect H$\alpha$, H$\beta$ and CaII H with a significance exceeding 3.5$\sigma$. After applying 0.1-0.5 Å masks to the cores of these lines to mitigate stellar contamination, all them still shows visible absorptions although not significant, suggesting at least partial planet contribution to them. Therefore, we are yet unable to confirm or reject the planetary origin of these spectral signals based on the current data set. Further investigation of WASP-77Ab's atmosphere, particularly in areas beyond the terminator region, is essential to illuminate the planet's two-dimensional atmospheric structure.

Super-Alfvenic turbulence is important for many astrophysical objects, particularly galaxy clusters. In this paper, we explore the accuracy of Synchrotron Intensity Gradients (SIGs) and X-ray intensity gradients to map magnetic fields in super-Alfvenic turbulence for a set of astrophysically relevant parameters of turbulent driving. Analyzing our synthetic observations, we report a good accuracy for both techniques. Our results are suggestive that other types of Gradient Technique (GT) can be successfully employed to trace magnetic fields within super-Alfvenic sub-sonic turbulence.

The Legacy Survey of Space and Time, to be conducted with the Vera C. Rubin Observatory, is poised to revolutionize our understanding of the Solar System by providing an unprecedented wealth of data on various objects, including the elusive interstellar objects (ISOs). Detecting and classifying ISOs is crucial for studying the composition and diversity of materials from other planetary systems. However, the rarity and brief observation windows of ISOs, coupled with the vast quantities of data to be generated by LSST, create significant challenges for their identification and classification. This study aims to address these challenges by exploring the application of machine learning algorithms to the automated classification of ISO tracklets in simulated LSST data. We employed various machine learning algorithms, including random forests (RFs), stochastic gradient descent (SGD), gradient boosting machines (GBMs), and neural networks (NNs), to classify ISO tracklets in simulated LSST data. We demonstrate that GBM and RF algorithms outperform SGD and NN algorithms in accurately distinguishing ISOs from other Solar System objects. RF analysis shows that many derived Digest2 values are more important than direct observables in classifying ISOs from the LSST tracklets. The GBM model achieves the highest precision, recall, and F1 score, with values of 0.9987, 0.9986, and 0.9987, respectively. These findings lay the foundation for the development of an efficient and robust automated system for ISO discovery using LSST data, paving the way for a deeper understanding of the materials and processes that shape planetary systems beyond our own. The integration of our proposed machine learning approach into the LSST data processing pipeline will optimize the survey's potential for identifying these rare and valuable objects, enabling timely follow-up observations and further characterization.

Accurate $N$-body simulations of multiple systems such as binaries and triples are essential for understanding the formation and evolution of interacting binaries and binary mergers, including gravitational wave sources, blue stragglers and X-ray binaries. The logarithmic time-transformed explicit symplectic integrator (LogH), also known as algorithmic regularization, is a state-of-the-art method for this this http URL, we show that this method is accurate for isolated Kepler orbits because of its ability to trace Keplerian trajectories, but much less accurate for hierarichal triple systems. The method can lead to an unphysical secular evolution of inner eccentricity in Kozal-Lidov triples, despite a small energy error. We demonstrate that hybrid methods, which apply LogH to the inner binary and alternative methods to the outer bodies, are significantly more effective, though not symplectic. Additionally, we introduce a more efficient hybrid method, BlogH, which eliminates the need for time synchronization and is time symmetric. The method is implemented in the few-body code SDAR. We explore suitable criteria for switching between the LogH and BlogH methods for general triple systems. These hybrid methods have the potential to enhance the integration performance of hierarchial triples.

We present the statistic results of GeV spectral breaks of bright gamma-ray flat-spectrum radio quasars (FSRQs) in the energy range of 0.1-10 GeV based on New Pass 8 date of the Large Area Telescope abroad on Fermi Gamma-ray Space Telescope. We have fitted the 15-year average gamma-ray spectra of 755 FSRQs by using both a broken power-law (BPL) and the Logarithmic Parabolas (LP) models, and obtained 87 bright gamma-ray FSRQs with their integrated photon fluxes greater than $2.16\times 10^{-8}$ cm$^{-2}$ s$^{-1}$. From our results, the FSRQ population shows similar preferences for both the BPL and LP models in gamma-ray spectral fitting and clustering analysis suggests that BPL-preferred and LP-preferred FSRQs belong to the same category. Our results indicate that GeV spectral breaks in bright gamma-ray FSRQs are located at $\rm 2.90\pm1.92$ GeV in the rest frame, and the observed change in photon index is $\rm \Delta \gamma =0.45 \pm 0.19$, which is consistent with the expected value for a cooling break of electrons scattering seed photons.

Current reverberation mapping (RM) studies primarily focus on single emission lines, particularly the \hb\ line, which may not fully reveal the geometry and kinematic properties of the broad-line region (BLR). To overcome this limitation, we conducted multiline RM observations on two highly variable active galactic nuclei (AGNs), KUG 1141+371 and UGC 3374, using the Lijiang 2.4 m telescope. Our goal was to investigate the detailed structure of different regions within the BLR. We measured the time lags of multiple broad emission lines (\ha, \hb, \hg, \hei, and \heii) and found clear evidence of radial ionization stratification in the BLRs of both AGNs. Velocity-resolved RM analysis revealed distinct geometry and kinematics between the inner and outer regions of the BLRs. Assuming that velocity-resolved lags reflect the kinematics of BLR, our observations indicate that: (1) in KUG 1141+371, the inner BLR exhibits outflow signatures, while the outer region is consistent with virialized motion; (2) in UGC 3374, the inner region displays virial motion, while the outer region shows inflow. Furthermore, we detected ``breathing" behavior in the outer BLR regions of both AGN, while the inner BLR regions show ``anti-breathing", which may be linked to intrinsic BLR properties. We discuss these findings in the context of various BLR formation models, highlighting importance of long-term, multiline RM campaigns in understanding of BLR structure and evolution. Additionally, our results suggest that the observed stratification in BLR geometry and kinematics may contribute to the scatter in black hole mass estimates and the rapid changes in velocity-resolved RM signatures reported in recent studies.

The rate coefficients of various isotopic variations of the H2+ + H2 and H3+ + H2 reactions in the 10-250 K temperature range were measured using a cryogenic 22 pole radio frequency ion trap. The processes involving diatomic ions were found to behave close to the Langevin rate, whereas temperature-dependent rate coefficients were obtained for the four isotopic exchange processes with triatomic ions. Fitting the experimental data using a chemical code allowed us in specific cases to constrain rate coefficients that were not directly measured in the ion trap. The reported rate coefficients suggest a more efficient hydrogenation of deuterated H3+ forms than usually assumed in astrochemical models, which might affect deuteration rates in warmer environments.

We conduct the first-ever Pulsar Polarization Array (PPA) analysis to detect the ultralight Axion-Like Dark Matter (ALDM) using the polarization data of 22 millisecond pulsars from the third data release of Parkes Pulsar Timing Array. As one of the major dark matter candidates, the ultralight ALDM exhibits a pronounced wave nature on astronomical scales and offers a promising solution to small-scale structure issues within local galaxies. While the linearly polarized pulsar light travels through the ALDM galactic halo, its position angle (PA) can be subject to an oscillation induced by the ALDM Chern-Simons coupling with electromagnetic field. The PPA is thus especially suited for detecting the ultralight ALDM by correlating polarization data across the arrayed pulsars. To accomplish this task, we develop an advanced Bayesian analysis framework that allows us to construct pulsar PA residual time series, model noise contributions properly and search for pulsar cross-correlations. We find that for an ALDM density of $\rho_0=0.4\,\textrm{GeV}/\textrm{cm}^3$, the Parkes PPA offers the best global limits on the ALDM Chern-Simons coupling, namely $\lesssim 10^{-13.5}-10^{-12.2}~{\rm GeV}^{-1}$, for the mass range of $10^{-22} - 10^{-21}~{\rm eV}$. The crucial role of pulsar cross-correlation in recognizing the nature of the derived limits is also highlighted.

Recent high-resolution observations at millimeter (mm) and sub-mm reveal a diverse spatial distribution for sub-pc scale dense cores within star-forming regions, ranging from clustered to aligned arrangements. To address the increasing volume of observational and simulation data, we introduce "alignment parameters" as a quantitative and reproducible method to automatically assess core alignment. We first demonstrate the effectiveness of these parameters by applying them to artificial test clumps and comparing the results with labels from visual inspection. A threshold value is then proposed to differentiate between "clustered" and "aligned" categories. Subsequently, we apply these parameters to dense cores identified from a sample of ALMA 1.3 mm dust continuum images in high-mass star-forming regions. Analysis exploring correlations between alignment parameters and clump properties rules out the presence of moderate or strong correlation, indicating that clump properties do not appear to strongly influence the outcome of fragmentation. One possible explanation for this is that the fragmentation process is chaotic, meaning that small variations in initial conditions can lead to significant differences in fragmentation outcomes, thus obscuring any direct link between clump properties and core alignment/distribution.

Considering radio-interferometric observations, we present a fast and efficient estimator to compute the binned angular bispectrum (ABS) from gridded visibility data. The estimator makes use of Fast Fourier Transform (FFT) techniques to compute the bispectrum covering all possible triangle shapes and sizes. Here, we present the formalism of the estimator and validate it using simulated visibility data for the Murchison Widefield Array (MWA) observations at $\nu=154.25$ MHz. We find that our estimator is able to faithfully recover the ABS of the simulated sky signal with $\approx 10-15 \%$ accuracy for a wide variety of triangle shapes and sizes across the range of angular multipoles $46 \le \ell \le 1320$. In future work, we plan to apply this to actual data and also generalize it to estimate the three-dimensional redshifted 21-cm bispectrum.

Bars are important in the secular evolution of galaxies. This study is aimed at exploring the reasons why some galaxies have bars at redshift $z=0$ while others do not. We use ellipse fitting to measure the properties and evolution of bars in the IllustrisTNG cosmological simulation. By using the K-S two-sample test and tracing their evolutionary changes, we analyze the parameter differences between barred and unbarred galaxies. The properties of galaxies with short bars are also studied. When tracing all disk galaxies at $z=0$ back to $z=1$, all of them show similar bar features at $z=1$. The fraction of bars increases in barred and short-bar galaxies but decreases in unbarred galaxies during $z=1-0$. In the case of disk galaxies with stellar mass log$(M_*/M_\odot)> 10.8$, nurture (mainly mergers) plays the most important role in suppressing or destroying bars. Bars are more likely to endure in galaxies that experience fewer mergers, which can be quantified by smaller stellar halos and ex-situ mass fractions. Approximately 60\% of the unbarred galaxies in the local Universe once had a bar. In contrast, the lack of responsiveness to bar instabilities (a larger Toomre-Q parameter) due to a less compact nature plays an important role in generating unbarred disk galaxies with stellar mass log$(M_*/M_\odot)<10.8$. Moreover, short bars generally form at a similar time to normal bars, during which they either grow mildly or contract significantly. The fact that IllustrisTNG simulations may produce too many short-bar galaxies indicates that the dynamical properties of the central regions in IllustrisTNG galaxies are less likely to be affected by external factors, such as mergers and gas inflows.

Rotation-powered pulsars represent the main class of identified gamma-ray sources in the Galaxy. The wealth of observational data collected by the AGILE and Fermi gamma-ray space telescopes in the GeV range, and by ground-based Cherenkov telescopes in the TeV band provide invaluable insights into how relativistic plasmas dissipate and accelerate particles. Decoding the information contained in the gamma-ray pulses profile is an important step to understand how pulsars work. In this study, we aim at putting an ab initio plasma model of pulsar magnetospheres to the test, in light of the most recent gamma-ray observations in the GeV and TeV bands. To this end, we present of a new series of global particle-in-cell simulations of an inclined pulsar magnetosphere. High-quality synthetic pulse profiles in the synchrotron and inverse Compton channels are reconstructed to study in greater details their morphology and their energy dependence. We also perform a fit of observed lightcurves with the model, using the third Fermi-LAT gamma-ray pulsar catalog. Reconnection in the wind current sheet powers synchrotron and inverse Compton emission. The modeled pulse profiles reproduce some of the salient features of observed gamma-ray pulsars, including the mysterious Vela-like lightcurves, such as: the generic double-peaked structure, the presence of a bridge or third peak in between the main pulses, the pulse narrowing with increasing energy. The bolometric synchrotron radiative efficiency is strictly limited by the reconnection rate. Our global kinetic simulations are able to match observed pulse profiles. Such direct comparisons will help drive and focus future simulation developments.

Making inferences about physical properties of the Universe requires knowledge of the data likelihood. A Gaussian distribution is commonly assumed for the uncertainties with a covariance matrix estimated from a set of simulations. The noise in such covariance estimates causes two problems: it distorts the width of the parameter contours, and it adds scatter to the location of those contours which is not captured by the widths themselves. For non-Gaussian likelihoods, an approximation may be derived via Simulation-Based Inference (SBI). It is often implicitly assumed that parameter constraints from SBI analyses, which do not use covariance matrices, are not affected by the same problems as parameter estimation with a covariance matrix estimated from simulations. We investigate whether SBI suffers from effects similar to those of covariance estimation in Gaussian likelihoods. We use Neural Posterior and Likelihood Estimation with continuous and masked autoregressive normalizing flows for density estimation. We fit our approximate posterior models to simulations drawn from a Gaussian linear model, so that the SBI result can be compared to the true posterior. We test linear and neural network based compression, demonstrating that neither methods circumvent the issues of covariance estimation. SBI suffers an inflation of posterior variance that is equal or greater than the analytical result in covariance estimation for Gaussian likelihoods for the same number of simulations. The assumption that SBI requires a smaller number of simulations than covariance estimation for a Gaussian likelihood analysis is inaccurate. The limitations of traditional likelihood analysis with simulation-based covariance remain for SBI with a finite simulation budget. Despite these issues, we show that SBI correctly draws the true posterior contour given enough simulations.

We compare two candidate nonlinearities for regulating the solar cycle within the Babcock-Leighton paradigm: tilt quenching (whereby the tilt of active regions is reduced in stronger cycles) and latitude quenching (whereby flux emerges at higher latitudes in stronger solar cycles). Digitized historical observations are used to build a database of individual magnetic plage regions from 1923 to 1985. The regions are selected by thresholding in Ca II K synoptic maps, with polarities constrained using Mount Wilson Observatory sunspot measurements. The resulting data show weak evidence for tilt quenching, but much stronger evidence for latitude-quenching. Further, we use proxy observations of the polar field from faculae to construct a best-fit surface flux transport model driven by our database of emerging regions. A better fit is obtained when the sunspot measurements are used, compared to a reference model where all polarities are filled using Hale's Law. The optimization suggests clearly that the "dynamo effectivity range" of the Sun during this period should be less than 10 degrees; this is also consistent with latitude quenching being dominant over tilt quenching.

We present an analysis of LAMOST J171013.211+532646.04 (hereafter J1710), a binary system comprising a hot subdwarf B star (sdB) and a white dwarf (WD) companion. Multi-epoch spectroscopy reveals an orbital period of 109.20279 minutes, consistent with TESS and ZTF photometric data, marking it as the sixth detached system known to harbor a WD companion with a period less than two hours. J1710 is remarkably close to Earth, situated at a distance of only \(350.68^{+4.20}_{-4.21} \, \mathrm{pc}\), with a GAIA G-band magnitude of 12.59, rendering it conducive for continuous observations. The spectral temperature is around 25164 K, in agreement with SED fitting results (\(25301^{+839}_{-743} \, \mathrm{K}\)). The TESS light curve displays ellipsoidal variation and Doppler beaming without eclipsing features. Through fitting the TESS light curve using the Wilson-Devinney code, we determined the masses for the sdB (\(M_1 = 0.44^{+0.06}_{-0.07} \, M_{\odot}\)) and the compact object (\(M_2 = 0.54^{+0.10}_{-0.07} \, M_{\odot}\)), with the compact object likely being a WD. Furthermore, MESA models suggest that the sdB, with a helium core mass of 0.431 \(M_{\odot}\) and a hydrogen envelope mass of \(1.3 \times 10^{-3}\, M_{\odot}\), is in the early helium main-sequence phase. The MESA binary evolution shows that the J1710 system is expected to evolve into a double white dwarf system, making it an important source of low-frequency gravitational waves.

After a merger of two massive black holes (MBHs), the remnant receives a gravitational wave (GW) recoil kick that can have a strong effect on its future evolution. The magnitude of the kick ($v_\mathrm{recoil}$) depends on the mass ratio and the alignment of the spins and orbital angular momenta, therefore on the previous evolution of the MBHs. We investigate the cosmic effect of GW recoil by running for the first time a high-resolution cosmological simulation including GW recoil that depends on the MBH spins (evolved through accretion and mergers), masses and dynamics computed self-consistently. We also run a twin simulation without GW recoil. The simulations are run down to $z=4.4$. We find that GW recoil reduces the growth of merger remnants, and can have a significant effect on the MBH-galaxy correlations and the merger rate. We find large recoil kicks across all galaxy masses in the simulation, up to a few $10^{11}\,\rm M_\odot$. The effect of recoil can be significant even if the MBHs are embedded in a rotationally supported gaseous structure. We investigate the dynamics of recoiling MBHs and find that MBHs remain in the centre of the host galaxy for low $v_\mathrm{recoil}/v_\mathrm{esc}$ and escape rapidly for high $v_\mathrm{recoil}/v_\mathrm{esc}$. Only if $v_\mathrm{recoil}$ is comparable to $v_\mathrm{esc}$ the MBHs escape the central region of the galaxy but might remain as wandering MBHs until the end of the simulation. Recoiling MBHs are a significant fraction of the wandering MBH population. Although the dynamics of recoiling MBHs may be complex, some retain their initial radial orbits but are difficult to discern from other wandering MBHs on radial orbits. Others scatter with the halo substructure or circularise in the asymmetric potential. Our work highlights the importance of including GW recoil in cosmological simulation models.

We present a catalogue of photometric redshifts for galaxies from DESI Legacy Imaging Surveys, which includes $\sim0.18$ billion sources covering 14,000 ${\rm deg}^2$. The photometric redshifts, along with their uncertainties, are estimated through galaxy images in three optical bands ($g$, $r$ and $z$) from DESI and two near-infrared bands ($W1$ and $W2$) from WISE using a Bayesian Neural Network (BNN). The training of BNN is performed by above images and their corresponding spectroscopic redshifts given in DESI Early Data Release (EDR). Our results show that categorizing galaxies into individual groups based on their inherent characteristics and estimating their photo-$z$s within their group separately can effectively improve the performance. Specifically, the galaxies are categorized into four distinct groups based on DESI's target selection criteria: Bright Galaxy Sample (BGS), Luminous Red Galaxies (LRG), Emission Line Galaxies (ELG) and a group comprising the remaining sources, referred to as NON. As measured by outliers of $|\Delta z| > 0.15 (1 + z_{\rm true})$, accuracy $\sigma_{\rm NMAD}$ and mean uncertainty $\overline{E}$ for BNN, we achieve low outlier percentage, high accuracy and low uncertainty: 0.14%, 0.018 and 0.0212 for BGS and 0.45%, 0.026 and 0.0293 for LRG respectively, surpassing results without categorization. However, the photo-$z$s for ELG cannot be reliably estimated, showing result of $>15\%$, $\sim0.1$ and $\sim0.1$ irrespective of training strategy. On the other hand, NON sources can reach 1.9%, 0.039 and 0.0445 when a magnitude cut of $z<21.3$ is applied. Our findings demonstrate that estimating photo-$z$s directly from galaxy images is significantly potential, and to achieve high-quality photo-$z$ measurement for ongoing and future large-scale imaging survey, it is sensible to implement categorization of sources based on their characteristics.

Deep-learning-based methods have been favored in astrophysics owing to their adaptability and remarkable performance and have been applied to the task of the classification of real and bogus transients. Different from most existing approaches which necessitate massive yet expensive annotated data, We aim to leverage training samples with only 1000 labels available to discover real sources that vary in brightness over time in the early stage of the WFST 6-year survey. Methods. We present a novel deep-learning method that combines active learning and semi-supervised learning to construct a competitive real/bogus classifier. Our method incorporates an active learning stage, where we actively select the most informative or uncertain samples for annotation. This stage aims to achieve higher model performance by leveraging fewer labeled samples, thus reducing annotation costs and improving the overall learning process efficiency. Furthermore, our approach involves a semi-supervised learning stage that exploits the unlabeled data to enhance the model's performance and achieve superior results compared to using only the limited labeled data.

The paper is devoted to the analysis of the influence of the galactic bar on the nature of the orbital motion (chaotic or regular) of globular clusters in the central region of the Galaxy with a radius of 3.5 kpc, which are subject to the greatest influence of the bar. The sample includes 45 globular clusters. To form the 6D phase space required for integrating the orbits, the most accurate astrometric data to date from the Gaia satellite (Vasiliev, Baumgardt, 2021) were used, as well as new refined average distances (Baumgardt, Vasiliev, 2021). The orbits of the globular clusters were obtained both in an axisymmetric potential and in a potential including the bar. The following, most realistic, bar parameters were adopted: mass $10^{10} M_\odot$, semi-major axis length 5 kpc, bar axis rotation angle 25$^o$, angular rotation velocity 40 km s$^{-1}$ kpc$^{-1}$. The analysis of the chaoticity/regularity of the orbital motion in both potentials was carried out using one of the most effective methods, namely, the frequency method, which consists in calculating the drift of fundamental frequencies. As a result, the influence of the bar on the dynamics of each GC of the sample was assessed. It is established that 8 GCs changed regular dynamics to chaotic under the influence of the bar, and 9 GCs changed chaotic dynamics to regular one.

The supersonic flow motions associated with infall of baryonic gas toward sheets and filaments, as well as cluster mergers, produces large-scale shock waves. The shocks associated with galaxy clusters can be classified mainly into two categories: internal shocks appear in the hot intracluster medium within the viral radius, and external accretion shocks form in the outer cold region well outside of the virial radius. Cosmic-ray (CR) electrons and/or protons accelerated by these shocks are expected to produce gamma-rays through inverse-Compton scatterings (ICS) or inelastic $pp$ collisions respectively. Recent studies have found a spatially extended GeV source within the virial radius, consistent with the internal shock origin. Here we report the detection of a new GeV source at a distance of about 2.8$^\circ$ from the center of the Coma cluster through the analysis of 16.2 years of Fermi-LAT data. The hard spectrum of the source, in agreement with the ICS origin, and its location in a large-scale filament of galaxies points to the external accretion shock origin. The gamma-ray ($0.1-10^3$ GeV) luminosity of the source, $1.4\times 10^{42}~ {\rm erg~s^{-1}}$, suggests that a fraction $\sim 10^{-3}$ of the kinetic energy flux through the shock-surface is transferred to relativistic CR electrons.

Understanding how galaxies form and evolve is at the heart of modern astronomy. With the advent of large-scale surveys and simulations, remarkable progress has been made in the last few decades. Despite this, the physical processes behind the phenomena, and particularly their importance, remain far from known, as correlations have primarily been established rather than the underlying causality. We address this challenge by applying the causal inference framework. Specifically, we tackle the fundamental open question of whether galaxy formation and evolution depends more on nature (i.e., internal processes) or nurture (i.e., external processes), by estimating the causal effect of environment on star-formation rate in the IllustrisTNG simulations. To do so, we develop a comprehensive causal model and employ cutting-edge techniques from epidemiology to overcome the long-standing problem of disentangling nature and nurture. We find that the causal effect is negative and substantial, with environment suppressing the SFR by a maximal factor of $\sim100$. While the overall effect at $z=0$ is negative, in the early universe, environment is discovered to have a positive impact, boosting star formation by a factor of $\sim10$ at $z\sim1$ and by even greater amounts at higher redshifts. Furthermore, we show that: (i) nature also plays an important role, as ignoring it underestimates the causal effect in intermediate-density environments by a factor of $\sim2$, (ii) controlling for the stellar mass at a snapshot in time, as is common in the literature, is not only insufficient to disentangle nature and nurture but actually has an adverse effect, though (iii) stellar mass is an adequate proxy of the effects of nature. Finally, this work may prove a useful blueprint for extracting causal insights in other fields that deal with dynamical systems with closed feedback loops, such as the Earth's climate.

The Taurus region is one of the most extensively studied star-forming regions. Surveys indicate that the young stars in this region are comprised of Young Stellar Objects (YSOs) that cluster in groups associated with the molecular cloud (Grouped Young Stellar Objects, GYSOs), and some older ones that are sparsely distributed throughout the region (Distributed Young Stellar Objects, DYSOs). To bridge the age gap between the GYSOs ($\le$5 Myr) and the DYSOs (10-20 Myr), we conducted a survey to search for new YSOs in this direction. Based on infrared color excesses and Li I absorption lines, we identified 145 new YSOs. Combining these with the previously identified GYSOs and DYSOs, we constructed a sample of 519 YSOs that encompasses the entire region. Subsequently, we calculated the ages of the sample based on their proximity to the local bubble. The age versus Distance to the Local Bubble ($D_{\rm LB}$) relationship for the DYSOs shows a clear trend: the farther they are from the Local Bubble, the younger they are, which is consistent with the supernovae-driven formation scenario of the Local Bubble. The GYSOs also exhibit a mild age versus $D_{\rm LB}$ trend. However, they are significantly younger and are mostly confined to distances of 120 to 220 pc. Considering their distribution in the age versus $D_{\rm LB}$ space is well separated from the older DYSOs, they may also be products of the Local Bubble but formed in more recent and localized events.

The stochastic gravitational-wave background originating from cosmic sources contains vital information about the early universe. In this work, we comprehensively study the cross-correlations between the energy-density anisotropies in scalar-induced gravitational waves (SIGWs) and the temperature anisotropies and polarization in the cosmic microwave background (CMB). In our analysis of the angular power spectra for these cross-correlations, we consider all contributions of the local-type primordial non-Gaussianity $f_{\mathrm{NL}}$ and $g_{\mathrm{NL}}$ that can lead to large anisotropies. We show that the angular power spectra are highly sensitive to primordial non-Gaussianity. Furthermore, we project the sensitivity of future gravitational-wave detectors to detect such signals and, consequently, measure the primordial non-Gaussianity.

Recent observations reveal a rapid dust build-up in high-redshift galaxies (z > 4), challenging current models of galaxy formation. While our understanding of dust production and destruction in the interstellar medium (ISM) is advancing, probing baryonic processes in the early Universe remains a complex task. We characterize the evolution of 98 z~5 star-forming galaxies observed as part of the ALPINE survey by constraining the physical processes underpinning the gas and dust production, consumption, and destruction in their ISM. We make use of chemical evolution models to simultaneously reproduce the observed dust and gas content. For each galaxy, we estimate initial gas mass, inflows and outflows, and efficiencies of dust growth and destruction. We test the models with the canonical Chabrier and top-heavy initial mass functions (IMFs), with the latter enabling rapid dust production on shorter timescales. Our models successfully reproduce gas and dust content in older galaxies (> 600 Myr) regardless of the IMF, with Type II SNe as the primary dust source and no dust growth in ISM with moderate inflow of primordial gas. In case of intermediate-age galaxies (300 - 600 Myr), we reproduce the gas and dust content through Type II SNe and dust growth in ISM, though we observe an over-prediction of dust mass in older galaxies, potentially indicating an unaccounted dust destruction mechanism and/or an overestimation of the observed dust masses. The number of young galaxies (< 300 Myr) reproduced, increases for models assuming top-heavy IMF but with maximal prescriptions of dust production. Galactic outflows are necessary to reproduce observed gas and dust masses. The Chabrier IMF models reproduce 65% of galaxies, while top-heavy IMF models improve this to 93%, easing tensions with observations. Upcoming JWST data will refine these models by resolving degeneracies in intrinsic galaxy properties.

The solar radiative zone spin down problem is treated without internal gravity waves, IGW. The shear-only model of the Reynolds stresses adopted thus far in all calculations, is substituted by a shear vorticity model. The latter component is shown to play a role analogous to the IGW and with a similar time scale of about 107years.

We study the formation of solitons inside scalar-field dark matter halos with a non-polynomial self-interaction potential. We consider a self-interaction potential that is quartic in the scalar field in the low-density regime but saturates at large densities. This mimics the behaviour of axion monodromy potentials. We concentrate on the semi-classical regime, where the de Broglie wavelength is much smaller than the size of the system. We find that depending on the strength and scale of the self-interactions, the system can form solitons of the Thomas-Fermi type (dominated by self-interactions) or of the Fuzzy Dark Matter type (dominated by the quantum pressure). The system can also display transitions from a Thomas-Fermi soliton to a Fuzzy Dark Matter soliton as the former becomes unstable. We show that these behaviours can be understood from a simple Gaussian ansatz. We find that even in cases where the self-interactions are always subdominant they can play a critical role, by providing a small density boost that is enough to generate the seed for the formation of a Fuzzy Dark Matter soliton at much later times. We also point out that the intuition derived from a hydrodynamical picture can be misleading in regimes where wave effects are important.

A plethora of inflationary models can produce interesting small-scale phenomenology, such as enhanced scalar fluctuations leading to primordial black hole (PBH) production and large scalar-induced GW. Nevertheless, good models must simultaneously explain current observations on all scales. In this work, we showcase our methodology to establish the small-scale phenomenology of inflationary models on firm grounds. We consider the case of hybrid $\alpha$-attractors, and focus on a reduced parameter space featuring the two potential parameters which roughly determine the position of the peak in the scalar power spectrum, $\mathcal{P}_\zeta$, and its amplitude. We first constrain the parameter space by comparing the large-scale predictions for $\mathcal{P}_\zeta$ with current CMB anisotropies measurements and upper limits on $\mu$-distortions. We take into account uncertainties due to the reheating phase, and observe that the parameter-space area compatible with large-scale constraints shrinks for extended reheating stages. We then move to smaller scales, where we find that non-Gaussianity at peak scales is of the local type and has amplitude $f_\text{NL}\sim \mathcal{O}(0.1)$. This ensures that non-linear effects are subdominant, motivating us to employ the tree-level $\mathcal{P}_\zeta$ to compute the abundance of PBHs and the spectrum of induced GWs for models consistent with large-scale tests. The former allows us to further constrain the parameter space, by excluding models which over-produce PBHs. We find that a subset of viable models can lead to significant production of PBHs, and a fraction of these is within reach for LISA, having a signal-to-noise ratio larger than that of astrophysical foregrounds. Our first-of-its-kind study systematically combines tests at different scales, and exploits the synergy between cosmological observations and theoretical consistency requirements.

The precise measurement of the sky-averaged HI absorption signal between 50 and 200 MHz is the primary goal of global 21-cm cosmology. This measurement has the potential to unravel the underlying physics of cosmic structure formation and evolution during the Cosmic Dawn. It is, however, hindered by various non-smooth, frequency-dependent effects, whose structures resemble those of the signal. One such effect is the leakage of polarised foregrounds into the measured intensity signal: polarised foreground emission undergoes Faraday rotation as it passes through the magnetic fields of the interstellar medium, imprinting a chromatic structure in the relevant frequency range which complicates the extraction of the cosmological HI absorption feature. We investigate the effect of polarised Galactic foregrounds on extracting the global 21-cm signal from simulated data using REACH's data analysis pipeline; the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) is an experiment designed to detect the sky-averaged 21-cm HI signal from the early Universe using physically informed models. Using the REACH pipeline, we successfully recover an injected global 21-cm signal with an amplitude of approximately 0.16 K, centred between 80 and 120 MHz, achieving a low root-mean-square error (less than 30\% of the injected signal strength) in all the tested cases. This includes scenarios with simulated polarised Galactic diffuse emissions and polarised point source emissions, provided the overall polarisation fraction is below $\sim 3\%$. The linear mixing of contamination, caused by the superposition of multiple patches with varying strengths of Faraday rotation, produces patterns that are more distinct from the global signal. This distinction makes global signal recovery easier compared to contamination resulting from a single, slow oscillation pattern.

[Abridged] Our view of dust in primordial galaxies is limited towards a few tens of z~7 galaxies, pre-selected from UV-optical observations, and are thus not necessarily representative of the bulk of the sources at these redshifts. In this work, we aim at constraining the dust properties of galaxies at 6<z<12 by making the most of the A3COSMOS database in the JADES/GOODS-South field. We stacked ALMA band 6 and 7 observations of 4464 galaxies covered by the A3COSMOS database and used the measurements as constraints to perform UV-to-FIR SED modelling. We obtain tentative signals for the brightest UV galaxies (MUV<-19mag) as well as for the most massive ones (log Mstar>9) at 6<z<7, and upper limits for fainter (MUV>-19mag), lower mass sources (log Mstar<9), and at higher redshift (z>7). Fitting these 6<z<7 galaxies with ALMA constraints results in lower star formation rates (-0.4dex) and FUV attenuation (-0.5mag) for galaxies with log Mstar>8, compared to the fit without FIR. We extend the LIR vs MUV relation down to MUV=-19mag and show a tentative breakdown of the relation at fainter UV magnitudes. The positions of the JADES z~6.5 sample on the IRX versus beta and IRX versus Mstar diagrams are consistent with those of the ALPINE (z~5.5) and REBELS (z~6.5) samples, suggesting that the dust composition and content of our mass-selected sample are similar to these UV-selected galaxies. Extending our analysis of the infrared properties to z>7 galaxies, we find a non-evolution of beta with redshift in the MUV range probed by our sample (-17.24+/-0.62) and highlight the fact that samples from the literature are not representative of the bulk of galaxy populations at z>6. We confirm a linear relation between AV and 1/sSFR with a flatter slope than previously reported due to the use of ALMA constraints. Our results suggest that rapid and significant dust production has already happened by z~7.

Many questions posed in the Astro2020 Decadal survey in both the New Messengers and New Physics and the Cosmic Ecosystems science themes require a gamma-ray mission with capabilities exceeding those of existing (e.g. Fermi, Swift) and planned (e.g. COSI) observatories. ComPair, the Compton Pair telescope, is a prototype of such a next-generation gamma-ray mission. It had its inaugural balloon flight from Ft. Sumner, New Mexico in August 2023. To continue the goals of the ComPair project to develop technologies that will enable a future gamma-ray mission, the next generation of ComPair (ComPair-2) will be upgraded to increase the sensitivity and low-energy transient capabilities of the instrument. These advancements are enabled by AstroPix, a silicon monolithic active pixel sensor, in the tracker and custom dual-gain silicon photomultipliers and front-end electronics in the calorimeter. This effort builds on design work for the All-sky Medium Energy Gamma-ray Observatory eXplorer (AMEGO-X) concept that was submitted the 2021 MIDEX Announcement of Opportunity. Here we describe the ComPair-2 prototype design and integration and testing plans to advance the readiness level of these novel technologies.

The theory of pebble accretion (PA) has become a popular planet formation theory during the past decade. However, PA studies generally rely on large planetary embryos, much larger than those expected from the streaming instability. This study analyses the formation of terrestrial planets around stars with masses ranging from 0.09 to 1.00 M$_\odot$ through PA, starting with small planetesimals with radii between 175 and 450 km. We perform numerical simulations, using a modified version of the N-body simulator SyMBA, which includes pebble accretion, type I and II migration, and eccentricity and inclination damping. Two prescriptions for the pebble accretion rate are analysed. We find that Earth-mass planets are consistently formed around 0.49, 0.70 and 1.00 M$_\odot$ stars, irrespective of the pebble accretion model. However, Earth-like planets seldom remain in the habitable zone, for they rapidly migrate to the inner edge of the disc. Furthermore, we find that pebble accretion onto small planetesimals cannot produce Earth-mass planets around 0.09 and 0.20 M$_\odot$, challenging the proposed narrative of the formation of the TRAPPIST-1 system. Further research is needed to determine if models with a lower pebble mass flux, or additional migration traps, could produce solar-system-like planetary systems in which Earth-like planets remain in the habitable zone.

Upcoming deep optical surveys, such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), will scan the sky to unprecedented depths, detecting billions of galaxies. However, this amount of detections will lead to the apparent superposition of galaxies in the images, a phenomenon known as blending, that can affect the accurate measurement of individual galaxy properties. In particular, galaxy shapes play a crucial role in estimating the masses of large-scale structures, such as galaxy clusters, through weak gravitational lensing. This proceeding introduces a new catalog matching algorithm, friendly, designed for detecting and characterizing blends in simulated LSST data for the Dark Energy Science Collaboration (DESC) Data Challenge 2. The aim of this algorithm is to combine several matching procedures, as well as a probabilistic method to quantify blended systems. By removing the resulted 27% of galaxies affected by blending from the dataset, we demonstrate that the amplitude of the excess surface mass density weak lensing profile - potentially biased low due to blending - may be partially corrected.

Robust inference of galaxy stellar masses from photometry is crucial for constraints on galaxy assembly across cosmic time. Here, we test a commonly-used Spectral Energy Distribution (SED) fitting code, using simulated galaxies from the SPHINX20 cosmological radiation hydrodynamics simulation, with JWST NIRCam photometry forward-modelled with radiative transfer. Fitting the synthetic photometry with various star formation history models, we show that recovered stellar masses are, encouragingly, generally robust to within a factor of ~3 for galaxies in the range M*~10^7-10^9M_sol at z=5-10. These results are in stark contrast to recent work claiming that stellar masses can be underestimated by as much as an order of magnitude in these mass and redshift ranges. However, while >90% of masses are recovered to within 0.5dex, there are notable systematic trends, with stellar masses typically overestimated for low-mass galaxies (M*<~10^8M_sol) and slightly underestimated for high-mass galaxies (M*>~10^9M_sol). We demonstrate that these trends arise due to the SED fitting code poorly modelling the impact of strong emission lines on broadband photometry. These systematic trends, which exist for all star formation history parametrisations tested, have a tilting effect on the inferred stellar mass function, with number densities of massive galaxies underestimated (particularly at the lowest redshifts studied) and number densities of lower-mass galaxies typically overestimated. Overall, this work suggests that we should be optimistic about our ability to infer the masses of high-z galaxies observed with JWST (notwithstanding contamination from AGN) but careful when modelling the impact of strong emission lines on broadband photometry.

A multi-wavelength, multi-scale study of the Mon R2 hub-filament system (HFS) reveals a spiral structure, with the central hub containing more mass than its filaments. ALMA C$^{18}$O(1-0) emission reveals several accreting filaments connected to a molecular ring (size $\sim$0.18 pc $\times$ 0.26 pc). The molecular ring surrounds the infrared (IR) ring (size $\sim$0.12 pc $\times$ 0.16 pc), which is not usually observed. The IR ring encircles IR dark regions and a population of embedded near-IR sources, including the massive stars IRS 1 and IRS 2. ALMA HNC(3-2) line data reveal a mirrored B-shaped feature (extent $\sim$19000 AU $\times$ 39000 AU) toward the eastern part of the molecular ring, suggesting expansion at $\sim$2.25 km s$^{-1}$. Distinct HNC sub-structures in both redshifted and blueshifted velocity components are investigated toward the B-shaped feature. The presence of these braid-like substructures in each velocity component strongly suggests instability in photon-dominated regions. A dusty shell-like feature (extent $\sim$0.04 pc $\times$ 0.07 pc; mass $\sim$7 M$_{\odot}$) hosting IRS 1 is identified in the ALMA 1.14 mm continuum map, centered toward the base of the B-shaped feature. The IR and dense molecular rings are likely shaped by feedback from massive stars, driven by high pressure values between 10$^{-8}$-10$^{-10}$ dynes cm$^{-2}$, observed within a 1 pc range of the B0 ZAMS star powering the ultracompact HII region. Overall, these outcomes support that the Mon R2 HFS transitioned from IR-quiet to IR-bright, driven by the interaction between gas accretion and feedback from massive stars.

GRB 221009A is one of the brightest transients ever observed with the highest peak gamma-ray flux for a gamma-ray burst (GRB). A type Ic-BL supernova (SN), SN 2022xiw, was definitively detected in late-time JWST spectroscopy (t = 195 days, observer-frame). However, photometric studies have found SN 2022xiw to be less luminous (10-70%) than the canonical GRB-SN, SN 1998bw. We present late-time Hubble Space Telescope (HST)/WFC3 and JWST/NIRCam imaging of the afterglow and host galaxy of GRB 221009A at t ~ 185, 277, and 345 days post-trigger. Our joint archival ground, HST, and JWST light curve fits show strong support for a break in the light curve decay slope at t = 50 +/- 10 days (observer-frame) and a supernova at $1.4^{+0.37}_{-0.40} \times$ the optical/NIR flux of SN 1998bw. This break is consistent with an interpretation as a jet break when requiring slow-cooling electrons in a wind medium with the electron energy spectral index, p > 2, and $\nu_m < \nu_c$. Our light curve and joint HST/JWST spectral energy distribution (SED) also show evidence for the late-time emergence of a bluer component in addition to the fading afterglow and supernova. We find consistency with the interpretations that this source is either a young, massive, low-metallicity star cluster or a scattered light echo of the afterglow with a SED shape of $f_{\nu} \propto \nu^{2.0\pm1.0}$.

Metallicity plays a crucial role in the evolution of massive stars and their final core-collapse supernova (CCSN) explosions. Integral-field-unit (IFU) spectroscopy can provide a spatially resolved view of SN host galaxies and serve as a powerful tool to study SN metallicities. Early transient surveys targeted bright galaxies with high star formation and SN rates; as a result, the discovered SNe are significantly biased toward high metallicities. More recently, the untargeted, wide-field transient surveys, such as ASAS-SN and ZTF, have discovered a large number of SNe without such a bias. In this work, we construct a large and unbiased sample of SNe discovered by (quasi-) untargted searches, consisting of 209 SNe of Types II (with unknown subtypes), IIP, IIn, IIb, Ib and Ic at z $\leq$ 0.02 with VLT/MUSE observations. This is currently the largest CCSN sample with IFU observations. With the strong-line method, we reveal the spatially-resolved metallicity maps of the SN host galaxies and acquire accurate metallicity measurements for the SN sites. Our results show that the SN metallicities range from 12 + log(O/H) = 8.1 to 8.7 dex, and the metallicity distributions for different SN types are very close to each other, with mean and median values of 8.4 - 8.5 dex. We carefully analysed the stochastic sampling effect, showing that our large sample size narrows the 1$\sigma$ uncertainty down to only 0.05 dex. The apparent metallicity differences among SN types are all within 3$\sigma$ uncertainties and the metallicity distributions for different SN types are all consistent with being randomly drawn from the same reference distribution. This suggests that metallicity plays a minor role in the origin of different CCSN types and some other metallicity-insensitive processes, such as binary interaction, dominate the distinction of CCSN types.

Are we alone? It is a compelling question that human beings have confronted for centuries. The search for extraterrestrial life is a broad range of quests for finding the simple forms of life up to intelligent beings in the Universe. The plausible assumption is that there is a chance that intelligent life will followed by advanced civilization equipped or even dominated by artificial intelligence (AI). In this work, we categorize the advanced civilizations (on an equal footing, an AI-dominated civilization) on the Kardashev scale. We propose a new scale known as space exploration distance to measure civilization advancement. We propose a relation between this length and the Kardashev scale. Then, we suggest the idea that advanced civilizations will use primordial black holes as sources of harvesting energy. We calculate the energy harvested by calculating the space exploration distance. Finally, we propose an observational method to detect the possibility of extraterrestrial AI using Dyson spheres-like structures around primordial black holes in the Milky Way and other galaxies.

Observational campaigns hunting the elusive reservoirs of cold gas in the host galaxies of quasars at the Epoch of Reionization (EoR) are crucial to study the formation and evolution of the first massive systems at early epochs. We present new Northern Extended Millimetre Array (NOEMA) observations tracing CO(6--5), CO(7--6) emission lines, and the underlying continuum in five of the eight quasars at redshift $z>7$ known to date, thus completing the survey of the cold molecular gas reservoir in the host galaxies of the first quasars. Combining NOEMA observations with archival Atacama Large Millimeter Array (ALMA) data available, we model the far-infrared spectral energy distribution with a modified blackbody to measure dust properties and star formation rates. We use CO and [CII] lines to derive molecular gas masses, which we compare with results from semi-analytical models and observations of galaxies at different epochs. No statistically significant detection of CO emission lines was reported for the five quasars in this sample, resulting in a relatively low amount of cold molecular gas in the host when compared with galaxies at later epochs. Nonetheless, gas-to-dust ratios are consistent with the local value, suggesting that the scaling relation between dust and cold gas holds up to $z>7$. Quasars at the EoR show star formation efficiencies which are among the highest observed so far, but comparable with that observed in luminous quasar at Cosmic Noon and that predicted for the brightest ($L_{bol}>3\times10^{46}$ erg s$^{-1}$) quasar objects drawn from the semi-analytical model GAEA. Quasar host galaxies at the EoR are undergoing an intense phase of star formation, which suggests a strong coupling between the luminous phase of the quasar and the rapid growth of the host.

Discoveries of giant planet candidates orbiting white dwarf stars and the demonstrated capabilities of the James Webb Space Telescope bring the possibility of detecting rocky planets in the habitable zones of white dwarfs into pertinent focus. We present simulations of an aqua planet with an Earth-like atmospheric composition and incident stellar insolation orbiting in the habitable zone of two different types of stars - a 5000 K white dwarf and main-sequence K-dwarf star Kepler-62 with a similar effective temperature - and identify the mechanisms responsible for the two differing planetary climates. The synchronously-rotating white dwarf planet's global mean surface temperature is 25 K higher than that of the synchronously-rotating planet orbiting Kepler-62, due to its much faster (10-hr) rotation and orbital period. This ultra-fast rotation generates strong zonal winds and meridional flux of zonal momentum, stretching out and homogenizing the scale of atmospheric circulation, and preventing an equivalent build-up of thick, liquid water clouds on the dayside of the planet compared to the synchronous planet orbiting Kepler-62, while also transporting heat equatorward from higher latitudes. White dwarfs may therefore present amenable environments for life on planets formed within or migrated to their habitable zones, generating warmer surface environments than those of planets with main-sequence hosts to compensate for an ever shrinking incident stellar flux.

We report forecasted constraints on warm inflation in the light of future cosmic microwave background (CMB) surveys, with data expected to be available in the coming decade. These observations could finally give us the missing information necessary to unveil the production of gravitational waves during inflation, reflected by detecting a non-zero tensor-to-scalar ratio crucial to the B-mode power spectrum of the CMB. We consider the impact of three future surveys, namely the CMB-S4, Simons Observatory, and the space-borne $\textit{LiteBIRD}$, in restricting the parameter space of four typical warm inflationary models in the context of a quartic potential, which is well motivated theoretically. We find that all three surveys significantly improve the models' parameter space, compared to recent results obtained with current $\textit{Planck}$+BICEP/Keck Array data. Moreover, the combination of ground-based and space-borne (CMB-S4+$\textit{LiteBIRD}$) tightens the constraints so that we expect to distinguish even better warm inflation scenarios. This result becomes clear when we compare the models' predictions with a $\Lambda$CDM+$r$ forecast, compatible with $r=0$, in which one of them already becomes excluded by data.

In this work, we analyze the characteristics of electromagnetic (EM) radiation associated with scalar and axion field oscillations in different background field setups. Because the scalar field and axion field have different parity and couple with the EM field in different forms, the EM signals generated by the scalar and axion can be used to distinguish them. More interestingly, resonance effects amplify the difference between the two fields and consequent EM signal strength, which helps us distinguish and detect them in future observations.

We discuss a recently proposed interpretation of the signal detected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) as due to a relic stochastic background of primordial gravitons, produced in the context of the string cosmology pre-big bang scenario. We show that such interpretation cannot be reconciled with a phenomenologically viable minimal version of such scenario, while it can be allowed if one considers an equally viable but generalised, non-minimal version of pre-big bang evolution. Maintaining the $S$-duality symmetry throughout the high-curvature string phase is possible although somewhat disfavoured. The implications of this non-minimal scenario for the power spectrum of curvature perturbations are also briefly discussed.

In the post-Newtonian era, the Eddington-inspired Born-Infeld (EiBI) theory, considered as an improved modification of the Einsteinian general relativity formalism in the weak field regime (non-relativistic), has enabled us to study the dynamics of dense astroobjects in light of the modified gravitational effects. This EiBI theory imparts a new shape to the usual gravitational Poisson equation through the addition of a cosmological correction factor, termed as the EiBI gravity parameter. A systematic inclusion of this gravity in the basic structure equation could lead to a realistic picture of the existing solar models free from any end-stage singularity. A theoretic model is accordingly proposed to investigate the effect of the EiBI gravity on the Gravito-Electrostatic Sheath (GES) formalism of the equilibrium solar plasma structure. This study shows that the GES-based solar plasma dynamics is noticeably modified against the previously reported Newtonian GES-model studies. An equilibrium bounded solution for the solar self-gravity shows the EiBI-modified solar surface boundary (SSB) to exist at a new helio-centric radial location $\xi = 4$ (on the Jeansean scale). It is found that the EiBI gravity shifts the present SSB outwards by 14.28% relative to the original Newtonian SSB. The EiBI-modified gravity effects on diverse relevant solar parameters, such as the gravito-electrostatic potentials, fields, and Mach numbers, are illustratively analyzed. It is anticipated that our analyses could be applied further to see the solar plasma equilibrium and fluctuation dynamics in realistically modified post-Newtonian gravity environments on both the bounded (interior) and unbounded (exterior) solar plasma scales.

We analyze the computation of $n$-point correlation functions in de Sitter spacetime, including loop corrections, using the wavefunction of the universe approach. This method consists of two stages employing distinct Feynman rules. First, one must compute the wavefunction coefficients using interactions as vertices. Then, in the second stage, one computes correlation functions using wavefunction coefficients as vertices. For massless fields, loop corrections in the first stage are free of infrared (IR) divergences, which leads to the question of how this matches the well-known IR behavior of correlators obtained via other methods. By considering a scalar field with an arbitrary potential, we compute $n$-point correlation functions to first order in the potential but to all orders in loops. We find that, although loop integrals in the first stage are indeed IR convergent, the second procedure reintroduces the IR divergence. We discuss how this induces renormalization of the interaction potential such that the final result combining both steps exactly matches the form of $n$-point functions previously calculated with other methods.

Superradiance, the process by which waves are amplified through energy and angular momentum transfer, can also occur in horizonless objects like boson stars, due to both the real space and internal field space rotations. In this work, we study superradiance in the frequency and time domains for static and spinning boson stars, constructed within general relativity and with a self-interacting complex scalar field as a matter source. Using linear perturbation theory and three dimensional nonlinear simulations, we calculate amplification factors and analyze energy and angular momentum transfer in scattering processes, with results showing consistency between approaches. Wave scattering inside a cavity containing a boson star is also examined, demonstrating the effects of confinement on amplification.

A large population of binary systems in the Universe emitting gravitational waves (GW) would produce a stochastic noise, known as the gravitational wave background (GWB). The properties of the GWB directly depend on the attributes of its constituents. If the binary systems are in eccentric orbits, it is well established that the GW power they radiate strongly depends on their instantaneous orbital phase. Consequently, their power spectrum varies over time, and the resulting GWB can appear nonstationary. In this work, we estimate the amplitude of time-dependent fluctuations in the GWB power spectrum as a function of the eccentricity of the binaries. Specifically, we focus on the GWB produced by a population of supermassive black hole binaries (SMBHB) that should be observable by pulsar timing arrays (PTA). We show that a large population of homogeneously distributed equal SMBHBs produces nonstationary features that are undetectable by current PTA datasets. However, using more realistic and astrophysically motivated populations of SMBHBs, we show that the nonstationarity might become very large and detectable, especially in the case of more massive and eccentric populations. In particular, when one binary is slightly brighter than the GWB, we demonstrate that time fluctuations can become significant. This is also true for individual binary systems with a low signal-to-noise ratio (SNR) relative to the GWB (SNR $\approx$ 1), which standard data analysis methods would struggle to detect. The detection of nonstationary features in the GWB could indicate the presence of some relatively bright GW sources in eccentric orbits, offering new insights into the origins of the signal.

Galactic dark matter (DM) particles, having non-gravitational interactions with nucleons, can interact with stellar constituents and eventually become captured within stars. Over the lifetime of the celestial body, these non-annihilating, heavy DM particles may accumulate and eventually form a comparable stellar mass black hole (BH), referred to as a Transmuted Black Hole (TBH). We investigate how current gravitational wave (GW) experiments could detect such particle DM through the presence of low-mass TBHs, which cannot form via standard stellar evolution. Different stellar objects (compact and non-compact) provide laboratories across DM-nucleon interaction regimes, offering insight into DM's mass and its non-gravitational properties.

Axion-like particles (ALPs) can decay into two photons with a rest-frame frequency given by half of the ALP mass. This implies that ultra-violet searches can be used to investigate ALPs in the multi-eV mass range. We use archival data from the Hubble Space Telescope between 110 and 170 nm to constrain ALPs with mass between 14.4-22.2 eV. We consider observations of a set of dwarf spheroidal galaxies and galaxy clusters and assume the ALP density in these objects to follow their dark matter density. The derived limit on the ALP-photon coupling $g_{a\gamma}$ excludes values above $10^{-12}~{\rm GeV}^{-1}$ over the whole mass range and surpasses previous limits by over one order of magnitude.

The singlet-doublet vector-like fermion dark matter model has been extensively studied in the literature over the past decade. An important parameter in this model is the singlet-doublet mixing angle ($\sin\theta$). All the previous studies have primarily focused on annihilation and co-annihilation processes for obtaining the correct dark matter relic density, assuming that the singlet and doublet components decouple at the same epoch. In this work, we demonstrate that this assumption holds only for larger mixing angles with a dependency on the mass of the dark matter. However, it badly fails for the mixing angle $\sin\theta<0.05$. We present a systematic study of the parameter space of the singlet-doublet dark matter relic, incorporating annihilation, co-annihilation, and, for the first time, co-scattering processes. Additionally, the freeze-in parameter space is also explored. We found that due to the inclusion of co-scattering processes, the correct relic density parameter space is shifted towards the detection sensitivity range of the LHC and MATHUSLA via displaced vertex signatures.

We study supercooled first-order phase transitions above the QCD scale in a wide class of conformal Majoron-like U(1)' models that explain the totality of active neutrino oscillation data and produce a detectable stochastic gravitational wave background (SGWB) at LIGO, LISA and ET. We place constraints on the U(1)' breaking scale and gauge coupling using current LIGO-Virgo-Kagra data. We find that strong supercooling can be ruled out in large regions of parameter space if a SGWB is not detected by these experiments. A null signal at LIGO and ET will disfavor a type-I seesaw scale above $10^{14}$ GeV, while a positive signal is a signature of heavy right-handed neutrinos. On the other hand, LISA will be sensitive to seesaw scales as low as a TeV, and could detect a SGWB even if the right-handed neutrinos are decoupled.

Bayesian inference with stochastic sampling has been widely used to obtain the properties of gravitational wave (GW) sources. Although computationally intensive, its cost remains manageable for current second-generation GW detectors because of the relatively low event rate and signal-to-noise ratio (SNR). The third-generation (3G) GW detectors are expected to detect hundreds of thousands of compact binary coalescence events every year with substantially higher SNR and longer signal duration, presenting significant computational challenges. In this study, we systematically evaluate the computational costs of source parameter estimation (PE) in the 3G era by modeling the PE time cost as a function of SNR and signal duration. We examine the standard PE method alongside acceleration methods including relative binning, multibanding, and reduced order quadrature. We predict that PE for a one-month-observation catalog with 3G detectors could require billions to quadrillions of CPU core hours with the standard PE method, whereas acceleration techniques can reduce this demand to millions of core hours. These findings highlight the necessity for more efficient PE methods to enable cost-effective and environmentally sustainable data analysis for 3G detectors. In addition, we assess the accuracy of accelerated PE methods, emphasizing the need for careful treatment in high-SNR scenarios.