Table of contents for issue 1, volume 947, The Astrophysical Journal

The cosmic dust particles found in space are mainly porous aggregates of smaller grains. Theoretically, these aggregates are replicated using fractal geometry, assuming a cluster of spheres. Although the light scattering response of cosmic dust aggregates has been thoroughly studied using clusters of spherical grains in the past few decades, the effect of irregularities on the surface of each grain in an entire aggregate has mostly been neglected. We introduce, for the first time, a visually realistic cosmic dust model that incorporates a mixture of rough fractal aggregates (RFA) and agglomerated debris (Solids) to replicate the unusual polarization–phase curve observed in the case of the interstellar comet 2I/Borisov at multiple wavelengths. The authenticity of the RFA structures has been verified by replicating light scattering results of circumstellar dust analogs from the Granada Amsterdam Light Scattering Database. We demonstrate that the light scattering response from the RFA structures has a very close resemblance to the experimental values. Finally, we model the observed polarization–phase curve of the interstellar comet 2I/Borisov using a mixture of RFA and solid particles. The best-fit data indicate the presence of a higher percentage of porous RFA structures (80%) owing to the fact that the comet carries a higher percentage of small and highly porous pristine cosmic dust particles. Further, the model indicates that the unusually steep polarimetric slope and the high dust-to-gas ratio in newer comets are mainly due to a higher porous-to-compact ratio.

The AGN STORM 2 Collaboration targeted the Seyfert 1 galaxy Mrk 817 for a year-long multiwavelength, coordinated reverberation mapping campaign including Hubble Space Telescope, Swift, XMM-Newton, NICER, and ground-based observatories. Early observations with NICER and XMM revealed an X-ray state 10 times fainter than historical observations, consistent with the presence of a new dust-free, ionized obscurer. The following analysis of NICER spectra attributes variability in the observed X-ray flux to changes in both the column density of the obscurer by at least one order of magnitude (N_H ranges from to ) and the intrinsic continuum brightness (the unobscured flux ranges from 10^−11.8 to 10^−10.5 erg s⁻¹ cm⁻²). While the X-ray flux generally remains in a faint state, there is one large flare during which Mrk 817 returns to its historical mean flux. The obscuring gas is still present at lower column density during the flare, but it also becomes highly ionized, increasing its transparency. Correlation between the column density of the X-ray obscurer and the strength of UV broad absorption lines suggests that the X-ray and UV continua are both affected by the same obscuration, consistent with a clumpy disk wind launched from the inner broad-line region.

We report the detection of gravity (g) and Rossby (r) mode pulsations of four short-period eclipsing binaries, KIC 5439790, KIC 7501230, KIC 9350889, and KIC 9453192 based on the 4 yr Kepler high-precision photometry. Light-curve modeling reveals that the four binaries are all contact systems with small mass ratios. We study the short-term variability of the light-curve residuals after the removal of the binary model and attribute them to the g- and r-mode pulsations that stem from the primary stars. By introducing a new criterion, we attempt to identify period-spacing patterns in the Fourier spectra, which refers to the determination of the rotation rates of stellar interior and envelope as well as the asymptotic period spacings of the pulsating stars. It is interesting to find that the rotation rates of the stellar envelopes are all nearly equal to the orbits. The near-core rotations, however, are significantly slower by about 10%. Based on the derived asymptotic period spacings, stellar parameters were constrained from asteroseismology models. The pulsators are revealed to be evolved main-sequence stars with high metallicities.

Coalescence of multiple magnetic islands is recognized as an effective process to energize particles during magnetic reconnection, while its energy conversion process still remains unclear. Here, a two-dimensional fully kinetic simulation of multiple island coalescence with a small reconnection guide field is studied. In the analysis of energy conversion within a magnetic island, the dot product of is a useful quantity to compare with j · E = w₂, since the average work done by the Lorentz force on the circulating particles is negligible in the island and . A bipolar pattern of w₁ is found at a secondary island when the electrons are in circular motion inside the island. Significant energy dynamo (w₃ < 0) resulting from j_∥E_∥ is found at the secondary island, which has not been reported before, where the parallel electric field E_∥ is highly correlated with w₃. Moreover, significant energy dissipation (w₃ > 0) due to is seen in the merging region between two coalescing islands. Both types of energy conversions are accompanied by enhancements in j_∥ and the parallel electron temperature T_e∥. Three ion-scale magnetic islands (FR1, FR2, and FR3) observed by the Magnetospheric Multiscale spacecraft are compared favorably with the simulated signatures of energy dynamo and dissipation of an evolving secondary island. In particular, FR1 displayed a similar energy dynamo signature as that simulated in an early stage of the secondary island. FR2 and FR3 showed a dominant energy conversion similar to that obtained in a later stage of the secondary island.

I report on validation and testing of a novel 3D reconstruction method than can obtain coronal plasma properties from a single snapshot perspective. I first reported on the method in 2021, and I have since named it the Coronal Reconstruction Onto B-Aligned Regions (CROBAR) method. The testing and validation are carried out with a cube from a MURaM 3D MHD simulation, which affords a coronal-like “ground truth” against which the reconstruction method can be applied and compared. I find that the method does quite well, recovering the “coronal veil”−like features recently reported from the MURaM simulations and allaying concerns that these features would thwart recovery of valid 3D coronal structure from a limited number of perspectives. I also find that a second perspective between ∼45° and 90° does significantly improve the reconstructions. Two distinct channels with soft-X-ray-like temperature response (peaking above 5 MK) would suffice for CROBAR’s optically thin observables A suite of AIA-like EUV passbads, with good coverage in the 3–8 MK range, is also well suited to CROBAR.

