Papers for Monday, Mar 11 2024

JWST observations of galaxies at $z\gtrsim 8$ suggest that they are more luminous and clumpier than predicted by most models, prompting several proposals on the physics of star formation and feedback in the first galaxies. In this paper, we focus on the role of ultraviolet (UV) radiation in regulating star formation by performing a set of cosmological radiation hydrodynamics simulations of one galaxy at sub-pc resolution with different radiative feedback models. We find that the suppression of cooling by far UV (FUV) radiation (i.e., $\mathrm{H_2}$ dissociating radiation) from Pop II stars is the main physical process triggering the formation of compact and massive star clusters and is responsible for the bursty star formation observed in metal-poor galaxies at $z\gtrsim 10$. Indeed, artificially suppressing FUV radiation leads to a less intense continuous mode of star formation distributed into numerous, but low-mass open star clusters. Due to the intense FUV field, low-metallicity clouds remain warm ($\sim 10^4\,\mathrm{K}$) until they reach a relatively high density ($\gtrsim 10^3\,\mathrm{cm^{-3}}$), before becoming self-shielded and transitioning to a colder ($\sim 100\,\mathrm{K}$), partially molecular phase. As a result, star formation is delayed until the clouds accumulate enough mass to become gravitationally unstable. At this point, the clouds undergo rapid star formation converting gas into stars with high efficiency. We, therefore, observe exceptionally bright galaxies (ten times brighter than for continuous star formation) and subsequent quenched "dead" galaxies that did not form stars for tens of Myrs.

JWST observations are leading to important new insights into exoplanetary atmospheres through transmission spectroscopy. In order to harness the full potential of the broad spectral range and high sensitivity of JWST, atmospheric retrievals of exoplanets require a high level of robustness and accuracy in the underlying models. We present the VIRA retrieval framework which implements a range of modelling and inference capabilities motivated by early JWST observations of exoplanet transmission spectra. This includes three complementary approaches to modelling atmospheric composition, three atmospheric aerosol models, including a physically-motivated Mie scattering approach, and consideration of correlated noise. VIRA enables a cascading retrieval architecture involving a sequence of retrievals with increasing sophistication. We demonstrate VIRA using a JWST transmission spectrum of the hot Saturn WASP-39 b in the $\sim$1-5 $\mu$m range. In addition to confirming prior chemical inferences, we retrieve molecular abundances for H$_2$O, CO, CO$_2$, SO$_2$ and H$_2$S, resulting in super-solar elemental abundances of log(O/H)=$-2.0\pm0.2$, log(C/H)=$-2.1\pm0.2$ and log(S/H)=$-3.6\pm0.2$, along with C/O and S/O ratios of $0.83^{+0.05}_{-0.07}$ and $0.029^{+0.012}_{-0.009}$, respectively, in the free chemistry case. The abundances correspond to $20.1^{+10.5}_{-8.1}\times$, $28.2^{+16.3}_{-12.1}\times$ and $20.8^{+10.3}_{-7.5}\times$ solar values for O/H, C/H and S/H, respectively, compared to C/H $= 8.67\pm0.35 \times$ solar for Saturn. Our results demonstrate how JWST transmission spectroscopy combined with retrieval frameworks like VIRA can measure multi-elemental abundances for giant exoplanets and enable comparative characterisation with solar system planets.

Many gravitational wave (GW) sources are expected to have non-negligible eccentricity in the millihertz band. These highly eccentric compact object binaries may commonly serve as a progenitor stage of GW mergers, particularly in dynamical channels where environmental perturbations bring a binary with large initial orbital separation into a close pericenter passage, leading to efficient GW emission and a final merger. This work examines the stochastic GW background from highly eccentric ($e\gtrsim 0.9$), stellar-mass sources in the mHz band. Our findings suggest that these binaries can contribute a substantial GW power spectrum, potentially exceeding the LISA instrumental noise at $\sim 3-7$~mHz. This stochastic background is likely to be dominated by eccentric sources within the Milky Way, thus introducing anisotropy and time dependence in LISA's detection. However, given efficient search strategies to identify GW transients from highly eccentric binaries, the unresolvable noise level can be substantially lower, approaching $\sim 2$ orders of magnitude below the LISA noise curve. Therefore, we highlight the importance of characterizing stellar-mass GW sources with extreme eccentricity, especially their transient GW signals in the millihertz band.

The first generations of stars left their chemical fingerprints on metal-poor stars in the Milky Way and its surrounding dwarf galaxies. While instantaneous and homogeneous enrichment implies that groups of co-natal stars should have the same element abundances, small amplitudes of abundance scatter are seen at fixed [Fe/H]. Measurements of intrinsic abundance scatter have been made with small, high-resolution spectroscopic datasets where measurement uncertainty is small compared to this scatter. In this work, we present a method to use mid-resolution survey data, which has larger errors, to make this measurement. Using APOGEE DR17, we calculate the intrinsic scatter of Al, O, Mg, Si, Ti, Ni, and Mn relative to Fe for 333 metal-poor stars across 6 classical dwarf galaxies around the Milky Way, and 1604 stars across 19 globular clusters. We first calibrate the reported abundance errors in bins of signal-to-noise and [Fe/H] using a high-fidelity halo dataset. We then apply these calibrated errors to the APOGEE data, and find small amplitudes of average intrinsic abundance scatter in dwarf galaxies ranging from 0.032 $-$ 0.14 dex with a median value of 0.043 dex. For the globular clusters, we find intrinsic scatters ranging from 0.018 $-$ 0.21 dex, with particularly high scatter for Al and O. Our measurements of intrinsic abundance scatter place important upper limits on the intrinsic scatter in these systems, as well as constraints on their underlying star formation history and mixing, that we can look to simulations to interpret.

S$H_0$ES $VI$-band photometry for classical Cepheids in the keystone galaxy NGC4258 yield discrepant absolute magnitudes. Specifically, the 2016 and 2022 published S$H_0$ES Cepheid data for NGC4258 exhibit a substantial offset of $\Delta W_{0,VI}\simeq0^{\rm m}.3$. That adds to a suite of existing concerns associated with the S$H_0$ES analysis, which in sum imply that their relatively non-changing Hubble constant for nearly twenty years should be viewed with caution.

During the epoch of reionization (EoR), the 21-cm signal allows direct observation of the neutral hydrogen (\hi{}) in the intergalactic medium (IGM). In the post-reionization era, this signal instead probes \hi{} in galaxies, which traces the dark matter density distribution. With new numerical simulations, we investigated the end stages of reionization to elucidate the transition of our Universe into the post-reionization era. Our models are consistent with the latest high-redshift measurements, including ultraviolet (UV) luminosity functions. Notably, these models consistently reproduced the evolution of the UV photon background, which is constrained from Lyman-$\alpha$ absorption spectra. We studied the dependence of this background on the nature of photon sinks in the IGM, requiring mean free path of UV photons to be $\sim$10 comoving-megaparsecs (cMpc) during the EoR that increases gradually with time during late stages ($z\lesssim 6$). Our models revealed that the reionization of the IGM transitioned from an \textit{inside-out} to an \textit{outside-in} process when the Universe is less than 0.01 per cent neutral. During this epoch, the 21-cm signal also shifted from probing predominantly the \hi{} in the IGM to that in galaxies. Furthermore, we identified a statistically significant number of large neutral islands (with sizes up to 40 cMpc) persisting until very late stages ($5 \lesssim z \lesssim 6$) that can imprint features in Lyman-$\alpha$ absorption spectra and also produce a knee-like feature in the 21-cm power spectrum.

