We present the discovery of Andromeda XXXV, the faintest Andromeda satellite galaxy discovered to date, identified as an overdensity of stars in the Pan-Andromeda Archaeological Survey and confirmed via Hubble Space Telescope imaging. Located at a heliocentric distance of $92{7}_{-63}^{+76}$ kpc and $15{8}_{-45}^{+57}$ kpc from Andromeda, Andromeda XXXV is an extended (${r}_{h}=5{3\,}_{-11}^{+13}$ pc), elliptical ( = 0.4 ± 0.2), metal-poor ([Fe/H] ∼ −1.9) system, and the least luminous (M_V = −5.2 ± 0.3) of Andromeda's dwarf satellites discovered so far. Andromeda XXXV's properties are consistent with the known population of dwarf galaxies around the Local Group, bearing close structural resemblance to the Canes Venatici II and Hydra II Milky Way (MW) dwarf satellite galaxies. Its stellar population, characterized by a red horizontal branch or a red clump feature, mirrors that of other Andromeda satellite galaxies in showing evidence for a spread in age and metallicity, with no signs of younger stellar generations. This age–metallicity spread is not observed in MW satellites of comparable stellar mass, and highlights the persistent differences between the satellite systems of Andromeda and the MW, extending even into the ultrafaint regime.

The search for habitable environments and biomarkers in exoplanetary atmospheres is the holy grail of exoplanet science. The detection of atmospheric signatures of habitable Earth-like exoplanets is challenging owing to their small planet–star size contrast and thin atmospheres with high mean molecular weight. Recently, a new class of habitable exoplanets, called Hycean worlds, has been proposed, defined as temperate ocean-covered worlds with H₂-rich atmospheres. Their large sizes and extended atmospheres, compared to rocky planets of the same mass, make Hycean worlds significantly more accessible to atmospheric spectroscopy with JWST. Here we report a transmission spectrum of the candidate Hycean world K2-18 b, observed with the JWST NIRISS and NIRSpec instruments in the 0.9–5.2 μm range. The spectrum reveals strong detections of methane (CH₄) and carbon dioxide (CO₂) at 5σ and 3σ confidence, respectively, with high volume mixing ratios of ∼1% each in a H₂-rich atmosphere. The abundant CH₄ and CO₂, along with the nondetection of ammonia (NH₃), are consistent with chemical predictions for an ocean under a temperate H₂-rich atmosphere on K2-18 b. The spectrum also suggests potential signs of dimethyl sulfide (DMS), which has been predicted to be an observable biomarker in Hycean worlds, motivating considerations of possible biological activity on the planet. The detection of CH₄resolves the long-standing missing methane problem for temperate exoplanets and the degeneracy in the atmospheric composition of K2-18 b from previous observations. We discuss possible implications of the findings, open questions, and future observations to explore this new regime in the search for life elsewhere.

Barnard's Star is an old, single M dwarf star that comprises the second-closest extrasolar system. It has a long history of claimed planet detections from both radial velocities and astrometry. However, none of these claimed detections have so far withstood further scrutiny. Continuing this story, extreme precision radial velocity measurements from the ESPRESSO instrument have recently been used to identify four new sub-Earth-mass planet candidates around Barnard's Star. We present here 112 radial velocities of Barnard's Star from the MAROON-X instrument that were obtained independently to search for planets around this compelling object. The data have a typical precision of 30 cm s⁻¹ and are contemporaneous with the published ESPRESSO measurements (2021–2023). The MAROON-X data on their own confirm planet b (P = 3.154 days) and planet candidates c and d (P = 4.124 and 2.340 days, respectively). Furthermore, adding the MAROON-X data to the ESPRESSO data strengthens the evidence for planet candidate e (P = 6.739 days), thus leading to its confirmation. The signals from all four planets are <50 cm s⁻¹, the minimum masses of the planets range from 0.19 to 0.34 M_⊕, and the system is among the most compact known among late M dwarfs hosting low-mass planets. The current data rule out planets with masses >0.57 M_⊕ (with a 99% detection probability) in Barnard's Star's habitable zone (P = 10–42 days).

