For the first time in nearly a decade, a new, bright transient was detected in the central parsec (pc) of the Galaxy. MAXI J1744-294 was never observed in outburst prior to January 2025. We present the results of a broadband, multi-wavelength study of this enigmatic source, including data from the NuSTAR, Chandra, XMM-Newton, Swift, and NICER X-ray telescopes, as well as complementary radio and near-infrared observations. We find that MAXI J1744-294 remained in the bright/soft state throughout the first months of 2025. Spectral hardening was observed in April 2025, followed by a decline in flux. Based on the spectral and temporal characteristics of the source, we identify MAXI J1744-294 as a candidate black hole (BH) low-mass X-ray binary (LMXB) $-$ the fourth candidate BH transient discovered within a (projected) distance of one pc from the Galactic supermassive black hole Sgr A*. This discovery provides further evidence for a cusp of BH-LMXBs in the central pc of our Galaxy, as argued for in previous observational work and suggested by analytical and theoretical work. Our multi wavelength study, involving a complementary range of observatories and spanning different outburst states, can serve as a model for future time domain astrophysics research.

With Cosmicflows-4, distances are compiled for 55,877 galaxies gathered into 38,065 groups. Eight methodologies are employed, with the largest numbers coming from the correlations between the photometric and kinematic properties of spiral galaxies (TF) and elliptical galaxies (FP). Supernovae that arise from degenerate progenitors (type Ia Sne) are an important overlapping component. Smaller contributions come from distance estimates from the surface brightness fluctuations of elliptical galaxies and the luminosities and expansion rates of core collapse supernovae (SNII). Cepheid period-luminosity relation and tip of the red giant branch observations founded on local stellar parallax measurements along with the geometric maser distance to NGC 4258 provide the absolute scaling of distances. The assembly of galaxies into groups is an important feature of the study in facilitating overlaps between methodologies. Merging between multiple contributions within a methodology and between methodologies is carried out with Bayesian Markov chain Monte Carlo procedures. The final assembly of distances is compatible with a value of the Hubble constant of $H_0=74.6$ km s$^{-1}$ Mpc$^{-1}$ with the small statistical error of $\pm 0.8$ km s$^{-1}$ Mpc$^{-1}$ but a large potential systematic error of ~3 km s$^{-1}$ Mpc$^{-1}$. Peculiar velocities can be inferred from the measured distances. The interpretation of the field of peculiar velocities is complex because of large errors on individual components and invites analyses beyond the scope of this study.

The Extragalactic Distance Database (EDD) was created as a repository for high quality, redshift-independent distances. A key component of EDD is the Color Magnitude Diagrams/Tip of the Red Giant Branch (CMDs/TRGB) catalog, which provides information on the stellar content of nearby galaxies observed with the Hubble Space Telescope (HST). Here we provide a decadal update to this catalog, which has now doubled in size to over 500 galaxies. We highlight the additions to our data reduction and analysis techniques, and provide examples of the science that has been made possible with this large data set. We find the TRGB to be a reliable measure for distance, and we aim to extend its distance coverage with HST to every galaxy within 10 Mpc. In the near-future, the combination of the James Webb Space Telescope and the Nancy Grace Roman Space Telescope will dramatically increase the number of targets within our grasp.

We present high cadence optical and ultraviolet observations of the Type II supernova (SN), SN 2022jox which exhibits early spectroscopic high ionization flash features of \ion{H}{1}, \ion{He}{2}, \ion{C}{4}, and \ion{N}{4} that disappear within the first few days after explosion. SN 2022jox was discovered by the Distance Less than 40 Mpc (DLT40) survey $\sim$0.75 days after explosion with followup spectra and UV photometry obtained within minutes of discovery. The SN reached a peak brightness of M$_V \sim$ $-$17.3 mag, and has an estimated $^{56}$Ni mass of 0.04 M$_{\odot}$, typical values for normal Type II SNe. The modeling of the early lightcurve and the strong flash signatures present in the optical spectra indicate interaction with circumstellar material (CSM) created from a progenitor with a mass loss rate of $\dot{M} \sim 10^{-3}-10^{-2}\ M_\odot\ \mathrm{yr}^{-1}$. There may also be some indication of late-time CSM interaction in the form of an emission line blueward of H$\alpha$ seen in spectra around 200 days. The mass-loss rate is much higher than the values typically associated with quiescent mass loss from red supergiants, the known progenitors of Type II SNe, but is comparable to inferred values from similar core collapse SNe with flash features, suggesting an eruptive event or a superwind in the progenitor in the months or years before explosion.