We study the fraction of the intra-cluster light (ICL) formed in-situ in the three most massive clusters of the TNG50 simulation, with virial masses $\sim 10^{14}$ M$_{\odot}$. We find that a significant fraction of ICL stars ($8\%$-$28\%$) are born in-situ. This amounts to a total stellar mass comparable to the central galaxy itself. Contrary to simple expectations, only a sub-dominant fraction of these in-situ ICL stars are born in the central regions and later re-distributed to more energetic orbits during mergers. Instead, many in-situ ICL stars form directly hundreds of kiloparsecs away from the central galaxy, in clouds condensing out of the circum-cluster medium. The simulations predict a present-date diffuse star formation rate of $\sim$1 $\mathrm{M}_{\odot}$/yr, with higher rates at higher redshifts. The diffuse star forming component of the ICL is filamentary in nature, extends for hundreds of kiloparsecs and traces the distribution of neutral gas in the cluster host halo. We discuss briefly how numerical details of the baryonic treatment in the simulation may play a role in this result and conclude that a sensitivity of $1.6 \times 10^{-19} - 2.6 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ in H$_\alpha$ flux (beyond current observational capabilities) would be necessary to detect this diffuse star-forming component in galaxy clusters.

Planetary-mass objects and brown dwarfs at the transition from relatively red L dwarfs to bluer mid-T dwarfs $\left(\rm{T}_{eff}\sim1300 \ \text{K} \right)$ show enhanced spectrophotometric variability. An open question is whether this variability is caused by atmospheric planetary-scale (Kelvin or Rossby) waves or by large spots associated with the precipitation of silicate and metal clouds. We applied both waves and spotted models to fit near-infrared (NIR), multi-band ($Y$/$J$/$H$/$K$) photometry of SIMP J013656.5+093347 (hereafter SIMP0136), collected at the Canada-France-Hawaii Telescope using the Wide-field InfraRed Camera. SIMP0136 is a planetary-mass object (12.7$\pm1.0 \ \rm{M_J}$) at the L/T transition (T2$\pm0.5$) known to exhibit light curve evolution over multiple rotational periods. We measure the maximum peak-to-peak variability of $6.17\pm0.46\%$, $6.45\pm0.33\%$, $6.51\pm0.42\%$, and $4.33\pm0.38\%$ in the $Y$, $J$, $H$, and $K$ bands respectively, and find evidence that wave models are preferred for all four NIR bands. Furthermore, we determine the spot size necessary to reproduce the observed variations is larger than the Rossby deformation radius and Rhines scale, which is unphysical. Through the correlation between light curves produced by the waves and associated color variability, we find evidence of planetary-scale, wave-induced cloud modulation and breakup, similar to Jupiter's atmosphere and supported by general circulation models. We also detect a $93.8^{\circ}\pm7.4^{\circ}$ ($12.7\sigma$) phase shift between the $H-K$ and $J-H$ color time series, providing evidence for complex vertical cloud structure in SIMP0136's atmosphere.

The prediction of a nearly scale-invariant spectrum of curvature and tensor fluctuations is among the main features of cosmic inflation. The current measurements of the primordial fluctuations in the cosmic microwave background (CMB) provide tight constraints on the amplitude of the scalar and tensor spectra, and the scalar tilt. However, the precise connection between these observables and a given inflationary model, depends on the expansion history between the end of inflation and the beginning of the radiation dominated era, which corresponds to the reheating epoch. This mapping between horizon exit and reentry of fluctuations, parametrized by the number of $e$-folds $N_*$, can therefore be affected by the presence of a transient epoch of non-perturbative particle production during reheating ({\em preheating}). Using a combination of perturbative and lattice computations, we quantify the impact of preheating in a non-equilibrated dark matter sector on the CMB observables, under the assumption of a simultaneous perturbative decay of the inflaton into Standard Model particles. Combined with structure formation constraints, this allows us to impose stringent bounds on the post-inflationary reheating temperature.

The large-scale mass distributions of strong-lensing galaxies have long been assumed to be well-described by a singular ellipsoidal power-law density profile with external shear. However, the inflexibility of this model could lead to systematic errors in astrophysical parameters inferred with gravitational lensing observables. Here, we present observations with the Atacama Large (sub-)Millimetre Array (ALMA) of three strongly lensed dusty star-forming galaxies at $\simeq30$ mas angular resolution and investigate the sensitivity of these data to angular structure in the lensing galaxies. We jointly infer the lensing mass distribution and the full surface brightness of the lensed sources with multipole expansions of the power-law density profile up to fourth order using a technique developed for interferometric data. All three data sets strongly favour third and fourth-order multipole amplitudes of $\approx1$ percent of the convergence. While the infrared stellar isophotes and isodensity shapes agree for one lens system, for the other two the isophotes disagree to varying extents, suggesting contributions to the angular structure from dark matter intrinsic or extrinsic to the lensing galaxy.

Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths difficult to access from the ground. Aims. We aim to image the circumstellar environment of PDS 70 at 1.87 $\mu$m and 4.83 $\mu$m, assess the presence of Pa-$\alpha$ emission due to accretion onto the protoplanets, and probe any IR excess indicative of heated circumplanetary material. Methods. We obtain non-coronagraphic JWST/NIRCam images of PDS 70 within the MINDS (MIRI mid-INfrared Disk Survey) program. We leverage the Vortex Image Processing (VIP) package for data reduction, and develop dedicated routines for optimal stellar PSF subtraction, unbiased imaging of the disk, and protoplanet flux measurement in this type of dataset. A radiative transfer model of the disk is used to disentangle the contributions from the disk and the protoplanets. Results. We re-detect both protoplanets and identify extended emission after subtracting a disk model, including a large-scale spiral-like feature. We interpret its signal in the direct vicinity of planet c as tracing the accretion stream feeding its circumplanetary disk, while the outer part of the feature may rather reflect asymmetric illumination of the outer disk. We also report a bright signal consistent with a previously proposed protoplanet candidate enshrouded in dust, near the 1:2:4 mean-motion resonance with planets b and c. The 1.87 $\mu$m flux of planet b is consistent with atmospheric model predictions, but not that of planet c. We discuss potential origins for this discrepancy, including significant Pa-$\alpha$ line emission. The 4.83 $\mu$m fluxes of planets b and c suggest enshrouding dust or heated CO emission from their circumplanetary environment.

Supermassive black holes have been known to disrupt passing stars, producing outbursts called tidal disruption events (TDEs). TDEs have recently gained attention due to their unique dynamics and emission processes, which are still not fully understood. Especially, the so-called optical TDEs, are of interest as they often exhibit delayed or obscured X-ray emission from the accretion disk, making the origin of the prompt emission unclear. In this paper, we present the optical photo-, spectro-, and polarimetry of a recent TDE candidate AT 2022fpx, in addition to accompanying monitoring observations in ultraviolet and X-rays. The optical spectra of AT 2022fpx show Bowen fluorescence as well as highly-ionized iron emission lines, which are characteristic of extreme coronal line emitters. Additionally, the source exhibits variable but low-polarized emission at the outburst peak, with a clear rotation of the polarization angle. X-ray emission observed later in the decay appear flare-like but is consistent with constant temperature black-body emission. The overall outburst decay is slower than for typical TDEs, and resembles more the ones seen from Bowen fluorescence flares. These observations suggest that AT 2022fpx could be a key source in linking different long-lived TDE scenarios. Its unique characteristics, such as the extreme coronal line emission, variable polarization, and the delayed X-ray flare, can be attributed to the outer shock scenario or a clumpy torus surrounding the supermassive black hole. Further studies, especially in the context of multi-wavelength observations, are crucial to fully understand the dynamics and emission mechanisms of these intriguing astrophysical events.

The most dynamically relaxed clusters of galaxies play a special role in cosmological studies as well as astrophysical studies of the intracluster medium (ICM) and active galactic nucleus feedback. While high spatial resolution imaging of the morphology of the ICM has long been the gold standard for establishing a cluster's dynamical state, such data are not available from current or planned surveys, and thus require separate, pointed follow-up observations. With optical and/or near-IR photometric imaging, and red-sequence cluster finding results from those data, expected to be ubiquitously available for clusters discovered in upcoming optical and mm-wavelength surveys, it is worth asking how effectively photometric data alone can identify relaxed cluster candidates, before investing in, e.g., high-resolution X-ray observations. Here we assess the ability of several simple photometric measurements, based on the redMaPPer cluster finder run on Sloan Digital Sky Survey data, to reproduce X-ray classifications of dynamical state for an X-ray selected sample of massive clusters. We find that two simple metrics contrasting the Bright Central Galaxy (BCG) to other cluster members can identify a complete sample of relaxed clusters with a purity of ~40 per cent in our data set. Including minimal ICM information in the form of a center position increases the purity to ~60 per cent. However, all three metrics depend critically on correctly identifying the BCG, which is presently a challenge for optical red-sequence cluster finders.

