Papers for Monday, Jan 13 2025

Planets orbiting one of the two stars in a binary are vulnerable to gravitational perturbations from the other star. Particularly, highly eccentric companion stars risk disrupting planetary orbits, such as in the extreme system TOI 4633 where close encounters between the companion and a gas giant planet in the habitable zone make it one of the most fragile systems discovered so far. Here, we report that TOI 4633's planet likely survived these encounters throughout the system's age by orbiting retrograde relative to the binary, stabilised by the Coriolis force. Using direct $N$-body simulations, we show it otherwise tends to collide with the binary stars or becomes free-floating after getting ejected. A retrograde planetary orbit has profound implications for TOI 4633's formation and evolution, suggesting an extraordinary history where its eccentric companion was likely randomly captured after planet formation in a single-star system. Alternatively, if stars and planet are born in situ from the same gas clump, we show the planet must have formed at sub-snow-line distances, contrary to the conventional core-accretion model. Our study highlights the importance of considering the long-term stability ($\gtrsim\rm Gyr$) of planets in eccentric binaries and demonstrates that the mere existence in such dynamically hostile environments places strong constraints on their orbital configuration and formation.

We present the discovery of PSR J1947-1120, a new huntsman millisecond pulsar with a red giant companion star in a 10.3 d orbit. This pulsar was found via optical, X-ray, and radio follow-up of the previously unassociated gamma-ray source 4FGL J1947.6-1121. PSR J1947-1120 is the second confirmed pulsar in the huntsman class and establishes this as a bona fide subclass of millisecond pulsar. We use MESA models to show that huntsman pulsars can be naturally explained as neutron star binaries whose secondaries are currently in the "red bump" region of the red giant branch, temporarily underfilling their Roche lobes and hence halting mass transfer. Huntsman pulsars offer a new view of the formation of typical millisecond pulsars, allowing novel constraints on the efficiency of mass transfer and recycling at an intermediate stage in the process.

We present the result from the April 2024 observation of the low-mass X-ray binary GX 13+1 with the Imaging X-ray Polarimetry Explorer (IXPE), together with NICER and Swift-XRT coordinated observations. Two light curve dips were observed; during them, the harder Comptonized spectral component was dominant and the polarization degree higher than in the softer, off-dip intervals. Through a joint analysis of the three IXPE observations, which also included the dip from the first observation, we demonstrate that the polarization properties varied in response to the intensity and spectral hardness changes associated with the dips. The polarization degree attained values up to ~4%. The polarization angle showed a swing of ~70° across the dip and off-dip states, comparable to the continuous rotation seen during the first IXPE observation. We discuss these results in the context of models for polarized emission from the accretion disk and the boundary/spreading layer on the neutron star surface. We also draw attention to the role that an extended accretion disk corona or disk wind can play in generating high polarization degrees and, possibly, swings of the polarization angle.

The extended narrow line region (NLR) of Active Galactic Nuclei (AGN) provides a valuable laboratory for exploring the relationship between AGN and their host galaxies, often appearing as an ''ionization cone'' that can extend out to the very edge of the galaxy. We use new James Webb Space Telescope (JWST) NIRCam imaging to study the morphologies and sizes of NLRs traced by [O III] at cosmic noon ($z\sim2-3$), measurements which were previously only well-studied at $z\sim0$ with IFU spectroscopy. To determine NIRCam's ability to probe the properties of the NLR in cosmic noon galaxies, we present simulated images of AGN at $z=2-3$ created with archival data cubes from the Multi Unit Spectroscopic Explorer (MUSE) of nine nearby ($z<0.05$) AGN host galaxies with previously confirmed extended NLRs. Our simulated images show that NIRCam is able to resolve the morphologies of NLRs at $z=2-3$ with narrow- and medium-band observations. We then search for extended NLRs with NIRCam medium-band observations targeting the [O III] emission in twenty-seven previously identified AGN at $z=2.4-3.4$ in the Great Observatories Origins Deep Survey South (GOODS-S) field. We detect six galaxies in our observed sample with [O III] morphologies consistent with AGN ionization cones with characteristic NLR sizes of $1-2.5$ kpc. Thanks to our simulated data, we can predict the effects of cosmological surface brightness dimming and instrument noise on the NLR size measurements at high redshift, which has the biasing effect of lowering the NLR size-AGN Luminosity trend that is observed at low redshift by a factor of $\sim 2$.

Transient astronomy in the early, high-redshift (z>3) Universe is an unexplored regime that offers the possibility of probing the first stars and the Epoch of Reionization. During Cycles 1 and 2 of the James Webb Space Telescope (JWST), the JWST Advanced Deep Extragalactic Survey (JADES) program enabled one of the first searches for transients in deep images (~ $30$ AB mag) over a relatively wide area ($25$ arcmin$^2$). One transient, AT 2023adsv, was discovered with an F200W magnitude of $28.04$ AB mag, and subsequent JWST observations revealed that the transient is a likely supernova (SN) in a host with $z_{\rm spec}=3.613\pm0.001$, a host mass of log(M$_{*}$/M$_\odot$) = $8.41^{+0.12}_{-0.12}$, and an inferred metallicity at the position of the SN of Z$_{*}$ = $0.3\pm0.1$ Z$_\odot$. At this redshift, the first detections in F115W and F150W show that AT 2023adsv had bright rest-frame ultraviolet flux at the time of discovery. The multi-band light curve of AT 2023adsv is best matched by a template of an SN IIP, with a peak absolute magnitude of $-18.3$ AB mag in the rest-frame $B$-band. We model AT 2023adsv's light curve and find a good match to a $20M_\odot$ red supergiant progenitor star with an explosion energy of $2.0\times 10^{51}$ ergs, likely higher than normally observed in the local Universe, but consistent with SNe IIP drawn from local, lower metallicity environments. AT 2023adsv is the most distant photometrically classified SN IIP yet discovered with a spectroscopic redshift measurement, and may represent a global shift in SNe IIP properties as a function of redshift. This discovery, and the ones sure to follow, demonstrate the continued need for facilities like JWST to build a statistical sample of core-collapse SNe to understand the evolution of their properties, and to constrain the poorly understood relationship between progenitor metallicity and massive star evolution.

We present XRISM Resolve observations of the core of the hot, relaxed galaxy cluster Abell 2029. We find that the line-of-sight bulk velocity of the intracluster medium (ICM) within the central 180 kpc is at rest with respect to the Brightest Cluster Galaxy, with a 3-sigma upper limit of |v_bulk| < 100 km/s. We robustly measure the field-integrated ICM velocity dispersion to be sigma_v = 169 +/- 10 km/s, obtaining similar results for both single-temperature and two-temperature plasma models to account for the cluster cool core. This result, if ascribed to isotropic turbulence, implies a subsonic ICM with Mach number M_3D ~ 0.21 and a non-thermal pressure fraction of 2%. The turbulent velocity is similar to what was measured in the core of the Perseus cluster by Hitomi, but here in a more massive cluster with an ICM temperature of 7 keV, the limit on non-thermal pressure fraction is even more stringent. Our result is consistent with expectations from simulations of relaxed clusters, but it is on the low end of the predicted distribution, indicating that Abell 2029 is an exceptionally relaxed cluster with no significant impacts from either a recent minor merger or AGN activity.

[Abridged] Gas kinematics is a new, unique way to study planet-forming environments by an accurate characterization of disk velocity fields. High angular resolution ALMA observations allow deep kinematical analysis of disks, by observing molecular line emission at high spectral resolution. In particular, rotation curves are key tools to study the disk pressure structure and estimate fundamental disk parameters, as mass and radius. In this work, we explore the potential of a multi-molecular approach to gas kinematics to provide a 2D characterization of the HD 163296 disk. From the high quality data of the MAPS Large Program we extract the rotation curves of rotational lines from seven distinct molecular species, spanning a wide range in the disk radial and vertical extents. To obtain reliable rotation curves for hyperfine lines, we extend standard methodologies to fit multiple-component line profiles. We then sample the likelihood of a thermally stratified model that reproduces all the rotation curves simultaneously, taking into account the molecular emitting layers and disk thermal structure. We obtain dynamical estimates of the stellar mass $M_\star=1.89$ M$_\odot$, the disk mass $M_\text{d}=0.12$ M$_\odot$ and scale radius $ R_\text{c}=143$ au. We also explore how rotation curves and the parameter estimates depend on the adopted emitting layers: the disk mass proves to be the most affected by these systematics, yet the main trends we find do not depend on the adopted parametrization. Finally, we investigate the impact of thermal structure on gas kinematics, showing that the thermal stratification can efficiently explain the measured rotation velocity discrepancies between tracers at different heights. Our results show that such a multi-molecular approach, tracing a large range of emission layers, can provide unique constraints on the ($R,z$) pressure structure of protoplanetary disks.

