The Panchromatic Hubble Andromeda Southern Treasury (PHAST) is a large 195-orbit Hubble Space Telescope program imaging ∼0.45 deg2 of the southern half of M31's star-forming disk at optical and near-ultraviolet (NUV) wavelengths. The PHAST survey area extends the northern coverage of the Panchromatic Hubble Andromeda Treasury (PHAT) down to the southern half of M31, covering out to a radius of ∼13 kpc along the southern major axis and in total ∼two-thirds of M31's star-forming disk. This new legacy imaging yields stellar photometry of over 90 million resolved stars using the Advanced Camera for Surveys in the optical (F475W and F814W), and the Wide Field Camera 3 (WFC3) in the NUV (F275W and F336W). The photometry is derived using all overlapping exposures across all bands, and achieves a 50% completeness-limited depth of F475W ∼ 27.7 in the lowest surface density regions of the outer disk and F475W ∼ 26.0 in the most crowded, high surface brightness regions near M31's bulge. We provide extensive analysis of the data quality, including artificial star tests to quantify completeness, photometric uncertainties, and flux biases, all of which vary due to the background source density and the number of overlapping exposures. We also present seamless population maps of the entire M31 disk, which show relatively well-mixed distributions for stellar populations older than 1–2 Gyr, and highly structured distributions for younger populations. The combined PHAST + PHAT photometry catalog of ∼0.2 billion stars is the largest ever produced for equidistant sources and is available for public download by the community.
The American Astronomical Society (AAS), established in 1899 and based in Washington, DC, is the major organization of professional astronomers in North America. Its membership of about 7,000 individuals also includes physicists, mathematicians, geologists, engineers, and others whose research and educational interests lie within the broad spectrum of subjects comprising contemporary astronomy. The mission of the AAS is to enhance and share humanity's scientific understanding of the universe.
The Institute of Physics (IOP) is a leading scientific society promoting physics and bringing physicists together for the benefit of all. It has a worldwide membership of around 50 000 comprising physicists from all sectors, as well as those with an interest in physics. It works to advance physics research, application and education; and engages with policy makers and the public to develop awareness and understanding of physics. Its publishing company, IOP Publishing, is a world leader in professional scientific communications.
A publishing partnership
ISSN: 1538-4357
The Astrophysical Journal is an open access journal devoted to recent developments, discoveries, and theories in astronomy and astrophysics. Publications in ApJ constitute significant new research that is directly relevant to astrophysical applications, whether based on observational results or on theoretical insights or modeling.
Remembering Judy Pipher (1940–2022)
GOLD OPEN ACCESS FROM 1 JANUARY 2022
Zhuo Chen et al 2025 ApJ 979 35
Thomas Williams et al 2020 ApJ 892 134
Following the success of the first mission, the High-Resolution Coronal Imager (Hi-C) was launched for a third time (Hi-C 2.1) on 2018 May 29 from the White Sands Missile Range, NM, USA. On this occasion, 329 s of 17.2 nm data of target active region AR 12712 were captured with a cadence of ≈4 s, and a plate scale of 0129 pixel−1. Using data captured by Hi-C 2.1 and co-aligned observations from SDO/AIA 17.1 nm, we investigate the widths of 49 coronal strands. We search for evidence of substructure within the strands that is not detected by AIA, and further consider whether these strands are fully resolved by Hi-C 2.1. With the aid of multi-scale Gaussian normalization, strands from a region of low emission that can only be visualized against the contrast of the darker, underlying moss are studied. A comparison is made between these low-emission strands and those from regions of higher emission within the target active region. It is found that Hi-C 2.1 can resolve individual strands as small as ≈202 km, though the more typical strand widths seen are ≈513 km. For coronal strands within the region of low emission, the most likely width is significantly narrower than the high-emission strands at ≈388 km. This places the low-emission coronal strands beneath the resolving capabilities of SDO/AIA, highlighting the need for a permanent solar observatory with the resolving power of Hi-C.
Ryo Tazaki et al 2025 ApJ 980 49
We present near- and mid-infrared (IR) broadband imaging observations of the edge-on protoplanetary disk around HH 30 with the James Webb Space Telescope/Near Infrared Camera and the Mid-Infrared Instrument (MIRI). We combine these observations with archival optical/near-IR scattered light images obtained with the Hubble Space Telescope and a millimeter-wavelength dust continuum image obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) with the highest spatial resolution ever obtained for this target. Our multiwavelength images clearly reveal the vertical and radial segregation of micron-sized and submillimeter-sized grains in the disk. In the near- and mid-IR, the images capture not only bireflection nebulae separated by a dark lane but also diverse dynamical processes occurring in the HH 30 disk, such as spiral- and tail-like structures, a conical outflow, and a collimated jet. In contrast, the ALMA image reveals a flat dust disk in the disk midplane. By performing radiative transfer simulations, we show that grains of about 3 μm in radius or larger are fully vertically mixed to explain the observed mid-IR scattered light flux and its morphology, whereas millimeter-sized grains are settled into a layer with a scale height of ≳1 au at 100 au from the central star. We also find a tension in the disk inclination angle inferred from optical/near-IR and millimeter observations, with the latter being closer to exactly edge-on. Finally, we report the first detection of the proper motion of an emission knot associated with the mid-IR collimated jet detected by combining two epochs of our MIRI 12.8 μm observations.
Peter G. Boorman et al 2025 ApJ 978 118
Hard X-ray-selected samples of active galactic nuclei (AGN) provide one of the cleanest views of supermassive black hole accretion but are biased against objects obscured by Compton-thick gas column densities of NH > 1024 cm−2. To tackle this issue, we present the NuSTAR Local AGN NH Distribution Survey (NuLANDS)—a legacy sample of 122 nearby (z < 0.044) AGN primarily selected to have warm infrared colors from IRAS between 25 and 60 μm. We show that optically classified Type 1 and 2 AGN in NuLANDS are indistinguishable in terms of optical [O iii] line flux and mid-to-far-infrared AGN continuum bolometric indicators, as expected from an isotropically selected AGN sample, while Type 2 AGN are deficient in terms of their observed hard X-ray flux. By testing many X-ray spectroscopic models, we show the measured line-of-sight column density varies on average by ∼1.4 orders of magnitude depending on the obscurer geometry. To circumvent such issues, we propagate the uncertainties per source into the parent column density distribution, finding a directly measured Compton-thick fraction of 35% ± 9%. By construction, our sample will miss sources affected by severe narrow-line reddening, and thus segregates sources dominated by small-scale nuclear obscuration from large-scale host-galaxy obscuration. This bias implies an even higher intrinsic obscured AGN fraction may be possible, although tests for additional biases arising from our infrared selection find no strong effects on the measured column density distribution. NuLANDS thus holds potential as an optimized sample for future follow-up with current and next-generation instruments aiming to study the local AGN population in an isotropic manner.
Aomawa L. Shields et al 2025 ApJ 979 45
Discoveries of giant planet candidates orbiting white dwarf (WD) stars and the demonstrated capabilities of the James Webb Space Telescope bring the possibility of detecting rocky planets in the habitable zones (HZs) of WDs into pertinent focus. We present simulations of an aqua planet with an Earth-like atmospheric composition and incident stellar insolation orbiting in the HZ of two different types of stars—a 5000 K WD and main-sequence K-dwarf star Kepler-62 (K62) with a similar effective temperature—and identify the mechanisms responsible for the two differing planetary climates. The synchronously rotating WD planet's global mean surface temperature is 25 K higher than that of the synchronously rotating planet orbiting K62, due to its much faster (10 hr) rotation and orbital period. This ultrafast rotation generates strong zonal winds and meridional flux of zonal momentum, stretching out and homogenizing the scale of atmospheric circulation, and preventing an equivalent buildup of thick, liquid water clouds on the dayside of the planet compared to the synchronous planet orbiting K62, while also transporting heat equatorward from higher latitudes. White dwarfs may therefore present amenable environments for life on planets formed within or migrated to their HZs, generating warmer surface environments than those of planets with main-sequence hosts to compensate for an ever shrinking incident stellar flux.
Qinyue Fei et al 2025 ApJ 980 84
We conduct a study of the gas kinematics of two quasar host galaxies at z ≳ 6 traced by the [C ii] emission line using the Atacama Large Millimeter/submillimeter Array. By combining deep observations at both low and high resolution, we recover the diffuse emission, resolve its structure, and measure the rotation curves from the inner region of the galaxy to its outskirts using DysmalPy and 3DBarolo. Assuming that both galaxies exhibit disk rotation driven by the gravitational potential of the galaxy, we find that the best-fit disk models have a Vrot/σ ≈ 2 and inferred circular velocities out to ∼6–8 kpc scales, well beyond the likely stellar distribution. We then determine the mass profiles of each component (stars, gas, dark matter) with priors on the baryon and dark matter properties. We find relatively large dark matter fractions within their effective radii (fDM(R < Re) = and
, respectively), which are significantly larger than those extrapolated from lower redshift studies and remain robust under different input parameters verified by Monte Carlo simulations. The large fDM(R < Re) corresponds to halo masses of ∼1012.5−1012.8 M⊙, thus representative of the most massive halos at these redshifts. Notably, while the masses of these supermassive black holes (SMBHs) are approximately 1 dex higher than the low-redshift relationship with stellar mass, the closer alignment of SMBH and halo masses with a local relationship may indicate that the early formation of these SMBHs is linked to their dark matter halos, providing insights into the coevolution of galaxies and black holes in the early Universe.
Thomas M. Do et al 2025 ApJ 979 50
Current multi-spacecraft in situ measurements allow for the investigation of the time evolution of energetic particles at interplanetary shocks (IPs) at small (≲0.1 au) heliocentric distances. The energy spectrum of accelerated particles at IPs was shown by a previous 1D transport model that includes both self-excited plus preexisting turbulence and a term representing the escape of particles from the system to gradually steepen as a result of a finite acceleration-to-escape timescales ratio; such a model was found in excellent agreement with the entire sample of the ground-level enhancement spectra of solar cycle 23. We solve the time-dependent case of such a model in the case of diffusion dominated by preexisting turbulence. The average timescale for particle acceleration at various heliocentric distances, from 1 au down to the inner heliosphere (<0.1 au), is shorter than in the no-escape case, as higher energy particles have a shorter time to accelerate before completely leaving the system into the upstream medium. A simple scaling with time of the time-dependent spectrum is provided. We compare the "nose" structure at a few ∼100s keV protons first measured in situ by Parker Solar Probe in crossing the very fast 2022 September 5 shock at 0.07 au; we find that the nose is reasonably well explained by a lack of the highest energy particles not yet produced by the young shock by both our model and the no-escape version.
Adam G. Riess et al 2024 ApJ 977 120
We cross-check the Hubble Space Telescope (HST) Cepheid/Type Ia supernova (SN Ia) distance ladder, which yields the most precise local H0, against early James Webb Space Telescope (JWST) subsamples (∼1/4 of the HST sample) from SH0ES and CCHP, calibrated only with NGC 4258. We find HST Cepheid distances agree well (∼1σ) with all combinations of methods, samples, and telescopes. The comparisons explicitly include the measurement uncertainty of each method in NGC 4258, an oft-neglected but dominant term. Mean differences are ∼0.03 mag, far smaller than the 0.18 mag "Hubble tension." Combining all measures produces the strongest constraint yet on the linearity of HST Cepheid distances, 0.994 ±0.010, ruling out distance-dependent bias or offset as the source of the tension at ∼7σ. However, current JWST subsamples produce large sampling differences in H0 whose size and direction we can directly estimate from the full HST set. We show that ΔH0 ∼ 2.5 km s−1 Mpc−1 between the CCHP JWST program and the full HST sample is entirely consistent with differences in sample selection. We combine all JWST samples into a new distance-limited set of 16 SNe Ia at D ≤ 25 Mpc. Using JWST Cepheids, JAGB, and tip of the red giant branch, we find 73.4 ± 2.1, 72.2 ± 2.2, and 72.1 ± 2.2 km s−1 Mpc−1, respectively. Explicitly accounting for common supernovae, the three-method JWST result is H0 = 72.6 ± 2.0, similar to H0 = 72.8 expected from HST Cepheids in the same galaxies. The small JWST sample trivially lowers the Hubble tension significance due to small-sample statistics and is not yet competitive with the HST set (42 SNe Ia and 4 anchors), which yields 73.2 ± 0.9. Still, the joint JWST sample provides important cross-checks that the HST data pass.
Hisashi Hayakawa et al 2025 ApJ 979 49
In 2024 May, the scientific community observed intense solar eruptions that resulted in a great geomagnetic storm and auroral extensions, highlighting the need to document and quantify these events. This study mainly focuses on their quantification. The source active region (AR; NOAA Active Region 13664) evolved from 113 to 2761 millionths of the solar hemisphere between May 4 and 14. NOAA AR 13664's magnetic free energy surpassed 1033 erg on May 7, triggering 12 X-class flares on May 8–15. Multiple interplanetary coronal mass ejections (ICMEs) were produced from this AR, accelerating solar energetic particles toward Earth. According to satellite and interplanetary scintillation data, at least four ICMEs erupted from AR 13664, eventually overcoming and combining each other. The shock arrival at 17:05 UT on May 10 significantly compressed the magnetosphere down to ≈5.04 RE and triggered a deep Forbush Decrease. GOES satellite data and ground-based neutron monitors confirmed a ground-level enhancement from 2 UT to 10 UT on 2024 May 11. The ICMEs induced exceptional geomagnetic storms, peaking at a provisional Dst index of −412 nT at 2 UT on May 11, marking the sixth-largest storm since 1957. The AE and AL indices showed great auroral extensions that located the AE/AL stations into the polar cap. We gathered auroral records at that time and reconstructed the equatorward boundary of the visual auroral oval to 298 invariant latitude. We compared naked-eye and camera auroral visibility, providing critical caveats on their difference. We also confirmed global disturbances of the storm-enhanced density of the ionosphere.
M. G. Aartsen et al 2017 ApJ 835 151
Since the recent detection of an astrophysical flux of high-energy neutrinos, the question of its origin has not yet fully been answered. Much of what is known about this flux comes from a small event sample of high neutrino purity, good energy resolution, but large angular uncertainties. In searches for point-like sources, on the other hand, the best performance is given by using large statistics and good angular reconstructions. Track-like muon events produced in neutrino interactions satisfy these requirements. We present here the results of searches for point-like sources with neutrinos using data acquired by the IceCube detector over 7 yr from 2008 to 2015. The discovery potential of the analysis in the northern sky is now significantly below = 10−12 TeV cm−2 s−1, on average 38% lower than the sensitivity of the previously published analysis of 4 yr exposure. No significant clustering of neutrinos above background expectation was observed, and implications for prominent neutrino source candidates are discussed.
Pinjian Chen et al 2025 ApJ 980 227
A consensus has been reached in recent years that binarity plays an important role in the formation and evolution of a significant fraction of planetary nebulae (PNe). Utilizing the archived photometric data from the Zwicky Transient Facility survey, we conducted a comprehensive data mining in search of brightness variations in a large sample of Galactic PNe. This effort led to identification of 39 PNe, whose central stars exhibit periodic variation in light curves. Among these objects, 20 are known binary central stars of PNe, while the remaining 19 are new discoveries. Additionally, we identified 14 PNe with central stars displaying anomalous variation in light curves, as well as eight variables based on the high-cadence photometric data. Among the new discoveries of periodicity, we found compelling evidence of binary systems at the centers of two archetypal quadrupolar PNe. We also report on very peculiar brightness variation observed in the central core of the compact PN NGC 6833. Several PNe in our sample deserve follow-up observations, both high-dispersion spectroscopy and high-precision photometry, to reveal the true nature of their central binarity or even multiplicity.
Emilio Tejeda and Jesús A. Toalá 2025 ApJ 980 226
The Bondi–Hoyle–Lyttleton (BHL) accretion model is widely used to describe how a compact object accretes material from a companion's stellar wind in binary systems. However, its standard implementation becomes inaccurate when the wind velocity (vw) is comparable to or less than the orbital velocity (vo), predicting nonphysical accretion efficiencies above unity. This limits its applicability to systems with low wind-to-orbital velocity ratios (w = vw/vo ≤ 1), such as symbiotic systems. We revisit the implementation of the BHL model and introduce a geometric correction factor that accounts for the varying orientation of the accretion cylinder relative to the wind direction. This correction ensures physically plausible accretion efficiencies (η ≤ 1) for all w in circular orbits. Our new implementation naturally predicts the flattening of the accretion efficiency observed in numerical simulations for w < 1, without the need for ad hoc adjustments. We also peer into the implications of our prescription for the less-explored case of eccentric orbits, highlighting the key role of the geometric correction factor in shaping the accretion process. We compare our predictions with numerical simulations, finding good agreement for a wide range of parameters. Applications to the symbiotic star R Aqr and the X-ray binary LS 5039 are presented. This improved implementation offers a more accurate description of wind accretion in binary systems, with implications for stellar evolution, population synthesis, and observational data interpretation.
Michael W. Topping et al 2025 ApJ 980 225
JWST has recently discovered a subset of reionization era galaxies with ionized gas that is metal-poor in oxygen and carbon but heavily enriched in nitrogen. This abundance pattern is almost never seen in lower-redshift galaxies but is commonly observed in globular cluster stars. We have recently demonstrated that this peculiar abundance pattern appears in a compact (≃20 pc) metal-poor galaxy undergoing a strong burst of star formation. This galaxy was originally selected based on strong C iv emission, indicating a hard radiation field rarely seen locally. In this paper, we present JWST/NIRSpec observations of another reionization-era galaxy known to power strong C iv emission, the z = 7.04 gravitationally lensed galaxy A1703-zd6. The emission-line spectrum reveals this is a metal-poor galaxy () dominated by a young stellar population (
Myr) that powers a very hard ionizing spectrum (C iv equivalent width, EW = 19.4 Å, He ii EW = 2.2 Å). The interstellar medium is highly enriched in nitrogen (
) with very high electron densities (8–19 × 104 cm−3) and extreme ionization conditions rarely seen at lower redshift. We also find intense CIV emission (EW ≳ 20 Å) in two new z ≳ 6 metal-poor galaxies. To put these results in context, we search for UV line emission in a sample of 737 z ≳ 4 galaxies with NIRSpec spectra, establishing that 40%(30%) of systems with [O iii]+Hβ EW > 2000 Å have N iv] (C iv) detections with EW > 5 Å(> 10 Å). These results suggest high N/O ratios, and hard ionizing sources appear in a brief phase following a burst of star formation in compact high-density stellar complexes.
Irin Babu Vathachira et al 2025 ApJ 980 224
We define two regimes of the parameter space of symbiotic systems based on the dominant mass transfer mechanism. A wide range of white dwarf (WD) mass, donor mass, and donor radius combinations are explored to determine the separation for each parameter combination, below which wind Roche-lobe (RL) overflow will be the dominant mass transfer mechanism. The underlying concept is the premise that the wind accelerates. If it reaches the RL before attaining sufficient velocity to escape, it will be trapped and gravitationally focused through the inner Lagrangian point toward the accreting WD. However, if the wind succeeds in attaining the required velocity to escape from the donor's RL, it will disperse isotropically, and the dominant mass transfer mechanism will be the Bondi–Hoyle–Lyttleton prescription, in which only a fraction of the wind will be accreted onto the WD. We present these two regimes of the 4D parameter space, covering 375 different parameter combinations.
Yuxuan Pang et al 2025 ApJ 980 223
The wide survey of the Chinese Space Station Telescope (CSST) will observe a large field of 17,500 deg2. The GU, GV, and GI grism observations of CSST will cover the wavelength range from 2550 to 10000 Å at a resolution of R ∼ 200 and a depth of about 22 AB magnitude for the continuum. In this paper, we present a pipeline to identify quasars and measure their physical properties with the CSST mock data. We simulate the raw images and extract the one-dimensional grism spectra for quasars, galaxies, and stars with r-band magnitudes of 18 < mr < 22 using the CSST Cycle 6 simulation code (https://csst-tb.bao.ac.cn/code/csst_sim/csst-simulation). Using a convolution neural network, we separate quasars from stars and galaxies. We measure the redshifts by identifying strong emission lines of quasars. We also fit the 1D slitless spectra with QSOFITMORE to estimate black hole masses and Eddington ratios. Our results show that CSST slitless spectroscopy can effectively separate quasars with redshifts z = 0−5 from other types of objects with an accuracy of 99%. Among those successfully classified quasars, 90% of them could have precise redshift measurements with σNMAD = 0.002. The scatters of black hole masses and Eddington ratios from the spectral fittings are 0.13 and 0.15 dex, respectively. The metallicity diagnosis line ratios have a scatter of 0.1–0.2 dex. Our results show that the CSST slitless spectroscopy survey has the potential to discover about 0.9 million new quasars and provide important contributions to the science of active galactic nuclei and cosmology.