We present the results of a search for Extended Emission-Line Regions (EELRs) ionized by extant or recently-faded active galactic nuclei (AGN), using [O III] narrowband imaging and spectroscopic followup. The sample includes 198 galaxies in 92 strongly interacting or merging galaxy systems in the range z=0.009-0.0285. Among these, three have EELRs extended beyond 10 kpc in projection from the nucleus detected in previous studies. We identify a single new distant emission region, projected 35 kpc from UGC 5941. Our optical spectrum does not detect He II, but its strong-line ratios put this in the same class as securely characterized EELR clouds. The nucleus of UGC 5941 is dominated by recent star formation, preventing detection of any weak ongoing AGN. Overall counts of distant EELRs in this and the previous TELPERION samples give incidence 2-5% depending on galaxy and AGN selection, 20-50 times higher than the Galaxy Zoo EELR survey with its higher surface-brightness threshold and much larger input sample. AGN in interacting and merging systems have an increased detection rate $12+/-6$%, while none are detected around noninteracting AGN. Some of these AGN are at luminosity low enough to require additional X-ray or far-infrared information to tell whether the EELR ionization level suggests long-term fading.

We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which have passed human inspection and morphological classification, and a larger sample of $\sim$100,000 classified by neural network models. Here, we report on the observed properties of these samples, and demonstrate that tremendous insight can be gained from just the observed properties of clusters, even in the absence of their transformation into physical quantities. In particular, we show the utility of the UBVI color-color diagram, and the three principal features revealed by the PHANGS-HST cluster sample: the young cluster locus, the middle-age plume, and the old globular cluster clump. We present an atlas of maps of the 2D spatial distribution of clusters and compact associations in the context of the molecular clouds from PHANGS-ALMA. We explore new ways of understanding this large dataset in a multi-scale context by bringing together once-separate techniques for the characterization of clusters (color-color diagrams and spatial distributions) and their parent galaxies (galaxy morphology and location relative to the galaxy main sequence). A companion paper presents the physical properties: ages, masses, and dust reddenings derived using improved spectral energy distribution (SED) fitting techniques.

Following the Pluto fly-by of the New Horizons spacecraft, the mission provided a unique opportunity to explore the Kuiper Belt in-situ. The possibility existed to fly-by a Kuiper Belt object (KBO) as well as to observe additional objects at distances closer than are feasible from earth-orbit facilities. However, at the time of launch no KBOs were known about that were accessible by the spacecraft. In this paper we present the results of 10 years of observations and three uniquely dedicated efforts -- two ground-based using the Subaru Suprime Camera, the Magellan MegaCam and IMACS Cameras, and one with the Hubble Space Telescope -- to find such KBOs for study. In this paper we overview the search criteria and strategies employed in our work and detail the analysis efforts to locate and track faint objects in the galactic plane. We also present a summary of all of the KBOs that were discovered as part of our efforts and how spacecraft targetability was assessed, including a detailed description of our astrometric analysis which included development of an extensive secondary calibration network. Overall, these efforts resulted in the discovery of 89 KBOs including 11 which became objects for distant observation by New Horizons and (486958) Arrokoth which became the first post-Pluto fly-by destination.

We test Deep-Learning (DL) techniques for the analysis of rotational core-collapse supernovae (CCSN) gravitational-wave (GW) signals by performing classification and parameter inference of the GW strain amplitude ($D \cdot \Delta h $) and the maximum (peak) frequency $(f_\textrm{peak})$, attained at core bounce. Our datasets are built from a catalog of numerically generated CCSN waveforms assembled by Richers 2017. Those waveforms are injected into noise from the Advanced LIGO and Advanced Virgo detectors corresponding to the O2 and O3a observing runs. For a signal-to-noise ratio (SNR) above 5, our classification network using time series detects Galactic CCSN GW signals buried in detector noise with a false positive rate (FPR) of 0.10% and a 98% accuracy, being able to detect all signals with SNR>10. The inference of $f_\textrm{peak}$ is more accurate than for $D \cdot \Delta h $, particularly for our datasets with the shortest time window (0.25 s) and for a minimum SNR=15. From the calibration plots of predicted versus true values of the two parameters, the standard deviation ($\sigma$) and the slope deviation with respect to the ideal value are computed. We find $\sigma_{D \cdot \Delta h}$ = 52.6cm and $\sigma_{f_\textrm{peak}}$ = 18.3Hz, with respective slope deviations of 11.6% and 8.3%. Our best model is also tested on waveforms from a recent CCSN catalog built by Mitra 2023, different from the one used for the training. For these new waveforms the true values of the two parameters are mostly within the 1$\sigma$ band around the network's predicted values. Our results show that DL techniques hold promise to infer physical parameters of Galactic rotational CCSN even in the presence of real (non-Gaussian) noise conditions from current GW detectors.

We examine the simple model put forth in a recent note by Loeb 2024 regarding the brightness of space debris in the size range of 1-10 cm and their impact on the Rubin Observatory LSST transient object searches. Their main conclusion was that "image contamination by untracked space debris might pose a bigger challenge [than large commercial satellite constellations in LEO]". Following corrections and improvements to this model, we calculate the apparent brightness of tumbling Low Earth Orbit (LEO) debris of various sizes, and we briefly discuss the likely impact and potential mitigations of glints from space debris in LSST. The largest difference from the Loeb 2024 estimates is that 1-10 cm debris in LEO pose no threat to LSST transient object alert generation because their signal-to-noise ratio (SNR) for detection will be much lower than estimated by Loeb 2024. Most of the difference in predicted SNR, about a factor of six, arises from defocus of LEO objects due to the large Simonyi Survey Telescope primary mirror. We find that only tumbling LEO debris larger than 10 cm or with significantly greater reflectivity, which give 1 ms glints, might be detected with high confidence (SNR>5). We estimate that only one in five LSST exposures low on the sky during twilight might be affected. More slowly tumbling objects of larger size can give flares in brightness that are easily detected, however these will not be cataloged by the LSST Science Pipelines Bosch et al. 2019 because of the resulting long streak.

Extreme-mass ratio inspirals (EMRIs) could be detected by space-borne gravitational-wave (GW) detectors, such as the Laser Interferometer Space Antenna (LISA), TianQin and Taiji. In general, locating EMRIs by space-borne GW detectors can help us select the candidate host galaxies which can be used to infer the cosmic expansion history. In this paper, we show that the localization information of EMRIs can also be used to select the candidate host AGNs which can be used to infer the formation channel of the EMRIs and extract more precisely the redshift probability distributions of the EMRIs. Using the EMRIs that are expected to be detected by TianQin and LISA and the galaxy catalog that can be provided by the Chinese Space Station Telescope for the analysis, we find that TianQin can constrain the Hubble-Lema\^itre constant $H_0$ to a precision of about $3\%-8\%$ and the dark energy equation of state parameter $w_0$ to about $10\%-40\%$, the network composed of TianQin and LISA can improve the precisions of $H_0$ and $w_0$ to about $0.4\% - 7\%$ and $4\%-20\%$, respectively, without considering the effects of AGNs. Furthermore, combining the detected EMRIs with AGN catalog by the statistical framework of likelihood-ratio-based, we find TianQin can establish the EMRI-AGN correlation with about $500$ EMRIs if all of the EMRIs indeed originate in AGNs, and the TianQin+LISA network can reduce this required number to about $30$. Once the EMRI-AGN correlation is confirmed, the combination of the EMRIs and AGN catalog can significantly improve the constraints on the cosmological parameters. These results demonstrate the potentials of using EMRIs as well as galaxy and AGN catalogs to constrain the cosmological parameters and astrophysical models.

