Papers for Wednesday, Sep 11 2024

On December 1980, supernova 1980N was discovered in NGC 1316, a galaxy of the Fornax cluster. Three months later, supernova 1981D was observed in the same galaxy. The light curves of these two supernovae Ia were found to be virtually identical, suggesting that they are images of the same event, the delay between them being due to strong gravitational lensing. If so, as anticipated by Sjur Refsdal, the distance to the lens can be determined accurately, namely, 90 $\pm$ 1 kpc, meaning that it belongs to the outer halo of the Milky Way. Interestingly, there is another pair of possible images in the Pantheon+ sample, namely, supernovae 2013aa and 2017cbv, the distance to the lens being 702 $\pm$ 1 kpc, that is, nearly the same as the distance to the Andromeda galaxy. In both cases, given the relatively large angle of deviation of the supernova light by the lens, namely, 271" and 325", respectively, the lens has to be a compact object, with a mass to radius ratio over 150 M$_\odot$ R$_\odot^{-1}$. It is likely to be an ultra massive white dwarf.

Schwarzschild dynamical models are now regularly employed in large surveys of galaxies in the local and distant Universe to derive information on galaxies' intrinsic properties such as their orbital structure and their (dark matter and stellar) mass distribution. Comparing the internal orbital structures and mass distributions of galaxies in the distant Universe with simulations is key to understanding what physical processes are responsible for shaping galaxy properties. However it is first crucial to understand whether observationally derived properties are directly comparable with intrinsic ones in simulations. To assess this, we build Schwarzschild dynamical models for MUSE-like IFS cubes (constructed to be like those obtained by the MAGPI survey) of 75 galaxies at z ~ 0.3 from the Eagle simulations. We compare the true particle-derived properties with the galaxies' model-derived properties. In general, we find that the models can recover the true galaxy properties qualitatively well, with the exception of the enclosed dark matter, where we find a median offset of 48%, which is due to the assumed NFW profile not being able to reproduce the dark matter distribution in the inner region of the galaxies. We then compare our model-derived properties with Schwarzschild models-derived properties of observed MAGPI galaxies and find good agreement between MAGPI and Eagle: the majority of our galaxies (57%) have non-oblate shapes within 1 effective radius. More triaxial galaxies show higher fractions of hot orbits in their inner regions and tend to be more radially anisotropic.

3C 264 is one of the few FRI radio galaxies with detected TeV emission. It is a low-luminosity AGN (LLAGN) and is generally associated with a radiatively inefficient accretion flow (RIAF). Earlier multiwavelength studies suggest that the X-ray emission originates from a jet. However, the possibility that the RIAF can significantly contribute to the X-rays cannot be ruled out. In particular, hard X-ray emission $\gtrsim$10 keV has never been detected, making it challenging to distinguish between X-ray models. Here we report a NuSTAR detection up to 25 keV from 3C 264. We also present subpixel deconvolved Chandra images to resolve jet emission down to ~0.2 arcsec from the center of the unresolved X-ray core. Together with a simultaneous Swift observation, we have constrained the dominant hard X-ray emission to be from its unresolved X-ray core, presumably in its quiescent state. We found evidence of a cutoff in the energy around 20 keV, indicating that at least some of the X-rays from the core can be attributed to the RIAF. The Comptonization model suggests an electron temperature of about 15 keV and an optical depth ranging between 4 and 7, following the universality of coronal properties of black hole accretion. The cutoff energy or electron temperature of 3C 264 is the lowest among those of other LLAGNs. The detected hard X-ray emission is at least an order of magnitude higher than that predicted by synchrotron self-Compton models introduced to explain $\gamma$-ray and TeV emission, suggesting that the synchrotron electrons might be accelerated to higher energies than previously thought.

Recent detections of gravitational waves from mergers of binary black holes (BBHs) with pre-merger source-frame individual masses in the so-called upper mass-gap, expected due to (pulsational) pair instability supernova ((P)PISN), have created immense interest in the astrophysical production of high-mass black holes (BHs). Previous studies show that high-mass BHs may be produced via repeated BBH mergers inside dense star clusters. Alternatively, inside dense star clusters, stars with unusually low core-to-envelope mass ratios can form via mergers of high-mass stars, which then can avoid (P)PISN, but produce high-mass BHs via mass fallback. We simulate detailed star-by-star multi-physics models of dense star clusters using the Monte Carlo cluster evolution code, CMC, to investigate the role of primordial binary fraction among high-mass stars (>=15 Msun) on the formation of high-mass BHs. We vary the high-mass stellar binary fraction (fb_15_prime) while keeping all other initial properties, including the population of high-mass stars, unchanged. We find that the number of high-mass BHs, as well as the mass of the most massive BH formed via stellar core-collapse are proportional to fb_15_prime. In contrast, there is no correlation between fb_15_prime and the number of high-mass BHs formed via BH-BH mergers. Since the total production of high-mass BHs is dominated by BH-BH mergers in old clusters, the overall number of high-mass BHs produced over the typical lifetime of globular clusters is insensitive to fb_15_prime. Furthermore, we study the differences in the demographics of BH-BH mergers as a function of fb_15_prime.

Stars form in dense cores within molecular clouds and newly formed stars influence their natal environments. How stellar feedback impacts core properties and evolution is subject to extensive investigation. We performed a hierarchical clustering (dendrogram) analysis of a STARFORGE simulation modelling a giant molecular cloud to identify gas overdensities (cores) and study changes in their radius, mass, velocity dispersion, and virial parameter with respect to stellar feedback. We binned these cores on the basis of the fraction of gas affected by protostellar outflows, stellar winds, and supernovae and analysed the property distributions for each feedback bin. We find that cores that experience more feedback influence are smaller. Feedback notably enhances the velocity dispersion and virial parameter of the cores, more so than it reduces their radius. This is also evident in the linewidth-size relation, where cores in higher feedback bins exhibit higher velocities than their similarly sized pristine counterparts. We conclude that stellar feedback mechanisms, which impart momentum to the molecular cloud, simultaneously compress and disperse the dense molecular gas.

High-resolution millimetre-imaging of protoplanetary discs has revealed many containing rings and gaps. These rings can contain large quantities of dust, often in excess of 10M$_\oplus$, providing prime sites for efficient and rapid planet formation. Rapid planet formation will produce high accretion luminosities, heating the surrounding disc. We investigate the importance of a planetary embryo's accretion luminosity by simulating the dynamics of the gas and dust in a dust ring, accounting for the energy liberated as a resident planetary embryo accretes. The resulting heating alters the flow structure near the planet, increasing the accretion rate of large, millimetre-to-centimetre-sized dust grains. We show how this process varies with the mass of dust in the ring and the local background gas temperature, demonstrating that the thermal feedback always acts to increase the planet's mass. This increase in planet mass is driven primarily by the formation of vortices, created by a baroclinic instability once the accreting planet heats the disc significantly outside its Hill radius. The vortices can then migrate with respect to the planet, resulting in a complex interplay between planetary growth, gap-opening, dust trapping and vortex dynamics. Planets formed within dust traps can have masses that exceed the classical pebble isolation mass, potentially providing massive seeds for the future formation of giant planets. Once pebble accretion ceases, the local dust size distribution is depleted in large grains, and much of the remaining dust mass is trapped in the system's L$_5$ Lagrange point, providing potentially observable signatures of this evolution.

A growing population of long-period radio transients has been discovered and their physical origin is still up to debate. Recently, a new such source named ILT J1101 + 5521 was discovered, which is in a white dwarf (WD) -- M dwarf (MD) binary system, with the observed 125.5 min period being identified as the orbital period and the radio emission phase coinciding with the conjunction configuration when the MD is at the far end. We suggest that the radio emission properties of the system can be well explained within the framework of the unipolar inductor magnetic interaction model between the magnetized WD and the MD with low magnetization, with the electron cyclotron maser being the most likely radiation mechanism. This mechanism is similar to that of Jupiter decametric emission due to Jupiter-Io interaction. We suggest that this mechanism can interpret at least some long-period radio transients, especially the ultra-long period sub-population.