Constraining L dwarf properties from their spectra is challenging. Near-infrared (NIR) spectra probe a limited range of pressures, while many species condense within their photospheres. Condensation creates two complexities: gas-phase species “rain out” (decreasing in abundances by many orders of magnitude) and clouds form. We designed tests using synthetic data to determine the best approach for retrieving L dwarf spectra, isolating the challenges in the absence of cloud opacity. We conducted atmospheric retrievals on synthetic cloud-free L dwarf spectra derived from the Sonora Bobcat models at SpeX resolution using a variety of thermal and chemical abundance profile parameterizations. For objects hotter than L5 (T_eff ∼ 1700 K), the limited pressure layers probed in the NIR are mostly convective; parameterized pressure–temperature (PT) profiles bias results and free, unsmoothed profiles should be used. Only when many layers both above and below the radiative-convective boundary are probed can parameterized profiles provide accurate results. Furthermore, a nonuniform abundance profile for FeH is needed to accurately retrieve bulk properties of early-to-mid L dwarfs. Nonuniform prescriptions for other gases in NIR retrievals may also be warranted near the L/T transition (CH₄) and early Y dwarfs (Na and K). We demonstrate the utility of using realistic, self-consistent models to benchmark retrievals and suggest how they can be used in the future.

The stellar cluster environment is expected to play a central role in the evolution of circumstellar disks. We use thermochemical modeling to constrain the dust and gas masses, disk sizes, UV and X-ray radiation fields, viewing geometries, and central stellar masses of 20 class II disks in the Orion Nebula Cluster (ONC). We fit a large grid of disk models to 350 GHz continuum, CO J = 3 − 2, and HCO⁺J = 4 − 3 Atacama Large Millimeter/submillimeter Array observations of each target, and we introduce a procedure for modeling interferometric observations of gas disks detected in absorption against a bright molecular cloud background. We find that the ONC disks are massive and compact, with typical radii <100 au, gas masses ≥10⁻³M_⊙, and gas-to-dust ratios ≥100. The interstellar‐medium‐like gas-to-dust ratios derived from our modeling suggest that compact, externally irradiated disks in the ONC are less prone to gas-phase CO depletion than the massive and extended gas disks that are commonly found in nearby low-mass star-forming regions. The presence of massive gas disks indicates that external photoevaporation may have only recently begun operating in the ONC; though it remains unclear whether other cluster members are older and more evaporated than the ones in our sample. Finally, we compare our dynamically derived stellar masses with the stellar masses predicted from evolutionary models and find excellent agreement. Our study has significantly increased the number of dynamical mass measurements in the mass range ≤0.5 M_⊙, demonstrating that the ONC is an ideal region for obtaining large samples of dynamical mass measurements toward low-mass M-dwarfs.

In dense stellar clusters like galactic nuclei and globular clusters, stellar densities are so high that stars might physically collide with each other. In galactic nuclei the energy and power output can be close to, and even exceed, those from supernovae events. We address the event rate and the electromagnetic characteristics of collisions of main-sequence stars (MS) and red giants (RGs). We also investigate the case in which the cores form a binary and emit gravitational waves. In the case of RGs, this is particularly interesting because the cores are degenerate. We find that MS event rate can be as high as tens per year, and that of RGs 1 order of magnitude larger. The collisions are powerful enough to mimic supernovae or tidal disruptions events. We find Zwicky Transient Facility observational data that seem to exhibit the features we describe. The cores embedded in the gaseous debris experience a friction force that has an impact on the chirping mass of the gravitational wave. As a consequence, the two small cores in principle mimic two supermassive black holes merging. However, their evolution in frequency along with the precedent electromagnetic burst and the ulterior afterglow are efficient tools to reveal the impostors. In the particular case of RGs, we derive the properties of the degenerate He cores and their H-burning shells to analyze the formation of the binaries. The merger is such that it can be misclassified with SN Ia events. Because the masses and densities of the cores are so dissimilar in values depending on their evolutionary stage, the argument about standard candles and cosmic ladder should be reevaluated.

The LIGO/Virgo detected a gravitational wave (GW) event, named GW200224_222234 (also known as S200224ca) and classified as a binary-black hole coalescence, on 2020 February 24. Given its relatively small localization skymap (71 deg² for a 90% credible region; revised to 50 deg² in GWTC-3), we performed target-of-opportunity observations using the Subaru/Hyper Suprime-Cam (HSC) in the r2 and z bands. Observations were conducted on 2020 February 25 and 28 and March 23, with the first epoch beginning 12.3 hr after the GW detection. The survey covered the highest-probability sky area of 56.6 deg², corresponding to a 91% probability. This was the first deep follow-up (m_r ≳ 24, m_z ≳ 23) for a binary-black hole merger covering >90% of the localization. By performing image subtraction and candidate screening including light-curve fitting with transient templates and examples, we found 22 off-nucleus transients that were not ruled out as the counterparts of GW200224_222234 with our Subaru/HSC data alone. We also performed GTC/OSIRIS spectroscopy of the probable host galaxies for five candidates; two are likely to be located within the 3D skymap, whereas the others are not. In conclusion, 19 transients remain as possible optical counterparts of GW200224_222234; but we could not identify a unique promising counterpart. If there are no counterparts in the remaining candidates, the upper limits of the optical luminosity are erg s⁻¹ and erg s⁻¹ in the r2 and z bands, respectively, at ∼12 hr after GW detection. We also discuss improvements in the strategies of optical follow-ups for future GW events.

New analyses of gravitational-wave events raise questions about the nature of some events. For example, LIGO–Virgo–KAGRA initially determined GW151226 to be a merger with a mass ratio q ≈ 0.5 and effective inspiral spin χ_eff ≈ 0.2. However, recent works offer an alternative picture: GW151226 is a lower mass ratio event q ≈ 0.3 with slightly higher spin χ_eff ≈ 0.3. This discrepancy has been challenging to resolve, as a wide range of differences are employed for each analysis. This work introduces a “deep follow-up” framework to efficiently compute the posterior odds between two different peaks in parameter space. In doing so, we aim to help resolve disputes about the true nature of gravitational-wave events associated with conflicting astrophysical interpretations. Our proposal is not a replacement for standard inference techniques; instead, our method provides a diagnostic tool to understand discrepancies between conflicting results. We demonstrate this method by studying three q–χ_eff peaks proposed for GW151226. We find that the (q ∼ 0.5, χ_eff ∼ 0.2) interpretation is only slightly preferred over the (q ∼ 0.3, χ_eff ∼ 0.3) hypothesis with posterior odds of ∼1.7 ± 0.4, suggesting that neither of the two peaks can be ruled out. We discuss strategies to produce more reliable parameter estimation studies in gravitational-wave astronomy.