We present the first results of JWST Cycle 1 and 2 observations of Sgr A* using NIRCam taken simultaneously at 2.1 and 4.8 μm for a total of ~48 hr over seven different epochs in 2023 and 2024. We find correlated variability at 2.1 and 4.8 μm in all epochs, continual short-timescale (a few seconds) variability, and epoch-to-epoch variable emission implying long-term (~days to months) variability of Sgr A*. A highlight of this analysis is the evidence for subminute, horizon-scale time variability of Sgr A*, probing inner accretion disk size scales. The power spectra of the light curves in each observing epoch also indicate long-term variable emission. With continuous observations, JWST data suggest that the flux of Sgr A* is fluctuating constantly. The flux density correlation exhibits a distinct break in the slope at ~3 mJy at 2.1 μm. The analysis indicates two different processes contributing to the variability of Sgr A*. Brighter emission trends toward shallower spectral indices than the fainter emission. Cross-correlation of the light curves indicates for the first time a time delay of 3–40 s in the 4.8 μm variability with respect to 2.1 μm. This phase shift leads to loops in plots of flux density versus spectral index as the emission rises and falls. Modeling suggests that the synchrotron emission from the evolving, age-stratified electron population reproduces the shape of the observed light curves with a direct estimate of the magnetic field strengths in the range between 40 and 90 G and an upper cutoff energy, E_c, between 420 and 720 MeV.

Isolated planetary-mass objects share their mass range with planets but do not orbit a star. They lack the necessary mass to support fusion in their cores and thermally radiate their heat from formation as they cool, primarily at infrared wavelengths. Many isolated planetary-mass objects show variations in their infrared brightness consistent with nonuniform atmospheric features modulated by their rotation. SIMP J013656.5+093347.3 is a rapidly rotating isolated planetary-mass object, and previous infrared monitoring suggests complex atmospheric features rotating in and out of view. The physical nature of these features is not well understood, with clouds, temperature variations, thermochemical instabilities, and infrared-emitting aurora all proposed as contributing mechanisms. Here we report JWST time-resolved low-resolution spectroscopy from 0.8 to 11 μm of SIMP J013656.5+093347.3, which supports the presence of three specific features in the atmosphere: clouds, hot spots, and changing carbon chemistry. We show that no single mechanism can explain the variations in the time-resolved spectra. When combined with previous studies of this object indicating patchy clouds and aurorae, these measurements reveal the rich complexity of the atmosphere of SIMP J013656.5+093347.3. Gas giant planets in the solar system, specifically Jupiter and Saturn, also have multiple cloud layers and high-altitude hot spots, suggesting these phenomena are also present in worlds both within and beyond our solar system.

The Earth sits inside a 300 pc-wide void that was carved by a series of supernova explosions that went off tens of millions of years ago, pushing away interstellar gas and creating a bubble-like structure. The ⁶⁰Fe peak deposits found in the deep-sea crust have been interpreted by the imprints left by the ejecta of supernova explosions occurring about 2–3 and 5–6 Myr ago. It is likely that the ⁶⁰Fe peak at about 2–3 Myr originated from a supernova occurring in the Upper Centaurus Lupus association in Scorpius Centaurus (≈140 pc) or the Tucana-Horologium association (≈70 pc), whereas the  ≈5–6 Myr peak is likely attributed to the solar system's entrance into the bubble. In this Letter, we show that the supernova source responsible for synthesizing the ⁶⁰Fe peak deposits  ≈2–3 Myr ago can consistently explain the cosmic-ray spectrum and the large-scale anisotropy between 100 TeV and 100 PeV. The cosmic-ray knee could then potentially be attributed entirely to a single nearby "PeVatron" source. Matching the intensity and shape of the cosmic-ray spectrum allows us to place stringent constraints on the cosmic-ray energy content from the supernova as well as on the cosmic-ray diffusion coefficient. Making use of such constraints, we provide a robust estimate of the temporal variation of terrestrial ionizing cosmic radiation levels and discuss their implications in the development of early life on Earth by plausibly influencing the mutation rate and, as such, conceivably assisting in the evolution of complex organisms.

