Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt–Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We describe the discovery of a solar neighborhood (d = 468 pc) binary system with a main-sequence sunlike star and a massive noninteracting black hole candidate. The spectral energy distribution of the visible star is described by a single stellar model. We derive stellar parameters from a high signal-to-noise Magellan/MIKE spectrum, classifying the star as a main-sequence star with T _eff = 5972 K, log g = 4.54 , and M = 0.91 M _⊙. The spectrum shows no indication of a second luminous component. To determine the spectroscopic orbit of the binary, we measured the radial velocities of this system with the Automated Planet Finder, Magellan, and Keck over four months. We show that the velocity data are consistent with the Gaia astrometric orbit and provide independent evidence for a massive dark companion. From a combined fit of our spectroscopic data and the astrometry, we derive a companion mass of 11.39 − 1.31 + 1.51 M _⊙. We conclude that this binary system harbors a massive black hole on an eccentric (e = 0.46 ± 0.02), 185.4 ± 0.1 day orbit. These conclusions are independent of El-Badry et al., who recently reported the discovery of the same system. A joint fit to all available data yields a comparable period solution but a lower companion mass of 9.32 − 0.21 + 0.22 M ⊙ . Radial velocity fits to all available data produce a unimodal solution for the period that is not possible with either data set alone. The combination of both data sets yields the most accurate orbit currently available.

We use compiled high-precision pulsar timing measurements to directly measure the Galactic acceleration of binary pulsars relative to the solar system barycenter. Given the vertical accelerations, we use the Poisson equation to derive the Oort limit, i.e., the total volume mass density in the Galactic mid-plane. Our best-fitting model gives an Oort limit of 0.08 − 0.02 0.05 M ⊙ pc − 3 , which is close to estimates from recent Jeans analyses. Given the accounting of the baryon budget from McKee et al., we obtain a local dark matter density of − 0.004 − 0.02 0.05 M ⊙ pc − 3 , which is slightly below other modern estimates but consistent within the current uncertainties of our method. The error bars are currently about five times larger than kinematical estimates, but should improve in the future for this novel dynamical method. We also constrain the oblateness of the potential, finding it consistent with that expected from the disk and inconsistent with a potential dominated by a spherical halo, as is appropriate for our sample that is within a ∼kpc of the Sun. We find that current measurements of binary pulsar accelerations lead to large uncertainties in the slope of the rotation curve. We give a fitting function for the vertical acceleration a _z : a _z = − α ₁ z; log 10 ( α 1 / Gyr − 2 ) = 3.69 − 0.12 0.19 . By analyzing interacting simulations of the Milky Way, we find that large asymmetric variations in da _z /dz as a function of vertical height may be a signature of sub-structure. We end by discussing the power of combining constraints from pulsar timing and high-precision radial velocity measurements toward lines-of-sight near pulsars, to test theories of gravity and constrain dark matter sub-structure.

A bright (m _F150W,AB = 24 mag), z = 1.95 supernova (SN) candidate was discovered in JWST/NIRCam imaging acquired on 2023 November 17. The SN is quintuply imaged as a result of strong gravitational lensing by a foreground galaxy cluster, detected in three locations, and remarkably is the second lensed SN found in the same host galaxy. The previous lensed SN was called “Requiem,” and therefore the new SN is named “Encore.” This makes the MACS J0138.0−2155 cluster the first known system to produce more than one multiply imaged SN. Moreover, both SN Requiem and SN Encore are Type Ia SNe (SNe Ia), making this the most distant case of a galaxy hosting two SNe Ia. Using parametric host fitting, we determine the probability of detecting two SNe Ia in this host galaxy over a ∼10 yr window to be ≈3%. These observations have the potential to yield a Hubble constant (H ₀) measurement with ∼10% precision, only the third lensed SN capable of such a result, using the three visible images of the SN. Both SN Requiem and SN Encore have a fourth image that is expected to appear within a few years of ∼2030, providing an unprecedented baseline for time-delay cosmography.

Supernovae (SNe) that have been multiply imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply imaged Type Ia SN 2022qmx (aka, “SN Zwicky”) at z = 0.3544. SN Zwicky was discovered by the Zwicky Transient Facility in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multicycle “LensWatch (www.lenswatch.org)” program. We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, and F814W) but unresolved with WFC3/IR F160W, and present an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between lens-model-predicted time delays (≲1 day), and delays estimated with the single epoch of Hubble Space Telescope colors (≲3.5 days), including the uncertainty from chromatic microlensing (∼1–1.5 days). Our lens models converge to an Einstein radius of θ E = ( 0.168 − 0.005 + 0.009 ) ″ , the smallest yet seen in a lensed SN system. The “standard candle” nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter than predicted by ∼ 1.7 − 0.6 + 0.8 mag and ∼ 0.9 − 0.6 + 0.8 mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models.

Supernova (SN) SN H0pe is a gravitationally lensed, triply imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply imaged SNe provide a rare opportunity to constrain the Hubble constant (H ₀), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is located at z = 1.783 and is the first SN Ia with sufficient light-curve sampling and long enough time delays for an H ₀ inference. Here we present photometric time-delay measurements and SN properties of SN H0pe. Using JWST/NIRCam photometry, we measure time delays of Δt _ab = − 116.6 − 9.3 + 10.8 observer-frame days and Δt _cb = − 48.6 − 4.0 + 3.6 observer-frame days relative to the last image to arrive (image 2b; all uncertainties are 1σ), which corresponds to a ∼5.6% uncertainty contribution for H ₀ assuming 70 km s⁻¹ Mpc⁻¹. We also constrain the absolute magnification of each image to μ _a = 4.3 − 1.8 + 1.6 , μ _b = 7.6 − 2.6 + 3.6 , μ _c = 6.4 − 1.5 + 1.6 by comparing the observed peak near-IR magnitude of SN H0pe to the nonlensed population of SNe Ia.