The Epoch of Reionization Spectrometer (EoR-Spec) is an upcoming Line Intensity Mapping (LIM) instrument designed to study the evolution of the early universe (z = 3.5 to 8) by probing the redshifted [CII] 158 $\mu$m fine-structure line from aggregates of galaxies. The [CII] emission is an excellent tracer of star formation since it is the dominant cooling line from neutral gas heated by OB star light and thus can be used to probe the reionization of the early Universe due to star formation. EoR-Spec will be deployed on Prime-Cam, a modular direct-detection receiver for the 6-meter Fred Young Submillimeter Telescope (FYST), currently under construction by CPI Vertex Antennentechnik GmbH and to be installed near the summit of Cerro Chajnantor in the Atacama Desert. This instrument features an image plane populated with more than 6500 Microwave Kinetic Inductance Detectors (MKIDs) that are illuminated by a 4-lens optical design with a cryogenic, scanning Fabry-Perot Interferometer (FPI) at the pupil of the optical system. The FPI is designed to provide a spectral resolving power of $R\sim100$ over the full spectral range of 210--420 GHz. EoR-Spec will tomographically survey the E-COSMOS and E-CDFS fields with a depth of about 4000 hours over a 5 year period. Here we give an update on EoR-Spec's final mechanical/optical design and the current status of fabrication, characterization and testing towards first light in 2026.

One of the primary science goals of the Habitable Worlds Observatory (HWO) as defined by the Astro2020 decadal survey is the imaging of the first Earth-like planet around a Sun-like star. A key technology gap towards reaching this goal are the development of ultra-low-noise photon counting detectors capable of measuring the incredibly low count rates coming from these planets which are at contrasts of $\sim 1 \times 10^{-10}$. Superconducting energy-resolving detectors (ERDs) are a promising technology for this purpose as, despite their technological challenges, needing to be cooled below their superconducting transition temperature ($< 1\mathrm{K}$), they have essentially zero read noise, dark current, or clock-induced charge, and can get the wavelength of each incident photon without the use of additional throughput-reducing filters or gratings that spread light over many pixels. The use of these detectors on HWO will not only impact the science of the mission by decreasing the required exposure times for exo-Earth detection and characterization, but also in a wavefront sensing and control context when used for starlight suppression to generate a dark zone. We show simulated results using both an EMCCD and an ERD to ``dig a dark zone'' demonstrating that ERDs can achieve the same final contrast as an EMCCD in about half of the total time. We also perform a simple case study using an exposure time calculator tool called the Error Budget Software (EBS) to determine the required integration times to detect water for HWO targets of interest using both EMCCDs and ERDs. This shows that once a dark zone is achieved, using an ERD can decrease these exposure times by factors of 1.5--2 depending on the specific host star properties.

The number density of galaxy clusters across mass and redshift has been established as a powerful cosmological probe. Cosmological analyses with galaxy clusters traditionally employ scaling relations. However, many challenges arise from this approach as the scaling relations are highly scattered, may be ill-calibrated, depend on the cosmology, and contain many nuisance parameters with low physical significance. In this paper, we use a simulation-based inference method utilizing artificial neural networks to optimally extract cosmological information from a shallow X-ray survey of galaxy clusters, solely using count rates (CR), hardness ratios (HR), and redshifts. This procedure enables us to conduct likelihood-free inference of cosmological parameters $\Omega_{\mathrm{m}}$ and $\sigma_8$. We analytically generate simulations of galaxy cluster distribution in a CR, HR space in multiple redshift bins based on totally random combinations of cosmological and scaling relation parameters. We train Convolutional Neural Networks (CNNs) to retrieve the cosmological parameters from these simulations. We then use neural density estimation (NDE) neural networks to predict the posterior probability distribution of $\Omega_{\mathrm{m}}$ and $\sigma_8$ given an input galaxy cluster sample. The 1 $\sigma$ errors of our density estimator on one of the target testing simulations are 1000 deg$^2$: 15.2% for $\Omega_{\mathrm{m}}$ and 10.0% for $\sigma_8$; 10000 deg$^2$: 9.6% for $\Omega_{\mathrm{m}}$ and 5.6% for $\sigma_8$. We also compare our results with Fisher analysis. We demonstrate, as a proof of concept, that it is possible to calculate cosmological predictions of $\Omega_{\mathrm{m}}$ and $\sigma_8$ from a galaxy cluster population without explicitly computing cluster masses and even, the scaling relation coefficients, thus avoiding potential biases resulting from such a procedure. [abridged]

In this paper, we obtain new measurements of the angular homogeneity scale ($\theta_H$) from the BOSS DR12 and eBOSS DR16 catalogs of Luminous Red Galaxies of the Sloan Digital Sky Survey. Considering the flat $\Lambda$CDM model, we use the $\theta_H(z)$ data to constrain the matter density parameter ($\Omega_{m0}$) and the Hubble constant ($H_{0}$). We find $H_0 = 65.7 \pm 7.0$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_{m0}>0.293$. By combining the $\theta_H$ measurements with current Baryon Acoustic Oscillations (BAO) and Type Ia Supernova (SN) data, we obtain $H_{0}= 69.9^{+4.9}_{-4.5}$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_{m0} = 0.294^{+0.013}_{-0.015}$ ($\theta_H$ + BAO) and $H_{0}=70.7^{+5.2}_{-4.1}$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_{m0}=0.299 \pm 0.022$ ($\theta_H$ + SN). We show that $\theta_H$ measurements help break the BAO and SN degeneracies concerning $H_0$, as they do not depend on the sound horizon scale at the drag epoch or the SN absolute magnitude value obtained from the distance ladder method. Hence, despite those constraints being loose compared to other probes, $\theta_H$ data may provide an independent cosmological probe of $H_0$ in light of the Hubble tension. For completeness, we also forecast the constraining power of future $\theta_H$ data via Monte Carlo simulations. Considering a relative error of the order of 1$\%$, we obtain competitive constraints on $\Omega_{m0}$ and $H_0$ ($\approx 5\%$ error) from the joint analysis with current SN and BAO measurements.

One of the most compelling pieces of evidence of the Hot Big Bang model is the realisation and confirmation that some nuclides were created shortly after the Big Bang. This process is referred to as Big Bang nucleosynthesis (or, sometimes, primordial nucleosynthesis), and is the end-product of putting neutrons and protons in a hot, expanding Universe. Big Bang nucleosynthesis currently provides our earliest test of cosmology, and it is the only experiment currently designed that is simultaneously sensitive to all four known fundamental forces: the gravitational force, the electromagnetic force, the strong force and the weak force. Our theoretical understanding of Big Bang nucleosynthesis and the measurement of the primordial abundances together represents one of the strongest pillars of the standard cosmological model. In this chapter, we will develop an intuitive understanding of Big Bang nucleosynthesis, discuss modern calculations of this process, and provide a summary of the current state-of-the-art measurements that have been made. Overall, Big Bang nucleosynthesis is in remarkable agreement with various cosmological probes, and it is this agreement that serves to strengthen our confidence in the general picture of cosmology that we have today.

Primordial black holes (PBHs) would induce non-Gaussianity in the cosmic microwave background (CMB) by sourcing recombination perturbations spatially modulated by relative velocities between PBHs and the baryons they accrete. The leading non-Gaussian signatures are non-vanishing connected 4-point correlation functions, or trispectra. Earlier, we computed the CMB temperature trispectrum, and forecasted Planck to be more sensitive to it than to changes in the CMB temperature power spectrum for light enough PBHs. Excitingly, accreting PBHs would also induce non-Gaussianity in CMB polarization, and source both E and B modes, which we compute in this paper. We first calculate linear-response perturbations to the tensor-valued photon distribution function sourced by a general spatially-varying ionization history, and apply our results to accreting PBHs. We then compute linear-order perturbations to the temperature and polarization 2-point functions sourced by inhomogeneities in recombination due to accreting PBHs; we find them to be negligible relative to their counterparts sourced by homogeneous perturbations to the ionization history. Lastly, we compute all CMB trispectra including temperature, E- and B-mode polarization at linear order in the PBH abundance. We forecast that including polarization data in a 4-point-function analysis would only increase Planck's sensitivity to accreting PBHs by a factor ~2 relative to using temperature alone. As a consequence, we find that a search for PBHs using all temperature and polarization trispectra with Planck data would mostly not be competitive with current bounds from temperature and polarization power spectra. In contrast, we forecast that a CMB Stage-4 experiment would gain significant sensitivity to accreting PBHs through a 4-point-function search, in particular through the contributions of parity-odd trispectra including one B-mode field.