We evaluate recent and upcoming low-redshift neutral hydrogen (HI) surveys as a cosmological probe of small scale structure with a goal of determining the survey criteria necessary to test ultra-light axion (ULA) dark matter models. Standard cold dark matter (CDM) models predict a large population of low-mass galactic halos, whereas ULA models demonstrate significant suppression in this small-scale regime, with halo mass cutoffs of $10^{12}\, \mathrm{M}_{\odot}$ to $10^{7}\, \mathrm{M}_{\odot}$ corresponding to ULA masses of $10^{-24}\,$eV to $10^{-20}\,$eV, respectively. We generate random, homogeneously populated mock universes with cosmological parameters adjusted to match CDM and ULA models. We simulate observations of these mock universes with hypothetical analogs of the mass-limited ALFALFA and WALLABY HI surveys and reconstruct the corresponding HI mass function (HIMF). We find that the ALFALFA HIMF can test for the presence of ULA DM with $m_{a}\lesssim 10^{-21.5}\,$eV, while WALLABY could reach the larger window $m_{a}\lesssim 10^{-20.9}\,$eV. These constraints are complementary to other probes of ULA dark matter, demonstrating the utility of local Universe HI surveys in testing dark matter models.

We report on the properties of pulsed X-ray emission from eight MeV pulsars using XMM-Newton, NICER, NuSTAR and HXMT data. For the five among eight MeV pulsars, the X-ray spectra can be fitted by a broken-power law model with a break energy of $\sim5-10$ keV. The photon index below and above break energy are $\sim 1$ and $\sim 1.5$, respectively. In comparison with the X-ray emission of the $Fermi$-LAT pulsars, the MeV pulsars have a harder spectrum and ahigher radiation efficiency in 0.3-10 keV energy bands. By assuming the isotropic emission, the emission efficiency in the keV-MeV bands is estimated to be $\eta_{MeV}\sim 0.01-0.1$, and it is similar to the efficiency of GeV emission of the $Fermi$-LAT pulsars that have similar spin-down power. To explain the observed efficiency of the MeV pulsars, we estimate the required pair multiplicity as $10^{4-7}$ that depends on the emission process (curvature radiation or synchrotron radiation) and the location in the magnetosphere. The large multiplicity indicates that the secondary pairs that are created by a pair-creation process of the GeV photons produce the X-ray/soft gamma-ray emissions of the MeV pulsars. We speculate that the difference between the MeV pulsars and $Fermi$-LAT pulsars is attributed to the difference in viewing angle measured from the spin-axis, if the emission originates from a region inside the light cylinder (canonical gap model) or the difference in the inclination angle of the magnetic axis, if the emission is produced from equatorial current sheet outside the light cylinder.

Leveraging a high-resolution 3D dust map of the solar neighborhood from Edenhofer et al. (2023), we derive a new 3D model for the dust-traced surface of the Local Bubble, the supernova-driven cavity surrounding the Sun. We find that the surface of the Local Bubble is highly irregular in shape, with its peak extinction surface falling at an average distance of 170 pc from the Sun (spanning 70-600+ pc) with a typical thickness of 35 pc and a total dust-traced mass of $(6.0 \pm 0.7) \times 10^5 \ \rm{M}_{\odot}$. The Local Bubble displays an extension in the Galactic Northern hemisphere that is morphologically consistent with representing a "Local Chimney." We argue this chimney was likely created by the "bursting" of this supernova-driven superbubble, leading to the funneling of interstellar medium ejecta into the lower Galactic halo. We find that many well-known dust features and molecular clouds fall on the surface of the Local Bubble and that several tunnels to other adjacent cavities in the interstellar medium may be present. Our new, parsec-resolution view of the Local Bubble may be used to inform future analysis of the evolution of nearby gas and young stars, the investigation of direct links between the solar neighborhood and the Milky Way's lower halo, and numerous other applications.

We revisit the proposal that an energy transfer from dark energy into dark matter can be described in field theory by a first order phase transition. We analyze the model proposed in Ref. Abdalla et al. (2013), using updated constraints on the decay time of a metastable dark energy from the work of Ref. Shafieloo et al. (2018). The results of our analysis show no prospects for potentially observable signals that could distinguish this scenario. We also show that such model would not drive a complete transition to a dark matter dominated phase even in a distant future. Nevertheless, the model is not excluded by the latest data and we confirm that the mass of the dark matter particle that would result from such a process corresponds to an axion-like particle, which is currently one of the best motivated dark matter candidates. We argue that extensions to this model, possibly with additional couplings, still deserve further attention as it could provide an interesting and viable description for an interacting dark sector scenario based in a single scalar field.

Neutron star$-$white dwarf (NSWD) binaries are one of the most abundant sources of gravitational waves (GW) in the Milky Way. These GW sources are the evolutionary products of primordial binaries that experienced many processes of binary interaction. We employ a binary population synthesis method to investigate the properties of Galactic NSWD binaries detectable by the Laser Interferometer Space Antenna (LISA). In this paper, only the NSWD systems with a COWD or ONeWD component are included. We consider various models related to mass transfer efficiencies during primordial binary evolution, supernova explosion mechanisms at NS formation, common envelope ejection efficiencies, and critical WD masses that determining the stability of mass transfer between WDs and NSs. Based on our calculations, we estimate that tens to hundreds of LISA NSWD binaries exist in the Milky Way. We find that the detection of LISA NSWD binaries is able to provide profound insights into mass transfer efficiencies during the evolution of primordial binaries and critical WD masses during mass transfer from a WD to an NS.

Primordial black holes (PBHs) are the plausible candidates for the cosmological dark matter. Theoretically, PBHs with masses $M_{\rm PBH}$ in the range of $4\times10^{14}\sim 10^{17}\,{\rm g}$ can emit sub-GeV electrons and positrons through Hawking radiation. Some of these particles could undergo diffusive reacceleration during their propagation in the Milky Way, resulting in acceleration up to the GeV level. This falls within the observation energy range of AMS-02. In this work, we utilize AMS-02 data to constrain the PBH abundance $f_{\rm PBH}$ employing the reacceleration mechanism. Our findings reveal that the limit is more stringent than those derived from Voyager 1 for nearly $\mathcal{O}(1)$ at $M_{\rm PBH}\sim10^{16}\,{\rm g}$ under the assumption of a monochromatic PBH mass distribution. The constraints are even more robust in a more realistic scenario involving a log-normal mass distribution of PBHs. Moreover, we explore the impact of varying propagation parameters and solar modulation potential within reasonable ranges, and find that such variations have minimal effects on the final results.

We investigate the morphological and structural evolutions of disk galaxies in simulations for a wide range of Toomre's $Q$ parameter. In addition to the inspection of conventional bar modes, we compute the concentration, asymmetry and clumpiness (CAS) parameters to enlarge the understanding of the galaxy evolution. These parameters are widely employed to analyze the light distribution of the observed galaxies, but the adaptation to numerical simulations is not much considered. While the bar formation takes place in a considerable range of $Q$ around $1$, barred galaxies originating from $Q>1$ and $Q<1$ disks yield the CAS values that differ significantly. Disks starting with $Q<1$ develop clumps due to local gravitational instabilities along with the bar and these clumps play a central role in enhancing the CAS values. That process is absent in $Q>1$ counterparts in which the evolution is dominated by linearly unstable two-armed modes that lead to lower CAS values. Likewise, unbarred galaxies that are obtainable from disks with $Q$ far below and far above $1$ exhibit greatly different CAS magnitudes. It turns out that the CAS parameters can serve as indicators of the initial kinematical state and the evolution history of a disk of any morphology. In addition, we find an alternative mechanism of the formation of the lopsided barred galaxy when $Q\lesssim 1$. Bars that evolve in the midst of the clumps can spontaneously become lopsided at the end.