We present high-definition observations with the James Webb Space Telescope (JWST) of >1000 Cepheids in a geometric anchor of the distance ladder, NGC 4258, and in five hosts of eight Type Ia supernovae, a far greater sample than previous studies with JWST. These galaxies individually contain the largest samples of Cepheids, an average of >150 each, producing the strongest statistical comparison to those previously measured with the Hubble Space Telescope (HST) in the near-infrared (NIR). They also span the distance range of those used to determine the Hubble constant with HST, allowing us to search for a distance-dependent bias in HST measurements. The superior resolution of JWST negates crowding noise, the largest source of variance in the NIR Cepheid period–luminosity relations (Leavitt laws) measured with HST. Together with the use of two epochs to constrain Cepheid phases and three filters to remove reddening, we reduce the dispersion in the Cepheid P–L relations by a factor of 2.5. We find no significant difference in the mean distance measurements determined from HST and JWST, with a formal difference of −0.01 ± 0.03 mag. This result is independent of zero-points and analysis variants including metallicity dependence, local crowding, choice of filters, and slope of the relations. We can reject the hypothesis of unrecognized crowding of Cepheid photometry from HST that grows with distance as the cause of the "Hubble tension" at 8.2σ, i.e., greater confidence than that of the Hubble tension itself. We conclude that errors in photometric measurements of Cepheids across the distance ladder do not significantly contribute to the tension.

We have examined the distribution of the position angle (PA) of the Galactic center filaments with lengths L > 66'' and <66'' as well as their length distribution as a function of PA. We find bimodal PA distributions of the filaments, and long and short populations of radio filaments. Our PA study shows the evidence for a distinct population of short filaments with PA close to the Galactic plane. Mainly thermal, short-radio filaments (<66'') have PAs concentrated close to the Galactic plane within 60° < PA < 120°. Remarkably, the short filament PAs are radial with respect to the Galactic center at l < 0° and extend in the direction toward Sgr A*. On a smaller scale, the prominent Sgr E H ii complex G358.7-0.0 provides a vivid example of the nearly radial distribution of short filaments. The bimodal PA distribution suggests a different origin for two distinct filament populations. We argue that the alignment of the short-filament population results from the ram pressure of a degree-scale outflow from Sgr A* that exceeds the internal filament pressure, and aligns them along the Galactic plane. The ram pressure is estimated to be 2 × 10⁶ cm⁻³ K at a distance of 300 pc, requiring biconical mass outflow rate 10⁻⁴M_⊙ yr⁻¹ with an opening angle of ∼40°. This outflow aligns not only the magnetized filaments along the Galactic plane but also accelerates thermal material associated with embedded or partially embedded clouds. This places an estimate of ∼6 Myr as the age of the outflow.

When surrounded by a transparent emission region, black holes are expected to reveal a dark shadow caused by gravitational light bending and photon capture at the event horizon. To image and study this phenomenon, we have assembled the Event Horizon Telescope, a global very long baseline interferometry array observing at a wavelength of 1.3 mm. This allows us to reconstruct event-horizon-scale images of the supermassive black hole candidate in the center of the giant elliptical galaxy M87. We have resolved the central compact radio source as an asymmetric bright emission ring with a diameter of 42 ± 3 μas, which is circular and encompasses a central depression in brightness with a flux ratio ≳10:1. The emission ring is recovered using different calibration and imaging schemes, with its diameter and width remaining stable over four different observations carried out in different days. Overall, the observed image is consistent with expectations for the shadow of a Kerr black hole as predicted by general relativity. The asymmetry in brightness in the ring can be explained in terms of relativistic beaming of the emission from a plasma rotating close to the speed of light around a black hole. We compare our images to an extensive library of ray-traced general-relativistic magnetohydrodynamic simulations of black holes and derive a central mass of M = (6.5 ± 0.7) × 10⁹ M_⊙. Our radio-wave observations thus provide powerful evidence for the presence of supermassive black holes in centers of galaxies and as the central engines of active galactic nuclei. They also present a new tool to explore gravity in its most extreme limit and on a mass scale that was so far not accessible.

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of $\sim 1.7\,{\rm{s}}$ with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg² at a luminosity distance of ${40}_{-8}^{+8}$ Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 $\,{M}_{\odot }$. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at $\sim 40\,{\rm{Mpc}}$) less than 11 hours after the merger by the One-Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ∼10 days. Following early non-detections, X-ray and radio emission were discovered at the transient's position $\sim 9$ and $\sim 16$ days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC 4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta.