Astrophysical jets are present in a range of environments, including young stellar objects, X-ray binaries, and active galactic nuclei, but their formation is still not fully understood. As one of the nearest symbiotic binary stars, R Aquarii (D ∼ 220 pc) offers a unique opportunity to study the inner region within ∼600 au of the jet source, which is particularly crucial to our understanding of nonrelativistic jet formation and origin. We present high-angular resolution UV and optical imaging from the Hubble Space Telescope in six emission-line regions of the inner jet. Using these observations to obtain a range of representative line ratios for our system and kinematic data derived from a comparison with previous studies, we model the shocked gas in order to determine the relative roles of shock heating and photoionization in the R Aquarii system. We find that our shock models suggest that a nonzero magnetic field is needed to describe the measured line ratios. We also find that the Mg iiλλ2795,2802 intensities are overpredicted by our models for most of the jet regions, perhaps because of depletion onto grains or to opacity in these resonance lines.

Close encounters between stars in star-forming regions are important as they can perturb or destroy protoplanetary disks, young planetary systems, and multiple-star systems. We simulate simple, virialized, equal-mass N-body star clusters and find that both the rate and the total number of encounters between stars vary by factors of several in statistically identical clusters due to the stochastic/chaotic details of the orbits and stellar dynamics. Encounters tend to “saturate” rapidly in the core of a cluster, with stars there each having many encounters, while more distant stars have none. However, we find that the fraction of stars that has had at least one encounter within a particular distance grows in the same way (scaling with the crossing time and half-mass radius) in all clusters, and we present a new (empirical) way of estimating the fraction of stars that has had at least one encounter at a particular distance.

We combine an unprecedented MaNGA sample of over 3000 passive galaxies in the stellar mass range 10⁹–10¹²M_⊙ with the Sloan Digital Sky Survey group catalog by Tinker to quantify how central and satellite formation, quantified by radial profiles in stellar age, [Fe/H], and [Mg/Fe], depends on the stellar mass of the galaxy (M_*) and the mass of the host halo (M_h). After controlling for M_* and M_h, the stacked spectra of centrals and satellites beyond the effective radius (r_e) show small, yet significant differences in multiple spectral features at the 1% level. According to spectral fitting with the code alf, a primary driver of these differences appears to be [Mg/Fe] variations, suggesting that stellar populations in the outskirts of satellites formed more rapidly than the outer populations of centrals. To probe the physical mechanisms that may be responsible for this signal, we examined how satellite stellar populations depend on M_h. We find that satellites in high-M_h halos show older stellar ages, lower [Fe/H], and higher [Mg/Fe] compared to satellites in low-M_h halos, especially for M_* = 10^9.5–10^10.5M_⊙. These signals lend support to environmentally driven processes that quench satellite galaxies, although variations in the merger histories of central and satellite galaxies also emerge as a viable explanation.

We study the gravitational lensing of acoustic charged black holes in strong and weak field limit approximations. For this purpose, we first numerically obtain the deflection limit coefficients and deflection angle in the strong field limit. We observe that the strong deflection angle α_D increases with increasing magnitude of the charged parameter Q and that the strong deflection angle α_D of an acoustic charged black hole with tuning parameter ξ = 4 is greater than that of a standard Reissner–Nordström black hole (ξ = 0). We also study the astrophysical consequences via strong gravitational lensing by taking the example of various supermassive black holes in the center of several galaxies and observe that the acoustic charged black hole could be quantitatively distinguished from standard Reissner–Nordström (ξ = 0) and standard Schwarzschild (ξ = 0, Q = 0) black holes. Furthermore, by using the Gauss–Bonnet theorem, we derive the weak deflection angle in the background of an acoustic charged black hole in the curved spacetime. We find that, for fixed values of the charged parameter Q and the tuning parameter (ξ = 0 or 4), the weak deflection angle σ_D decreases with the impact parameter b. We also observe that the weak deflection angle σ_D decreases with increasing magnitude of the charged parameter Q for a fixed value of the tuning parameter (ξ = 0 or 4). Our results suggest that the observational test for an acoustic charged black hole is indeed feasible, and it is generalized to the cases of acoustic Schwarzschild (Q = 0), standard Reissner–Nordström (ξ = 0), and standard Schwarzschild (ξ = 0, Q = 0) black holes.

Among the several candidate models for chondrule formation, the lighting model has been recognized to be less likely than the other two major models, shock-wave heating and planetesimal collision. It might be because we have believed that the lightning model predicts cooling rates of chondrules that are too fast to reproduce their textures with the assumption that the discharge channels must be optically thin. However, the previous works revealed that the buildup of a strong electric field to generate the lightning in protoplanetary disks requires the enhancement of the solid density. Moreover, some properties of chondrules indicate their formation in environments with such a high solid density. Therefore, the discharge channels may be optically thick, and the lightning model can potentially predict the proper cooling rates of chondrules. In this study, we reinvestigate the cooling rates of chondrules produced by the lightning in the solid-rich environments considering the radiative transfer and the expansion of the hot channel. Chondrules must interact dynamically with the surrounding gas and dust via the drag force. We consider two limiting cases for the dynamics of chondrules: the drag force is ignored in the first case, and chondrules are completely coupled with their surroundings in the second case. In both cases, the lightning model predicts the proper cooling rates of chondrules under the optically thick conditions with high solid enhancement. Therefore, the lightning model is worth further investigation to judge its reliability as the source of chondrule formation.

