Reanalyzing DESI DR1: 4. Percent-Level Cosmological Constraints from Combined Probes and Robust Evidence for the Normal Neutrino Mass Hierarchy

We present cosmological parameter measurements from the full combination of DESI DR1 galaxy clustering data, described with large-scale structure effective field theory. By incorporating photometric galaxies and CMB lensing cross-correlations, and extending the bispectrum likelihood to smaller scales with a consistent one-loop computation, we achieve substantial gains in constraining power. Combined with the latest DESI baryon acoustic oscillation (BAO) data and cosmic microwave background (CMB) priors on the spectral tilt and baryon density, we find, in Λ CDM, H 0=69.08 plus or minus 0.37 km, s-1 Mpc-1, Ω m=0.2974 plus or minus 0.0050, and σ 8 = 0.838 plus or minus 0.017 (S 8 = σ 8 sqrt Ω m/0.3 =0.834 plus or minus 0.018). Adding the Pantheon+ supernovae (SNe), we find a 2.2σ preference for the w 0w a dynamical dark energy model from low-redshift data alone, rising to 2.7σ when exchanging the SNe for Planck CMB data. Combining the full-shape, BAO, CMB, and SNe likelihoods improves the dark energy figure-of-merit by 15% and bounds the neutrino mass sum to M ν<0.049 eV (Λ CDM) and M ν<0.077 eV (w 0w a CDM) at 95% CL. This is the strongest w 0w a CDM bound to date, 37% tighter than from the background expansion data alone. The preference for the normal neutrino mass ordering thus holds regardless of the background model: the inverted hierarchy is disfavored at about,3.5σ in Λ CDM and about,2.4σ in w 0w a CDM, with the latter constraint free of the geometric tension between CMB and BAO that is known to sharpen the Λ CDM bound.