Using the Square Kilometre Array (SKA) mid precursor MeerKAT, we acquired broadband spectro-polarimetric data in the context of the MeerKAT Fornax Survey to study the Fornax cluster's magnetic fields in detail by building the densest rotation measure (RM) grid to date. Here, we present the survey, the analysis, and a discussion of the RM grid properties. We analyzed a circular region centered on the Fornax cluster center with a radius of $\sim1.4^\circ$; that is, $\rm\sim 0.73 R_{vir}$. The mosaics have a resolution of 13arcsec and cover the frequencies between 900\,MHz and 1.4\,GHz, reaching an average noise of 16$\mu$Jy beam$^{-1}$ in total intensity and 3$\mu$Jy beam$^{-1}$ in the Q and U Stokes images. With these data, we detected 508 polarized sources over an area of $\sim$6.35 deg$^2$ corresponding to a density of $\sim$80 polarized sources/deg$^2$. This is the densest RM grid ever built. Of the polarized sources, five are cluster sources. Excluding the cluster sources, we built the Euclidean-normalized differential source counts in polarization and we went a factor of ten deeper than previous surveys. We tentatively detect for the first time an increment in the differential source counts at low polarized flux densities; that is, $\sim$9\,$\mu$Jy at 1.4\,GHz. The average degree of polarization of about 3--4\% suggests that the sub$-\mu$Jansky population is not dominated by star-forming galaxies, typically showing a degree of polarization lower than 1\%. The majority of the polarized sources are Faraday simple; in other words, their polarization plane rotates linearly with the wavelength squared. The RM shows the typical decrement going from the center to the outskirts of the Fornax cluster. However, interesting features are observed both in the RM grid and in the RM radial profiles across different directions. A combination of the ...

State-of-the-art JWST observations are unveiling unprecedented views into the atmospheres of substellar objects in the infrared, further highlighting the importance of clouds. Current forward models struggle to fit the silicate clouds absorption feature at ~$10\,\mu$m observed in substellar atmospheres. In the MSG model, we aim to couple the MARCS 1D radiative-convective equilibrium atmosphere model with the 1D kinetic, stationary, non-equilibrium, cloud formation model DRIFT, to create a new grid of self-consistent cloudy substellar atmosphere models with microphysical cloud formation. We aim to test if this new grid is able to reproduce the silicate cloud absorption feature at ~$10\,\mu$m. We model substellar atmospheres with effective temperatures in the range 1200-2500 K and with $\log(g)=4.0$. We compute atmospheric structures that self-consistently account for condensate cloud opacities based on microphysical properties. We present an algorithm based on control theory to help converge such self-consistent models. Synthetic atmosphere spectra are computed for each model to explore the observable impact of the cloud microphysics. We additionally explore the impact of choosing different nucleation species (TiO$_2$ or SiO) and the effect of less efficient atmospheric mixing on these spectra. The new MSG cloudy grid using TiO$_2$ nucleation shows spectra which are redder in the near-infrared compared to the currently known population of substellar atmospheres. We find the models with SiO nucleation, and models with reduced mixing efficiency are less red in the near-infrared. The grid is unable to reproduce the silicate features similar to those found in recent JWST observations and Spitzer archival data. We thoroughly discuss further work that may better approximate the impact of convection in cloud-forming regions and steps that may help resolve the silicate cloud feature.

Upon their evaporation via Hawking radiation, primordial black holes (PBHs) may deposit energy in the ambient plasma on scales smaller than the typical distance between two black holes, leading to the formation of hot spots around them. We investigate how the corresponding rise of the local temperature during the evaporation may act as a shield against the release of low-energy photons, affecting PBH's capacity to dissociate light nuclei after Big-Bang Nucleosynthesis through photo-dissociation. We study the different ways PBH hot spots affect the flux of low-energy photons expected from PBH evaporation, and we find that such effects can be particularly relevant to the physics of photo-dissociation during Big-Bang Nucleosynthesis for PBHs with masses between $10^{11}$g and $3\times 10^{12}$g. We emphasize that the magnitude of this effect is highly dependent on the specific shape of the temperature profile around PBHs and its time evolution. This underscores the necessity for a comprehensive study of PBH hot spots and their dynamics in the future.

Space-based experiments, either orbiting the Earth or from scientific balloon altitudes, measure high-energy cosmic rays by measuring from above the atmosphere the optical and radio signals generated by extensive air showers (EAS). These experiments are designed to have a large field-of-view (FoV) for observing EAS which translates to monitoring the atmosphere over a large, $\sim 10^6$ km$^2$ area on the ground. Ultra-high energy cosmic rays (UHECRs, $E_{CR} > ~\sim 1$ EeV) are measured by using the isotropic near-UV air fluorescence signal to finely sample the EAS development and to efficiently use the atmosphere as a vast calorimeter. At UHE, these immense EAS particle cascades have sufficient charged particle content to generate the relatively dim fluorescence light that propagates to the space-based instrument. Additionally, the beamed Cherenkov light and geomagnetic radio emission from EAS arrive with $> ~\sim 10$ ns impulse and are measured at small angles away from the cosmic ray trajectories. In particular for optical Cherenkov measurements, the energy thresholds can be $> ~\sim 1$ PeV, i.e. for very-high energy cosmic rays (VHECRs). ... In this chapter, the nature of observing the UHECR-induced shower development from orbiting and balloon-borne experiments is detailed, both for missions that have been flown and those currently in development. This will be accomplished by discussing experimental performance in terms of measuring the UHECR spectrum, UHECR nuclear composition, and UHECR arrival direction. The ability for these experiments to also perform VHECR and VHE and UHE cosmic neutrino measurements will also be discussed.

We present an analysis of four Chandra observations of the 45 Myr old DS Tuc binary system. We observed X-ray variability of both stars on timescales from hours to months, including two strong X-ray flares from star A. The implied flaring rates are in agreement with past observations made with XMM-Newton, though these rates remain imprecise due to the relatively short total observation time. We find a clear, monotonic decline in the quiescent level of the star by a factor 1.8 across eight months, suggesting stellar variability that might be due to an activity cycle. If proven through future observations, DS Tuc A would be the youngest star for which a coronal activity cycle has been confirmed. The variation in our flux measurements across the four visits is also consistent with the scatter in empirical stellar X-ray relationships with Rossby number. In simulations of the possible evolution of the currently super-Neptune-sized planet DS Tuc Ab, we find a range of scenarios for the planet once it reaches a typical field age of 5 Gyr, from Neptune-size down to a completely stripped super-Earth. Improved constraints on the planet's mass in the future would significantly narrow these possibilities. We advocate for further Chandra observations to better constrain the variability of this important system.

The reionization of helium is thought to occur at $2.5\lesssim z\lesssim4$, marking the last phase transition and final global heating event of the intergalactic medium (IGM). Since it is driven by rare quasars, helium reionization should give rise to strong temperature fluctuations in the IGM between neutral and recently-ionized regions of order $\sigma (\ln T) \sim \Delta T/T = 20-50\%$. We introduce a novel method to search for reionization-induced temperature fluctuations in the IGM by using the effective optical depths of the Lyman-$\alpha$ forest towards a large number of background quasars. Higher IGM temperatures give rise to lower effective optical depths in the Lyman-$\alpha$ forest, implying that temperature fluctuations will broaden the observed optical depth distribution. We measured the distributions of effective Lyman-$\alpha$ forest optical depths across $71$ X-Shooter spectra from the XQ-100 survey in four redshift bins from $z=3.76$ to $z=4.19$ and compared them to a large-volume cosmological hydrodynamical simulation. A good agreement is found between the observations and the simulation, which does not include temperature fluctuations; therefore, we do not detect a signature of helium reionization. We then post-process the simulations to include an increasing amount of temperature fluctuations until the model becomes inconsistent with the observations. We obtain tight constraints on $\sigma (\ln T) < 0.29 \ (<0.40)$ at $2 \sigma\ (3 \sigma)$ at $z=3.76$ when averaging over scales of $100$ comoving Mpc, and weaker constraints for higher redshifts and smaller scales. Our constraints are the tightest to date, and imply that either the IGM temperature contrast caused by helium reionization is less than $\sim30\%$, or that the process has not yet significantly started at $z=3.76$.

We report the analysis of ~1Ms of XMM-Newton observations of the rapidly accreting active galactic nucleus REJ1034+396. The 0.3-9 keV EPIC-pn spectra are well described by a model consisting of steep continuum emission from the corona accompanied by relativistically-blurred reflection from a highly ionized accretion disk. The source is known to exhibit strong excess soft X-ray emission, which we show is well represented by thermal disk photons Comptonized by a warm plasma spanning the inner accretion flow. Additionally, the EPIC-pn data provide compelling evidence ($\Delta C\sim60$ for 4 additional parameters) for the presence of an ultrafast outflow (UFO) with a line-of-sight velocity $v/c=0.307^{+0.001}_{-0.005}$, and an emission signature consistent with reflection of the corona from modestly ionized, outflowing gas. The simultaneous $0.5-2.5$\,keV RGS spectra show clear absorption lines. Modelling of these data confirms the presence of the UFO and constrains its equivalent hydrogen column density, log $N_\mathrm{H}/$(atom cm$^{-2}$) = $21.7^{+0.1}_{-0.2}$. The RGS data also reveal at least two warm absorber components with a modest outflow velocity ($1680^{+40}_{-50}$ km/s). The measured properties and time evolution of the UFO in REJ1034+396 suggest that it is formed from collisionally ionized plasma, launched from the disk surface and accelerated by radiation pressure. The high terminal velocity and substantial absorbing column density imply that the outflow carries sufficient momentum and energy to transform its environment, being capable of driving out essentially all dust and gas it interacts with along the line of sight, even if the AGN were initially surrounded by a Compton thick absorber.