Earth's geodynamo has operated for over 3.5 billion years. The magnetic field is currently powered by thermocompositional convection in the outer core, which involves the release of light elements and latent heat as the inner core solidifies. However, since the inner core nucleated no more than 1.5 billion years ago, the early dynamo could not rely on these buoyancy sources. Given recent estimates of the thermal conductivity of the outer core, an alternative mechanism may be required to sustain the geodynamo prior to nucleation of the inner core. One possibility is a silicate dynamo operating in a long-lived basal magma ocean. Here, we investigate the structural, thermal, buoyancy, and magnetic evolution of an Earth-like terrestrial planet. Using modern equations of state and melting curves, we include a time-dependent parameterization of the compositional evolution of an iron-rich basal magma ocean. We combine an internal structure integration of the planet with energy budgets in a coupled core, basal magma ocean, and mantle system. We determine the thermocompositional convective stability of the core and the basal magma ocean, and assess their respective dynamo activity using entropy budgets and magnetic Reynolds numbers. Our conservative nominal model predicts a transient basal magma ocean dynamo followed by a core dynamo after 1 billion years. The model is sensitive to several parameters, including the initial temperature of the core-mantle boundary, the parameterization of mantle convection, the composition of the basal magma ocean, the radiogenic content of the planet, as well as convective velocity and magnetic scaling laws. We use the nominal model to constrain the range of basal magma ocean electrical conductivity and core thermal conductivity that sustain a dynamo.

As mass-losing asymptotic giant branch (AGB) stars evolve to planetary nebulae (PNe), the mass outflow geometries transform from nearly spherical to extreme aspherical. The physical mechanisms governing this transformation are widely believed to be linked to binarity and the associated production of disks and fast jets during transitional (post-AGB) evolutionary stages. We are carrying out a systematic ALMA survey ($P$re-planet$A$ry $N$ebulae high-angular-res$O$lution su$R$vey with $A$L$MA$ or PANORAMA) of a representative sample of bipolar and multipolar post-AGB objects. We have obtained high angular-resolution (0".1-0".4) observations of the CO(3--2) and/or 6--5 emission in order to probe the spatio-kinematic structure of the collimated outflows and the central disk/torii. The results are remarkable, generally showing the presence of bipolar or multipolar high-velocity outflows, dense toroidal waists, and in one case, a geometrically-thin circular ring around the central bipolar nebula. A high degree of point-symmetry characterizes the morphology of the mass ejecta. In this contribution, we present these and other highlights from our survey. We aim to use 2D/3D radiative transfer modeling in order to derive accurate outflow momenta, masses and mass-loss rates for our sample, and build hydrodynamical models that can explain the observed spatio-kinematic structures. These results will then be used to distinguish between different classes of PN-shaping binary interaction models.

The active galactic nucleus (AGN) feeding and feedback process in the centers of galaxy clusters and groups is still not well understood. NGC5044 is the ideal system in which to study AGN feedback. It hosts the largest known reservoir of cold gas in any cool-core galaxy group, and features several past epochs of AGN feedback imprinted as cavities in the X-ray bright intragroup medium (IGrM), as well as parsec scale jets. We present Submillimeter Array (SMA) and Karl G. Jansky Very Large Array (VLA) high frequency observations of NGC5044 to assess the time variability of the mm-waveband emission from the accretion disk, and quantify the Spectral Energy Distribution (SED) from the radio to sub-millimeter band. The SED is well described by advection dominated accretion flow (ADAF) model and self-absorbed jet emission from an aging plasma with tau ~1kyr. We find a characteristic variability timescale of 150 days, which constrains the ADAF emission region to about 0.1pc, and the magnetic field to 4.7mG in the jets and and 870G in the accretion disk. A longer monitoring/sampling will allow to understand if the underlying process is truly periodic in nature.

Over the past decade and a half, adoption of Bayesian inference in pulsar timing analysis has led to increasingly sophisticated models. The recent announcement of evidence for a stochastic background of gravitational waves by various pulsar timing array projects highlighted Bayesian inference as a central tool for parameter estimation and model selection. Despite its success, Bayesian inference is occasionally misused in the pulsar timing community. A common workflow is that the data is analyzed in multiple steps: a first analysis of single pulsars individually, and a subsequent analysis of the whole array of pulsars. A mistake that is then sometimes introduced stems from using the posterior distribution to craft the prior for the analysis of the same data in a second step, a practice referred to in the statistics literature as ``circular analysis.'' This is done to prune the model for computational efficiency. Multiple recent high-profile searches for gravitational waves by pulsar timing array (PTA) projects have this workflow. This letter highlights this error and suggests that Spike and Slab priors can be used to carry out model averaging instead of model selection in a single pass. Spike and Slab priors are proved to be equal to Log-Uniform priors.

Ganymede's aurora are the product of complex interactions between its intrinsic magnetosphere and the surrounding Jovian plasma environment and can be used to derive both atmospheric composition and density. In this study, we analyzed a time-series of Ganymede's optical aurora taken with Keck I/HIRES during eclipse by Jupiter on 2021-06-08 UTC, one day after the Juno flyby of Ganymede. The data had sufficient signal-to-noise in individual 5-minute observations to allow for the first high cadence analysis of the spatial distribution of the aurora brightness and the ratio between the 630.0 and 557.7 nm disk-integrated auroral brightnesses -- a quantity diagnostic of the relative abundances of O, O$_2$ and H$_2$O in Ganymede's atmosphere. We found that the hemisphere closer to the centrifugal equator of Jupiter's magnetosphere (where electron number density is highest) was up to twice as bright as the opposing hemisphere. The dusk (trailing) hemisphere, subjected to the highest flux of charged particles from Jupiter's magnetosphere, was also consistently almost twice as bright as the dawn (leading) hemisphere. We modeled emission from simulated O$_2$ and H$_2$O atmospheres during eclipse and found that if Ganymede hosts an H$_2$O sublimation atmosphere in sunlight, it must collapse on a faster timescale than expected to explain its absence in our data given our current understanding of Ganymede's surface properties.

We present the detection and analysis of GRB 231115A, a candidate extragalactic magnetar giant flare (MGF) observed by Fermi/GBM and localized by INTEGRAL to the starburst galaxy M82. This burst exhibits distinctive temporal and spectral characteristics that align with known MGFs, including a short duration and a high peak energy. Gamma-ray analyses reveal significant insights into this burst, supporting conclusions already established in the literature: our time-resolved spectral studies provide further evidence that GRB 231115A is indeed a MGF. Significance calculations also suggest a robust association with M82, further supported by a high Bayes factor that minimizes the probability of chance alignment with a neutron star merger. Despite extensive follow-up efforts, no contemporaneous gravitational wave or radio emissions were detected. The lack of radio emission sets stringent upper limits on possible radio luminosity. Constraints from our analysis show no fast radio bursts (FRBs) associated with two MGFs. X-ray observations conducted post-burst by Swift/XRT and XMM/Newton provided additional data, though no persistent counterparts were identified. Our study underscores the importance of coordinated multi-wavelength follow-up and highlights the potential of MGFs to enhance our understanding of short GRBs and magnetar activities in the cosmos. Current MGF identification and follow-up implementation are insufficient for detecting expected counterparts; however, improvements in these areas may allow for the recovery of follow-up signals with existing instruments. Future advancements in observational technologies and methodologies will be crucial in furthering these studies.

Apophis' current trajectory takes it safely past our planet at a distance of several Earth radii on 2029 April 13. Here the possibility is considered that Apophis could collide with a small asteroid, like the ones that frequently and unpredictably strike Earth, and the resulting perturbation of its trajectory. The probability of an impact that could significantly displace Apophis relative to its keyholes is found to be less than 1 in $10^6$, requiring a delta-v greater than 0.3 mm/s, while for an impact that could significantly displace Apophis compared to its miss distance in 2029 it is less than 1 in $10^9$, requiring a delta-v greater than 5 cm/s. These probabilities are below the usual thresholds considered by asteroid impact warning systems. Apophis is in the daytime sky and unobservable from mid-2021 to 2027. It will be challenging to determine from single night observations in 2027 if Apophis has moved on the target plane enough to enter a dangerous keyhole, as the deviation from the nominal ephemeris might be only a few tenths of an arcsecond. An impending Earth impact would, however, be signalled clearly in most cases by deviations of tens of arcseconds of Apophis from its nominal ephemeris in 2027. Thus most of the impact risk could be retired by a single observation of Apophis in 2027, though a minority of cases present some ambiguity and are discussed in more detail. Charts of the on-sky position of Apophis under different scenarios are presented for quick assessment by observers.