The W. M. Keck Observatory Archive (KOA) has released the Observers Data Access Portal (ODAP), a web-application that delivers astronomical data from the W. M. Keck Observatory to the scheduled program's principal investigator and their collaborators anywhere in the world in near real-time. Data files and their associated metadata are streamed to a user's desktop machine moments after they are written to disk and archived in KOA. The ODAP User Interface is built in React and uses the WebSocket protocol to stream data between KOA and the user. This document describes the design of the tool, challenges encountered, shows how ODAP is integrated into the Keck observing model, and provides an analysis of usage metrics.

We analyze the \ovi\ content and kinematics for 126 {\hi}-selected absorbers at $0.14 \lesssim z \lesssim 0.73$ for which the metallicities of their cool photoionized phase have been determined. We separate the absorbers into 100 strong {\lya} forest systems (SLFSs with $15 \la$\,{\colden}\,< 16.2) and 26 partial Lyman Limit systems (pLLSs with $16.2\le$\,{\colden}\,$\le 17.2$). The sample is drawn from the COS CGM Compendium (CCC) and has \ovi\ coverage in $\sn \geq 8$ {\it HST}/COS G130M/G160M QSO spectra, yielding a $2\sigma$ completeness level of {\coldenovi}$\,\geq 13.6$. The \ovi\ detection rates differ substantially between low-metallicity (LM; {\xh} $\leq -1.4$) and high-metallicity (HM; {\xh} > -1.4) SLFSs, with 20\% and 60\% detection rates, respectively. The \ovi\ detection frequency for the HM and LM pLLSs is, however, similar at $\sim$60\%. The SLFSs and pLLSs without detected \ovi\ are consistent with the absorbing gas being in a single phase, while those with \ovi\ trace multiphase gas. We show that the \ovi\ velocity widths and column densities have different distributions in LM and HM gas. We find a strong correlation between \ovi\ column density and metallicity. The strongest (\coldenovi$\,\ga 14$) and broadest {\ovi} absorbers are nearly always associated with HM absorbers, while weaker \ovi\ absorbers are found in both LM and HM absorbers. From comparisons with galaxy-selected and blind \ovi\ surveys, we conclude absorbers with \coldenovi$\,\ga 14$ most likely arise in the circumgalactic medium (CGM) of star-forming galaxies. Absorbers with weak \ovi\ likely trace the extended CGM or intergalactic medium (IGM), while those without \ovi\ likely originate in the IGM.

We present the analysis of an ancient galaxy at $z=2.675$ which we dub ``Eridu.'' Simultaneously modeling the JWST/NIRSpec G140M and G235M spectra from the SMILES program and $0.4-25\ \mu\mathrm{m}$ HST, JWST/NIRCam, and JWST/MIRI photometry from the the JADES+SMILES photometric catalogs shows that Eridu is massive and quiescent with stellar mass $\log(M_*/\mathrm{M_\odot})=10.96^{+0.01}_{-0.01}$ and average star formation rate <1\ \mathrm{M_\odot\ yr^{-1}} over the last 100 Myr. Star formation histories inferred from various models produce disconcertingly early and fast formation within $\sim300$ Myr of the Big Bang and quenching 2 Gyr prior to observation ($z\sim10$). This stellar mass assembly implies that the progenitor of Eridu had $M_*\approx10^{11}\ \mathrm{M_\odot}$ at z>10, nearly two orders of magnitude more than the most massive current high redshift observations. From Eridu's spectrum we infer $\mathrm{[Mg/Fe]} =+0.65^{+0.20}_{-0.19}$, indicating its stellar population is extremely $\alpha$-enhanced, which is consistent with the rapid formation timescale inferred from its star formation history. Eridu inhabits a massive protostructure which offers additional explanations for rapid mass assembly and quenching via environmental mechanisms, e.g. major mergers. Though its inferred formation is at odds with observations of the brightest cosmic dawn galaxies, we anticipate that future high-redshift galaxy formation models and sophisticated stellar population modeling codes will unearth how Eridu formed at the dawn of time.