In protoplanetary disks, atomic carbon is expected to originate from the PDR at the disk surface where CO is dissociated by UV photons coming from the stellar, or external interstellar, radiation field. Even though atomic carbon has been detected in several protoplanetary disks, there is a lack of spatially resolved observations of it. For HD 163296 protoplanetary disk, we aim to obtain both radial and vertical structure of [CI]$={^3}{P_1} - {^3}{P_0}$ line emission and perform the first direct comparison of this tracer with optically thick line emission $^{12}$CO $J=2-1$. We used archival ALMA data for [CI]$={^3}{P_1} - {^3}{P_0}$ and previously published ${^{12}}$ CO $J=2-1$ data in HD 163296. Through the software of disksurf we extracted the vertical structure, meanwhile radial profiles were obtained directly from imaging. Brand new DALI modelling was employed to perform direct comparison to the data. We found that these tracers are collocated radially but not vertically, where $^{12}$CO $J=2-1$ emission is, on average, located at higher altitudes, as it is also the case for other tracers in the same disk. Due to this difference in vertical height of the emission, the optically thick $^{12}$CO $J=2-1$ emission line appears to trace the highest altitudes, despite the expected formation mechanism of [CI] in the disk. The latter phenomena may be due to efficient mixing of the upper layers of the disk, or UV photons penetrating deeper than we expected.

The unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 micron are ubiquitously seen in a wide variety of astrophysical regions and commonly attributed to polycyclic aromatic hydrocarbon (PAH) molecules. However, the unambiguous identification of any individual, specific PAH molecules has proven elusive until very recently two isomers of cyanonapthalene, which consists of two fused benzene rings and substitutes a nitrile (-CN) group for a hydrogen atom, were discovered in the Taurus Molecular Cloud based on their rotational transitions at radio frequencies. To facilitate the James Webb Space Telescope (JWST) to search for cyanonapthalenes in astrophysical regions, we model the vibrational excitation of cyanonapthalenes and calculate their infrared emission spectra in a number of representative astrophysical regions. The model emission spectra and intensities will allow JWST to quantitatively determine or place an upper limit on the abundance of cyanonapthalenes.

A large fraction of Fermi-LAT sources in the 4th Fermi LAT 14-year Catalog (4FGL) still remain unidentified (unIDed). We continued to improve our machine learning pipeline (MUWCLASS) and used it to classify 1206 X-ray sources located within the extent of 73 unIDed 4FGL sources with Chandra X-ray Observatory (CXO) observations included in the Chandra Source Catalog 2.0. Recent improvements to our pipeline include astrometric corrections, probabilistic cross-matching to lower-frequency counterparts, and a more realistic oversampling method. X-ray sources are classified into 8 broad pre-determined astrophysical classes defined in the updated training dataset which we also release. The classifications give 115 plausible X-ray counterparts to 45 GeV sources. We present details of the machine learning classification, describe the pipeline improvements, and perform additional spectral and variability analysis for brighter sources. We identify 3 GeV sources as isolated neutron star candidates, 1 as an X-ray binary candidate, 16 as active galactic nucleus candidates, 7 sources are associated with star-forming regions, and 9 are ambiguous cases. For the remaining 37 unIDed 4FGL sources, we could not identify any plausible counterpart in X-rays or they are too close to the Galactic Center. Finally, we outline the observational strategies and further improvements in the pipeline that can lead to more accurate classifications.

Although the MeV gamma-ray band is a promising energy-band window in astrophysics, the current situation of MeV gamma-ray astronomy significantly lags behind those of the other energy bands in angular resolution and sensitivity. An electron-tracking Compton camera (ETCC), a next-generation MeV detector, is expected to revolutionize the situation. An ETCC tracks each Compton-recoil electron with a gaseous electron tracker and determines the incoming direction of each gamma-ray photon; thus, it has a strong background rejection power and yields a better angular resolution than classical Compton cameras. Here, we study ETCC events in which the Compton-recoil electrons do not deposit all energies to the electron tracker but escape and hit the surrounding pixel scintillator array (PSA). We developed an analysis method for this untapped class of events and applied it to laboratory and simulation data. We found that the energy spectrum obtained from the simulation agreed with that of the actual data within a factor of 1.2. We then evaluated the detector performance using the simulation data. The angular resolution for the new-class events was found to be twice as good as in the previous study at the energy range 1.0--2.0~MeV, where both analyses overlap. We also found that the total effective area is dominated by the contribution of the double-hit events above an energy of 1.5~MeV. Notably, applying this new method extends the sensitive energy range with the ETCC from 0.2--2.1 MeV in the previous studies to up to 3.5~MeV. Adjusting the PSA dynamic range should improve the sensitivity in even higher energy gamma-rays. The development of this new analysis method would pave the way for future observations by ETCC to fill the MeV-band sensitivity gap in astronomy.

The census of open clusters has exploded in size thanks to data from the Gaia satellite. However, it is likely that many of these reported clusters are not gravitationally bound, making the open cluster census impractical for many scientific applications. We test different physically motivated methods for distinguishing between bound and unbound clusters, using them to create a cleaned cluster catalogue. We derived completeness-corrected photometric masses for 6956 clusters from our earlier work. Then, we used these masses to compute the size of the Roche surface of these clusters (their Jacobi radius) and distinguish between bound and unbound clusters. We find that only 5647 (79%) of the clusters from our previous catalogue are compatible with bound open clusters, dropping to just 11% of clusters within 250 pc. 3530 open clusters are in a strongly cut high quality sample. The moving groups in our sample show different trends in their size as a function of age and mass, suggesting that they are unbound and undergoing different dynamical processes. Our cluster mass measurements constitute the largest catalogue of Milky Way cluster masses to date, which we also use for further science. Firstly, we inferred the mass-dependent completeness limit of the open cluster census, showing that the census is complete within 1.8 kpc only for objects heavier than 230 M$_\odot$. Next, we derived a completeness-corrected age and mass function for our open cluster catalogue, including estimating that the Milky Way contains a total of $1.3 \times 10^5$ open clusters, only ~4% of which are currently known. Finally, we show that most open clusters have mass functions compatible with the Kroupa initial mass function. We demonstrate Jacobi radii for distinguishing between bound and unbound star clusters, and publish an updated star cluster catalogue with masses and improved cluster classifications. (abridged)

The development of low earth orbit (LEO) mega-constellation fundamentally threatens ground-based optical astronomical observations. To study the photometric properties of the LEO mega-constellations, we used the Xinglong 50 cm telescope to conduct a large-sample, high-precision, and multicolor target-tracking photometry of two typical LEO Mega-constellations: Starlink and OneWeb. Over a three-month observation period starting on 2022 January 1st, we collected 1,447 light curves of 404 satellites in four typical versions: Starlink v1.0, DarkSat, VisorSat, Starlink v1.5, and OneWeb. According to data statistics, Starlink v1.0 has the smallest median magnitude at clear and SDSS $gri$ band, and OneWeb is the dimmest bus. The brightness of Starlink v1.5 is slightly brighter than VisorSat. We construct a detailed photometric model with solar phase angle variations by calculating the illumination-visibility geometry based on the orbital parameters. Our data analysis shows that the solar phase angle is the significant characteristic which influencing Starlink satellites' brightness, but it is not sensitive to OneWeb satellites. VisorSat and Starlink v1.5 version, which are equipped with deployable visors, have significantly reduced scattered light compared to the previous Starlink v1.0 version. The multiband LOWESS and color-index are analyzed in characterizing the energy and color features of LEO mega-constellation satellites. This work found that the proportion of scattered sunlight mitigation achieved with VisorSat and Starlink v1.5 was 55.1\% and 40.4\%, respectively. The color index of different buses shows an evident clustering feature. Our observation and analysis could provide valuable quantitative data and photometric models, which can contribute to assessing the impact of LEO mega-constellations on astronomical observations.