We present a detection of 21 cm emission from large-scale structure (LSS) between redshift 0.78 and 1.43 made with the Canadian Hydrogen Intensity Mapping Experiment. Radio observations acquired over 102 nights are used to construct maps that are foreground filtered and stacked on the angular and spectral locations of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasars (QSOs) from the eBOSS clustering catalogs. We find decisive evidence for a detection when stacking on all three tracers of LSS, with the logarithm of the Bayes factor equal to 18.9 (LRG), 10.8 (ELG), and 56.3 (QSO). An alternative frequentist interpretation, based on the likelihood ratio test, yields a detection significance of 7.1σ (LRG), 5.7σ (ELG), and 11.1σ (QSO). These are the first 21 cm intensity mapping measurements made with an interferometer. We constrain the effective clustering amplitude of neutral hydrogen (H i), defined as , where Ω_{H i} is the cosmic abundance of H i, b_{H i} is the linear bias of H i, and 〈fμ²〉 = 0.552 encodes the effect of redshift-space distortions at linear order. We find for LRGs (z = 0.84), for ELGs (z = 0.96), and for QSOs (z = 1.20), with constraints limited by modeling uncertainties at nonlinear scales. We are also sensitive to bias in the spectroscopic redshifts of each tracer, and we find a nonzero bias Δ v = − 66 ± 20 km s⁻¹ for the QSOs. We split the QSO catalog into three redshift bins and have a decisive detection in each, with the upper bin at z = 1.30 producing the highest-redshift 21 cm intensity mapping measurement thus far.

We explore how the fraction of quenched galaxies changes in groups of galaxies with respect to the distance to the center of the group, redshift, and stellar mass to determine the dominant process of environmental quenching in 0.2 < z < 0.8 groups. We use new UV data from the UVCANDELS project in addition to existing multiband photometry to derive new galaxy physical properties of the group galaxies from the zCOSMOS 20 k group catalog. Limiting our analysis to a complete sample of log (M_*/M_⊙) > 10.56 group galaxies, we find that the probability of being quenched increases slowly with decreasing redshift, diverging from the stagnant field galaxy population. A corresponding analysis on how the probability of being quenched increases with time within groups suggests that the dominant environmental quenching process is characterized by slow (∼Gyr) timescales. We find a quenching time of approximately Gyr, consistent with the slow processes of strangulation and delayed-then-rapid quenching although more data are needed to confirm this result.

We analyzed state-of-the-art observations of the solar atmosphere to investigate the dependence of the Ca ii K brightness of several solar features on spectral bandwidth and spatial resolution of the data. In particular, we study data obtained at the Swedish Solar Telescope with the Crisp Imaging Spectropolarimeter and Chromospheric Imaging Spectrometer instruments. The analyzed data, which are characterized by a spectral bandwidth of 0.12 Å and a spatial resolution of 0078, were acquired close to the disk center by targeting a quiet-Sun area and an active region. We convolved the original observations with Gaussian kernels to degrade their spectral bandwidth and spatial resolution to the instrumental characteristics of the most prominent series of Ca ii K observations available to date. We then studied the effect of data degradation on the observed regions and on parameters derived from Ca ii K line measurements that are largely employed as diagnostics of the solar and stellar chromospheres. We find that the effect of degrading the spectral resolution of Ca ii K observations and line profiles depends on both the employed bandwidth and observed solar region. Besides, we found that the spatial degradation impacts the data characterized by a broad bandwidth to a larger extent compared to those acquired with a narrow band. However, the appearance of the observed solar regions is only slightly affected by the spatial resolution of data with bandwidths up to 1 Å and in the range [3,10] Å. Finally, we derived relationships that can be used to intercalibrate results from observations taken with different instruments in diverse regions of the solar atmosphere.

Based on independent shear measurements using the Dark Energy Camera Legacy Survey/DR8 imaging data, we measure the weak lensing signals around isolated central galaxies (ICGs) from Sloan Digital Sky Survey/DR7 at z ∼ 0.1. The projected stellar mass density profiles of satellite galaxies are further deduced, using photometric sources from the Hyper Suprime-cam survey. The signals of ICGs + their extended stellar halos are taken from Wang et al. All measurements are compared with predictions by the IllustrisTNG300-1 simulation. We find, overall, a good agreement between observation and TNG300. In particular, a correction to the stellar mass of massive observed ICGs is applied based on the calibration of He et al., which brings a much better agreement with TNG300 predicted lensing signals at . In real observation, red ICGs are hosted by more massive dark matter halos and have more satellites and more extended stellar halos than blue ICGs at fixed stellar mass. However, in TNG300 there are more satellites around blue ICGs at fixed stellar mass, and the outer stellar halos of red and blue ICGs are similar. The stellar halos of TNG galaxies are more extended compared with real observed galaxies, especially for blue ICGs with . We find the same trend for TNG100 galaxies and for true halo central galaxies. The tensions between TNG and real galaxies indicate that satellite disruptions are stronger in TNG. In both TNG300 and observation, satellites approximately trace the underlying dark matter distribution beyond 0.1R₂₀₀, but the fraction of total stellar mass in TNG300 does not show the same radial distribution as real galaxies.

We present the results of a systematic search for candidate quiescent galaxies in the distant universe in 11 JWST fields with publicly available observations collected during the first 3 months of operations and covering an effective sky area of ∼145 arcmin². We homogeneously reduce the new JWST data and combine them with existing observations from the Hubble Space Telescope. We select a robust sample of ∼80 candidate quiescent and quenching galaxies at 3 < z < 5 using two methods: (1) based on their rest-frame UVJ colors, and (2) a novel quantitative approach based on Gaussian mixture modeling of the near-UV − U, U − V, and V − J rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive (M_⋆ ≥ 10^10.6 M_⊙) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of 2–3, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at z > 3, as could be expected for highly clustered massive systems. Importantly, JWST enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available.

It is not yet settled how the combination of secular processes and merging gives rise to the bulges and pseudobulges of galaxies. The nearby (D ∼ 4.2 Mpc) disk galaxy M94 (NGC 4736) has the largest pseudobulge in the local universe, and offers a unique opportunity for investigating the role of merging in the formation of its pseudobulge. We present a first ever look at M94's stellar halo, which we expect to contain a fossil record of M94's past mergers. Using Subaru's Hyper Suprime-Cam, we resolve and identify red giant branch (RGB) stars in M94's halo, finding two distinct populations. After correcting for completeness through artificial star tests, we can measure the radial profile of each RGB population. The metal-rich RGB stars show an unbroken exponential profile to a radius of 30 kpc that is a clear continuation of M94's outer disk. M94's metal-poor stellar halo is detectable over a wider area and clearly separates from its metal-rich disk. By integrating the halo density profile, we infer a total accreted stellar mass of ∼2.8 × 10⁸ M_⊙, with a median metallicity of [M/H] = −1.4. This indicates that M94's most-massive past merger was with a galaxy similar to, or less massive than, the Small Magellanic Cloud. Few nearby galaxies have had such a low-mass dominant merger; therefore we suggest that M94's pseudobulge was not significantly impacted by merging.

