2025-12-04 | | Total: 28
Dark matter (DM) annihilation can power the first generation of stars as long lived dark stars (DSs) that grow to supermassive scales $M_{\rm DS}\gtrsim 10^{5} M_{\odot}$ and eventually collapse into heavy black holes that could seed the supermassive black holes observed at high redshifts. We compute the diffuse electromagnetic emission from a cosmological population of such supermassive DSs and their black hole remnants, tracking the entire DS history and including thermal surface radiation, DM annihilation in adiabatically contracted halos as well as late-time emission from DM overdensity spikes around the resulting black holes. After accounting for photon attenuation, we find that DS related contributions can exceed the Fermi-LAT extragalactic $γ$-ray background for thermal relic annihilation cross-sections and DM masses below $\sim 1$ TeV. Our results constitute the first population integrated diffuse multimessenger constraints on supermassive DSs as progenitors of early black holes and demonstrate that diffuse photon and neutrino backgrounds offer a powerful and complementary avenue for probing the role of DM in the formation of the earliest massive structures.
We investigate the impact of shear viscosity on the propagation of primordial gravitational waves (pGW) after inflation. Without assuming a specific inflationary scenario we focus on the evolution of pGWs after they re-enter the horizon during a cosmological epoch characterized by the presence of shear viscosity. We show that shear viscosity introduces an additional damping term in the tensor equation, modifying both the transfer function and the energy density power spectrum. For a constant shear viscosity-to-Hubble ratio the transfer function acquires an extra red tilt, while a time-dependent viscosity leads to a running spectral index $Ω_\text{GW}\sim k^{n_\text{eff}(k)}$ controlled by the time evolution of the mean free path of the viscous fluid. Our analysis provides a general framework to analytically quantify how shear viscosity can alter the primordial gravitational wave background in standard and non-standard post-inflationary scenarios. As a case study we evaluate the effect of viscosity of the electron-photon-baryon plasma, on both the transfer function and the normalized energy density, finding a $k$-dependent blue tilt due to gravitational wave freeze-out from the viscous phase. This effect corresponds to a fractional difference of order $10^{-3}$.
[Abridged] The primary gravitational waves (PGWs) are generated in the early universe from the quantum vacuum during inflation. In slow roll inflation, the power spectrum (PS) of PGWs over large scales, which leave the Hubble radius during inflation, is nearly scale-invariant. However, over very small scales, which never leave the Hubble radius, the PS of PGWs behaves as k^2, where k denotes the wave number. We examine the PS of PGWs at such high wave numbers or frequencies when the PGWs are evolved post-inflation, through the epochs of radiation and matter domination. Firstly, we argue that the PS has to be regularized in order to truncate the unphysical k^2 rise at high frequencies. Assuming instantaneous transitions from inflation to the epochs of radiation and matter domination, we carry out the method of adiabatic regularization to arrive at the PS of PGWs over a wide range of frequencies. We show that the process of regularization truncates the k^2 rise and the PS of PGWs oscillates with a fixed amplitude about a vanishing mean value over small scales or, equivalently, at high frequencies. Secondly, we smooth the transition from inflation to radiation domination (to be precise, we smooth the 'effective potential' governing the equation of motion of PGWs) and examine the impact of the smoothing on the regularized PS of PGWs. With the help of a linear smoothing function, we explicitly show that the smoother transition leads to a power-law suppression in the amplitude of the oscillations (about the zero mean value) of the regularized PS of PGWs over small scales that never leave the Hubble radius during inflation. Our analysis indicates that, when transitions are involved, regularization as well as smooth transitions seem essential to ensure that the correlation functions of the PGWs in real space are well behaved. We discuss the directions in which our results need to be extended.
There is tentative evidence for a mismatch between the rest frames of matter and the cosmic microwave background, the "quasar dipole anomaly". We consider such a dipole in tilted anisotropic models, for a range of scenarios and sources: spatial curvature, cosmic heat flux, large scale electromagnetic fields and a Khronon field. Crucially, we determine the ancillary effects on other cosmological observables in each of these models and we show that, apart from the case of the Khronon field, it is unlikely that one can obtain a dipole with the amplitude that is being observed unless one considers additional exotica.
Current and forthcoming galaxy surveys will map the observable Universe with unprecedented depth, sky coverage, and precision. These maps are affected by relativistic redshift-space distortions (RSDs), which become increasingly relevant on ultra-large scales. Accurate modelling of these relativistic RSDs is essential to avoid systematic biases in key cosmological measurements, such as primordial non-Gaussianity (PNG). To address this, we introduce an updated implementation of the LIGER method, LIGER4GAL, which incorporates all linear-order relativistic RSDs directly at the tracer level of high-resolution N-body simulations. We demonstrate that LIGER4GAL improves upon previous iterations of the LIGER method by reproducing the expected non-linear clustering while maintaining accuracy for relativistic RSDs on large scales. We use the updated code to generate a DESI-like sample of luminous red galaxies from the Huge MultiDark Planck simulation. By measuring the power spectrum multipoles of this sample with and without the imprint of relativistic RSDs, we assess the impact of relativistic effects on measurements of the local PNG signal ($f_\mathrm{nl}$). We find that the omission of the''finger-of-the-observer'' (sourced by the peculiar velocity of the observer) effect in the power spectrum modelling can bias measurements of $f_{\rm nl}$ by more than $1$ ($0.25$) $ σ_{f_{\rm nl}}$ in 40% (80%) of the possible realizations of the universe if scales down to $k_\mathrm{min} = 0.0015\,h/\mathrm{Mpc}$ are included.
Type Ia supernovae (SNe Ia) provide the most robust means of measuring extragalactic distances. While most of the effort has focused on increasing the number of SNe Ia observed in the optical, near-infrared (NIR) observations remain scarce despite their advantages, that is, reduced dust extinction and a more intrinsic standard candle behavior, requiring little to no empirical corrections. Here, we present ASNOS (ANDICAM-SOFI Near-infrared and Optical type Ia Supernova), a dataset with sample size of 1,482 epochs in the $BVRIYJH$ filters from the ANDICAM instrument on the 1.3-meter SMARTS telescope at Cerro Tololo Inter-American Observatory, along with 125 $JHK$ epochs from the SOFI instrument on the 3.58-meter New Technology Telescope on the La Silla Observatory. Additionally, we incorporate optical forced photometry from the Zwicky Transient Facility and the Asteroid Terrestrial-impact Last Alert System. The sample comprises 41 SNe Ia in total, including 29 normal events, eight 1991T-like objects, and four peculiar subtypes, all located at redshifts $z < 0.085$. This paper provides a detailed overview of the ASNOS sample selection, data reduction, SN photometry, host-galaxy spectral energy distribution construction, both global and local, and SN light-curve fitting using three methods: SALT3-NIR, SNooPy, and BayeSN. A companion paper will present the cosmological analysis.
Recent results from the DESI collaboration suggest a preference for an evolving dark energy (DE) component rather than a cosmological constant, motivating the exploration of alternative models for the background expansion. These data also reveal tension in the inferred matter density parameter -lower in DESI and higher in Planck- as well as a neutrino mass posterior that approaches the lower bounds permitted by oscillation experiments. In this work, we propose and test a conformal holographic DE (CHDE) model in which the DE density depends on a power law of the conformal time, characterized by an exponent (n). This formulation introduces a single additional parameter relative to LambdaCDM and reduces to it in the limit n = 0. We confront the CHDE model with BAO, CMB, and supernova datasets, following the same combinations used by DESI, and perform parameter inference under both flat and non-flat cosmologies. Our analyses show that LambdaCDM is not favored as the best-fit model when using CMB data alone or in joint analyses including BAO and SNla, and it is disfavored at the 4.4 sigma level for non-flat model and 4.5 sigma for the flat model. We obtain consistent values of n= -0.28 to -0.32 with uncertainties less than +-0.1 across multiple data combinations. Similar to LambdaCDM, the CHDE model predicts a lower matter density when employing DESI data instead of Planck data. This, in turn, influences the neutrino mass constraints, yielding values close to the minimal allowed range. Despite these dataset-dependent tensions, both the flat and curved CHDE models remain compatible with neutrino mass constraints from terrestrial experiments and yield posterior distributions that peaks at positive values. This behavior avoids the issue encountered in the LambdaCDM model, where the posterior peaks at negative mass values.
Ultra-light bosonic dark matter (ULDM) is an interesting and promising dark matter candidate. While the wave-like nature of ULDM has been widely studied in the literature, we explore another distinctive feature of ULDM as Bose-Einstein Condensate (BEC) in this paper: the emergence of vortices in a rotating BEC-ULDM halos. Using numerical solution of the GPP equation, we demonstrate that a lattice of vortices ,underdensity columns that carry angular momentum, naturally forms in a ULDM halo under conditions similar to those of the Milky Way. Furthermore, we study the gravitational lensing by these vortices as a possible observational signature of BEC-ULDM. If the vortices are large enough and the halo's rotational axis align with the line of sight, regularly separated brightness anomalies can be produced, providing strong evidence for BEC-ULDM.
In this work we forecast constraints on models of decaying dark matter (DCDM) by using redshift-space-distortion (RSD) measurements implemented in a Fisher information matrix. In particular, we focus on the fraction of unstable dark matter, $α_{\mathrm{dcdm}}$ respect to the ordinary CDM component, and the decay rate, $Γ_{\mathrm{dcdm}}$ as the key parameters of the model. Fiducial values are derived from a MontePython MCMC analysis. The derivatives of the growth-related observable, $fσ_8(z)$ with respect to the parameters are numerically around the fiducial model. For the Fisher analysis, we employ mock data designed for upcoming surveys, particularly Euclid and BOSS, where RSD measurements yield constraints on $fσ_8.$ Our results show that when both stable and unstable components are allowed, constrains on the DCDM lifetime remain weak, with $τ_{\mathrm{dcdm}}>1.18$ Gyr. In the limiting case of fully unstable dark matter ($α_{\mathrm{dcdm}}=1$), the uncertainty improves to $τ_{\mathrm{dcdm}}>235.89$ Gyr. Our findings highlight the potential of RSD probes in testing and complementing decaying dark matter scenarios.
The effective field theory (EFT) of dark energy provides a model-independent framework for studying cosmology within scalar-tensor theories. In this work, we explore how the time evolution of the cosmological background, inferred from cosmic chronometer measurements of the Hubble parameter, can be used to reconstruct the relevant EFT functions. Our approach enables the direct determination of these EFT functions from observational data without assuming any specific cosmological model. This makes it possible to test the background evolution of a wide range of dark energy models, including the $Λ$CDM model. We further demonstrate how the reconstructed EFT functions can be applied to constrain concrete theories, such as the quintessence model.
The Dark Energy Spectroscopic Instrument (DESI) in its first Data Release (DR1) already provides more than 100,000 galaxies with relative distance measurements. The primary purpose of this paper is to perform the calibration of the zero-point for the DESI Fundamental Plane and Tully-Fisher relations, which allows us to measure the Hubble constant, $H_0$. This sample has a lower statistical uncertainty than any previously used to measure $H_0$, and we investigate the systematic uncertainties in absolute calibration that could limit the accuracy of that measurement. We improve upon the DESI Early Data Release Fundamental Plane $H_0$ measurement by a) using a group catalog to increase the number of calibrator galaxies and b) investigating alternative calibrators in the nearby universe. Our baseline measurement calibrates to the SH0ES/Pantheon+ type Ia supernovae, and finds $H_0=73.7\pm 0.06\;(\text{stat.})\pm 1.1\;(\text{syst.})$ km s$^{-1}$ Mpc$^{-1}$. Calibrating to surface brightness fluctuation (SBF) distances yields a similar $H_0$. We explore measurements using other calibrators, but these are currently less precise since the overlap with DESI peculiar velocity tracers is much smaller. In future data releases with an even larger peculiar velocity sample, we plan to calibrate directly to Cepheids and the tip of the red giant branch, which will enable the uncertainty to decrease towards a percent-level measurement of $H_0$. This will provide an alternative to supernovae as the Hubble flow sample for $H_0$ measurements.
The large-scale structure of the Universe and its evolution encapsulate a wealth of cosmological information. A powerful means of unlocking this knowledge lies in measuring the auto-power spectrum and/or the cross-power spectrum of the galaxy density and momentum fields, followed by the estimation of cosmological parameters based on these spectrum measurements. In this study, we generalize the cross-power spectrum model to accommodate scenarios where the density and momentum fields are derived from distinct galaxy surveys. The growth rate of the large-scale structures of the Universe, commonly represented as $fσ_8$, is extracted by jointly fitting the monopole and quadrupole moments of the auto-density power spectrum, the monopole of the auto-momentum power spectrum, and the dipole of the cross-power spectrum. Our estimators, theoretical models and parameter-fitting framework have been tested using mocks, confirming their robustness and accuracy in retrieving the fiducial growth rate from simulation. These techniques are then applied to analyze the power spectrum of the DESI Bright Galaxy Survey and Peculiar Velocity Survey, and the fit result of the growth rate is $fσ_8=0.440^{+0.080}_{-0.096}$ at effective redshift $z_{\rm eff}=0.07$. By synthesizing the fitting outcomes from correlation functions, maximum likelihood estimation and power spectrum, yields a consensus value of $fσ_8(z_{\rm eff}=0.07) = 0.450 ^{+0.055}_{-0.055}$, and correspondingly we obtain $γ=0.580^{+0.110}_{-0.110}$, $Ω_\mathrm{m}=0.301^{+0.011}_{-0.011}$ and $σ_8=0.834^{+0.032}_{-0.032}$. The measured $fσ_8$ and $γ$ are consistent with the prediction of the $Λ$ Cold Dark Matter Model and General Relativity.
Joint analysis of the local peculiar velocity and galaxy density fields offers a promising route to testing cosmological models of gravity. We present a measurement of the normalised growth rate of structure, $fσ_8$, from the two-point correlations of velocity and density tracers from the DESI DR1 Peculiar Velocity and Bright Galaxy Surveys, the largest catalogues of their kind assembled to date. We fit the two-point correlation measurements with non-linear correlation function models, constructed from density and momentum power spectra generated using 1-loop Eulerian perturbation theory, and validate our methodology using representative mock catalogues. We find $fσ_8 = 0.391^{+0.080}_{-0.081}$, consistent to within $1σ$ with accompanying analyses of the same datasets using power spectrum and maximum-likelihood fields methods. Combining these growth rate results from different methods including appropriate correlations, we find a consensus determination $fσ_8(z = 0.07) = 0.4497 \pm 0.0548$, consistent with predictions from \textit{Planck}$+Λ$CDM cosmology. Jointly fitting to this consensus low-redshift growth rate and the DESI DR1 full-shape clustering dataset, we measure gravitational growth index $γ_{\rm L} = 0.580^{+0.110}_{-0.110}$, consistent with the prediction of general relativity.
We present the constraint on the growth rate of structure from the combination of DESI DR1 BGS sample, Fundamental Plane, and Tully-Fisher peculiar velocity catalogues using the maximum likelihood fields method. The combined catalogue contains 415,523 galaxy redshifts and 76,616 peculiar velocity measurements. To handle the large amount of data in the DESI DR1 peculiar velocity catalogue, we significantly improve the computational efficiency by rewriting the algorithm with JAX. After removing outliers and Tully-Fisher galaxies that are affected by systematics, we find $fσ_8 = 0.483_{-0.043}^{+0.080}(\mathrm{stat}) \pm 0.018(\mathrm{sys})$, consistent within $1σ$ with the power spectrum and correlation function analysis using the same dataset. Combining all three measurements with appropriate correlations, the consensus measurement is $fσ_8 (z_{\mathrm{eff}}=0.07) = 0.450\pm0.055$, consistent with Planck $+Λ$CDM cosmology $(fσ_8 = 0.449 \pm 0.008)$. Combining with the high redshift growth rate of structure measurements from DESI ShapeFit, the constraint on the growth index is $γ= 0.58\pm0.11$, consistent with GR.
We describe the production of the official set of mock catalogs for the Dark Energy Spectroscopic Instrument Peculiar Velocity Survey (DESI-PV) Data Release 1 (DR1). Our mock catalogs reproduce the Bright Galaxy Survey number density and clustering at low redshift $(z<0.1)$ and the DESI PV samples of Fundamental plane and Tully-Fisher distances, from which we derive peculiar velocities. We carefully match mock and data properties and we mimic measurements of distance indicators and peculiar velocities, which follow the same statistical properties as real data. Mock samples of type-Ia supernovae also complement the other two distance indicators. Our 675 available mock realizations were used consistently by our three different methodologies described in our companion papers that measure the growth rate of structure $fσ_8$ with DESI PV DR1. Those mocks allow us to perform precise tests of clustering models and uncertainty estimation to an unprecedented level of accuracy, and compute correlations between methodologies. The consensus value for the DESI DR1 PV growth rate measurement is $fσ_8 = 0.450 \pm 0.055$. This sample of mock catalogs represents the largest and most realistic set for cosmological measurements with peculiar velocities to date.
We calibrate the Tully-Fisher relation (TFR) using observations of spiral galaxies taken during the first year (DR1) of the DESI galaxy redshift survey. The rotational velocities of 10,262 galaxies are measured at 0.4 R26 by comparing the redshifts at 0.4 R26 with those at the galaxy centers of spatially- resolved galaxies targeted as part of the DESI Peculiar Velocity Survey. The DESI DR1 TFR slope is calibrated by separating the spiral galaxies into redshift bins of width dz = 0.005 from 0.03 < z < 0.1 and jointly fitting the TFR across all bins. We find a slope of -7.22+/-0.01 AB mag in the r-band for the TFR, with an intrinsic scatter of 0.466+/-0.001 AB mag. We present a catalog of the distances and peculiar velocities to these 10,262 galaxies using our calibrated TFR. For cosmological analyses, we also present a clustering catalog and associated random catalogs using a subset of 6807 of the DESI DR1 TF galaxies.
Measurements of peculiar velocities in the local Universe are a powerful tool to study the nature of dark energy at low ($z < 0.1$) redshifts. Here we present the largest single set of $z<0.1$ peculiar velocity measurements to date, obtained using the Fundamental Plane (FP) of galaxies in the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI). We describe the photometric and spectroscopic selection criteria used to define the sample, as well as extensive quality control checks on the photometry and velocity dispersion measurements. Additionally, we perform detailed systematics checks for the many analysis parameters in our pipeline. Our DESI DR1 catalogue contains FP-based distances and peculiar velocities for $98,292$ unique early-type galaxies, increasing the total number of $z < 0.1$ FP distances ever measured by a factor of $\sim2$. We achieve a precision of $26\%$ random error in our distance measurements which is comparable to previous surveys. A series of companion DESI papers use the distances and peculiar velocities presented in this paper to measure cosmological parameters.
The extragalactic magnetic field could be detected by searching for signatures of the electromagnetic cascade initiated by high-energy photons on the intergalactic radiation and deflected by the field. This process produces a time delay and an extended gamma-ray halo around the source, which are looked for. We propose a new signature of electromagnetic echoes: the asymmetry of the gamma-ray distribution around blazars. As a measure of asymmetry, we use the offset of the gamma-ray distribution to the location of the blazar. This offset is due to the tilt of the jet of the blazar relative to the line of sight. Using a subsample of the 10 brightest BL Lacs, we exclude the range of extragalactic magnetic fields from $10^{-16}$ to $10^{-14}$ G, assuming that these objects have maintained a constant average luminosity over hundreds of thousands of years.
We analyze the origin of the large scale bulk flow using the CosmicFlows 4 (CF4) peculiar velocity catalog. We decompose the observed motions into internal components, generated by mass fluctuations within 200Mpc/h, and external ones arising from structures beyond this volume. A weighted average technique is developed to test the model's self consistency while minimizing the impact of non Gaussian distance errors. The CF4 velocities show excellent agreement with the predicted internal field, yielding beta = 0.31 pm 0.01. We also determine that the value of the Hubble constant that should be used for calculating peculiar velocities from the CF4 to be H0 = 75.9 pm 0.1 km/s/ Mpc, consistent with CF4 calibrations. Using the minimum variance formalism, we further separate the bulk flow into its internal and external contributions and find that the observed large scale bulk flow is dominated by sources beyond 200Mpc/h. The amplitude of this externally driven flow increases monotonically with scale, consistent with the influence of a distant, massive overdensity. These findings reinforce the reliability of the CF4 velocity field while calling into question the assumption of a spatially uniform flow generated by external sources. Our results challenge the commonly made assumption that the flow in our local volume due to external mass concentrations can be modeled as being spatially uniform.
One possible way to explain the observed effective dark energy equation of state crossing $w=-1$ (the phantom divide) is through modified gravity. A key point is to not view the expansion history in isolation but to take into account the other gravitational impacts on growth of large scale structure, lensing, etc. Within shift symmetric Horndeski gravity this implies three main paths for the late time cosmic expansion. All require unusual kinetic structure and we analyze their various implications for how $w$ should behave after phantom crossing.
The correlation between galaxy stellar mass and gas-phase metallicity, known as the mass-metallicity relation (MZR), gives key insights into the processes that govern galaxy evolution. However, unquantified observational and selection biases can result in systematic errors in attempts to recover the intrinsic MZR, particularly at higher redshifts. We characterize the MZR at z~3-6 within a fully Bayesian framework using JWST NIRSpec spectra of 193 galaxies from the RUBIES survey. We forward model the observed mass-metallicity surface using prospector-generated spectra to account for two selection biases: the survey selection function and success in observing high signal-to-noise emission lines. We demonstrate that the RUBIES selection function, based on F444W magnitude and F150W-F444W color, has a negligible effect on our measured MZR. A correct treatment of the non-Gaussian metallicity uncertainties from strong-line calibrations lowers the derived MZR normalization by 0.2 dex and flattens the slope by ~20%; forward-modeling the effect of emission line observability steepens the slope by ~15%. Both of these biases must be taken into account in order to properly measure the intrinsic MZR. This novel forward modeling process motivates careful consideration of selection functions in future surveys, and paves the way for robust, high-redshift chemical enrichment studies that trace the evolution of the mass-metallicity relation across cosmic time.
Dark photons with masses in the 1-100 eV range can be produced in the Sun and subsequently absorbed in terrestrial detectors, offering a promising avenue for probing hidden-sector physics beyond the Standard Model. In this work, we develop a theoretically grounded strategy to identify optimal detector materials for solar dark photon absorption. Our strategy builds on a material-independent upper limit on the absorption rate, which we derive from Kramers-Kronig relations applied separately to the longitudinal and transverse dark photon modes. We show how the optimal material properties depend on the dark photon mass relative to the detector's plasma frequency, identifying the conditions under which a detector can saturate the theoretical upper limit. We then assess the performance of commonly used detector materials in light of these criteria and comment on the prospects of metamaterials featuring tunable plasma frequencies. Our results provide a general and model-independent framework to effectively guide the design of next-generation experiments targeting solar dark photons.
We discuss graviton dark radiation produced by the decay of a secluded scalar field that couples to the Standard Model (SM) only through gravity. Such scalar fields are long-lived, and their decay channels generically include gravitons. If such particles existed and dominated the early universe, a sizable branching ratio into gravitons would yield non-negligible dark radiation that significantly alters the subsequent thermal history of the universe. In this work, we focus on the dark glueball as a representative secluded hidden scalar and compare the decay rates into SM particles via a non-minimal coupling to gravity with those into gravitons, paying attention to how the breaking of conformal invariance affects the amount of graviton dark radiation. We find that decays into the SM are dominated by two-body decay channels into Higgs bosons and gluons. In particular, when the Higgs field has a large non-minimal coupling to gravity, the production of graviton dark radiation is naturally suppressed in the metric formalism, and the SM sector is preferentially reheated and energy transfer to other hidden sectors is suppressed. Finally, we present the expected gravitational-wave spectrum resulting from dark glueball domination.
Covariant $f(Q)$ gravity is a viable extension of General Relativity, however its strong-field predictions remain largely untested. Using the static, spherically symmetric black-hole solutions of the theory, we confront it with the most stringent probes available: black-hole shadows, Event Horizon Telescope (EHT) measurements, S2-star precession, and strong gravitational lensing. We show that the two admissible solution branches behave very differently: Case~I produces negligible deviations from Schwarzschild solution, whereas Case~II yields significant, potentially observable corrections to the photon sphere and shadow size. From the EHT shadow diameters of M87* and Sgr~A*, we obtain tight bounds, which are further strengthened by strong-lensing coefficients. These results provide the sharpest strong-field constraints on covariant $f(Q)$ gravity to date, and point toward future tests using next-generation horizon-scale imaging and precision Galactic-center astrometry.
Time-delay cosmography leverages strongly lensed quasars to measure the Universe's current expansion rate, H_0, independently from other methods. While the latest TDCOSMO results relied mainly on quadruply lensed quasars, doubly lensed systems are far more common and offer precise time delays, potentially enlarging the usable sample by a factor of five and enabling percent-level constraints on H_0. We present the first TDCOSMO analysis of a doubly imaged source, HE1104-1805, including the measurement of the four necessary ingredients. First, by combining 17 years of data from the SMARTS, Euler and WFI telescopes, we measure a time delay of 176.3\pm 10.8 days. Second, using MUSE data, we extract stellar velocity dispersion measurements in three radial bins with up to 5% precision. Third, employing F160W HST imaging for lens modelling and marginalising over various modelling choices, we measure the Fermat potential difference between the images. Fourth, using wide-field imaging, we measure the convergence added by objects not included in the lens modelling. Hence, we measure the time delay distance and the angular diameter distance to the deflector, favouring a power-law mass model over a baryonic and dark matter composite model. The measurement was performed blindly and yielded H_0 = 64.2^{+5.8}_{-5.0} x $λ_{int} km s^{-1} Mpc^{-1}, where λ_{int} is the internal mass sheet degeneracy parameter. This is in agreement with the TDCOSMO-2025 milestone and its precision for λ_{int}=1 is comparable to that obtained with the best-observed quadruply lensed quasars (4-6%). This work is a stepping stone towards a precise measurement of H_0 using a large sample of doubly lensed quasars, supplementing the current sample. The next TDCOSMO milestone paper will include this system in its hierarchical analysis, constraining λ_{int} and H_0 jointly with multiple lenses.