Spectral energy distribution (SED) modeling is one of the main methods to estimate galaxy properties, such as photometric redshifts, z, and stellar masses, M_*, for extragalactic imaging surveys. SEDs are currently modeled as light from a composite stellar population attenuated by a geometrically homogeneous foreground dust screen. This is despite evidence from simulations and observations that find large spatial variations in dust attenuation due to the detailed geometry of stars and gas within galaxies. In this work, we examine the impact of this simplistic dust assumption on inferred z and M_* for Rubin LSST. We first construct synthetic LSST-like observations (ugrizy magnitudes) from the Numerical Investigation of Hundred Astrophysical Objects (NIHAO)-SKIRT catalog, which provides SEDs from high-resolution hydrodynamic simulations using 3D Monte Carlo radiative transfer. We then infer z and M_* from the synthetic observations using the PROVABGS Bayesian SED modeling framework. Overall, the homogeneous dust screen assumption biases both z and M_* in galaxies, consistently and significantly for galaxies with dust attenuation A_V ≳ 0.5, and likely below. The biases depend on the orientation in which the galaxies are observed. At z = 0.4, z is overestimated by  ∼0.02 for face-on galaxies and M_* is underestimated by  ∼0.4 dex for edge-on galaxies. The bias in photo-z is equivalent to the desired redshift precision level of the LSST "gold sample" and will be larger at higher redshifts. Our results underscore the need for SED models with additional flexibility in the dust parameterization to mitigate significant systematic biases in cosmological analyses with LSST.

The Pleiades is a young open cluster that has not yet dynamically relaxed, making it an ideal target to observe various internal dynamical effects. By employing a well-defined sample of main-sequence (MS) cluster members, including both MS single stars and unresolved MS+MS binaries, we revisited their individual masses and mass functions and quantified the mass dependence of their radial distributions. We found that the mass function of binaries is more top-heavy than that of single stars. Significant mass segregation is observed for both single and binary populations, respectively, with more massive objects concentrated toward the cluster center. Notably, within given mass ranges, binaries are distributed more scattered than single stars, providing direct evidence for more efficient dynamical disruption of binaries in the inner region. The radial distribution of the binary fraction, expressed as the f_b–R relation can be characterized by a bimodal shape, with higher f_b values in both innermost and outermost regions of the cluster. The lower-mass subsample exhibits a monotonic increase in f_b with radius, reflecting the impact of binary disruption. Conversely, for the higher-mass subsample, f_b decreases with radius. It can be explained that these massive cluster members, which possess higher binary probabilities, have already undergone significant mass segregation. All these observational evidence and analyses related to the radial mass distribution imply that the Pleiades is currently undergoing a complicated interplay of various internal dynamical effects, of which the modulation between mass segregation and binary disruption is particularly pronounced.