While solar-like oscillations in red giants have been observed at massive scales by the Kepler mission, few features of these oscillation mode frequencies, other than their global properties, have been exploited for stellar characterization. The signatures of acoustic glitches in mode frequencies have been used for studying main-sequence stars, but the validity of applying such techniques to evolved red giants, particularly pertaining to the inclusion of nonradial modes, has been less well examined. Making use of new theoretical developments, we characterize glitches using the π modes associated with red giant stellar models, and use our procedure to examine for the first time how the properties of the He ii acoustic glitch—specifically its amplitude and associated acoustic depth—vary over the course of evolution up the red giant branch, and with respect to other fundamental stellar properties. We find that the acoustic depths of these glitches, in conjunction with other spectroscopic information, discriminate between red giants in the first-ascent and core-helium-burning phases. We critically reexamine previous attempts to constrain acoustic glitches from nonradial (in particular dipole) modes in red giants. Finally, we apply our fitting procedure to Kepler data, to evaluate its robustness to noise and other observational systematics.

We construct a sample of 644 carbon-enhanced metal-poor (CEMP) stars with abundance analyses based on moderate- to high-resolution spectroscopic studies. Dynamical parameters for these stars are estimated based on radial velocities, Bayesian parallax-based distance estimates, and proper motions from Gaia EDR3 and DR3, supplemented by additional available information where needed. After separating our sample into the different CEMP morphological groups in the Yoon–Beers diagram of absolute carbon abundance versus metallicity, we used the derived specific energies and actions (E, J_r, J_ϕ, J_z) to cluster them into Chemodynamically Tagged Groups (CDTGs). We then analyzed the elemental-abundance dispersions within these clusters by comparing them to the dispersion of clusters that were generated at random. We find that, for the Group I (primarily CEMP-s and CEMP-r/s) clustered stars, there exist statistically insignificant intracluster dispersions in [Fe/H], [C/Fe]_c (evolution corrected carbon), and [Mg/Fe] when compared to the intracluster dispersions of randomly clustered Group I CEMP stars. In contrast, the Group II (primarily CEMP-no) stars exhibit clear similarities in their intracluster abundances, with very low, statistically significant, dispersions in [C/Fe]_c and marginally significant results in [Mg/Fe]. These results strongly indicate that Group I CEMP stars received their carbon enhancements from local phenomena, such as mass transfer from an evolved binary companion in regions with extended star formation histories, while the CDTGs of Group II CEMP stars formed in low-metallicity environments that had already been enriched in carbon, likely from massive rapidly rotating ultra- and hyper-metal-poor stars and/or supernovae associated with high-mass early-generation stars.

Microlensing events have historically been discovered throughout the Galactic bulge and plane by surveys designed solely for that purpose. We conduct the first multiyear search for microlensing events on the Zwicky Transient Facility (ZTF), an all-sky optical synoptic survey that observes the entire visible northern sky every few nights. We discover 60 high-quality microlensing events in the 3 yr of ZTF-I using the bulk lightcurves in the ZTF Public Data Release 5.19 of our events are found outside of the Galactic plane (∣b∣ ≥ 10°), nearly doubling the number of previously discovered events in the stellar halo from surveys pointed toward the Magellanic Clouds and the Andromeda galaxy. We also record 1558 ongoing candidate events as potential microlensing that can continue to be observed by ZTF-II for identification. The scalable and computationally efficient methods developed in this work can be applied to future synoptic surveys, such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope, as they attempt to find microlensing events in even larger and deeper data sets.

We present Atacama Large Millimeter/submillimeter Array observations of the ∼10,000 au environment surrounding 21 protostars in the Orion A molecular cloud tracing outflows. Our sample is composed of Class 0 to flat-spectrum protostars, spanning the full ∼1 Myr lifetime. We derive the angular distribution of outflow momentum and energy profiles and obtain the first two-dimensional instantaneous mass, momentum, and energy ejection rate maps using our new approach: the pixel flux-tracing technique. Our results indicate that by the end of the protostellar phase, outflows will remove ∼2–4 M_⊙ from the surrounding ∼1 M_⊙ low-mass core. These high values indicate that outflows remove a significant amount of gas from their parent cores and continuous core accretion from larger scales is needed to replenish core material for star formation. This poses serious challenges to the concept of cores as well-defined mass reservoirs, and hence to the simplified core-to-star conversion prescriptions. Furthermore, we show that cavity opening angles, and momentum and energy distributions all increase with protostar evolutionary stage. This is clear evidence that even garden-variety protostellar outflows: (a) effectively inject energy and momentum into their environments on 10,000 au scales, and (b) significantly disrupt their natal cores, ejecting a large fraction of the mass that would have otherwise fed the nascent star. Our results support the conclusion that protostellar outflows have a direct impact on how stars get their mass, and that the natal sites of individual low-mass star formation are far more dynamic than commonly accepted theoretical paradigms.