The interstellar medium of galaxies is filled with holes, bubbles, and shells, typically interpreted as remnants of stellar evolution. There is growing interest in the study of their properties to investigate stellar and supernova feedback. So far, the detection of cavities in observational and numerical data is mostly done visually and, hence, is prone to biases. Therefore, we present an automated, objective method for discovering cavities in particle simulations, with demonstrations using hydrodynamical simulations of a dwarf galaxy. The suggested technique extracts holes based on the persistent homology of particle positions and identifies tight boundary points around each. With a synthetic ground-truth analysis, we investigate the relationship between data density and the detection radius, demonstrating that higher data density also allows for the robust detection of smaller cavities. By tracking the boundary points, we can measure the shape and physical properties of the cavity, such as its temperature. In this contribution, we detect 808 holes in 21 simulation snapshots. We classified the holes into supernova-blown bubbles and cavities unrelated to stellar feedback activity based on their temperature profile and expansion behaviour during the 100 million years covered by the simulation snapshots analysed for this work. Surprisingly, less than 40% of the detected cavities can unequivocally be linked to stellar evolution. Moreover, about 36% of the cavities are contracting, while 59% are expanding. The rest do not change for a few million years. Clearly, it is erroneous to interpret observational data based on the premise that all cavities are supernova-related and expanding. This study reveals that supernova-driven bubbles typically exhibit smaller diameters, larger expansion velocities, and lower kinetic ages (with a maximum of 220 million years) compared to other cavities.

The detection of periodic signals in irregularly-sampled time series is a problem commonly encountered in astronomy. Traditional tools used for periodic searches, such as the periodogram, have poorly defined statistical properties under irregular sampling, which complicate inferring the underlying aperiodic variability used for hypothesis testing. The problem is exacerbated in the presence of stochastic variability, which can be easily mistaken by genuine periodic behaviour, particularly in the case of poorly sampled lightcurves. Here we present a method based on Gaussian Processes (GPs) modelling for period searches and characterization, specifically developed to overcome these problems. We argue that in cases of irregularly-sampled time series, GPs offer an appealing alternative to traditional periodograms, because the known distribution of the data (correlated Gaussian) allows a well-defined likelihood to be constructed. We exploit this property and draw from existing statistical methods to perform traditional likelihood ratio tests for an additional, (quasi-)periodic component, using the aperiodic variability inferred from the data as the null hypothesis. Inferring the noise from the data allows the method to be fully generalizable, with the only condition that the data can be described as a Gaussian process. We demonstrate the method by applying it to a variety of objects showing varying levels of noise and data quality. Limitations of the method are discussed and a package implementing the proposed methodology is made publicly available.

This paper aims to describe the Pointing Reconstruction Model (PRM) and the prototype Star Tracker, which will be mounted on LSPE-Strip, a microwave Q- and W-band CMB telescope planned for installation at the "Observatorio del Teide" in Tenerife. The PRM integrates information on the instantaneous attitude provided by the telescope control system to determine the actual pointing direction and focal plane orientation of the telescope. It accounts for various non-idealities in the telescope setup, represented by eight configuration angles, which will be calibrated using the Star Tracker. Following the derivation of the PRM formalism and its implementation, we investigate the pointing errors caused by incorrect calibration of these configuration angles to validate the required 1 arcminute maximum systematic pointing error for the LSPE-Strip survey. This paper also describes the main structure and operations of the Star Tracker and presents the results of a campaign of actual sky observations conducted with a prototype. The results demonstrate a Star Tracker RMS accuracy of approximately 3 arcseconds, while systematic errors remain below 10 arcseconds. Based on these results, we analyzed the problem of reconstructing the PRM configuration angles. Two methods for intercalibrating the Star Tracker's pointing direction with respect to the focal plane's pointing direction were examined: (1) observations of planets and (2) observations of a drone carrying both an optical beacon and a radio beacon. In the first case, an intercalibration accuracy between 1/3 arcminute and 1 arcminute is achievable. In the second case, the expected intercalibration accuracy ranges from 0.25 arcminute to 1 arcminute.

Irradiated brown dwarfs offer a unique opportunity to bridge the gap between stellar and planetary atmospheres. We present high-quality $\mathit{HST}$/WFC3/G141 phase-resolved spectra of the white dwarf + brown dwarf binary GD 1400, covering more than one full rotation of the brown dwarf. Accounting for brightness variations caused by ZZ Ceti pulsations, we revealed weak ($\sim$1\%) phase curve amplitude modulations originating from the brown dwarf. Sub-band light curve exploration in various bands showed no significant wavelength dependence on amplitude or phase shift. Extracted day- and night-side spectra indicated chemically similar hemispheres, with slightly higher day-side temperatures, suggesting efficient heat redistribution or dominance of radiative escape over atmospheric circulation. A simple radiative and energy redistribution model reproduced observed temperatures well. Cloud-inclusive models fit the day and night spectra better than cloudless models, indicating global cloud coverage. We also begin qualitatively exploring atmospheric trends across six irradiated brown dwarfs, from the now complete "Dancing with Dwarfs" WD$-$BD sample. The trend we find in the day-side/night-side temperature and irradiation levels is consistent with efficient heat redistribution for irradiation levels less than $\sim$10$^9$ ergs/s/cm$^2$ and decreasing efficiency above that level.

The Rossiter-McLaughlin effect allows us to measure the projected stellar obliquity of exoplanets. From the spin-orbit alignment, planet formation and migration theories can be tested to improve our understanding of the currently observed exoplanetary population. Despite having the spin-orbit measurements for more than 200 planets, the stellar obliquity distribution is still not fully understood, warranting additional measurements to sample the full parameter space. We analyze archival HARPS and HARPS-N spectroscopic transit time series of eight gas giant exoplanets on short orbits and derive their projected stellar obliquity $\lambda$. We report a prograde, but misaligned orbit for HAT-P-50b ($\lambda =41^\circ\ ^{+10}_{-9}$), possibly hinting at previous high-eccentricity migration given the presence of a close stellar companion. We measured sky-projected obliquities that are consistent with aligned orbits for the rest of the planets: WASP-48b ($\lambda =-4^\circ\ \pm$ 4), WASP-59b ($\lambda =-1^\circ\ ^{+20}_{-21}$), WASP-140 Ab ($\lambda =-1^\circ\ \pm$ 3), WASP-173 Ab ($\lambda =9^\circ\ \pm$ 5), TOI-2046b ($\lambda =1^\circ\ \pm$ 6), HAT-P-41 Ab ($\lambda =-4^\circ\ ^{+5}_{-6}$), and Qatar-4b ($\lambda =-13^\circ\ ^{+15}_{-19}$). We measure the true stellar obliquity $\psi$ for four systems. We infer a prograde, but misaligned, orbit for TOI-2046b with $\psi=$42$^{+10}_{-8}$\,deg. Additionally, $\psi = 30^\circ\ ^{+18}_{-15}$ for WASP-140 Ab, $\psi = 21^\circ\ ^{+9}_{-10}$ for WASP-173 Ab, and $\psi = 32^\circ\ ^{+14}_{-13}$ for Qatar-4b. The aligned orbits are consistent with slow disk migration, ruling out violent events that would excite the orbits over the history of these systems. Finally, we provide a new age estimate for TOI-2046 of at least 700 Myr and for Qatar-4 of at least 350-500 Myr, contradicting previous results.

We report our findings on a spectroscopic survey of seven unresolved DA+DB binary white dwarf candidates. We have discovered extreme spectroscopic variations in one of these candidates, SDSS J084716.21+484220.40. Previous analysis failed to reproduce the optical spectrum using a single object with a homogeneous atmosphere. Our time-resolved spectroscopy reveals a double-faced white dwarf that switches between a DBA and DA spectral type over 6.5 or 8.9 hours due to varying surface abundances. We also provide time-series spectroscopy of the magnetic DBA, SDSS J085618.94+161103.6 (LB 8915), and confirm an inhomogeneous atmosphere. We employ an atmosphere model with hydrogen caps and a helium belt that yields excellent fits to our time-resolved spectra. We use the oblique rotator model to derive the system geometry for both targets. With the addition of these two objects, the emerging class of double-faced white dwarfs now consists of seven members. We summarize the properties of this new class of objects, and discuss how magnetism impacts the convective processes and leads to the formation of double-faced white dwarfs. We identify cooler versions of white dwarfs with inhomogeneous atmospheres among the cool magnetic DA white dwarf sample, where the H$\alpha$ line is shallower than expected based on pure hydrogen atmosphere models.

The short-term X-ray variability of Cyg X-1 can be interpreted as random occurrence of mini-flares known as the shots, whose physical nature is still unclear. We propose a new algorithm for shot identification in the X-ray light curve, based on baseline detection and template fitting. Compared with previous techniques, our algorithm allows us to detect shots with lower amplitudes and shorter time separations. With NICER observations, we find that, after correction for detection sensitivity, both the shot amplitude and recurrence rate are positively scaled with the mean count rate, while the recurrence rate has a much higher dependence on the count rate. These suggest that a higher mass accretion rate will drive more and slightly larger shots. We also find that the abrupt hardening near the shot peak found in previous studies is attributed to different shot profiles in different energy bands; there is no need to involve a rapid physical process to suddenly harden the emitting spectrum.

Very long baseline interferometry (VLBI) enables high-angular-resolution observations in astronomy and geodesy by synthesizing a virtual telescope with baselines spanning hundreds to thousands of kilometres. Achieving high instrumental phase stability in VLBI relies on the generation of high-quality, atomic-referenced RF local oscillator (LO) and RF-comb signals for the effective downconversion of celestial RF signals and precise phase calibration, respectively. As observing frequencies move into higher ranges with wider bandwidth, conventional electronic methods face significant challenges in maintaining the quality of these signals. Here, we demonstrate that an optical frequency comb (OFC) can be used as a versatile tool to generate and distribute low-noise and atomic-referenced RF-comb and RF-LO signals in the VLBI telescope. Hydrogen maser-stabilized optical pulses are transmitted over a timing-stabilized fibre link from the observatory building to the VLBI receiver system at the telescope, where photodetection converts them into the required RF-comb and RF-LO signals. In VLBI test observation, we successfully detected VLBI fringes and extracted the RF-combs characteristics in a format suitable for VLBI instrumental phase calibration. These results highlight the high potential of OFC-based technology for enhancing next-generation broadband VLBI measurements, advancing astrophysical research and facilitating intercontinental clock comparison.