I have used the third data release of the Gaia mission to improve the reliability and completeness of membership samples in the beta Pic moving group (BPMG) and other nearby associations with ages of 20-50 Myr (Sco Body, Carina, Columba, chi1 For, Tuc-Hor, IC 2602, IC 2391, NGC 2547). I find that Carina, Columba, and chi1 For are physically related and coeval, and that Carina is the closest fringe of a much larger association. Similarly, Tuc-Hor and IC 2602 form a coeval population that is spatially and kinematically continuous. Both results agree with hypotheses from Gagne et al. (2021). I have used the new catalogs to study the associations in terms of their initial mass functions, X-ray emission, ages, and circumstellar disks. For instance, using the model for Li depletion from Jeffries et al. (2023), I have derived an age of 24.7+0.9/-0.6 Myr for BPMG, which is similar to estimates from previous studies. In addition, I have used infrared photometry from the Wide-field Infrared Survey Explorer to check for excess emission from circumstellar disks among the members of the associations, which has resulted in a dramatic increase in the number of known disks around M stars at ages of 30-50 Myr and a significant improvement in measurements of excess fractions for those spectral types and ages. Most notably, I find that the W3 excess fraction for M0-M6 initially declines with age to a minimum in BPMG <0.015), increases to a maximum in Carina/Columba chi1 For (0.041+0.009/-0.007, 34 Myr), and declines again in the oldest two associations (40-50 Myr). The origin of that peak and the nature of the M dwarf disks at >20 Myr are unclear.

We analyze $\textit{JWST}$ CEERS NIRCam images to present {the first estimate} of the observed fraction and properties of bars out to $z \sim 4$. We analyze a sample of 1770 galaxies with stellar mass $M_\star > 10^{10} M_\odot$ at $0.5 \leq z \leq 4$ and identify barred galaxies via ellipse fits and visual classification of both F200W and F444W images. Our results apply mainly to bars with projected semi-major axis $a_{\rm bar}$ $> 1.5 $ kpc ($\sim$ 2 $\times$ PSF in F200W images) that can be robustly traced by ellipse fits. For such bars, the {observed} bar fraction at $z\sim$ 2-4 is low ($\lesssim 10\%$), and they appear to be emerging at least as early as $z\sim 4$ when the Universe was $\sim$ 13\% of its present age. At $z\sim$ 2-4, compared to our results, TNG50 simulations {predict} a significantly larger bar fraction due to a large population of small bars with $a_{\rm bar}$ $< 1.5$ kpc {that we cannot robustly detect}. If such a population exists, the true bar fraction may be significantly higher than our results. At $z \ge 1.5$, many barred galaxies show nearby neighbors, suggesting bars may be tidally triggered. {From $z \sim 4$ to $z \sim 0.5$, the observed bar fraction, average projected bar length, and projected bar strength rise.} Our results highlight the early emergence and evolution of barred galaxies and the rising importance of bar-driven secular evolution from $z \sim$4 to today.

We have used the Five-hundred-meter Aperture Spherical radio Telescope (FAST) to make a blind ultra-deep survey for neutral hydrogen (HI). We present the complete results from the first of six fields (FUDS0). This observation of 95 hours allowed us to achieve a high sensitivity ($\sim 50~\mu$Jy beam$^{-1}$) and a high frequency resolution (22.9 kHz) over an area of 0.72 deg$^2$. We detected 128 galaxies in HI distributed over the redshift range of $0<z<0.4$ with HI masses in the range of $6.67 \leq \log(M_{\rm HI}/h_{70}^{-2} \rm M_\odot) \leq 10.92$, and three faint high-velocity clouds (HVCs) with peak column density of $N_{\rm HI} \leq 3.1 \times 10^{17}$ cm$^{-2}$. Of the galaxies, 95 are new detections and six have $z > 0.38$, where no unlensed HI emission has previously been directly detected. Estimates of completeness and reliability are presented for the catalog. Consistency of continuum and HI flux estimates with NVSS and AUDS, respectively, confirms the accuracy of calibration method and data reduction pipeline developed for the full FUDS survey.

Stars form in the Galaxy with a wide range in mass. If the mass is below 7% of the Sun's, then the object does not become hot enough for stable hydrogen burning. These substellar objects are called brown dwarfs. Maps of the sky at infrared wavelengths have found large numbers of brown dwarfs. However only 24 objects have been found (as of April 2017) that are cold enough to be classified as "Y dwarfs": these have atmospheres that are cooler than 500 K (or 200 C, 400 F) and have masses only 5 - 20 times that of Jupiter. The coolest Y dwarf currently known, discovered in 2014, has a temperature around freezing, has a mass of about 5 Jupiter masses, and is only 2 pc away from the Sun. These small and cold objects are faint and difficult to find. This chapter describes the discovery and characterization of the Y dwarfs. Finding more of these very cold planet-like brown dwarfs will require an as-yet unplanned space mission mapping large areas of sky at wavelengths around 5 microns.

The Fe II/Mg II emission line flux ratio in quasar spectra serves as a proxy for the relative Fe to alpha-element abundances in the broad line regions of quasars. Due to the expected different enrichment timescales of the two elements, they can be used as a cosmic clock in the early Universe. We present a study of the Fe II/Mg II ratios in a sample of luminous quasars exploiting high-quality near-IR spectra taken primarily by the XQR-30 program with VLT XSHOOTER. These quasars have a median bolometric luminosity of log(L_bol[erg s^-1])~47.3 and cover a redshift range of z=6.0-6.6. The median value of the measured Fe II/Mg II ratios is ~7.9 with a normalized median absolute deviation of ~2.2. In order to trace the cosmic evolution of Fe II/Mg II in an unbiased manner, we select two comparison samples of quasars with similar luminosities and high-quality spectra from the literature, one at intermediate redshifts (z=3.5-4.8) and the other at low redshifts (z=1.0-2.0). We perform the same spectral analysis for all these quasars, including the usage of the same iron template, the same spectral fitting method, and the same wavelength fitting windows. We find no significant redshift evolution in the Fe II/Mg II ratio over the wide redshift range from z=1 to 6.6. The result is consistent with previous studies and supports the scenario of a rapid iron enrichment in the vicinity of accreting supermassive black holes at high redshift.

We explore the co-evolution of dark matter halos, their central galaxies, and central supermassive black holes (SMBHs) using the IllustrisTNG (TNG) simulation. We find that the evolutionary histories of individual galaxies in the $M_{\rm BH}$-$M_*$ plane can be decomposed into four distinct phases, separated by three transition points. We identify the driving processes of galaxy evolution within each phase and derive the conditions necessary and sufficient for transitions to subsequent phases. The first phase is dominated by star formation, with its duration primarily determined by the mass of the SMBH seed and the surrounding gas environment. The second phase is characterized by rapid SMBH growth, and the transition to the next phase occurs when the thermal-mode feedback of active galactic nucleus (AGN) can unbind gas from the galaxy. The third phase involves self-regulation of the SMBH, and the transition to the quenched phase occurs when the kinetic-mode feedback of AGN counterbalances gas cooling within the subhalo. The final phase is dominated by mergers. We investigate the use of scaling relations among different mass components and evolutionary phases to understand processes implemented in TNG and other simulations, and discuss how current and forthcoming observations can be used to constrain models.

We used the Five-hundred-meter Aperture Spherical radio Telescope (FAST) to search for the molecular emissions in the L-band between 1.0 and 1.5 GHz toward four comets, C/2020 F3 (NEOWISE), C/2020 R4 (ATLAS), C/2021 A1 (Leonard), and 67P/Churyumov-Gerasimenko during or after their perihelion passages. Thousands of molecular transition lines fall in this low-frequency range, many attributed to complex organic or prebiotic molecules. We conducted a blind search for the possible molecular lines in this frequency range in those comets and could not identify clear signals of molecular emissions in the data. Although several molecules have been detected at high frequencies of great than 100 GHz in comets, our results confirm that it is challenging to detect molecular transitions in the L-band frequency ranges. The non-detection of L-band molecular lines in the cometary environment could rule out the possibility of unusually strong lines, which could be caused by the masers or non-LTE effects. Although the line strengths are predicted to be weak, for FAST, using the ultra-wide bandwidth receiver and improving the radio frequency interference environments would enhance the detectability of those molecular transitions at low frequencies in the future.