A detailed modeling of simultaneous UV-photochemical and thermochemical processes in exoplanet atmosphere-like conditions is essential for the analysis and interpretation of a vast amount of current and future spectral data from exoplanets. However, a detailed reaction kinetic model that incorporates both UV photochemistry and thermal chemistry is challenging due to the massive size of the chemical system as well as the lack of understanding of photochemistry compared to thermal-only chemistry. Here, we utilize an automatic chemical reaction mechanism generator to build a high-fidelity thermochemical reaction kinetic model later then incorporated with UV photochemistry enhanced by metastable triplet-state carbon monoxide (a³Π). Our model results show that two different photochemical reactions driven by Lyα photons (i.e., H₂ + CO(a³Π) → H + HCO and CO(X¹Σ⁺) + CO(a³Π) → C(³P) + CO₂) can enhance thermal chemistry resulting in significant increases in the formation of CH₄, H₂O, and CO₂ in H₂-dominated systems with trace amounts of CO, which qualitatively matches the observations from previous experimental studies. Our model also suggests that at temperatures above 2000 K, thermal chemistry becomes the dominant process. Finally, the chemistry simulated up to 2500 K does not produce any larger species such as C₃ species, benzene, or larger (i.e., PAHs). This might indicate that the photochemistry of C₂ species such as C₂H₂ might play a key role in the formation of organic aerosols observed in a previous experimental study.

We present 18 yr of OGLE photometry together with spectra obtained over 12 yr revealing that the early Oe star AzV 493 shows strong photometric (ΔI < 1.2 mag) and spectroscopic variability with a dominant, 14.6 yr pattern and ∼40 day oscillations. We estimate the stellar parameters T_eff = 42,000 K, , M/M_⊙ = 50 ± 9, and v sin i = 370 ± 40 km s⁻¹. Direct spectroscopic evidence shows episodes of both gas ejection and infall. There is no X-ray detection, and it is likely a runaway star. The star AzV 493 may have an unseen companion on a highly eccentric (e > 0.93) orbit. We propose that close interaction at periastron excites ejection of the decretion disk, whose variable emission-line spectrum suggests separate inner and outer components, with an optically thick outer component obscuring both the stellar photosphere and the emission-line spectrum of the inner disk at early phases in the photometric cycle. It is plausible that AzV 493’s mass and rotation have been enhanced by binary interaction followed by the core-collapse supernova explosion of the companion, which now could be either a black hole or a neutron star. This system in the Small Magellanic Cloud can potentially shed light on OBe decretion disk formation and evolution, massive binary evolution, and compact binary progenitors.

The JWST early release data show unexpected high stellar mass densities of massive galaxies at 7 < z < 11. A high star formation efficiency is probably needed to explain this. However, such a high star formation efficiency would greatly increase the number of ionizing photons, which would be in serious conflict with current cosmic microwave background (CMB) and other measurements of cosmic reionization history. To solve this problem, we explore fuzzy dark matter (FDM), which is composed of ultra-light scalar particles, e.g., ultra-light axions, and calculate its halo mass function and stellar mass density for different axion masses. We find that a FDM model with m_a ≃ 5 × 10⁻²³ eV and a possible uncertainty range ∼3 × 10⁻²³–10⁻²² eV can effectively suppress the formation of small halos and galaxies, so that with higher star formation efficiency both the JWST data at z ∼ 8 and the reionization history measurements from optical depth of CMB scattering and ionization fraction can be simultaneously matched. We also find that the JWST data at z ∼ 10 are still too high to fit in this scenario. We note that the estimated mean redshift of the sample may have large uncertainty, that it can be as low as z ∼ 9 depending on adopted spectral energy distribution templates and photometric-redshift code. In addition, warm dark matter with ∼keV mass can also be an alternative choice, since it should have similar effects on halo formation as FDM.

Primordial B-mode detection is one of the main goals of current and future cosmic microwave background (CMB) experiments. However, the weak B-mode signal is overshadowed by several Galactic polarized emissions, such as thermal dust emission and synchrotron radiation. Subtracting foreground components from CMB observations is one of the key challenges in searching for the primordial B-mode signal. Here, we construct a deep convolutional neural network (CNN) model, called CMBFSCNN (Cosmic Microwave Background Foreground Subtraction with CNN), which can cleanly remove various foreground components from simulated CMB observational maps at the sensitivity of the CMB-S4 experiment. Noisy CMB Q (or U) maps are recovered with a mean absolute difference of 0.018 ± 0.023 μK (or 0.021 ± 0.028 μK). To remove the residual instrumental noise from the foreground-cleaned map, inspired by the needlet internal linear combination method, we divide the whole data set into two “half-split maps,” which share the same sky signal, but have uncorrelated noise, and perform a cross-correlation technique to reduce the instrumental noise effects at the power spectrum level. We find that the CMB EE and BB power spectra can be precisely recovered with significantly reduced noise effects. Finally, we apply this pipeline to current Planck observations. As expected, various foregrounds are cleanly removed from the Planck observational maps, with the recovered EE and BB power spectra being in good agreement with the official Planck results.

The conditions under which galactic nuclear regions become active are largely unknown, although it has been hypothesized that secular processes related to galaxy morphology could play a significant role. We investigate this question using optical i-band images of 3096 SDSS quasars and galaxies at 0.3 < z < 0.6 from the Hyper Suprime-Cam Subaru Strategic Program, which possesses a unique combination of area, depth, and resolution, allowing the use of residual images, after removal of the quasar and smooth galaxy model, to investigate internal structural features. We employ a variational auto-encoder, which is a generative model that acts as a form of dimensionality reduction. We analyze the lower-dimensional latent space in search of features that correlate with nuclear activity. We find that the latent space does separate images based on the presence of nuclear activity, which appears to be associated with more pronounced components (i.e., arcs, rings, and bars) as compared to a matched control sample of inactive galaxies. These results suggest the importance of secular processes and possibly mergers (by their remnant features) in activating or sustaining black hole growth. Our study highlights the breadth of information available in ground-based imaging taken under optimal seeing conditions and having an accurate characterization of the point-spread function (PSF), thus demonstrating future science to come from the Rubin Observatory.

We characterize the optical counterparts to the compact X-ray source population within the nearby spiral galaxy M81 using multiband Hubble Space Telescope (HST) imaging data. By comparing the optical luminosities and colors measured for candidate donor stars and host clusters to stellar and cluster evolutionary models, respectively, we estimate the likely masses and upper age limits of the field and cluster X-ray binaries. We identify 15 low-mass X-ray binaries (i.e., donor star mass ≲ 3 M_⊙) within ancient globular clusters, as well as 42 candidate high-mass X-ray binaries (i.e., donor star mass ≳ 8 M_⊙). To estimate the likelihood of misclassifications, we inject 4000 artificial sources into the HST mosaic image and conclude that our classifications of globular clusters and high-mass X-ray binaries are reliable at the >90% level. We find that globular clusters that host X-ray binaries are on average more massive and more compact than globular clusters that do not. However, there is no apparent correlation between the X-ray brightness of the clusters and their masses or densities, nor are X-ray binary hosts more X-ray luminous than the general field population of low-mass X-ray binaries. This work represents one of the first in-depth analyses of the population of X-ray binaries within globular clusters in a spiral galaxy.