The intersection of art and science offers novel ways to interpret and represent complex phenomena. This project explores the convergence of high-energy astrophysics, concrete poetry, and natural language processing (NLP) by proposing a Markov chain-based algorithm that generates poetic texts from gamma-ray emission maps of the universe. Gamma rays, the most energetic form of electromagnetic radiation, arise from violent astrophysical processes such as supernovae, pulsars, and black hole accretion, observable through instruments like the \textit{Fermi Large Area Telescope} (FermiLAT). These high-energy events are mapped and processed into matrices, treated as Markov chains, where each state's transition probability is determined by the intensity of gamma-ray sources in a region of interest.

The ExoEcho project is designed to study the photodynamics of exoplanets by leveraging high-precision transit timing data from ground- and space-based telescopes. Some exoplanets are experiencing orbital decay, and transit timing variation (TTV) is a useful technique to study their orbital period variations. In this study, we have obtained transit middle-time data from the Hubble Space Telescope (HST) observations for 37 short-period exoplanets, most of which are hot Jupiters. To search for potential long- and short-term orbital period variations within the sample, we conduct TTV model fitting using both linear and quadratic ephemeris models. Our analysis identifies two hot Jupiters experiencing strong periodic decays. Given the old age of the host stars of the hot Jupiter population, our findings call for a scenario where HJs are continuously being destructed and created. Our study demonstrates the importance of incorporating high-precision transit timing data to TTV study in the future.

Despite the broad successes of the flat $\Lambda$CDM model and its fitness to the various cosmological observations, it confronts challenges stemming from anomalies in the measurements of the Hubble constant ($H_0$) and the amplitude of matter fluctuations ($\sigma_8$). These inconsistencies have necessitated a reassessment of the model parameters, with a particular focus on their potential dependence on redshift. This study pioneers a new investigation to probe this redshift dependency by generating mock data simulated from observational data of Type Ia supernovae (SNIa) and cosmic chronometers (CC), thereby increasing the data density in this field. By sorting the data into high-redshift and low-redshift bins, we aim to refine the cosmological constraints on the parameters of the $\Lambda$CDM model and determine whether the noted dependence on redshift is due to a lack of high-redshift observational data or if they signify intrinsic issues within the model itself. Our approach employs the Markov Chain Monte Carlo (MCMC) algorithm to minimize the $\chi^2$ function, thus tightening the cosmological constraints. Our findings within the mock analysis reveal discrepancies between the values of $\Omega_{m0}$ and $H_0$ derived from the mock data bins with high redshift and low redshift, indicating the potential deviation of the standard $\Lambda$ CDM cosmology from the high-redshift SNIa and CC data. If this deviation proposes a new physics beyond the standard model, then with better quality future data tracking the new physics, these discrepancies will be statistically significant.

We performed a detailed spectroscopic analysis of three extremely metal-poor RR Lyrae stars, exploring uncharted territories at these low metallicities for this class of stars. Using high-resolution spectra acquired with HARPS-N at TNG, UVES at VLT, and PEPSI at LBT, and employing Non-Local Thermodynamic Equilibrium (NLTE) spectral synthesis calculations, we provide abundance measurements for Fe, Al, Mg, Ca, Ti, Mn, and Sr. Our findings indicate that the stars have metallicities of [Fe/H] = -3.40 \pm 0.05, -3.28 \pm 0.02, and -2.77 \pm 0.05 for HD 331986, DO Hya, and BPS CS 30317-056, respectively. Additionally, we derived their kinematic and dynamical properties to gain insights into their origins. Interestingly, the kinematics of one star (HD 331986) is consistent with the Galactic disc, while the others exhibit Galactic halo kinematics, albeit with distinct chemical signatures. We compared the [Al/Fe] and [Mg/Mn] ratios of the current targets with recent literature estimates to determine whether these stars were either accreted or formed in situ, finding that the adopted chemical diagnostics are ineffective at low metallicities ([Fe/H] $\lesssim -$1.5). Finally, the established horizontal branch evolutionary models, indicating that these stars arrive at hotter temperatures on the Zero-Age Horizontal Branch (ZAHB) and then transition into RR Lyrae stars as they evolve, fully support the existence of such low-metallicity RR Lyrae stars. As a consequence, we can anticipate detecting more of them when larger samples of spectra become available from upcoming extensive observational campaigns.

We present the statistical properties of Ly$\alpha$ emission in 586 galaxies at $z=4.5-14.2$, observed by multiple JWST/NIRSpec spectroscopy projects, including JADES, GLASS, CEERS, and GO/DDT programs. We obtain Ly$\alpha$ equivalent width (EW), Ly$\alpha$ escape fraction, and ionizing photon production efficiency measurements or upper limits for these galaxies, and confirm that the Ly$\alpha$ emitting galaxy fraction decreases towards higher redshifts. We derive Ly$\alpha$ luminosity functions from $z\sim 5$ to $z\sim 10-14$ with the observed Ly$\alpha$ EW distributions and galaxy UV luminosity functions, and find a $\sim3$ dex decrease in number density at $L_\mathrm{Ly\alpha}=10^{42}-10^{43}$ erg s$^{-1}$ over the redshift range. We obtain the neutral hydrogen fractions of $x_\mathrm{HI}=0.17_{-0.16}^{+0.23}$, $0.63_{-0.28}^{+0.18}$, $0.79_{-0.21}^{+0.13}$, and $0.88_{-0.13}^{+0.11}$ at $z\sim6$, $7$, $8-9$, and $10-14$, respectively, via comparisons of the reionization models developed by semi-numerical simulations with 21cmFAST explaining the observations of Ly$\alpha$, UV continuum, and Planck electron optical depth. The high $x_\mathrm{HI}$ values over $z\sim 7-14$ suggest a late and sharp reionization, with the primary reionization process occurring at $z\sim 6-7$. Such a late and sharp reionization is not easily explained by either a clumpy inter-galactic medium or sources of reionization in a classical faint-galaxy or a bright-galaxy/AGN scenario, unless a very high escape fraction or AGN duty cycle is assumed at $z\sim 6-7$.

The accretion of material onto a black hole modifies the properties of that hole owing to the capture of both matter and radiation. Adding matter to the hole through an accretion disc generally acts to increase the hole's spin parameter, while the capture of radiation generally provides a retarding torque. The balance between the torques provided by adding matter and radiation leads to a maximum spin parameter that can be obtained by a black hole which grows by accretion, known as the Thorne limit. In the simplest theory of thin disc accretion this Thorne limit has the value $a_{\bullet, {\rm lim}} \simeq 0.998$. The purpose of this paper is to highlight that any modification to theories of accretion flows also modify this limiting value, and to compute precisely the modification arising from a particular extension of accretion theory: the inclusion of bright emission from within the plunging region which is sourced from the magnetohydrodynamic stresses ubiquitously observed in simulations. This extra emission further suppresses black hole spin-up and results in new, lower, limits on the final black hole spin. These limits depend on the details of the magnetic stresses acting within the plunging region, but typical values seen in simulations and observations would lower the limit to $a_{\bullet, {\rm lim}} \simeq 0.99$, a subtle but not negligible deviation.