Recently Loeb et al. (2024, "Recovery and Classification of Spherules from the Pacific Ocean Site of the CNEOS 2014 January 8 (IM1) Bolide", Res. Notes. Amer. Astron. Soc. 8, 39) reported the magnetic collection of millimeter-sized spherules from the seafloor near Papua New Guinea. About 22% had Mg/Si < 1/3 and were identified as a new "differentiated" variety of cosmic spherule ("D-type"). In a subset of 26 of these "D-type" spherules, 12 "BeLaU" spherules were found to be dominated by Fe and Al, marked by low Si and even lower Mg content, depletions of volatile species like Pb and Cs, and remarkable enrichments of Be, La, U, Ba, and other elements. Loeb et al. claimed these have exotic compositions different from other Solar System materials. We show that in fact samples with these compositions are not just found on Earth, they are from Earth; specifically, we identify them as microtektites of terrestrial lateritic sandstone. Based on the location of the sample site, we associate them with the Australasian tektite strewn field, generated 788 kyr ago by an impactor that melted and ejected ~10^8 tons of sandstone, including a lateritic layer, from Indochina. A tektite origin for the spherules is corroborated by their terrestrial Fe isotopic compositions and the compound, non-spherical nature of many of them, which preclude formation as ablation spherules from a bolide. Due to the restriction of laterites to the tropics, iron-rich tektites may be uncommon, but we predict they should comprise ~3% of the Australasian microtektites.

The supernova remnant SNR 0540-69.3, twin of the Crab Nebula, offers an excellent opportunity to study the continuum emission from a young pulsar and pulsar-wind nebula (PWN). We present observations taken with the VLT instruments MUSE and X-shooter in the wavelength range 3000-25,000 \r{A}, which allow us to study spatial variations of the optical spectra, along with the first near-infrared (NIR) spectrum of the source. We model the optical spectra with a power law (PL) $F_\nu\propto\nu^{-\alpha}$ and find clear spatial variations (including a torus-jet structure) in the spectral index across the PWN. Generally, we find spectral hardening toward the outer parts, from $\alpha\sim1.1$ to $\sim0.1$, which may indicate particle reacceleration by the PWN shock at the inner edge of the ejecta or alternatively time variability of the pulsar wind. The optical-NIR spectrum of the PWN is best described by a broken PL, confirming that several breaks are needed to model the full spectral energy distribution of the PWN, suggesting the presence of more than one particle population. Finally, subtracting the PWN contribution from the pulsar spectrum we find that the spectrum is best described with a broken-PL model with a flat and a positive spectral index, in contrast to the Crab pulsar that has a negative spectral index and no break in the optical. This might imply that pulsar differences propagate to the PWN spectra.

We determined shock speed in the Galactic supernova remnant Cygnus Loop, using proper motion of its optical filaments and the latest estimate for its distance. The proper motion was measured by comparing H$\alpha$ images of the remnant observed in two epochs: 1993 (Kitt Peak National Observatory) and 2018/2019 (National Astronomical Observatory Rozhen and Astronomical station Vidojevica). We derived shock speed for 35 locations along different filaments, which is twice as much as in earlier studies of north-eastern part of Cygnus Loop. For the first time, we have measured shock speed of radiative filaments in this region. Three of the analyzed locations where we measured proper motion of filaments are radiative, based on their presence in [SII] images from the second epoch. The other filaments are non-radiative. The speed we obtained for the non-radiative filaments is in the range of 240-650 $\mathrm{km\ s^{-1}}$, with an estimate for the uncertainty of 70 $\mathrm{km\ s^{-1}}$. These values are mostly in agreement with previous studies. The radiative filaments have lower speed of 100-160$\pm${70} $\mathrm{km\ s^{-1}}$, which is in agreement with the assumption that they are older in evolutionary terms. This clear distinction between the speed of the two types of filaments proves that the [SII] emission can be used for identifying radiative filaments in supernova remnants.

During the process of star formation, the dense gas undergoes significant chemical evolution leading to the emergence of a rich variety of molecules associated with hot cores and hot corinos. However, the physical and chemical conditions involved in this evolution are poorly constrained. We provide here a full inventory of the emission from complex organic molecules (COMs) to investigate the physical structure and chemical composition of six high-mass protostellar envelopes. We aim to investigate the conditions for the emergence of COMs in hot cores. We performed an unbiased spectral survey towards six infrared-quiet massive clumps between 159 GHz and 374 GHz with the APEX 12 m telescope. We detect up to 11 COMs, of which at least five COMs are detected towards all sources. Towards all the objects, most of the COM emission is found to be cold, with respect to the typical temperatures at which COMs are found, with a temperature of 30 K and extended with a size of ~0.3 pc. Although for our sample of young massive clumps the bulk of the gas has a cold temperature, we also detect emission from COMs originating from the immediate vicinity of the protostar revealing a compact and hot component of the envelope. Only three out of the six sources exhibit a hot gas component. We find a gradual emergence of the warm component in terms of size and temperature, together with an increasing molecular complexity, allowing us to establish an evolutionary sequence for our sample based on COMs. Our findings confirm that our sample of infrared-quiet massive clumps are in an early evolutionary stage during which the bulk of the gas is cold. The presence of COMs is found to be characteristic of these early evolutionary stages. We suggest that the emergence of hot cores is preceded by a phase in which mostly O-bearing COMs appear first with similar abundances to hot corinos albeit with larger source sizes.

We study the effects due to mismatches in passbands, polarization angles, and temperature and polarization calibrations in the context of the upcoming cosmic microwave background experiment Simons Observatory (SO). Using the SO multi-frequency likelihood, we estimate the bias and the degradation of constraining power in cosmological and astrophysical foreground parameters assuming different levels of knowledge of the instrumental effects. We find that incorrect but reasonable assumptions on the values of all the systematics examined here can have important effects in cosmological analyses, hence requiring marginalization approaches at likelihood level. When doing so, we find that the most relevant effect is due to bandpass shifts. When marginalizing over them, the posteriors of parameters describing astrophysical microwave foregrounds (such as radio point sources or dust) get degraded, while cosmological parameters constraints are not significantly affected. Marginalization over polarization angles with up to 0.25$^\circ$ uncertainty causes an irrelevant bias $\lesssim 0.05 \sigma$ in all parameters. Marginalization over calibration factors in polarization broadens the constraints on the effective number of relativistic degrees of freedom $N_\mathrm{eff}$ by a factor 1.2, interpreted here as a proxy parameter for non standard model physics targeted by high-resolution CMB measurements.

Clusters and their progenitors (protoclusters) at z = 2-4, the peak epoch of star formation, are ideal laboratories to study the formation process of both the clusters themselves and their member galaxies. However, a complete census of their member galaxies has been challenging due to observational difficulties. Here we present new JWST/NIRCam observations targeting the distant cluster CLJ1001 at z = 2.51 from the COSMOS-Webb program, which, in combination with previous narrow-band imaging targeting H-alpha emitters and deep millimeter surveys of CO-emitters, provide a complete view of massive galaxy assembly in CLJ1001. In particular, JWST reveals a population of massive, extremely red cluster members in the long-wavelength bands that were invisible in previous HST/F160W imaging (HST-dark). Based on this highly complete spectroscopic sample of member galaxies, we show that the spatial distribution of galaxies in CLJ1001 exhibits a strong central concentration, with the central galaxy density already resembling that of low-z clusters. Moreover, we reveal a ``top-heavy" stellar mass function for the star-forming galaxies (SFGs), with an overabundance of massive SFGs piled up in the cluster core. These features strongly suggest that CLJ1001 is caught in a rapid transition, with many of its massive SFGs likely soon becoming quiescent. In the context of cluster formation, these findings suggest that the earliest clusters form from the inside out and top to bottom, with the massive galaxies in the core assembling first, followed by the less-massive ones in the outskirts.

We estimate the neutrino flux from different kinds of galactic sources and compare it with the recently diffuse neutrino flux detected by IceCube. We find that the flux from these sources may contribute to ~ 20% of the IceCube neutrino flux. Most of the sources selected in this work populate the southern hemisphere, therefore a detector like KM3NeT could help in resolving the sources out of the observed diffused galactic neutrino flux.