During the inspiralling of a white dwarf (WD) into an intermediate-mass black hole (∼10²⁻⁵M_⊙), both gravitational waves (GWs) and electromagnetic (EM) radiation are emitted. Once the eccentric orbit’s pericenter radius approaches the tidal radius, the WD would be tidally stripped upon each pericenter passage. The accretion of this stripped mass would produce EM radiation. It is suspected that the recently discovered new types of transients, namely the quasiperiodic eruptions and the fast ultraluminous x-ray bursts, might originate from such systems. Modeling these flares requires a prediction of the amount of stripped mass from the WD and the details of the mass supply to the accretion disk. We run hydrodynamical simulations to study the orbital parameter dependence of the stripped mass. We find that our results match the analytical estimate that the stripped mass is proportional to z^5/2, where z is the excess depth by which the WD overfills its instantaneous Roche lobe at the pericenter. The corresponding fallback rate of the stripped mass is calculated, which may be useful in interpreting the individual flaring light curve in candidate EM sources. We further calculate the long-term mass-loss evolution of a WD during its inspiral and the detectability of the GW and EM signals. The EM signal from the mass-loss stage can be easily detected: the limiting distance is ∼320(M_h/10⁴M_⊙) Mpc for the Einstein Probe. The GW signal, for space-borne detectors such as Laser Interferometer Space Antenna or TianQin, can be detected only within the Local Supercluster (∼33 Mpc).

The Makani galaxy hosts the poster child of a galactic wind on scales of the circumgalactic medium. It consists of a two-episode wind in which the slow, outer wind originated 400 Myr ago (Episode I; R_I = 20 − 50 kpc) and the fast, inner wind is 7 Myr old (Episode II; R_II = 0 − 20 kpc). While this wind contains ionized, neutral, and molecular gas, the physical state and mass of the most extended phase—the warm, ionized gas—are unknown. Here we present Keck optical spectra of the Makani outflow. These allow us to detect hydrogen lines out to r = 30–40 kpc and thus constrain the mass, momentum, and energy in the wind. Many collisionally excited lines are detected throughout the wind, and their line ratios are consistent with 200–400 km s⁻¹ shocks that power the ionized gas, with v_shock = σ_wind. Combining shock models, density-sensitive line ratios, and mass and velocity measurements, we estimate that the ionized mass and outflow rate in the Episode II wind could be as high as those of the molecular gas: and yr⁻¹. The outer wind has slowed, so that yr⁻¹, but it contains more ionized gas, M_⊙. The momentum and energy in the recent Episode II wind imply a momentum-driven flow (p “boost” ∼7) driven by the hot ejecta and radiation pressure from the Eddington-limited, compact starburst. Much of the energy and momentum in the older Episode I wind may reside in a hotter phase, or lie further into the circumgalactic medium.

We present new radial velocity measurements from the Magellan and the Anglo-Australian Telescopes for 175 previously known and 121 newly confirmed globular clusters (GCs) around NGC 5128, the nearest accessible massive early-type galaxy at D = 3.8 Mpc. Remarkably, 28 of these newly confirmed GCs are at projected radii (≳54 kpc), extending to ∼130 kpc, in the outer halo where few GCs had been confirmed in previous work. We identify several subsets of GCs that spatially trace halo substructures that are visible in red giant branch star maps of the galaxy. In some cases, these subsets of GCs are kinematically cold, and may be directly associated with and originate from these specific stellar substructures. From a combined kinematic sample of 645 GCs, we see evidence for coherent rotation at all radii, with a higher rotation amplitude for the metal-rich GC subpopulation. Using the tracer mass estimator, we measure a total enclosed mass of 2.5 ± 0.3 × 10¹² M_⊙ within ∼120 kpc, an estimate that will be sharpened with forthcoming dynamical modeling. The combined power of stellar mapping and GC kinematics makes NGC 5128 an ongoing keystone for understanding galaxy assembly at mass scales inaccessible in the Local Group.

We compare the core-collapse evolution of a pair of 15.8 M_☉ stars with significantly different internal structures, a consequence of the bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79 M_☉ progenitors have core masses (masses interior to an entropy of 4 k_B baryon⁻¹) of 1.47 and 1.78 M_☉ and compactness parameters ξ_1.75 of 0.302 and 0.604, respectively. The core-collapse simulations are carried out in 2D to nearly 3 s postbounce and show substantial differences in the times of shock revival and explosion energies. The 15.78 M_☉ model begins exploding promptly at 120 ms postbounce when a strong density decrement at the Si–Si/O shell interface, not present in the 15.79 M_☉ progenitor, encounters the stalled shock. The 15.79 M_☉ model takes 100 ms longer to explode but ultimately produces a more powerful explosion. Both the larger mass accretion rate and the more massive core of the 15.79 M_☉ model during the first 0.8 s postbounce time result in larger ν_e/ luminosities and RMS energies along with a flatter and higher-density heating region. The more-energetic explosion of the 15.79 M_☉ model resulted in the ejection of twice as much ⁵⁶Ni. Most of the ejecta in both models are moderately proton rich, though counterintuitively the highest electron fraction (Y_e = 0.61) ejecta in either model are in the less-energetic 15.78 M_☉ model, while the lowest electron fraction (Y_e = 0.45) ejecta in either model are in the 15.79 M_☉ model.

