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

Blazars exhibit multiwavelength variability, a phenomenon whose underlying mechanisms remain elusive. This study investigates the origin of such variability through leptonic blazar emission simulations, focusing on stochastic fluctuations in environmental parameters. By analyzing the spectral indices of the power spectral densities of the variability, we assess their relationship with the underlying fluctuations. Our findings reveal that the variability spectral indices remain almost independent of the variations responsible for their emergence. This suggests a complex interplay of factors contributing to the observed multiwavelength variability in blazars.

We present a focused X-ray and multiwavelength study of the ultraluminous weak-line quasar (WLQ) SDSS J1521+5202, one of the few X-ray weak WLQs that is amenable to basic X-ray spectral and variability investigations. J1521+5202 shows striking X-ray variability during 2006–2023, by up to a factor of ≈32 in 0.5–2 keV flux, and our new 2023 Chandra observation caught it in its brightest X-ray flux state to date. Concurrent infrared/optical observations show only mild variability. The 2023 Chandra spectrum can be acceptably described by a power law with intrinsic X-ray absorption, and it reveals a nominal intrinsic level of X-ray emission relative to its optical/ultraviolet emission. In contrast, an earlier Chandra spectrum from 2013 shows apparent spectral complexity that is not well fit by a variety of models, including ionized absorption or standard Compton-reflection models. Overall, the observations are consistent with the thick-disk plus outflow model previously advanced for WLQs, where a nominal level of underlying X-ray emission plus variable absorption leads to the remarkable observed X-ray variability. In the case of J1521+5202, it appears likely that the outflow, and not the thick disk itself, lies along our line of sight and causes the X-ray absorption.

The ability to accurately discern active massive black holes (BHs) in nearby dwarf galaxies is paramount to understanding the origins and processes of “seed” BHs in the early Universe. We present Chandra X-ray Observatory observations of a sample of three local dwarf galaxies (M_* ≤ 3 × 10⁹M_⊙, z ≤ 0.15) previously identified as candidates for hosting active galactic nuclei (AGN). The galaxies were selected from the NASA-Sloan Atlas with spatially coincident X-ray detections in the eROSITA Final Equatorial Depth Survey. Our new Chandra data reveal three X-ray point sources in two of the target galaxies with luminosities between log(L_{2−10 keV} [erg s⁻¹]) = 39.1 and 40.4. Our results support the presence of an AGN in these two galaxies and an ultraluminous X-ray source (ULX) in one of them. For the AGNs, we estimate BH masses of M_BH ∼ 10⁵⁻⁶M_⊙ and Eddington ratios on the order of ∼10⁻³.

Subhalo dynamics in galaxy cluster host halos govern the observed distribution and properties of cluster member galaxies. We use the IllustrisTNG simulation to investigate the accretion and orbits of subhalos found in cluster-size halos. We find that the median change in the major axis direction of cluster-size host halos is approximately 80° between a ∼ 0.1 and the present day. We identify coherent regions in the angular distribution of subhalo accretion, and ∼68% of accreted subhalos enter their host halo through ∼38% of the surface area at the virial radius. The majority of galaxy clusters in the sample have ∼2 such coherent regions. We further measure angular orbits of subhalos with respect to the host major axis and use a clustering algorithm to identify distinct orbit modes with varying oscillation timescales. The orbit modes correlate with subhalo accretion conditions. Subhalos in orbit modes with shorter oscillations tend to have lower peak masses and accretion directions somewhat more aligned with the major axis. One orbit mode, exhibiting the least oscillatory behavior, largely consists of subhalos that accrete near the plane perpendicular to the host halo major axis. Our findings are consistent with expectations from inflow from major filament structures and internal dynamical friction: most subhalos accrete through coherent regions, and more massive subhalos experience fewer orbits after accretion. Our work offers a unique quantification of subhalo dynamics that can be connected to how the intracluster medium strips and quenches cluster galaxies.

The central region of M87 has been observed by the Atacama Large Millimeter/submillimeter Array (ALMA) several times to resolve its molecular cloud distribution. We present the ALMA CO(1–0), CO(2–1), and CO(3–2) observations of the M87 circumnuclear region. The CO(1–0) shows emission as well as molecular absorption features toward its nucleus, indicating the existence of cold molecular clouds/gas. ALMA imaging of the CO(1–0) observations show the distribution of molecular clouds toward the nucleus of M87 with a radius of ∼105 pc. The clouds seen in the CO(1–0) are moving with very high velocity and have an apparent motion in the range of −1000 to 3600 km s⁻¹ relative to the central supermassive black hole, which has a systemic velocity of ∼1266 km s⁻¹. The molecular gas mass was estimated to be about 1.08 ± 0.64 × 10⁷M_⊙ in the central 210 pc of M87. We calculated the optical depths using the absorption lines of the CO(1–0) and CO(2–1). The excitation temperature was estimated to be ∼7.133 ± 0.03 K using the optical depth ratio of the CO(2-1) to CO(1-0). We also found the line-of-sight column density to be ∼2.25 × 10¹⁵ cm⁻², which corresponds to a hydrogen column density of ∼0.28 × 10²⁰ cm⁻².

Meridional flow is crucial in generating the solar poloidal magnetic field by facilitating poleward transport of the field from decayed bipolar magnetic regions (BMRs). As the meridional circulation changes with the stellar rotation rate, the properties of stellar magnetic cycles are expected to be influenced by this flow. In this study, we explore the role of meridional flow in generating magnetic fields in the Sun and Sun-like stars using the STABLE (surface flux transport and Babcock–Leighton) dynamo model. We find that a moderate meridional flow increases the polar field by efficiently driving the trailing polarity flux toward the pole, while a strong flow tends to transport both polarities of BMRs poleward, potentially reducing the polar field. Our findings are in perfect agreement with what one can expect from the surface flux transport model. Similarly, the toroidal field initially increases with moderate flow speeds and then decreases beyond a certain value. This trend is due to the competitive effects of shearing and diffusion. Furthermore, our study highlights the impact of meridional flow on the strength and duration of stellar cycles. By including the meridional flow from a mean-field hydrodynamics model in STABLE, we show that the magnetic field strength initially increases with the stellar rotation rate and then declines in rapidly rotating stars, offering an explanation of the observed variation of stellar magnetic field with rotation rate.

In this article we confront a class of f(Q) gravity models with observational data of galaxy–galaxy lensing. Specifically, we consider f(Q) gravity models containing a small quadratic correction when compared with general relativity (GR), and quantify this correction by a model parameter, α. To derive the observational constraints, we start by extracting spherically symmetric solutions, which correspond to deviations from the Schwarzschild solution that depends on the model parameters in a twofold way, i.e., a renormalized mass and a new term proportional to r⁻². Then, we calculate the effective lensing potential, the deflection angle, the shear component, and the effective excess surface density profile. After that, we employ a group catalog and a shape catalog from the Sloan Digital Sky Survey Data Release 7 for the lens and source samples, respectively. Moreover, we handle the off-center radius as a free parameter and constrain it using a Markov Chain Monte Carlo method. Concerning the deviation parameter from GR we derive at the 1σ confidence level, and then compare the fitting efficiency with the standard Λ cold dark matter paradigm by applying the Akaike information criterion and Bayesian information criterion. Our results indicate that the f(Q) corrections alongside the off-center effect yield a scenario that is slightly favored.

The two main competing theories proposed to explain the formation of massive (>10 M_⊙) stars—competitive accretion and monolithic core collapse—make different observable predictions for the environment of the massive stars during, and immediately after, their formation. Proponents of competitive accretion have long predicted that the most massive stars should have a different spatial distribution to lower-mass stars, through the stars being either mass segregated or being in areas of higher relative densities or sitting deeper in gravitational potential wells. We test these predictions by analyzing a suite of smoothed-particle hydrodynamics simulations where star clusters form massive stars via competitive accretion with and without feedback. We find that the most massive stars have higher relative densities, and sit in deeper potential wells, only in simulations in which feedback is not present. When feedback is included, only half of the simulations have the massive stars residing in deeper potential wells, and there are no other distinguishing signals in their spatial distributions. Intriguingly, in our simple models for monolithic core collapse, the massive stars may also end up in deeper potential wells because if massive cores fragment then the stars that form are also massive, and dominate their local environs. We find no robust diagnostic test in the spatial distributions of massive stars that can distinguish their formation mechanisms, and so other predictions for distinguishing between competitive accretion and monolithic collapse are required.

Constraining the formation processes of small solar system bodies is crucial for gaining insights into planetesimal formation. Their bulk densities, determined by their compressive strengths, offer valuable information about their formation history. In this paper, we utilize a formulation of the compressive strength of dust aggregates obtained from dust N-body simulations to establish the relation between the bulk density and diameter. We find that this relation can be effectively approximated by a polytrope with an index of 0.5, coupled with a formulation of the compressive strength of dust aggregates. The lowest-density trans-Neptunian objects (TNOs) and main-belt asteroids (MBAs) are well reproduced by dust aggregates composed of 0.1 μm sized grains. However, most TNOs, MBAs, comets, and near-Earth asteroids (NEAs) exhibit higher densities, suggesting the influence of compaction mechanisms such as collision, dust grain disruption, sintering, or melting, leading to further growth. We speculate that there are two potential formation paths for small solar system bodies. One involves the direct coagulation of primordial dust grains, resulting in the formation of first-generation planetesimals, including the lowest-density TNOs, MBAs, and the parent bodies of comets and NEAs. In this case, comets and NEAs are fragments or rubble piles of first-generation planetesimals, and the objects themselves or the rubble are composed of 0.1 μm sized grains. The other path involves the further potential fragmentation of first-generation planetesimals into the compact dust aggregates observed in protoplanetary disks, resulting in the formation of second-generation planetesimals composed of compact dust aggregates, which may contribute to explaining another formation process of comets and NEAs.

Recent radio surveys have revealed pulsars with dispersion and scattering delays induced by ionized gas that are larger than the rest of the observed pulsar population, in some cases with electron column densities (or dispersion measures, DMs) larger than the maximum predictions of Galactic electron density models. By cross-matching the observed pulsar population against H ii region catalogs, we show that the majority of pulsars with DM > 600 pc cm⁻³ and scattering delays τ(1 GHz) > 10 ms lie behind H ii regions, and that H ii region intersections may be relevant to as much as a third of the observed pulsar population. The fraction of the full pulsar population with sightlines intersecting H ii regions is likely larger. Accounting for H ii regions resolves apparent discrepancies where Galactic electron density models place high-DM pulsars beyond the Galactic disk. By comparing emission measures inferred from recombination line observations to pulsar DMs, we show that H ii regions can contribute tens to hundreds of parsecs per cubic centimeter in electron column density along a pulsar line of sight. We find that nearly all pulsars with significant excess (and deficit) scattering from the mean τ–DM relation are spatially coincident with known discrete ionized gas structures, including H ii regions. Accounting for H ii regions is critical to the interpretation of radio dispersion and scattering measurements as electron density tracers, both in the Milky Way and in other galaxies.

We perform 3D3V hybrid-Vlasov simulations of turbulence with quasi-isotropic, compressible injection near ion scales to mimic the Earth’s magnetosheath plasma, and investigate the novel electron-only reconnection, recently observed by NASA’s Magnetospheric Multiscale mission, and its impact on ion heating. Retaining electron inertia in the generalized Ohm's law enables collisionless magnetic reconnection. Spectral analysis shows a shift from kinetic Alfvén waves to inertial kinetic Alfvén and inertial whistler waves near electron scales. To distinguish the roles of inertial scale and gyroradius (d_i and ρ_i), three ion beta (β_i = 0.25, 1, 4) values are studied. Ion-electron decoupling increases with β_i, as ions become less mobile when the injection scale is closer to ρ_i than d_i, highlighting the role of ρ_i in achieving an electron magnetohydrodynamic regime at sub-ion scales. This regime promotes electron-only reconnection in turbulence with small-scale injection at β_i ≳ 1. We observe significant ion heating even at large β_i, with Q_i/ ≈ 69%, 91%, and 96% at β_i = 0.25, 1, and 4, respectively. While ion heating is anisotropic at β_i ≤ 1 (T_i,⊥ > T_i,∥), it is marginally anisotropic at β_i > 1 (T_i,⊥ ≳ T_i,∥). Our results show ion turbulent heating in collisionless plasmas is sensitive to the separation between injection scales (λ_inj) and ρ_i, β_i, and finite-k_∥ effects, necessitating further investigation for accurate modeling. These findings have implications for other collisionless astrophysical environments, like high-β plasmas in intracluster medium, where processes such as microinstabilities or shocks may inject energy near ion-kinetic scales.

We present a detailed model atmosphere analysis of massive white dwarfs with M > 0.9 M_⊙ and T_eff ≥ 11,000 K in the Montreal White Dwarf Database 100 pc sample and the Pan-STARRS footprint. We obtained follow-up optical spectroscopy of 109 objects with no previous spectral classification in the literature. Our spectroscopic follow-up is now complete for all 204 objects in the sample. We find 118 normal DA white dwarfs, including 45 massive DAs near the ZZ Ceti instability strip. There are no normal massive DBs: the six DBs in the sample are strongly magnetic and/or rapidly rotating. There are 20 massive DQ white dwarfs in our sample, and all are found in the crystallization sequence. In addition, 66 targets are magnetic (32% of the sample). We use magnetic white dwarf atmosphere models to constrain the field strength and geometry using offset dipole models. We also use magnetism, kinematics, and rotation measurements to constrain the fraction of merger remnant candidates among this population. The merger fraction of this sample increases from 25% for 0.9–1 M_⊙ white dwarfs to 49% for 1.2–1.3 M_⊙. However, this fraction is as high as % for 1.1–1.2 M_⊙ white dwarfs. Previous works have demonstrated that 5%–9% of high-mass white dwarfs stop cooling for ∼8 Gyr due to the ²²Ne distillation process, which leads to an overdensity of Q-branch stars in the solar neighborhood. We demonstrate that the overabundance of the merger remnant candidates in our sample is likely due to the same process.

Type Ia supernovae (SNe Ia) discovered at redshift z ≲ 2.5 are presumed to be produced from Population I/II stars. In this work, we investigate the production of SNe Ia from Population III binaries in the cosmological framework. We derive the SN Ia rate as a function of redshift in a theoretical context for the production of first-generation stars and evaluate the likelihood of their detection by the James Webb Space Telescope (JWST). Assuming the initial stellar mass function (IMF) favors low-mass stars, as from recent numerical simulations, we found that Population III stars may give rise to a considerable number of SNe Ia at high redshift, and Population III stars may even be the dominant SN Ia producer at z ≳ 6. In an optimistic scenario, we expect ∼1(2) SNe Ia from Population III stars at z ≈ 4(5) for a survey of area during a 3 yr period with JWST. The same survey may record more than ∼400 SNe Ia at lower redshift (z ≲ 2.5) but with only about one of them from Population III progenitors. There will be ∼six Population III SNe Ia in the same field of view at redshifts of 5–10. Observational constraints on SN Ia rates at the redshift range of 5–10 can place crucial constraints on the IMF of Population III stars.

Although supermassive black holes (SMBHs) reside in the heart of virtually every massive galaxy, it remains debated whether dwarf galaxies commonly host SMBHs. Because low-mass galaxies may retain memory of the assembly history of their black holes (BHs), probing the BH occupation fraction of local dwarf galaxies might offer insights into the growth and seeding mechanisms of the first BHs. In this work, we exploit the Western half of the eROSITA all-sky survey (covering 20,000 deg²) and compile a catalog of accreting SMBHs in local (D < 200 Mpc) dwarf galaxies. Cleaning our sample from X-ray background sources, X-ray binaries, and ultraluminous X-ray sources, we identify 74 active galactic nucleus (AGN)–dwarf galaxy pairs. Using this large and uniform sample, we derive the luminosity function of the dwarf galaxy AGN, fitting it with a power-law function and obtaining . Measuring the offset between the dwarf galaxies' centroids and the X-ray sources, we find that ≈50% of the AGN are likely off-nuclear, in agreement with theoretical predictions. We compare the BH-to-stellar mass relation of our sample with the local and high-redshift relations, finding that our sources better adhere to the former, which suggests that local AGN across different mass scales undergo similar growth histories. Finally, we compare our sources with semianalytical models: while our sample’s shallowness prevents distinguishing between different seeding models, we find that the data favor models that keep SMBHs in dwarf galaxies active at a moderate rate, motivating model improvement by comparison to AGN in the dwarf galaxy regime.

We present a numerical analysis investigating the reliability of Type Ia supernova (SN Ia) delay-time distributions recovered from individual host galaxy star formation histories. We utilize star formation histories of mock samples of galaxies generated from the IllustrisTNG simulation at two redshifts to recover delay-time distributions. The delay-time distributions are constructed through piecewise constants as opposed to typically employed parametric forms such as power laws or Gaussian or skew/lognormal functions. The SN Ia delay-time distributions are recovered through a Markov Chain Monte Carlo exploration of the likelihood space by comparing the expected SN Ia rate within each mock galaxy to the observed rate. We show that a reduced representative sample of nonhost galaxies is sufficient to reliably recover delay-time distributions while simultaneously reducing the computational load. We also highlight a potential systematic between recovered delay-time distributions and the mass-weighted ages of the underlying host galaxy stellar population.

We present infrared photometry and Hubble Space Telescope imaging and spectroscopy of the Sab galaxy NGC 4826. Schwarzschild dynamical modeling is used to measure its central black hole mass M. Photometric decomposition is used to enable a comparison of M to published scaling relations between black hole masses and properties of host bulges. This decomposition implies that NGC 4826 contains classical and pseudobulges of approximately equal mass. The classical bulge has best-fit Sérsic index n = 3.27. The pseudobulge is made up of three parts, an inner lens (n = 0.18 at r ≲ 4″), an outer lens (n = 0.17 at r ≲ 45″), and a n = 0.58 Sérsic component required to match the surface brightness between the lens components. The total V-band luminosity of the galaxy is M_VT = −21.07, the ratio of classical bulge to total light is B/T ≃ 0.12, and the ratio of pseudobulge to total light is PB/T ≃ 0.13. The outer disk is exponential (n = 1.07) and makes up D/T = 0.75 of the light of the galaxy. Our best-fit Schwarzschild model has a black hole mass with 1σ uncertainties of and a stellar population with a K-band mass-to-light ratio of at the assumed distance of 7.27 Mpc. Our modeling is marginally consistent with M = 0 at the 3σ limit. These best-fit parameters were calculated assuming the black hole is located where the velocity dispersion is largest; this is offset from the maximum surface brightness, probably because of dust absorption. The black hole mass—one of the smallest measured by modeling stellar dynamics—satisfies the well known correlations of M with the K-band luminosity, stellar mass, and velocity dispersion of the classical bulge only. In contrast, the black hole is undermassive with respect to the correlation of M with total (classical plus pseudo) bulge luminosity. Thus the composite (classical bulge plus pseudobulge) galaxy NGC 4826 is consistent with previous results on black hole scaling relations and helps to strengthen these results at low black hole masses.