We present near-infrared spectroscopic confirmations of a sample of 16 photometrically-selected galaxies with stellar masses log(M_star/M_sun) > 11 at redshift z > 3 from the XMM-VIDEO and COSMOS-UltraVISTA fields using Keck/MOSFIRE as part of the MAGAZ3NE survey. Eight of the ultra-massive galaxies (UMGs) have specific star formation rates (sSFR) < 0.03 Gyr-1, with negligible emission lines. Another seven UMGs show emission lines consistent with active galactic nuclei and/or star formation, while only one UMG has sSFR > 1 Gyr-1. Model star formation histories of these galaxies describe systems that formed the majority of their stars in vigorous bursts of several hundred Myr duration around 4 < z < 6during which hundreds to thousands of solar masses were formed per year. These formation ages of < 1 Gyr prior to observation are consistent with ages derived from measurements of Dn(4000) and EW0(H\delta). Rapid quenching followed these bursty star-forming periods, generally occurring less than 350 Myr before observation, resulting in post-starburst SEDs and spectra for half the sample. The rapid formation timescales are consistent with the extreme star formation rates observed in 4 < z < 7 dusty starbursts observed with ALMA, suggesting that such dusty galaxies are progenitors of these UMGs. While such formation histories have been suggested in previous studies, the large sample introduced here presents the most compelling evidence yet that vigorous star formation followed by rapid quenching is almost certainly the norm for high mass galaxies in the early universe. The UMGs presented here were selected to be brighter than Ks = 21.7 raising the intriguing possibility that even (fainter) older quiescent UMGs could exist at this epoch.

The inner circumgalactic medium (CGM) of galaxies, where disk and halo processes intersect, remains poorly characterized despite its critical role in regulating galaxy evolution. We present results from Project AMIGA Insider, mapping Andromeda's (M31) inner CGM within 0.25 R_vir (~75 kpc) using 11 QSO sightlines, bringing our total sample to 54 sightlines from the disk to 2 R_vir. We detect a clear transition between M31's thick disk and CGM at R < 30 kpc, where low/intermediate ions show thick-disk corotating components with higher column densities than the CGM ones, while high ions exhibit similar column densities in both the CGM and thick disk. Beyond this region, all ion column densities decrease with impact parameter, with steeper gradients for low ions than high ions. The inner CGM (R < 100 kpc) shows more complex gas phases and multi-component absorption compared to the predominantly single-component outer CGM. We find no significant azimuthal dependence for any observed ions, suggesting M31's CGM is shaped by radial processes (e.g., cooling flows, precipitation) rather than disk-aligned outflows. We estimate the total metal mass in M31's cool (SiII, SiIII, SiIV) CGM within R_vir to be (1.9+/-0.3_stat+/-0.7_sys)x10^7 M_sun, leading to a cool gas mass of approximately 6x10^9 (Z/0.3 Z_sun)^-1 M_sun. The warmer OVI gas may contain at least 10 times more metal and gas mass. Compared to the COS-Halos L* galaxies, M31's cool CGM shows lower Si column densities at R < 0.4 R_200 and lower cool CGM masses, possibly resulting from M31's higher halo mass or different environments.

HST-10 is one of the largest proplyds in the Orion Nebula and is located approximately 1' SE of the Trapezium. Unlike other proplyds in Orion, however, the long-axis of HST-10 does not align with theta 1 C Ori, but is instead aligned with the rotational axis of the HST-10 disk. This cannot be easily explained using current photo-evaporation models. In this letter, we present high spatial resolution near-infrared images of the Orion proplyd HST-10 using Keck/NIRC2 with the Laser Guide Star Adaptive Optics system, along with multi-epoch analysis of HH objects near HST-10 using Hubble Space Telescope WFPC2 and ACS cameras. Our narrow-band near-IR images resolve the proplyd ionization front (IF) and circumstellar disk down to 23 AU at the distance to Orion in Br gamma, He I, H_2, and PAH emission. Br gamma and He I emission primarily trace the IF (with the disk showing prominently in silhouette), while the H_2 and PAH emission trace the surface of the disk itself. PAH emission also traces small dust grains within the proplyd envelope which is asymmetric and does not coincide with the IF. The curious morphology of the PAH emission may be due to UV-heating by both theta 1C Ori and theta 2A Ori. Multi-epoch HST images of the HST-10 field show proper motion of 3 knots associated with HH 517, clearly indicating that HST-10 has a jet. We postulate that the orientation of HST-10 is determined by the combined ram-pressure of this jet and the FUV-powered photo-ablation flow from the disk surface.