We explore the stability and fate of gravitational triple systems comprising a central massive body and a tight binary of less massive pairs. Our present purpose is two fold; (1) to improve the Hill-type stability criterion for the binary in those systems, and (2) to examine the fate of the triple systems after the binary break-up, with particular attention to the effects of the eccentricities of the inner and outer orbits. We perform direct Newtonian N-body simulations over much longer integration times than previous studies, which is essential to determine the eventual fate of those systems statistically in a reliable fashion. We obtain an empirical fitting formula of the binary stability boundary that incorporates effects of the inner and outer eccentricities, the mutual inclination of the inner and outer orbits, and the mass ratios of the three bodies. We also find that those triple systems are stable for a much longer timescale after the binary break-up, and that their final fates (ejection of the outer body, merger to the central massive body, and collision of two less massive bodies) are very sensitive to the initial outer eccentricity.

Temperate exoplanets between the sizes of Earth and Neptune, known as "sub-Neptunes", have emerged as intriguing targets for astrobiology. It is unknown whether these planets resemble Earth-like terrestrial worlds with a habitable surface, Neptune-like giant planets with deep atmospheres and no habitable surface, or something exotic in between. Recent JWST transmission spectroscopy observations of the canonical sub-Neptune K2-18 b revealed ~1% CH4, ~1% CO2, and a non-detection of CO in the atmosphere. While previous studies have proposed that the observed atmospheric composition could help constrain the lower atmosphere conditions and determine the interior structure of sub-Neptunes like K2-18 b, the possible interactions between the atmosphere and a hot, supercritical water ocean at its base remain unexplored. In this work, we investigate whether a global supercritical water ocean, resembling a planetary-scale hydrothermal system, can explain these observations on K2-18 b-like sub-Neptunes through equilibrium aqueous geochemical calculations. We find that the observed atmospheric CH4/CO2 ratio implies a minimum ocean temperature of ~715 K, whereas the corresponding CO/CO2 ratio allows ocean temperatures up to ~1060 K. These results indicate that a global supercritical water ocean on K2-18 b is plausible. While life cannot survive in this ocean, this work represents the first step towards understanding how a global supercritical water ocean may influence observable atmospheric characteristics on volatile-rich sub-Neptunes. Future observations with better constrained NH3 and CO mixing ratios could further help distinguish between possible interior compositions of K2-18 b.

High-precision and high-cadence photometric surveys such as Kepler or TESS are making huge progress not only in the detection of new extrasolar planets but also in the study of a great number of variable stars. This is the case for central stars of planetary nebulae (PNe), which have similarly benefited from the capabilities of these missions, increasing the number of known binary central stars and helping us to constrain the relationship between binarity and the complex morphologies of their host PNe. In this paper, we analyse the TESS light curves of a large sample of central stars of PNe with the aim of detecting signs of variability that may hint at the presence of short-period binary nuclei. We analysed 62 central stars of true, likely, or possible PNe and modelled the detected variability through an MCMC approach accounting for three effects: reflection, ellipsoidal modulations, and Doppler beaming. Among the 62 central stars, only 38 are amenable for this study. The remaining 24 show large contamination from nearby sources preventing an optimal analysis. Also, eight targets are already known binary central stars, which we revisit here with the new high precision of the TESS data. In addition, we find that 18 further central stars show clear signs of periodic variability in the TESS data, probably resulting from different physical effects compatible with the binary scenario. We propose them as new candidate binary central stars. We also discuss the origin of the detected variability in each particular case by using the TESS_localize algorithm. Finally, 12 targets show no or only weak evidence of variability at the sensitivity of TESS. Our study demonstrates the power of space-based photometric surveys in searching for close binary companions of central stars of PNe.

We present a plausible and coherent view of the evolution of the protosolar disk that is consistent with the cosmochemical constraints and compatible with observations of other protoplanetary disks and sophisticated numerical simulations. The evidence that high-temperature condensates, CAIs and AOAs, formed near the protosun before being transported to the outer disk can be explained by either an early phase of vigorous radial spreading of the disk, or fast transport of these condensates from the vicinity of the protosun towards large disk radii via the protostellar outflow. The assumption that the material accreted towards the end of the infall phase was isotopically distinct allows us to explain the observed dichotomy in nucleosynthetic isotopic anomalies of meteorites and leads to intriguing predictions on the isotopic composition of refractory elements in comets. When the infall of material waned, the disk started to evolve as an accretion disk. Initially, dust drifted inwards, shrinking the radius of the dust component to ~ 45 au, probably about 1/2 of the width of the gas component. Then structures must have emerged, producing a series of pressure maxima in the disk which trapped the dust on My timescales. This allowed planetesimals to form at radically distinct times without changing significantly of isotopic properties. There was no late accretion of material onto the disk via streamers. The disk disappeared in ~5 Myr, as indicated by paleomagnetic data in meteorites. In conclusion, the evolution of the protosolar disk seems to have been quite typical in terms of size, lifetime, and dust behavior, suggesting that the peculiarities of the Solar system with respect to extrasolar planetary system probably originate from the chaotic nature of planet formation and not at the level of the parental disk.

Mass distribution of black holes identified through X-ray emission suggests a paucity of black holes in the mass range of 3 to 5 solar masses. Modified theories have been devised to explain this mass gap, and it is suggested that natal kicks during supernova explosion can more easily disrupt binaries with lower mass black holes. Although recent LIGO observations reveal the existence of compact remnants within this mass gap, the question of whether low-mass black holes can exist in binaries remains a matter of debate. Such a system is expected to be noninteracting without X-ray emission, and can be searched for using radial velocity and astrometric methods. Here we report Gaia DR3 3425577610762832384, a wide binary system including a red giant star and an unseen object, exhibiting an orbital period of approximately 880 days and near-zero eccentricity. Through the combination of radial velocity measurements from LAMOST and astrometric data from Gaia DR2 and DR3 catalogs, we determine a mass of $3.6^{+0.8}_{-0.5}$ $M_{\odot}$ of the unseen component. This places the unseen companion within the mass gap, strongly suggesting the existence of binary systems containing low-mass black holes. More notably, the formation of its surprisingly wide circular orbit challenges current binary evolution and supernova explosion theories.

Quantum evaporation of a black hole is conventionally studied semiclassically by assuming self-similarity of the black hole throughout the evaporation process. However, its validity was recently questioned, and the lifetime of a black hole is conjectured to be much extended by the memory burden effect. It gives rise to the possibility that the primordial black holes (PBHs) lighter than $10^{10}$ grams are the dark matter in the Universe. To probe such PBH dark matter, we study gravitational waves (GWs) induced by primordial curvature perturbations that produced the PBHs. We find $\Omega_\text{GW}(f_\text{peak})h^2 = 7 \times 10^{-9}$ with the peak frequency $f_\text{peak} = 1\times 10^{3} \, (M_\text{PBH}/(10^{10}\,\mathrm{g}))^{-1/2}\, \mathrm{Hz}$, and the induced GWs associated with the PBH dark matter whose initial mass is greater than about $10^7$ grams can be tested by future observations such as Cosmic Explorer. Furthermore, the scenario can be in principle confirmed by detecting another GW signal from the mergers of PBHs, which leads to high-frequency GWs with $f_\text{peak} = 2 \times 10^{27}\, (M_\text{PBH, ini}/(10^{10}\, \mathrm{g}))^{-1} \, \mathrm{Hz} $. On the other hand, the induced GW signals stronger than expected would contradict the dark matter abundance and exclude the memory burden effect.