WASP-49Ab, a low-density, Saturn-like planet in a tight orbit around a Sun-like star within a wide binary system, is a compelling candidate for hosting a volcanic moon, as suggested by the detection of Doppler-shifted this http URL study evaluates the stability of potential satellites around WASP-49Ab under the influence of planetary oblateness, relativistic effects, and perturbations from a close companion star, focusing on their impact on light curve parameters such as transit duration and impact parameter variations, driven by the evolution of the planet's orbit in this extreme environment. Using N-body simulations and semi-analytical methods, we analysed moon's dynamics across varied initial conditions and gravitational frameworks including the potential of an oblate planet and the effects of the general relativity. We find that `selenity', a moon survival indicator, is high in close orbits with low eccentricity, near the Roche limit, especially for masses greater than Io's. Stability decreases as eccentricity or distance from the planet increases. Additionally, We observe a strong destabilising resonance near 1.4 planetary radii when planetary eccentricities are slightly greater than zero. This study confirms the potential for stable exomoons around WASP-49Ab despite its hostile environment, emphasizing the importance of incorporating diverse physical effects in stability analyses, aiding future detection efforts.

Infrared spectroscopy, e.g., with JWST, provides a glimpse into the chemical inventory of the innermost region of protoplanetary discs, where terrestrial planets eventually form. The chemical make-up of regions inside snowlines is connected to the material drifting from the outer regions, which can be modeled with dust evolution models. However, infrared observations are limited by the high dust extinction in the inner disc, and only probes the abundances of gaseous species in the disc surface layers. As a result, the bulk mass of delivered volatiles is not directly relatable to what is measured through infrared spectra. In this paper, we investigate how the delivery of dust and ice after prolonged pebble drift affects the observable reservoir of water vapor in the inner disc. We develop a 1+1D approach based on dust evolution models to determine the delivery and distribution of vapor compared to the height of the $\tau = 1$ surface in the dust continuum. We find that the observable column density of water vapor at wavelengths probed by JWST spans many orders of magnitude over time, exhibiting different radial profiles depending on dust properties, drift rate, and local processing. In the presence of a traffic-jam effect inside the snowline, the observable vapor reservoir appears constant in time despite the ongoing delivery by pebble drift, such that water is effectively smuggled unnoticed. Differences in measured column densities then originate not only from variations in bulk vapor content, but also from differences in the properties and distribution of dust particles.

We investigate the spectral characteristics of the Be/X-ray binary system, EXO 051910+3737.7, in which Be/X-ray systems are the largest sub-class of high-mass X-ray binaries. Spectroscopic observations are taken by the Thai National Telescope (TNT) with a Medium-RESolution spectrograph (MRES) instrument for seven nights spanning from 2020 to 2021. Our primary focus is directed towards the analysis of two Balmer lines, namely H$\alpha$ and H$\beta$, given that Be stars typically exhibit emission features in at least one of these hydrogen Balmer lines during certain phases. Our observations reveal split Balmer emission lines throughout the entire duration of our monitoring. Double Gaussian profiles were employed for line fitting to characterize these lines. The presence of double peaks in the Balmer lines indicates the presence of asymmetries within the circumstellar disc. We then analyze V/R variations and the changes in H$\beta$ spectra. Our analysis of V/R variation which Violet (V) and Red (R) peak intensity components, revealed rapid fluctuations occurring within a single day, although determining the precise periodicity was constrained by instrumental limitations and the duration of observability. Furthermore, employing observed the wavelength differences ($\Delta\lambda$) in conjunction with typical Be star parameters allowed us to estimate the radius ($r_{\beta}$) of the H$\beta$ emitting envelope. The average value was calculated to be 2.585$r_*$, with a standard deviation of 0.050$r_*$.

A common feature of Active Galactic Nuclei (AGN) is their random variations in brightness across the whole emission spectrum, from radio to $\gamma$-rays. Studying the nature and origin of these fluctuations is critical to characterising the underlying variability process of the accretion flow that powers AGN. Random timing fluctuations are often studied with the power spectrum; this quantifies how the amplitude of variations is distributed over temporal frequencies. Red noise variability -- when the power spectrum increases smoothly towards low frequencies -- is ubiquitous in AGN. The commonly used Fourier analysis methods, have significant challenges when applied to arbitrarily sampled light curves of red noise variability. Several time-domain methods exist to infer the power spectral shape in the case of irregular sampling but they suffer from biases which can be difficult to mitigate, or are computationally expensive. In this paper, we demonstrate a method infer the shape of broad-band power spectra for irregular time series, using a Gaussian process regression method scalable to large datasets. The power spectrum is modelled as a power-law model with one or two bends with flexible slopes. The method is fully Bayesian and we demonstrate its utility using simulated light curves. Finally, Ark 564, a well-known variable Seyfert 1 galaxy, is used as a test case and we find consistent results with the literature using independent X-ray data from XMM-Newton and Swift. We provide publicly available, documented and tested implementations in Python and Julia.

Dust grains in protoplanetary disks are the building blocks of planets. Investigating the dust composition and size, and their variation over time, is crucial for understanding the planet formation process. The PDS 70 disk is so far the only protoplanetary disk with concrete evidence for the presence of young planets. Mid-infrared spectra were obtained for PDS 70 by the Infrared Spectrograph (IRS) on the Spitzer Space Telescope (SST) and the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2007 and 2022, respectively. In this work, we investigate the dust mineralogy through a detailed decomposition of the observed mid-infrared spectra. The results show that both the dust size and crystallinity increased by a factor of about two during the two epochs of observation, indicating evident dust processing in the terrestrial planet-forming region of the PDS 70 disk. The dust size (~0.8 micron) and crystallinity (~6%) in the PDS 70 disk are similar to those of other disks, which implies that the two nascent planets, PDS 70b and PDS 70c located at radial distances of ~22AU and ~34AU, do not have a significant impact on the dust processing in the inner disk. The flux densities at wavelengths longer than ~16 micron measured by JWST/MIRI are only ~60% of those obtained by Spitzer/IRS. Based on self-consistent radiative transfer modeling, we found that such a strong variability in mid-infrared fluxes can be produced by adjustments to the dust density distribution and structure of the inner disk probably induced by planet-disk interaction.

We present axially symmetric analytical potential-density pairs with surface density similar to the Miyamoto-Nagai model, but with more realistic vertical structure. Our models closely approximate an exponential, a sech$^2$, or a cored exponential vertical density profile. The latter profile has a density core of adjustable width, which provides more flexibility when modelling galaxy discs.

Multiple studies on radial migration in disc galaxies have proven the importance of the effect of resonances with non-axisymmetric components on the evolution of galactic discs. However, the dynamical effects of classic Newtonian dynamics with dark matter (DM) differ from MOdified Newtonian Dynamics (MOND) and might trigger different radial migration. A thorough analysis of radial migration considering these two gravitational regimes might shed some light on different predictions of DM and MOND theories. We aim to quantitatively and qualitatively compare the effects of resonances and stellar radial migration (churning) in a Milky Way-like (MW-like) galaxy in the DM and MOND regimes. We performed simulations of a MW-like galaxy to analyse the effect of non-axisymmetric structures (galactic bar and spiral arms) considering various parameters of the spiral structure. We conducted a two-dimensional numerical simulation consisting of the integration of $2 \cdot 10^6$ stars in a static rotating galactic potential for $6~\mbox{Gyr}$. We analysed the change in the star's position, the guiding radius, as well as the frequency phase space. We investigated DM and MOND approaches. The outcome of the simulation shows that the radial migration is much more pronounced in the MOND regime compared to the DM one. Compared to the DM approach, in the MOND regime, we observe up to five times as many stars with a maximum change in the guiding radius of more than $1.5~\mbox{kpc}$ during the time interval from $2-6~\mbox{Gyr}$.Analysis of the frequency phase space reveals that the most prominent resonances in all DM and MOND configurations are the co-rotation resonance with the spiral arms ($m=p=1$), outer Lindblad resonance with the galactic bar and spiral arms, and the co-rotation resonance ($m=2$, $p=1$) with the superposition of the galactic bar and spiral arms, $2 \Omega = \Omega_b + \Omega_{sp}$.

The CIELO project introduces a novel set of chemo-dynamical zoom-in simulations designed to simultaneously resolve galaxies and their nearby environments. The initial conditions include a diverse range of cosmic structures, such as local groups, filaments, voids, and walls, allowing for a detailed exploration of galaxies within the broader context of the cosmic web. This study presents the initial conditions and characterizes the global properties of CIELO galaxies and their environments. It focuses on galaxies with stellar masses ranging from log [8,11] solar masses and examines key scaling relations, including the mass-size relation, the Tully-Fisher relation, and the mass-metallicity relation for both stars and star-forming gas. The DisPerSe algorithm was used to determine the positions of CIELO galaxies within the cosmic web, with a specific focus on the Pehuen haloes. The selection of local group volumes was guided by criteria based on the relative positions and velocities of the two primary galaxies. The Pehuen regions were chosen to map walls, filaments, and voids. Synthetic images in the SDSS i, r, and g bands were generated using the SKIRT radiative transfer code. Additionally, a dynamical decomposition was performed to classify galaxy morphologies into bulge, disc, and stellar halo components (abridged).

This paper discusses a fast algorithm for analyzing the phases of a multi-tone phase calibration signal. The multi-tone method is superior to its single-tone counterpart. The PCal signal is a wide-band frequency comb derived from an ultra-stable frequency standard which is then coupled into the input stage of the radio astronomy receiver chain. The frequency comb allows to determine the signal delay, frequency response and phase shifts induced by the radio telescope instrumentation. The instrumentation delay is one of the several delay components that have to be compensated for in multi-telescope aperture synthesis. The frequency response also helps to monitor the instrumentation. To extract the information from the digitized baseband output of the receiver chain calls for fast numerical methods. A simple fast algorithm is presented in this paper. Two implementations are provided. One is based on precomputed look-up tables. The other is based on fast digital sinewave generators. Both lend themselves to parallel software processing on multi-core and SIMD-capable platforms as well as processing with FPGA gateware. We present optimized reference implementations for Intel/AMD, IBM Cell and nVidia GPU platforms.

The spectrum of cosmic-ray electrons depends sensitively on the history and spatial distribution of nearby sources. Given our limited observational handle on cosmic-ray sources, any model remains necessarily probabilistic. Previously, predictions were performed in a Monte Carlo fashion, summing the contributions from individual, simulated sources to generate samples from the statistical ensemble of possible electron spectra. Such simulations need to be re-run if the cosmic-ray transport parameters (e.g. diffusion coefficient, maximum energy) are changed, rendering any parameter study computationally expensive. In addition, a proper statistical analysis of observations and comparison with such probabilistic models requires the joint probability distribution of the full spectrum instead of only samples. Note that parametrising this joint distribution is rendered difficult by the non-Gaussian statistics of the cosmic-ray fluxes. Here, we employ machine learning to compute the joint probability distribution of cosmic-ray electron fluxes. Specifically, we employ masked autoregressive density estimation (MADE) for a representation of the high-dimensional joint probability distribution. In a first step, we train the network on a Monte Carlo simulation for a fixed set of transport parameters, thus significantly accelerating the generation of samples. In a second step, we extend this setup to SECRET (Stochasticity Emulator for Cosmic Ray Electrons), allowing to reliably interpolate over the space of transport parameters. We make the MADE and SECRET codes available at this https URL .

We report the X-ray spectral and timing analysis of the high mass X-ray binary EXO 2030+375 during the 2021 type-II outburst based on the Insight-HXMT observations. Pulsations can be detected in the energy band of 1-150 keV. The pulse profile shows energy and luminosity dependence and variability. We observed transitions in the pulse profile shape during the rising and the decaying phase of the outburst. The pulse fraction exhibits an anti-correlation with luminosity and a non-monotonic energy dependence, with a possible dip near 30 keV during the outburst peak. The hardness-intensity diagrams (7-10 keV/4-7 keV) suggest state transitions during the early and late phases of the outburst. These transitions are consistent with the luminosity at which the pulse profile shape changes occur, revealing the source reaching the critical luminosity and transitioning between super-critical and sub-critical accretion regimes. We performed the average and phase-resolved spectral analysis, where the flux-resolved average spectra show a stable spectral evolution with luminosity. The phase-resolved spectral analysis reveals that the dependence of spectral parameters on the pulse phase varies with different luminosities.

X-ray emission from wind-driven bow shocks is both difficult to measure and predict, but may give important insights into the energy budget of the hot phase of the ISM by quantifying mixing at the interface between hot and warm gas phases. We investigate the effect of magnetic fields and numerical resolution on predicted X-ray emission and other observable properties of bow shocks, to study convergence properties and assess robustness of predicted observables from simulations. A suite of 2D and 3D HD and MHD simulations of bow shocks were run and analysed to generate synthetic emission maps and lightcurves in X-ray and infrared emission. Resolving the Kelvin-Helmholtz (KH) instability at the wind-ISM contact discontinuity is crucial for obtaining converged results and for predicting X-ray emission and the properties of the hot shocked wind. When sufficient spatial resolution is used, we measure time variation of X-ray emission of at least an order of magnitude on a timescale comparable to the advection timescale of the wake downstream from the bow shock. Good correspondence is found between 2D and 3D simulations with comparable resolution, and 3D simulations can achieve the required resolution with reasonable computing resources. Development of the KH instability is inhibited for shear flows parallel to the ISM magnetic field, compared with what is seen in the perpendicular direction, resulting in synthetic IR emission maps of bow shocks that are smooth when seen from one perspective but show strong distortions from another. Measuring the X-ray morphology and luminosity in bow shocks may be useful for constraining mixing and energy-transfer rates between hot and warm gas phases of the ISM. Dynamical instabilities at the wind-ISM interface are a crucial ingredient in determining the properties of the hot-gas phase in stellar bow-shocks, in particular to capture its time dependence.

We present a comprehensive comparison of different Markov Chain Monte Carlo (MCMC) sampling methods, evaluating their performance on both standard test problems and cosmological parameter estimation. Our analysis includes traditional Metropolis-Hastings MCMC, Hamiltonian Monte Carlo (HMC), slice sampling, nested sampling as implemented in dynesty, and PolyChord. We examine samplers through multiple metrics including runtime, memory usage, effective sample size, and parameter accuracy, testing their scaling with dimension and response to different probability distributions. While all samplers perform well with simple Gaussian distributions, we find that HMC and nested sampling show advantages for more complex distributions typical of cosmological problems. Traditional MCMC and slice sampling become less efficient in higher dimensions, while nested methods maintain accuracy but at higher computational cost. In cosmological applications using BAO data, we observe similar patterns, with particular challenges arising from parameter degeneracies and poorly constrained parameters.

This paper investigates the formation, populations, and evolutionary paths of UltraLuminous X-ray Sources (ULXs) within Globular Clusters (GCs). ULXs, characterised by their extreme X-ray luminosities, present a challenge to our understanding of accretion physics and compact object formation. While previous studies have largely focused on field populations, this research examines the unique environment of GCs, where dynamical interactions play a significant role. Using the MOCCA Monte Carlo code, we explore how dynamics influences ULX populations within these dense stellar clusters. Our findings reveal that dynamical processes, such as binary hardening and exchanges, can both facilitate and impede ULX formation in GCs. The study explores the impact of parameters including the initial binary fraction, tidal filling, and multiple stellar populations on the evolution of ULXs. We find that non-tidally filling clusters exhibit significantly larger ULX populations compared to tidally filling ones. The results indicate that the apparent scarcity of ULXs in GCs may be related to the older stellar populations of GCs relative to the field. Furthermore, the study identifies a population of "escaper" ULXs, which originate in GCs but are ejected and emit X-rays outside the cluster. These escapers may significantly contribute to the observed field ULX population.

The formation time-scales of quiescent galaxies can be estimated in two different ways, by their star formation history and by their chemistry. Previously, both methods yielded conflicting results, especially when considering $\alpha$-enhanced objects. This is primarily due to the time resolution limitations of very old stellar populations that prevent us from accurately constraining their star formation histories. We analyse the JWST observations of the extremely massive galaxy ZF-UDS-7329 at z$\sim$3.2 and show that we can achieve the higher time resolution necessary to match the chemical formation time-scales using stellar population synthesis by studying galaxies at high redshift. We compare it to the well known relic galaxy NGC 1277, arguing that ZF-UDS-7329 is an early Universe example of the cores of present day massive elliptical galaxies or, if left untouched, a relic galaxy.

A consensus has been reached in recent years that binarity plays an important role in the formation and evolution of a significant fraction of planetary nebulae (PNe). Utilizing the archived photometric data from the Zwicky Transient Facility survey, we conducted a comprehensive data mining in search for brightness variations in a large sample of Galactic PNe. This effort leads to identification of 39 PNe, whose central stars exhibit periodic variation in light curves. Among these objects, 20 are known binary central stars of PNe, while the remaining 19 are new discoveries. Additionally, we identified 14 PNe with central stars displaying anomalous variation in light curves, as well as eight variables based on the high-cadence photometric data. Among the new discoveries of periodicity, we found compelling evidence of binary systems at the centres of two archetypal quadrupolar PNe. We also report on very peculiar brightness variation observed in the central core of the compact PN NGC6833. Several PNe in our sample deserve follow-up observations, both high-dispersion spectroscopy and high-precision photometry, to reveal the true nature of their central binarity or even multiplicity.

We combine ultraviolet imaging of the 13H survey field, taken with the XMM-Newton Optical Monitor telescope (XMM-OM) and the Neil Gehrels Swift Observatory Ultraviolet and Optical Telescope (UVOT) in the UVM2 band, to measure rest-frame ultraviolet 1500A luminosity functions of star-forming galaxies with redshifts between 0.4 and 0.6. In total the UVM2 imaging covers a sky area of 641 square arcmin, and we detect 273 galaxies in the UVM2 image with 0.4<z<0.6. The luminosity function is fit by a Schechter function with best-fit values for the faint end slope alpha = -1.8 +0.4 -0.3 and characteristic absolute magnitude M* = -19.1 +0.3 -0.4. In common with XMM-OM based studies at higher redshifts, our best-fitting value for M* is fainter than previous measurements. We argue that the purging of active galactic nuclei from the sample, facilitated by the co-spatial X-ray survey carried out with XMM-Newton is important for the determination of M*. At the brightest absolute magnitudes (M1500<-18.5) the average UV colour of our galaxies is consistent with that of minimal-extinction local analogues, but the average UV colour is redder for galaxies at fainter absolute magnitudes, suggesting that higher levels of dust attenuation enter the sample at absolute magnitudes somewhat fainter than M*.

We report the stellar mass functions obtained from 20 radiation hydrodynamical simulations of star cluster formation in 500 M$_\odot$ molecular clouds with metallicities of 3, 1, 1/10 and 1/100 of the solar value, with the clouds subjected to levels of the cosmic microwave background radiation that are appropriate for star formation at redshifts z=0, 3, 5, 7, and 10. The calculations include a thermochemical model of the diffuse interstellar medium and treat dust and gas temperatures separately. We find that the stellar mass distributions obtained become increasingly bottom light as the redshift and/or metallicity are increased. Mass functions that are similar to a typical Galactic initial mass function are obtained for present-day star formation (z=0) independent of metallicity, and also for the lowest-metallicity (1/100 solar) at all redshifts up to z=10, but for higher metallicities there is a larger deficit of brown dwarfs and low-mass stars as the metallicity and redshift are increased. These effects are a result of metal-rich gas being unable to cool to as lower temperatures at higher redshift due to the warmer cosmic microwave background radiation. Based on the numerical results we provide a parameterisation that may be used to vary the stellar initial mass function with redshift and metallicity; this could be used in simulations of galaxy formation. For example, a bottom-light mass function reduces the mass-to-light ratio compared to a typical Galactic stellar initial mass function, which may reduce the estimated masses of high-redshift galaxies.

We introduce a novel neural network architecture, SkyReconNet, which is tailored to reconstruct missing information in the masked regions of the full-sky Cosmic Microwave Background (CMB) map. Inpainting the CMB maps is a particularly formidable challenge when dealing with extensive and irregular masks, such as galactic masks which can obscure substantial fractions of the sky. Our network addresses this challenge by combining the strengths of the dilated convolutional layers, with their expanded receptive fields and standard convolutional layers, to capture both the global and local features in a CMB map. This hybrid design enables our network to accurately predict CMB fluctuations in the masked regions by effectively leveraging the information from the surrounding unmasked areas. The network is trained on simulated CMB maps generated using the Code for Anisotropies in the Microwave Background (CAMB) package. The galactic region of these simulated maps are excluded using the Planck 2018 common mask provided by the Planck CMB mission. During training, the network optimizes its weights by minimizing a composite loss function that combines the Structural Similarity Index Measure (SSIM) and mean squared error (MSE). SSIM preserves the essential structural features of the CMB, ensuring an accurate and coherent reconstruction of the full-sky CMB map, while MSE minimizes the pixel-wise deviations, enhancing the overall accuracy of the predictions. The predicted CMB maps and their corresponding angular power spectra align closely with the targets, achieving the performance limited only by the fundamental uncertainty of cosmic variance. Our results demonstrate that this approach effectively addresses the challenges posed by large irregular masks, offering a significant inpainting tool in CMB and other cosmological analyses.

Large spectroscopic surveys aim to consistently compute stellar parameters of very diverse stars while minimizing systematic errors. We explore the use of stellar clusters as benchmarks to verify the precision of spectroscopic parameters in the fourth data release (DR4) of the GALAH survey. We examine 58 open and globular clusters and associations to validate measurements of temperature, gravity, chemical abundances, and stellar ages. We focus on identifying systematic errors and understanding trends between stellar parameters, particularly temperature and chemical abundances. We identify trends by stacking measurements of chemical abundances against effective temperature and modelling them with splines. We also refitted spectra in three clusters with the Spectroscopy Made Easy and Korg packages to reproduce the trends in DR4 and to search for their origin by varying temperature and gravity priors, linelists, and spectral continuum. Trends are consistent between clusters of different ages and metallicities, can reach amplitudes of ~0.5 dex and differ for dwarfs and giants. We use the derived trends to correct the DR4 abundances of 24 and 31 chemical elements for dwarfs and giants, respectively, and publish a detrended catalogue. While we couldn't pinpoint the trends' origins, we found that: i) photometric priors affect derived abundances, ii) temperature, metallicity, and continuum levels are degenerate in spectral fitting, and it is hard to break the degeneracy even by using independent measurements, iii) the extent of the linelist used in spectral synthesis is essential for cool stars, and iv) different spectral fitting codes produce significantly different iron abundances for stars of all temperatures. We conclude that clusters can be used to characterise the systematic errors of parameters produced in large surveys, but further research is needed to explain the origin of the trends.

Atmospheric escape is now considered the major contributing factor in shaping the demographic of detected exoplanets. However, inferences about the exoplanet populations strongly depend on the accuracy of the models. Direct observational tests of atmospheric models are still in their infancy. Helium escape from planetary atmospheres has rapidly become the primary observational probe, already observed in $\gtrsim$20 exoplanets. Grounding our understanding in the basic physics of atmospheric escape, we present a new theoretical model to predict the excess absorption from the helium absorption line. We constrain the atmosphere properties, such as mass-loss rates and outflow temperatures, by implementing a Parker wind solution with an energy limited evaporating outflow. Importantly, we self-consistently link the mass-loss rates and outflow temperatures, which are critical to understanding helium absorption as the triplet-level population is typically exponentially sensitive to temperature. Furthermore, helium absorption is typically optically thin and the absorption is dominated far from the planet. Therefore, the absorption depth is not a measure of the size of the helium outflow. Our results indicate that for planets with a detectable signal, typically the helium triplet population in the atmosphere rapidly approaches a statistical equilibrium between populations by recombination and depopulation caused by electron collisions. We suggest that excess helium absorption can be quantified by a scaled equivalent width, which is positively correlated with the mass loss rate. We also show that the helium absorption scales with incident radiation, particularly with the XEUV to FUV flux ratios.

The chemistry within a protoplanetary disk is greatly affected by external radiation from the local stellar environment. Previous work has focused on extreme radiation fields, representative of the center of something like the Orion Nebula Cluster. However, even in such environments, many disks exist at the edges of a cluster where the lower stellar density leads to radiation fields weaker by orders of magnitude compared to the center. We present new chemical models of a T-Tauri disk in the presence of a moderately increased interstellar radiation field (ISRF). Such an environment has a background UV strength of 10 to 100 times higher than the galactic average ISRF. Moderate radiation fields are among the most prevalent disk-harboring environments and have interesting implications for the chemistry of the outer disk radii. We find that the external UV radiation creates an outer ionization front that impacts the cold disk chemistry to varying degrees, depending on outer disk structure. Certain molecules like C$^+$, N$_2$H$^+$, C, and CS are more strongly impacted by the ISRF in their abundance, column density, and observable emission. Other abundant species like HCO$^+$ and CO are less affected by the external UV flux in the outer disk under such moderate UV conditions. Further, we demonstrate that the chemistry occurring in the inner tens of au is relatively unchanged, which suggests that even in moderately externally irradiated disks, the inner disk chemistry may be more similar to isolated disks like those in, e.g., the Taurus and Lupus star-forming regions.

Models of the galaxy-halo connection are needed to understand both galaxy clusters and large scale structure. To make said models, we need a robust method that assigns galaxies to halos and matches the observed and simulated stellar-halo mass relation. We employ an empirical Subhalo Abundance Matching (SHAM) model implemented in the halos module of SkyPy which assigns blue and red galaxies based on the Peng et al. (2010) (arXiv:1003.4747v2) model containing three parameters: $M_\mu$ (halo mass where half the galaxies assigned should be quenched), $\sigma$ (transition width from star forming to quenched) and $b$ (baseline quenched fraction at low mass). We test two sets of galaxy stellar mass functions for four populations of galaxies (central/satellite, blue/red) and run parameter estimation using Approximate Bayesian Computation over each model when compared to a set of applicable literature models. For the Weigel et al. (2016) (arXiv:1604.00008v1) galaxies we find best fit values of log $M_\mu = 11.94^{+0.02}_{-0.02}$, $\sigma = 0.49^{+0.04}_{-0.04}$ and $b = 0.31^{+0.01}_{-0.01}$. For the Birrer et al. (2014) (arXiv:1401.3162v2) galaxies we find best fit values of log $M_\mu = 11.93^{+0.01}_{-0.01}$, $\sigma = 0.53^{+0.04}_{-0.04}$ and $b = 0.51^{+0.05}_{-0.04}$. Overall, we demonstrate that these constraints produce a model that is consistent with literature models for the central galaxies. Future research will focus on the normalisation of the satellite galaxies in order to better constrain the $b$ parameter.

Quasi-periodic eruptions (QPEs) are a recently identified class of X-ray transient associated with tidal disruption events by supermassive black holes, and for which there are multiple possible explanations. In this paper we present a simple model which requires the black hole be spinning, be misaligned with the accretion flow (both conditions of which are almost certainly met) and that the accretion rate is a few times the Eddington limit. We speculate that the resulting Lense-Thirring torques force the disc and entrained outflows to precess, leading to increased X-ray flux when the wind-cone is oriented at lower inclinations to the observer. We test the range of parameters for which this model could explain the period and brightness of the QPE events discovered thus far, and make qualitative comparisons between the observed X-ray spectra and lightcurves to those extracted from GR-RMHD simulations. Overall, we find some areas of promising concordance, and identify challenges related to the details of current simulations.

We measure the radio continuum fluxes at the locations of all Gaia${-}$confirmed members of Taurus${-}$Auriga using Karl G. Jansky Very Large Array Sky Survey data (VLASS; 2${-}$4 GHz, $\sigma_{\rm{VLASS}}{\sim}110{-}140 \mu$Jy, $2.5''$ resolution) spanning 3 VLASS epochs (2019, 2021, and 2023). We present 35 detections coincident with young Taurus${-}$Auriga stars (29 in individual VLASS images, 6 via stacking). We find a strong dependence on spectral type, wherein the fractional detection rate of radio emission coincident with early-type young stellar objects (YSOs) is systematically higher than late-type YSOs, ranging from 25% $\pm$ 13% for B${-}$F YSOs, 21% $\pm$ 11% for G YSOs, 18.4% $\pm$ 6.3% for K0${-}$K4 YSOs, 15.5% $\pm$ 5.4% for K5${-}$K9 YSOs, 7.0% $\pm$ 2.7% for M0${-}$M2 YSOs, 2.3% $\pm$ 0.9% for M3${-}$M6 YSOs, and 1.9% $\pm$ 1.9% for YSOs with SpTs later than M7. We present cumulative density distributions of radio luminosity densities that demonstrate a significant luminosity enhancement for early- versus late-type YSOs. We find 25% of the detected sources to be significantly variable. We discuss possible interpretations of this dependence, which may reflect stellar magnetic activity, binary interactions, or stellar flaring. We find that mid-infrared YSO class is a strong indicator of radio detectability consistent with higher frequency Taurus-Auriga VLA surveys, with class III stars detected at a rate of 8.8% $\pm$ 1.6%, class IIs at 2.0% $\pm$ 1.2%, and combined class 0s, Is and Fs at 8.0% $\pm$ 5.4%.

Identifying useful flat-space limits for cosmological correlators, where they can be expressed in terms of observables in Minkowski space is nontrivial due to their scale-invariant nature. In recent years, it has been shown that momentum-space correlators encode flat-space amplitudes at specific singularities that emerge in the complex plane of their kinematics after analytical continuation. This flat-space limit is massless in the sense that the amplitude corresponds to the ultraviolet regime of the associated flat-space process, where the masses of the internal propagators are effectively zero. In this paper, we introduce a novel massive flat-space (MFS) limit, in which the internal masses in the corresponding flat-space Feynman graph remain finite. Our proposal applies to arbitrary graphs with light external legs and heavy internal lines, using a double-scaling limit. In this limit, the external energies, treated as independent variables, approach zero in inverse proportion to the propagator masses, which are sent to infinity. We present a general reduction formula that expresses diagrams in this limit in terms of amputated Feynman graphs in flat space. Our findings underscore the deep connections between the rich structure of massive Feynman integrals and the properties of cosmological correlators involving the exchange of heavy fields. Using this reduction formula, we compute sample one-loop contributions from heavy particles to inflationary correlators in the small sound-speed regime, revealing novel bispectrum shapes. The non-Gaussian signals we uncover, which are especially pronounced around the equilateral configuration, cannot be reproduced by adding local terms to the effective field theory of single-field inflation. Instead, they are captured by incorporating prescribed spatially non-local operators into the EFT.

Many body gravity (MBG) is an alternate theory of gravity, which has been able to explain the galaxy rotation curves, the radial acceleration relation (RAR) and the wide binary stars (WBS). The genesis of MBG is a novel theory, which models systems with thermal gradients, by recasting the variation in the temperature as a variation in the metric. Merging the above concept with Einstein's gravity, leads to the theory of thermal gravity in 5-D space-time-temperature. Thermal gravity when generalized for partially thermalized systems, results in the theory of many body gravity. The bullet cluster is supposed to be a smoking gun evidence for the presence of dark matter. However, this work demonstrates that the MBG theory can explain the weak gravitational lensing effect of the bullet cluster, without the need for yet undiscovered baryonic matter or dark matter.

Understanding astrophysical and cosmological processes can be challenging due to their complexity and lack of intuitive analogies. To address this, we present \texttt{AstronomyCalc}, a Python package specifically designed to aid university-level teaching by integrating theoretical physics with practical astronomical data analysis methods. The package enables students to solve key cosmological calculations, such as the Friedmann equations, and explore various models while visualizing how parameter variations affect cosmic dynamics. It includes tools for generating synthetic astronomical data, such as Type Ia supernova measurements, and supports analysis of publicly available datasets, including Pantheon+ and the SPARC galaxy database. Simplified implementations of advanced algorithms, such as Monte Carlo Markov Chains, allow students to engage with data analysis techniques used in contemporary research. Additionally, \texttt{AstronomyCalc} will be consistently updated with more tools and user-friendly Jupyter notebooks, making it a continually evolving educational resource for developing conceptual understanding and practical skills in astrophysics and cosmology.

The kinetically mixed dark photon is a simple, testable dark matter candidate with strong theoretical motivation. Detecting the feeble electric field dark photon dark matter produces requires extremely sensitive detectors. Bulk acoustic resonators (BARs), with their exceptionally high-quality phonon modes, are capable of achieving incredible sensitivity to gravitational waves in the MHz to GHz frequency range. The BAR phonons are typically read out by detecting the electric field generated by the BAR materials' piezoelectricity. Here we show that this piezoelectricity also rewards such detectors sensitivity to dark photon dark matter, as the dark electric field can resonantly excite BAR phonons. A single 10 g piezoelectric BAR in a large, cold, environment can be orders of magnitude more sensitive to the kinetic mixing parameter than any current experiment, with only a month-long exposure and thermally-limited backgrounds.

The study of regular black holes and black hole mimickers as alternatives to standard black holes has recently gained significant attention, driven both by the need to extend general relativity to describe black hole interiors, and by recent advances in observational technologies. Despite considerable progress in this field, significant challenges remain in identifying and characterizing physically well-motivated classes of regular black holes and black hole mimickers. This report provides an overview of these challenges, and outlines some of the promising research directions -- as discussed during a week-long focus programme held at the Institute for Fundamental Physics of the Universe (IFPU) in Trieste from November 11th to 15th, 2024.

Further bright sirens - gravitational wave events with electromagnetic counterparts - are keenly awaited, but proving elusive. The exceptional event GW170817 had a profound impact on the landscape of viable cosmological extensions of General Relativity (GR); can we expect this kind of shift to be repeated in the next decade? In this work we will assess the potential constraints from bright sirens in the LIGO-Virgo-KAGRA O5 era and third generation detector era. We set up the statistical formalism for our constraints, and generate and analyse simulated data in the context of general scalar-tensor theories. We highlight the important role that gamma-ray burst detection has in breaking key parameter degeneracies. We find that the next ten bright sirens alone will not competitively constrain cosmological gravity, but that one year of third generation observations could confidently detect mild departures from GR, e.g. the Horndeski parameter $\alpha_{\rm M}\neq 0$ is detected at greater than $3\sigma$. This justifies investment in a broad range of methods for gravitational wave cosmology (dark sirens, bright sirens and cross-correlation with large-scale structure) to ensure tests of cosmological gravity advance in both the short-term and the long-term.

LiteBIRD is an upcoming JAXA-led mission that aims to measure primordial gravitational waves in the B-mode polarization of the cosmic microwave background. It is set to launch in 2032. The LiteBIRD detector array consists of around 5000 TES detectors which are read out using digital frequency multiplexing over a bandwidth of 1-6 MHz. The multiplexing factor ranges from 58x to 68x. We are presently developing single-stage SQUID array amplifiers for LiteBIRD readout. Due to the reduced complexity and cost, and greater heritage from ground-based experiments such as the South Pole Telescope and Simons Array, single-stage SQUID array amplification is preferable for the first-stage amplification, as long as it can meet the requirements. The LiteBIRD single-stage SQUID Array is required to have high transimpedance amplification while maintaining a low input inductance and low dynamic resistance. In addition, the input-referred current noise must be very low, and the power dissipation must remain below about 100 nW. These requirements have non-trivial interactions. To maximize performance within these requirements we have performed lumped-element SQUID simulation. We find that by optimizing SQUID internal damping elements and inductive loading, good single-stage SQUID array performance can be obtained for LiteBIRD, including significant engineering margin.

We study the energy-momentum tensor of a bubble wall beyond the approximation of an infinitely thin wall. To this end, we discuss the proper decomposition into wall and bulk contributions, and we use a systematic method to calculate the energy-momentum tensor at any order in the wall width. We consider the specific examples of spherical bubbles with different initial configurations, and we compare our approximations with a numerical computation.

We review the conservation laws of magnetohydrodynamics (MHD) in an expanding homogeneous and isotropic Universe that can be applied to the study of early Universe physics during the epoch of radiation domination. The conservation laws for a conducting perfect fluid with relativistic bulk velocities in an expanding background are presented, extending previous results that apply in the limit of subrelativistic bulk motion. Furthermore, it is shown that the subrelativistic limit presents new corrections that have not been considered in previous work. Imperfect relativistic fluids are briefly described but their detailed study is not included in this work. We review the propagation of sound waves, Alfvén waves, and magnetosonic waves, as well as the Boris correction for relativistic Alfvén speeds. This review is an extension, including new results, of part of the lectures presented at the minicourse ``Simulations of Early Universe Magnetohydrodynamics'' lectured by A. Roper Pol and J. Schober at EPFL, as part of the six-week program ``Generation, evolution, and observations of cosmological magnetic fields'' at the Bernoulli Center in May 2024.

Stable $^{205}$Tl ions have the lowest known energy threshold for capturing electron neutrinos ($\nu_e$) of ${ E}_{\nu_e}\ge50.6$\,keV. The Lorandite Experiment (LOREX), proposed in the 1980s, aims at obtaining the longtime averaged solar neutrino flux by utilizing natural deposits of Tl-bearing lorandite ores. To determine the $\nu_e$ capture cross section, it is required to know the strength of the weak transition connecting the ground state of $^{205}$Tl and the 2.3 keV first excited state in $^{205}$Pb. The only way to experimentally address this transition is to measure the bound-state beta decay ($\beta_{b}$) of fully ionized $\mathrm{^{205}Tl^{81+}}$ ions. After three decades of meticulous preparation, the half-life of the $\beta_{b}$ decay of $\mathrm{^{205}Tl^{81+}}$ has been measured to be $291_{-27}^{+33}$ days using the Experimental Storage Ring (ESR) at GSI, Darmstadt. The longer measured half-life compared to theoretical estimates reduces the expected signal-to-noise ratio in the LOREX, thus challenging its feasibility.

The role of memory time fluctuations for instabilities in random media is considered. It is shown that fluctuations can result in infinitely fast growth of statistical moments. The effect is demonstrated in the framework of light propagation in the Universe, which contains curvature fluctuations while remaining homogeneous and isotropic on average.

The scientific status of physical cosmology has been the subject of philosophical debate ever since detailed mathematical models of the Universe emerged from Einstein's general theory of relativity. Such debates revolve around whether and to what extent cosmology meets established demarcation criteria for a discipline to be scientific, as well as determining how to best characterize cosmology as a science, given the unique challenges and limitations faced by a discipline which aims to study the origin, composition, and fate of the Universe itself. The present article revisits, in light of the dramatic progress in cosmology in recent decades, an earlier debate held in the 1950s between Herman Bondi and Gerald Whitrow regarding the scientific status of cosmology. We analyse cosmology's transition from an emerging science to a cornerstone of modern physics, highlighting its empirical successes in establishing the $\Lambda$-Cold Dark Matter ($\Lambda$CDM) model and in its delivery of various successful novel predictions. Despite this remarkable scientific success and progress, we argue that modern cosmology faces a further profound challenge: the permanent underdetermination of the microphysical nature of its exotic energy components: inflation, dark matter, and dark energy. Drawing historical parallels with the role of spectroscopy in revealing the microphysical nature of atomic physics, we argue that the epistemic barriers obstructing us from ascertaining the microphysical nature of these exotic energy components are significant, in turn casting doubt upon whether cosmology can ever transcend these particular epistemic challenges. We conclude by reflecting on the prospects for future breakthroughs and/or non-empirical arguments which could decide this issue conclusively.