We present deep optical observations of the stellar halo of NGC 300, an LMC-mass galaxy, acquired with the DEEP subcomponent of the DECam Local Volume Exploration survey using the 4 m Blanco Telescope. Our resolved star analysis reveals a large, low surface brightness stellar stream (M_V ∼ −8.5; [Fe/H] = −1.4 ± 0.15) extending more than 40 kpc north from the galaxy's center. We also find other halo structures, including potentially an additional stream wrap to the south, which may be associated with the main stream. The morphology and derived low metallicities of the streams and shells discovered surrounding NGC 300 are highly suggestive of a past accretion event. Assuming a single progenitor, the accreted system is approximately Fornax-like in luminosity, with an inferred mass ratio to NGC 300 of approximately 1:15. We also present the discovery of a metal-poor globular cluster (GC) (R_proj = 23.3 kpc; M_V = −8.99 ± 0.16; [Fe/H] ≈ −1.6 ± 0.6) in the halo of NGC 300, the furthest identified GC associated with NGC 300. The stellar structures around NGC 300 represent the richest features observed in a Magellanic Cloud analog to date, strongly supporting the idea that accretion and subsequent disruption is an important mechanism in the assembly of dwarf galaxy stellar halos.

GJ 436b is the archetype warm Neptune exoplanet. The planet's thermal emission spectrum was previously observed via intensive secondary eclipse campaigns with Spitzer. The atmosphere has long been interpreted to be extremely metal-rich, out of chemical equilibrium, and potentially tidally heated. We present the first panchromatic emission spectrum of GJ 436b observed with JWST's NIRCAM (F322W2 and F444W) and MIRI (LRS) instruments between 2.4 and 11.9 μm. Surprisingly, the JWST spectrum appears significantly fainter around 3.6 μm than that implied by Spitzer photometry. The molecular absorption features in the spectrum are relatively weak, and we only find tentative evidence of CO₂ absorption at 2σ. Under the assumption of a dayside blackbody, we find T_day = 662.8 ± 5.0 K, which is similar to the zero Bond albedo equilibrium temperature. We use it to obtain a 3σ upper limit on the Bond albedo of A_B ≤ 0.66. To understand the spectrum, we employ 1D radiative–convective models but find that atmospheric constraints depend strongly on model assumptions. If thermochemical equilibrium is assumed, we find a cloudy metal-enriched atmosphere (metallicity ≥300× solar). We employ 1D photochemical modeling to show that the observed spectrum is also consistent with a cloud-free, relatively lower metallicity atmosphere (metallicity ≥80× solar) with a cold internal temperature (T_int ∼ 60 K). These are much lower metallicities and internal temperatures than inferences from Spitzer photometry. The low T_day and nondetection of transmission features at high spectral resolution do suggest a role for cloud opacity, but this is not definitive.

The study of the atmosphere of exoplanets orbiting white dwarfs is a largely unexplored field. With WD 0806-661 b, we present the first deep dive into the atmospheric physics and chemistry of a cold exoplanet around a white dwarf. We observed WD 0806-661 b using JWST's Mid-InfraRed Instrument Low-Resolution Spectrometer, covering the wavelength range from 5 to 12 μm, and the Imager, providing us with 12.8, 15, 18, and 21 μm photometric measurements. We carried the data reduction of those data sets, tackling second-order effects to ensure a reliable retrieval analysis. Using the TauREx retrieval code, we inferred the pressure–temperature structure, atmospheric chemistry, mass, and radius of the planet. The spectrum of WD 0806-661 b is shaped by molecular absorption of water, ammonia, and methane, consistent with a cold Jupiter atmosphere, allowing us to retrieve their abundances. From the mixing ratio of water, ammonia, and methane we derive C/O = 0.34 ± 0.06, ${\rm{C}}/{\rm{N}}=14.{4}_{-1.8}^{+2.5}$, and N/O = 0.023 ± 0.004 and the ratio of detected metals as a proxy for metallicity. We also derive upper limits for the abundance of CO and CO₂ (1.2 × 10⁻⁶ and 1.6 × 10⁻⁷, respectively), which were not detected by our retrieval models. While our interpretation of WD 0806-661 b's atmosphere is mostly consistent with our theoretical understanding, some results—such as the lack of evidence for water clouds, an apparent increase in the mixing ratio of ammonia at low pressure, or the retrieved mass at odds with the supposed age—remain surprising and require follow-up observational and theoretical studies to be confirmed.