Supernova remnants (SNRs) carry vast amounts of mechanical and radiative energy that heavily influence the structural, dynamical, and chemical evolution of galaxies. To this day, more than 300 SNRs have been discovered in the Milky Way, exhibiting a wide variety of observational features. However, existing classification schemes are mainly based on their radio morphology. In this work, we introduce a novel unsupervised deep learning pipeline to analyse a representative subsample of the Galactic SNR population ($\sim$ 50% of the total) with the aim of finding a connection between their multi-wavelength features and their physical properties. The pipeline involves two stages: (1) a representation learning stage, consisting of a convolutional autoencoder that feeds on imagery from infrared and radio continuum surveys (WISE 22$\mu$m, Hi-GAL 70 $\mu$m and SMGPS 30 cm) and produces a compact representation in a lower-dimensionality latent space; and (2) a clustering stage that seeks meaningful clusters in the latent space that can be linked to the physical properties of the SNRs and their surroundings. Our results suggest that this approach, when combined with an intermediate uniform manifold approximation and projection (UMAP) reprojection of the autoencoded embeddings into a more clusterable manifold, enables us to find reliable clusters. Despite a large number of sources being classified as outliers, most clusters relate to the presence of distinctive features, such as the distribution of infrared emission, the presence of radio shells and pulsar wind nebulae, and the existence of dust filaments.

An ongoing challenge for radio-based detectors of high-energy cosmic particles is the accurate description of radio signal propagation in natural non-uniform media. For radio signals originating from extensive air showers, the current state of the art simulations often implicitly assume straight-line signal propagation. This while the refraction due to a non-uniform atmosphere is expected to have an effect on the received signal and associated reconstruction that is currently not completely understood for the most inclined geometries. Here, we present a study regarding the validity of assuming straight-line signal propagation when simulating radio emission associated with very inclined air shower geometries. To this end, the calculation of the electric field based on the end-point formalism used in CoREAS was improved by use of tabulated ray tracing data. We find that including ray curvature effects into the end-point formalism calculation introduces changes of up to a few percent in fluence for frequencies up to 1.2 GHz and zenith angles up to 88°.

Ly$\alpha$ is the transition to the ground state from the first excited state of hydrogen (the most common element). Resonant scattering of this line by neutral hydrogen greatly impedes its emergence from galaxies, so the fraction of galaxies which show Ly$\alpha$ is a tracer of the neutral fraction of the intergalactic medium (IGM), and thus the history of reionization. In previous works, we used early JWST/NIRSpec data from the JWST Advanced Deep Extragalactic Survey (JADES) to classify and characterise Ly$\alpha$ emitting galaxies (LAEs). This survey is now approaching completion, and the current sample is nearly an order of magnitude larger. From a sample of 784 galaxies in JADES at $4.0<z<14.3$, we find evidence for Ly$\alpha$ emission in 145 sources. We reproduce the previously found correlation between Ly$\alpha$ escape fraction (\fesc) - Ly$\alpha$ rest-frame equivalent width (\rew) and the negative correlation between Ly$\alpha$ velocity offset - \fesc. Both \fesc and \rew decrease with redshift ($z\gtrsim5.5$), indicating the progression of reionization on a population scale. Our data are used to demonstrate an increasing IGM transmission of Ly$\alpha$ from $z\sim14-6$. We measure the completeness-corrected fraction of LAEs ($X_{Ly\alpha}$) from $z=4-9.5$. An application of these $X_{Ly\alpha}$ values to the results of cosmological models suggests a high neutral fraction at $z=7$ ($\rm X_{HI}=0.81_{-0.10}^{+0.07}$), likely suggesting the need for models with updated \rew distributions (based on comparison to other works). This large sample of LAEs and the completeness correction we have detailed will be paramount for unbiased population studies of galaxies in the EoR.

Energy-dependence of X-ray Fourier power spectral states and the characteristic frequencies of the Band-Limited Noise (BLN) components have been seen in the hard state and intermediate states of black hole X-ray binaries. Here we report our analysis of the \emph{Insight}-HXMT observations of the black hole transient MAXI J1820$+$070 during its 2018 outburst when the source was brightest in hard X-rays. We found opposite trends of the low-frequency ($<$ 0.1 Hz) and the high-frequency ($>$ 10 Hz) BLN components, i.e., decreasing vs. increasing in frequency with increasing photon energy up to beyond 200 keV, respectively. This establishes an apparent two-way broadening of the power plateau formed by multiple BLNs in the power spectra towards higher photon energies. The trend of increasing characteristic frequency of the highest BLN component with increasing photon energy has been interpreted as due to that the corresponding seed photons which are up-scatted to relatively higher energies originate in a region relatively more central in the corona previously. Following the same framework, the decreasing trend of the characteristic frequency of the low-frequency BLN component with increasing photon energy can be interpreted as due to that the corresponding seed photons which are up-scattered to higher photon energies originate from further out in the disk flow but on the opposite side of the central corona as to the observer. The opposite trends then implies that the the plateau in the power spectra formed by the multiple BLNs represents the radial extension of the accretion disk that contributes seed photons which produce the observed BLNs; the higher the photon energy is, the wider the power plateau and the smaller the fractional variability are, probably approaching to a Power-Law Noise (PLN) seen in the soft state.

We introduce the Contour Analysis Tool (CAT), a Python toolkit aimed at identifying and analyzing structural elements in density maps. CAT employs various contouring techniques, including the lowest-closed contour (LCC), linear and logarithmic Otsu thresholding, and average gradient thresholding. These contours can aid in foreground and background segmentation, providing natural limits for both, as well as edge detection and structure identification. Additionally, CAT provides image processing methods such as smoothing, background removal, and image masking. The toolkit features an interactive suite of controls designed for Jupyter environments, enabling users to promptly visualize the effects of different methods and parameters. We describe, test, and demonstrate the performance of CAT, highlighting its potential use cases. CAT is publicly available on GitHub, promoting accessibility and collaboration.

We have developed a new set of Application-Specific Integrated Circuits (ASICs) of the TARGET family (CTC and CT5TEA), designed for the readout of signals from photosensors in cameras of Imaging Atmospheric Cherenkov Telescopes (IACTs) for ground-based gamma-ray astronomy. We present the performance and design details. Both ASICs feature 16 channels, with CTC being a Switched-Capacitor Array (SCA) sampler at 0.5 to 1 GSa/s with a 16,384 sample deep storage buffer, including the functionality to digitize full waveforms at arbitrary times. CT5TEA is its companion trigger ASIC (though may be used on its own), which provides trigger information for the analog sum of four (and 16) adjacent channels. Since sampling and triggering takes place in two separate ASICs, the noise due to interference from the SCA is suppressed, and allows a minimal trigger threshold of $\leq$ 2.5 mV (0.74 photo electrons (p.e.)) with a trigger noise of $\leq$ 0.5 mV (0.15 p.e.). For CTC, a maximal input voltage range from $-$0.5 V up to 1.7 V is achieved with an effective bit range of $>$ 11.6 bits and a baseline noise of 0.7 mV. The cross-talk improved to $\leq$ 1% over the whole $-$3 dB bandwidth of 220 MHz and even down to 0.2% for 1.5 V pulses of 10 ns width. Not only is the performance presented, but a temperature-stable calibration routine for pulse mode operation is introduced and validated. The resolution is found to be $\sim$ 2.5% at 33.7 mV (10 p.e.) and $\leq$ 0.3% at 337 mV (100 p.e.) with an integrated non-linearity of $<$ 1.6 mV. Developed for the Small-Sized Telescope (SST) and Schwarzschild-Couder Telescope (SCT) cameras of the Cherenkov Telescope Array Observatory (CTAO), CTC and CT5TEA are deployed for both prototypes and shall be integrated into the final versions.

Context. (4015) Wilson-Harrington (hereafter, WH) was discovered as a comet in 1949 but has a dynamical property consistent with that of a near-Earth asteroid. Although there is a report that the 1949 activity is associated with an ion tail, the cause of the activity has not yet been identified. Aims. This work aims to reveal the mysterious comet-like activity of the near-Earth asteroid. Methods. We conducted new polarimetric observations of WH from May 2022 to January 2023, reanalyses of the photographic plate images taken at the time of its discovery in 1949, and dust tail simulation modelings, where the dust terminal velocity and ejection epoch are taken into account. Results. We found that this object shows polarization characteristics similar to those of low-albedo asteroids. We derived the geometric albedo ranging from pV = 0.076 +- 0.010 to pV = 0.094 +- 0.018 from our polarimetry (the values vary depending on the data used for fitting and the slope-albedo relationship coefficients). In addition, the 1949 image showed an increase in brightness around the nucleus. Furthermore, we found that the color of the tail is consistent with sunlight, suggesting that the 1949 activity is associated with dust ejection. From the dust tail analysis, ~9 x 10^5 kg of material was ejected episodically at a low velocity equivalent to or even slower than the escape velocity. Conclusions. We conclude that WH is most likely an active asteroid of main belt origin and that the activity in 1949 was likely triggered by mass shedding due to fast rotation.