Reservoirs of dense atomic gas (primarily hydrogen), contain approximately 90 percent of the neutral gas at a redshift of 3, and contribute to 2-3 percent of the total baryons in the Universe. These damped Lyman-${\alpha}$ systems (so called because they absorb Lyman-${\alpha}$ photons from within and from background sources) have been studied for decades, but only through absorption lines present in the spectra of background quasars and gamma-ray bursts. Such pencil beams do not constrain the physical extent of the systems. Here, we report integral-field spectroscopy of a bright, gravitationally lensed galaxy at a redshift of 2.7 with two foreground damped Lyman-${\alpha}$ systems. These systems are &gt; 238 $kpc^2$ in extent, with column densities of neutral hydrogen varying by more than an order of magnitude on &lt; 3 kpc-scales. The mean column densities are $10^{20.46}$ - $10^{20.84} cm^{-2}$ and the total masses are &gt; 5.5 \times 10^{8} - $1.4 \times 10^{9} M_{\odot}$, showing that they contain the necessary fuel for the next generation of star formation, consistent with relatively massive, low-luminosity primeval galaxies at redshifts &gt; 2.

The prioritization and improvement of ethics, planetary protection, and safety standards in the astro-sciences is the most critical priority as our scientific and exploratory capabilities progress, both within government agencies and the private sector. These priorities lie in the belief that every single science mission - crewed or non-crewed, ground-based or not - should heed strict ethical and safety standards starting at the very beginning of a mission. Given the inevitability of the private sector in influencing future crewed missions both in and beyond low-Earth orbit, it is essential to the science community to agree on universal standards of safety, mission assurance, planetary protection, and especially anti-colonization. These issues will impact all areas of space science. Examples that are particularly relevant to the Astro2020 Decadal Survey include but are not limited to: light pollution from satellites, the voices and rights of Native people when constructing telescopes on their lands, and the need to be cognizant of contamination when searching for and exploring habitable environments beyond Earth. Ultimately, moving international space law and domestic space policy from a reactive nature to a proactive one will ensure the future of space exploration is one that is safe, transparent, and anti-imperialist.

We present spectroscopic confirmation of an ultra-massive galaxy (UMG) with $\log(M_\star/M_\odot) = 10.98 \pm 0.07$ at $z_\mathrm{spec} = 4.8947$ in the Extended Groth Strip (EGS), based on deep observations of Ly$\alpha$ emission with Keck/DEIMOS. The ultra-massive galaxy (UMG-28740) is the most massive member in one of the most significant overdensities in the EGS, with four additional photometric members with $\log(M_\star/M_\odot) > 10.5$ within $R_\mathrm{proj} \sim 1$ cMpc. The Ly$\alpha$ profile is highly asymmetric ($A_f = 3.56$), suggesting the presence of neutral gas within the interstellar medium, circumgalactic medium, or via AGN-driven outflows. Spectral energy distribution (SED) fitting using a large suite of star formation histories and two sets of high-quality photometry from ground- and space-based facilities consistently estimates the stellar mass of UMG-28740 to be $\log(M_\star/M_\odot) \sim 11$ with a small standard deviation between measurements ($\sigma = 0.07$). While the best-fit SED models agree on stellar mass, we find discrepancies in the estimated star formation rate for UMG-28740, resulting in either a star-forming or quiescent system. JWST/NIRCam photometry of UMG-28740 strongly favors a quiescent scenario, demonstrating the need for high-quality mid-IR observations. Assuming the galaxy to be quiescent, UMG-28740 formed the bulk of its stars at $z > 10$ and is quenching at $z \sim 8$, resulting in a high star formation efficiency at high redshift ($\epsilon \sim 0.2$ at $z \sim 5$ and $\epsilon \gtrsim 1$ at $z \gtrsim 8$). As the most massive galaxy in its protocluster environment, UMG-28740 is a unique example of the impossibly early galaxy problem.

The ASTrometric and phase-Referenced Astronomy (ASTRA) project will provide phase referencing and astrometric observations at the Keck Interferometer, leading to enhanced sensitivity and the ability to monitor orbits at an accuracy level of 30-100 microarcseconds. Here we discuss recent scientific results from ASTRA, and describe new scientific programs that will begin in 2010-2011. We begin with results from the "self phase referencing" (SPR) mode of ASTRA, which uses continuum light to correct atmospheric phase variations and produce a phase-stabilized channel for spectroscopy. We have observed a number of protoplanetary disks using SPR and a grism providing a spectral dispersion of ~2000. In our data we spatially resolve emission from dust as well as gas. Hydrogen line emission is spectrally resolved, allowing differential phase measurements across the emission line that constrain the relative centroids of different velocity components at the 10 microarcsecond level. In the upcoming year, we will begin dual-field phase referencing (DFPR) measurements of the Galactic Center and a number of exoplanet systems. These observations will, in part, serve as precursors to astrometric monitoring of stellar orbits in the Galactic Center and stellar wobbles of exoplanet host stars. We describe the design of several scientific investigations capitalizing on the upcoming phase-referencing and astrometric capabilities of ASTRA.