The old nova and intermediate polar (IP) GK Persei underwent one of its recurrent dwarf nova (DN) outbursts in 2018. We proposed monitoring it in UV and X-rays with the Neil Gehrels Swift Observatory, starting less than six days after the eruption, until 16 days after the eruption ended. For the first time we could follow the decay to minimum light UV and X-rays. We present the timing and spectral analysis, comparing the results with the previous outbursts and with the quiescent status. We confirm the spin modulation in X-rays with a period 351.325(9) s, only in the 2-10 keV range. The period was not detected in the 0.3-2 keV range and in the UV band, suggesting that the soft portion of the X-ray spectrum in GK Per does not originate near the poles, but in a wind or circumstellar material. The amplitude of the modulation was less prominent than in 2015, a fact that seems correlated with a lower average mass accretion rate. The spectral fits are consistent with a mass accretion rate increasing by a factor of 2 from rise to maximum and decreasing during the return to minimum, following the trend of the modulation amplitude. The maximum plasma temperature is higher than the Swift XRT energy range of 0.3-10 keV, thus it is not well constrained, but our spectral fits indicate that it may have varied irregularly during the outburst.

We present an implicit method for solving the diffusion equation for the evolution of the dust fraction in the terminal velocity approximation using dust-as-mixture smoothed particle hydrodynamics (SPH). The numerical scheme involves casting the dust diffusion equation into implicit form, rearranging into its resolvent cubic equation and solving analytically. This method is relevant for small grains that are tightly coupled to the gas, such as sub-micron dust grains in the interstellar medium or millimetre-sized dust grains in protoplanetary discs. The method avoids problems with the variable used to evolve the dust fraction becoming negative when evolved explicitly and is fast and accurate, avoiding the need for dust stopping time limiters and significantly reducing computational expense. Whilst this method is an improvement over using the explicit terminal velocity approximation method, as with any dust-as-mixture method it still fails to give accurate solutions in the limit of large (weakly coupled) grains.

A key feature of active galactic nuclei (AGN) is their variability across all wavelengths. Typically, AGN vary by a few tenths of a magnitude or more over periods lasting from hours to years. By contrast, extreme variability of AGN -- large luminosity changes that are a significant departure from the baseline variability -- are known as AGN flares. These events are rare and their timescales poorly constrained, with most of the literature focusing on individual events. It has been suggested that extreme AGN variability including flares can provide insights into the accretion processes in the disk. With surveys such as the Legacy Survey of Space and Time (LSST) promising millions of transient detections per night in the coming decade, there is a need for fast and efficient classification of AGN flares. The problem with the systematic detection of AGN flares is the requirement to detect them against a stochastically variable baseline; the ability to define a signal as a significant departure from the ever-present variability is a statistical challenge. Recently, Gaussian Processes (GPs) have revolutionised the analysis of time-series data in many areas of astronomical research. They have, however, seen limited uptake within the field of transient detection and classification. Here we investigate the efficacy of Gaussian Processes to detect AGN flares in both simulated and real optical light curves. We show that GP analysis can successfully detect AGN flares with a false-positive rate of less than seven per cent, and we present examples of AGN light curves that show extreme variability.

Extreme emission line galaxies (EELGs) are typically characterized by high equivalent widths (EWs) which are driven by elevated specific star formation rates (sSFR) in low-mass galaxies with subsolar metallicities and little dust. Such extreme systems are rare in the local universe, but the number density of EELGs increases with redshift. Such starburst galaxies are currently presumed to be the main drivers of hydrogen reionization over 5.5<z<15, which serves to motivate many of the searches for high-z EELGs. We aim to characterize the physical properties of a sample of ~730 EELGs at 4<z<9 photometrically selected from the CEERS survey using JWST/NIRCam. We validate our method and demonstrate the main physical properties of a subset of EELGs using NIRSpec spectra. We create synthetic NIRCam observations of EELGs using empirical templates based on ~2000 local metal-poor starbursts to select EELGs based on color-color criteria. We study their properties based on SED fitting and flux excess from emission lines in the photometric filters. Our sample has a mean stellar mass of $10^{7.84}$Msun with high sSFRs with a mean value of $10^{-7.03}$ yr$^{-1}$. We consider a delayed-$\tau$ model for the star formation history and find our sample of EELGs are young with a mean value of the time after the onset of star formation of 45Myr. We find that they have similar line ratios to local metal-poor starbursts with high log([OIII]/H$\beta$)>0.4-1 which indicates that star formation may be the dominant source of ionization. Based on the photometric fluxes, we find an increase of EW([OIII]+H$\beta$) with sSFR and $\Sigma_{SFR}$, and a decrease with age and stellar mass. The sample of EELGs can reach $\Sigma_{SFR}>$10Msun yr$^{-1}$kpc$^{-2}$ which indicate they are strong candidates of LyC leakers. Another indirect indicator is the high values of O32>5 that can be reached for some galaxies in the sample.

The temperature structure of a giant planet was traditionally thought to be an adiabat assuming convective mixing homogenizes entropy. The only in-situ measurement made by the Galileo Probe detected a near-adiabatic temperature structure within one of Jupiter's 5$\mu$m hot spots with small but definite local departures from adiabaticity. We analyze Juno's microwave observations near Jupiter's equator (0 ~ 5$^o$N) and find that the equatorial temperature structure is best characterized by a stable super-adiabatic temperature profile rather than an adiabatic one. Water is the only substance with sufficient abundance to alter the atmosphere's mean molecular weight and prevent dynamic instability if a super-adiabatic temperature gradient exists. Thus, from the super-adiabaticity, our results indicate a water concentration (or the oxygen to hydrogen ratio) of about 4.9 times solar with a possible range of 1.5 ~ 8.3 times solar in Jupiter's equatorial region.

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate

Radiation transport plays a crucial role in star formation models, as certain questions within this field cannot be accurately addressed without taking it into account. Given the high complexity of the interstellar medium from which stars form, numerical simulations are frequently employed to model the star formation process. This study reviews recent methods for incorporating radiation transport into star formation simulations, discussing them in terms of the used algorithms, treatment of radiation frequency dependence, the interaction of radiation with the gas, and the parallelization of methods for deployment on supercomputers. Broadly, the algorithms fall into two categories: (i) moment-based methods, encompassing the flux-limited diffusion approximation, M1 closure, and variable Eddington tensor methods, and (ii) methods directly solving the radiation transport equation, including forward and reverse ray tracing, characteristics-based methods, and Monte Carlo techniques. Beyond discussing advantages and disadvantages of these methods, the review also lists recent radiation hydrodynamic codes implemented the described methods.

Radio-interferometric (RI) imaging entails solving high-resolution high-dynamic range inverse problems from large data volumes. Recent image reconstruction techniques grounded in optimization theory have demonstrated remarkable capability for imaging precision, well beyond CLEAN's capability. These range from advanced proximal algorithms propelled by handcrafted regularization operators, such as the SARA family, to hybrid plug-and-play (PnP) algorithms propelled by learned regularization denoisers, such as AIRI. Optimization and PnP structures are however highly iterative, which hinders their ability to handle the extreme data sizes expected from future instruments. To address this scalability challenge, we introduce a novel deep learning approach, dubbed ``Residual-to-Residual DNN series for high-Dynamic range imaging'. R2D2's reconstruction is formed as a series of residual images, iteratively estimated as outputs of Deep Neural Networks (DNNs) taking the previous iteration's image estimate and associated data residual as inputs. It thus takes a hybrid structure between a PnP algorithm and a learned version of the matching pursuit algorithm that underpins CLEAN. We present a comprehensive study of our approach, featuring its multiple incarnations distinguished by their DNN architectures. We provide a detailed description of its training process, targeting a telescope-specific approach. R2D2's capability to deliver high precision is demonstrated in simulation, across a variety of image and observation settings using the Very Large Array (VLA). Its reconstruction speed is also demonstrated: with only few iterations required to clean data residuals at dynamic ranges up to 105, R2D2 opens the door to fast precision imaging. R2D2 codes are available in the BASPLib library on GitHub.