Third generation and future upgrades of current gravitational-wave detectors will present exquisite sensitivities which will allow to detect a plethora of gravitational wave signals. Hence, a new problem to be solved arises: the detection and parameter estimation of overlapped signals. The problem of separating and identifying two signals that overlap in time, space or frequency is something well known in other fields (e.g. medicine and telecommunication). Blind source separation techniques are all those methods that aim at separating two or more unknown signals. This article provides a methodological review of the most common blind source separation techniques and it analyses whether they can be successfully applied to overlapped gravitational wave signals or not, while comparing the limits and advantages of each method.

In the instant preheating scenario efficient particle production occurs immediately following the period of inflationary expansion in the early Universe. We demonstrate that instant preheating predicts unique gravitational wave (GW) signals arising from two distinct origins. One source is the bremsstrahlung GWs produced through the decay of superheavy particles, an inevitable consequence of instant preheating. The other is GWs generated from the nonlinear dynamics of the inflaton and coupled scalar fields. Using numerical simulations, we show that the peak of the GW spectrum shifts depending on the coupling constants of the theory. The detection of these dual GW signatures, characteristic of instant preheating, provides novel opportunities for probing the dynamics of the early Universe.

Context. Supernova remnants (SNRs) are the late stages of supernovae before their merging into the surrounding medium. Oxygen-rich supernova remnants represent a rare subtype with strong visible light oxygen emission. Aims. We present a new method to detect SNRs exploiting the capabilities of modern visible-light integral-field units based on the shapes of the SNR emission lines. Methods. We search for unresolved shocked regions with broadened emission lines using the medium-resolution integral-field spectrograph MUSE on the Very Large Telescope. The spectral resolving power allows shocked emission sources to be differentiated from photoionised sources based on the linewidths. Results. We find 307 supernova remnants, including seven O-rich SNRs. For all O-rich SNRs, we observe the [O III]{\lambda}{\lambda}4959,5007 emission doublet. In addition, we observe emissions from [O I]{\lambda}{\lambda}6300,6364, [O II]{\lambda}{\lambda}7320,7330, H{\alpha}+[N II]{\lambda}6583 and [S II]{\lambda}{\lambda}6717,6731 to varying degrees. The linewidths for the O-rich SNRs are generally broader than the rest of the SNRs in the sample of this article. The oxygen emission complexes are reminiscient of SNR 4449-1 and some long-lasting SNe. For the O-rich SNRs, we also search for counterparts in archival data of other telescopes; we detect X-ray and mid-IR counterparts for a number of remnants. Conclusions. We have shown efficacy of the method to detect SNRs presented in this article. In addition, the method is also effective in detecting the rare O-rich SNRs, doubling the sample size in the literature. The origin of O-rich SNRs and their link to specific SN types or environments is still unclear, but further work into this new sample will unquestionably help us shed light on these rare remnants.

This work describes the instrumental error budget for space-based measurements of the absolute flux of the sky synchrotron spectrum at frequencies below the ionospheric cutoff (<20 MHz). We focus on an architecture using electrically short dipoles onboard a small satellite. The error budget combines the contributions of the dipole dimensions, plasma noise, stray capacitance, and front-end amplifier input impedance. We treat the errors using both a Monte Carlo error propagation model and an analytical method. This error budget can be applied to a variety of experiments and used to ultimately improve the sensing capabilities of space-based electrically short dipole instruments. The impact of individual uncertainty components, particularly stray capacitance, is explored in more detail.

Massive black hole binary (MBHB) mergers will be detectable in large numbers by the Lisa Interferometer Space Antenna (LISA), which will thus provide new insights on how they form via repeated dark matter (DM) halo and galaxy mergers. Here we present a simple analytical model to generate a population of MBHB mergers based on a theoretical prescription that connects them to DM halo mergers. The high flexibility of our approach allows us to explore the broad and uncertain range of MBH seeding and growth mechanisms, as well as the different effects behind the interplay between MBH and galactic astrophysics. Such a flexibility is fundamental for the successful implementation and optimisation of the hierarchical Bayesian parameter estimation approach that here we apply to the MBHB population of LISA for the first time. Our inferred population hyper-parameters are chosen as proxies to characterise the MBH--DM halo mass scaling relation, the occupation fraction of MBHs in DM halos and the delay between halo and MBHB mergers. We find that LISA will provide tight constraints at the lower-end of the MBH-halo scaling relation, well complementing EM observations which are biased towards large masses. Furthermore, our results suggest that LISA will constrain some features of the MBH occupation fraction at high redshift, as well as merger time delays of the order of a few hundreds of Myr, opening the possibility to constrain dynamical evolution time scales such as the dynamical friction. The analysis presented here constitutes a first attempt at developing a hierarchical Bayesian inference approach to the LISA MBHB population, opening the way for several further improvements and investigations.

We study the perturbed-from-synchronous librational state of a double asteroid, modeled by the Full Two Rigid Body Problem (F2RBP), with primary emphasis on deriving analytical formulas which describe the system's evolution after deflection by a kinetic impactor. To this end, both a linear and nonlinear (canonical) theory are developed. We make the simplifying approximations (to be relaxed in a forthcoming paper) of planar binary orbit and axisymmetric shape of the primary body. To study the effect of a DART-like hit on the secondary body, the momentum transfer enhancement parameter $\beta$ is introduced and retained as a symbolic variable throughout all formulas derived, either by linear or nonlinear theory. Our approach can be of use in the context of the analysis of the post impact data from kinetic impactor missions, by providing a precise modeling of the impactor's effect on the seconadry's librational state as a function of $\beta$.

Stage IV large scale structure surveys are promising probes of gravity on cosmological scales. Due to the vast model-space in the modified gravity literature, model-independent parameterisations represent useful and scalable ways to test extensions of $\Lambda$CDM. In this work we use a recently validated approach of computing the non-linear $3\times 2$pt observables in modified gravity models with a time-varying effective gravitational constant $\mu$ and a gravitational slip $\eta$ that is binned in redshift to produce Fisher forecasts for an LSST Y10-like survey. We also include in our modelling an effective nulling scheme for weak-lensing by applying the BNT transformation that localises the weak-lensing kernel enabling well-informed scale cuts. We show that the combination of improved non-linear modelling and better control of the scales that are modelled/cut yields high precision constraints on the cosmological and modified gravity parameters. We find that 4 redshift bins for $\mu$ of width corresponding to equal incremental $\Lambda$CDM growth is optimal given the state-of-the-art modelling and show how the BNT transformation can be used to mitigate the impact of small-scale systematic effects, such as baryonic feedback.

Exploiting IRAM 30 m CO spectroscopy, we find that SDSS post-merger galaxies display gas fractions and depletion times enhanced by 25-50%, a mildly higher CO excitation, and standard molecular-to-atomic gas ratios, compared to non-interacting galaxies with similar redshift, stellar mass ($M_{\star}$) and star-formation rate (SFR). To place these results in context, we compile further samples of interacting or starbursting galaxies, from pre-coalescence kinematic pairs to post-starbursts, carefully homogenising gas mass, $M_{\star}$ and SFR measurements in the process. We explore systematics by duplicating our analysis for different SFR and $M_{\star}$ estimators, finding good qualitative agreement in general. Molecular gas fractions and depletion times are enhanced in interacting pairs, albeit less than for post-mergers. Among all samples studied, gas fraction and depletion time enhancements appear largest in young (a few 100 Myr) post-starbursts. While there is only partial overlap between post-mergers and post-starbursts, this suggests that molecular gas reservoirs are boosted throughout most stages of galaxy interactions, plausibly due to torque-driven inflows of halo gas and gas compression. The gas fraction and depletion time offsets of mergers and post-starbursts anti-correlate with their distance from the galaxy main sequence $\Delta({\rm MS})$, evidencing the role of SFE in driving the high SFRs of the strongest starbursts. Post-starbursts display the steepest dependency of gas fraction and SFE-offsets on $\Delta({\rm MS})$, with an evolving normalisation that reflects gas reservoir depletion over time. Our multi-sample analysis paints a coherent picture of the starburst-merger throughout the low-z merger sequence. It reconciles contradictory literature findings by highlighting that gas fraction enhancements and SFE variations both play their part in merger-driven star formation.