We present a comprehensive multi-epoch photometric and spectroscopic study of SN 2024bch, a nearby (19.9 Mpc) Type II supernova (SN) with prominent early high ionization emission lines. Optical spectra from 2.9 days after the estimated explosion reveal narrow lines of H I, He II, C IV, and N IV that disappear by day 6. High cadence photometry from the ground and TESS show that the SN brightened quickly and reached a peak M$_V \sim$ $-$17.8 mag within a week of explosion, and late-time photometry suggests a $^{56}$Ni mass of 0.050 M$_{\odot}$. High-resolution spectra from day 8 and 43 trace the unshocked circumstellar medium (CSM) and indicate a wind velocity of 30--40 km s$^{-1}$, a value consistent with a red supergiant (RSG) progenitor. Comparisons between models and the early spectra suggest a pre-SN mass-loss rate of $\dot{M} \sim 10^{-3}-10^{-2}\ M_\odot\ \mathrm{yr}^{-1}$, which is too high to be explained by quiescent mass loss from RSGs, but is consistent with some recent measurements of similar SNe. Persistent blueshifted H I and [O I] emission lines seen in the optical and NIR spectra could be produced by asymmetries in the SN ejecta, while the multi-component H$\alpha$ may indicate continued interaction with an asymmetric CSM well into the nebular phase. SN 2024bch provides another clue to the complex environments and mass-loss histories around massive stars.

We present deep imaging of Sirius B, the closest and brightest white dwarf, to constrain post-main-sequence planetary evolution in the Sirius system. We use Keck/NIRC2 in L'-band (3.776 $\mu$m) across three epochs in 2020 using the technique of angular differential imaging. Our observations are speckle-limited out to 1 AU and background-limited beyond. The 5$\sigma$ detection limits from our best performing epoch are 17 to 20.4 L' absolute magnitude. We consider multiple planetary formation pathways in the context of Sirius B's evolution to derive mass sensitivity limits, and achieve sub-Jupiter sensitivities at sub-AU separations, reaching 1.6 $\mathrm{M_J}$ to 2.4 $\mathrm{M_J}$ at 0.5 AU down to a sensitivity of 0.7 $\mathrm{M_J}$ to 1.2 $\mathrm{M_J}$ at >1 AU. Consistent with previous results, we do not detect any companions around Sirius B. Our strong detection limits demonstrate the potential of using high-contrast imaging to characterize nearby white dwarfs.

This study is a part of the Cosmicflows-4 project with the aim of measuring the distances of more than ~10,000 spiral galaxies in the local universe up to ~15,000 km/s. New HI linewidth information has come primarily from the Arecibo Legacy Fast ALFA Survey. Photometry of our sample galaxies has been carried out in optical (SDSS u, g, r, i, z) and infrared (WISE W1 and W2) bands. Inclinations have been determined using an online graphical interface accessible to a collaboration of citizen scientists. Galaxy distances are measured based on the correlation between the rotation rate of spirals and their absolute luminosity, known as the Tully-Fisher Relation (TFR). In this study, we present the calibration of the TFR using a subsample of ~600 spirals located in 20 galaxy clusters. Correlations among such observables as color, surface brightness, and relative HI content are explored in an attempt to reduce the scatter about the TFR with the goal of obtaining more accurate distances. A preliminary determination of the Hubble constant from the distances and velocities of the calibrator clusters is H0=76.0+-1.1(stat.)+-2.3(sys.) km/s/Mpc.