We update the field-level inference code KARMMA to enable tomographic forward-modelling of shear maps. Our code assumes a lognormal prior on the convergence field, and properly accounts for the cross-covariance in the lensing signal across tomographic source bins. We use mock weak lensing data from N-body simulations to validate our mass-mapping forward model by comparing our posterior maps to the input convergence fields. We find that KARMMA produces more accurate reconstructions than traditional mass-mapping algorithms. More-over, the KARMMA posteriors reproduce all statistical properties of the input density field we tested -- one- and two-point functions, and the peak and void number counts -- with $\lesssim~10\%$ accuracy. Our posteriors exhibit a small bias that increases with decreasing source redshift, but these biases are small compared to the statistical uncertainties of current (DES) cosmic shear surveys. Finally, we apply KARMMA to Dark Energy Survey Year 3 (DES-Y3) weak lensing data, and verify that the two point shear correlation function $\xi_+$ is well fit by the correlation function of the reconstructed convergence field. This is a non-trivial test that traditional mass mapping algorithms fail. The code is publicly available at https://github.com/Supranta/KaRMMa.git. KARMMA DES-Y3 mass maps are publicly available at https://zenodo.org/records/10672062.

Some cosmic microwave background (CMB) data allow a cosmological scenario in which the free streaming of neutrinos is delayed until close to matter-radiation equality. Interestingly, recent analyses have revealed that large-scale structure (LSS) data also align with this scenario, discarding the possibility of an accidental feature in the CMB sky and calling for further investigation into the free-streaming nature of neutrinos. By assuming a simple representation of self-interacting neutrinos, we investigate whether this nonstandard scenario can accommodate a consistent cosmology for both the CMB power spectra and the large-scale distribution of galaxies simultaneously. Employing three different approaches - a profile likelihood exploration, a nested sampling method, and a heuristic Metropolis-Hasting approximation - we exhaustively explore the parameter space and demonstrate that galaxy data exacerbates the challenge already posed by the Planck polarization data for this nonstandard scenario. We find that the most conservative value of the Bayes factor disfavors the interactions among neutrinos over a $\Lambda$CDM + $N_\mathrm{eff}$ + $\sum m_\nu$ model with odds of $23:1000$ and that the difficulty of simultaneously fitting the galaxy and CMB data relates to the so-called $S_8$ discrepancy. Our analysis not only emphasizes the need to consider a broader range of phenomenologies in the early Universe but also highlights significant numerical and theoretical challenges ahead in uncovering the exact nature of the feature observed in the data or, ultimately, confirming the standard chronological evolution of the Universe.

The Near-infrared Spectrograph (NIRSpec) on the James Webb Space Telescope is uniquely suited to studying galaxies in the distant Universe with its combination of multi-object capabilities and sensitivity over a large range in wavelength (0.6-5.3 microns). Here we present the NIRSpec Wide survey, part of the NIRSpec Instrument Science Team's Guaranteed Time Observations, using NIRSpec's microshutter array to obtain spectra of more than 3200 galaxies at $z>1$ at both low- and high-resolution ($R\approx100$ and 2700) for a total of 105 hours. With 31 pointings covering $\approx$320 arcmin$^2$ across the five CANDELS fields with exquisite ancillary photometry from the Hubble Space Telescope, the NIRSpec Wide survey represents a fast and efficient way of using JWST to probe galaxies in the early Universe. Pointing centers are determined to maximize the observability of the rarest, high-value sources. Subsequently, the microshutter configurations are optimized to observe the maximum number of "census" galaxies with a selection function based primarily on HST/F160W magnitude, photometric/slitless grism redshift, and predicted \ha\ flux tracing the bulk of the galaxy population at cosmic noon ($z_{\rm med}=2.0$). We present details on the survey strategy, the target selection, an outline of the motivating science cases, and discuss upcoming public data releases to the community.

(abridged) Ultra-high energy cosmic rays (UHECRs) are extremely energetic charged particles with energies surpassing $10^{18}$ eV. Their sources remain elusive, obscured by deflections caused by the Galactic magnetic field (GMF). This challenge is further complicated by our limited understanding of the three-dimensional structure of the GMF, as current GMF observations consist primarily of quantities integrated along the line-of-sight (LOS). Nevertheless, data from upcoming stellar polarisation surveys along with Gaia's stellar parallax data are expected to yield local GMF measurements.. In this work, we employ methods of Bayesian statistical inference in order to sample the posterior distribution of the GMF within part of the Galaxy. By assuming a known rigidity and arrival direction of an UHECR, we backtrack its trajectory through various GMF configurations drawn from the posterior distribution. Our objective is to rigorously evaluate our algorithm's performance in scenarios that closely mirror the setting of expected future applications. In pursuit of this, we condition the posterior to synthetic integrated LOS measurements of the GMF, in addition to synthetic local POS-component measurements. In this proof of concept work, we assume the ground truth to be a magnetic field produced by a dynamo simulation of the Galactic ISM. Our results demonstrate that for all locations of the observed arrival direction on the POS, our algorithm is able to substantially update our knowledge on the original arrival direction of UHECRs with rigidity $E/Z = 5 \times 10^{19}$ eV, even in the case of complete absence of LOS information. If integrated data is included in the inference, then the regions of the celestial sphere where the maximum error occurs diminishes greatly. Even in those regions the maximum error is diminished by a factor of about $3$ in the specific setting studied.

We present a new model of string inflation driven by a blow-up K\"ahler modulus of type IIb compactifications with a potential generated by string loops. Slow-roll is naturally realized thanks to the fact that the blow-up mode is a leading-order flat direction lifted by string loops which are unavoidable and generate a plateau at large field values. We check that throughout the whole inflationary dynamics the effective field theory is under control. We perform a phenomenological analysis determining the exact number of efoldings by studying the post-inflationary evolution. We determine the values of the microscopic parameters which lead to agreement with CMB data, together with the prediction of a tensor-to-scalar ratio of order $r\sim 10^{-5}$.

Astrophysical uncertainties in dark matter direct detection experiments are typically addressed by parametrizing the velocity distribution in terms of a few uncertain parameters that vary around some central values. Here we propose a method to optimize over all velocity distributions lying within a given distance measure from a central distribution. We discretize the dark matter velocity distribution as a superposition of streams, and use a variety of information divergences to parametrize its uncertainties. With this, we bracket the limits on the dark matter-nucleon and dark matter-electron scattering cross sections, when the true dark matter velocity distribution deviates from the commonly assumed Maxwell-Boltzmann form. The methodology pursued is general and could be applied to other physics scenarios where a given physical observable depends on a function that is uncertain.

In this work we present a new approach to produce spectroscopic constants and model first-principles synthetic spectra for all molecules of astrophysical interest. We have generalized our previous diatomic molecule simulation framework, employing Transition-Optimised Shifted Hermite (TOSH) theory, thereby enabling the modelling of polyatomic rotational constants for molecules with three or more atoms. These capabilities, are now provided by our new code Epimetheus. As a first validation of our approach, we confront our predictions and assess their accuracy against the well-studied triatomic molecule, ozone 666 ($^{16}$O$_3$), in addition to eight of its potential isotopomers: ozone 668 ($^{16}$O$^{16}$O$^{18}$O), 686 ($^{16}$O$^{18}$O$^{16}$O), 667 ($^{16}$O$^{16}$O$^{17}$O), 676 ($^{16}$O$^{17}$O$^{16}$O), 688 ($^{16}$O$^{18}$O$^{18}$O), 868 ($^{18}$O$^{16}$O$^{18}$O), 888 ($^{18}$O$_3$), and 777 ($^{17}$O$_3$). We then assess the accuracy of these rotational constants using the Epimetheus data in our code Pandora, and generate synthetic molecular spectra. The ozone spectra presented here are purely infrared and not Raman. Epimetheus builds upon the work from our previous code Prometheus, which used the TOSH theory to account for anharmonicity for the fundamental $\nu=0 \rightarrow \nu=1$ band, going further to now account for triatomic molecules. This is combined with thermal profile modeling for the rotational transitions. We have found that this extended method performs promisingly, typically approximating the spectroscopic constants and spectra well. Some issues do arise depending on the symmetry group of the ozone isotopomer. In general, we show that Epimetheus can provide the data to produce appreciable molecular spectra, to help drive future high-resolution studies.