Intrinsic alignments (IAs) of galaxies/halos observed via galaxy imaging survey, combined with redshift information, offer a novel probe of cosmology as a tracer of tidal force field of large-scale structure. In this paper, we present a perturbation theory based model for the redshift-space power spectra of galaxy/halo IAs that can keep the impact of Finger-of-God damping effect, known as a nonlinear systematics of redshift-space distortions, under control. Focusing particularly on galaxy/halo density and IA cross power spectrum, we derive analytically the explicit expressions for the next-to-leading order corrections. Comparing the model predictions with $N$-body simulations, we show that these corrections indeed play an important role for an unbiased determination of the growth-rate parameter, and hence the model proposed here can be used for a precision test of gravity on cosmological scales.

The observed spectra from black hole (BH) X-ray binaries (XRBs) typically consist of two primary components: multitemperature blackbody (BB) originating from the accretion disk in soft X-ray, and a power-law like component in hard X-ray due to Comptonization of soft photons by the hot corona. Illumination of the disk by the corona gives rise to another key component known as reflection. A fraction of the incident hard X-ray radiation is naturally absorbed and re-emitted as a BB at lower energies, referred to as reprocessed BB. For densities relevant to XRBs and typical ionization values, the reprocessed BB may become significant in the soft X-ray region and should be noticeable in the observed spectra as a consequence of reflection. The absence of any BB component in the low/hard state of BH XRB may not be consistent with reflection of high irradiating flux observed as power-law from appropriately dense disk of XRB. We focus on the low/hard state of the BH XRB MAXI J1820+070. We simultaneously fit the shape and flux of the reflection spectra, allowing us to estimate the correct density and ionization of the slab and, correspondingly, the reprocessed BB. Our fitting suggests that the disk in principle may extend close to the BH and still the reprocessed BB due to disk illumination remains small enough to be consistent with the data as opposed to earlier study. The inner reflection component is highly ionized and its fit is primarily driven by its contribution to the continuum. The reprocessed BB cannot resolve whether the disk is extended close to the BH or not in the hard state. For this specific observation, the flux in inner reflection component turns out to be quite low with respect to outer reflection or power-law. Outflowing slab corona covering the inner region of the disk could be the possible geometry of the source with the underlying disk reaching close to the BH. (shortened)

We examine 4 years of Kepler 30-min data, and 5 Sectors of TESS 2-min data for the dM3 star KIC-8507979/TIC-272272592. This rapidly rotating (P=1.2 day) star has previously been identified as flare active, with a possible long-term decline in its flare output. Such slow changes in surface magnetic activity are potential indicators of Solar-like activity cycles, which can yield important information about the structure of the stellar dynamo. We find that while TIC-272272592 shows evidence for both short and long timescale variations in its flare activity, it is unlikely physically motivated. Only a handful of stars have been subjected to such long baseline point-in-time flare studies, and we urge caution in comparing results between telescopes due to differences in bandpass, signal to noise, and cadence. In this work, we develop an approach to measure variations in the flare frequency distributions over time, which is quantified as a function of the observing baseline. For TIC-272272592, we find a $2.7\sigma$ detection of a Sector which has a flare deficit, therefore indicating the short term variation could be a result of sampling statistics. This quantifiable approach to describing flare rate variation is a powerful new method for measuring the months-to-years changes in surface magnetic activity, and provides important constraints on activity cycles and dynamo models for low mass stars.

Stellar magnetic fields have a major impact on space weather around exoplanets orbiting low-mass stars. From an analysis of Zeeman-broadened Fe I lines measured in near-infrared SDSS/APOGEE spectra, mean magnetic fields are determined for a sample of 29 M dwarf stars that host closely orbiting small exoplanets. The calculations employed the radiative transfer code Synmast and MARCS stellar model atmospheres. The sample M dwarfs are found to have measurable mean magnetic fields ranging between $\sim$0.2 to $\sim$1.5 kG, falling in the unsaturated regime on the $<$B$>$ vs P$_{\rm rot}$ plane. The sample systems contain 43 exoplanets, which include 23 from Kepler, nine from K2, and nine from TESS. We evaluated their equilibrium temperatures, insolation, and stellar habitable zones and found that only Kepler-186f and TOI-700d are inside the habitable zones of their stars. Using the derived values of $<$B$>$ for the stars Kepler-186 and TOI-700 we evaluated the minimum planetary magnetic field that would be necessary to shield the exoplanets Kepler-186f and TOI-700d from their host star's winds, considering reference magnetospheres with sizes equal to those of the present-day and young Earth, respectively. Assuming a ratio of 5$\%$ between large-to-small scale B-fields, and a young-Earth magnetosphere, Kepler-186f and TOI-700d would need minimum planetary magnetic fields of, respectively, 0.05 and 0.24 G. These values are considerably smaller than Earth's magnetic field of 0.25 G$\lesssim$B$\lesssim$0.65 G, which suggests that these two exoplanets might have magnetic fields sufficiently strong to protect their atmospheres and surfaces from stellar magnetic fields.

Core accretion is the standard scenario of planet formation, wherein planets are formed by sequential accretion of gas and solids, and is widely used to interpret exoplanet observations. However, no direct probes of the scenario have been discussed yet. Here, we introduce an onion-like model as one idealization of sequential accretion and propose that bulk and atmospheric metallicities of exoplanets can be used as direct probes of the process. Our analytical calculations, coupled with observational data, demonstrate that the trend of observed exoplanets supports the sequential accretion hypothesis. In particular, accretion of planetesimals that are $\gtrsim $ 100 km in size is most favored to consistently explain the observed trends. The importance of opening gaps in both planetesimal and gas disks following planetary growth is also identified. New classification is proposed, wherein most observed planets are classified into two interior statuses: globally mixed and locally (well-)mixed. Explicit identification of the locally (well-)mixed status enables reliable verification of sequential accretion. During the JWST era, the quality and volume of observational data will increase drastically and improve exoplanet characterization. This work provides one key reference of how both the bulk and atmospheric metallicities can be used to constrain gas and solid accretion mechanisms of planets.

Aims. We present MAMBO, a flexible and efficient workflow to build empirical galaxy and Active Galactic Nuclei (AGN) mock catalogues that reproduce the physical and observational properties of these sources. Methods. We start from simulated dark matter (DM) haloes, to preserve the link with the cosmic web, and we populate them with galaxies and AGN using abundance matching techniques. We follow an empirical methodology, using stellar mass functions (SMF), host galaxy AGN mass functions and AGN accretion rate distribution functions studied at different redshifts to assign, among other properties, stellar masses, the fraction of quenched galaxies, or the AGN activity (demography, obscuration, multiwavelength emission, etc.). Results. As a proof test, we apply the method to a Millennium DM lightcone of 3.14 $\rm deg^2$ up to redshift $z=10$ and down to stellar masses $\mathcal{M} \gtrsim 10^{7.5} \, M_\odot$. We show that the AGN population from the mock lightcone here presented reproduces with good accuracy various observables, such as state-of-the-art luminosity functions in the X-ray up to $z \sim 7$ and in the ultraviolet up to $z \sim 5$, optical/NIR colour-colour diagrams, and narrow emission line diagnostic diagrams. Finally, we demonstrate how this catalogue can be used to make useful predictions for large surveys. Using \textit{Euclid} as a case example, we compute, among other forecasts, the expected surface densities of galaxies and AGN detectable in the \textit{Euclid} $H_{\rm E}$ band. We find that \textit{Euclid} might observe (on $H_{\rm E}$ only) about $10^{7}$ and $8 \times 10^{7}$ Type 1 and 2 AGN respectively, and $2 \times 10^{9}$ galaxies at the end of its 14 679 $\rm deg^2$ Wide survey, in good agreement with other published forecasts.