Advanced AO systems will likely utilise Pyramid wave-front sensors (PWFS) over the traditional Shack-Hartmann sensor in the quest for increased sensitivity, peak performance and ultimate contrast. Here, we wish to bring knowledge and quantify the PWFS theoretical limits as a means to highlight its properties and use cases. We explore forward models for the PWFS in the spatial-frequency domain for they prove quite useful since a) they emanate directly from physical-optics (Fourier) diffraction theory; b) provide a straightforward path to meaningful error breakdowns, c) allow for reconstruction algorithms with $O (n\,log(n))$ complexity for large-scale systems and d) tie in seamlessly with decoupled (distributed) optimal predictive dynamic control for performance and contrast optimisation. All these aspects are dealt with here. We focus on recent analytical PWFS developments and demonstrate the performance using both analytic and end-to-end simulations. We anchor our estimates with observed on-sky contrast on existing systems and then show very good agreement between analytical and Monte-Carlo estimates for the PWFS. For a potential upgrade of existing high-contrast imagers on 10\,m-class telescopes with visible or near-infrared PWFS, we show under median conditions at Paranal a contrast improvement (limited by chromatic and scintillation effects) of 2x-5x by replacing the wave-front sensor alone at large separations close to the AO control radius where aliasing dominates, and factors in excess of 10x by coupling distributed control with the PWFS over most of the AO control region, from small separations starting with the Inner Working Angle of typically 1-2 $\lambda/D$ to the AO correction edge (here 20 $\lambda/D$).

We present a census of neutral gas in the Milky Way disk and halo down to limiting column densities of $N$(HI)$\sim10^{14}$ cm$^{-2}$ using measurements of HI Lyman-series absorption from the Far Ultraviolet Spectroscopic Explorer (FUSE). Our results are drawn from an analysis of 25 AGN sightlines spread evenly across the sky with Galactic latitude |b|$\gtrsim 20^{\circ}$. By simultaneously fitting multi-component Voigt profiles to 11 Lyman-series absorption transitions covered by FUSE (Ly$\beta$-Ly$\mu$) plus HST measurements of Ly$\alpha$, we derive the kinematics and column densities of a sample of 152 HI absorption components. While saturation prevents accurate measurements of many components with column densities 17$\lesssim$log$N$(HI)$\lesssim$19, we derive robust measurements at log$N$(HI)$\lesssim$17 and log$N$(HI)$\gtrsim$19. We derive the first ultraviolet HI column density distribution function (CDDF) of the Milky Way, both globally and for low-velocity (ISM), intermediate-velocity clouds (IVCs), and high-velocity clouds (HVCs). We find that IVCs and HVCs show statistically indistinguishable CDDF slopes, with $\beta_{\rm IVC}=$ $-1.01_{-0.14}^{+0.15}$ and $\beta_{\rm HVC}=$ $-1.05_{-0.06}^{+0.07}$. Overall, the CDDF of the Galactic disk and halo appears shallower than that found by comparable extragalactic surveys, suggesting a relative abundance of high-column density gas in the Galactic halo. We derive the sky covering fractions as a function of HI column density, finding an enhancement of IVC gas in the northern hemisphere compared to the south. We also find evidence for an excess of inflowing HI over outflowing HI, with $-$0.88$\pm$0.40 M$_\odot$ yr$^{-1}$ of HVC inflow versus 0.20$\pm$0.10 M$_\odot$ yr$^{-1}$ of HVC outflow, confirming an excess of inflowing HVCs seen in UV metal lines.

Context. The recent estimates of the 3D shape of the M/Xe-type triple asteroid system (216) Kleopatra indicated a density of 5 g.cm$^{-3}$. Such a high density implies a high metal content and a low porosity which is not easy to reconcile with its peculiar dumbbell shape. Aims. Given the unprecedented angular resolution of the VLT/SPHERE/ZIMPOL camera, we aim to constrain the mass and the shape of Kleopatra with high accuracy, hence its density. Methods. We combined our new VLT/SPHERE observations of Kleopatra recorded in 2017 and 2018 with archival data, as well as lightcurve, occultation, and delay-Doppler images, to derive its 3D shape model using two different algorithms (ADAM, MPCD). Furthermore, an N-body dynamical model allowed us to retrieve the orbital elements of the two moons as explained in the accompanying paper. Results. The shape of Kleopatra is very close to an equilibrium dumbbell figure with two lobes and a thick neck. Its volume equivalent diameter (118.75$\pm$1.40) km and mass (2.97$\pm$0.32) 10$^{18}$ kg imply a bulk density of (3.38$\pm$0.50) g cm$^{-3}$. Such a low density for a supposedly metal-rich body indicates a substantial porosity within the primary. This porous structure along with its near-equilibrium shape is compatible with a formation scenario including a giant impact followed by reaccumulation. Kleopatra's current rotation period and dumbbell shape imply that it is in a critically rotating state. The low effective gravity along the equator of the body, together with the equatorial orbits of the moons and possibly rubble-pile structure, opens the possibility that the moons formed via mass shedding. Conclusions. Kleopatra is a puzzling multiple system due to the unique characteristics of the primary. It deserves particular attention in the future, with the Extremely Large Telescopes and possibly a dedicated space mission.