While cosmic rays (E ≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales of L > 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.

We present detailed radio observations of the tidal disruption event (TDE) ASASSN-19bt/AT 2019ahk, obtained with the Australia Telescope Compact Array, the Atacama Large Millimeter/submillimeter Array, and the MeerKAT radio telescopes, spanning 40–1464 days after the onset of the optical flare. We find that ASASSN-19bt displays unusual radio evolution compared to other TDEs, as the peak brightness of its radio emission increases rapidly until 457 days post-optical discovery and then plateaus. Using a generalized approach to standard equipartition techniques, we estimate the energy and corresponding physical parameters for two possible emission geometries: a nonrelativistic spherical outflow and a relativistic outflow observed from a range of viewing angles. We find that the nonrelativistic solution implies a continuous energy rise in the outflow from E ∼ 10⁴⁶ to E ∼ 10⁴⁹ erg with outflow speed β ≈ 0.05, while the off-axis relativistic jet solution instead suggests E ≈ 10⁵² erg with Lorentz factor Γ ∼ 10 at late times in the maximally off-axis case. We find that neither model provides a holistic explanation for the origin and evolution of the radio emission, emphasizing the need for more complex models. ASASSN-19bt joins the population of TDEs that display unusual radio emission at late times. Conducting long-term radio observations of these TDEs, especially during the later phases, will be crucial for understanding how these types of radio emission in TDEs are produced.

Transient current sheets (CSs) are common magnetic structures in the solar wind that can significantly disturb the planetary space environment and cause space weather phenomena. This study focuses on their properties in the Martian space environment and investigates their transformations caused by the Martian bow shock. Based on 14 months of data from the Mars Atmosphere and Volatile Evolution mission during similar solar activity, 786 CSs with directional changes greater than 90° and magnetic field depression were automatically identified, including 256 events in the upstream solar wind and 530 events in the Martian magnetosheath. After crossing the bow shock, the duration and thickness of the CSs decrease, while their current density and occurrence rate significantly increase. In the Martian magnetosheath, CSs located on the dayside exhibit stronger current density, a relatively lower occurrence rate, and a smaller than CSs located on the nightside. From the dayside to the nightside, the predominant magnetic field component of CSs changes to due to the draping process. Moreover, a clear dawn–dusk asymmetry in CS properties emerges from the foreshock to the downstream region of the Martian bow shock. Our results reveal the properties and evolution of solar wind CSs in the Martian space environment.

We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at subsolar metallicity (i.e., the PAH–metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features from 5–18 μm as metallicity decreases continuously from solar (Z_⊙) to 0.2 Z_⊙. The total PAH-to-dust luminosity ratio remains relatively constant until reaching a threshold of ∼ 2/3 Z_⊙, below which it declines smoothly but rapidly. The PZR has been attributed to preferential destruction of the smallest grains in the hard radiation environments found at low metallicity. In this scenario, a decrease in emission from the shortest-wavelength PAH features is expected. In contrast, we find a steep decline in long-wavelength power below Z_⊙, especially in the 17 μm feature, with the shorter-wavelength PAH bands carrying an increasingly large fraction of power at low metallicity. We use newly developed grain models to reproduce the observed PZR trends, including these variations in fractional PAH feature strengths. The model that best reproduces the data employs an evolving grain size distribution that shifts to smaller sizes as metallicity declines. We interpret this as a result of inhibited grain growth at low metallicity, suggesting continuous replenishment in the interstellar medium is the dominant process shaping the PAH grain population in galaxies.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary Class 0 protostellar system BHR 71 IRS1 and IRS2 as part of the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program. We describe the ¹²CO (J = 2–1), ¹³CO (J = 2–1), C¹⁸O (J = 2–1), H₂CO (J = 3_2,1–2_2,0), and SiO (J = 5–4) molecular lines along with the 1.3 mm continuum at high spatial resolution (∼008 or ∼5 au). Dust continuum emission is detected toward BHR 71 IRS1 and IRS2, with a central compact component and extended continuum emission. The compact components are smooth and show no sign of substructures such as spirals, rings, or gaps. However, there is a brightness asymmetry along the minor axis of the presumed disk in IRS1, possibly indicative of an inclined geometrically and optically thick disk-like component. Using a position–velocity diagram analysis of the C¹⁸O line, clear Keplerian motions were not detected toward either source. If Keplerian rotationally supported disks are present, they are likely deeply embedded in their envelope. However, we can set upper limits of the central protostellar mass of 0.46 M_⊙ and 0.26 M_⊙ for BHR 71 IRS1 and BHR 71 IRS2, respectively. Outflows traced by ¹²CO and SiO are detected in both sources. The outflows can be divided into two components, a wide-angle outflow and a jet. In IRS1, the jet exhibits a double helical structure, reflecting the removal of angular momentum from the system. In IRS2, the jet is very collimated and shows a chain of knots, suggesting episodic accretion events.

The Large Magellanic Cloud (LMC) is home to many H ii regions, which may lead to significant outflows. We examine the LMC’s multiphase gas (T∼10^4-5 K) in H i, S ii, Si iv, and C iv using 110 stellar sight lines from the Hubble Space Telescope’s Ultraviolet Legacy Library of Young Stars as Essential Standards program. We develop a continuum fitting algorithm based on the concept of Gaussian process regression and identify reliable LMC interstellar absorption over v_helio = 175–375 km s⁻¹. Our analyses show disk-wide ionized outflows in Si iv and C iv across the LMC with bulk velocities of ∣v_{out, bulk}∣ ∼ 20–60 km s⁻¹, which indicates that most of the outflowing mass is gravitationally bound. The outflows’ column densities correlate with the LMC’s star formation rate surface densities (Σ_SFR), and the outflows with higher Σ_SFR tend to be more ionized. Considering outflows from both sides of the LMC as traced by C iv, we conservatively estimate a total outflow rate of and a mass-loading factor of η ≳ 0.15. We compare the LMC’s outflows with those detected in starburst galaxies and simulation predictions, and find a universal scaling relation of over a wide range of star-forming conditions (Σ_SFR ∼ 10^−4.5–10²M_⊙ yr⁻¹ kpc⁻²). Lastly, we find that the outflows are corotating with the LMC’s young stellar disk and the velocity field does not seem to be significantly impacted by external forces; we thus speculate on the existence of a bow shock leading the LMC, which may have shielded the outflows from ram pressure as the LMC orbits the Milky Way.

Several large JWST blank field observing programs have not yet discovered the first galaxies expected to form at 15 ≤ z ≤ 20. This has motivated the search for more effective survey strategies that will be able to effectively probe this redshift range. Here, we explore the use of gravitationally lensed cluster fields, which have historically been the most effective discovery tool with the Hubble Space Telescope. In this paper, we analyze the effectiveness of the most massive galaxy clusters that provide the highest median magnification factor within a single JWST NIRCam module in uncovering this population. The results of exploiting these lensing clusters to break the z > 15 barrier are compared against the results from large-area, blank-field surveys such as JADES and CEERS in order to determine the most effective survey strategy for JWST. We report that the fields containing massive foreground galaxy clusters specifically chosen to occupy the largest fraction of a single NIRCam module with high magnification factors in the source plane while containing all multiple images in the image plane within a single module provide the highest probability of both probing the 15 ≤ z ≤ 20 regime as well as discovering the highest-redshift galaxy possible with JWST. We also find that using multiple massive clusters in exchange for shallower survey depths is a more time-efficient method of probing the z > 15 regime.

The origin of X-ray emission from the resolved kiloparsec-scale jets and hotspots of many active galactic nuclei remains uncertain, particularly where the X-ray emission is separate from the radio-optical synchrotron component. Possible explanations include synchrotron emission from a second electron population and external Compton or synchrotron self-Compton processes—alternatives which imply very different physical conditions within the jet. Until recently, X-ray studies of resolved jets and hotspots have been restricted to below ∼10 keV, often showing a hard spectral index indicating a spectral peak beyond this energy range. Here we present NuSTAR observations of the nearby powerful radio galaxy Pictor A, in which we clearly detect the western hotspot at approximately 4′ from the host galaxy, the most significant detection of hotspot emission above 10 keV to date. The NuSTAR spectrum is best fit by a single power law of index Γ = 2.03 ± 0.04; an exponential cutoff gives a 1σ lower limit on the cutoff energy of 40.7 keV. We confirm previous findings of variations in the soft X-ray flux detected by Chandra over the 2000 to 2015 period, at a significance of 6.5σ. This rises to >8σ in the common 3–8 keV band using the combined 22 yr span of Chandra and NuSTAR observations. The variability of the western Pictor A hotspot strongly confirms the previously argued synchrotron nature of the X-ray emission for the hotspot, while the lower bound to the spectral cutoff energy implies electron energies in the hotspot reach up to at least a few TeV.

We have discovered a triply eclipsing triple-star system, TIC 290061484, with the shortest known outer period, P_out, of only 24.5 days. This “eclipses” the previous record set by λ Tauri at 33.02 days, which held for 68 yr. The inner binary, with an orbital period of P_in = 1.8 days, produces primary and secondary eclipses and exhibits prominent eclipse timing variations with the same periodicity as the outer orbit. The tertiary star eclipses, and is eclipsed by, the inner binary with pronounced asymmetric profiles. The inclinations of both orbits evolve on observable timescales such that the third-body eclipses exhibit dramatic depth variations in TESS data. A photodynamical model provides a complete solution for all orbital and physical parameters of the triple system, showing that the three stars have masses of 6.85, 6.11, and 7.90 M_⊙, radii near those corresponding to the main sequence, and T_eff in the range of 21,000–23,700 K. Remarkably, the model shows that the triple is in fact a subsystem of a hierarchical 2+1+1 quadruple with a distant fourth star. The outermost star has a period of ∼3200 days and a mass comparable to the stars in the inner triple. In ∼20 Myr, all three components of the triple subsystem will merge, undergo a Type II supernova explosion, and leave a single remnant neutron star. At the time of writing, TIC 290061484 is the most compact triple system and one of the tighter known compact triples (i.e., P_out/P_in = 13.7).

We derive the central stellar velocity dispersion function (VDF) for quiescent galaxies in 280 massive clusters with in IllustrisTNG300. The VDF is an independent tracer of the dark matter mass distribution of subhalos in galaxy clusters. Based on the IllustrisTNG cluster catalog, we select quiescent member subhalos with a specific star formation rate <2 × 10⁻¹¹ yr⁻¹ and stellar mass . We then simulate fiber spectroscopy to measure the stellar velocity dispersion of the simulated galaxies; we compute the line-of-sight velocity dispersions of star particles within a cylindrical volume that penetrates the core of each subhalo. We construct the VDFs for quiescent subhalos within R₂₀₀. The simulated cluster VDF exceeds the simulated field VDF for , indicating the preferential formation of large velocity dispersion galaxies in dense environments. The excess is similar in simulations and in the observations. We also compare the simulated VDF for the three most massive clusters with with the observed VDF for the two most massive clusters in the local Universe, Coma and A2029. Intriguingly, the simulated VDFs are significantly lower for . This discrepancy results from (1) a smaller number of subhalos with in TNG300 compared to the observed clusters, and (2) a significant offset between the observed and simulated M_*–σ_* relations. The consistency in the overall shape of the observed and simulated VDFs offers a unique window into galaxy and structure formation in simulations.

The evolution of halos with masses around M_h ≈ 10¹¹M_⊙ and M_h ≈ 10¹²M_⊙ at redshifts z > 9 is examined using constrained N-body simulations. The average specific mass accretion rates, , exhibit minimal mass dependence and generally agree with existing literature. Individual halo accretion histories, however, vary substantially. About one-third of simulations reveal an increase in around z ≈ 13. Comparing simulated halos with observed galaxies having spectroscopic redshifts, we find that for galaxies at z ≳ 9, the ratio between observed star formation rate and is approximately 2%. This ratio remains consistent for the stellar-to-halo mass ratio (SHMR) but only for z ≳ 10. At z ≃ 9, the SHMR is notably lower by a factor of a few. At z ≳ 10, there is an agreement between specific star formation rates (sSFRs) and . However, at z ≃ 9, observed sSFRs exceed simulated values by a factor of 2. It is argued that the mildly elevated SHMR in high-z halos with M_h ≈ 10¹¹M_⊙ can be achieved by assuming the applicability of the local Kennicutt–Schmidt law and a reduced effectiveness of stellar feedback due to deeper gravitational potential of high-z halos of a fixed mass.

Recent observations by the Parker Solar Probe (PSP) suggest that protons and heavier ions are accelerated to high energies by magnetic reconnection at the heliospheric current sheet (HCS). By solving the energetic particle transport equation in large-scale MHD simulations, we study the compression acceleration of protons and heavier ions in the reconnecting HCS. We find that the acceleration of multispecies ions results in nonthermal power-law distributions with a spectral index consistent with the PSP observations. Our study shows that the high-energy cutoff of protons can reach –1 MeV depending on the particle diffusion coefficients. We also study how the high-energy cutoff of different ion species scales with the charge-to-mass ratio . When determining the diffusion coefficients from the quasi-linear theory with a Kolmogorov magnetic power spectrum, we find that α ∼ 0.4, which is somewhat smaller than α ∼ 0.7 observed by PSP.

We report evidence of galaxy assembly bias—the correlation between galaxy properties and biased secondary halo properties at fixed halo mass (M_H)—in the stellar-to-halo mass relation for red central galaxies from the Sloan Digital Sky Survey. In the M_H = 10^11.5–10^13.5h⁻¹M_⊙ range, central galaxy stellar mass (M_*) is correlated with the number density of galaxies within 10 h⁻¹ Mpc (δ₁₀), a common proxy for halo formation time. This galaxy assembly bias signal is also present when M_H, M_*, and δ₁₀ are substituted with group luminosity, galaxy luminosity, and metrics of the large-scale density field. To associate differences in δ₁₀ with variations in halo formation time, we fitted a model that accounts for (1) errors in the M_H measured by the J. L. Tinker group catalog and (2) the level of correlation between halo formation time and M_* at fixed M_H. Fitting of this model yields that (1) errors in M_H are ∼0.15 dex and (2) halo formation time and M_* are strongly correlated (Spearman’s rank correlation coefficient ∼0.85). At fixed M_H, variations of ∼0.4 dex in M_* are associated with ∼1–3 Gyr variations in halo formation time and galaxy formation time (from stellar population fitting). These results are indicative that halo properties other than M_H can impact central galaxy assembly.

Despite its biogenic and astrochemical importance, sulfur (S), the 10th most abundant element in the interstellar medium (ISM) with a total abundance of S/H ≈ 2.2 × 10⁻⁵, largely remains undetected in molecular clouds. Even in the diffuse ISM where S was previously often believed to be fully in the gas phase, in recent years, observational evidence has suggested that S may also be appreciably depleted from the gas. What might be the dominant S reservoir in the ISM remains unknown. Solid sulfides like MgS, FeS, and SiS₂ are excluded as major S reservoirs due to the nondetection of their expected infrared spectral bands in the ISM. In this work, we explore the potential role of sulfurated polycyclic aromatic hydrocarbon (PAH) molecules—PAHs with sulfur heterocycles (PASHs)—as a sink for the missing S. Utilizing density function theory, we compute the vibrational spectra of 18 representative PASH molecules. It is found that these molecules exhibit a prominent C–S stretching band at ∼10 μm and two relatively weak C–S deformation bands at 15 and 25 μm that are not mixed with the nominal PAH bands at 6.2, 7.7, 8.6, 11.3, and 12.7 μm. If several parts per million of S (relative to H) are locked up in PAHs, the 10 μm C–S band would be detectable by Spitzer and the James Webb Space Telescope (JWST). To quantitatively explore the amount of S/H depleted in PASHs, a detailed comparison of the infrared emission spectra of PASHs with the Spitzer and JWST observations is needed.

The solar-type subgiant β Hyi has long been studied as an old analog of the Sun. Although the rotation period has never been measured directly, it was estimated to be near 27 days. As a Southern Hemisphere target, it was not monitored by long-term stellar activity surveys, but archival International Ultraviolet Explorer data revealed a 12 yr activity cycle. Previous ground-based asteroseismology suggested that the star is slightly more massive and substantially larger and older than the Sun, so the similarity of both the rotation rate and the activity cycle period to solar values is perplexing. We use two months of precise time-series photometry from the Transiting Exoplanet Survey Satellite to detect solar-like oscillations in β Hyi and determine the fundamental stellar properties from asteroseismic modeling. We also obtain a direct measurement of the rotation period, which was previously estimated from an ultraviolet activity–rotation relation. We then use rotational evolution modeling to predict the rotation period expected from either standard spin-down or weakened magnetic braking (WMB). We conclude that the rotation period of β Hyi is consistent with WMB and that changes in stellar structure on the subgiant branch can reinvigorate the large-scale dynamo and briefly sustain magnetic activity cycles. Our results support the existence of a “born-again” dynamo in evolved subgiants—previously suggested to explain the cycle in 94 Aqr Aa—which can best be understood within the WMB scenario.

Coronal mass ejections (CMEs) expel multithermal, magnetized plasma from the Sun, and when directed toward Earth, can cause extensive damage to space and ground-based electronics. To better understand the triggering, acceleration, and evolution of CMEs, it is critical to study CME plasma properties close to the Sun. High-resolution ultraviolet and extreme ultraviolet (UV-EUV) spectroscopy can give the most detailed plasma diagnostics of CMEs in the low solar corona. Unfortunately, very few spectrally resolved observations of CMEs in the low solar corona exist. However, with the recent launch of the Spectral Imaging of the Coronal Environment on board Solar Orbiter and the upcoming missions, including the EUV High-Throughput Solar Telescope (EUVST) on Solar-C and the Multi-slit Solar Explorer (MUSE), we will have the opportunity to obtain unprecedented, spectrally resolved CME observations. Using the only full EUV spectral observation of a CME by the Hinode/EUV Imaging Spectrometer, we predict the spectra that SPICE, EUVST, and MUSE are expected to observe during an off-limb CME eruption to investigate the diagnostic capabilities of each instrument. Finally, we provide a list of density-sensitive and temperature-sensitive ratios for CME plasma diagnostics along with the expected spectral atlas for each instrument to facilitate observing sequence planning.

Abundance measurements for the volatile element phosphorus are important for measuring the metallicity in interstellar and circumgalactic gas, where their accuracies are limited by uncertainties in the oscillator strengths. We report updated oscillator strength values for two resonant transitions of the dominant ion P II, the transitions at 961.041 and 963.801 Å, which have historically shown large uncertainties. Using a combination of observational measurements and highly accurate quasi-relativistic Hartree–Fock theoretical calculations, we present an updated oscillator strength of f = 0.147 ± 0.021 for the poorly constrained P II resonant transition at 961.401 Å, which arises from the ground electronic state 3s²3p²³P₀ to the excited level 3s²3p3d³D. This result utilizes archival optical spectra obtained with the Very Large Telescope for the quasar PKS 0528–250, which has a damped Lyα absorber at z = 2.811. We calculate a theoretical f-value = 0.153 for 961.401 Å consistent with our empirically derived value, and calculate a theoretical f-value = 1.79 for 963.801 Å. We also present theoretical oscillator strengths for the P II resonant transitions at 972.779, 1124.945, 1152.818, 1301.874, and 1532.533 Å, as well as for multiple P II fine-structure and excited-level transitions. The updated f-value for the P II 961 Å transition will be useful in future studies of P abundances, especially in sight lines where the 1152 Å line is saturated.

We employ observations of the 205 μm [N ii] fine structure (FS) line and radio recombination line (RRL) emission to derive the electron density in 10 well-known H ii regions. The combination of these two spectral lines (the RRL–FS line method) provides a sensitive probe of electron density in regions with n(e) ≥ 30 cm⁻³ without requiring knowledge of the size of the ionized region. By using H54α data from the Green Bank Telescope and 205 μm data from the SOFIA Airborne Observatory, we have almost identical 18″ beamwidths, removing a significant source of error for observations of H ii regions due to nonuniform density across the sources observed. The electron densities vary widely among the sources observed, from 2600 to 36,000 cm⁻³, with two low-density outliers at 94 and 520 cm⁻³. On average, these densities are a factor of 4 greater than the highest-resolution single-antenna data and a factor of almost 13 greater than the 182″ angular resolution single-antenna data having more sources in common. The total 1σ fractional uncertainties in n(e) are in the range 0.15–0.29. In the RRL–FS line method, the observationally determined quantity is proportional to ∫n²(z)dz / ∫n(z)dz. For a Gaussian density distribution much more extended than its 1/e radius, this is equal to , where n₀ is the peak electron density. The high values of electron density found are plausibly the result of the RRL–FS line technique sampling the peak of a centrally condensed density distribution.

Ultrahigh-energy cosmic rays are often characterized indirectly by analyzing the properties of secondary cosmic ray particles produced in the collisions with air nuclei. The particle number N_μ of muon and the depth of shower maximum X_max after air shower cascade are mostly studied to infer the energy and mass of the incident cosmic rays. Research has shown that there is a significant excess in the observed number of muons arriving at the ground from extensive air showers compared to the simulations using the existing cosmic ray hadronic interaction model. To explain this muon excess phenomenon, a new theoretical model, the gluon condensation model, is introduced in this paper and simulated by using the AIRES engine. We assume that the gluon condensation (GC) effect appears mainly in the first collision of the cascade, leading to a significant increase in the strangeness production, consequently, the production rate of kaons is increased, and n_K/n_π is greater than the value of the usual hadronic interaction process. In the calculation, the model assumes that only pions and kaons are produced in the GC state. The increase of strange particle yield would mean that the energy transferred from the hadronic cascade to the electromagnetic cascade through π⁰ → 2γ decay is reduced. This would in turn increase the number of muons at the ground level due to meson decays. Our model provides a new theoretical possibility to explain the muon excess puzzle.

Active galactic nuclei (AGN) feedback and its impact on their host galaxies are critical to our understanding of galaxy evolution. Here, we present a combined analysis of new high resolution ultraviolet (UV) data from the Ultraviolet Imaging Telescope (UVIT) on AstroSat and archival optical spectroscopic data from the Very Large Telescope/MUSE, for the Seyfert galaxy, NGC 1365. Concentrating on the central 5 kpc region, the UVIT images in the far- and near-UV show bright star-forming knots in the circumnuclear ring as well as a faint central source. After correcting for extinction, we found the star formation rate (SFR) surface density of the circumnuclear 2 kpc ring to be similar to other starbursts, despite the presence of an AGN outflow, as seen in [O iii] 5007 Å. On the other hand, we found fainter UV and thus lower SFR in the direction southeast of the AGN relative to northwest in agreement with observations at other wavelengths from JWST and Atacama Large Millimeter/submillimeter Array. The AGN outflow velocity is found to be lesser than the escape velocity, suggesting that the outflowing gas will rain back into the galaxy. The deep UV data have also revealed diffuse UV emission in the direction of the AGN outflow. By combining [O iii] and UV data, we found the diffuse emission to be of AGN origin.

Collisionless low-Mach-number shocks are abundant in astrophysical and space plasma environments, exhibiting complex wave activity and wave–particle interactions. In this paper, we present 2D Particle-in-Cell (PIC) simulations of quasi-perpendicular nonrelativistic (v_sh ≈ (5500–22000) km s⁻¹) low-Mach-number shocks, with a specific focus on studying electrostatic waves in the shock ramp and precursor regions. In these shocks, an ion-scale oblique whistler wave creates a configuration with two hot counterstreaming electron beams, which drive unstable electron acoustic waves (EAWs) that can turn into electrostatic solitary waves (ESWs) at the late stage of their evolution. By conducting simulations with periodic boundaries, we show that the EAW properties agree with linear dispersion analysis. The characteristics of ESWs in shock simulations, including their wavelength and amplitude, depend on the shock velocity. When extrapolated to shocks with realistic velocities (v_sh ≈ 300 km s⁻¹), the ESW wavelength is reduced to one-tenth of the electron skin depth and the ESW amplitude is anticipated to surpass that of the quasi-static electric field by more than a factor of 100. These theoretical predictions may explain a discrepancy, between PIC and satellite measurements, in the relative amplitude of high- and low-frequency electric field fluctuations.

Blazars can be detected from very large distances due to their high luminosity. However, the detection of γ-ray emission of blazars beyond z = 3 has only been confirmed for a small number of sources. Such observations probe the growth of supermassive black holes close to the peak of star formation in the history of galaxy evolution. As a result from a continuous monitoring of a sample of 80 z > 3 blazars with the Fermi Large Area Telescope (Fermi-LAT), we present the first detection of a γ-ray flare from the z = 4.31 blazar TXS 1508+572. This source showed high γ-ray activity from 2022 February to August, reaching a peak luminosity comparable to the most luminous flares ever detected with Fermi-LAT. We conducted a multiwavelength observing campaign involving XMM-Newton, the Neil Gehrels Swift Observatory, the Effelsberg 100 m radio telescope, and the Very Long Baseline Array. In addition, we make use of the monitoring programs by the Zwicky Transient Facility and the Near-Earth Object Wide-field Infrared Survey Explorer at optical and infrared wavelengths, respectively. We find that the source is particularly variable in the infrared band on daily timescales. The spectral energy distribution collected during our campaign is well described by a one-zone leptonic model, with the γ-ray flare originating from an increase of external Compton emission as a result of a fresh injection of accelerated electrons.

We report nucleosynthetic results for both ⁴⁴Ti and nickel isotopes for 18 three-dimensional (3D) core-collapse supernova (CCSN) simulations extended to ∼20 s after bounce. We find that many of our long-term models are able to achieve ⁴⁴Ti/⁵⁶Ni ratios similar to that observed in Cassiopeia A, and modern supernova models can synthesize up to 2 × 10⁻⁴M_⊙ of ⁴⁴Ti. Neutrino-driven winds and the fact that there can be simultaneous accretion and explosion in 3D models of CCSNe play central roles in its production. We conclude that the ⁴⁴Ti underproduction problem in previous CCSN models is no longer an issue. In addition, we discuss the production of both ⁵⁷Ni and stable nickel/iron ratios and compare our results to observations of SN 1987A and the Crab.

In this study, galaxy samples have been generated using mock observation techniques based on the results of TNG100-1 simulations to investigate three forms of intrinsic alignment: satellite-central alignment between the orientation of the brightest group galaxies (BGG) and the spatial distribution of their satellites, radial alignment between the satellites’ orientation and the direction toward their BGG, as well as direct alignment between the orientation of BGG and that of its satellites. Overall, the predictions of galaxy alignment generally align with observations, although minor discrepancies have been identified. For satellite-central alignment, the alignment strength and color-dependence trends are well replicated by the mock observations. Regarding radial alignment, the signals are weak but discernible, with no apparent color dependence. As for direct alignment, no signal is detected, nor is there any color dependence. We also investigate the alignment dependencies on halo or the BGG properties, and proximity effect. For satellite-central alignment, the predicted alignment signal shows a positive correlation with halo and BGG mass, consistent with observations and previous predictions. Similar correlations have also been observed with the BGG age and metallicity, which merit future observational analysis for confirmation. Proximity effects have been observed for all three types of alignment, with satellites closer to the BGG exhibiting stronger alignment signals. The influence of galaxy definition and shape determination on alignment studies is also analyzed. This study underscores the importance of employing mock observation techniques for a fair comparison between predictions and observations.

We conduct a systematic search for high-redshift galaxy overdensities at 4.9 < z_spec < 8.9 in both the Great Observatories Origins Deep Survey (GOODS)-N and GOODS-S fields using James Webb Space Telescope/Near-Infrared Camera (JWST/NIRCam) imaging from the JWST Advanced Deep Extragalactic Survey and JWST Extragalactic Medium-band Survey in addition to JWST/NIRCam wide field slitless spectroscopy from the First Reionization Epoch Spectroscopic Complete Survey. High-redshift galaxy candidates are identified using Hubble Space Telescope + JWST photometry spanning λ = 0.4–5.0 μm. We confirmed the redshifts for roughly a third of these galaxies using JWST spectroscopy over λ = 3.9–5.0 μm through identification of either Hα or around the best-fit photometric redshift. The rest-ultraviolet magnitudes and continuum slopes of these galaxies were inferred from the photometry: the brightest and reddest objects appear in more dense environments and thus are surrounded by more galaxy neighbors than their fainter and bluer counterparts, suggesting accelerated galaxy evolution within overdense environments. We find 17 significant (δ_gal ≥ 3.04, N_gal ≥ 4) galaxy overdensities across both fields (seven in GOODS-N and 10 in GOODS-S), including the two highest redshift spectroscopically confirmed galaxy overdensities to date at and (representing densities around ∼6 and ∼12 times that of a random volume). We estimate the total halo mass of these large-scale structures to be using an empirical stellar mass-to-halo mass relation, which are likely underestimates as a result of incompleteness. These protocluster candidates are expected to evolve into massive galaxy clusters with by z = 0.

We present a sample of 1165 extreme emission-line galaxies (EELGs) at 4 < z < 9 selected using James Webb Space Telescope (JWST) NIRCam photometry in the Cosmic Evolution Early Release Science (CEERS) program. We use a simple method to photometrically identify EELGs with Hβ + [O iii] (combined) or Hα emission of observed-frame equivalent width (EW) > 5000 Å. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: All spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the spectral-energy-distribution-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between Hβ + [O iii] and Hα in the neighboring broadband filters. We present examples of EELGs that could be incorrectly classified as ultrahigh redshift (z > 12) as a result of extreme Hβ + [O iii] emission blended across the reddest photometric filters. We compare the EELGs to the broader (subextreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high-EW tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. The photometrically inferred emission-line ratios are consistent with interstellar medium conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies.

Precise spectroscopic classification of planet hosts is an important tool of exoplanet research at both the population and individual system level. In the era of large-scale surveys, data-driven methods offer an efficient approach to spectroscopic classification that leverages the fact that a subset of stars in any given survey has stellar properties that are known with high fidelity. Here, we use The Cannon, a data-driven framework for modeling stellar spectra, to train a generative model of spectra from the Gaia Data Release 3 Radial Velocity Spectrometer (RVS). Our model derives stellar labels with precisions of 72 K in T_eff, 0.09 dex in logg, 0.06 dex in [Fe/H], 0.05 dex in [α/Fe], and 1.9 km s⁻¹ in v_broad for main-sequence stars observed by Gaia DR3 by transferring GALAH labels, and is publicly available at https://github.com/isabelangelo/gaiaspec. We validate our model performance on planet hosts with available Gaia RVS spectra at SNR>50 by showing that our model is able to recover stellar parameters at ≥20% improved accuracy over the existing Gaia stellar parameter catalogs, measured by the agreement with high-fidelity labels from the Spectroscopic Observations of Cool Stars survey. We also provide metrics to test for stellar activity, binarity, and reliability of our model outputs and provide instructions for interpreting these metrics. Finally, we publish updated stellar labels and metrics that flag suspected binaries and active stars for Kepler Input Catalog objects with published Gaia RVS spectra.

We present detailed multiband photometric and spectroscopic observations and analysis of a rare core-collapse supernova, SN 2021wvw, that includes photometric evolution up to 250 days and spectroscopic coverage up to 100 days postexplosion. A unique event that does not fit well within the general trends observed for Type IIP supernovae, SN 2021wvw shows an intermediate luminosity with a short plateau phase of just about 75 days, followed by a very sharp (∼10 days) transition to the tail phase. Even in the velocity space, it lies at a lower velocity compared to a larger Type II sample. The observed peak absolute magnitude is −16.1 mag in r-band, and the nickel mass is well constrained to 0.020 ± 0.006 M_⊙. Detailed hydrodynamical modeling using MESA+STELLA suggests a radially compact, low-metallicity, high-mass red supergiant progenitor (M_ZAMS = 18 M_⊙), which exploded with ∼0.2 × 10⁵¹ erg s⁻¹ leaving an ejecta mass of M_ej ≈ 5 M_⊙. Significant late-time fallback during the shock propagation phase is also seen in progenitor+explosion models consistent with the light-curve properties. As the faintest short-plateau supernova characterized to date, this event adds to the growing diversity of transitional events between the canonical ∼100 days plateau Type IIP and stripped-envelope events.

At the extreme densities reached in the core of neutron stars, it is possible that deconfined quark matter is produced. The formation of this new phase of strongly interacting matter is likely to occur via a first-order phase transition for the typical temperatures reached in astrophysical processes. The first seeds of quark matter would then form through a process of nucleation within the metastable hadronic phase. Here, we address the role of the thermal fluctuations in the hadronic composition on the nucleation of two-flavor quark matter. At finite temperature, the thermodynamic quantities in a system fluctuate around average values. Nucleation being a local process, it is possible that it occurs in a subsystem whose composition makes the nucleation easier. We will consider the total probability of the nucleation as the product between the probability that a subsystem has a certain hadronic composition different from the average in the bulk, and the nucleation probability in that subsystem. We will show how those fluctuations of the hadronic composition can increase the efficiency of nucleation already for temperatures ∼(0.1−1) keV. However, for temperatures ≲(1−10) MeV, the needed overpressure exceeds the maximum pressure reached in compact stars. Finally, for even larger temperatures the process of nucleation can take place, even taking into account finite-size effects.

Pulsars have been primarily detected by their narrow pulses or periodicity in time domain data. Interferometric surveys for pulsars are challenging due to the trade-off between beam sensitivity and beam size and the corresponding tradeoff between survey sensitivity (depth), sky coverage, and computational efforts. The detection sensitivity of time domain searches for pulsars is affected by dispersion smearing, scattering, and rapid orbital motion of pulsars in binaries. We have developed a new technique to select pulsar candidates in interferometric radio images by identifying scintillating sources and measuring their scintillation bandwidth and timescale. Identifying likely candidates allows sensitive, focused time-domain searches, saving computational effort. Pulsar scintillation is independent of its timing properties and hence offers a different selection of pulsars compared to time-domain searches. Candidates identified from this method could allow us to find hard-to-detect pulsars, such as submillisecond pulsars and pulsars in very compact, highly accelerated binary orbits. We use upgraded Giant Metrewave Radio Telescope (uGMRT) observations in the fields of PSR B1508+55, PSR J0437−4715, and PSR B0031−07 as test cases for our technique. We demonstrate that the technique correctly differentiates between the pulsar and other nonscintillating point sources and show that the extracted dynamic spectrum of the pulsar is equivalent to that extracted from the uGMRT phased array beam. We show the results from our analysis of known pulsar fields and discuss challenges in dealing with interference and instrumental effects.

Understanding plasma dynamics and nonthermal particle acceleration in 3D magnetic reconnection has been a long-standing challenge. In this paper, we explore these problems by performing large-scale fully kinetic simulations of multi-X-line plasmoid reconnection with various parameters in both the weak- and strong-guide-field regimes. In each regime, we have identified its unique 3D dynamics that lead to field-line chaos and efficient acceleration, and we have achieved nonthermal acceleration of both electrons and protons into power-law spectra. The spectral indices agree well with a simple Fermi acceleration theory that includes guide-field dependence. In the low-guide-field regime, the flux rope kink instability governs the 3D dynamics for efficient acceleration. The weak dependence of the spectra on the ion-to-electron mass ratio and β (≪1) implies that the particles are sufficiently magnetized for Fermi acceleration in our simulations. While both electrons and protons are injected at reconnection exhausts, protons are primarily injected by perpendicular electric fields through Fermi reflections and electrons are injected by a combination of perpendicular and parallel electric fields. The magnetic power spectra agree with in situ magnetotail observations, and the spectral index may reflect a reconnection-driven size distribution of plasmoids instead of the Goldreich–Sridhar vortex cascade. As the guide field becomes stronger, the oblique flux ropes of large sizes capture the main 3D dynamics for efficient acceleration. Intriguingly, the oblique flux ropes can also experience flux rope kink instability, to drive extra 3D dynamics. This work has broad implications for 3D reconnection dynamics and particle acceleration in heliophysics and astrophysics.

One of the more surprising findings after the first year of JWST observations is the large number of spatially extended galaxies (ultrared flattened objects, or UFOs) among the optically faint galaxy (OFG) population otherwise thought to be compact. Leveraging the depth and survey area of the JWST Advanced Deep Extragalactic Survey, we extend observations of the OFG population to an additional 112 objects, 56 of which are well-resolved in F444W with effective sizes, R_e > 025, more than tripling previous UFO counts. These galaxies have redshifts around 2 < z < 4, high stellar masses (), and star formation rates around ∼100–1000 M_⊙ yr⁻¹. Surprisingly, UFOs are red across their entire extents, which spatially resolved analysis of their stellar populations shows is due to large values of dust attenuation (typically A_V > 2 mag even at large radii). Morphologically, the majority of our UFO sample tends to have low Sérsic indices (n ∼ 1) suggesting that these large, massive, OFGs have little contribution from a bulge in F444W. Further, a majority have axis ratios between 0.2 < q < 0.4, which Bayesian modeling suggests that their intrinsic shapes are consistent with being a mixture of inclined disks and prolate objects with little to no contribution from spheroids. While kinematic constraints will be needed to determine the true intrinsic shapes of UFOs, it is clear that an unexpected population of large, disky or prolate objects contributes significantly to the population of OFGs.

ACT-CL J0034.4+0225 is a previously unrecognized merging galaxy cluster at z = 0.38588 ± 0.00068. Our primary evidence is provided by a 21 ks Chandra image that shows two surface brightness peaks separated by ∼49″ (259 kpc) surrounded by an extended cluster gas distribution. Each gas peak contains a brightest cluster galaxy, offset from the gas peak. We collect new South African Large Telescope optical spectra that, when augmented by archival data, yield redshifts for the two BGCs and 58 other cluster members. Archival Giant Metrewave Radio Telescope and MeerKAT data reveal a radio halo that encompasses the X-ray peaks. We provide and compare three X-ray-based mass estimates (5.0 × 10¹⁴M_⊙, 6.4 × 10¹⁴M_⊙, and 8.6 × 10¹⁴M_⊙). The Planck and ACT Sunyaev–Zel’dovich masses are ≈5.8 × 10¹⁴M_⊙. We constrain the merger state and properties by comparing them to an existing suite of N-body/hydrodynamical models using the measured gas peak separation (259 kpc, projected) and radial velocity difference (0–1000 km s⁻¹). This constrains the epoch of the merger to be within ∼100 Myr of first pericenter passage. A strong lensing analysis constrains the mass ratio to be in the range 1:1–1:20, while the cluster morphology prefers values near the equal-mass range.

We present multiwavelength polarization measurements of the luminous blazar Mrk 501 over a 14 month period. The 2–8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100 ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optical–infrared polarization measurements were also available in the B, V, R, I, and J bands, as were radio polarization measurements from 4.85 GHz to 225.5 GHz. Among the first five IXPE observations, we did not find significant variability in the X-ray polarization degree and angle with IXPE. However, the most recent sixth observation found an elevated polarization degree at >3σ above the average of the other five observations. The optical and radio measurements show no apparent correlations with the X-ray polarization properties. Throughout the six IXPE observations, the X-ray polarization degree remained higher than, or similar to, the R-band optical polarization degree, which remained higher than the radio value. This is consistent with the energy-stratified shock scenario proposed to explain the first two IXPE observations, in which the polarized X-ray, optical, and radio emission arises from different regions.

We present the discovery of a luminous X-ray active galactic nucleus (AGN) in the dwarf galaxy merger RGG 66. The black hole is predicted to have a mass of M_BH ∼ 10^5.4 M_⊙ and to be radiating close to its Eddington limit (L_bol/L_Edd ∼ 0.75). The AGN in RGG 66 is notable both for its presence in a late-stage dwarf–dwarf merger and for its luminosity of L_{2–10 keV} = 10^42.2 erg s⁻¹, which is among the most powerful AGNs known in nearby dwarf galaxies. The X-ray spectrum has a best-fit photon index of Γ = 2.4 and an intrinsic absorption of N_H ∼ 10²¹ cm⁻². These results come from a follow-up Chandra X-ray Observatory study of four irregular/disturbed dwarf galaxies with evidence for hosting AGNs based on optical spectroscopy. The remaining three dwarf galaxies do not have detectable X-ray sources with upper limits of L_{2–10 keV} ≲ 10⁴⁰ erg s⁻¹. Taken at face value, our results on RGG 66 suggest that mergers may trigger the most luminous of AGNs in the dwarf galaxy regime, just as they are suspected to do in more massive galaxy mergers.

We report the temporal evolution of electron pitch angle distributions behind the dipolarization front (DF) in the Earth's magnetotail with observations of the Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft. Taking advantage of multipoint observations from the THEMIS mission lined up in space, we study pitch angle distributions of energetic electrons behind the DF during two typical events. Pancake, rolling-pin, and cigar distributions are observed sequentially during the acceleration process. Based on Liouville's theorem, it is revealed that pancake distribution is dominantly formed by betatron acceleration in the early stage, and rolling-pin distribution is generated by both dominant Fermi and weak betatron acceleration in the transition stage, while cigar distribution is formed by Fermi acceleration finally. Our results provide comprehensive in situ observational evidence of the temporal evolution of electron acceleration behind the DF during propagation.

YO bands are conspicuous in the spectra of S stars. The laboratory spectrum of the A²Π–X²Σ⁺ electronic transition has been recorded with a Fourier transform spectrometer using a composite-wall hollow cathode lamp source. The A²Π–X²Σ⁺ transition with and v″ ≤ 4 has been rotationally analyzed using the modern spectral fitting program PGOPHER. Vibronic band strengths were calculated using an ab initio transition dipole moment function. A line list for the A²Π–X²Σ⁺ transition is provided and can be utilized in the modeling of S stars.

We have surveyed ³He-rich events on the Solar Orbiter mission from 2020 April to 2024 April, selecting isolated injections whose rollover ³He spectral shape is presumed to represent the initial acceleration state, unprocessed by subsequent activity such as coronal mass ejections or jets. A main goal has been to find relationships between the spectra of ³He and heavy ions C–Fe, in order to explore a common acceleration mechanism in spite of the fact that these events show ³He enrichments of up to ∼10⁴, while the heavy-ion enrichment is rarely larger than ∼10. Selecting 34 ³He injections, we find that heavy ions are always present, and arrive at the same time as the ³He signaling a common origin. Concentrating on Fe since it is a minor ion but with higher abundance than many others, we find its spectral shape and intensity is similar to ³He. In ∼two-thirds of the cases, if the ³He spectrum is shifted to lower energy by a factor 3.0 ± 1.3, it nearly coincides with the Fe spectrum, illustrating their close connection. Several plasma wave turbulence models have calculated spectra that also show the ion rollovers around 1 MeV nucleon⁻¹. The unique mass-to-charge ratio of ³He allows it to interact more efficiently with the turbulence, thereby gaining several times more energy per nucleon than the other heavy ions. In the spectral rollover region this can lead to the observed enormous enhancements of ³He. The acceleration appears to be associated with magnetic reconnection in emerging flux regions on the Sun.

Classification of gamma-ray bursts (GRBs) has been a long-standing puzzle in high-energy astrophysics. Recent observations challenge the traditional short versus long viewpoint, where long GRBs are thought to originate from the collapse of massive stars and short GRBs from compact binary mergers. Machine learning (ML) algorithms have been instrumental in addressing this problem, revealing five distinct GRB groups within the Swift Burst Alert Telescope (BAT) light-curve data, two of which are associated with kilonovae (KNe). In this work, we extend our analysis to the Fermi Gamma-ray Burst Monitor catalog and identify five clusters using unsupervised ML techniques, consistent with the Swift/BAT results. These five clusters are well separated in the fluence-duration plane, hinting at a potential link between fluence, duration, and complexities (or structures) in the light curves of GRBs. Further, we confirm two distinct classes of KN-associated GRBs. The presence of GRB 170817A in one of the two KN-associated clusters lends evidence to the hypothesis that this class of GRBs could potentially be produced by binary neutron star mergers. The second KN-associated GRB cluster could potentially originate from neutron star–black hole mergers. Future multimessenger observations of compact binaries in gravitational waves and electromagnetic waves can be paramount in understanding these clusters better.

The Large High Altitude Air Shower Observatory (LHAASO) has detected very-high-energy gamma rays from the low-ionization nuclear emission-line region galaxy NGC 4278, which has a low-luminosity active galactic nucleus (LLAGN) and symmetric, mildly relativistic S-shaped twin jets detected by radio observations. Few LLAGNs have been detected in gamma rays due to their faintness. Earlier, several radio-emitting components were detected in the jets of NGC 4278. We model their radio emission with synchrotron emission of ultra-relativistic electrons to estimate the strength of the magnetic field inside these components within a time-dependent framework after including the ages of the different components. We show that the synchrotron and synchrotron self-Compton emission by these components cannot explain the Swift X-ray data and the LHAASO gamma-ray data from NGC 4278. We suggest that a separate component in one of the jets is responsible for the high-energy emission, whose age, size, magnetic field, and the spectrum of the ultra-relativistic electrons inside it have been estimated after fitting the multiwavelength data of NGC 4278 with the sum of the spectral energy distributions from the radio components and the high-energy component. We note that the radio components of NGC 4278 are larger than the high-energy component, which has also been observed in several high-luminosity active galactic nuclei.

Shock drift acceleration (SDA) plays an important role in generating high-energy electrons at quasi-perpendicular shocks, but its efficiency in low-beta plasmas is questionable. In this article, we perform a two-dimensional particle-in-cell simulation of a low-Mach-number, low-plasma-beta quasi-perpendicular shock, and find that the electron cyclotron drift instability is unstable at the leading edge of the shock foot, which is excited by the relative drift between the shock-reflected ions and the incident electrons. The electrostatic waves triggered by the electron cyclotron drift instability can scatter and heat the incident electrons, which facilitates their escape from the shock’s loss cone. These electrons are then reflected by the shock and energized by SDA. In this way, the acceleration efficiency of SDA at low-plasma-beta quasi-perpendicular shocks is highly enhanced.

PG 1448+273 is a luminous, nearby (z = 0.0645), narrow-line Seyfert 1 galaxy, which likely accretes close to the Eddington limit. Previous X-ray observations of PG 1448+273 with XMM-Newton in 2017 and NuSTAR in 2022 revealed the presence of an ultrafast outflow, as seen through its blueshifted iron (Fe) K absorption profile, where the outflow velocity appeared to vary in the range 0.1−0.3c. In this work, new X-ray observations of PG 1448+273 are presented, in the form of four simultaneous XMM-Newton and NuSTAR observations performed in 2023 July and August. The X-ray spectra appeared at a similar flux in each observation, making it possible to analyze the mean 2023 X-ray spectrum at high signal-to-noise ratio. A broad (σ = 1 keV) and highly blueshifted (E = 9.8 ± 0.4 keV) Fe K absorption profile is revealed in the mean spectrum. The profile can be modeled by a fast, geometrically thick accretion disk wind, which reveals a maximum terminal velocity of v_∞ = −0.43 ± 0.03c, one of the fastest known winds in a nearby active galactic nucleus. As a result, the inferred mass outflow rate of the wind may reach a significant fraction of the Eddington accretion rate.

We present a sample of 335 Mira variables, extracted from DR9 of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey. These variables are characterized by the Balmer emissions (Hδ, Hγ, Hβ, and Hα) and the metal emissions (Fe iλλ 4202, 4308, 4376, and Mg iλ4571) observed in M giant spectra. We distinguished oxygen-rich stars from carbon-rich stars through the identification of carbon molecular bands present in the optical spectra. For the oxygen-rich stars we examined multiple attributes, such as the link between line strength and bolometric luminosity, and the connection between atmospheric parameters and their periods. We observed that Fe iλλ 4202, 4308 showed a significantly gradual progression, which can be postulated to trace the fluorescent emission resulting from pulsation shocks. Regarding the correlation between T_eff and the period, T_eff remains relatively constant over varying periods, with no clear trend, while both log g and [Fe/H] show a decreasing trend within a period range of 450 days, and this decreasing of the log g trend is consistent with the results in the literature. To shed more light on the variations of Balmer lines, we showcase time-series spectra for two objects, demonstrating that the Balmer lines reach their peak intensity during the brightest phase of the stellar cycle.

Due to observational challenges, the mass function of black holes (BH) at lower masses is poorly constrained in the local universe. Understanding the occupation fraction of BHs in low-mass galaxies is crucial for constraining the origins of supermassive BH seeds. Compact stellar systems (CSSs), including ultracompact dwarf galaxies (UCDs) and compact elliptical galaxies (cEs), are potential intermediate-mass BH hosts. Despite the difficulties posed by their limited spheres of influence, stellar dynamical modeling has been effective in estimating central BH masses in CSSs. Some CSSs may harbor a BH constituting up to 20% of their host stellar mass, while others might not have a central BH. In support of our ongoing efforts to determine the BH masses in select CSSs in the Virgo cluster using JWST/NIRSpec IFU observations and orbit-superposition dynamical models, we create mock kinematic data mimicking the characteristics of observed cEs/UCDs in the Virgo cluster with different BH masses. We then construct a series of dynamical models using the orbit-superposition code FORSTAND and explore the accuracy of recovering the BH mass. We find that the mass of BHs comprising 1% or more of the total host stellar mass can be accurately determined through kinematic maps featuring higher-order velocity moments. We also assess how BH mass measurement is affected by deprojection methods, regularization factors, anisotropy parameters, orbit initial conditions, the absence of higher-order velocity moments, the spatial resolution, and the signal-to-noise ratio.

This paper investigates the origin of the γ-ray emission from MGRO J1908+06 in the GeV–TeV energy band. By analyzing the data collected by the Fermi Large Area Telescope, the Very Energetic Radiation Imaging Telescope Array System, and High Altitude Water Cherenkov, with the addition of spectral data previously reported by LHAASO, a multiwavelength study of the morphological and spectral features of MGRO J1908+06 provides insight into the origin of the γ-ray emission. The mechanism behind the bright TeV emission is studied by constraining the magnetic field strength, the source age, and the distance through detailed broadband modeling. Both spectral shape and energy-dependent morphology support the scenario that inverse Compton emission of an evolved pulsar wind nebula associated with PSR J1907+0602 is responsible for the MGRO J1908+06 γ-ray emission with a best-fit true age of T = 22 ± 9 kyr and a magnetic field of B = 5.4 ± 0.8 μG, assuming the distance to the pulsar d_PSR = 3.2 kpc.

The Kelvin–Helmholtz instability (KHI), characterized by vortices forming at a perturbed velocity shear layer, is a prominent candidate mechanism for mass, momentum, and energy transport across boundaries with velocity shear in various space plasma environments. It is of particular interest at the flanks of Earth’s magnetopause, which separates the plasma of the magnetosphere from the adjacent shocked solar wind flow in the magnetosheath. In the present study, we use local hybrid-Vlasov simulations to investigate the ion velocity distribution functions (VDFs) associated with KHI in a magnetopause-like, transverse velocity shear layer setting (magnetic field perpendicular to the shear plane). We look for signatures of ion finite Larmor radius (FLR) effects, which could be utilized in spacecraft measurements to recognize when such effects are active, influencing KHI evolution and driving plasma mixing. We show that when a density/temperature asymmetry exists across the shear layer, FLR effects produce a heat flux along the vortex edges. With a magnitude (≳0.1 mW m⁻²) that is a significant fraction of the total magnetosheath energy flux, the heat flux provides a distinct signature that could be measured with a single spacecraft. During the late nonlinear stage of KHI, mixed non-Maxwellian ion VDFs are additionally found within the vortices. Our results are also valid in the presence of a small magnetic shear across the magnetopause.

Spectropolarimetric results of Fraunhofer lines between 516.3 and 532.6 nm are presented in local upper solar chromosphere and inner corona below a height of about 0.04 solar radius above the solar limb. The data were acquired on 2013 November 3 during a total solar eclipse in Gabon by the prototype Fiber Arrayed Solar Optical Telescope. It is found that the linear polarizations of the Fraunhofer lines in these layers depend strongly on specific spectral lines and positions. A Fraunhofer line at Mg ib₁518.4 nm can have a polarization amplitude up to 0.36% with respect to the continuum polarization level, while polarizations of lines like Fe i/Cr i524.7 nm are often merged in the noise level of 6.0 × 10⁻⁴. The polarizations of the Fraunhofer lines, like the emission ones and the continuum, increase with height as a whole trend, and their amplitudes can be close to those of emission ones yielded in close positions, and generally larger than those of the continuum. Rotations of the polarization directions of the Fraunhofer lines are often accompanied by variations in their polarization amplitudes and profile shapes. It is also judged from these polarimetric properties, along with other evidence, that neutral metal atoms exist in these atmospheric layers.

Using Swedish 1 m Solar Telescope Crisp Imaging Spectro-Polarimeter 6563 Å (Hα) observations and Mancha3D simulations, we analyze the formation and evolution of falling knots beneath a hedgerow prominence. By comparing the observed knot widths and kinematics to those of a parametric survey of simulations, we estimate the range of magnetic field values and characteristic wavelengths to test if the magnetic Rayleigh–Taylor instability (MRTI) can provide a physically meaningful explanation. We recover observational parameters using a novel semiautomated method and find knot velocities with a mean of −9.68 km s⁻¹ and a mean width of 614 km. Our simulations survey a range of critical wavelengths, λ_c, of 100 to 500 km, and magnetic field strengths, B₀, of 1 to 20 G, finding the closest match to observations around λ_c = 300 km, and B₀ = 2 to 6 G. As both the observational and simulated values match expected values, we conclude that the MRTI can provide a physically meaningful explanation of this observation. Additionally, we also predict that the Daniel K. Inouye Solar Telescope will be able to observationally recover secondary instabilities on the leading edge of the falling mass through applying a point-spread function to an example from the simulated results.

We present a combined Stratospheric Observatory for Infrared Astronomy (SOFIA) and two-epoch Hubble Space Telescope (HST) study of the region NGC 2071 IR containing the outflow from the protostar HOPS 361-A (NGC 2071 IRS 1). New [O i] 63 μm spectra, taken with the German REciever for Astronomy at Terahertz frequencies (GREAT) instrument aboard SOFIA, trace the outflow, spectrally resolve line profiles, and give radial velocities of the far-IR emission to yield 3D kinematics, shock properties, and feedback properties. Proper motions and extinction values from HST give a more complete quasi-3D picture that we use to model shock speeds and densities. We identify 14 distinct features that show proper motions. Combined with the GREAT data, we estimate that HOPS 361-A has a mass outflow rate of 3.16 × 10⁻⁶M_⊙ yr⁻¹ with shock speeds up to 51 km s⁻¹. The outflow force summed over all knots per area is measured to be 1.45 × 10⁻⁴M_⊙ km⁻¹ s⁻¹ yr⁻¹ and a power summed over all knotsper area of 2.37 × 10⁻²M_⊙ s⁻³ dissipated assuming constant , with large uncertainties on the feedback values. The total [O i] luminosity using all SOFIA apertures encompassing the HOPS 361 clump is 0.7 L_⊙, corresponding to a mass outflow rate of 7 × 10⁻⁵M_⊙ yr⁻¹. This is the second 3D kinematic study of an outflow in NGC 2071 IR in the near- and far-IR regimes following Rubinstein et al. (2023). The SOFIA [O i] spectra add unique focus to the entrained gas connecting the recent episode of outflows from HOPS 361-A to the outflow cavity and the surrounding NGC 2071 IR star-forming environment.

The supermassive black holes (M_BH ∼ 10⁶–10¹⁰M_⊙) that power luminous active galactic nuclei (AGNs), i.e., quasars, generally show a correlation between thermal disk emission in the ultraviolet (UV) and coronal emission in hard X-rays. In contrast, some “massive” black holes (mBHs; M_BH ∼ 10⁵–10⁶M_⊙) in low-mass galaxies present curious X-ray properties with coronal radiative output up to 100× weaker than expected. To examine this issue, we present a pilot study incorporating Very Large Array radio observations of a sample of 18 high-accretion-rate (Eddington ratios L_bol/L_Edd > 0.1), mBH-powered AGNs (M_BH ∼ 10⁶M_⊙) with Chandra X-ray coverage. Empirical correlations previously revealed in samples of radio-quiet, high-Eddington AGNs indicate that the radio–X-ray luminosity ratio, L_R/L_X, is approximately constant. Through multiwavelength analysis, we instead find that the X-ray-weaker mBHs in our sample tend toward larger values of L_R/L_X even though they remain radio-quiet per their optical–UV properties. This trend results in a tentative but highly intriguing correlation between L_R/L_X and X-ray weakness, which we argue is consistent with a scenario in which X-rays may be preferentially obscured from our line of sight by a “slim” accretion disk. We compare this observation to weak emission-line quasars (AGNs with exceptionally weak broad-line emission and a significant X-ray-weak fraction) and conclude by suggesting that our results may offer a new observational signature for finding high-accretion-rate AGNs.

We consider Roche lobe overflow (RLO) from a low-mass star on a nearly circular orbit, onto a supermassive black hole (SMBH). If mass transfer is unstable, its rate accelerates in a runaway process, resulting in highly super-Eddington mass accretion rates, accompanied by an optically thick outflow emanating from the SMBH vicinity. This produces a 1–4 week long, bright optical/ultraviolet flare, accompanied by a 1–10 year long X-ray precursor and post-cursor emitted from the accretion flow onto the SMBH. Such “Circular Tidal Disruption Events” (TDEs) represent a new class of nuclear transients, occurring at up to 1%–10% of the canonical parabolic tidal disruption event rate. Near-breakup rotation and strong tidal deformation of the star prior to disruption could lead to strong magnetic fields, making circular TDEs possible progenitors of jetted TDEs. Outflows prior to the final stellar disruption produce a circumnuclear environment (CNM) with ∼10⁻² M_⊙ at distances of ∼0.01–0.1 pc, likely leading to bright radio emission, and also similar to the CNM inferred for jetted TDEs. We discuss broader connections between circular TDEs and other recently identified classes of transients associated with galactic nuclei, such as repeating TDEs and quasi-periodic X-ray eruptions, as well as possible connections to luminous fast blue optical transients such as AT2018cow. We also discuss observational signatures of the analogous RLO of a white dwarf around an intermediate-mass BH, which may be a multimessenger source in the LISA era.

Multispacecraft observations of ³He-rich solar energetic particle (SEP) events are scarce, but much needed in order to understand and properly constrain the source and transport of these remarkably enriched ³He SEP events. In this paper, we report 15 ³He-rich SEP events that were detected by the Advanced Composition Explorer, the Solar Terrestrial Relations Observatory, and Solar Orbiter near 1 au during Solar Orbiter’s aphelion pass at the end of 2022 and early 2023. Three (five) of these events were detected simultaneously by at least two (three) spacecraft at up to ∼40° longitudinal separation, while seven events were detected by only a single spacecraft, even though an adjacent spacecraft was less than 20° apart. Using a magnetic connectivity tool, we show statistically that there is a >50% probability of detection when the spacecraft-modeled footpoints have an angular separation angle of <24° to the potential source region back at the Sun. This supports previous studies suggesting that the source of these ³He-rich SEP events is narrow in longitudinal extent. On the other hand, the magnetic connectivity due to the presence of coronal mass ejections, footpoint motion, and/or field-line meandering may also lead to difference in a detection at 1 au.

The Local Volume Complete Cluster Survey is an ongoing program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts 0.03 < z < 0.12 and X-ray luminosities in the 0.1–2.4 keV band L_X500c > 10⁴⁴ erg s⁻¹) with the Dark Energy Camera, capturing data in the u, g, r, i, z bands with a 5σ point source depth of approximately 25th–26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19° within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32° within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales.

We have used the Dark Energy Camera (DECam) on the CTIO Blanco 4 m telescope to perform a new emission-line survey of the Large Magellanic Cloud (LMC) using narrowband Hα and [S ii] filters in addition to a continuum band to create pure emission-line images. We refer to this new survey as DeMCELS, to distinguish it from the earlier Magellanic Cloud Emission-line Survey (MCELS). DeMCELS covers ∼54 deg², encompassing most of the bright optical disk of the LMC. With DECam's pixel size of only 027, our DeMCELS survey provides a seeing-limited improvement of 3–5 times over MCELS and is comparable in depth, with surface brightness limits of and in Hα and [S ii], respectively. DeMCELS provides detailed morphological information on nebulae of all scales, from the largest supershells to individual H ii regions and supernova remnants, to bubbles of emission surrounding individual stars, and even to faint structures in the diffuse ionized gas of the LMC. Many complex regions of emission show significant variations in the ratio of [S ii] to Hα—a sign of a mixture of shocks from stellar winds and/or supernovae with photoionization by embedded hot, young stars. We present the details of the observing strategy and data processing for this survey, and show selected results in comparison with previous data. A companion project for the Small Magellanic Cloud is in progress and will be reported separately. We are making these new data available to the community at large via NOIRLab’s Data Lab site.

Version 11 of the chianti database and software package is presented. Advanced ionization equilibrium models have been added for low charge states of seven elements (C, N, O, Ne, Mg, Si, and S), and represent a significant improvement especially when modeling the solar transition region. The models include the effects of higher electron density and charge transfer on ionization and recombination rates. As an illustration of the difference these models make, a synthetic spectrum is calculated for an electron pressure of 7 × 10¹⁵ cm⁻³ K and compared with an active region observation from HRTS. Increases are seen in factors of 2–5 in the predicted radiances of the strongest lines in the UV from Si iv, C iv, and N v, compared to the previous modeling using the coronal approximation. Much better agreement (within 20%) with the observations is found for the majority of the lines. The new atomic models better equip both those who are studying the transition region and those who are interpreting the emission from higher-density astrophysical and laboratory plasma. In addition to the advanced models, several ion data sets have been added or updated, and data for the radiative recombination energy loss rate have been updated.

We report on the detection of radio bursts from the Galactic bulge using the real-time transient detection and localization system, realfast. The pulses were detected commensally on the Karl G. Jansky Very Large Array during a survey of unidentified Fermi γ-ray sources. The bursts were localized to subarcsecond precision using realfast fast-sampled imaging. Follow-up observations with the Green Bank Telescope detected additional bursts from the same source. The bursts do not exhibit periodicity in a search up to periods of 480 s, assuming a duty cycle of <20%. The pulses are nearly 100% linearly polarized, show circular polarization up to 12%, and exhibit variable scattering on timescales of months. The arcsecond-level realfast localization links the source confidently with the Fermi γ-ray source and places it nearby (though not coincident with) an XMM-Newton X-ray source. Based on the source’s overall properties, we discuss various options for the nature of this object and propose that it could be a young pulsar, a magnetar, or a binary pulsar system.

Cyclopropenethione falls into the class of complex organic molecules but has not yet been observed in the interstellar medium or any circumstellar disks. However, its existence is very likely, and thus this study provides high-level ab initio predictions, which may serve as reference data for future observations or experimental work. Rovibrational configuration interaction theory based on multidimensional potential energy surfaces being obtained from explicitly correlated coupled-cluster calculations has been used to predict the fundamental vibrational modes, the microwave spectrum, and low-lying rovibrational transitions. Rotational constants as well as quartic and sextic centrifugal distortion constants were obtained from vibrational perturbation theory.

The results of simulations of magnetic reconnection accompanied by electron and proton heating and energization in a macroscale system are presented. Both species form extended power-law distributions that extend nearly three decades in energy. The primary drive mechanism for the production of these nonthermal particles is Fermi reflection within evolving and coalescing magnetic flux ropes. While the power-law indices of the two species are comparable, the protons overall gain more energy than electrons, and their power law extends to higher energy. The power laws roll into a hot thermal distribution at low energy with the transition energy occurring at lower energy for electrons compared with protons. A strong guide field diminishes the production of nonthermal particles by reducing the Fermi drive mechanism. In solar flares, proton power laws should extend down to tens of keV, far below the energies that can be directly probed via gamma-ray emission. Thus, protons should carry much more of the released magnetic energy than expected from direct observations.

The cores of active galactic nuclei are potential accelerators of 10–100 TeV cosmic rays, in turn producing high-energy neutrinos. This picture was confirmed by the compelling evidence of a TeV neutrino signal from the nearby active galaxy NGC 1068, leaving open the question of what is the site and mechanism of cosmic-ray acceleration. One candidate is the magnetized turbulence surrounding the central supermassive black hole. Recent particle-in-cell simulations of magnetized turbulence indicate that stochastic cosmic-ray acceleration is nonresonant, in contrast to the assumptions of previous studies. We show that this has important consequences on a self-consistent theory of neutrino production in the corona, leading to a more rapid cosmic-ray acceleration than previously considered. The turbulent magnetic-field fluctuations needed to explain the neutrino signal are consistent with a magnetically powered corona. We find that strong turbulence, with turbulent magnetic energy density higher than 1% of the rest-mass energy density, naturally explains the normalization of the IceCube neutrino flux, in addition to the neutrino spectral shape. Only a fraction of the protons in the corona, which can be directly inferred from the neutrino signal, are accelerated to high energies. Thus, in this framework, the neutrino signal from NGC 1068 provides a testbed for particle acceleration in magnetized turbulence.

Crush curves are of fundamental importance to numerical modeling of small and porous astrophysical bodies. The empirical literature often measures them for silica grains, and different studies have used various methods, sizes, textures, and pressure conditions. Here, we review past studies and supplement further experiments in order to develop a full and overarching understanding of the silica crush curve behavior. We suggest a new power-law function that can be used in impact simulations of analog materials similar to microgranular silica. We perform a benchmarking study to compare this new crush curve to the parametric quadratic crush curve often used in other studies, based on the study case of the DART impact onto the asteroid Dimorphos. We find that the typical quadratic crush curve parameters do not closely follow the silica crushing experiments, and as a consequence, they under (over) estimate compression close (far) from the impact site. The new crush curve presented here, applicable to pressures between a few hundred Pa and up to 1.1 GPa, might therefore be more precise. Additionally, it is not calibrated by case-specific parameters, and can be used universally for comet- or asteroid-like bodies, given an assumed composition similar to microgranular silica.

Asteroseismic modeling is a powerful way to derive stellar properties. However, the derived quantities are limited by built-in assumptions used in stellar models. This work presents a detailed characterization of stellar model uncertainties in asteroseismic red giants, focusing on the mixing-length parameter α_MLT, the initial helium fraction Y_init, the solar abundance scale, and the overshoot parameters. First, we estimate error floors due to model uncertainties to be ≈0.4% in mass, ≈0.2% in radius, and ≈17% in age, primarily due to the uncertain state of α_MLT and Y_init. The systematic uncertainties in age exceed typical statistical uncertainties, suggesting the importance of their evaluation in asteroseismic applications. Second, we demonstrate that the uncertainties from α_MLT can be entirely mitigated by direct radius measurements or partially through . Utilizing radii from Kepler eclipsing binaries, we determined the α_MLT values and calibrated the α_MLT–[M/H] relation. The correlation observed between the two variables is positive, consistent with previous studies using 1D stellar models, but in contrast with outcomes from 3D simulations. Third, we explore the implications of using asteroseismic modeling to test the scaling relation. We found that a perceived dependency of on [M/H] from individual frequency modeling can be largely removed by incorporating the calibrated α_MLT–[M/H] relation. Variations in Y_init can also affect predictions. These findings suggest that conveys information not fully captured by individual frequencies, and that it should be carefully considered as an important observable for asteroseismic modeling.

Barium stars are peculiar stars with enhanced slow neutron capture process (s-process) elements. Abundance analysis of them aids in better understanding the chemical evolution of the Milky Way. In this paper, we introduce a data-driven method named the memory-enhanced adaptive spectral network (MEASNet) to search for barium candidates in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution survey (LRS) and estimate the abundance of five s-process elements: Sr, Y, Ba, Ce, and Nd. MEASNet, trained using spectra from common stars in both LAMOST and the Galactic Archaeology with HERMES survey, showcases notable performance: for the classification task, precision = 98.22% and recall = 94.12%; in prediction, the mean absolute error for the seven elements range between 0.07 and 0.15 dex. After training, we apply the model to 4,083,003 stellar spectra from LAMOST DR10 LRS, successfully identifying 1,803,670 spectra of barium candidates ([Ba/Fe] ≥ 0.25 dex) along with their five s-process elemental abundances. The catalog enlarges the sample size, providing a wealth of data for further statistical analysis of the formation and evolution of barium stars. Meanwhile, this work highlights the potential value of MEASNet in star classification and abundance estimation, offering a strong reference for future data-driven models.

We present a comprehensive study based on multiwavelength observations from the NuSTAR, NICER, Swift, Fermi, NEOWISE, and ATCA telescopes during the 2022 outburst of the black-hole X-ray binary IGR J17091–3624. Our investigation concentrates on the heartbeat-like variability in the X-ray emission, with the aim of using it as a tool to unravel the origin of the nonthermal emission during the heartbeat state. Through X-ray timing and spectral analysis, we observe that the heartbeat-like variability correlates with changes in the disk temperature, supporting the disk radiation pressure instability scenario. Moreover, in addition to a Comptonization component, our time-averaged and phase-resolved spectroscopy reveal the presence of a power-law component that varies independently from the disk component. Combined with the radio–X-ray spectral energy distribution fitting, our results suggest that the power-law component could originate from synchrotron self-Compton radiation in the jet, which requires a strong magnetic field of about B = (0.3–3.5) × 10⁶ G. Additionally, assuming that IGR J17091-3624 and GRS 1915 + 105 share the same radio–X-ray correlation coefficient during both the hard and the heartbeat states, we obtain a distance of 13.7 ± 2.3 kpc for IGR J17091–3624.

Some electromagnetic outbursts from the nuclei of distant galaxies have been found to repeat on months-to-years timescales, and each of these sources can putatively arise from the accretion flares generated through the repeated tidal stripping of a star on a bound orbit about a supermassive black hole (SMBH), i.e., a repeating partial tidal disruption event (rpTDE). Here, we test the rpTDE model through analytical estimates and hydrodynamical simulations of the interaction between a range of stars, which differ from one another in mass and age, and an SMBH. We show that higher-mass (≳1M_⊙), evolved stars can survive many (≳10−100) encounters with an SMBH while simultaneously losing few × 0.01M_⊙, resulting in accretion flares that are approximately evenly spaced in time with nearly the same amplitude, quantitatively reproducing ASASSN-14ko. We also show that the energy imparted to the star via tides can lead to a change in its orbital period that is comparable to the observed decay in the recurrence time of ASASSN-14ko’s flares, . Contrarily, lower-mass and less-evolved stars lose progressively more mass and produce brighter accretion flares on subsequent encounters for the same pericenter distances, leading to the rapid destruction of the star and cessation of flares. Such systems cannot reproduce ASASSN-14ko-like transients, but are promising candidates for recreating events such as AT2020vdq, which displayed a second and much brighter outburst compared to the first. Our results imply that the lightcurves of repeating transients are tightly coupled with stellar type.

We present a suite of high-resolution numerical simulations to study the evolution and survival of dust in hot galactic winds. We implement a novel dust framework in the Cholla hydrodynamics code and use wind tunnel simulations of cool, dusty clouds to understand how thermal sputtering affects the dust content of galactic winds. Our simulations illustrate how various regimes of cloud evolution impact dust survival, dependent on cloud size, wind properties, and dust grain size. We find that significant amounts of dust can survive in winds in all scenarios, even without shielding from the cool phase of outflows. We present an analytic framework that explains this result, along with an analysis of the impact of cloud evolution on the total fraction of dust survival. Using these results, we estimate that 60% of 0.1 μm dust that enters a starburst-driven wind could survive to populate both the hot and cool phases of the halo, based on a simulated distribution of cloud properties. We also investigate how these conclusions depend on grain size, exploring grains from 0.1 μm to 10 Å. Under most circumstances, grains smaller than 0.01 μm cannot withstand hot-phase exposure, suggesting that the small grains observed in the circumgalactic medium (CGM) are either formed in situ due to the shattering of larger grains, or must be carried there in the cool phase of outflows. Finally, we show that the dust-to-gas ratio of clouds declines as a function of distance from the galaxy due to cloud–wind mixing and condensation. These results provide an explanation for the vast amounts of dust observed in the CGMs of galaxies and beyond.

3C 264 is one of the few FRI radio galaxies with detected TeV emission. It is a low-luminosity active galactic nucleus (LLAGN) and is generally associated with a radiatively inefficient accretion flow (RIAF). Earlier multiwavelength studies suggest that the X-ray emission originates from a jet. However, the possibility that the RIAF can significantly contribute to the X-rays cannot be ruled out. In particular, hard X-ray emission ≳10 keV has never been detected, making it challenging to distinguish between X-ray models. Here we report a NuSTAR detection up to 25 keV from 3C 264. We also present subpixel deconvolved Chandra images to resolve jet emission down to ∼02 from the center of the unresolved X-ray core. Together with a simultaneous Swift observation, we have constrained the dominant hard X-ray emission to be from its unresolved X-ray core, presumably in its quiescent state. We found evidence of a cutoff in the energy around 20 keV, indicating that at least some of the X-rays from the core can be attributed to the RIAF. The Comptonization model suggests an electron temperature of about 15 keV and an optical depth ranging between 4 and 7, following the universality of coronal properties of black hole accretion. The cutoff energy or electron temperature of 3C 264 is the lowest among those of other LLAGNs. The detected hard X-ray emission is at least an order of magnitude higher than that predicted by synchrotron self-Compton models introduced to explain γ-ray and TeV emission, suggesting that the synchrotron electrons might be accelerated to higher energies than previously thought.

Characterizing the chemistry of complex organic molecules (COMs) at the epoch of planet formation provides insights into the chemical evolution of the interstellar medium (ISM) and the origin of organic materials in our solar system. We report a detection of dimethyl ether (CH₃OCH₃) in the disk around the Herbig Ae star MWC 480 with sensitive Atacama Large Millimeter/submillimeter Array observations. This is the first detection of CH₃OCH₃ in a nontransitional Class II disk. The spatially unresolved, compact (≲25 au in radius) nature, broad line width (∼30 km s⁻¹), and high excitation temperature (∼200 K) indicate the sublimation of COMs in the warm inner disk. Despite the detection of CH₃OCH₃, methanol (CH₃OH), the most abundant COM in the ISM, has not been detected, from which we constrain the column density ratio of CH₃OCH₃/CH₃OH ≳ 7. This high ratio may indicate the reprocessing of COMs during the disk phase, as well as the effect of the physical structure in the inner disk. We also find that this ratio is higher than in COM-rich transition disks recently discovered. This may indicate that in the full disk of MWC 480, COMs have experienced substantial chemical reprocessing in the innermost region, while the COM emission in the transition disks predominantly traces the inherited ice sublimating at the dust cavity edge located at larger radii (≳20 au).

In the past, high-z active galactic nuclei (AGNs) were given a minor role as possible drivers of reionization, despite initial evidence in favor of their large space densities at low luminosities by Chandra and the Hubble Space Telescope. Recent observations from JWST are finding relatively large numbers of faint AGNs at z > 4, convincingly confirming these early results. We present a sample of z ∼ 5 AGNs, both from wide, shallow ground-based surveys and from deep, pencil-beam observations from JWST, allowing us to estimate their space densities with unprecedented accuracy. The bright end (M₁₄₅₀ < −26) of the z ∼ 5 AGN luminosity function is well constrained, with a rather steep slope. The faint end (M₁₄₅₀ ≥ −22) indicates a high space density, the scatter is significant, and the knee (M₁₄₅₀ ∼ −24) is mostly undetermined. Comparisons with state-of-the-art models find reasonable agreement with the observed AGN luminosity function at z = 5, while the predicted space density evolution at higher redshifts appears to be too fast with respect to observational constraints. Given the large variance at the faint end, we consider different options in fitting the luminosity functions and deriving the ionizing emissivity. Even in the most conservative scenario, the photoionization rate produced by z ∼ 5 AGNs is consistent with the ultraviolet background measurements. A slow evolution of the space density of faint AGNs is observed, indicating that active SMBHs are probably producing large amounts of ionizing photons at z > 6, well into the Epoch of Reionization. This is an important indication that high-z AGNs could be major contributors to the reionization of the Universe.

Metal-poor massive stars drive the evolution of low-mass galaxies, both locally and at high redshift. However, quantifying the feedback they impart to their local surroundings remains uncertain because models of stellar evolution, mass loss, and ionizing spectra are unconstrained by observations below 20% solar metallicity (Z_⊙). We present new Keck Cosmic Web Imager optical spectroscopy of three O-type stars in the nearby dwarf galaxies Leo P, Sextans A, and WLM, which have gas-phase oxygen abundances of 3%–14% Z_⊙. To characterize their fundamental stellar properties and radiation-driven winds, we fit PoWR atmosphere models to the optical spectra simultaneously with Hubble Space Telescope far-ultraviolet (FUV) spectra and multiwavelength photometry. We find that all three stars have effective temperatures consistent with their spectral types and surface gravities typical of main-sequence dwarf stars. Yet, the combination of those inferred parameters and luminosity for the two lower-Z stars is not reproduced by stellar evolution models, even those that include rotation or binary interactions. The scenario of multiple-star systems is difficult to reconcile with all available data, suggesting that these observations pose a challenge to current evolution models. We highlight the importance of validating the relationship between stellar mass, temperature, and luminosity at very low Z for accurate estimates of ionizing photon production and spectral hardness. Finally, all three stars’ FUV wind profiles reveal low mass-loss rates and terminal wind velocities in tension with expectations from widely adopted radiation-driven wind models. These results provide empirical benchmarks for future development of mass-loss and evolution models for metal-poor stellar populations.

We conduct a reverberation mapping (RM) campaign to spectroscopically monitor a sample of selected bright active galactic nuclei with large anticipated broad-line region (BLR) sizes adequate for spectroastrometric (SA) observations by the GRAVITY instrument on the Very Large Telescope Interferometer. We report the first results for five objects, IC 4329A, Mrk 335, Mrk 509, Mrk 1239, and PDS 456, among which Mrk 1239 and PDS 456 are for the first time spectroscopically monitored. We obtain multiyear monitoring data and perform multicomponent spectral decomposition to extract the broad Hβ profiles. We detect significant time lags between the Hβ and continuum variations, generally obeying the previously established BLR size–luminosity relation. Velocity-resolved Hβ time lags illustrate diverse, possibly evolving, BLR kinematics. We further measure the Hβ line widths from mean and rms spectra and the resulting virial products show good consistency among different observing seasons. Adopting a unity virial factor and the FWHM of the broad Hβ line from the mean spectra as a measure of velocity, the obtained black hole masses averaged over seasons are , , , , and for the five objects, respectively. Black hole mass estimations using other line width measures are also reported. For objects with previous RM campaigns, our mass estimates are in agreement with earlier results. In a companion paper, we will employ BLR dynamical modeling to directly infer the black hole mass and thereby determine the virial factors.

The size, shape, and composition of the heliosphere are affected by its interaction with the local interstellar medium. A means of quantifying this interaction is by measuring the relative motion of the two media. We determine the flow direction, , of the interstellar wind through the heliosphere and its trend over an 11 yr solar cycle using Advanced Composition Explorer/Solar Wind Ion Composition Spectrometer He⁺ double-coincidence pickup ion (PUI) measurements from a data set spanning 13 full orbits, with focusing cone crossings occurring between the years 1998 and 2010. We measure the flow direction by fitting a kappa function to the focusing cone signature in the count rates of He⁺ at the PUI cutoff measured as a function of ecliptic longitude. We determine for each three-orbit boxcar centered on the years 1999 through 2009 and find the trend of the resulting linear fit. We account for solar transients by removing measurements associated with CMEs. However, there still may be effects from compression regions, such as stream–stream and corotating interaction regions. We additionally make considerations to ensure the PUI torus is in view, to account for varying interplanetary magnetic field angles, to ensure that we analyze newly injected interstellar PUIs largely unaffected by transport effects, and to account for general solar activity. We measure a slope in the flow direction to be 0.00° ± 0.51°​​​​ yr⁻¹, indicating that we do not measure a varying trend over this period. We repeat the focusing cone analysis for the combined data set and determine an overall flow direction of 75.37° ± 0.43°.

We perform a thorough analysis of the projected shapes of nearby galaxies in both observations and cosmological simulations. We implement a forward-modeling approach to overcome the limitations in previous studies, which hinder accurate comparisons between observations and simulations. We measure axis ratios of z = 0 (snapshot 99) TNG50 galaxies from their synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images and compare them with those obtained from real HSC-SSP images of a matched galaxy sample. Remarkably, the comparison shows excellent agreement between the observations and the TNG50 simulation, challenging previous claims that ΛCDM models underproduced the abundance of thin galaxies. Specifically, for galaxies with stellar masses , we find ≲0.1σ tensions between the observations and the simulation, a stark contrast to the previously reported ≳10σ tensions. We reveal that low-mass galaxies in TNG50 are thicker than their observed counterparts in HSC-SSP and attribute this to the spurious dynamical heating effects that artificially puff up galaxies. We also find that, despite the overall broad agreement, TNG50 galaxies are more concentrated than the HSC-SSP ones at the low- and high-mass end of the stellar mass range of and are less concentrated at intermediate stellar masses. But we argue that the higher concentrations of the low-mass TNG50 galaxies are not likely the cause of their thicker/rounder appearances. Our study underscores the critical importance of conducting mock observations of simulations and applying consistent measurement methodologies to facilitate proper comparison with observations.

It is generally recognized that the electromagnetic multipolar emission from magnetars can be used to explain radiation from soft gamma repeaters or anomalous X-ray pulsars, but they have little impact on the spin-down of magnetars. We here present an analytical solution for the neutron star multipolar electromagnetic fields and their associated expected luminosities. We find that for newborn millisecond magnetars, the spin-down luminosity from higher multipolar components can match or even exceed that from the dipole component. Such high-intensity radiation will undoubtedly affect related astrophysical phenomena at the birth of a magnetar. We show that the spin-down luminosity from multipoles can well explain the majority of gamma-ray burst (GRB) afterglows, from the plateau starting at several hundred seconds until the normal decay phase lasting for many years. The fitted magnetar parameters for GRB afterglows are all typical values, with spins in the millisecond range and magnetic field strengths on the order of 10¹⁴–10¹⁵ G. Our results, in turn, provide support for the hypothesis that GRBs originate from the birth of magnetars with a period of a few milliseconds, thus deepening our understanding of the complex magnetic field structure and the equation of state of magnetars.

In this study, we undertake a spectral-timing analysis of the black hole X-ray binary source GRS 1915+105 using simultaneous observations carried out by AstroSat (Large Area X-ray Proportional Counter, LAXPC and Soft X-ray Telescope, SXT) and Neutron Star Interior Composition Explorer (NICER) in 2017. The source showed two flux levels (high and low), whose energy spectra can be described by the thermal comptonization of disk photons. The spectral parameters obtained by the joint fitting of SXT/LAXPC and NICER/LAXPC were consistent. The power density spectra from LAXPC and NICER revealed a broad, prominent feature at ∼2 Hz. The energy dependence of the fractional rms and time lag of this feature cannot be explained by only variations of coronal spectral parameters. Instead, a model where the coronal heating rate varies first and induces a change in the disk temperature and inner radius can explain the variation. Our results underline the importance of simultaneous observations by AstroSat and NICER and highlight the need for more sophisticated models to explain the spectral-temporal behavior of black hole systems.

We present the results of the XMM-Newton and NuSTAR observations taken as part of the ongoing, intensive multiwavelength monitoring program of the Seyfert 1 galaxy Mrk 817 by the AGN Space Telescope and Optical Reverberation Mapping 2 (AGN STORM 2) Project. The campaign revealed an unexpected and transient obscuring outflow, never before seen in this source. Of our four XMM-Newton/NuSTAR epochs, one fortuitously taken during a bright X-ray state has strong narrow absorption lines in the high-resolution grating spectra. From these absorption features, we determine that the obscurer is in fact a multiphase ionized wind with an outflow velocity of ∼5200 km s⁻¹, and for the first time find evidence for a lower ionization component with the same velocity observed in absorption features in the contemporaneous Hubble Space Telescope spectra. This indicates that the UV absorption troughs may be due to dense clumps embedded in diffuse, higher ionization gas responsible for the X-ray absorption lines of the same velocity. We observe variability in the shape of the absorption lines on timescales of hours, placing the variable component at roughly 1000 R_g if attributed to transverse motion along the line of sight. This estimate aligns with independent UV measurements of the distance to the obscurer suggesting an accretion disk wind at the inner broad line region. We estimate that it takes roughly 200 days for the outflow to travel from the disk to our line of sight, consistent with the timescale of the outflow's column density variations throughout the campaign.

In this paper we describe the survey design for the Ultradeep NIRSpec and NIRCam Observations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program, which executed its early imaging component in 2022 November. The UNCOVER survey includes ultradeep (∼29–30AB) imaging of ∼45 arcmin² on and around the well-studied A2744 galaxy cluster at z = 0.308 and will follow up ∼500 galaxies with extremely deep low-resolution spectroscopy with the NIRSpec/PRISM during the summer of 2023, with repeat visits in summer 2024. We describe the science goals, survey design, target selection, and planned data releases. We also present and characterize the depths of the first NIRCam imaging mosaic, highlighting previously unparalleled resolved and ultradeep 2–4 μm imaging of known objects in the field. The UNCOVER primary NIRCam mosaic spans 28.8 arcmin² in seven filters (F115W, F150W, F200W, F277W, F356W, F410M, and F444W) and 16.8 arcmin² in our NIRISS parallel (F115W, F150W, F200W, F356W, and F444W). To maximize early community use of the Treasury data set, we publicly release the full reduced mosaics of public JWST imaging including 45 arcmin² NIRCam and 17 arcmin² NIRISS mosaics on and around the A2744 cluster, including the Hubble Frontier Field primary and parallel footprints.

Measuring the properties of the cold neutral medium (CNM) in low-metallicity galaxies provides insights into heating and cooling mechanisms in early Universe-like environments. We report detections of two localized atomic neutral hydrogen (H i) absorption features in NGC 6822, a low-metallicity (0.2 Z_⊙) dwarf galaxy in the Local Group. These are the first unambiguous CNM detections in a low-metallicity dwarf galaxy outside the Magellanic Clouds. The Local Group L-band Survey (LGLBS) enabled these detections, due to its high spatial (15 pc for H i emission) and spectral (0.4 km s⁻¹) resolution. We introduce LGLBS and describe a custom pipeline for searching for H i absorption at high angular resolution and extracting associated H i emission. A detailed Gaussian decomposition and radiative transfer analysis of the NGC 6822 detections reveals five CNM components, with key properties: a mean spin temperature of 32 ± 6 K, a mean CNM column density of 3.1 × 10²⁰ cm⁻², and CNM mass fractions of 0.33 and 0.12 for the two sightlines. Stacking nondetections does not reveal low-level signals below our median optical depth sensitivity of 0.05. One detection intercepts a star-forming region, with the H i absorption profile encompassing the CO (2−1) emission, indicating coincident molecular gas and a depression in high-resolution H i emission. We also analyze a nearby sightline with deep, narrow H i self-absorption dips, where the background warm neutral medium is attenuated by intervening CNM. The association of CNM, CO, and Hα emissions suggests a close link between the colder, denser H i phase and star formation in NGC 6822.

We examine Ulysses magnetic field observations from 1993 to 1996 as the spacecraft made its first fast-latitude scan from the southern to the northern hemisphere. Most of the observations we use are representative of high-latitude solar minimum conditions. We examine magnetic field power spectra characteristics of interplanetary turbulence at high frequencies, where the spectrum breaks from an inertial range into the ion dissipation range. The onset and spectral index of the dissipation spectrum are consistent with low-latitude observations at 1 au. Both ranges have a ratio of power in perpendicular magnetic field components to parallel components near 3. The power spectrum ratio test developed by Bieber et al. for single-spacecraft analyses that determines the underlying anisotropy of the wave vectors yields only marginally more energy associated with field-aligned wave vectors than perpendicular wave vectors when comparing the inertial and dissipation-range spectra. The lack of significant change in the anisotropies between the inertial and dissipation ranges contrasts strongly with the turbulence found typically for 1 au near-ecliptic observations, where significant differences in both anisotropies are observed.

Fragmentation contributes to the formation and evolution of stars. Observationally, high-mass stars are known to form multiple-star systems, preferentially in cluster environments. Theoretically, Jeans instability has been suggested to determine characteristic fragmentation scales, and thermal or turbulent motion in the parental gas clump mainly contributes to the instability. To search for such a characteristic fragmentation scale, we have analyzed Atacama Large Millimeter/submillimeter Array (ALMA) 1.33 mm continuum observations toward 30 high-mass star-forming clumps taken by the Digging into the Interior of Hot Cores with ALMA survey. We have identified 573 cores using the dendrogram algorithm and measured the separation of cores by using the Minimum Spanning Tree technique. The core separation corrected by projection effects has a distribution peaked around 5800 au. In order to remove biases produced by different distances and sensitivities, we further smooth the images to a common physical scale and perform completeness tests. Our careful analysis finds a characteristic fragmentation scale of ∼7000 au, comparable to the thermal Jeans length of the clumps. We conclude that thermal Jeans fragmentation plays a dominant role in determining the clump fragmentation in high-mass star-forming regions, without the need to invoke turbulent Jeans fragmentation.

With nearly two billion stars observed and their corresponding astrometric parameters evaluated in the recent Gaia mission, the number of astrometric binary candidates has risen significantly. Due to the surplus of astrometric data, the current computational methods employed to inspect these astrometric binary candidates are both computationally expensive and cannot be executed in a reasonable time frame. In light of this, a machine learning (ML) technique to automatically classify whether a set of stars belongs to an astrometric binary pair via an artificial neural network (ANN) is proposed. Using data from Gaia Data Release 3, the ANN was trained and tested on 1.5 million highly probable true and visual binaries, considering the proper motions, parallaxes, and angular and physical separations as features. The ANN achieves high classification scores, with an accuracy of 99.3%, a precision rate of 0.988, a recall rate of 0.991, and an area under the curve of 0.999, indicating that the utilized ML technique is a highly effective method for classifying astrometric binaries. Thus, the proposed ANN is a promising alternative to the existing methods for the classification of astrometric binaries.

We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) observations of warm CO and H₂O gas in emission toward the low-mass protostar IRAS 15398-3359, observed as part of the CORINOS program. The CO is detected via the rovibrational fundamental band and hot band near 5 μm, whereas the H₂O is detected in the rovibrational bending mode at 6–8 μm. Rotational analysis indicates that the CO originates in a hot reservoir with an excitation temperature of 1598 ± 118 K, while the water is much cooler at 204 ± 7 K. Neither the CO nor the H₂O line images are significantly spatially extended, constraining the emission to within ∼40 au of the protostar. The compactness and high temperature of the CO are consistent with an origin in the embedded protostellar disk, or in a compact disk wind. In contrast, the water must arise from a cooler region and requires a larger emitting area (compared to the CO) to produce the observed fluxes. The water may arise from a more extended part of the disk, or from the inner portion of the outflow cavity. Thus, the origin of the molecular emission observed with JWST remains ambiguous. Better constraints on the overall extinction, comparison with realistic disk models, and future kinematically resolved observations may all help to pinpoint the true emitting reservoirs.

The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test models of radiatively inefficient accretion flows (RIAFs). Previous fits to the quiescent Chandra ACIS-S spectrum found that a RIAF model with an equal inflow–outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Program, which displays features suggestive of temperature and velocity structures within the plasma. A comprehensive forward model analysis accounting for the accretion flow geometry and HETG-S instrumental effects is required for a full interpretation of the quiescent Chandra HETG-S spectrum. We present a RIAF model that takes these effects into account. Our fits to the high-resolution grating spectrum indicate an inflow balanced by an outflow (s ∼ 1) alongside a temperature profile that appears shallower than what would be expected from a gravitational potential following 1/r. The data require that the abundance of iron relative to solar is Z_Fe < 0.32 Z_⊙ (90% credible interval), much lower than the 2 Z_⊙ metallicity measured in nearby late-type giants. While future missions like NewAthena will provide higher spectral resolution, source separation will continue to be a problem. Leveraging Chandra’s unparalleled spatial resolution, which is not expected to be surpassed for decades, remains essential for detailed investigations of the densely populated Galactic center in X-rays.

Hydrodynamic simulations of the stellar winds from Wolf–Rayet stars within the Galactic center can provide predictions for the X-ray spectrum of the supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. and Russell et al., and find that a “cold” (10⁴ K) gas disk forms around Sgr A* with a simulation runtime of 3500 yr. This result is consistent with previous grid-based simulations, demonstrating that a cold disk can form regardless of numerical method. We examine the plasma scenarios arising from an environment with and without this cold disk, by generating synthetic spectra for comparison to the quiescent Fe Kα Sgr A* spectrum from Chandra HETGS, taken through the Chandra X-ray Visionary Program. We find that current and future X-ray missions are unlikely to distinguish between the kinematic signatures in the plasma in these two scenarios. Nonetheless, the stellar wind plasma model presents a good fit to the dispersed Chandra spectra within 15 of Sgr A*. We compare our results to the radiatively inefficient accretion flow (RIAF) model fit to the HETGS spectrum presented in Paper I and find that the Bayesian model evidence does not strongly favor either model. With 9″ angular resolution and high spectral resolution of the X-IFU, NewAthena will offer a clearer differentiation between the RIAF plasma model and hydrodynamic simulations, but only a future X-ray mission with arcsecond resolution will significantly advance our understanding of Sgr A*’s accretion flow in X-rays.

About 25%–50% of white dwarfs (WDs) are found to be polluted by heavy elements. It has been argued that the pollution could be caused by the tidal disruption of an approaching planet around the WD, during which a large number of clumps would be produced and would finally fall onto the WD. The reason that the planet approaches the WD is usually believed to be due to gravitational perturbations from another distant planet or stellar companion. However, the dynamics of the perturbations and the detailed partial disruption process are still poorly understood. In this study, we present an in-depth investigation of these issues. A triple system composed of a WD, an inner orbit planet, and an outer orbit planet is considered. The inner planet would be partially disrupted periodically over its long-term evolution. Fragments generated in the process are affected by gravitational perturbations from the remnant planet, facilitating their fall toward the WD. The mass-loss rate of the inner planet depends on both its internal structure and also on the orbital configuration of the planetary system.

X-ray polarization provides a new way to probe accretion geometry in black hole systems. If the accretion geometry of black holes is similar regardless of mass, we should expect the same to be true of their polarization properties. We compare the polarimetric properties of all nonblazar black holes observed with the Imaging X-ray Polarimetry Explorer. We find that their polarization properties are very similar, particularly in the hard state, where the corona dominates. This tentatively supports the idea that stellar and supermassive black holes share a common coronal geometry.

Detecting planet signatures in protoplanetary disks is fundamental to understanding how and where planets form. In this work, we report dust and gas observational hints of planet formation in the disk around 2MASS J16120668-301027, as part of the Atacama Large Millimeter/submillimeter Array (ALMA) Large Program “AGE-PRO: ALMA survey of Gas Evolution in Protoplanetary disks.” The disk was imaged with the ALMA at Band 6 (1.3 mm) in dust continuum emission and four molecular lines: ¹²CO(J = 2–1), ¹³CO(J = 2–1), C¹⁸O(J = 2–1), and H₂CO(J = 3_(3,0)–2_(2,0)). Resolved observations of the dust continuum emission (angular resolution of ∼150 mas, 20 au) show a ring-like structure with a peak at 057 (75 au), a deep gap with a minimum at 024 (31 au), an inner disk, a bridge connecting the inner disk and the outer ring, along with a spiral arm structure, and a tentative detection (to 3σ) of a compact emission at the center of the disk gap, with an estimated dust mass of ∼2.7−12.9 Lunar masses. We also detected a kinematic kink (not coincident with any dust substructure) through several ¹²CO channel maps (angular resolution ∼200 mas, 30 au), located at a radius of ∼0875 (115.6 au). After modeling the ¹²CO velocity rotation around the protostar, we identified a purple tentative rotating-like structure at the kink location with a geometry similar to that of the disk. We discuss potential explanations for the dust and gas substructures observed in the disk and their potential connection to signatures of planet formation.

The first detailed photometric and spectroscopic analysis of the G-type eclipsing binary KM UMa is presented, which indicates that the system is a short-period detached eclipsing binary. The radial velocity curves were calculated using the cross-correlation function method based on Large Sky Area Multi-Object Fiber Spectroscopic Telescope, Sloan Digital Sky Survey, and our observations, which determined the mass ratio as q = 0.45 (±0.04). Based on the light curves from the Transiting Exoplanet Survey Satellite, other survey data, and our multiband observations, the positive and negative O’Connell effects have been detected evolving gradually and alternately over the last 20 yr, which can be explained by the presence of spots on the primary component. A superflare event was detected in the SuperWASP data on 2007 February 28, further indicating that KM UMa is a very active system. We calculated its energy to be 5 × 10³⁴ erg by assuming it occurred on the primary star. Utilizing hundreds of medium-resolution spectra and one low-resolution spectrum, the equivalent width variations of the H_α line were calculated, indicating the presence of a 5.21 (±0.67) yr magnetic activity cycle. The orbital period variations were analyzed using the O–C method, detecting a long-term decrease superimposed with a periodic variation. The amplitude of the cyclic variation is 0.01124 (±0.00004) day, with a period of 33.66 (±0.0012) yr, which exceeds the 5.21 yr activity cycle, suggesting that this is more likely attributable to the light travel time effect of a third body. Simultaneously, a visual companion has been detected based on the Gaia astrometric data, indicating that KM UMa is actually in a 2+1+1 hierarchical quadruple system.

Based on observations from the Solar Dynamics Observatory (SDO) and the Solar and Heliospheric Observatory (SOHO), we investigate 30 apparent current sheets during 1999–2021, including 10 on-disk and 8 limb ones from the SDO, as well as 12 limb ones from the SOHO. Each on-disk current sheet is formed among an X-type configuration consisting of two sets of atmospheric structures, and each limb one is involved in a flare–coronal mass ejection event. During magnetic reconnection period, the on-disk apparent current sheet evolves from a bright point to an elongated line-like structure, and the structure becomes thin in the late stage of the reconnection. Subsequently, the plasma distribution within the current sheet manifests as a plasmoid chain. For the limb apparent current sheet, the length elongation is faster than that of the on-disk one, and the thinning process is also detected. Although the aspect ratios of the limb cases are comparable to the value for the occurrence of tearing mode instability from simulation research, no obvious plasmoid chain is detected within these limb current sheets, and the density distribution is locally uniform. We suggest that due to the rapid extension of limb cases, the tearing mode instability is very fast, resulting in the formation of tiny plasmoids that are smaller than the instrument resolution. Moreover, there is another possible scenario. The observed limb apparent current sheet is just a bright ray, and the actual current sheet is only a small segment of the ray.

The origin of our Galaxy’s high-velocity clouds (HVCs) remains a mystery after many decades of effort. In this paper, we use the TNG50 simulation of the IllustrisTNG project to identify cool, dense clouds that match observations of Galactic H i HVCs. We track these clouds back in time to determine their origin. For a TNG50 Milky Way−like galaxy, we find that only 17% of HVCs can be tracked directly to the disk and 21% to material stripped out of satellites. The majority of HVCs (62%) arise from warm and hot circumgalactic gas that cools through thermal instability. They then obtain their anomalous velocities through interactions with the turbulent circumgalactic medium. At TNG50 resolution, we do not see evidence for HVCs forming out of very low metallicity intergalactic material. Instead, low-metallicity HVCs are most likely associated with satellites. These results suggest that Galactic HVCs are highly heterogeneous in their origin and can provide insight into the physical processes that shape the circumgalactic medium, such as disk outflows, satellite accretion, and thermal instabilities.

Massive stars may form in or be captured into active galactic nuclei (AGN) disks. Recent 1D studies employing stellar-evolution codes have demonstrated the potential for rapid growth of such stars through accretion up to a few hundred solar masses. We perform 3D radiation hydrodynamic simulations of moderately massive stars’ envelopes in order to determine the rate and critical radius R_crit of their accretion process in an isotropic gas-rich environment in the absence of luminosity-driven mass loss. We find that in the “fast-diffusion” regime where characteristic radiative diffusion speed c/τ is faster than the gas sound speed c_s, the accretion rate is suppressed by feedback from gravitational and radiative advection energy flux, in addition to the stellar luminosity. Alternatively, in the “slow-diffusion” regime where c/τ < c_s, due to adiabatic accretion, the stellar envelope expands quickly to become hydrostatic and further net accretion occurs on thermal timescales in the absence of self-gravity. When the radiation entropy of the medium is less than that of the star, however, this hydrostatic envelope can become more massive than the star itself. Within this subregime, the self-gravity of the envelope excites runaway growth. Applying our results to realistic environments, moderately massive stars (≲100M_⊙) embedded in AGN disks typically accrete in the fast-diffusion regime, leading to a reduction of steady-state accretion rate 1–2 orders of magnitudes lower than expected by previous 1D calculations and R_crit smaller than the disk scale height, except in the opacity window at temperature T ∼ 2000 K. Accretion in slow diffusion regime occurs in regions with very high density ρ ≳ 10⁻⁹ g cm⁻³, and needs to be treated with caution in 1D long-term calculations.

We use continuous wavelet transform techniques to construct the global and environment-dependent wavelet statistics, such as energy spectrum and kurtosis, to study the fluctuation and intermittency of the turbulent motion in the cosmic fluid velocity field with the IllustrisTNG simulation data. We find that the peak scale of the energy spectrum defines a characteristic scale, which can be regarded as the integral scale of turbulence, and the Nyquist wavenumber can be regarded as the dissipation scale. With these two characteristic scales, the energy spectrum can be divided into the energy-containing range, the inertial range, and the dissipation range of turbulence. The wavelet kurtosis is an increasing function of the wavenumber k, which first grows rapidly then slowly with k, indicating that the cosmic fluid becomes increasingly intermittent with k. In the energy-containing range, the energy spectrum increases significantly from z = 2 to 1, but remains almost unchanged from z = 1 to 0. We find that both the environment-dependent spectrum and kurtosis are similar to the global ones, and the magnitude of the spectrum is smallest in the lowest-density and largest in the highest-density environment, suggesting that the cosmic fluid is more turbulent in a high-density than in a low-density environment. In the inertial range, the energy spectrum’s exponent is steeper than both the Kolmogorov and Burgers exponents, indicating more efficient energy transfer compared to Kolmogorov or Burgers turbulence.

The Vera C. Rubin Observatory is a unique facility for survey astronomy that will soon be commissioned and begin operations. Crucial to many of its scientific goals is the achievement of sustained high image quality, limited only by the seeing at the site. This will be maintained through an active optics system that controls optical element misalignments and corrects mirror figure error to minimize aberrations caused by both thermal and gravitational effects. However, the large number of adjustment degrees of freedom available on the Rubin Observatory introduces a range of degeneracies, including many that are induced by noise due to imperfect measurement of the wave-front errors. We present a structured methodology for identifying these degeneracies through an analysis of image noise level. We also present a novel scaling strategy based on truncated singular value decomposition that mitigates the degeneracy and optimally distributes the adjustment over the available degrees of freedom. Our approach ensures the attainment of optimal image quality, while avoiding excursions around the noise-induced subspace of degeneracies, marking a significant improvement over the previous techniques adopted for Rubin, which were based on an optimal integral controller. This new approach is likely to also yield significant benefits for all telescopes that incorporate large numbers of degrees of freedom of adjustment.

The solar wind dynamics are thought to be governed by counter-propagating Alfvén waves. Alfvén waves generate a turbulent cascade through nonlinear couplings between shearing wave packets. However, imbalances, structures, and intermittency complicate the solar wind’s magnetohydrodynamic (MHD) turbulence. Simulations and theories have suggested that dynamic alignment, or the tendency for velocity and magnetic fluctuations to become more correlated as a function of scale, may dictate the turbulent cascade. Observations have hinted at dynamic alignment at large scales but remain inconclusive at small scales. To investigate the nature of dynamically aligning fluctuations in the solar wind, we examine slow wind intervals from the Helios 2 spacecraft. We develop a two-component model of solar wind turbulence consisting of Alfvénic and non-Alfvénic contributions. We assume that only counter-propagating Alfvén waves experience dynamic alignment, and the non-Alfvénic structures do not participate in the alignment process. Our model allows us to constrain the relative contribution of inward Alfvén waves and non-Alfvénic structures with respect to the amplitude of outward Alfvén waves in the slow solar wind as a function of scale under the assumption of dynamic alignment. We show that the power ratio of inward to outward Alfvén fluctuations decreases as a function of scale. At the same time, the non-Alfvénic structures to outward Alfvén fluctuations increase with decreasing scale. Increasing structures provide a possible explanation for no dynamic alignment at small scales. Our study implies the need for new theoretical models to fully account for the solar wind’s compressibility, intermittency, and imbalanced nature.

We postprocess a three-dimensional, general relativistic, full transport neutrino radiation magnetohydrodynamics simulation of the black-hole-accretion disk-wind system thought to be a potential outcome of the GW170817 merger to investigate the presence of electron lepton number (ELN-XLN) crossings in the neutrino angular distribution. Neutrinos are evolved with an explicit Monte Carlo method and can interact with matter via emission, absorption, or scattering. Within the postprocessing framework, we find ubiquitous occurrence of ELN-XLN crossings at early times (∼11 ms), but this does not hold for later times in the simulation. At postmerger times of ∼60 ms and beyond, ELN-XLN crossings are only present near the equator. We provide a detailed analysis of the neutrino radiation field to investigate the origin and time evolution of these crossings. Previous reports have suggested ubiquitous flavor crossings persisting throughout the simulation lifetime, albeit for different sets of conditions for the merger remnant, the treatment of hydrodynamics, and neutrino transport. Even though we do not perform a direct comparison with other published works, we qualitatively assess the reasons for the difference with our results. The geometric structure and evolution of the ELN-XLN crossings found in our analysis, and by extension, fast flavor instabilities, have important implications for heavy element nucleosynthesis in neutron star mergers.

We present a multifrequency polarimetric study for the quasar 1604+159. The source was observed at the L band with the American Very Long Baseline Array and the L, X, and U bands with the Very Large Array. These observations provide different resolutions from mas to arcsec, enabling us to probe the morphology and magnetic field from tens of parsec to hundreds of kiloparsec scale. We detect a symmetrical Fanaroff–Riley Class I–like structure. The source has several lobes and bulges, forming a cocoon shape. The polarization is normal to the edges of the structure with high fractional polarization up to ∼60%. Two hotspots are observed at the eastern and western sides of the source, located symmetrically relative to the core. The flux density ratio (>1.5) between the two hotspots suggests the Doppler beaming effect exists at a large scale. The polarized emission in the hotspots also shows a symmetrical structure with an oblique direction from the jet direction. In general, the jet propagates in a collimating structure with several bends. Polarization is also detected perpendicular to the local jet from ∼100 mas to ∼1″. The jet shows strong polarized intensity and high fractional polarization at the bending edges. We discuss the possible origins of the observed structure and magnetic field.

We investigate the bright CO fundamental emission in the central regions of five protostars in their primary mass assembly phase using new observations from JWST’s Near-Infrared Spectrograph and Mid-Infrared Instrument. CO line emission images and fluxes are extracted for a forest of ∼150 rovibrational transitions from two vibrational bands, v = 1−0 and v = 2−1. However, ¹³CO is undetected, indicating that ¹²CO emission is optically thin. We use H₂ emission lines to correct fluxes for extinction and then construct rotation diagrams for the CO lines with the highest spectral resolution and sensitivity to estimate rotational temperatures and numbers of CO molecules. Two distinct rotational temperature components are required for v = 1 (∼600 to 1000 K and 2000 to ∼10⁴ K), while one hotter component is required for v = 2 (≳3500 K). ¹³CO is depleted compared to the abundances found in the interstellar medium, indicating selective UV photodissociation of ¹³CO; therefore, UV radiative pumping may explain the higher rotational temperatures in v = 2. The average vibrational temperature is ∼1000 K for our sources and is similar to the lowest rotational temperature components. Using the measured rotational and vibrational temperatures to infer a total number of CO molecules, we find that the total gas masses range from lower limits of ∼10²² g for the lowest mass protostars to ∼10²⁶ g for the highest mass protostars. Our gas mass lower limits are compatible with those in more evolved systems, which suggest the lowest rotational temperature component comes from the inner disk, scattered into our line of sight, but we also cannot exclude the contribution to the CO emission from disk winds for higher mass targets.

Massive evolved stars such as red supergiants and hypergiants are potential progenitors of Type II supernovae, and they are known for ejecting substantial amounts of matter, up to half their initial mass, during their final evolutionary phases. The rate and mechanism of this mass loss play a crucial role in determining their ultimate fate and the likelihood of their progression to supernovae. However, the exact mechanisms driving this mass ejection have long been a subject of research. Recent observations, such as the Great Dimming of Betelgeuse, have suggested that the activity of large convective cells, combined with pulsation, could be a plausible explanation for such mass-loss events. In this context, we conducted interferometric observations of the famous yellow hypergiant, ρ Cassiopeiae using the CHARA Array in H- and K-band wavelengths. ρ Cas is well known for its recurrent eruptions, characterized by periods of visual dimming (∼1.5–2 mag) followed by recovery. From our observations, we derived the diameter of the limb-darkened disk and found that this star has a radius of 1.04 ± 0.01 mas, or 564–700 R_⊙. We performed image reconstructions with three different image reconstruction software packages, and they unveiled the presence of giant hot and cold spots on the stellar surface. We interpret these prominent hot spots as giant convection cells, suggesting a possible connection to mass ejections from the star’s envelope. Furthermore, we detected spectral CO emission lines in the K band (λ = 2.31–2.38 μm), and the image reconstructions in these spectral lines revealed an extended circumstellar envelope with a radius of 1.45 ± 0.10 mas.

In the upcoming gravitational-wave (GW) observing runs, identifying host galaxies is crucial as it provides essential redshift information and enables the use of GW events as standard sirens. However, pinpointing host galaxies remains challenging due to the large localization uncertainties and the rapidly fading nature of their optical counterparts. Analyzing the host galaxies of short gamma-ray bursts (sGRBs) offers an alternative approach to deepen our understanding of the environments where binary neutron stars (BNS) primarily merge. This study compiles archival photometric data for the host galaxies of 76 sGRBs and four hybrid GRBs that are long GRBs accompanied by kilonova-like signals. We use this data to evaluate their physical properties through spectral energy distribution fitting. To assess the characteristics of the host galaxies, we utilized a volume-limited sample (z < 0.5) from the COSMOS field as a control group. Contrary to expectations that the BNS merger rate is proportional to host stellar mass, the short and hybrid GRB population appears less massive than the mass-weighted distribution of the control sample. Instead, we propose a formulation for the expected BNS merger rate from a galaxy as + , which optimally explains the deviation between the stellar mass distributions of the GRB host galaxies and the control sample. These insights provide a strategic framework for targeted GW follow-ups and enhance our ability to identify potential host galaxies for future GW events.

One of the most serious limitations of current astrochemical models with the rate equation (RE) approach is that only a single type of binding site is considered in grain surface chemistry, although laboratory and quantum chemical studies have found that surfaces contain various binding sites with different potential energy depths. When various sites exist, adsorbed species can be trapped in deep potential sites, increasing the resident time on the surface. On the other hand, adsorbed species can be populated in shallow sites, activating thermal hopping and thus two-body reactions even at low temperatures, where the thermal hopping from deeper sites is not activated. Such behavior cannot be described by the conventional RE approach. In this work, I present a framework for incorporating various binding sites (i.e., binding energy distribution) in gas-ice astrochemical models as an extension of the conventional RE approach. I propose a simple method to estimate the probability density function (pdf) for the occupation of various sites by adsorbed species, assuming a quasi-steady state. By using thermal desorption and hopping rates weighted by the pdfs, the effect of binding energy distribution is incorporated into the RE approach without increasing the number of ordinary differential equations to be solved. This method is found to be accurate and computationally efficient, and enables us to consider binding energy distribution even for a large gas-ice chemical network which contains hundreds of icy species. The impact of the binding energy distribution on interstellar ice composition is discussed quantitatively for the first time.

Event Horizon Telescope (EHT) images of the horizon-scale emission around the Galactic center supermassive black hole Sagittarius A* (Sgr A*) favor accretion flow models with a jet component. However, this jet has not been conclusively detected. Using the “best-bet” models of Sgr A* from the EHT Collaboration, we assess whether this nondetection is expected for current facilities and explore the prospects of detecting a jet with very-long-baseline interferometry (VLBI) at four frequencies: 86, 115, 230, and 345 GHz. We produce synthetic image reconstructions for current and next-generation VLBI arrays at these frequencies that include the effects of interstellar scattering, optical depth, and time variability. We find that no existing VLBI arrays are expected to detect the jet in these best-bet models, consistent with observations to date. We show that next-generation VLBI arrays at 86 and 115 GHz—in particular, the EHT after upgrades through the ngEHT program and the ngVLA—successfully capture the jet in our tests due to improvements in instrument sensitivity and (u, v) coverage at spatial scales critical to jet detection. These results highlight the potential of enhanced VLBI capabilities in the coming decade to reveal the crucial properties of Sgr A* and its interaction with the Galactic center environment.

We present results from our deep far-ultraviolet (FUV) survey using AstroSat/UVIT of a filamentary structure at z ∼0.072. A total of four filaments comprising 58 galaxies were probed in our study. We detect 18 filament galaxies in our FUV observation. All filament galaxies are further classified based on their photometric color, nuclear activity, and morphology. The filaments contain galaxies with mixed stellar population types and structures. We do not detect galaxies in our UVIT survey up to a distance of 0.4 Mpc h⁻¹ from the filament axis, implying a lack of recent star formation in the inner region of filaments. The FUV star formation rate (SFR) for star-forming galaxies agrees well with the SFR_144MHz calculated using Low-Frequency Array radio-continuum observations. We witness an increase in the FUV specific-SFR (sSFR) of filament galaxies with increasing distance from the filament spine (D_fil). The intermediate-to-high stellar mass filament galaxies were more star-forming than cluster galaxies in a fixed stellar mass bin. The FUV morphology of some filament galaxies detected in the filament outskirts (D_fil ≳ 0.7 Mpc h⁻¹) is comparable to or slightly extended than their optical counterpart. The mass assembly of galaxies examined by estimating (FUV − r) color gradients shows that more “red-cored’ galaxies reside in the outer region of the filaments. Our results prove that the likelihood of merger interaction and gas starvation increases when approaching the filament spine. We report a definitive and inhomogeneous impact of filaments on the galaxies residing inside them.

Polarization of starlight and thermal dust emission caused by aligned dust grains is a valuable tool to characterize magnetic fields (B-fields) and constrain dust properties. However, the grain alignment physics is not fully understood. To test the radiative torque (RAT) paradigm, including RAT alignment (RAT-A) and disruption (RAT-D), we use dust polarization observed by Planck and Stratosphere Observatory for Infrared Astronomy/High-resolution Airborne Wideband Camera Plus toward two filaments, Musca and OMC-1, with contrasting physical conditions. Musca, a quiescent filament, is ideal for testing RAT-A, while OMC-1, an active star-forming region OMC-1, is most suitable for testing RAT-D. We found that polarization fraction, P, decreases with increasing polarization angle dispersion function, , and increasing column density, N(H₂), consistent with RAT-A. However, P increases with increasing dust temperature, T_d, but decreases when T_d reaches a certain high value. We compute the polarization fraction for the ideal models with B-fields in the plane of sky based on the RAT paradigm, accounting for the depolarization effect by B-field tangling. We then compare the realistic polarization model with observations. For Musca with well-ordered B-fields, our numerical model successfully reproduces the decline of P toward the filament spine (aka polarization hole), having higher N(H₂) and lower T_d, indicating the loss of grain alignment efficiency due to RAT-A. For OMC-1, with stronger B-field variations and higher temperatures, our model can reproduce the observed P–T_d and P−N(H₂) relations only if the B-field tangling and RAT-D effect are incorporated. Our results provide more robust observational evidence for the RAT paradigm, particularly the recently discovered RAT-D.

X-ray observations provide important and valuable insights into the acceleration and propagation of nonthermal electrons during solar flares. Improved X-ray spectral analysis requires a deeper understanding of the dynamics of energetic electrons. Previous studies have demonstrated that the dynamics of accelerated electrons with a few thermal speeds are more complex than those with significantly higher speeds. To better describe the energetic electrons after injection, a model considering energy diffusion and thermalization effects in flare conditions (the warm-target model) has recently been developed for spectral analysis of hard X-rays. This model has demonstrated how the low-energy cutoff, which can hardly be constrained in cold-target modeling, can be determined. However, the power-law form may not be the most suitable representation of injected electrons. The kappa distribution, which is proposed as a physical consequence of electron acceleration, has been applied successfully in RHESSI spectral analysis. In this study, we employ the kappa-form injected electrons in the warm-target model to analyze two M-class flares, observed by RHESSI and the Spectrometer/Telescope for Imaging X-rays, respectively. The best-fit results show that the kappa-form energetic electron spectrum generates lower nonthermal energy than the power-law form when producing a similar photon spectrum in the fit range. We also demonstrated that the fit parameters associated with the kappa-form electron spectrum can be well determined with small uncertainty. Further, the kappa distribution, which covers the entire electron energy range, enables the determination of key electron properties such as total electron number density and average energy at the flare site, providing valuable information on electron acceleration processes.

Using an established classification technique, we leverage standard observations and analyses to predict the progenitors of gamma-ray bursts (GRBs). This technique, utilizing support vector machine statistics, provides a more nuanced prediction than the previous two-component Gaussian mixture in duration of the prompt gamma-ray emission. Based on further covariance testing from Fermi/Gamma Ray Burst Monitor, Swift/Burst Alert Telescope, and Swift/X-Ray Telescope data, we find that our classification based only on prompt emission properties gives perspective on the recent evidence that mergers and collapsars exist in both “long” and “short” GRB populations.

We present an analysis searching for dual active galactic nuclei (AGN) among 62 high-redshift (2.5 < z < 3.5) X-ray sources selected from the X-UDS, AEGIS-XD, CDF-S, and COSMOS-Legacy Chandra surveys. We aim to quantify the frequency of dual AGN in the high-redshift Universe, which holds implications for black hole merger timescales and low-frequency gravitational wave detection rates. We analyze each X-ray source using BAYMAX, an analysis tool that calculates the Bayes factor for whether a given archival Chandra AGN is more likely a single or dual point source. We find no strong evidence for dual AGN in any individual source in our sample. We increase our sensitivity to search for dual AGN across the sample by comparing our measured distribution of Bayes factors to that expected from a sample composed entirely of single point sources and find no evidence for dual AGN in the sample distribution. Although our analysis utilizes one of the largest Chandra catalogs of high-z X-ray point sources available to study, the findings remain limited by the modest number of sources observed at the highest spatial resolution with Chandra and the typical count rates of the detected sources. Our nondetection allows us to place an upper limit on the X-ray dual AGN fraction at 2.5 < z < 3.5 of 4.8% at the 95% confidence level. Expanding substantially on these results at X-ray wavelengths will require future surveys spanning larger sky areas and extending to fainter fluxes than has been possible with Chandra. We illustrate the potential of the AXIS mission concept in this regard.

Over the last decade, studies of large samples of binary systems have identified chemical anomalies and shown that they might be attributed to planet formation or planet engulfment. However, both scenarios have primarily been tested in pairs without known exoplanets. In this work, we explore these scenarios in the newly detected planet-hosting wide binary TOI-1173 A/B (projected separation ∼11,400 au), using high-resolution MAROON-X and ARCES spectra. We determined photospheric stellar parameters both by fitting stellar models and via the spectroscopic equilibrium approach. Both analyses agree and suggest that they are cool main-sequence stars located in the thin disk. A line-by-line differential analysis between the components (B−A) displays an abundance pattern in the condensation temperature plane, where the planet-hosting star TOI-1173 A is enhanced in refractory elements such as iron by more than 0.05 dex. This suggests the engulfment of ∼18 M_⊕ of rocky material in star A. Our hypothesis is supported by the dynamics of the system (detailed in our companion paper), which suggest that the super-Neptune TOI-1173 A b might have been delivered to its current short period (∼7 days) through circularization and von Zeipel–Lidov–Kozai mechanisms, thereby triggering the engulfment of inner rocky exoplanets.

Minifilament eruptions producing small jets and microflares have mostly been studied based on coronal observations at extreme-ultraviolet and X-ray wavelengths. This study presents chromospheric plasma diagnostics of a quiet-Sun minifilament of size ∼ 2″ × 5″ with a sigmoidal shape and an associated microflare observed on 2021 August 7 17:00 UT using high temporal and spatial resolution spectroscopy from the Fast Imaging Solar Spectrograph (FISS) and high-resolution magnetograms from the Near InfraRed Imaging Spectropolarimeter (NIRIS) installed on the 1.6 m Goode Solar Telescope at Big Bear Solar Observatory. Using FISS Hα and Ca ii 8542 Å line spectra at the time of the minifilament activation we determined a temperature of 8600 K and a nonthermal speed of 7.9 km s⁻¹. During the eruption, the minifilament was no longer visible in the Ca ii 8542 Å line, and only the Hα line spectra were used to find the temperature of the minifilament, which reached 1.2 × 10⁴ K and decreased afterward. We estimated thermal energy of 3.6 × 10²⁴ erg from the maximum temperature and kinetic energy of 2.6 × 10²⁴ erg from the rising speed (18 km s⁻¹) of the minifilament. From the NIRIS magnetograms we found small-scale flux emergence and cancellation coincident with the minifilament eruption, and the magnetic energy change across the conjugate footpoints reaches 7.2 × 10²⁵ erg. Such spectroscopic diagnostics of the chromospheric minifilament complement earlier studies of minifilament eruptions made using coronal images.

Giant eruptions (GEs) in luminous blue variables are years-to-decades-long episodes of enhanced mass loss from the outer layers of the star during which the star undergoes major changes in its physical and observed properties. We use the Modules for Experiments in Stellar Astrophysics stellar evolution code to model the evolution of a 70 M_⊙ star that undergoes a GE. We let the star evolve to the termination of the main sequence, and when it reaches T ≃ 19,400 K we emulate a GE by removing mass from its outer layers at a rate of 0.15 M_⊙ yr⁻¹ for 20 yr. As mass is being lost, the star contracts and releases a substantial amount of gravitational energy. The star undergoes an initial ≃3 days of expansion followed by years of contraction. During that time the star tries to reach an equilibrium state, and as a result of loss in gravitational energy, its luminosity drops about 1 order of magnitude. As the GE terminates, we let the star continue to evolve without any further mass loss and track its recovery as it regains its equilibrium by adjusting its internal structure. After ≃87 yr it reaches a state very close to the one where the GE was first initiated. We suggest that at this point another GE or a cycle of GEs may occur.

Magnetar outbursts are powered by an intense magnetic field. The phenomenon has recently drawn significant attention because of a connection to some fast radio bursts that has been reported. Understanding magnetar outbursts may provide the key to mysterious transient events. The elastic deformation of the solid crust due to magnetic field evolution accumulates over a secular timescale. Eventually, the crust fractures or responds plastically beyond a particular threshold. Determination of the critical limit is required to obtain the shear strain tensor in response to magnetic stress. In some studies, the tensor was substituted with an approximate expression determined algebraically from the magnetic stress. This study evaluated the validity of the approximation by comparing it with the strain tensor obtained through appropriate calculations. The differential equations for the elastic deformation driven by the magnetic field were solved. The results indicated that the approximation did not represent the correct strain tensor value, in both magnitude and spatial profile. Previous evolutionary calculations based on spurious criteria are likely to overestimate the magnitude of the strain tensor, and crustal failure occurs on a shorter timescale. Therefore, revisiting evolutionary calculations using the correct approach is necessary. This study is essential for developing the dynamics of crustal fractures and the magnetic field evolution in a magnetar.

The statistical analysis of cosmic large-scale structure is most often based on simple two-point summary statistics, like the power spectrum or the two-point correlation function of a sample of galaxies or other types of tracers. In contrast, topological measures of clustering are also sensitive to higher-order correlations and thus offer the prospect to access additional information that may harbor important constraining power. We here revisit one such geometric measure of the cosmic web in the form of the so-called percolation analysis, using the recent MillenniumTNG simulation suite of the ΛCDM paradigm. We analyze continuum percolation statistics both for high-resolution dark matter particle distributions and for galaxy mock catalogs from a semianalytic galaxy formation model within a periodic simulation volume of 3000 Mpc on a side. For comparison, we also investigate the percolation statistics of random particle sets and neutrino distributions with two different summed particle masses. We find that the percolation statistics of the dark matter distribution evolves strongly with redshift and thus clustering strength, yielding a progressively lower percolation threshold toward later times. However, there is a sizable residual dependence on numerical resolution, which we interpret as a residual influence of different levels of shot noise. This is corroborated by our analysis of galaxy mock catalogs, whose results depend on sampling density more strongly than on galaxy selection criteria. While this limits the discriminative power of percolation statistics, our results suggest that it still remains useful as a complementary cosmological test when controlled for sampling density.

We use the James Webb Space Telescope Mid-Infrared Instrument medium-resolution spectrometer observations of the radio-loud active galactic nucleus (AGN) host UGC 8782 to map the warm molecular and ionized gas kinematics. The data reveal spatially resolved outflows in the inner 2 kpc, seen in low ionization (traced by the [Ar ii] 6.99 μm emission) and in warm molecular gas (traced by the H₂ rotational transitions). We find a maximum mass-outflow rate of 4.90 ± 2.04 M_⊙ yr⁻¹ at ∼900 pc from the nucleus for the warm outflow (198 K ≤ T ≤ 1000 K) and estimate an outflow rate of up to 1.22 ± 0.51 M_⊙ yr⁻¹ for the hotter gas phase (T > 1000 K). These outflows can clear the entire nuclear reservoir of warm molecular gas in about 1 Myr. The derived kinetic power of the molecular outflows leads to coupling efficiencies of 2%–5% of the AGN luminosity, way above the minimum expected for the AGN feedback to be effective in quenching the star formation.

The supermassive black holes (SMBHs) that power active galactic nuclei found at z ≥ 6 were formed during the Epoch of Reionization. Because reionization is an inhomogeneous process, the physical properties of SMBH host-galaxy environments will vary spatially during reionization. We construct a semi-analytic model to estimate the impact of reionization on SMBH growth. Using a series of merger trees, reionization models, and black hole (BH) growth models, we find that early reionization can reduce a SMBH’s mass by up to [50, 70, 90]% within dark matter halos of mass [10¹², 10¹¹, 10¹⁰] M_⊙ by z = 6. Our findings also suggest that the redshift range in which BH growth is impacted by reionization strongly depends on whether the Eddington accretion rate can be exceeded. If so, we find that BH masses are significantly suppressed principally during the early phases of reionization (z ≳ 10), while they are more readily suppressed across the full redshift range if super-Eddington growth is not allowed. We find that the global average impact of reionization may be to reduce the masses of BHs residing in ≲10¹¹M_⊙ halos by a factor of ≳2. The census of SMBHs being uncovered by the JWST may offer a means to test the basic prediction that more massive BHs reside in cosmological volumes that are reionized at later times.

In recent years, the continued detection of complex organic molecules of prebiotic interest has refueled the interest in a panspermic origin of life. The prebiotic molecule glyceraldehyde is proposed to be formed from (Z)-1,2-ethenediol, a molecule recently detected toward the G+0.693-0.027 molecular cloud at the galactic center. In this work, we computationally simulate the formation of (Z)-1,2-ethenediol from vinyl alcohol on the surface of amorphous solid water in a two-step synthesis involving an OH addition and an H abstraction reaction. In total, we considered all reaction possibilities of the 1,1- and 1,2-OH addition to vinyl alcohol followed by H abstraction or H addition reactions on the resulting radicals. The combination of these reactions is capable of explaining the formation of (Z)-1,2-ethenediol provided a suprathermal diffusion of OH. We also conclude that our proposed formation pathway is not selective and also yields other abstraction and addition products. Key in our findings is the connection between the adsorption modes of the reactants and intermediates and the stereoselectivity of the reactions.

We present a photometric study of EX Dra, a dwarf nova that has been extensively observed by the Transiting Exoplanet Survey Satellite. The data reveal the occurrence of 20 complete outbursts, exhibiting several intriguing and rare characteristics. The light curves exhibit a distinct superorbital signal with a period of approximately P_sor ∼ 4.39(7) days, along with a negative superhump showing an approximate period of P_nsh ∼ 4.805(1) hr, indicating that the accretion disk is tilted and undergoing precession with the period of P_sor. In addition, the time-varying nature of P_sor suggests that the precession period is fluctuating.The eclipsing light minima O – C analysis during quiescence shows an oscillation with period of 3.9(5) days, which is a little shorter than the superorbital period. We contend that this is unlikely to be a sudden alteration of the orbital period, but rather, it is influenced by the tilt and precession of the accretion disk. Notably, we found an amplitude shift in the outburst behavior from 3.5 mag with a periodicity of about 26 days to an amplitude of around 2.5 mag with a periodicity of about 12 days, which persisted for 14 yr before reverting. Furthermore, we have extracted quasiperiodic oscillations in the plateau at the noneclipsed phases, characterized by periods ranging between 37 and 40 minutes.

The prediction of solar energetic particle (SEP) events garners increasing interest as space missions extend beyond Earth’s protective magnetosphere. These events, which are, in most cases, products of magnetic-reconnection-driven processes during solar flares or fast coronal-mass-ejection-driven shock waves, pose significant radiation hazards to aviation, space-based electronics, and particularly space exploration. In this work, we utilize the recently developed data set that combines the Solar Dynamics Observatory/Space-weather Helioseismic and Magnetic Imager Active Region Patches and the Solar and Heliospheric Observatory/Space-weather Michelson Doppler Imager Active Region Patches. We employ a suite of machine learning strategies, including support vector machines (SVMs) and regression models, to evaluate the predictive potential of this new data product for a forecast of post-solar flare SEP events. Our study indicates that despite the augmented volume of data, the prediction accuracy reaches 0.7 ± 0.1 (experimental setting), which aligns with but does not exceed these published benchmarks. A linear SVM model with training and testing configurations that mimic an operational setting (positive–negative imbalance) reveals a slight increase (+0.04 ± 0.05) in the accuracy of a 14 hr SEP forecast compared to previous studies. This outcome emphasizes the imperative for more sophisticated, physics-informed models to better understand the underlying processes leading to SEP events.

Tilted disk precession exists in different objects. Negative superhumps (NSHs) in cataclysmic variable stars are believed to arise from the interaction between the reverse precession of a tilted disk and the streams from the secondary star. Utilizing Transiting Exoplanet Survey Satellite photometry, we present a comprehensive investigation into the tilted disk precession and NSHs in the dwarf nova (DN) HS 2325+8205, employing eclipse minima, eclipse depths, NSH frequencies, and NSH amplitudes and the correlation between them as the windows. We identified NSHs with a period of 0.185671(17) day in HS 2325+8205. The NSH frequency exhibits variability with a period of 3.943(9) days, akin to the tilted disk precession period validated in nova-like stars (SDSS J0812) and intermediate polars (IPs; TV Col). The O − C of the eclipse minima were similarly found to vary cyclically in a period of 4.135(5) days, characterized by a faster rise than fall. Furthermore, the NSH amplitude exhibits complex and diverse variations, which may be linked to changes in the disk radius, the mass transfer rate, and the apparent area of the hot spot. For the first time in DNe, we observe biperiodic variations in eclipse depth (P₁ = 4.131(4) days and P₂ = 2.065(2) days ≈ P_prec/2) resembling those seen in IPs, suggesting that variations with P₂ are not attributable to an accretion curtain, as previously suspected. Moreover, NSH amplitude and eclipse depth decrease with increasing NSH frequency, while NSH amplitude correlates positively with eclipse depth. These complex variations observed across multiple observational windows provide substantial evidence for the understanding of tilted disk precession and NSHs.

Recently, the detection of a coherent radio flash associated with short-duration GRB 201006A, occurring 76.6 minutes after the burst, has attracted great attention. However, the physical origin of the coherent radio flash remains under debate. By reanalyzing its data observed by Fermi and Swift, we find that an early radio afterglow as the physical origin of the radio flash can be ruled out, but the coherent radio emission seems to be consistent with the hypothesis of a supramassive magnetar as the central engine collapsing into a black hole. Within this scenario, the derived magnetar surface magnetic field (B_p) and the initial spin period (P₀) fall into a reasonable range but require a preferably low value of η_R = 10⁻⁷ or 10⁻⁶. Moreover, the calculated low-ε value and E_γ,iso–E_p correlation of GRB 201006A also supports the progenitor which is from the merger of compact stars. We also discuss the non-detected kilonova emission associated with GRB 201006A, and then compare with its upper limits of optical observations.

NGC 4278, a low-luminosity active galactic nucleus, is generally classified as a low-ionization nuclear emission-line region (LINER). Recently, it has been reported to be associated with a very high energy γ-ray source 1LHAASO J1219+2915 in the first Large High Altitude Air Shower Observatory source catalog. However, no associated counterpart has been detected by analyzing the data collected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. By analyzing its X-ray observation data from Swift-XRT, we find that NGC 4278 is in a high-flux state on MJD 59546, with the X-ray flux more than one order of magnitude higher than that observed ∼11.7 yr earlier by Chandra. Interestingly, this Swift-XRT observation was conducted during the active phase of the γ-ray source 1LHAASO J1219+2915. We propose that the detection of very high energy γ-rays from NGC 4278 may be attributed to the presence of an active nucleus in its center. To reproduce the spectral energy distribution (SED) of NGC 4278, we employ a one-zone leptonic model, typically used for fitting broadband SEDs of BL Lacs, and find that a smaller magnetic field strength is required than that of typical TeV BL Lacs. Furthermore, NGC 4278 exhibits significantly lower luminosity in both radio and TeV bands when compared with typical TeV BL Lacs. In the radio luminosity versus Eddington ratio plane, NGC 4278 shows greater similarity to Seyfert galaxies and LINERs than to BL Lacs; however, it still roughly follows the extension toward lower luminosity seen in BL Lacs.

We present a spatially resolved study of stellar populations in six galaxies with stellar masses M_* ∼ 10¹⁰M_☉ at z ∼ 3.7 using 14-filter James Webb Space Telescope (JWST)/NIRCam imaging from the JADES and JEMS surveys. The six galaxies are visually selected to have clumpy substructures with distinct colors over rest frame 3600−4100 Å, including a red, dominant stellar core that is close to their stellar-light centroids. With 23-filter photometry from the Hubble Space Telescope to JWST, we measure the stellar-population properties of individual structural components via spectral energy distribution fitting using Prospector. We find that the central stellar cores are ≳2 times more massive than the Toomre mass, indicating they may not form via single in situ fragmentation. The stellar cores have stellar ages of 0.4−0.7 Gyr that are similar to the timescale of clump inward migration due to dynamical friction, suggesting that they likely instead formed through the coalescence of giant stellar clumps. While they have not yet quenched, the six galaxies are below the star-forming main sequence by 0.2−0.7 dex. Within each galaxy, we find that the specific star formation rate is lower in the central stellar core, and the stellar-mass surface density of the core is already similar to quenched galaxies of the same masses and redshifts. Meanwhile, the stellar ages of the cores are either comparable to or younger than the extended, smooth parts of the galaxies. Our findings are consistent with model predictions of the gas-rich compaction scenario for the buildup of galaxies’ central regions at high redshifts. We are likely witnessing the coeval formation of dense central cores, along with the onset of galaxy-wide quenching at z > 3.

Stars form within dense cores composed of both gas and dust within molecular clouds. However, despite the crucial role that dust plays in the star formation process, its dynamics is frequently overlooked, with the common assumption being a constant, spatially uniform dust-to-gas ratio and grain size spectrum. In this study, we introduce a set of radiation-dust-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project. These simulations expand upon the earlier radiation MHD models, which included cooling, individual star formation, and feedback. Notably, they explicitly address the dynamics of dust grains, considering radiation, drag, and Lorentz forces acting on a diverse size spectrum of live dust grains. We find that once stars exceed a certain mass threshold (∼2 M_⊙), their emitted radiation can evacuate dust grains from their vicinity, giving rise to a dust-suppressed zone of size ∼100 au. This removal of dust, which interacts with gas through cooling, chemistry, drag, and radiative transfer, alters the gas properties in the region. Commencing during the early accretion stages and preceding the main-sequence phase, this process results in a mass-dependent depletion in the accreted dust-to-gas (ADG) mass ratio within both the circumstellar disk and the star. We predict that massive stars (≳10 M_⊙) would exhibit ADG ratios that are approximately 1 order of magnitude lower than that of their parent clouds. Consequently, stars, their disks, and circumstellar environments would display notable deviations in the abundances of elements commonly associated with dust grains, such as carbon and oxygen.

Cosmic dawn represents a critical juncture in cosmic history when the first population of stars emerged. The astrophysical processes that govern this transformation need to be better understood. The detection of redshifted 21 cm radiation emitted from neutral hydrogen during this era offers a direct window into the thermal and ionization state of the Universe. This emission manifests as differential brightness between spin temperature and the cosmic microwave background. The SARAS experiment aims to detect the sky-averaged signal in the frequency range 40–200 MHz. SARAS’s unique design and operational strategy to float the antenna over a water body minimizes spectral features that may arise due to stratified ground beneath the antenna. However, the antenna environment can be prone to configuration changes due to variations in critical design parameters such as conductivity and antenna tilt. In this paper, we connect the variations in antenna properties to signal detection prospects. By using realistic simulations of a direction- and frequency-dependent radiation pattern of the SARAS antenna and its transfer function, we establish critical parameters and estimate bias in the detectability of different models of the global 21 cm signal. We find a correlation between the nature of chromaticity in antenna properties and the bias in the recovered spectral profiles of 21 cm signals. We also find stringent requirements for transfer function corrections, which can otherwise make detection prospects prohibitive. We finally explore a range of critical parameters that allow robust signal detection.

Intrinsic colors (ICs) of stars are essential for studies on both stellar physics and dust reddening. In this work, we developed an XGBoost model to predict ICs with the atmospheric parameters T_eff, , and [M/H]. The model was trained and tested for three colors at Gaia and Two Micron All Sky Survey bands with 1,040,446 low-reddening sources. The atmospheric parameters were determined by the Gaia DR3 General Stellar Parameterizer from Photometry (GSP-phot) module and were validated by comparing with the Apache Point Observatory Galactic Evolution Experiment and Large Sky Area Multi-Object fiber Spectroscopic Telescope. We further confirmed that the biases in GSP-phot parameters, especially for [M/H], do not present a significant impact on the IC prediction. The generalization error of the model estimated by the test set is 0.014 mag for , 0.050 mag for , and 0.040 mag for . The model was applied to a sample containing 5,714,528 reddened stars with stellar parameters from R. Andrae et al. to calculate ICs and reddenings. The high consistency in the comparison of E(J − K_S) between our results and literature values further validates the accuracy of the XGBoost model. The variation of E(G_BP − K_S)/E(G_BP − G_RP), a representation of the extinction law, with Galactic longitude is found on large scales. This work preliminarily presents the feasibility and accuracy of the machine learning approach for IC and dust reddening calculation, whose products could be widely applied to spectrophotometric data. The data sets and trained model can be accessed via Zenodo, doi:10.5281/zenodo.12787594. Models for more bands will be completed in following work.

Interpretation of the ongoing efforts to simulate the atmospheres of potentially habitable terrestrial exoplanets requires that we understand the underlying dynamics and chemistry of such objects to a much greater degree than 1D or even simple 3D models enable. Here, for the tidally locked habitable-zone planet TRAPPIST-1e, we explore one effect which can shape the dynamics and chemistry of terrestrial planets: the inclusion of an Earth-like land–ocean distribution with orography. To do this we use the Earth-system model WACCM6/CESM2 to run a pair of TRAPPIST-1e models with N₂–O₂ atmospheres and with the substellar point fixed over either land or ocean. The presence of orography shapes atmospheric transport, and in the case of Earth-like orography, breaks the symmetry between the Northern and Southern Hemispheres which was previously found in slab ocean models. For example, peak zonal jet speeds in the Southern Hemisphere are 50%–100% faster than similar jets in the Northern Hemisphere. This also affects the meridional circulation, transporting equatorial material toward the south pole. As a result we also find significant changes in the atmospheric chemistry, including the accumulation of potentially lethal quantities of ozone at both the south pole and the surface. Future studies which investigate the effects of landmass distribution on the dynamics of exoplanetary atmospheres should pay close attention to both the dayside land fraction as well as the orography of the land. Simply modeling a flat landmass will not give a complete picture of its dynamical impact.

We analyze the interaction between an interplanetary coronal mass ejection (ICME) detected in situ at the L1 Lagrange point on 2016 October 12 with a trailing high-speed stream (HSS). We aim to estimate the region in the interplanetary (IP) space where the interaction happened/started using a combined observational-modeling approach. We use minimum variance analysis (MVA) and the Walen test to analyze possible reconnection exhaust at the interface of ICME and HSS. We perform a graduated cylindrical shell reconstruction of the CME to estimate the geometry and source location of the CME. Finally, we use a two-step drag-based model (DBM) model to estimate the region in IP space where the interaction took place. The magnetic obstacle observed in situ shows a fairly symmetric and undisturbed structure and shows the magnetic flux, helicity, and expansion profile/speed of a typical ICME. The MVA together with the Walen test, however, confirms reconnection exhaust at the ICME–HSS boundary. Thus, in situ signatures are in favor of a scenario where the interaction is fairly recent. The trailing HSS shows a distinct velocity profile which first reaches a semi-saturated plateau with an average velocity of 500 km s⁻¹ and then saturates at a maximum speed of 710 km s⁻¹. We find that the HSS's interaction with the ICME is influenced only by this initial plateau. The results of the two-step DBM suggest that the ICME has started interacting with the HSS close to Earth (∼0.81 au), which compares well with the deductions from in situ signatures.

The Minkowski functionals (MFs), a set of topological summary statistics, have emerged as a powerful tool for extracting non-Gaussian information. We investigate the prospect of constraining the reionization parameters using the MFs of the 21 cm brightness temperature field from the epoch of reionization (EOR). Realistic effects, including thermal noise, synthesized beam, and foreground avoidance, are applied to the mock observations from radio interferometric array experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the Square Kilometre Array (SKA). We demonstrate that the MFs of the 21 cm signal measured with SKA-Low can be used to distinguish different reionization models, whereas the MF measurement with a HERA-like array cannot be made accurately enough. We further forecast the accuracies with which the MF measurements can place constraints on reionization parameters, using the standard Markov Chain Monte Carlo analysis for parameter inference based on forward modeling. We find that for SKA-Low observation, MFs provide unbiased estimations of the reionization parameters with accuracies comparable to the power spectrum (PS) analysis. Furthermore, joint constraints using both MFs and PS can improve the constraint accuracies by up to 30% compared to those with the PS alone. Nevertheless, the constraint accuracies can be degraded if the EOR window is shrunk with strong foreground avoidance. Our analysis demonstrates the promise of MFs as a set of summary statistics that extract complementary information from the 21 cm EOR field to the two-point statistics, which suggests a strong motivation for incorporating the MFs into the data analysis of future 21 cm observations.

The circumgalactic medium (CGM) of star-forming dwarf galaxies plays a key role in regulating the galactic baryonic cycle. We investigate how susceptible the CGM of dwarf satellite galaxies is to ram pressure stripping in Milky Way–like environments. In a suite of hydrodynamical wind tunnel simulations, we model an intermediate-mass dwarf satellite galaxy (M_* = 10^7.2M_⊙) with a multiphase interstellar medium (ISM; M_ISM = 10^7.9M_⊙) and CGM (M_CGM,vir = 10^8.5M_⊙) along two first-infall orbits to more than 500 Myr past pericenter of a Milky Way–like host. The spatial resolution is ∼79 pc in the star-forming ISM and 316−632 pc in the CGM. Our simulations show that the dwarf satellite CGM removal is fast and effective: more than 95% of the CGM mass is ram pressure stripped within a few hundred megayears, even under a weak ram pressure orbit where the ISM stripping is negligible. The conditions for CGM survival are consistent with the analytical halo gas stripping predictions in McCarthy et al. We also find that including the satellite CGM does not effectively shield its galaxy, and therefore the ISM stripping rate is unaffected. Our results imply that a dwarf galaxy CGM is unlikely to be detected in satellite galaxies; and that the star formation of gaseous dwarf satellites is likely devoid of replenishment from a CGM.

We propose an analytic dual-cone accretion model for horizon-scale images of the cores of low-luminosity active galactic nuclei, including those observed by the Event Horizon Telescope (EHT). Our model is of synchrotron emission from an axisymmetric, magnetized plasma, constrained to flow within two oppositely oriented cones that are aligned with the black hole’s spin axis. We show this model can accurately reproduce images of a variety of time-averaged general relativistic magnetohydrodynamic simulations and that it accurately recovers the black hole spin, orientation, emission scale height, peak emission radius, and fluid flow direction from these simulations within a Bayesian inference framework using radio interferometric data. We show that nontrivial topologies in the images of relativistic accretion flows around black holes can result in nontrivial multimodal solutions when applied to observations with a sparse array, such as the EHT 2017 observations of M87*. The presence of these degeneracies underscores the importance of employing Bayesian techniques to adequately sample the posterior space for the interpretation of EHT measurements. We fit our model to the EHT observations of M87* and find a 95% highest posterior density interval for the mass-to-distance ratio of θ_g ∈ (2.84, 3.75) μas, and give an inclination of θ_o ∈ (11°, 24°). These new measurements are consistent with mass measurements from the EHT and stellar dynamical estimates and with the spin axis inclination inferred from properties of the M87* jet.

Millisecond pulsars (MSPs) are abundant in globular clusters (GCs), which offer favorable environments for their creation. While the advent of recent, powerful facilities led to a rapid increase in MSP discoveries in GCs through pulsation searches, detection biases persist. In this work, we investigate the ability of current and future detections in GCs to constrain the parameters of the MSP population in GCs through a careful study of their luminosity function. Parameters of interest are the number of MSPs hosted by a GC, as well as the mean and the width of their luminosity function, which are typically affected by large uncertainties. While, as we show, likelihood-based studies can lead to ill-behaved posteriors on the size of the MSP population, we introduce a novel, likelihood-free analysis, based on marginal neural ratio estimation, which consistently produces well-behaved posteriors. We focus on the GC Terzan 5 (or Ter 5), which currently counts 48 detected MSPs. We find that 158 MSPs should be hosted in this GC, but the uncertainty on this number remains large. We explore the performance of our new method on simulated Terzan 5-like data sets mimicking possible future observational outcomes. We find that significant improvement on the posteriors can be obtained by adding a reliable measurement of the diffuse radio emission of the GC to the analysis or by improving the detection threshold of current radio pulsation surveys by at least a factor of 2.

We report the spectroscopic confirmation of a massive (), Hubble Space Telescope–dark (m_F150W − m_F444W = 3.6) quiescent galaxy at z_spec = 3.97 in the UNCOVER survey. NIRSpec/PRISM spectroscopy and a nondetection in deep Atacama Large Millimeter/submillimeter Array imaging surprisingly reveals that the galaxy is consistent with a low (<10 M_⊙ yr⁻¹) star formation rate (SFR) despite evidence for moderate dust attenuation. The F444W image is well modeled with a two-component Sérsic fit that favors a compact, r_e ∼ 200 pc, n ∼ 2.9 component and a more extended, r_e ∼ 1.6 kpc, n ∼ 1.7 component. The galaxy exhibits strong color gradients: the inner regions are significantly redder than the outskirts. Spectral energy distribution models that reproduce both the red colors and low SFR in the center of UNCOVER 18407 require both significant (A_v ∼ 1.4 mag) dust attenuation and a stellar mass-weighted age of 900 Myr, implying 50% of the stars in the core already formed by z = 7.5. Using spatially resolved annular mass-to-light measurements enabled by the galaxy’s moderate magnification () to reconstruct a radial mass profile from the best-fitting two-component Sérsic model, we infer a total mass-weighted kpc and log . The early formation of a dense, low SFR, and dusty core embedded in a less attenuated stellar envelope suggests an evolutionary link between the earliest-forming massive galaxies and their elliptical descendants. Furthermore, the disparity between the global, integrated dust properties and the spatially resolved gradients highlights the importance of accounting for radially varying stellar populations when characterizing the early growth of galaxy structure.

The atmospheric response for MeV γ rays (∼0.1–10 MeV) can be characterized in terms of two observed components. The first component is due to photons that reach the detector without scattering. The second component is due to photons that reach the detector after scattering one or more times. While the former can be determined in a straightforward manner, the latter is much more complex to quantify, as it requires tracking the transport of all source photons that are incident on Earth’s atmosphere. The scattered component can cause a significant energy-dependent distortion in the measured spectrum, which is important to account for when making balloon-borne observations. In this work, we simulate the full response for γ-ray transport in the atmosphere. We find that the scattered component becomes increasingly more significant toward lower energies, and at 0.1 MeV, it may increase the measured flux by as much as a factor of ∼2–4, depending on the photon index and off-axis angle of the source. This is particularly important for diffuse sources, whereas the effect from scattering can be significantly reduced for point sources observed with an imaging telescope.

Low-luminosity active galactic nuclei (AGNs) with low-mass black holes (BHs) in the early universe are fundamental to understanding the BH growth and their coevolution with the host galaxies. Utilizing JWST NIRCam Wide Field Slitless Spectroscopy, we perform a systematic search for broad-line Hα emitters (BHAEs) at z ≈ 4–5 in 25 fields of the A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) project, covering a total area of 275 arcmin². We identify 16 BHAEs with FWHM of the broad components spanning from ∼1000 to 3000 km s⁻¹. Assuming that the broad line widths arise as a result of Doppler broadening around BHs, the implied BH masses range from 10⁷ to 10⁸M_⊙, with broad Hα-converted bolometric luminosities of 10^44.5–10^45.5 erg s⁻¹ and Eddington ratios of 0.07–0.47. The spatially extended structure of the F200W stacked image may trace the stellar light from the host galaxies. The Hα luminosity function indicates an increasing AGN fraction toward the higher Hα luminosities. We find possible evidence for clustering of BHAEs: two sources are at the same redshift with a projected separation of 519 kpc; one BHAE appears as a composite system residing in an overdense region with three close companion Hα emitters. Three BHAEs exhibit blueshifted absorption troughs indicative of the presence of high column density gas. We find that the broad-line-selected and photometrically selected BHAE samples exhibit different distributions in the optical continuum slopes, which can be attributed to their different selection methods. The ASPIRE broad-line Hα sample provides a good database for future studies of faint AGN populations at high redshift.