Many of the M-dwarf stars, though they are tiny and dim, are observed to possess strong surface magnetic fields and exhibit remarkably intense flaring. Such magnetism may severely impact habitability on the exoplanets now discovered nearby. The origin of the magnetism must rest with dynamo action achieved by turbulent convection coupled to rotation within the M-dwarfs. To further explore the nature and diversity of the magnetism that can result, we turn here to an extensive set of 45 global MHD simulations to explore dynamos operating within deep convective envelopes of rapidly rotating M2 (0.4 M_⊙) stars. We observe a wide range of cycle periods present in the convection zones, whose durations we find to scale with the Rossby number as Ro^{−1.66±0.07} in concurrence with scalings identified in simulations of more massive stars. We find a unifying relationship between the ratio of magnetic to convective kinetic energy (ME/CKE) and the degree to which the differential rotation is quenched by magnetic fields. We show that the presence of a tachocline in these model stars enhances their axisymmetric magnetic field components, leading to a surface dipole fraction on average 78% greater than an equivalent star with only a CZ, potentially shedding light on the nature of the tachocline divide through resultant effects on the spin-down rate.

We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4 m telescope at the Kitt Peak National Observatory. Over the next 5 yr DESI MWS will observe approximately seven million stars at Galactic latitudes ∣b∣ > 20°, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also include several high-completeness samples of rare stellar types, including white dwarfs, low-mass stars within 100 pc of the Sun, and horizontal branch stars. We summarize the potential of DESI to advance understanding of the Galactic structure and stellar evolution. We introduce the final definitions of the main MWS target classes and estimate the number of stars in each class that will be observed. We describe our pipelines for deriving radial velocities, atmospheric parameters, and chemical abundances. We use ≃500,000 spectra of unique stellar targets from the DESI Survey Validation program (SV) to demonstrate that our pipelines can measure radial velocities to ≃1 km s⁻¹ and [Fe/H] accurate to ≃0.2 dex for typical stars in our main sample. We find the stellar parameter distributions from ≈100 deg² of SV observations with ≳90% completeness on our main sample are in good agreement with expectations from mock catalogs and previous surveys.

A number of neutron stars have been observed within the remnants of the core-collapse supernova explosions that created them. In contrast, black holes are not yet clearly associated with supernova remnants (SNRs). Indeed, some observations suggest that black holes are “born in the dark,” i.e., without a supernova explosion. Herein, we present a multiwavelength analysis of the X-ray transient Swift J1728.9−3613, based on observations made with Chandra, ESO-VISTA, MeerKAT, NICER, NuSTAR, Swift, and XMM-Newton. Three independent diagnostics indicate that the system likely harbors a black hole primary. Infrared imaging signals a massive companion star that is broadly consistent with an A or B spectral type. Most importantly, the X-ray binary lies within the central region of the cataloged SNR G351.9−0.9. Our deep MeerKAT image at 1.28 GHz signals that the remnant is in the Sedov phase; this fact and the nondetection of the soft X-ray emission expected from such a remnant argue that it lies at a distance that could coincide with the black hole. Utilizing a formal measurement of the distance to Swift J1728.9−3613 (d = 8.4 ± 0.8 kpc), a lower limit on the distance to G351.9−0.9 (d ≥ 7.5 kpc), and the number and distribution of black holes and SNRs within the Milky Way, extensive simulations suggest that the probability of a chance superposition is <1.7% (99.7% credible interval). The discovery of a black hole within an SNR would support numerical simulations that produce black holes and remnants, and thus provide clear observational evidence of distinct black hole formation channels. We discuss the robustness of our analysis and some challenges to this interpretation.

The origin and distribution of stellar-mass black hole spins are a rare window into the progenitor stars and supernova events that create them. Swift J1728.9-3613 is an X-ray binary, likely associated with the supernova remnant (SNR) G351.9-0.9. An NuSTAR X-ray spectrum of this source during its 2019 outburst reveals reflection from an accretion disk extending to the innermost stable circular orbit. Modeling of the relativistic Doppler shifts and gravitational redshifts imprinted on the spectrum measures a dimensionless spin parameter of a = 0.86 ± 0.02 (1σ confidence), a small inclination angle of the inner accretion disk θ < 10°, and a subsolar iron abundance in the disk A_Fe < 0.84. This high spin value rules out a neutron star primary at the 5σ level of confidence. If the black hole is located in a still visible SNR, it must be young. Therefore, we place a lower limit on the natal black hole spin of a > 0.82, concluding that the black hole must have formed with a high spin. This demonstrates that black hole formation channels that leave an SNR, and those that do not (e.g., Cyg X-1), can both lead to high natal spin with no requirement for subsequent accretion within the binary system. Emerging disparities between the population of high-spin black holes in X-ray binaries and the low-spin black holes that merge in gravitational wave events may therefore be explained in terms of different stellar conditions prior to collapse, rather than different environmental factors after formation.

Star formation is ubiquitously associated with the ejection of accretion-powered outflows that carve bipolar cavities through the infalling envelope. This feedback is expected to be important for regulating the efficiency of star formation from a natal prestellar core. These low-extinction outflow cavities greatly affect the appearance of a protostar by allowing the escape of shorter-wavelength photons. Doppler-shifted CO line emission from outflows is also often the most prominent manifestation of deeply embedded early-stage star formation. Here, we present 3D magnetohydrodynamic simulations of a disk wind outflow from a protostar forming from an initially 60 M_⊙ core embedded in a high-pressure environment typical of massive star-forming regions. We simulate the growth of the protostar from m_* = 1 M_⊙ to 26 M_⊙ over a period of ∼100,000 yr. The outflow quickly excavates a cavity with a half opening angle of ∼10° through the core. This angle remains relatively constant until the star reaches 4 M_⊙. It then grows steadily in time, reaching a value of ∼50° by the end of the simulation. We estimate a lower limit to the star formation efficiency (SFE) of 0.43. However, accounting for continued accretion from a massive disk and residual infall envelope, we estimate that the final SFE may be as high as ∼0.7. We examine observable properties of the outflow, especially the evolution of the cavity's opening angle, total mass, and momentum flux, and the velocity distributions of the outflowing gas, and compare with the massive protostars G35.20-0.74N and G339.88-1.26 observed by the Atacama Large Millimeter/submillimeter Array (ALMA), yielding constraints on their intrinsic properties.