Table of contents

Carbon is an essential element for life on Earth, and the relative abundances of major carbon species (CO₂, CO, and CH₄) in the atmosphere exert fundamental controls on planetary climate and biogeochemistry. Here we employed a theoretical model of atmospheric chemistry to investigate diversity in the atmospheric abundances of CO₂, CO, and CH₄ on Earth-like lifeless planets orbiting Sun-like (F-, G-, and K-type) stars. We focused on the conditions for the formation of a CO-rich atmosphere, which would be favorable for the origin of life. Results demonstrated that elevated atmospheric CO₂ levels trigger photochemical instability of the CO budget in the atmosphere (i.e., CO runaway) owing to enhanced CO₂ photolysis relative to H₂O photolysis. Higher volcanic outgassing fluxes of reduced C (CO and CH₄) also tend to initiate CO runaway. Our systematic examinations revealed that anoxic atmospheres of Earth-like lifeless planets could be classified in the phase space of CH₄/CO₂ versus CO/CO₂, where a distinct gap in atmospheric carbon chemistry is expected to be observed. Our findings indicate that the gap structure is a general feature of Earth-like lifeless planets with reducing atmospheres orbiting Sun-like (F-, G-, and K-type) stars.

The PASSAGES (Planck All-Sky Survey to Analyze Gravitationally-lensed Extreme Starbursts) collaboration has recently defined a sample of 30 gravitationally lensed dusty star-forming galaxies (DSFGs). These rare, submillimeter-selected objects enable high-resolution views of the most extreme sites of star formation in galaxies at cosmic noon. Here, we present the first major compilation of strong lensing analyses using lenstool for PASSAGES, including 15 objects spanning z = 1.1–3.3, using complementary information from 06-resolution 1.1 mm Atacama Large Millimeter/submillimeter Array and 04 5 cm Jansky Very Large Array continuum imaging, in tandem with 1.6 μm Hubble and optical imaging with Gemini-S. Magnifications range from μ = 2 to 28 (median μ = 7), yielding intrinsic infrared luminosities of L_IR = 0.2–5.9 × 10¹³L_⊙ (median 1.4 × 10¹³L_⊙) and inferred star formation rates of 170–6300 M_⊙ yr⁻¹ (median 1500 M_⊙ yr⁻¹). These results suggest that the PASSAGES objects comprise some of the most extreme known starbursts, rivaling the luminosities of even the brightest unlensed objects, further amplified by lensing. The intrinsic sizes of far-infrared continuum regions are large (R_e = 1.7–4.3 kpc; median 3.0 kpc) but consistent with L_IR–R_e scaling relations for z > 1 DSFGs, suggesting a widespread spatial distribution of star formation. With modestly high angular resolution, we explore if these objects might be maximal starbursts. Instead of approaching Eddington-limited surface densities, above which radiation pressure will disrupt further star formation, they are safely sub-Eddington—at least on global, galaxy-integrated scales.

In this study, we report on small-scale magnetic flux ropes (SMFRs) observed as a compact series in a narrow Carrington longitudinal range during three Parker Solar Probe (PSP) encounters. First, during ∼1.5 days of PSP's inbound part of Encounter 4, we identified a series of 11 SMFRs within 14 in longitude over the radial distance of ∼8.4 R_☉ (from ∼44 to 35 R_☉). The identified SMFRs lasted from ∼0.5 to 1.8 hr, and adjacent events were separated mostly by a few hours and up to ∼6.5 hr at the longest, but some events were very closely spaced with intervals of a few ∼tens of minutes or less apart. Most of the identified SMFRs are successfully fitted to the force-free model. The SMFRs are clearly distinguished from the surroundings by a notable reduction in plasma β, which itself was comparably low (less than unity) in the background plasma. Furthermore, the magnetic field and plasma flow within the SMFRs fluctuated significantly less than the more turbulent surroundings. The fluctuations in the surrounding medium exhibited occasional Br polarity reversal (possibly switchbacks) and were Alfvénic to a large extent with far weaker compressional components. The majority of these key features with some differences have also been found in the series of SMFRs and their surroundings identified within 13 or less in longitude during Encounters 1 and 5. We speculate that these SMFRs were repetitively generated by successive reconnection within a very narrow angular zone located close to the Sun but not necessarily at the same radial position.

We present 500 and 700 au resolution 1 and 3 mm Atacama Large Millimeter/submillimeter Array observations, respectively, of protostellar cores in protoclusters Sagittarius B2 (Sgr B2) North (N) and Main (M), parts of the most actively star-forming cloud in our Galaxy. Previous lower-resolution (5000 au) 3 mm observations of this region detected ∼150 sources inferred to be young stellar objects (YSOs) with M > 8 M_⊙. With a 10-fold increase in resolution, we detect 371 sources at 3 mm and 218 sources in the smaller field of view at 1 mm. The sources seen at low resolution are observed to fragment into an average of two objects. About one-third of the observed sources fragment. Most of the sources we report are marginally resolved and are at least partially optically thick. We determine that the observed sources are most consistent with Stage 0/I YSOs, i.e., rotationally supported disks with an active protostar and an envelope, that are warmer than those observed in the solar neighborhood. We report source-counting-based inferred stellar mass and the star formation rate of the cloud: 2800 M_⊙ and 0.0038 M_⊙ yr⁻¹ for Sgr B2 N and 6900 M_⊙ and 0.0093 M_⊙ yr⁻¹ for Sgr B2 M, respectively.

The variability of active galactic nuclei (AGNs) is ubiquitous, but has not yet been understood. Measuring the optical variation properties of AGNs, such as the variation timescale and amplitude, and then correlating them with their fundamental physical parameters has long served as a critical way of exploring the origin of AGN variability and the associated physics of the accretion process in AGNs. Obtaining accurate variation properties of AGNs is thus essential. It has been found that the damped random walk process can describe the AGN optical variation well, but there is a controversy over how long a minimal monitoring baseline is required to obtain unbiased variation properties. In this work, we settle the controversy by exhaustively scrutinizing the complex combination of assumed priors, adopted best-fit values, ensemble averaging methods, and fitting methods. The new proposal is then an optimized solution where unbiased variation properties of an AGN sample possessing the same variation timescale can be obtained with a minimal baseline of about 10 times their variation timescale. Finally, the new optimized solution is used to demonstrate the positive role of the time-domain surveys to be conducted by the Wide Field Survey Telescope in improving constraints on AGN variation properties.

The transport of energetic charged particles through magnetized plasmas is ubiquitous in interplanetary space and astrophysics, and the important physical quantities are the parallel and perpendicular diffusion coefficients of energetic charged particles. In this paper, the influence of solar wind on particle transport is investigated. Using the focusing equation, we obtain parallel and perpendicular diffusion coefficients, accounting for the solar wind effect. For different conditions, the relative importance of the solar wind effect to diffusion is investigated. It is shown that, when energetic charged particles are close to the Sun, for parallel diffusion, the solar wind effect needs to be taken into account. These results are important for studying energetic charged particle transport processes in the vicinity of the Sun.

Magnetic fields are now widely recognized as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Using imposed magnetic fields cannot reliably model the interstellar medium's (ISM) dynamical structure nor phase interactions. Dynamos must be modeled. ISM models exist of turbulent magnetic fields using small-scale dynamo (SSD). Others model the large-scale dynamo (LSD) organizing magnetic fields at the scale of the disk or spiral arms. Separately, neither can fully describe the galactic magnetic field dynamics nor topology. We model the LSD and SSD together at a sufficient resolution to use the low explicit Lagrangian resistivity required. The galactic SSD saturates within 20 Myr. We show that the SSD is quite insensitive to the presence of an LSD and is even stronger in the presence of a large-scale shear flow. The LSD grows more slowly in the presence of SSD, saturating after 5 Gyr versus 1–2 Gyr in studies where the SSD is weak or absent. The LSD primarily grows in warm gas in the galactic midplane. Saturation of the LSD occurs due to α-quenching near the midplane as the growing mean-field produces a magnetic α that opposes the kinetic α. The magnetic energy in our models of the LSD shows a slightly sublinear response to increasing resolution, indicating that we are converging toward the physical solution at 1 pc resolution. Clustering supernovae in OB associations increases the growth rates for both the SSD and the LSD, compared to a horizontally uniform supernova distribution.

Primordial black holes (PBHs) are known as one of the potential candidates for dark matter. They are expected to have formed due to the direct gravitational collapse of density fluctuations in the early Universe. In this regard, examining the merger rate of PBHs within modified theories of gravity can offer a deeper insight into their abundance. In this work, we delve into the calculation of the merger rate of PBHs within the theoretical framework of f(R) gravity. Our analysis reveals an enhancement in the merger rate of PBHs compared to that obtained from general relativity. Additionally, modulating the field strength f_R0 induces shifts in the PBH merger rate, presenting a potential observational signature of modified gravity. We also explore the upper bounds on the abundance of PBHs obtained from f(R) gravity models by comparing the results with gravitational-wave and observational data. The results indicate that in certain regions not excluded by benchmarking data, the parameter space for these upper bounds may be considered reliable.

The electromagnetic emission from the nonrelativistic ejecta launched in neutron star mergers (either dynamically or through a disk wind) has the potential to probe both the total mass and composition of this ejecta. These observations are crucial in understanding the role of these mergers in the production of r-process elements in the Universe. However, many properties of the ejecta can alter the light curves and we must both identify which properties play a role in shaping this emission and understand the effects these properties have on the emission before we can use observations to place strong constraints on the amount of r-process elements produced in the merger. This paper focuses on understanding the effect of the velocity distribution (amount of mass moving at different velocities) for lanthanide-rich ejecta on the light curves and spectra. The simulations use distributions guided by recent calculations of disk outflows and compare the velocity-distribution effects to those of ejecta mass, velocity, and composition. Our comparisons show that uncertainties in the velocity distribution can lead to a factor of 2–4 uncertainties in the inferred ejecta mass based on peak infrared luminosities. We also show that early-time UV or optical observations may be able to constrain the velocity distribution, reducing the uncertainty in the ejecta mass.

In recent years, the Canadian Hydrogen Intensity Mapping Experiment (CHIME) interferometer has revealed a large number of fast radio bursts (FRBs), including a sizable population that demonstrates repeating behavior. This transit facility, employing a real-time FRB search pipeline, continually scans the sky with declinations between −10° and 90° for events with fluences ⪆0.4 Jy ms. We simulate a population of repeating FRBs by performing Monte Carlo simulations of underlying source populations processed through a mock CHIME/FRB observing pipeline. Assuming intrinsic repeater rates follow a Poisson distribution, we test assumptions about the burst populations of the repeater sample, and construct models of the FRB sample assuming various cosmological distributions. We infer the completeness of CHIME/FRB observations as a function of observing cadence and redshifts out to 0.5. We find that, if all simulated bursts have a fixed Poisson probability of repetition over their integrated time of observation, repeating burst detections across comoving volume should continue to grow near linearly on the order of decades. We predict that around 170 of the current CHIME/FRB one-off sources will ultimately repeat. We also make projections for FRB repeaters by future facilities and demonstrate that the number of repeaters they find could saturate on a ∼3 yr timescale.

We present a magnetic configuration of a compound solar eruption observed on 2012 March 10, from NOAA AR 11429 near the disk center, which displayed a soft X-ray sigmoid before the eruption. We constructed a series of magnetic field models, including double-decker flux rope configurations, using the flux rope insertion method. This produces three-dimensional coronal magnetic field models constrained by the photospheric magnetogram and observed EUV coronal structures. We used different combinations of flux rope paths. We found that two flux ropes sharing the same path at different heights quickly experience a partial merging in the initial iteration of the magnetofrictional relaxation process. Different paths with less than 30% overlap allowed us to construct stable double-decker structures. The high spatial and temporal resolution of the Solar Dynamics Observatory/Atmospheric Imaging Assembly facilitated the selection of a best-fit model that matches the observations best. Moreover, by varying fluxes in this validated nonlinear force-free field double-decker configuration, we successfully reproduce all three scenarios of eruptions of double-decker configurations: (i) eruption due to the instability of higher flux rope; (ii) eruption due to rising lower flux rope and merging with higher flux rope; and (iii) eruption due to the instability of both flux ropes. This demonstrates that magnetofrictional simulation can capture the large-scale magnetic structure of eruptions for a realistic field configuration at eruption onset.

We report on the X-ray spectral and spatial evolution of the symbiotic star R Aqr. Through a multiepoch observational campaign performed with Chandra between 2017 and 2022, we study the X-ray emission of this binary system, composed of an evolved red giant star and a white dwarf (WD). This analysis is particularly timely as the WD approached the periastron in late 2018/early 2019; thus, mass transfer, jet emission, and outburst phenomena are to be expected. Through detailed spectral analysis, we detect a significant rise in the soft X-ray (0.5–2 keV) emission of R Aqr, likely linked to jet emission, followed by a decay toward the previous quiescent state. The hard X-ray emission (5–8 keV) is not immediately affected by the periastron passage; the hard component, after maintaining the same flux level between 2017 and 2021, rapidly decays after 2022. Possible explanations for this are a change in the reflection properties of the medium surrounding the binary, obscuration of the central region by material ejected during the periastron passage, or even the partial/complete destruction of the inner regions of the accretion disk surrounding the WD. In addition to this activity in the central region, extended emission is also detected, likely linked to a hot spot in a pre-outburst-emitted jet, which can be observed moving away from the system's central region.

The Galactic global magnetic field is thought to play a vital role in shaping Galactic structures such as spiral arms and giant molecular clouds. However, our knowledge of magnetic field structures in the Galactic plane at different distances is limited, as measurements used to map the magnetic field are the integrated effect along the line of sight. In this study, we present the first ever tomographic imaging of magnetic field structures in a Galactic spiral arm. Using optical stellar polarimetry over a $17^{\prime} \times 10^{\prime}$ field of view, we probe the Sagittarius spiral arm. Combining these data with stellar distances from the Gaia mission, we can isolate the contributions of five individual clouds along the line of sight by analyzing the polarimetry data as a function of distance. The observed clouds include a foreground cloud (d < 200 pc) and four clouds in the Sagittarius arm at 1.23, 1.47, 1.63, and 2.23 kpc. The column densities of these clouds range from 0.5 to 2.8 × 10²¹ cm⁻². The magnetic fields associated with each cloud show smooth spatial distributions within their observed regions on scales smaller than 10 pc and display distinct orientations. The position angles projected on the plane of the sky, measured from the Galactic north to the east, for the clouds in increasing order of distance are 135°, 46°, 58°, 150°, and 40°, with uncertainties of a few degrees. Notably, these position angles deviate significantly from the direction parallel to the Galactic plane.

We have targeted the dusty symbiotic Mira system HM Sge with four instruments from the IR to the UV. We have used these observations along with archival observations to study how the system has been evolving after its 1975 nova-like outburst. We have detected rovibrational water emission in a symbiotic system for the first time using new EXES high-spectral-resolution infrared spectroscopy. The features, detected in emission, have velocities consistent with the systemic velocity but do not show any clear evidence of high-velocity outflows. Mid-infrared photometry and grism spectroscopy show that the oxygen-rich asymptotic giant branch dust and dust output have shown little to no change over the past 39 years. In the optical/UV, we detect three main [N ii] nebular features that were detected 22 years ago. Two of these features show a small amount of movement, corresponding to average outflows speeds of 38 and 78 km s⁻¹ since they were previously observed; some previously detected [N ii] features are no longer visible. New UV spectroscopy has shown that the nebular environment continues to steadily relax after the system's 1975 outburst. The data suggest, however, that the temperature of the hot component has increased from 200,000 K in 1989 to greater than 250,000 K now. Our new and archival observations suggest that the evolution of the system after its outburst is swift with little to no major changes after a period of a couple of years.

Recent radio observations with the Low Frequency Array (LOFAR) discovered diffuse emission extending beyond the scale of classical radio halos. The presence of such megahalos indicates that the amplification of the magnetic field and acceleration of relativistic particles are working in the cluster outskirts, presumably due to the combination of shocks and turbulence that dissipate energy in these regions. Cosmological magnetohydrodynamical (MHD) simulations of galaxy clusters suggest that solenoidal turbulence has a significant energy budget in the outskirts of galaxy clusters. In this paper, we explore the possibility that this turbulence contributes to the emission observed in megahalos through second-order Fermi acceleration of relativistic particles and magnetic field amplification by the dynamo. We focus on the case of A2255 and find that this scenario can explain the basic properties of the diffuse emission component that is observed under assumptions that are used in previous literature. More specifically, we conduct a numerical follow-up, solving the Fokker–Planck equation by using a snapshot of an MHD simulation and deducing the synchrotron brightness integrated along the lines of sight. We find that a volume-filling emission, ranging between 30% and almost 100% of the projected area, depending on our assumptions on the particle diffusion and transport, can be detected at LOFAR sensitivities. Assuming a magnetic field B ∼ 0.2 μG, as derived from a dynamo model applied to the emitting region, we find that the observed brightness can be matched when ∼1% of the solenoidal turbulent energy flux is channeled into particle acceleration.

We present the first star formation history (SFH) and age–metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color–magnitude diagram (CMD) reaches below the oldest main sequence turnoff with a signal-to-noise ratio = 10 at M_F090W = + 4.6 mag. This is the deepest CMD for any galaxy that is not a satellite of the Milky Way. We use Hubble Space Telescope (HST) optical imaging that overlaps with the NIRCam observations to directly evaluate the SFHs derived based on data from the two great observatories. The JWST and HST-based SFHs are in excellent agreement. We use the metallicity distribution function measured from stellar spectra to confirm the trends in the AMRs based on the JWST data. Together, these results confirm the efficacy of recovering an SFH and AMR with the NIRCam F090W−F150W filter combination, and validate the sensitivity and accuracy of stellar evolution libraries in the near-infrared relative to the optical for SFH recovery work. From the JWST data, WLM shows an early onset to star formation, followed by an extended pause post-reionization before star formation reignites, which is qualitatively similar to what has been observed in the isolated galaxies Leo A and Aquarius. Quantitatively, 15% of the stellar mass formed in the first Gyr, while only 10% formed over the next ∼5 Gyr. The stellar mass then rapidly doubled in ∼2.5 Gyr, followed by constant star formation over the last ∼5 Gyr.

We present the first test of coasting cosmological models with gravitational-wave (GW) standard sirens observed in the first three observing runs of the LIGO–Virgo–KAGRA detector network. We apply the statistical galaxy catalog method adapted to coasting cosmologies and infer constraints on the H₀ Hubble constant for the three fixed values of the curvature parameter $k=\left\{-1,0,+1\right\}$ in ${H}_{0}^{2}{c}^{-2}$ units. The maximum posteriors and 68.3% highest density intervals we obtained from a combined analysis of 46 dark siren detections and a single bright siren detection are ${H}_{0}=\left\{{68.1}_{-5.6}^{+8.5},{67.5}_{-5.2}^{+8.3},{67.1}_{-5.8}^{+6.6}\right\}\,\mathrm{km}\ {{\rm{s}}}^{-1}\ {\mathrm{Mpc}}^{-1}$, respectively. All our constraints on H₀ are consistent within 1σ with the H₀ measured with the differential age method, which provides a constraint on H₀ in coasting cosmologies independently from k. Our results constrain all cosmological models with a(t) ∝ t linear expansion in the luminosity distance and redshift range of the 47 LIGO–Virgo detections, i.e., d_L ≲ 5Gpc and z ≲ 0.8, which practically include all (both strictly linear and quasi-linear) models in the coasting model family. As we have found, the coasting models and the Lambda cold dark matter (or ΛCDM) model fit equally well to the applied set of GW detections.

We present a spatially resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the Atacama Large Millimeter/submillimeter Array Comprehensive High-resolution Extragalactic Molecular Inventory, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3–7 transitions from nine molecular species to delineate the position–position–velocity distributions of column density (${N}_{{{\rm{H}}}_{2}}$), volume density (${n}_{{{\rm{H}}}_{2}}$), and temperature (T_kin) at 27 pc resolution. Two distinct components are detected: a low-density component with $({n}_{{{\rm{H}}}_{2}},{T}_{\mathrm{kin}})\sim ({10}^{3.3}\ {\mathrm{cm}}^{-3},85\ {\rm{K}})$ and a high-density component with $({n}_{{{\rm{H}}}_{2}},{T}_{\mathrm{kin}})\,\sim ({10}^{4.4}\ {\mathrm{cm}}^{-3},110\ {\rm{K}})$, separated at ${n}_{{{\rm{H}}}_{2}}\sim {10}^{3.8}\ {\mathrm{cm}}^{-3}$. NGC 253 has ∼10 times the high-density gas mass and ∼3 times the dense-gas mass fraction of the GC. These properties are consistent with their HCN/CO ratio but cannot alone explain the factor of ∼30 difference in their star formation efficiencies (SFEs), contradicting the dense-gas mass to star formation rate scaling law. The ${n}_{{{\rm{H}}}_{2}}$ histogram toward NGC 253 exhibits a shallow declining slope up to ${n}_{{{\rm{H}}}_{2}}\sim {10}^{6}\ {\mathrm{cm}}^{-3}$, while that of the GC steeply drops in ${n}_{{{\rm{H}}}_{2}}\gtrsim {10}^{4.5}\ {\mathrm{cm}}^{-3}$ and vanishes at 10⁵ cm⁻³. Their dense-gas mass fraction ratio becomes consistent with their SFEs when the threshold ${n}_{{{\rm{H}}}_{2}}$ for the dense gas is taken at ∼10^4.2−4.6 cm⁻³. The rich abundance of gas above this density range in the NGC 253 CMZ, or its scarcity in the GC, is likely to be the critical difference characterizing the contrasting star formation in the centers of the two galaxies.

The galaxy CID-42 (CXOC J100043.1+020637.2) at z = 0.359 has been proposed to contain a promising candidate for a gravitational-wave recoiling supermassive black hole (SMBH), a slingshot SMBH from a triple-SMBH interaction, or a kiloparsec-scale dual active galactic nucleus (AGN). These claims were primarily based on a pair of bright cores separated by ∼05 resolved in optical Hubble Space Telescope (HST) imaging. Existing HST, Chandra, and ground-based imaging and spectroscopy are unable to confirm either scenario. With improved spatial resolution, depth, and IR wavelength coverage, NIRCam multiband imaging from the COSMOS-Web JWST treasury program well resolved the two cores in CID-42, revealing a significant stellar bulge for both cores (with stellar masses of ∼10¹⁰M_⊙ for both). JWST imaging further revealed that only the SE core contains an unobscured AGN point source, based on both image decomposition and spectral energy distribution fitting. There is no evidence for AGN activity in the NW core based on the current data. These new observations unambiguously rule out the gravitational-wave-recoiling and slingshot-SMBH scenarios and establish CID-42 as a low-redshift merging pair of galaxies, with likely only one active AGN in the system. These results demonstrate the unparalleled capabilities of JWST (even with imaging alone) in studying the galactic-scale environment of merging galaxies and SMBHs.

The (sub)millimeter wavelengths (86–690 GHz) very long baseline interferometry will provide ∼5–40 μas angular resolution, ∼10 mJy baseline sensitivity, and ∼1 μas yr⁻¹ proper-motion precision, which can directly detect supermassive black hole binary (SMBHB) systems by imaging the two visible sources and tracking their relative motions. Such a way exhibits an advantage compared to indirect detect methods of observing periodic signals in motion and light curves, which are difficult to confirm from competing models. Moreover, tracking relative motion at (sub)millimeter wavelengths is more reliable, as there is a negligible offset between the emission region and the black hole center. In this way, it is unnecessary to correct the black hole location by a prior of jet morphology as it would be required at longer wavelengths. We extend the formalism developed in D'Orazio & Loeb (2018) to link the observations with the orbital evolution of SMBHBs from the ≲10 kpc dynamical friction stages to the ≲0.01 pc gravitational radiation stages, and estimate the detectable numbers of SMBHBs. By assuming 5% of active galactic nuclei holding SMBHBs, we find that the number of detectable SMBHBs with redshift z ≤ 0.5 and mass M ≤ 10¹¹M_⊙ is about 20. Such a detection relies heavily on proper-motion precision and sensitivity. Furthermore, we propose that the simultaneous multifrequency technique plays a key role in meeting the observational requirements.

We continue the numerical modeling of a corotating interaction region (CIR) and the effects it has on solar-rotational recurrent variations of galactic cosmic rays (GCRs). A magnetohydrodynamic model is adapted to simulate the background solar wind plasma with a CIR structure in the inner heliosphere, which is incorporated into a comprehensive Parker-type transport model. The focus is on the simulation of the effects of a CIR on GCR protons and the two helium isotopes as a function of heliolongitude. This is to establish whether the difference in composition affects how they are modulated by the CIR in terms of their distribution in longitude. It is demonstrated that particle diffusion and drift influence the effects of the CIR with increasing rigidity from 100 MV up to 15 GV. It is found that protons and helium isotopes are modulated differently with longitude by the CIR and that particle drift influences the modulation effects in longitude. These differences dissipate with increasing rigidity. The final results are focused on the simulated amplitude of these GCR flux variations as a function of rigidity. The amplitude displays a power-law behavior above ∼1 GV with an index similar to the power index of the rigidity dependence of the assumed diffusion coefficients. The simulations further show that below this rigidity, the amplitude at first flattens off, displaying a plateau-like profile, but it then increases systematically with decreasing rigidity below ∼0.3 GV. Again, a power-law behavior is displayed, but it is completely different from that above 1 GV.

Tidal heating on Io due to its finite eccentricity was predicted to drive surface volcanic activity, which was subsequently confirmed by the Voyager spacecraft. Although the volcanic activity in Io is more complex, in theory volcanism can be driven by runaway melting in which the tidal heating increases as the mantle thickness decreases. We show that this runaway melting mechanism is generic for a composite planetary body with liquid core and solid mantle, provided that (i) the mantle rigidity, μ, is comparable to the central pressure, i.e., μ/(ρgR_P) ≳ 0.1 for a body with density ρ, surface gravitational acceleration g, and radius R_P; (ii) the surface is not molten; (iii) tides deposit sufficient energy; and (iv) the planet has nonzero eccentricity. We calculate the approximate liquid core radius as a function of μ/(ρgR_P), and find that more than 90% of the core will melt due to this runaway for μ/(ρgR_P) ≳ 1. From all currently confirmed exoplanets, we find that the terrestrial planets in the L 98-59 system are the most promising candidates for sustaining active volcanism. However, uncertainties regarding the quality factors and the details of tidal heating and cooling mechanisms prohibit definitive claims of volcanism on any of these planets. We generate synthetic transmission spectra of these planets assuming Venus-like atmospheric compositions with an additional 5%, 50%, and 98% SO₂ component, which is a tracer of volcanic activity. We find a ≳3σ preference for a model with SO₂ with 5–10 transits with JWST for L 98-59bcd.

Young solar-type stars frequently produce superflares, serving as a unique window into the young Sun-Earth environments. Large solar flares are closely linked to coronal mass ejections (CMEs) associated with filament/prominence eruptions, but observational evidence for stellar superflares remains scarce. Here, we present a 12-day, multiwavelength campaign observation of young solar-type star EK Draconis (G1.5V, 50–120 Myr age) utilizing the Transiting Exoplanet Survey Satellite, the Neutron star Interior Composition ExploreR, and the Seimei telescope. The star has previously exhibited blueshifted Hα absorptions as evidence for a filament eruption associated with a superflare. Our simultaneous optical and X-ray observations identified three superflares of 1.5 × 10³³–1.2 × 10³⁴ erg. We report the first discovery of two prominence eruptions on a solar-type star, observed as blueshifted Hα emissions at speeds of 690 and 430 km s⁻¹ and masses of 1.1 × 10¹⁹ and 3.2 × 10¹⁷ g, respectively. The faster, massive event shows a candidate of post-flare X-ray dimming with the amplitude of up to ∼10%. Several observational aspects consistently point to the occurrence of a fast CME associated with this event. The comparative analysis of the estimated length scales of flare loops, prominences, possible dimming region, and starspots provides the overall picture of the eruptive phenomena. Furthermore, the energy partition of the observed superflares in the optical and X-ray bands is consistent with flares from the Sun, M-dwarfs, and close binaries, yielding the unified empirical relations. These discoveries provide profound implications of the impact of these eruptive events on early Venus, Earth, and Mars and young exoplanets.

To better understand the mixing and mass loss experienced by low-mass stars as they ascend the asymptotic giant branch (AGB), I have gathered from the literature the abundances of CNO and s-process elements in post-AGB stars in Galactic globular clusters. These species are mixed to the surface during third dredge-up (3DU) events, so their abundance should increase as the star ascends the AGB. Of the 17 stars in this sample, CNO abundances are available for 11. Of these, four are enhanced in CNO relative to the red giant branch stars from which they descended, which I take as evidence of 3DU on the AGB. The enhancement is mainly in the form of carbon. Of the six stars for which only heavy-element abundances are available, one shows s-process enhancements that previous authors have interpreted as evidence of 3DU. Combining these 17 stars with other recent samples reveals that most globular-cluster post-AGB stars have luminosities $\mathrm{log}(L/{L}_{\odot })\sim 3.25$. They are the progeny of blue horizontal-branch (HB) stars in clusters with intermediate metallicity ([Fe/H] ∼ −1.5). A second group consists of sub-luminous stars associated with high-metallicity clusters ([Fe/H] ∼ −1.0) with red HBs. They may be burning helium, rather than hydrogen. A third group of hot, super-luminous stars is evolving quickly across the Hertzsprung–Russell diagram. Some of them may be merger remnants.

To investigate how magnetic reconnection (MR) accelerates electrons to a power-law energy spectrum in solar flares, we explore the scaling of a kinetic model proposed by Che & Zank (CZ) and compare it to observations. Focusing on thin current sheet MR particle-in-cell (PIC) simulations, we analyze the impact of domain size on the evolution of the electron Kelvin–Helmholtz instability (EKHI). We find that the duration of the growth stage of the EKHI (${t}_{G}\sim {{\rm{\Omega }}}_{e}^{-1}$) is short and remains nearly unchanged because the electron gyrofrequency Ω_e is independent of domain size. The quasi-steady stage of the EKHI (t_MR) dominates the electron acceleration process and scales linearly with the size of the simulations as L/v_A0, where v_A0 is the Alfvén speed. We use the analytical results obtained by CZ to calculate the continuous temporal evolution of the electron energy spectra from PIC simulations and linearly scale them to solar flare observational scales. For the first time, an electron acceleration model predicts the sharp two-stage transition observed in typical soft–hard–harder electron energy spectra, implying that the electron acceleration model must be efficient with an acceleration timescale that is a small fraction of the duration of solar flares. Our results suggest that we can use PIC MR simulations to investigate the observational electron energy spectral evolution of solar flares if the ratio t_MR/t_G is sufficiently small, i.e., ≲10%.

In the latest data release from the Fermi Gamma-ray Space Telescope (the 4th Fermi LAT 14 yr Catalog, or 4FGL), more than 50% of the Galactic sources are yet to be identified. We observed 13 unidentified Fermi LAT sources with the Chandra X-Ray Observatory to explore their nature. We report the results of the classification of X-ray sources in the fields of these γ-ray sources and discuss the implications for their nature. We use multiwavelength (MW) data for a machine-learning classification, accompanied by a more detailed spectral/variability analysis for brighter sources. Eight 4FGL sources have γ-ray pulsars within their position error ellipses. We consider three of these pulsars (PSR J1906+0722, PSR J1105–6037, and PSR J1358–6025) to be detected in X-rays, while PSR J1203–6242 shows a hint of X-ray emission. Within the positional uncertainties of three of the 4FGL sources, we detect X-ray sources that may be yet unknown pulsars, depending on the MW association. In addition to point sources, we discovered two extended sources, one of which is likely to be a bow-shock pulsar-wind nebula associated with PSR J1358–6025. Finally, we classify other X-ray sources detected in these observations and report the most interesting classifications.

The connection between galaxies and dark matter halos is often quantified using the stellar mass–halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at z ≳ 2 can also constrain the SMHM relation via the galaxy autocorrelation function and through the cross-correlation between galaxies and Lyα absorption measured in transverse sight lines; however, such studies are very few and have produced some unexpected or inconclusive results. We use ∼3000 spectra of z ∼ 2.5 galaxies from the Lyα Tomography IMACS Survey (LATIS) to measure the galaxy–galaxy and galaxy–Lyα correlation functions in four bins of stellar mass spanning 10^9.2 ≲ M_*/M_⊙ ≲ 10^10.5. Parallel analyses of the MultiDark N-body and ASTRID hydrodynamic cosmological simulations allow us to model the correlation functions, estimate covariance matrices, and infer halo masses. We find that results of the two methods are mutually consistent and broadly accord with standard SMHM relations. This consistency demonstrates that we are able to measure and model Lyα transmission fluctuations δ_F in LATIS accurately. We also show that the galaxy–Lyα cross-correlation, a free by-product of optical spectroscopic galaxy surveys at these redshifts, can constrain halo masses with similar precision to galaxy–galaxy clustering.

While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial-resolution observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A.

Filamentary structures in neutral hydrogen (H i) emission are well aligned with the interstellar magnetic field, so H i emission morphology can be used to construct templates that strongly correlate with measurements of polarized thermal dust emission. We explore how the quantification of filament morphology affects this correlation. We introduce a new implementation of the Rolling Hough Transform (RHT) using spherical harmonic convolutions, which enables efficient quantification of filamentary structure on the sphere. We use this Spherical RHT algorithm along with a Hessian-based method to construct H i-based polarization templates. We discuss improvements to each algorithm relative to similar implementations in the literature and compare their outputs. By exploring the parameter space of filament morphologies with the Spherical RHT, we find that the most informative H i structures for modeling the magnetic field structure are the thinnest resolved filaments. For this reason, we find a ∼10% enhancement in the B-mode correlation with polarized dust emission with higher-resolution H i observations. We demonstrate that certain interstellar morphologies can produce parity-violating signatures, i.e., nonzero TB and EB, even under the assumption that filaments are locally aligned with the magnetic field. Finally, we demonstrate that B modes from interstellar dust filaments are mostly affected by the topology of the filaments with respect to one another and their relative polarized intensities, whereas E modes are mostly sensitive to the shapes of individual filaments.

Current endeavours in exoplanet characterization rely on atmospheric retrieval to quantify crucial physical properties of remote exoplanets from observations. However, the scalability and efficiency of said technique are under strain with increasing spectroscopic resolution and forward model complexity. The situation has become more acute with the recent launch of the James Webb Space Telescope and other upcoming missions. Recent advances in machine learning provide optimization-based variational inference as an alternative approach to perform approximate Bayesian posterior inference. In this investigation we developed a normalizing-flow-based neural network, combined with our newly developed differentiable forward model, Diff-τ, to perform Bayesian inference in the context of atmospheric retrievals. Using examples from real and simulated spectroscopic data, we demonstrate the advantages of our proposed framework: (1) training our neural network does not require a large precomputed training set and can be trained with only a single observation; (2) it produces high-fidelity posterior distributions in excellent agreement with sampling-based retrievals; (3) it requires up to 75% fewer forward model calls to converge to the same result; and (4) this approach allows formal Bayesian model selection. We discuss the computational efficiencies of Diff-τ in relation to TauREx3's nominal forward model and provide a "lessons learned" account of developing radiative transfer models in differentiable languages. Our proposed framework contributes toward the latest development of neural network–powered atmospheric retrieval. Its flexibility and significant reduction in forward model calls required for convergence holds the potential to be an important addition to the retrieval tool box for large and complex data sets along with sampling-based approaches.

We performed synthetic observations of the Ulrich, Cassen, and Moosman (UCM) model to understand the relation between the physical structures of the infalling envelope around a protostar and their observational features in molecular lines, adopting L1527 as an example. We also compared the physical structure and synthetic position–velocity (P–V) diagrams of the UCM model and a simple ballistic (SB) model. There are multiple ways to compare synthetic data with observational data. We first calculated the correlation coefficient. The UCM model and the SB model show similarly good correlation with the observational data. While the correlation reflects the overall similarity between the cube datasets, we can alternatively compare specific local features, such as the centrifugal barrier in the SB model or the centrifugal radius in the UCM model. We evaluated systematic uncertainties in these methods. In the case of L1527, the stellar mass values estimated using these methods are all lower than the value derived from previous Keplerian analysis of the disk. This may indicate that the gas infall motion in the envelope is retarded by, e.g., magnetic fields. We also showed analytically that, in the UCM model, the spin-up feature of the P–V diagram is due to the infall velocity rather than the rotation. The line-of-sight velocity V is thus ∝x^−0.5, where x is the offset. If the infall is retarded, rotational velocity should dominate so that V is proportional to x⁻¹, as is often observed in the protostellar envelope.

Among Type Ia supernova remnants (SNRs), Tycho's SNR has been considered as a typical object from the viewpoints of its spectroscopic, morphological, and environmental properties. A recent reanalysis of Chandra data showed that its forward shock is experiencing a substantial deceleration since around 2007, which suggests recent shock interactions with a dense medium as a consequence of a cavity-wall environment inside a molecular cloud. Such a nonuniform environment can be linked back to the nature and activities of its progenitor. In this study, we perform hydrodynamic simulations to characterize Tycho's cavity-wall environment using the latest multiepoch proper motion measurements of the forward shock. A range of parameters for the environment is explored in the hydrodynamic models to fit with the observational data for each azimuthal region. Our results show that a wind-like cavity with ρ(r) ∝ r⁻² reconciles with the latest data better than a uniform medium with a constant density. In addition, our best-fit model favors an anisotropic wind with an azimuthally varying wind parameter. The overall result indicates a mass-loss rate which is unusually high for the conventional single-degenerate explosion scenario.

We present a semi-analytic model for the growth, drift, desorption, and fragmentation of millimeter- to meter-sized particles in protoplanetary disks. Fragmentation occurs where particle collision velocities exceed critical fragmentation velocities. Using this criterion, we produce fragmentation regions in disk orbital radius–particle size phase space for particles with a range of material properties, structures, and compositions (including SiO₂, Mg₂SiO₄, H₂O, CO₂, and CO). For reasonable disk conditions, compact aggregate H₂O, CO₂, and CO ice particles do not reach destructive relative velocities and are thus not likely to undergo collisional fragmentation. Uncoated silicate particles are more susceptible to collisional destruction and are expected to fragment in the inner disk, consistent with previous work. We then calculate the growth, drift, and sublimation of small particles, initially located in the outer disk. We find that ice-coated particles can avoid fragmentation as they grow and drift inward under a substantial range of disk conditions, as long as the particles are aggregates composed of 0.1 μm-sized monomers. Such particles may undergo runaway growth in disk regions abundant in H₂O or CO₂ ice, depending on the assumed disk temperature structure. These results indicate that icy collisional growth to planetesimally relevant sizes may happen efficiently throughout a disk's lifetime, and is particularly robust at early times when the disk's dust-to-gas ratio is comparable to that of the interstellar medium.

Low-frequency gravitational-wave experiments such as the Laser Interferometer Space Antenna and pulsar timing arrays are expected to detect individual massive black hole (MBH) binaries and mergers. However, secure methods of identifying the exact host galaxy of each MBH merger among the large number of galaxies in the gravitational-wave localization region are currently lacking. We investigate the distinct morphological signatures of MBH merger host galaxies, using the Romulus25 cosmological simulation. We produce mock telescope images of 201 simulated galaxies in Romulus25 hosting recent MBH mergers through stellar population synthesis and dust radiative transfer. Based on comparisons to mass- and redshift-matched control samples, we show that combining multiple morphological statistics via a linear discriminant analysis enables identification of the host galaxies of MBH mergers, with accuracies that increase with chirp mass and mass ratio. For mergers with high chirp masses (≳10^8.2M_⊙) and high mass ratios (≳0.5), the accuracy of this approach reaches ≳80%, and does not decline for at least ∼1 Gyr after numerical merger. We argue that these trends arise because the most distinctive morphological characteristics of MBH merger and binary host galaxies are prominent classical bulges, rather than relatively short-lived morphological disturbances from their preceding galaxy mergers. Since these bulges are formed though major mergers of massive galaxies, they lead to (and become permanent signposts for) MBH binaries and mergers that have high chirp masses and mass ratios. Our results suggest that galaxy morphology can aid in identifying the host galaxies of future MBH binaries and mergers.

Compositional convection is atmospheric mixing driven by density variations caused by compositional gradients. Previous studies have suggested that compositional gradients of atmospheric trace species within planetary atmospheres can impact convection and the final atmospheric temperature profile. In this work, we employ 3D convection-resolving simulations using Cloud Model 1 (CM1) to gain a fundamental understanding of how compositional variation influences convection and the final atmospheric state of exoplanet atmospheres. We perform 3D initial value problem simulations of noncondensing compositional convection for Earth-air, H₂, and CO₂ atmospheres. Conventionally, atmospheric convection is assumed to mix the atmosphere to a final, marginally stable state defined by a unique temperature profile. However, when there is compositional variation within an atmosphere, a continuous family of stable end states is possible, differing in the final state composition profile. Our CM1 simulations are used to determine which of the family of possible compositional end states is selected. Leveraging the results from our CM1 simulations, we develop a dry convective adjustment scheme for use in general circulation models (GCMs). This scheme relies on an energy analysis to determine the final adjusted atmospheric state. Our convection scheme produces results that agree with our CM1 simulations and can easily be implemented in GCMs to improve modeling of compositional convection in exoplanet atmospheres.

The current sample of 12 radio-quiet isolated neutron stars that emit strongly in X-rays (XINSs) is both small and heterogeneous, limiting its usefulness for understanding the physics of neutron star atmospheres and cooling rates and for constraining the equation of state of neutron degenerate matter. Utilizing the ROSAT 1RXS and 2RXS data sets, in conjunction with the Sloan Digital Sky Survey Data Release 17 and other companion multiwavelength surveys, we have extended previous searches for blank-field X-ray source candidate XINSs, ultimately recovering two known XINSs while identifying 46 new, unstudied candidate fields devoid of likely multiwavelength counterparts. In this publication, we describe our selection approach and provide detailed information regarding our sample of new candidate XINSs. Future opportunities to verify or to refute these X-ray sources as isolated neutron stars by obtaining more accurate X-ray source positions, quality X-ray spectra, or deeper optical imaging are also discussed.

Galaxies that are invisible in deep optical–near-infrared imaging but detected at longer wavelengths have been the focus of several recent observational studies, with speculation that they could constitute a substantial missing population and even dominate the cosmic star formation rate density at z ≳ 4. The depths now achievable with JWST at the longest wavelengths probed by the Hubble Space Telescope (HST), coupled with the transformative resolution at longer wavelengths, are already enabling detailed, spatially resolved characterization of sources that were invisible to HST, often known as "HST-dark" galaxies. However, until now, there has been little theoretical work to compare against. We present the first simulation-based study of this population, using highly resolved galaxies from the Feedback in Realistic Environments project, with multiwavelength images along several lines of sight forward-modeled using radiative transfer. We naturally recover a population of modeled sources that meet commonly used selection criteria (H_AB > 27 mag and H_AB − F444W > 2.3). These simulated HST-dark galaxies lie at high redshifts (z = 4–7), have high levels of dust attenuation (A_V = 2–4), and display compact recent star formation (R_{1/2,4.4 μm} ≲ 1 kpc). Orientation is very important: for all but one of the 17 simulated galaxy snapshots with HST-dark sight lines, there exist other sight lines that do not meet the criteria. This result has important implications for comparisons between observations and models that do not resolve the detailed star-dust geometry, such as semianalytic models or coarsely resolved hydrodynamical simulations. Critically, we demonstrate that HST-dark sources are not an unexpected or exotic population, but a subset of high-redshift, highly dust-attenuated sources viewed along certain lines of sight.

Here we report an optical quasiperiodic oscillation (QPO) with a period of ∼134 days detected in g- and r-band light curves of the narrow-line Seyfert 1 galaxy TXS 1206+549 at a redshift of 1.34 with data from observations at the Zwicky Transient Facility (ZTF). After considering the trial factor, the significance levels in the two bands are 3.1σ and 2.6σ, respectively. The QPO signal presents about 10 cycles ranging from 2018 March to 2021 December, thus lasting ∼4 yr. A nearly sinusoidal profile also appears in the folded light curves by using a phase-resolved analysis. Interestingly, in the simultaneous light curve with the timescale of ZTF observations, a potential periodic signal with a similar period is detected in the o-band light curve from Asteroid Terrestrial-impact Last Alert System data; additionally, a weak peak is also detected at a similar period in the γ-ray light curve obtained from Fermi Gamma-ray Space Telescope data. Some potential origins of periodicities in active galactic nuclei are discussed for the QPO reported here.

We present a systematic study of the environmental impact on star formation activities of galaxies using a mass-complete sample of ∼170k galaxies at z < 4 from the latest COSMOS2020 catalog. At z < 1, we find that the mean star formation rate (SFR) of all galaxies decreases with increasing density of the environment. However, when we only consider star-forming galaxies, the mean SFR becomes independent of the environment at z < 1. At z > 2, we observe a clear positive correlation between the SFR and the density of the environment for all the galaxies. On the other hand, the stellar mass of the galaxies increases significantly with the environment at all redshifts except for star-forming galaxies at z < 1. The fraction of quiescent galaxies increases with increasing density of the environment at z < 2, and the morphology–density relation is confirmed to be present up to z ∼ 1. We also find that environmental quenching is negligible at z > 1, whereas mass quenching is the dominant quenching mechanism for massive galaxies at all redshifts. Based on these results, we argue that stellar mass-regulated physical processes might be the major driving force for star formation activities of galaxies. At low redshift (z < 1) massive galaxies are quenched primarily due to their high mass, resulting in a normal SFR–density relation. At high redshift (z > 2) most of the galaxies are star-forming ones tightly following the star-forming main sequence, and the difference in their stellar mass in different environments naturally leads to a reversal of the SFR–density relation.

Chromospheric differential rotation is a key component in comprehending the atmospheric coupling between the chromosphere and the photosphere at different phases of the solar cycle. In this study, we therefore utilize the newly calibrated multidecadal Ca ii K spectroheliograms (1907–2007) from the Kodaikanal Solar Observatory (KoSO) to investigate the differential rotation of the solar chromosphere using the technique of image cross-correlation. Our analysis yields the chromospheric differential rotation rate Ω(θ) = (14.61 ± 0.04–2.18 ± 0.37 ${\sin }^{2}\theta -1.10\pm 0.61{\sin }^{4}\theta )^\circ$ day⁻¹. These results suggest the chromospheric plages exhibit an equatorial rotation rate 1.59% faster than the photosphere when compared with the differential rotation rate measured using sunspots and also a smaller latitudinal gradient compared to the same. To compare our results to those from other observatories, we have applied our method on a small sample of Ca ii K data from Rome, Meudon, and Mount Wilson observatories, which support our findings from KoSO data. Additionally, we have not found any significant north–south asymmetry or any systematic variation in chromospheric differential rotation over the last century.

Solar wind magnetic field fluctuations exhibit a complex multiscale nature, often encompassing ion-scale discontinuities and MHD-scale Alfvénic fluctuations. Both of these types of structures are thought to play a critical role in plasma heating and turbulence dissipation. Here we comparatively analyze the plasma pressure anisotropies within discontinuities and adjacent Alfvénic fluctuations, leveraging unique solar wind observations from orbit conjunctions between the ARTEMIS and WIND missions, along the same flow streamline, though about 150 Earth radii apart. Based on 11 cases of such observations, we compare direct measurements of plasma anisotropy from particle instruments with its estimates from anisotropic MHD theory using the ratios of correlated ion velocity and Alfvén speed variations Δv_i/Δv_A. We find that (1) sporadically observed discontinuities associated with bifurcated reconnection current sheets harbor significant parallel electron anisotropies of >0.2; (2) direct electron measurements in all events reveal a median anisotropy of ∼0.07 for Alfvénic fluctuations and ∼0.17 for discontinuities; (3) anisotropic MHD predicts even more disparate total anisotropies within Alfvénic fluctuations and discontinuities, with a median value of ∼0.15 for the former and ∼0.57 for the latter; (4) the differences between theory-predicted and directly measured anisotropies imply that the ion contribution to anisotropy is significant and likely dominant within both types of structures, an assertion which we partly verify using simultaneous ion measurements from WIND. Our observations confirm that such discontinuities play a uniquely important role in producing solar wind plasma heating and anisotropy.

Determining how the galactic environment, especially the high gas densities and complex dynamics in bar-fed galaxy centers, alters the star formation efficiency (SFE) of molecular gas is critical to understanding galaxy evolution. However, these same physical or dynamical effects also alter the emissivity properties of CO, leading to variations in the CO-to-H₂ conversion factor (α_CO) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence of α_CO on the local CO velocity dispersion at 150 pc scales using a new set of dust-based α_CO measurements and propose a new α_CO prescription that accounts for CO emissivity variations across galaxies. Based on this prescription, we estimate the SFE in a sample of 65 galaxies from the PHANGS–Atacama Large Millimeter/submillimeter Array survey. We find increasing SFE toward high-surface-density regions like galaxy centers, while using a constant or metallicity-based α_CO results in a more homogeneous SFE throughout the centers and disks. Our prescription further reveals a mean molecular gas depletion time of 700 Myr in the centers of barred galaxies, which is overall three to four times shorter than in nonbarred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2–3 Gyr, regardless of the choice of α_CO prescription. All together, our results suggest that the high level of star formation activity in barred centers is not simply due to an increased amount of molecular gas, but also to an enhanced SFE compared to nonbarred centers or disk regions.

We calculate cross sections for fine-structure transitions of Ne⁺, Ar⁺, Ne²⁺, and Ar²⁺ in collisions with atomic hydrogen by using quantum-mechanical methods. Relaxation rate coefficients are calculated for temperatures up to 10,000 K. The temperature-dependent critical densities for the relaxation of Ne⁺, Ar⁺, Ne²⁺, and Ar²⁺ in collisions with H have been determined and compared to the critical densities for collisions with electrons. The present calculations will be useful for studies utilizing the infrared lines [Ne ii] 12.8, [Ne iii] 15.6, [Ne iii] 36.0, [Ar ii] 6.99, [Ar iii] 8.99, and [Ar iii] 21.8 μm as diagnostics of, for example, planetary nebulae and star formation.

We present high-resolution 1.5–6 GHz Karl G. Jansky Very Large Array and Hubble Space Telescope (HST) optical and infrared observations of the extremely active repeating fast radio burst (FRB) FRB 20201124A and its barred spiral host galaxy. We constrain the location and morphology of star formation in the host and search for a persistent radio source (PRS) coincident with FRB 20201124A. We resolve the morphology of the radio emission across all frequency bands and measure a star formation rate (SFR) ≈ 8.9 M_⊙ yr⁻¹, approximately ≈2.5–6 times larger than optically inferred SFRs, demonstrating dust-obscured star formation throughout the host. Compared to a sample of all known FRB hosts with radio emission, the host of FRB 20201124A has the most significantly obscured star formation. While HST observations show the FRB to be offset from the bar or spiral arms, the radio emission extends to the FRB location. We propose that the FRB progenitor could have formed in situ (e.g., a magnetar born from a massive star explosion). It is still plausible, although less likely, that the progenitor of FRB 20201124A migrated from the central bar of the host. We further place a limit on the luminosity of a putative PRS at the FRB position of L_6.0GHz ≲ 1.8 ×10²⁷ erg s⁻¹ Hz⁻¹, among the deepest PRS luminosity limits to date. However, this limit is still broadly consistent with both magnetar nebulae and hypernebulae models assuming a constant energy injection rate of the magnetar and an age of ≳10⁵ yr in each model, respectively.

The physical origin of fast radio bursts (FRBs) is still unclear. However, young magnetars associated with short-duration gamma-ray bursts (SGRBs) have been thought to be possible central engines for some FRBs. In this paper, we perform a systematic search for SGRBs that are associated with FRBs in a sample including 623 FRBs (601 one-off bursts and 22 repeaters) and 168 SGRBs with precise localizations. We find that FRB 190309A is spatially associated with GRB 060502B, with a chance probability of 0.05 when temporal and redshift information is taken into account. Considering the high chance probability (the statistical significance is <3σ), we examine other observational properties such as the host galaxy, the dispersion measure, and the energy budget of the central engine to check the possibility of their association. Although the available observational information is insufficient to determine whether they are physically associated, it does not rule out such a possibility. As the only pair of FRB and GRB that are spatially associated, it remains an interesting case worthy of further attention.

Galaxies are observed to be lopsided, meaning that they are more massive and more extended along one side than the opposite side. In this work, we provide a statistical analysis of the lopsided morphology of 1780 isolated satellite galaxies generated by the TNG50-1 simulation, incorporating the effect of tidal fields from halo centers. The isolated satellites are galaxies without nearby substructures whose mass is over 1% of the satellites within their virial radii. We study the radial alignment (RA) between the major axes of satellites and the radial direction of their halo centers in radial ranges of 0–2, 2–5, and 5–10 R_h, with R_h being the stellar half-mass radius. According to our results, the RA is virtually undetectable in inner and intermediate regions, yet it is significantly evident in outer regions. We also calculate the far-to-near-side semiaxial ratios of the major axes, denoted by a₋/a₊, which measure the semiaxial ratios of the major axes in the hemispheres between those facing away from (far side) and facing toward (nearside) halo centers. In all the radial bins of the satellites, the numbers of satellites with longer semiaxes on the far side are found to be almost equal to those with longer semiaxes on the near side. Therefore, the tidal fields from halo centers play a minor role in the generation of lopsided satellites. The long semimajor-axes radial alignment (LRA), i.e., an alignment between the long semimajor axes of satellite galaxies and the radial directions to their halo centers, is further studied. No clear evidence of LRA is found in our sample within the framework of ΛCDM Newtonian dynamics. Finally, we briefly discuss the possible origins of the asymmetry of galaxies in TNG50-1.

Some supernovae (SNe) are powered by the collision of the SN ejecta with dense circumstellar matter (CSM). Their emission spectra show characteristic line shapes of combined broad emission and narrow P Cygni lines, which should closely relate to the CSM structure and the mass-loss mechanism that creates the dense CSM. We quantitatively investigate the relationship between the line shape and the CSM structure by Monte Carlo radiative transfer simulations, considering two representative cases of dense CSM formed by steady and eruptive mass loss. Comparing the Hα emission between the two cases, we find that a narrow P Cygni line appears in the eruptive case but does not appear in the steady case due to the difference in the velocity gradient in the dense CSM. We also reproduce the blueshifted photon excess observed in some Type IIn SNe, which is formed by photon transport across the shock wave, and find the relationship between the velocity of the shocked matter and the amount of blueshift of the photon excess. We conclude that the presence or absence of narrow P Cygni lines can distinguish the mass-loss mechanism and suggest high-resolution spectroscopic observations with λ/Δλ ≳ 10⁴ after the light-curve peak for applying this diagnostic method.

We report detection results of nine millisecond pulsars (MSPs) at 8600 MHz using simultaneous 2250 and 8600 MHz observations conducted with the Shanghai Tian Ma Radio Telescope. Mainly benefiting from updated ephemerids with 2250 MHz observations, integrated profiles of all nine MSPs at 8600 MHz are successfully obtained by coherently adding multi-epoch (3–83 epochs) observation data spanning from 19–1210 days, which significantly increases the number of MSPs with published profiles (from 4 to 11) above 8000 MHz, as seven of our target MPSs had no related results previously. Combining our new flux density and pulse width measurements with previous low-frequency results, we study their integrated profile evolution and spectral behaviors in a wider frequency range. We find their component separations and pulse widths remain almost constant, which is consistent with previous findings. While dramatic evolution of integrated profiles exists at low frequencies, we observe a potential end of the related evolution around 5 GHz in eight MSPs. The spectra of four MSPs are found to deviate from a single power law, and we fit them with a broken power law. The change in the profile of PSR J1713+0747, which started around MJD 59320−59321, seems to be more prominent as the observation frequency increases. Compared with the effects of the interstellar medium, we prefer to explain this event as some changes in the magnetosphere. We also find its integrated profile possibly had not recovered to the pre-event state until MJD 59842–59857.

We  present the demography of the dynamics and gas mass fraction of 33 extremely metal-poor galaxies (EMPGs) with metallicities of 0.015–0.195 Z_⊙ and low stellar masses of 10⁴–10⁸M_⊙ in the local universe. We conduct deep optical integral field spectroscopy (IFS) for the low-mass EMPGs with the medium-high resolution (R = 7500) grism of the 8 m Subaru FOCAS IFU instrument by the EMPRESS 3D survey, and investigate the Hα emission of the EMPGs. Exploiting the resolution high enough for the low-mass galaxies, we derive gas dynamics with the Hα lines by the fitting of three-dimensional disk models. We obtain an average maximum rotation velocity (v_rot) of 15 ± 3 km s⁻¹ and an average intrinsic velocity dispersion (σ₀) of 27 ± 10 km s⁻¹ for 15 spatially resolved EMPGs out of 33 EMPGs, and find that all 15 EMPGs have v_rot/σ₀ < 1 suggesting dispersion-dominated systems. There is a clear decreasing trend of v_rot/σ₀ with the decreasing stellar mass and metallicity. We derive the gas mass fraction (f_gas) for all 33 EMPGs, and find no clear dependence on stellar mass and metallicity. These v_rot/σ₀ and f_gas trends should be compared with young high-z galaxies observed by the forthcoming JWST IFS programs to understand the physical origins of the EMPGs in the local universe.

A recent ultraviolet luminosity function (UVLF) analysis in the Hubble Frontier Fields, behind foreground lensing clusters, has helped solidify estimates of the faint-end of the z ∼ 5–9 UVLF at up to 5 mag fainter than in the field. These measurements provide valuable information regarding the role of low-luminosity galaxies in reionizing the universe and can help in calibrating expectations for JWST observations. We fit a semiempirical model to the lensed and previous UVLF data from Hubble. This fit constrains the average star formation efficiency (SFE) during reionization, with the lensed UVLF measurements probing halo mass scales as small as M ∼ 2 × 10⁹M_⊙. The implied trend of SFE with halo mass is broadly consistent with an extrapolation from previous inferences at M ≳ 10¹⁰M_⊙, although the joint data prefer a shallower SFE. This preference, however, is partly subject to systematic uncertainties in the lensed measurements. Near z ∼ 6, we find that the SFE peaks at ∼20% between ∼10¹¹ and 10¹²M_⊙. Our best-fit model is consistent with the Planck 2020 determinations of the electron scattering optical depth, and most current reionization history measurements, provided the escape fraction of ionizing photons is f_esc ∼ 10%–20%. The joint UVLF accounts for nearly 80% of the ionizing photon budget at z ∼ 8. Finally, we show that recent JWST UVLF estimates at z ≳ 11 require strong departures from the redshift evolution suggested by the Hubble data.

Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at z ≥ 3. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with the physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample (M_*/M_⊙ > 10⁹) at 3 ≤ z ≤ 6 from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not exactly populate the same manifold in the representation space from contrastive learning. We also find that half the galaxies classified as disks—either convolutional neural network-based or visually—populate a similar region of the representation space as TNG50 galaxies with low stellar specific angular momentum and nonoblate structure. Although our data-driven study does not allow us to firmly conclude on the true nature of these galaxies, it suggests that the disk fraction at z ≥ 3 remains uncertain and possibly overestimated by traditional supervised classifications. Deeper imaging and spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies, and establish more robust constraints on the emergence of disks at very high redshift.

Correlating the Shack–Hartmann wave front sensor (SH-WFS) with extended targets is widely used in solar adaptive optics systems. This paper aims to introduce a theoretical analysis that evaluates the accuracy of the SH-WFS on extended sources, with a specific focus on the implementation of the Normalized Cross-correlation (NCC) algorithm. To obtain an accurate error description, we utilized the calculation formula of the NCC algorithm to directly express the coordinates of the maximum value in the correlation function matrix. Furthermore, we determined the variance of the centroid position through the error transfer function, which quantifies the measurement error. In comparison with the previous findings of Michau et al., our result exhibits a coefficient disparity, specifically obtaining results 1.5 times higher than their work. The extensive solar granulation simulation and experimental results validate the theoretical error formulas. These error formulas can effectively estimate the accuracy of the SH-WFS, providing a theoretical foundation for the design of optical systems.

Galaxy formation and evolution involve a variety of effectively stochastic processes that operate over different timescales. The extended regulator model provides an analytic framework for the resulting variability (or "burstiness") in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow, equilibrium, and dynamical processes influencing giant molecular clouds' creation and destruction using the power spectral density (PSD) formalism. We use the connection between the PSD and autocovariance function for general stochastic processes to reformulate this model as an autocovariance function, which we use to model variability in galaxy star formation histories (SFHs) using physically motivated Gaussian processes in log star formation rate (SFR) space. Using stellar population synthesis models, we then explore how changes in model stochasticity can affect spectral signatures across galaxy populations with properties similar to the Milky Way and present-day dwarfs, as well as at higher redshifts. We find that, even at fixed scatter, perturbations to the stochasticity model (changing timescales vs. overall variability) leave unique spectral signatures across both idealized and more realistic galaxy populations. Distributions of spectral features including Hα and UV-based SFR indicators, Hδ and Ca H and K absorption-line strengths, D_n(4000), and broadband colors provide testable predictions for galaxy populations from present and upcoming surveys with the Hubble Space Telescope, James Webb Space Telescope, and Nancy Grace Roman Space Telescope. The Gaussian process SFH framework provides a fast, flexible implementation of physical covariance models for the next generation of spectral energy distribution modeling tools. Code to reproduce our results can be found at https://github.com/kartheikiyer/GP-SFH.

We present the results of an ultradeep radio continuum survey, containing ∼480 hr of observations, of the Galactic globular cluster 47 Tucanae with the Australia Telescope Compact Array. This comprehensive coverage of the cluster allows us to reach rms noise levels of 1.19 μJy beam⁻¹ at 5.5 GHz, 940 nJy beam⁻¹ at 9 GHz, and 790 nJy beam⁻¹ in a stacked 7.25 GHz image. This is the deepest radio image of a globular cluster and the deepest image ever made with the Australia Telescope Compact Array. We identify ATCA J002405.702-720452.361, a faint (6.3 ± 1.2 μJy at 5.5 GHz, 5.4 ± 0.9 μJy at 9 GHz), flat-spectrum (α = −0.31 ± 0.54) radio source that is positionally coincident with the cluster center and potentially associated with a faint X-ray source. No convincing optical counterpart was identified. We use radio, X-ray, optical, and UV data to show that explanations involving a background active galactic nucleus, a chromospherically active binary, or a binary involving a white dwarf are unlikely. The most plausible explanations are that the source is an undiscovered millisecond pulsar or a weakly accreting black hole. If the X-ray source is associated with the radio source, the fundamental plane of black-hole activity suggests a black hole mass of ∼54–6000 M_⊙, indicating an intermediate-mass black hole or a heavy stellar-mass black hole.

We report the discovery of three Transiting Exoplanet Survey Satellite Objects of Interest (TOI) with signatures of pulsation, observed in more than one sector. Our main goal is to explore how large is the variety of classical pulsators such as δ Sct, γ Dor, RR Lyrae and Cepheid among TOI pulsators. The analysis reveals two stars with signatures of δ Sct and one of γ Dor, out of a sample of 3901 TOIs with available light curves (LCs). To date, there is a very scarce number of known pulsating stars hosting planets. The present finding also emerges as an exciting laboratory for studying different astrophysical phenomena, including the effects of star–planet interaction on pulsation and timing detection of planetary companions. We have also identified 16 TOI stars with periodicities and LCs morphology compatible with different classical pulsating classes, but for most of them, the dominant frequency signals originate from contaminating sources.

The Mapper of the IGM Spin Temperature (MIST) is a new ground-based, single-antenna, radio experiment attempting to detect the global 21 cm signal from the Dark Ages and Cosmic Dawn. A significant challenge in this measurement is the frequency dependence, or chromaticity, of the antenna beam directivity. MIST observes with the antenna above the soil and without a metal ground plane, and the beam directivity is sensitive to the electrical characteristics of the soil. In this paper, we use simulated observations with MIST to study how the detection of the global 21 cm signal from Cosmic Dawn is affected by the soil and the MIST beam directivity. We simulate observations using electromagnetic models of the directivity computed for single- and two-layer models of the soil. We test the recovery of the Cosmic Dawn signal with and without beam chromaticity correction applied to the simulated data. We find that our single-layer soil models enable a straightforward recovery of the signal even without chromaticity correction. Two-layer models increase the beam chromaticity and make the recovery more challenging. However, for the model in which the bottom soil layer has a lower electrical conductivity than the top layer, the signal can be recovered even without chromaticity correction. For the other two-layer models, chromaticity correction is necessary for the recovery of the signal, and the accuracy requirements for the soil parameters vary between models. These results will be used as a guideline to select observation sites that are favorable for the detection of the Cosmic Dawn signal.

In this paper we explore the idea of using multi-spacecraft observations of Jovian electrons to measure the 3D distribution of these particles in the inner heliosphere. We present simulations of Jovian electron intensities along selected spacecraft trajectories for 2021 and compare these, admittedly qualitatively, to these measurements. Using the data–model comparison we emphasize how such a study can be used to constrain the transport parameters in the inner heliosphere, and how this can lead to additional insight into energetic particle transport. Model results are also shown along the expected trajectories of selected spacecraft, including the off-ecliptic phase of the Solar Orbiter mission from 2025 onward. Lastly, we revisit the use of historical data and discuss upcoming missions that may contribute to Jovian electron measurements.

Peptide-like molecules, which have a close connection with the origin of life, have been detected in the Universe. Mapping observations of HCONH₂ and CH₃CONH₂, two of the simplest peptide-like molecules, are performed toward the Sagittarius B2 (Sgr B2) complex with the IRAM 30 m telescope. Seven transitions of HCONH₂ and five transitions of CH₃CONH₂ are used in the analysis. The spatial distributions of the excitation temperature and column density of HCONH₂ in the molecular envelope of Sgr B2 are obtained by rotation diagrams. Assuming the same excitation temperature of HCONH₂, the column densities of CH₃CONH₂ are also calculated. The results show that the excitation temperature ranges from 6 to 46 K in the molecular envelope of Sgr B2. The abundance ratios between HCONH₂ and CH₃CONH₂ are calculated to explore the relationship between them, as are those between HCONH₂ and HNCO. The abundance ratio of CH₃CONH₂/HCONH₂ varies from 10% to 20%, while that of HCONH₂/HNCO ranges from 1.5% to 10%. CH₃CONH₂ is enhanced with respect to HCONH₂ in the northwest region of Sgr B2. One transition of H¹³CONH₂ is detected toward 12 positions of Sgr B2, from which a ¹²C/¹³C ratio of 28.7 is obtained. A time-dependent chemical model with a short-duration X-ray burst is used to explain the observed abundances of HCONH₂ and CH₃CONH₂, with the best-fitting result at T_dust = 53–56 K. More chemical reactions are required to be included in the model since the modeled abundance is lower than the observed one at the observed T_dust.

We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including the Dark Energy Survey (DES), DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of z = 2.35 to a depth of m_r = 27. Compared to the Buzzard mocks, the Cardinal mocks include an updated subhalo abundance matching model that considers orphan galaxies and includes mass-dependent scatter between galaxy luminosity and halo properties. This model can simultaneously fit galaxy clustering and group–galaxy cross-correlations measured in three different luminosity threshold samples. The Cardinal mocks also feature a new color assignment model that can simultaneously fit color-dependent galaxy clustering in three different luminosity bins. We have developed an algorithm that uses photometric data to further improve the color assignment model and have also developed a novel method to improve small-scale lensing below the ray-tracing resolution. These improvements enable the Cardinal mocks to accurately reproduce the abundance of galaxy clusters and the properties of lens galaxies in the DES data. As such, these simulations will be a valuable tool for future cosmological analyses based on large sky surveys.

The coronal heating problem has been a major challenge in solar physics, and a tremendous amount of effort has been made over the past several decades to solve it. In this paper, we aim at answering how the physical processes behind the Alfvén wave turbulent heating adopted in the Alfvén Wave Solar atmosphere Model (AWSoM) unfold in individual plasma loops in an active region (AR). We perform comprehensive investigations in a statistical manner on the wave dissipation and reflection, temperature distribution, heating scaling laws, and energy balance along the loops, providing in-depth insights into the energy allocation in the lower solar atmosphere. We demonstrate that our 3D global model with a physics-based phenomenological formulation for the Alfvén wave turbulent heating yields a heating rate exponentially decreasing from loop footpoints to top, which had been empirically assumed in the past literature. A detailed differential emission measure (DEM) analysis of the AR is also performed, and the simulation compares favorably with DEM curves obtained from Hinode/Extreme-ultraviolet Imaging Spectrometer observations. This is the first work to examine the detailed AR energetics of our AWSoM model with high numerical resolution and further demonstrates the capabilities of low-frequency Alfvén wave turbulent heating in producing realistic plasma properties and energetics in an AR.

The very-high-energy γ-ray source HESS J1809-193 has been detected by the LHAASO and HAWC observatory beyond 100 TeV energy. It is an interesting candidate for exploring the underlying mechanisms of γ-ray production due to the presence of supernova remnants, pulsars, and molecular clouds close to it. We have considered the injection of the energetic cosmic rays from a past explosion, whose reminiscent may be SNR G011.0-00.0, which is located within the extended γ-ray source HESS J1809-193. We explain the multiwavelength data from the region of HESS J1809-193 with synchrotron, inverse Compton, and bremsstrahlung emission of cosmic-ray electrons and secondary γ-ray production in interactions of cosmic-ray protons with the cold protons in the local molecular clouds within a time-dependent framework including the diffusion loss of cosmic rays. The observational data have been modeled with the secondary photons produced by the time-evolved cosmic-ray spectrum assuming the age of the explosion is 4500 yr.

In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state characterizing dense matter. Two independent analyses found a mass of ∼1.3–1.4 M_⊙ and a radius of ∼13 km. They also both found that the hot spots were all located on the same hemisphere, opposite to the observer, and that at least one of them had a significantly elongated shape. Here we reanalyze, in greater detail, the same NICER data set, incorporating the effects of an updated NICER response matrix and using an upgraded analysis framework. We expand the adopted models and also jointly analyze XMM-Newton data, which enables us to better constrain the fraction of observed counts coming from PSR J0030+0451. Adopting the same models used in previous publications, we find consistent results, although with more stringent inference requirements. We also find a multimodal structure in the posterior surface. This becomes crucial when XMM-Newton data is accounted for. Including the corresponding constraints disfavors the main solutions found previously, in favor of the new and more complex models. These have inferred masses and radii of ∼[1.4 M_⊙, 11.5 km] and ∼[1.7 M_⊙, 14.5 km], depending on the assumed model. They display configurations that do not require the two hot spots generating the observed X-rays to be on the same hemisphere, nor to show very elongated features, and point instead to the presence of temperature gradients and the need to account for them.

In this paper, we present a Lyα halo (LAH) identified by stacking ∼3300 Lyα emitters (LAEs) at z = 2.2–2.3. We carry out imaging observations and data reduction with Subaru/Hyper Suprime-Cam. Our total survey area is ∼12 deg² and the imaging depths are 25.5–27.0 mag. Using the imaging data, we select 1240 and 2101 LAE candidates at z = 2.2 and 2.3, respectively. We carry out spectroscopic observations of our LAE candidates and data reduction with Magellan/IMACS to estimate the contamination rate of our LAE candidates. We find that the contamination rate of our sample is low (8%). We stack our LAE candidates with a median stacking method to identify the LAH at z = 2. We show that our LAH is detected until ∼100 kpc at the 2σ significance level and likely extended to ∼200 kpc at a surface brightness level of ∼10⁻²⁰ erg s⁻¹ cm⁻² arcsec⁻². Compared to those of previous studies, our LAH is brighter at radii of ∼25–100 kpc, which is not likely caused by the contamination in our sample but by the different redshifts, fields, and selection methods instead. To investigate how central galaxies affect surrounding LAHs, we divide our LAEs into subsamples based on the Lyα luminosity (L_Lyα), rest-frame Lyα equivalent width (EW₀), and UV magnitude (M_uv). We stack the subsamples and find that higher L_Lyα, smaller EW₀, and brighter M_uv cause more extended halos. Our results suggest that more massive LAEs generally have more extended LAHs.

Previous work has established an empirical relationship between densities gained from coronal rotational tomography near the ecliptic plane with solar wind outflow speeds at heliocentric distance r₀ = 8R_⊙. This work aims to include solar wind acceleration, and thus velocity profiles out to 1 au. Inner boundary velocities are given as a function of normalized tomographic densities, ρ_N, as ${V}_{0}=\left(75\ast {e}^{-\left[5.2* {\rho }_{N}\right]}+108\right)$, and typically range from 100 to 180 km s⁻¹. The subsequent acceleration is defined as $V(r)={V}_{0}\left(1+{\alpha }_{\mathrm{IP}}\left[1-{e}^{\left(-\left[r-{r}_{0}\right]/{r}_{{\rm{H}}}\right)}\right]\right)$, with α_IP ranging between 1.75 and 2.7, and r_H between 50 and 35 R_⊙ dependent on V₀. These acceleration profiles approximate the distribution of in situ measurements by Parker Solar Probe (PSP) and other measurements at 1 au. Between 2018 November and 2021 September these constraints are applied using the HUXt model and give good agreement with in situ observations at PSP, with a ∼6% improvement compared with using a simpler constant acceleration model previously considered. Given the known tomographical densities at 8 R_⊙, we extrapolate density to 1 au using the model velocities and assuming mass flux conservation. Extrapolated densities agree well with OMNI measurements. Thus coronagraph-based estimates of densities define the ambient solar wind outflow speed, acceleration, and density from 8 R_⊙ to at least 1 au. This sets a constraint on more advanced models, and a framework for forecasting that provides a valid alternative to the use of velocities derived from magnetic field extrapolations.

In the epoch of deep photometric surveys, a large number of substructures—e.g., overdensities and streams—have been identified. With the help of astrometry and spectroscopy, the community has revealed a complex picture of our Milky Way (MW) after investigating their origins. The off-plane substructures the Anticenter Stream (ACS) and Monoceros Ring (MNC), once considered as dissolving dwarf galaxies, were later found to share similar kinematics and metallicity with the Galactic outer thin disk. In this work, we aim to chemically tag ACS and MNC with high-accuracy abundances from the APOGEE survey. By extrapolating chemical abundance trends in the outer thin-disk region (10 < R_GC < 18 kpc, 0 < ∣Z_GC∣ < 3 kpc), we found that ACS and MNC stars show consistent chemical abundances as the extrapolating values for 12 elements, including C, N, O, Mg, Al, Si, K, Ca, Cr, Mn, Co, and Ni. The similar chemical patterns indicate that ACS and MNC have a similar star formation history as the MW outer thin disk, while we also excluded their dwarf galaxy association, as they are distinctive in multiple chemical spaces. The ages of ACS and MNC stars are consistent with the time of the first Sagittarius dSph passage, indicating their possible connection.

Inhomogeneous cloud formation and wavelength-dependent phenomena are expected to shape hot Jupiter atmospheres. We present a general circulation model with multiwavelength "picket fence" radiative transfer and radiatively active, temperature-dependent clouds, and compare the results to those of a double gray routine. The double gray method inherently fails to model polychromatic effects in hot Jupiter atmospheres, while picket fence captures these non-gray aspects and performs well compared to fully wavelength-dependent methods. We compare both methods with radiatively active clouds and cloud-free models, assessing the limitations of the double gray method. Although there are broad similarities, the picket fence models have larger dayside–nightside temperature differences, nonisothermal upper atmospheres, and multiwavelength effects in the presence of radiatively active clouds. We model the well-known hot Jupiters HD 189733 b and HD 209458 b. For the hotter HD 209458 b, the picket fence method prevents clouds from thermostating dayside temperatures, resulting in hotter upper atmospheres and the dissipation of dayside clouds. Differences in the temperature structures are then associated with nuanced differences in the circulation patterns and clouds. Models of the cooler HD 189733 b have global cloud coverage, regardless of the radiative transfer scheme, whereas there are larger differences in the models of HD 209458 b, particularly in the extent of the partial cloud coverage on its dayside. This results in minor changes to the thermal and reflected light phase curves of HD 189733 b, but more significant differences for the picket fence and double gray versions of HD 209458 b.

We present the updated open-source code Complete History of Interaction-Powered Supernovae (CHIPS) that can be applied to modeling supernovae (SNe) arising from an interaction with the massive circumstellar medium (CSM) as well as the formation process of the CSM. Our update mainly concerns extensions to hydrogen-poor SNe from stripped progenitors, targeting the modeling of interaction-powered SNe Ibc, such as Type Ibn and Icn SNe. We successfully reproduce the basic properties of the light curves of these types of SNe that occur after the partial eruption of the outermost layer with a mass of 0.01–0.1 M_⊙ at ≲1 year before explosion. We also find that the luminosity of the observed precursors can be naturally explained by the outburst that creates the dense CSM, given that the energy of the outburst is efficiently dissipated by collision with an external material, possibly generated by a previous mass eruption. We discuss possible scenarios causing eruptive mass loss based on our results.

In this study, we examine photoionization outflows during the late stages of galaxy mergers, with a specific focus on the relation between the observed velocity of outflowing gas and the apparent effects of dust extinction. We used the N-body/smoothed particle hydrodynamics code ASURA for galaxy merger simulations. These simulations concentrated on identical galaxy mergers featuring supermassive black holes of 10⁸M_⊙ and gas fractions of 30% and 10%. From the simulation data, we derived velocity and velocity dispersion diagrams for the active galactic nuclei (AGN)-driven ionized outflowing gas. Our findings show that high-velocity outflows with velocity dispersions of 500 km s⁻¹ or greater can be observed in the late stages of galactic mergers. Particularly, in buried AGNs, both the luminosity-weighted outflow velocity and velocity dispersion increase owing to the apparent effects of dust extinction. Owing to these effects, velocity–velocity dispersion diagrams display a noticeable blue-shifted tilt in models with higher gas fractions. Crucially, this tilt is not influenced by the AGN luminosity but emerges from the observational impacts of dust extinction. Our results imply that the observed high-velocity [O iii] λ5007 outflow exceeding 1000 km s⁻¹ in buried AGNs may be linked to the dust extinction that occurs during the late stages of gas-rich galaxy mergers.

HDF850.1 is the brightest submillimeter galaxy (SMG) in the Hubble Deep Field. It is known as a heavily dust-obscured star-forming galaxy embedded in an overdense environment at z = 5.18. With nine-band NIRCam images at 0.8–5.0 μm obtained through the JWST Advanced Deep Extragalactic Survey, we detect and resolve the rest-frame UV–optical counterpart of HDF850.1, which splits into two components because of heavy dust obscuration in the center. The southern component leaks UV and Hα photons, bringing the galaxy ∼100 times above the empirical relation between infrared excess and UV continuum slope (IRX–β_UV). The northern component is higher in dust attenuation and thus fainter in UV and Hα surface brightness. We construct a spatially resolved dust-attenuation map from the NIRCam images, well matched with the dust continuum emission obtained through millimeter interferometry. The whole system hosts a stellar mass of 10^10.8±0.1M_⊙ and star formation rate (SFR) of 10^2.8±0.2M_⊙ yr⁻¹, placing the galaxy at the massive end of the star-forming main sequence at this epoch. We further confirm that HDF850.1 resides in a complex overdense environment at z = 5.17–5.30, which hosts another luminous SMG at z = 5.30 (GN10). The filamentary structures of the overdensity are characterized by 109 Hα-emitting galaxies confirmed through NIRCam slitless spectroscopy at 3.9–5 μm, of which only eight were known before the JWST observations. Given the existence of a similar galaxy overdensity in the GOODS-S field, our results suggest that 50% ± 20% of the cosmic star formation at z = 5.1–5.5 occur in protocluster environments.

The discovery and timing follow up of millisecond pulsars (MSPs) are necessary not just for their usefulness in pulsar timing arrays (PTAs) but also for investigating their own intriguing properties. In this work, we provide the findings of the decade-long timing of four MSPs discovered by the Giant Meterwave Radio Telescope (GMRT), including their timing precision, model parameters, and newly detected proper motions. We compare the timing results for these MSPs before and after the GMRT upgrade in 2017 and characterize the improvement in timing precision due to the bandwidth upgrade. We discuss the suitability of these four GMRT MSPs as well as the usefulness of the decade-long timing data for PTA experiments. These data may aid in the global effort to improve the signal-to-noise ratios of recently detected signature of gravitational waves in cross-correlation statistics of residuals of MSPs.

We present results of [C ii] 158 μm emission line observations, and report the spectroscopic redshift confirmation of a strongly lensed (μ ∼ 20) star-forming galaxy, MACS0308-zD1 at z = 6.2078 ± 0.0002. The [C ii] emission line is detected with a signal-to-noise ratio >6 within the rest-frame UV-bright clump of the lensed galaxy (zD1.1) and exhibits multiple velocity components; the narrow [C ii] has a velocity full width half maximum (FWHM) of 110 ± 20 km s⁻¹, while broader [C ii] is seen with an FWHM of 230 ± 50 km s⁻¹. The broader [C ii] component is blueshifted (−80 ± 20 km s⁻¹) with respect to the narrow [C ii] component, and has a morphology that extends beyond the UV-bright clump. We find that, while the narrow [C ii] emission is most likely associated with zD1.1, the broader component is possibly associated with a physically distinct gas component from zD1.1 (e.g., outflowing or inflowing gas). Based on the nondetection of λ_158μm dust continuum, we find that MACS0308-zD1's star formation activity occurs in a dust-free environment indicated by a strong upper limit of infrared luminosity ≲9 × 10⁸L_⊙. Targeting this strongly lensed faint galaxy for follow-up Atacama Large Millimeter/submillimeter Array and JWST observations will be crucial to characterize the details of typical galaxy growth in the early Universe.

HD 148937 is a peculiar massive star (Of?p) with a strong magnetic field (1 kG). The hourglass-shaped emission nebula NGC 6164/5 surrounds this star. This nebula is presumed to originate from episodic mass-loss events of the central O-type star, but the detailed formation mechanism is not yet well understood. Grasping its three-dimensional structure is essential to uncovering the origin of this nebula. Here we report the high-resolution multiobject spectroscopic observations of NGC 6164/5 using the GIRAFFE on the 8.2 m Very Large Telescope. Integrated intensity maps constructed from several spectral lines delineate well the overall shape of this nebula, such as the two bright lobes and the inner gas region. The position–velocity diagrams show that the two bright lobes are found to be redshifted and blueshifted, respectively, while the inner region has multiple layers. We consider a geometric model composed of a bilateral outflow harboring nitrogen-enriched knots and expanding inner shells. Its spectral features are then simulated by using a Monte Carlo radiative transfer technique for different sets of velocities. Some position–velocity diagrams from simulations are very similar to the observed ones. According to the model that best reproduces the observational data, the two bright lobes and the nitrogen-enriched knots are moving away from HD 148937 at about 120 km s⁻¹. Their minimum kinematic age is estimated to be about 7500 yr. We discuss possible formation mechanisms of this nebula in the context of binary interaction.

We demonstrate that the inference of galaxy stellar masses via spectral energy distribution (SED) fitting techniques for galaxies formed in the first billion years after the Big Bang carries fundamental uncertainties owing to the loss of star formation history (SFH) information from the very first episodes of star formation in the integrated spectra of galaxies. While this early star formation can contribute substantially to the total stellar mass of high-redshift systems, ongoing star formation at the time of detection outshines the residual light from earlier bursts, hampering the determination of accurate stellar masses. As a result, order-of-magnitude uncertainties in stellar masses can be expected. We demonstrate this potential problem via direct numerical simulation of galaxy formation in a cosmological context. In detail, we carry out two cosmological simulations with significantly different stellar feedback models, which span a significant range in SFH burstiness. We compute the mock SEDs for these model galaxies at z = 7 via calculations of 3D dust radiative transfer, and then backward fit these SEDs with prospector SED fitting software. The uncertainties in derived stellar masses that we find for z > 7 galaxies motivate the development of new techniques and/or priors for SFH to model star formation in the early Universe.

Galactic conformity is the phenomenon whereby a galaxy of a certain physical property is correlated with its neighbors of the same property, implying a possible causal relationship. The observed auto correlations of emission-line galaxies (ELGs) from the highly complete DESI One-Percent Survey exhibit a strong clustering signal on small scales, providing clear evidence for the conformity effect of ELGs. Building upon the original subhalo abundance-matching (SHAM) method developed by Gao et al., we propose a concise conformity model to improve the ELG–halo connection. In this model, the number of satellite ELGs is boosted by a factor of ∼5 in the halos whose central galaxies are ELGs. We show that the mean ELG satellite number in such central halos is still smaller than 1 and that the model does not significantly increase the overall satellite fraction. With this model, we can well recover the ELG auto correlations to the smallest scales explored with the current data (i.e., r_p > 0.03 Mpc h⁻¹ in real space and at s > 0.3 Mpc h⁻¹ in redshift space), while the cross correlations between luminous red galaxies and ELGs are nearly unchanged. Although our SHAM model has only eight parameters, we further verify that it can accurately describe the ELG clustering in the entire redshift range from z = 0.8 to 1.6. We therefore expect that this method can be used to generate high-quality ELG lightcone mocks for DESI.

Magnetic reconnection, a fundamental plasma process transforming magnetic field energy into particle energy, is ubiquitous in space and responsible for many explosive phenomena, such as solar flares and gamma-ray bursts. Recent numerical theories have predicted that reconnection fronts far from the primary reconnection region can host secondary reconnection in three-dimensional scenarios, different from the conventional two-dimensional diagram where only one X-line stands to sustain reconnection. In this study, we provide direct observational evidence for ongoing secondary reconnection in the reconnection front via the unprecedentedly high-cadence data from NASA's MMS mission. The secondary reconnection is identified by the presence of an X-line, a super-Alfvénic electron jet, and nonideal energy dissipation. Different from the primary ion–electron reconnection, the secondary reconnection is electron-only, with its X-line quasi-perpendicular to the primary X-line. Hence reconnection, when evolving from local to global scales, becomes essentially three-dimensional with different patterns developed. These results provide crucial insights into understanding cross-scale energy transport driven by reconnection in space plasmas.

We present the results of the first fully cosmological hydrodynamical simulations studying the merger-driven model for massive black hole (BH) seed formation via direct collapse. Using the zoom-in technique as well as particle splitting, we achieve a final spatial resolution of 2 pc. We show that the major merger of two massive galaxies at redshift z ∼ 8 results in the formation of a nuclear supermassive disk (SMD) of only 4 pc in radius, owing to a prodigious gas inflow sustained at 100–1000 M_⊙ yr⁻¹. The core of the merger remnant is metal-rich, well above solar abundance, and the SMD reaches a gaseous mass of 3 × 10⁸M_⊙ in less than a million years after the merger, despite a concurrent prominent nuclear starburst. Dynamical heating as gas falls into the deepest part of the potential well, and heating and stirring by supernova blastwaves, generate a turbulent multiphase interstellar medium, with a gas velocity dispersion exceeding 100 km s⁻¹. As a result, only moderate fragmentation occurs in the inner 10–20 pc, despite the temperature falling below 1000 K. The SMD is Jeans-unstable as well as bar-unstable and will collapse further adiabatically, becoming warm and ionized. We show that the SMD, following inevitable contraction, will become general-relativistic-unstable and directly form a supermassive BH of mass in the range 10⁶–10⁸M_⊙, essentially skipping the stage of BH seed formation. These results confirm that mergers between the most massive galaxies at z ∼ 8–10 can naturally explain the rapid emergence of bright high-redshift quasars.

We present resolved images of the inner disk component around HD 141569A using the Magellan adaptive optics system with the Clio2 1–5 μm camera, offering a glimpse of a complex system thought to be in a short evolutionary phase between protoplanetary and debris disk stages. We use a reference star along with the Karhunen–Loéve image projection (KLIP) algorithm for point-spread function subtraction to detect the disk inward to about 024 (∼25 au assuming a distance of 111 pc) at high signal-to-noise ratios at $L^{\prime}$ (3.8 μm), Ls (3.3 μm), and narrowband Ice (3.1 μm). We identify an arc or spiral arm structure at the southeast extremity, consistent with previous studies. We implement forward modeling with a simple disk model within the framework of a Markov Chain Monte Carlo sampler to better constrain the geometrical attributes and photometry using our KLIP-reduced disk images. We then leverage these modeling results to facilitate a comparison of the measured brightness in each passband to find a reduction in scattered light from the disk in the Ice filter, implying significant absorption due to water ice in the dust. Additionally, our best-fit disk models exhibit peak brightness in the southwestern, back-scattering region of the disk, which we suggest to be possible evidence of 3.3 μm polycyclic aromatic hydrocarbon emission. However, we point out the need for additional observations with bluer filters and more complex modeling to confirm these hypotheses.

Meridional circulation regulates the Sun's interior dynamics and magnetism. While it is well accepted that meridional flows are poleward at the Sun's surface, helioseismic observations have yet to provide a definitive answer for the depth at which those flows return to the equator, or the number of circulation cells in depth. Here, we explore the observability of multiple circulation cells stacked in radius. Specifically, we examine the seismic signature of several meridional flow profiles by convolving time–distance averaging kernels with mean flows obtained from a suite of 3D hydrodynamic simulations. At mid and high latitudes, we find that weak flow structures in the deep convection zone can be obscured by signals from the much stronger surface flows. This contamination of 1–2 m s⁻¹ is caused by extended side lobes in the averaging kernels, which produce a spurious equatorward signal with flow speeds that are 1 order of magnitude stronger than the original flow speeds in the simulations. At low latitudes, the flows in the deep layers of the simulations are stronger (>2 m s⁻¹) and multiple cells across the convection zone can produce a sufficiently strong signal to survive the convolution process. Now that meridional flows can be measured over two decades of data, the uncertainties arising from convective noise have fallen to a level where they are comparable in magnitude to the systematic biases caused by nonlocal features in the averaging kernels. Hence, these systematic errors are beginning to influence current helioseismic deductions and need broader consideration.

We used 29 high-resolution line-of-sight magnetograms acquired with the Goode Solar Telescope (GST) in a quiet-Sun area to extrapolate a series of potential field configurations and study their time variations. The study showed that there are regions that consistently exhibit changes in loop connectivity, whereas other vast areas do not show such changes. Analysis of the topological features of the potential fields indicates that the photospheric footprint of the separatrix between open- and closed-loop systems closely matches the roots of rapid blue- and redshifted excursions, which are disk counterparts of type II spicules. There is a tendency for the footpoints of the observed H_α features to be cospatial with the footpoints of the loops that most frequently change their connectivity, while the area occupied by the open fields that did not show any significant and persistent connectivity changes is void of prominent jet and spicular activity. We also detected and tracked magnetic elements using the Southwest Automatic Magnetic Identification Suite and GST magnetograms, which allowed us to construct artificial magnetograms and calculate the corresponding potential field configurations. Analysis of the artificial data showed tendencies similar to those found for the observed data. The present study suggests that a significant amount of chromospheric activity observed in the far wings of the H_α spectral line may be generated by reconnecting closed-loop systems and canopy fields consisting of "open" field lines.

We numerically study the diffusion and scattering of cosmic rays (CRs) together with their acceleration processes in the framework of the modern understanding of magnetohydrodynamic (MHD) turbulence. Based on the properties of compressible MHD turbulence obtained from observations and numerical experiments, we investigate the interaction of CRs with plasma modes. We find that (1) the gyroradius of particles exponentially increases with the acceleration timescale; (2) the momentum diffusion presents the power-law relationship with the gyroradius in the strong turbulence regime, and shows a plateau in the weak turbulence regime implying a stochastic acceleration process; (3) the spatial diffusion is dominated by the parallel diffusion in the sub-Alfvénic regime, while it is dominated by the perpendicular diffusion in the super-Alfvénic one; (4) as for the interaction of CRs with plasma modes, the particle acceleration is dominated by the fast mode in the high β case, while in the low β case, it is dominated by the fast and slow modes; and (5) in the presence of acceleration, magnetosonic modes still play a critical role in the diffusion and scattering processes of CRs, which is in good agreement with earlier theoretical predictions.

Interplanetary magnetic flux ropes (MFRs) are commonly observed structures in the solar wind, categorized as magnetic clouds (MCs) and small-scale MFRs (SMFRs) depending on whether they are associated with coronal mass ejections. We apply machine learning to systematically compare SMFRs, MCs, and ambient solar wind plasma properties. We construct a data set of 3-minute averaged sequential data points of the solar wind's instantaneous bulk fluid plasma properties using about 20 years of measurements from Wind. We label samples by the presence and type of MFRs containing them using a catalog based on Grad–Shafranov (GS) automated detection for SMFRs and NASA's catalog for MCs (with samples in neither labeled non-MFRs). We apply the random forest machine learning algorithm to find which categories can be more easily distinguished and by what features. MCs were distinguished from non-MFRs with an area under the receiver-operator curve (AUC) of 94% and SMFRs with an AUC of 89%, and had distinctive plasma properties. In contrast, while SMFRs were distinguished from non-MFRs with an AUC of 86%, this appears to rely solely on the 〈B〉 > 5 nT threshold applied by the GS catalog. The results indicate that SMFRs have virtually the same plasma properties as the ambient solar wind, unlike the distinct plasma regimes of MCs. We interpret our findings as additional evidence that most SMFRs at 1 au are generated within the solar wind. We also suggest that they should be considered a salient feature of the solar wind's magnetic structure rather than transient events.

Strong evidence is reported to support unobscured broad-line regions (BLRs) in Type-1.9 active galactic nucleus (AGN) SDSS J1241+2602 with reliable broad Hα but no broad Hβ. Commonly, the disappearance of broad Hβ can be explained by the AGN unified model, in which heavily obscured BLRs are expected in Type-1.9 AGNs. Here, based on properties of two kinds of BH masses, the virial BH mass and the BH mass obtained through the M_BH–σ relation, an independent method is proposed to test whether there are unobscured central BLRs in a Type-1.9 AGN. By the reliable measurement of stellar velocity dispersion of about 110 ± 12 km s ⁻¹ through the host galaxy absorption features in SDSS J1241+2602, the BH mass obtained through the M_BH–σ relation is consistent with the virial BH mass (3.43 ± 1.25) × 10⁷ M_⊙ determined through properties of the observed broad Hα without considering the effects of obscurations in SDSS J1241+2602. Meanwhile, if considering heavily obscured BLRs in SDSS J1241+2602, the reddening-corrected virial BH mass is tens of times larger than the value expected from M_BH–σ, leading SDSS J1241+2602 to be an outlier in the M_BH–σ space with a confidence level higher than 5σ. Therefore, unobscured BLRs are preferred in the Type-1.9 AGN SDSS J1241+2602. The results indicate that it is necessary to check whether unobscured central BLRs are common in Type-1.9 AGNs when testing the unified model of AGNs through properties of Type-1.9 AGNs.

In the magnetohydrodynamic (MHD) perspective, the planet's bow shock would disappear when the fast-mode Mach number (M_F) of the solar wind is less than one. Compared to Earth, Mercury is subject to a lower M_F solar wind due to its proximity to the Sun, resulting in a higher possibility of the disappearance of its bow shock. To examine the variability of Mercury's bow shock in response to the solar wind properties, analyses of the observations by the Helios spacecraft at 0.30–0.50 au during 1975–1983, covering solar cycle 21, together with the theoretical solutions and MHD simulations are conducted in this study. Our observational analyses show that more solar wind data with extremely low fast-mode Mach numbers (say, M_F ≤ 1.5) are observed during the rising and maximum phases and are characterized by a significantly low proton number density. It is also found that approximately 35% of the extremely low fast-mode Mach number solar wind events (M_F ≤ 1.5) occur within the main body of interplanetary coronal mass ejections (ICMEs), while about 58% of them are unrelated to ICMEs. Three of these events are selected to demonstrate that the occurrences of the solar wind with M_F ≤ 1.5 may not be necessarily affected by ICMEs. Our theoretical and numerical results indicate that when Mercury encounters the solar wind with M_F ≤ 1.5, its bow shock would move farther away, become flattened, and even disappear. Furthermore, our calculations suggest that Mercury's bow shock would become a slow-mode shock with a concave-upward structure under such extreme solar wind conditions.

The IceCube Neutrino Observatory has been continuously taking data to search for ${ \mathcal O }(0.5\mbox{--}10)$ s long neutrino bursts since 2007. Even if a Galactic core-collapse supernova is optically obscured or collapses to a black hole instead of exploding, it will be detectable via the ${ \mathcal O }(10)$ MeV neutrino burst emitted during the collapse. We discuss a search for such events covering the time between 2008 April 17 and 2019 December 31. Considering the average data taking and analysis uptime of 91.7% after all selection cuts, this is equivalent to 10.735 yr of continuous data taking. In order to test the most conservative neutrino production scenario, the selection cuts were optimized for a model based on an 8.8 solar mass progenitor collapsing to an O–Ne–Mg core. Conservative assumptions on the effects of neutrino oscillations in the exploding star were made. The final selection cut was set to ensure that the probability to detect such a supernova within the Milky Way exceeds 99%. No such neutrino burst was found in the data after performing a blind analysis. Hence, a 90% C.L. upper limit on the rate of core-collapse supernovae out to distances of ≈25 kpc was determined to be 0.23 yr⁻¹. For the more distant Magellanic Clouds, only high neutrino luminosity supernovae will be detectable by IceCube, unless external information on the burst time is available. We determined a model-independent limit by parameterizing the dependence on the neutrino luminosity and the energy spectrum.

Surveys of exoplanet host stars are valuable tools for assessing population level trends in exoplanets, and their outputs can include stellar ages, activity, and rotation periods. We extracted chromospheric activity measurements from the California-Kepler Survey Gaia survey spectra in order to probe connections between stellar activity and fundamental stellar properties. Building on the California Kepler Survey's legacy of 1189 planet host star stellar properties including temperature, surface gravity metallicity, and isochronal age, we add measurements of the Ca ii H and K lines as a proxy for chromospheric activity for 879 planet hosting stars. We used these chromospheric activity measurements to derive stellar rotation periods. We find a discrepancy between photometrically derived and activity-derived rotation periods for stars on the Rossby Ridge. These results support the theory of weakened magnetic braking. We find no evidence for metallicity-dependent activity relations, within the metallicity range of −0.2 to +0.3 dex. With our single epoch spectra we identify stars that are potentially in Maunder minimum–like state using a combination of log(${R}_{\mathrm{HK}}^{{\prime} }$) and position below the main sequence. We do not yet have the multiyear time series needed to verify stars in Maunder minimum–like states. These results can help inform future theoretical studies that explore the relationship between stellar activity, stellar rotation, and magnetic dynamos.

Planets can excite density waves and open annular gas gaps in protoplanetary disks. The depth of gaps is influenced by the evolving angular momentum carried by density waves. While the impact of radiative cooling on the evolution of density waves has been studied, a quantitative correlation to connect gap depth with the cooling timescale is lacking. To address this knowledge gap, we employ the grid-based code Athena++ to simulate disk-planet interactions, treating cooling as a thermal relaxation process. We establish quantitative dependencies of steady-state gap depth (Equation 36) and width (Equation 41) on planetary mass, Shakura–Sunyaev viscosity, disk scale height, and thermal relaxation timescale (β). We confirm previous results that gap opening is the weakest when the thermal relaxation timescale is comparable to the local dynamical timescale. Significant variations in gap depth, up to an order of magnitude, are found with different β. In terms of width, a gap is at its narrowest around β = 1, approximately 10%–20% narrower compared to the isothermal case. When β ∼ 100, it can be ∼20% wider, and higher viscosity enhances this effect. We derive possible masses of the gas gap-opening planets in AS 209, HD 163296, MWC 480, and HL Tau, accounting for the uncertainties in the local thermal relaxation timescale.

GRAPES-3 is a mid-altitude (2200 m) and near-equatorial (114N) air shower array, overlapping in its field of view for cosmic-ray observations with experiments that are located in the Northern and Southern Hemispheres. We analyze a sample of 3.7 × 10⁹ cosmic-ray events collected by the GRAPES-3 experiment between 2013 January 1 and 2016 December 31 with a median energy of ∼16 TeV for study of small-scale (<60°) angular-scale anisotropies. We observed two structures, labeled A and B, that deviate from the expected isotropic distribution of cosmic rays in a statistically significant manner. Structure A spans 50°–80° in R.A. and from −15° to 30° in decl. The relative excess observed in structure A is at the level of (6.5 ± 1.3) × 10⁻⁴ with a statistical significance of 6.8 standard deviations. Structure B is observed in the R.A. range 110°–140° and at decl. from −10° to 30°. The relative excess observed in this region is at the level of (4.9 ± 1.4) × 10⁻⁴ with a statistical significance of 4.7 standard deviations. These structures are consistent with those reported by Milagro, ARGO-YBJ, and HAWC. These observations could provide a better understanding of the sources of cosmic rays, their propagation, and the magnetic structures in our Galaxy.

We present a comprehensive study of type III radio bursts and their association with solar flares of magnitude M1.0 and larger, as observed by four widely separated spacecraft (Parker Solar Probe, Solar Orbiter, STEREO-A, and Wind). Our main focus is the introduction and validation of two methods for localizing radio bursts using the available multispacecraft data. The first method utilizes intensity fitting with a circular Gaussian distribution, while the second method is based on the time arrival of radio bursts. We demonstrate the effectiveness of these methods through the analysis of a single type III burst event and compare their results with the traditional radio triangulation technique. Furthermore, we conduct a statistical study of 17 type III bursts associated with M- and X-class solar flares in years 2020–2022. Our findings suggest a possible correlation between solar flare intensities and longitudes, with east limb flares tending to be weaker than west limb flares. We also observe a systematic drift of radio burst longitudes toward the east, potentially explained by a poleward component of the local density gradient. Our results suggest a strong correlation between solar flare intensities and radio burst properties, enhancing our understanding of the relationship between solar flares and type III radio bursts.

We calculate new evolutionary models of rotating primordial very massive stars, with initial mass from 100 M_⊙ to 200 M_⊙, for two values of the initial metallicity Z = 0 and Z = 0.0002. For the first time in this mass range, we consider stellar rotation and pulsation-driven mass loss, along with radiative winds. The models evolve from the zero-age main sequence until the onset of pair-instability. We discuss the main properties of the models during their evolution and then focus on the final fate and the possible progenitors of jet-driven events. All tracks that undergo pulsational-pair instability produce successful gamma-ray bursts (GRB) in the collapsar framework, while those that collapse directly to black holes (BH) produce jet-driven supernova events. In these latter cases, the expected black hole mass changes due to the jet propagation inside the progenitor, resulting in different models that should produce BH within the pair-instability black hole mass gap. Successful GRBs predicted here from zero metallicity, and very metal-poor progenitors, may be bright enough to be detected even up to redshift ∼20 using current telescopes such as the Swift-BAT X-ray detector and the JWST.

Low dust opacity spectral indices (β < 1) measured in the inner envelopes of class 0/I young stellar objects (age ∼10^4–5 yr) have been interpreted as the presence of (sub-)millimeter dust grains in these environments. The density conditions and the lifetimes of collapsing envelopes have proven unfavorable for the growth of solids up to millimeter sizes. As an alternative, magnetohydrodynamical simulations suggest that protostellar jets and outflows might lift grains from circumstellar disks and diffuse them in the envelope. We reframe available data for the CALYPSO sample of Class 0/I sources and show tentative evidence for an anticorrelation between the value of β_{1–3 mm} measured in the inner envelope and the mass-loss rate of their jets and outflows, supporting a connection between the two. We discuss the implications that dust transport from the disk to the inner envelope might have for several aspects of planet formation. Finally, we urge for more accurate measurements of both correlated quantities and the extension of this work to larger samples, necessary to further test the transport scenario.

In the forthcoming era of big astronomical data, it is a burden to find target sources from ground-based and space-based telescopes. Although machine-learning methods have been extensively utilized to address this issue, the incorporation of in-depth data analysis can significantly enhance the efficiency of identifying target sources when dealing with massive volumes of astronomical data. In this work, we focused on the task of finding active galactic nucleus (AGN) candidates and identifying BL Lacertae objects (BL Lac) or flat spectrum radio quasar (FSRQ) candidates from the 4FGL_DR3 uncertain sources. We studied the correlations among the attributes of the 4FGL_DR3 catalog and proposed a novel method, named fractal dimension–inverse discrete wavelet transform (FDIDWT), to transform the original data. The transformed data set is characterized as low-dimensional and feature-highlighted, with the estimation of correlation features by fractal dimension theory and the multi-resolution analysis by inverse discrete wavelet transform (IDWT). Combining the FDIDWT method with an improved lightweight MatchboxConv1D model, we accomplished two missions: (1) to distinguish the AGNs from others (non-AGNs) in the 4FGL_DR3 uncertain sources with an accuracy of 96.65% ± 1.32%, namely Mission A; and (2) to classify blazar candidates of uncertain type into BL Lacs or FSRQs with an accuracy of 92.03% ± 2.2%, namely Mission B. There are 1354 AGN candidates in Mission A, and 482 BL Lacs candidates and 128 FSRQ candidates were found in Mission B. The results show a high consistency of greater than 98% with the results in previous works. In addition, our method has the advantage of finding less variable and relatively faint sources than ordinary methods.

We present the discovery of Ursa Major III/UNIONS 1, the least luminous known satellite of the Milky Way, which is estimated to have an absolute V-band magnitude of $+{2.2}_{-0.3}^{+0.4}$ mag, equivalent to a total stellar mass of ${16}_{-5}^{+6}$M_⊙. Ursa Major III/UNIONS 1 was uncovered in the deep, wide-field Ultraviolet Near Infrared Optical Northern Survey (UNIONS) and is consistent with an old (τ > 11 Gyr), metal-poor ([Fe/H] ∼ −2.2) stellar population at a heliocentric distance of ∼10 kpc. Despite its being compact (r_h = 3 ± 1 pc) and composed of few stars, we confirm the reality of Ursa Major III/UNIONS 1 with Keck II/DEIMOS follow-up spectroscopy and identify 11 radial velocity members, eight of which have full astrometric data from Gaia and are co-moving based on their proper motions. Based on these 11 radial velocity members, we derive an intrinsic velocity dispersion of ${3.7}_{-1.0}^{+1.4}$ km s⁻¹ but some caveats preclude this value from being interpreted as a direct indicator of the underlying gravitational potential at this time. Primarily, the exclusion of the largest velocity outlier from the member list drops the velocity dispersion to ${1.9}_{-1.1}^{+1.4}$ km s⁻¹, and the subsequent removal of an additional outlier star produces an unresolved velocity dispersion. While the presence of binary stars may be inflating the measurement, the possibility of a significant velocity dispersion makes Ursa Major III/UNIONS 1 a high-priority candidate for multi-epoch spectroscopic follow-ups to deduce the true nature of this incredibly faint satellite.

This work studies the relationship between accretion-disk size and quasar properties, using a sample of 95 quasars from the Sloan Digital Sky Survey Reverberation Mapping Project with measured lags between the g and i photometric bands. Our sample includes disk lags that are both longer and shorter than predicted by the Shakura and Sunyaev model, requiring explanations that satisfy both cases. Although our quasars each have one lag measurement, we explore the wavelength-dependent effects of diffuse broad-line region (BLR) contamination through our sample's broad redshift range, 0.1 < z < 1.2. We do not find significant evidence of variable diffuse Fe ii and Balmer nebular emission in the rms spectra, nor from Anderson–Darling tests of quasars in redshift ranges with and without diffuse nebular emission falling in the observed-frame filters. Contrary to previous work, we do not detect a significant correlation between the measured continuum and BLR lags in our luminous quasar sample, similarly suggesting that our continuum lags are not dominated by diffuse nebular emission. Similar to other studies, we find that quasars with larger-than-expected continuum lags have lower 3000 Å luminosities, and we additionally find longer continuum lags with lower X-ray luminosities and black hole masses. Our lack of evidence for diffuse BLR contribution to the lags indicates that the anticorrelation between continuum lag and luminosity is not likely to be due to the Baldwin effect. Instead, these anticorrelations favor models in which the continuum lag increases in lower-luminosity active galactic nuclei, including scenarios featuring magnetic coupling between the accretion disk and X-ray corona, and/or ripples or rims in the disk.

High-velocity clouds (HVCs) are multiphase gas structures whose velocities (∣v_LSR∣ ≥ 100 km s⁻¹) are too high to be explained by Galactic disk rotation. While large HVCs are well characterized, compact and small HVCs (with H i angular sizes of a few degrees) are poorly understood. Possible origins for such small clouds include Milky Way (MW) halo gas or fragments of the Magellanic System, but neither their origin nor their connection to the MW halo has been confirmed. We use new Hubble Space Telescope/Cosmic Origins Spectrograph UV spectra and Green Bank Telescope H i spectra to measure the metallicities of five small HVCs in the southern Galactic sky projected near the Magellanic System. We build a set of distance-dependent Cloudy photoionization models for each cloud and calculate their ionization-corrected metallicities. All five small HVCs have oxygen metallicities ≤0.17 Z_⊙, indicating they do not originate in the disk of the MW. Two of the five have metallicities of 0.16–0.17 Z_⊙, similar to the Magellanic Stream, suggesting these clouds are fragments of the Magellanic System. The remaining three clouds have much lower metallicities of 0.02–0.04 Z_⊙. While the origin of these low-metallicity clouds is unclear, they could be gaseous minihalos or gas stripped from dwarf galaxies by ram pressure or tidal interactions. These results suggest that small HVCs do not all reside in the inner MW halo or the Magellanic System, but instead can trace more distant structures.

We present James Webb Space Telescope imaging from 2 to 21 μm of the edge-on protoplanetary disk around the embedded young star IRAS04302+2247. The structure of the source shows two reflection nebulae separated by a dark lane. The source extent is dominated by the extended filamentary envelope at ∼4.4 μm and shorter wavelengths, transitioning at 7.7 μm and longer wavelengths to more compact lobes of scattered light from the disk itself. The dark lane thickness does not vary significantly with wavelength, which we interpret as an indication for intermediate-sized (∼10 μm) grains in the upper layers of the disk. Intriguingly, we find that the brightest nebula of IRAS40302 switches side between 12.8 and 21 μm. We explore the effect of a tilted inner region on the general appearance of edge-on disks. We find that radiative transfer models of a disk including a tilted inner region can reproduce an inversion in the brightest nebula. In addition, for specific orientations, the model predicts strong lateral asymmetries, which can occur for more than half possible viewing azimuths. A large number of edge-on protoplanetary disks observed in scattered light show such lateral asymmetries (15/20), which suggests that a large fraction of protoplanetary disks might host a tilted inner region. Stellar spots may also induce lateral asymmetries, which are expected to vary over a significantly shorter timescale. Variability studies of edge-on disks would allow us to test the dominant scenario for the origin of these asymmetries.

Solar radio emissions provide several unique diagnostics to estimate different physical parameters of the solar corona, which are otherwise simply inaccessible. However, imaging the highly dynamic solar coronal emissions spanning a large range of angular scales at radio wavelengths is extremely challenging. At gigahertz frequencies, MeerKAT radio telescope is possibly globally the best-suited instrument at present for providing high-fidelity spectroscopic snapshot solar images. Here, we present the first published spectroscopic images of the Sun made using the observations with MeerKAT in the 880–1670 MHz band. This work demonstrates the high fidelity of spectroscopic snapshot MeerKAT solar images through a comparison with simulated radio images at MeerKAT frequencies. The observed images show extremely good morphological similarities with the simulated images. Our analysis shows that below ∼900 MHz MeerKAT images can recover essentially the entire flux density from the large angular-scale solar disk. Not surprisingly, at higher frequencies, the missing flux density can be as large as ∼50%. However, it can potentially be estimated and corrected for. We believe once solar observation with MeerKAT is commissioned, it will enable a host of novel studies, open the door to a large unexplored phase space with significant discovery potential, and also pave the way for solar science with the upcoming Square Kilometre Array-Mid telescope, of which MeerKAT is a precursor.

We present an effective strategy for extensive analysis of eclipsing time variations (ETVs) using modern and sophisticated optimization methods that comprise a set of tools to investigate period variability mechanisms in eclipsing binary stars such as the light-time effect, the Applegate mechanism, and mass transfer. We implement these methods for the first time assuming that the above mechanisms can act simultaneously in the puzzling W UMa–type binary star TZ Bootis by using archival and new TESS data spanning 75 yr and reexamining the up-to-date ETVs. Preliminary analysis of the TESS data revealed for the first time the presence of a second binary in agreement with previous spectroscopic data and astrometric results from Gaia DR3. We consider the most credible scenario for the ETV: two stellar circumbinary companions of minimum masses M₃ = 0.5 M_☉ and M₄ = 0.14 M_☉ in highly eccentric orbits e₃ = 0.70 and e₄ = 0.82 with periods P₃ = 38 yr and P₄ = 20 yr along with a 24 yr magnetic activity of the secondary component and a long-term period increase (dP/dt = 1.2 × 10⁻⁸ days yr⁻¹), interpreted as a conservative mass transfer from the secondary to the primary component at a rate of dM₁/dt = 3.7 × 10⁻⁹ days yr⁻¹. Further spectroscopic observations, analytical modeling of the second pair, and ETV analysis of both pairs are needed to investigate the quadruple nature of the system.

Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we create N-body models that reproduce the ratio of single-to-binary pulsars in Milky Way–like GCs. We focus especially on NGC 6752, a typical core-collapsed cluster with many observed millisecond pulsars. Previous studies suggested that an increased rate of neutron star binary disruption in the densest clusters could explain the overabundance of single pulsars in these systems. Here, we demonstrate that binary disruption is ineffective and instead we propose that two additional dynamical processes play dominant roles: (1) tidal disruption of main-sequence stars by neutron stars and (2) gravitational collapse of heavy white dwarf binary merger remnants. Neutron stars formed through these processes may also be associated with fast radio bursts similar to those observed recently in an extragalactic GC.

We present a search for host galaxy associations for the third set of repeating fast radio burst (FRB) sources discovered by the CHIME/FRB Collaboration. Using the ∼1' CHIME/FRB baseband localizations and probabilistic methods, we identify potential host galaxies of two FRBs, 20200223B and 20190110C at redshifts of 0.06024(2) and 0.12244(6), respectively. We also discuss the properties of a third marginal candidate host galaxy association for FRB 20191106C with a host redshift of 0.10775(1). The three putative host galaxies are all relatively massive, fall on the standard mass–metallicity relationship for nearby galaxies, and show evidence of ongoing star formation. They also all show signatures of being in a transitional regime, falling in the green valley, which is between the bulk of star-forming and quiescent galaxies. The plausible host galaxies identified by our analysis are consistent with the overall population of repeating and nonrepeating FRB hosts while increasing the fraction of massive and bright galaxies. Coupled with these previous host associations, we identify a possible excess of FRB repeaters whose host galaxies have M_u − M_r colors redder than the bulk of star-forming galaxies. Additional precise localizations are required to confirm this trend.

The Herschel Gould Belt Survey showed that stars form in dense filaments in nearby molecular clouds. Recent studies suggest that massive filaments are bound by the slow shocks caused by accretion flows onto the filaments. The slow shocks are known to be unstable to corrugation deformation of the shock front. Corrugation instability could convert the accretion flow's ram pressure into turbulent pressure that influences the width of the filament, which, according to theory, determines the self-gravitational fragmentation scale and core mass. In spite of its importance, the effect of slow-shock instability on star-forming filaments has not been investigated. In addition, the linear dispersion relation obtained from ideal magnetohydrodynamics (MHD) analysis shows that the most unstable wavelength of shock corrugation is infinitesimally small. In the scale of dense filaments, the effect of ambipolar diffusion can suppress the instability at small scales. This study investigates the influence of ambipolar diffusion on the instability of the slow shock. We perform two-dimensional MHD simulations to examine the linear growth of the slow-shock instability, considering the effect of ambipolar diffusion. The results demonstrate that the most unstable scale of slow-shock instability is approximately 5 times the length scale of ambipolar diffusion ℓ_AD calculated using post-shock variables, where ℓ_AD corresponds to the scale where the magnetic Reynolds number for ambipolar diffusivity is unity.

We use the upgraded Giant Metrewave Radio Telescope (uGMRT) to measure scintillation arc properties in six bright canonical pulsars with simultaneous dual-frequency coverage. These observations, at frequencies from 300 to 750 MHz, allowed for detailed analysis of arc evolution across frequency and epoch. We perform more robust determinations of frequency dependence for arc curvature, scintillation bandwidth, and scintillation timescale, and comparison between arc curvature and pseudo-curvature than allowed by single-frequency-band-per-epoch measurements, which we find to agree with theory and previous literature. We find a strong correlation between arc asymmetry and arc curvature, which we have replicated using simulations, and attribute to a bias in the Hough transform approach to scintillation arc analysis. Possible evidence for an approximately week-long timescale over which a given scattering screen dominates signal propagation was found by tracking visible scintillation arcs in each epoch in PSR J1136+1551. The inclusion of a 155-minute observation allowed us to resolve the scale of scintillation variations on short timescales, which we find to be directly tied to the amount of interstellar medium sampled over the observation. Some of our pulsars showed either consistent or emerging asymmetries in arc curvature, indicating instances of refraction across their lines of sight. Significant features in various pulsars, such as multiple scintillation arcs in PSR J1136+1551 and flat arclets in PSR J1509+5531, that have been found in previous works, were also detected. The simultaneous multiple-band observing capability of the upgraded GMRT shows excellent promise for future pulsar scintillation work.

Cosmic variance introduces significant uncertainties into galaxy number density properties when surveying the high-redshift Universe with a small volume. Such uncertainties produce the field-to-field variance σ_g of galaxy numbers in observational astronomy, which significantly affects the luminosity function (LF) measurement of Lyα emitters (LAEs). For most previous Lyα LF studies, σ_g is often adopted from predictions by cosmological simulations, but barely confirmed by observations. Measuring cosmic variance requires a huge sample over a large volume, exceeding the capabilities of most astronomical instruments. In this study, we demonstrate an observational approach for measuring the cosmic variance contribution for z ≈ 2.2 Lyα LFs. The LAE candidates are observed using the narrow band and broad band of the Subaru/Hyper Suprime-Cam in eight independent fields, making the total survey area ≃11.62 deg² and a comoving volume of ≃8.71 × 10⁶ Mpc³. We report a best-fit Schechter function with parameters α = −1.75 (fixed), ${L}_{\mathrm{Ly}\alpha }^{* }={5.95}_{-0.96}^{+1.22}\times {10}^{42}$ erg s⁻¹, and ${\phi }_{\mathrm{Ly}\alpha }^{* }={5.26}_{-1.27}^{+1.65}$ × 10⁻⁴ Mpc⁻³ for the overall Lyα LFs. After clipping out the regions that may bias the cosmic variance measurements, we calculate σ_g by sampling LAEs within multiple pointings on the field image. We investigate the relation between σ_g and survey volume V, and fit a simple power-law ${\sigma }_{g}=k\times {\left(\tfrac{V}{{10}^{5}\,{\mathrm{Mpc}}^{3}}\right)}^{\beta }$. We find best-fit values of $-{1.399}_{-0.156}^{+0.160}$ for β and ${1.249}_{-0.193}^{+0.213}$ for k. We compare our measurements with predictions from simulations and find that the cosmic variance of LAEs is likely larger than that of general star-forming galaxies.

Aligned interstellar grains produce polarized extinction (observed at wavelengths from the far-ultraviolet to the mid-infrared) and polarized thermal emission (observed at far-infrared and submm wavelengths). The grains must be quite nonspherical, but the actual shapes are unknown. The relative efficacy for aligned grains to produce polarization at optical versus infrared wavelengths depends on particle shape. The discrete dipole approximation is used to calculate polarization cross sections for 20 different convex shapes, for wavelengths from 0.1 to 100 μm, and grain sizes a_eff from 0.05 to 0.3 μm. Spheroids, cylinders, square prisms, and triaxial ellipsoids are considered. Minimum aspect ratios required by the observed starlight polarization are determined. Some shapes can also be ruled out because they provide too little or too much polarization at far-infrared and submm wavelengths. The ratio of 10 μm polarization to integrated optical polarization is almost independent of grain shape, varying by only ±8% among the viable convex shapes; thus, at least for convex grains, uncertainties in grain shape cannot account for the discrepancy between predicted and observed 10 μm polarization toward Cyg OB2-12.

Galactic gamma-ray diffuse emission (GDE) is emitted by cosmic rays (CRs), ultra-relativistic protons, and electrons, interacting with gas and electromagnetic radiation fields in the interstellar medium. Here we present the analysis of teraelectronvolt diffuse emission from a region of the Galactic plane over the range in longitude of l ∈ [43°, 73°], using data collected with the High Altitude Water Cherenkov (HAWC) detector. Spectral, longitudinal, and latitudinal distributions of the teraelectronvolt diffuse emission are shown. The radiation spectrum is compatible with the spectrum of the emission arising from a CR population with an index similar to that of the observed CRs. When comparing with the DRAGONbase model, the HAWC GDE flux is higher by about a factor of 2. Unresolved sources such as pulsar wind nebulae and teraelectronvolt halos could explain the excess emission. Finally, deviations of the Galactic CR flux from the locally measured CR flux may additionally explain the difference between the predicted and measured diffuse fluxes.

The presence of supermassive black holes with M ∼ 10⁹M_⊙ hosted by the luminous quasars at cosmological redshift z ≥ 6 is still an open problem in astrophysical cosmology. Here we study the formation of massive black holes at high redshift (z ≫ 7) through Hoyle–Lyttleton–Bondi accretion of self-interacting dark matter (SIDM) onto a 20 M_⊙ seed black hole moving with a velocity ∼100 km s⁻¹ inside the short, mean-free path region of an SIDM halo. We consider observational constraints on a specific SIDM cross section, σ/m_dm = (0.1–5) cm² g⁻¹. Formation timescale of massive black holes with M = (10³–10⁸) M_⊙ is calculated for the universal Navarro–Frenk–White (NFW) profile, singular isothermal sphere (SIS), other power-law profiles with a cusp index 2.19 ≤ γ ≤ 2.5 of accreted dark matter, and modified-core isothermal profiles. The ambient sound speed is taken as C_s = (10–100) km s⁻¹. It is found that an NFW profile with halo concentration C = 4.75–32.58 estimated at z = 20 and 30 for halo masses M₂₀₀ = (10¹²–10¹⁴) M_⊙ favors formation of massive black holes with M = (10³–10⁸) M_⊙ at high redshift, well before quasar epoch. In this profile, these black holes grow within timescales (0.1–69) Myr at z = 16–20. For the SIS profile, the black hole formation timescales are short compared to NFW. For power-law profiles, massive black holes with M = (10⁶–10⁸) M_⊙ grow within a few tens to 100 Myr (z = 5–30). For modified-core isothermal profiles, the timescale of massive black hole formation is in the range (0.79–464.08) Myr (z = 8–30).

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are believed to be associated with magnetars, which have extremely strong magnetic fields. Recently, with the operation of the Imaging X-ray Polarimetry Explorer (IXPE), the polarization information of two AXPs and one SGR have been investigated. In this work, we report the observational results of the fourth magnetar, 1E 2259+586, with IXPE, and perform a joint analysis with observations from Neutron Star Interior Composition Explorer. We find that the emission from 1E 2259+586 is linearly polarized, with a polarization degree (5.3% ± 1.3%) and a polarization angle −77° ± 7° in the 2–8 keV energy range. Additionally, both the polarization degree and polarization angle exhibit variability with the pulse phase, and there is a hint of anticorrelation between the polarization degree and the flux, which is similar to AXP 1RXS J170849.0-400910. The phase-dependent polarization angle displays a sinusoidal profile and can be well fitted with the rotating vector model, indicating that the magnetic dipole field dominated the magnetic structure of the pulsar, and the variation in the polarization angle was modulated by the pulsar's rotation.

The unidentified infrared emission features at 3.3, 6.2, 7.7, 8.6, 11.3, and 12.7 μm are ubiquitously seen in a wide variety of astrophysical regions and commonly attributed to polycyclic aromatic hydrocarbon (PAH) molecules. However, the unambiguous identification of any individual, specific PAH molecules has proven elusive until very recently, when two isomers of cyanonapthalene, which consists of two fused benzene rings and substitutes a nitrile (–CN) group for a hydrogen atom, were discovered in the Taurus Molecular Cloud, based on their rotational transitions at radio frequencies. To facilitate the James Webb Space Telescope (JWST) to search for cyanonapthalenes in astrophysical regions, we model the vibrational excitation of cyanonapthalenes and calculate their infrared emission spectra in a number of representative astrophysical regions. The model emission spectra and intensities will allow JWST to quantitatively determine or place an upper limit on the abundances of cyanonapthalenes.

Protostellar outflows often present a knotty appearance, providing evidence of sporadic accretion in stellar mass growth. To understand the direct relation between mass accretion and ejection, we analyze the contemporaneous accretion activity and associated ejection components in B335. B335 has brightened in the mid-IR by 2.5 mag since 2010, indicating increased luminosity, presumably due to an increased mass accretion rate onto the protostar. Atacama Large Millimeter/submillimeter Array (ALMA) observations of ¹²CO emission in the outflow reveal high-velocity emission, estimated to have been ejected 4.6–2 yr before the ALMA observations and consistent with the jump in mid-IR brightness. The consistency in timing suggests that the detected high-velocity ejection components are directly linked to the most recent accretion activity. We calculated the kinetic energy, momentum, and force for the ejection component associated with the most recent accretion activity and found that, at least, about 1.0% of the accreted mass has been ejected. More accurate information on the jet inclination and the temperature of the ejected gas components will better constrain the ejected mass induced by the recently enhanced accretion event.

We present results from the Very Long Baseline Array multifrequency (1.6, 4.4, 8.6, and 22 GHz), high-sensitivity (∼25 μJy beam⁻¹), subparsec-scale (<1 pc) observations and spectral energy distributions for a sample of 12 local active galactic nuclei (AGNs), a subset from our previous volume-complete sample with hard-X-ray (14–195 keV) luminosities above 10⁴² erg s⁻¹, out to a distance of 40 Mpc. All 12 of the sources presented here were detected in the C (4.4 GHz) and X (8.6 GHz) bands, 75% in the L band (1.6 GHz), and 50% in the K band (22 GHz). Most sources showed compact, resolved/slightly resolved, central subparsec-scale radio morphology, except for a few with extended outflow-like features. A couple of sources have an additional component that may indicate the presence of a dual-core, single or double-sided jet or a more intricate feature, such as radio emission resulting from interaction with the nearby interstellar medium. The spectral slopes are mostly gigahertz-peaked or curved, with a few showing steep, flat, or inverted spectra. We found that at the subparsec scale, the gigahertz-peaked spectra belong to the low-accreting, radio-loud AGNs, with a tendency to produce strong outflows, possibly small-scale jets, and/or have a coronal origin. In contrast, flat/inverted spectra suggest compact radio emission from the central regions of highly accreting AGNs, possibly associated with radio-quiet AGNs producing winds/shocks or nuclear star formation in the vicinity of black holes.

Black hole (BH) ultracompact X-ray binaries (UCXBs) are potential Galactic low-frequency gravitational wave (GW) sources. As an alternative channel, BH UCXBs can evolve from BH+He star binaries. In this work, we perform a detailed stellar evolution model for the formation and evolution of BH UCXBs evolving from the He star channel to diagnose their detectability as low-frequency GW sources. Our calculations found that some nascent BH+He star binaries after the common-envelope (CE) phase could evolve into UCXB-LISA sources with a maximum GW frequency of ∼5 mHz, which can be detected in a distance of 10 kpc (or 100 kpc). Once BH+He star systems become UCXBs through mass transfer, they would emit X-ray luminosities of ∼10³⁸ erg s⁻¹, making them ideal multimessenger objects. If the initial He-star masses are ≥0.7 M_⊙, those systems are likely to experience two Roche lobe overflows, and the X-ray luminosity can reach a maximum of 3.5 × 10³⁹ erg s⁻¹ in the second mass-transfer stage. The initial He-star masses and initial orbital periods of progenitors of Galactic BH UCXB-LISA sources are in the range of 0.32–2.9 M_⊙ and 0.02–0.19 days, respectively. Nearly all BH+He star binaries in the above parameter space can evolve into GW sources whose chirp masses can be accurately measured. Employing a population synthesis simulation, we predict the birthrate and detection number of Galactic BH UCXB-LISA sources evolving from the He star channel are R = 2.2 × 10⁻⁶ yr⁻¹ and 33 for an optimistic CE parameter, respectively.

We suggest a physically motivated model of the uncorrelated background, which can be used to improve the accuracy of helioseismic frequency measurements when the background contributes significantly to the formation of spectral lines of acoustic resonances. The basic assumption of our model is that the correlation length of the convective motions is small compared with the horizontal wavelength R_⊙/ℓ of the observations, where ℓ is the degree of the spherical harmonic Y_ℓm(θ, φ). When applied to solar power spectra at frequencies below acoustic resonances, the model reveals a distinct sensitivity to solar rotation: advection of the convective velocity pattern brings spatial correlations in the apparent stochastic velocity field (temporal correlations in the corotating frame induce spatial correlations in the inertial frame). The induced spatiotemporal correlations manifest themselves as an antisymmetric component in the dependence of the convective noise power on azimuthal order m, which allows us to address the solar differential rotation. With 360 days of data obtained by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, we measure three components of the rotation rate as a function of latitude using only ℓ = 300. This result indicates that the model suggests a new way of measuring solar subsurface rotation. This approach can complement traditional measurements based on correlation tracking.

Extreme scattering events (ESEs) are observed as dramatic (>50%) drops in flux density that occur over an extended period of weeks to months. Discrete plasma lensing structures are theorized to scatter the radio waves produced by distant sources such as pulsars, causing the signature decrease in flux density and characteristic caustic spikes in ESE light curves. While plasma lens models in the extant literature have reproduced key features of ESE light curves, they have all faced the problem of being highly overdense and overpressured relative to the surrounding interstellar medium by orders of magnitude. We model ESEs by numerically ray tracing through analytic, volumetric plasma lens models by solving the eikonal equation. Delaunay triangulation connecting the rays approximates the wave front, generating a mapping from the observer plane to the source plane to account for multiple imaging. This eikonal method of ray tracing is tested against known analytic solutions and is then applied to a three-dimensional Gaussian-distributed electron volume density lens and a filament model inspired by Grafton et al. We find convergence of our numerical results with established analytic solutions, validating our numerical method, and reproduce ESE-like light curves. Our numerical ray-tracing method lends itself well to exploring the lensing effects of volumetric turbulence as well as sheet-like lenses, which is currently in progress.

The orbit distribution of young stars in the Galactic disk is highly structured, from well-defined clusters to streams of stars that may be widely dispersed across the sky, but are compact in orbital action-angle space. The age distribution of such groups can constrain the timescales over which conatal groups of stars disperse into the "field." Gaia data have proven powerful in identifying such groups in action-angle space, but the resulting member samples are often too small and have too narrow a color–magnitude diagram (CMD) coverage to allow robust age determinations. Here, we develop and illustrate a new approach that can estimate robust stellar population ages for such groups of stars. This first entails projecting the predetermined action-angle distribution into the 5D space of positions, parallaxes, and proper motions, where much larger samples of likely members can be identified over a much wider range of the CMD. It then entails isochrone fitting that accounts for: (a) widely varying distances and reddenings; (b) outliers and binaries; (c) sparsely populated main-sequence turnoffs, by incorporating the age information of the low-mass main sequence; and (d) the possible presence of an intrinsic age spread in the stellar population. When we apply this approach to 92 nearby stellar groups identified in 6D orbit space, we find that they are predominantly young (≲1 Gyr), mono-age populations. Many groups are established (known) localized clusters with possible tidal tails, while others tend to be widely dispersed and manifestly unbound. This new age-dating tool offers a stringent approach to understanding on which orbits stars form in the solar neighborhood and how quickly they disperse into the field.

We have carried out a detailed study of individual pulse emission from the pulsar J2022+5154 (B2021+51), observed at 2250 MHz using the Jiamusi 66 m radio telescope. We have investigated the modulations in single-pulse behavior using fluctuation spectral analysis, which shows the presence of two prominent periodicities, around 5 and 40 rotation periods, respectively. The shorter periodicity is associated with the phenomenon of subpulse drifting. In the absence of aliasing, the emission pattern is demonstrated to consist of eight subbeams, which rotate around the magnetic axis in about 45 periods. In addition to subpulse drifting, the pulsar also shows the presence of periodic amplitude modulation with a longer periodicity in the single-pulse sequence. The pulsar joins a select group that shows the presence of periodic phase-modulated drifting as well as amplitude-modulated drifting. This provides further evidence for the two phenomena being distinct from each other with different physical origins.

A critically important process affecting the climate evolution and potential habitability of an exoplanet is atmospheric escape, in which high-energy radiation from a star drives the escape of hydrogen atoms and other light elements from a planet's atmosphere. L 98-59 is a benchmark system for studying such atmospheric processes, with three transiting terrestrial-sized planets receiving Venus-like instellations (4–25 S_⊕) from their M3 host star. We use the VPLanet model to simulate the evolution of the L 98-59 system and the atmospheric escape of its inner three small planets, given different assumed initial water quantities. We find that, regardless of their initial water content, all three planets accumulate significant quantities of oxygen due to efficient water photolysis and hydrogen loss. All three planets also receive enough strong X-ray and extreme-ultraviolet flux to drive rapid water loss, which considerably affects their developing climates and atmospheres. Even in scenarios of low initial water content, our results suggest that the JWST will be sensitive to observations of retained oxygen on the L 98-59 planets in its future scheduled observations, with planets b and c being the most likely targets to possess an extended atmosphere. Our results constrain the atmospheric evolution of these small rocky planets, and they provide context for current and future observations of the L 98-59 system to generalize our understanding of multiterrestrial planet systems.

While general relativity predicts only two tensor modes for gravitational-wave (GW) polarization, general metric theories of gravity allow for up to four additional modes, including two vector and two scalar modes. Observing the polarization modes of GWs could provide a direct test of the modified gravity. The stochastic GW background (SGWB), which can be detected by space-based laser-interferometric detectors at design sensitivity, will provide an opportunity to directly measure alternative polarization. In this paper, we investigate the performance of the LISA-TianQin network for detecting alternative polarizations of stochastic backgrounds, and propose a method to separate different polarization modes. First, we generalize the small antenna approximation to compute the overlap reduction functions for the SGWB with arbitrary polarization, which is suitable for any time-delay interferometry combination. Then we analyze the detection capability of LISA-TianQin for the SGWB with different polarizations. Based on the orbital characteristics of LISA-TianQin, we propose a method to distinguish different polarization modes from their mixed data. Finally, simulation tests are performed to verify the effectiveness of the method. The results of the simulations demonstrate that LISA-TianQin, when employing our proposed method, has the ability to differentiate between various polarization modes, with a specific emphasis on the ability to distinguish between the breathing and longitudinal modes.

We present JCMT POL-2 850 μm dust polarization observations and Mimir H-band stellar polarization observations toward the starless core L 1512. We detect the highly ordered core-scale magnetic field traced by the POL-2 data, of which the field orientation is consistent with the parsec-scale magnetic fields traced by Planck data, suggesting the large-scale fields thread from the low-density region to the dense core region in this cloud. The surrounding magnetic field traced by the Mimir data shows a wider variation in the field orientation, suggesting there could be a transition of magnetic field morphology at the envelope-scale. L 1512 was suggested to be presumably older than 1.4 Myr in a previous study via time-dependent chemical analysis, hinting that the magnetic field could be strong enough to slow the collapse of L 1512. In this study, we use the Davis–Chandrasekhar–Fermi method to derive a plane-of-sky magnetic field strength (B_pos) of 18 ± 7 μG and an observed mass-to-flux ratio (λ_obs) of 3.5 ± 2.4, suggesting that L 1512 is magnetically supercritical. However, the absence of significant infall motion and the presence of an oscillating envelope are inconsistent with the magnetically supercritical condition. Using a virial analysis, we suggest the presence of a hitherto hidden line-of-sight magnetic field strength of ∼27 μG with a mass-to-flux ratio (λ_tot) of ∼1.6, in which case both magnetic and kinetic pressures are important in supporting the L 1512 core. On the other hand, L 1512 may have just reached supercriticality and will collapse at any time.

We present individual star formation histories (SFHs) of ∼3000 massive galaxies (log(M_*/M_⊙) > 10.5) from the Large Early Galaxy Astrophysics Census spectroscopic survey at a lookback time of ∼7 billion yr and quantify the population trends leveraging 20 hr deep-integrated spectra of these ∼1800 star-forming and ∼1200 quiescent galaxies at 0.6 < z < 1.0. Essentially all galaxies at this epoch contain stars of age <3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectrophotometric analysis using parametric and nonparametric Bayesian stellar population synthesis modeling tools—Bagpipes and Prospector—to constrain the median SFHs of this mass complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified as t₅₀ and t₉₀, corresponding to the age of the Universe when galaxies formed 50% and 90% of their total stellar mass, although the two methods disagree at the earliest formation times (e.g., t₁₀). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local Universe—low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star formation (τ_SF,late = t₉₀–t₅₀) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectrophotometric studies of the more distant Universe.

Kilonovae are likely a key site of heavy r-process element production in the Universe, and their optical/infrared spectra contain insights into both the properties of the ejecta and the conditions of the r-process. However, the event GW170817/AT2017gfo is the only kilonova so far with well-observed spectra. To understand the diversity of absorption features that might be observed in future kilonovae spectra, we use the TARDIS Monte Carlo radiative transfer code to simulate a suite of optical spectra spanning a wide range of kilonova ejecta properties and r-process abundance patterns. To identify the most common and prominent absorption lines, we perform dimensionality reduction using an autoencoder, and we find spectra clusters in the latent space representation using a Bayesian Gaussian Mixture model. Our synthetic kilonovae spectra commonly display strong absorption by strontium ₃₈Sr ii, yttrium ₃₈Y ii, and zirconium ₄₀Zr i–ii, with strong lanthanide contributions at low electron fractions (Y_e ≲ 0.25). When a new kilonova is observed, our machine-learning framework will provide context on the dominant absorption lines and key ejecta properties, helping to determine where this event falls within the larger "zoo" of kilonovae spectra.

We have observed the z = 4.3 protocluster SPT2349−56 with the Australia Telescope Compact Array (ATCA) with the aim of detecting radio-loud active galactic nuclei (AGNs) among the ∼30 submillimeter (submm) galaxies (SMGs) identified in the structure. We detect the central complex of submm sources at 2.2 GHz with a luminosity of L_2.2 = (4.42 ± 0.56) × 10²⁵ W Hz⁻¹. MeerKAT and the Australian Square Kilometre Array Pathfinder also detect the source at 816 MHz and 888 MHz, respectively, constraining the radio spectral index to α = −1.45 ± 0.16, implying L_1.4,rest = (2.2 ± 0.2) × 10²⁶ W Hz⁻¹. The radio observations do not have sufficient spatial resolution to uniquely identify one of the three Atacama Large Millimeter/submillimeter Array (ALMA) galaxies as the AGN, however the ALMA source properties themselves suggest a likely host. This radio luminosity is ∼100× higher than expected from star formation, assuming the usual far-infrared–radio correlation, indicating an AGN driven by a forming brightest cluster galaxy. None of the SMGs in SPT2349−56 show signs of AGNs in any other diagnostics available to us, highlighting the radio continuum as a powerful probe of obscured AGNs. We compare these results to field samples of radio sources and SMGs, along with the 22 gravitationally lensed SPT-SMGs also observed in the ATCA program, as well as powerful radio galaxies at high redshifts. The (3.3 ± 0.7) × 10³⁸ W of power from the radio-loud AGN sustained over 100 Myr is comparable to the binding energy of the gas mass of the central halo, and similar to the instantaneous energy injection from supernova feedback from the SMGs in the core region. The SPT2349−56 radio-loud AGNs may be providing strong feedback on a nascent intracluster medium.

We present a uniform forward-modeling analysis of 90 late-M and L dwarfs in nearby young (∼10–200 Myr) moving groups, the Pleiades, and the Hyades using low-resolution (R ≈ 150) near-infrared (0.9–2.4 μm) spectra and the BT-Settl model atmospheres. We derive the objects' effective temperatures, surface gravities, radii, and masses by comparing our spectra to the models using a Bayesian framework with nested sampling and calculate the same parameters using evolutionary models. Assuming the evolutionary-based parameters are more robust, our spectroscopically inferred parameters from BT-Settl exhibit two types of systematic behavior for objects near the M-L spectral type boundary. Several objects are clustered around T_eff ≈ 1800 K and $\mathrm{log}g\approx 5.5$ dex, implying impossibly large masses (150–1400 M_Jup), while others are clustered around T_eff ≳ 3000 K and $\mathrm{log}g\lesssim 3.0$ dex, implying unphysically low masses and unreasonably young ages. We find the fitted BT-Settl model spectra tend to overpredict the peak J- and H-band flux for objects located near the M-L boundary, suggesting the dust content included in the model atmospheres is insufficient to match the observations. By adding an interstellar medium–like reddening law to the BT-Settl model spectra, we find the fits between models and observed spectra are greatly improved, with the largest reddening coefficients occurring at the M-L transition. This work delivers a systematic examination of the BT-Settl model atmospheres and constitutes the largest spectral analysis of benchmark late-M- and L-type brown dwarfs to date.

Emission in forbidden lines of oxygen, neon, and other species are commonly used to trace winds from protoplanetary disks. Using Cloudy, we calculate such emission for parametrized wind models of the magnetothermal type, following Bai et al. These models share characteristics with both photoevaporative and magnetocentrifugal winds, which can be regarded as end members, and are favored by recent theoretical research. Both broad and narrow low-velocity components of the lines can be produced with plausible wind parameters, something that traditional wind models have difficulty with. Line luminosities, blueshifts, and widths, as well as trends of these with accretion luminosity and disk inclination, are in general accordance with observations.

The internal structure of the prestellar core G208.68-19.02-N2 (G208-N2) in the Orion Molecular Cloud 3 (OMC-3) region has been studied with the Atacama Large Millimeter/submillimeter Array. The dust continuum emission revealed a filamentary structure with a length of ∼5000 au and an average H₂ volume density of ∼6 × 10⁷ cm⁻³. At the tip of this filamentary structure, there is a compact object, which we call a nucleus, with a radius of ∼150–200 au and a mass of ∼0.1 M_⊙. The nucleus has a central density of ∼2 × 10⁹ cm⁻³ with a radial density profile of r^{−1.87±0.11}. The density scaling of the nucleus is ∼3.7 times higher than that of the singular isothermal sphere (SIS). This as well as the very low virial parameter of 0.39 suggests that the gravity is dominant over the pressure everywhere in the nucleus. However, there is no sign of CO outflow localized to this nucleus. The filamentary structure is traced by the N₂D⁺ 3–2 emission, but not by the C¹⁸O 2–1 emission, implying the significant CO depletion due to high density and cold temperature. Toward the nucleus, the N₂D⁺ also shows the signature of depletion. This could imply either the depletion of the parent molecule, N₂, or the presence of the embedded very-low luminosity central source that could sublimate the CO in the very small area. The nucleus in G208-N2 is considered to be a prestellar core on the verge of first hydrostatic core (FHSC) formation or a candidate for the FHSC.

We study the magnetic field structures in six giant filaments associated with the spiral arms of the Milky Way by applying the velocity gradient technique (VGT) to the ¹³CO spectroscopic data from the GRS, FUGIN, and SEDIGSM surveys. Unlike dust-polarized emission, the VGT allows us to separate the foreground and background using the velocity information, from which the orientation of the magnetic field can be reliably determined. We find that in most cases the magnetic fields stay aligned with the filament bodies, which are parallel to the disk midplane. Among these, G29, G47, and G51 exhibit smooth magnetic fields, and G24, G339, and G349 exhibit discontinuities. The fact that most filaments have magnetic fields that stay aligned with the Galactic disk midplane suggests that Galactic shear may be responsible for shaping the filaments. The fact that the magnetic field can stay regular at the resolution of our analysis (≲10 pc), where the turbulence crossing time is short compared to the shear time, suggests that turbulent motion cannot effectively disrupt the regular orientation of the magnetic field. The discontinuities found in some filaments can be caused by processes including filament reassembly, gravitational collapse, and stellar feedback.

As the plasma boundary between two distinct plasma populations, dipolarization fronts (DFs) host abundant kinetic-scale substructures that change their normal directions and thus cause their deformation. However, studies on such deformation caused by an electron vortex have been lacking. Here, we present novel observations of a subion-scale magnetic hump (MHu) associated with an oblique electron vortex at a DF through strengthening three components of the magnetic field. A radial electric field in the MHu, showing bipolar variation, is also associated with the electron vortex as it is mainly ascribed to the electron convection term. There is apparent energy conversion ($\overrightarrow{J}\cdot \overrightarrow{E}$∼−0.3 nw m⁻³) from the particles to the electromagnetic field in the MHu's leading part, which is accompanied by inflow and outflow of electromagnetic energy (nonzero ${\rm{\nabla }}\cdot \overrightarrow{S}$). The other regions of the DF host opposite energy conversion ($\overrightarrow{J}\cdot \overrightarrow{E}$ > 0). Broadband parallel electrostatic waves are also observed in the MHu. Our study provides insights into the kinetic-scale processes at DFs.

We report results from a systematic wide-area search for faint dwarf galaxies at heliocentric distances from 0.3 to 2 Mpc using the full 6 yr of data from the Dark Energy Survey (DES). Unlike previous searches over the DES data, this search specifically targeted a field population of faint galaxies located beyond the Milky Way virial radius. We derive our detection efficiency for faint, resolved dwarf galaxies in the Local Volume with a set of synthetic galaxies and expect our search to be complete to M_V ∼ (−7, −10) mag for galaxies at D = (0.3, 2.0) Mpc. We find no new field dwarfs in the DES footprint, but we report the discovery of one high-significance candidate dwarf galaxy at a distance of ${2.2}_{-0.12}^{+0.05}\,\mathrm{Mpc}$, a potential satellite of the Local Volume galaxy NGC 55, separated by 47' (physical separation as small as 30 kpc). We estimate this dwarf galaxy to have an absolute V-band magnitude of $-{8.0}_{-0.3}^{+0.5}\,\mathrm{mag}$ and an azimuthally averaged physical half-light radius of ${2.2}_{-0.4}^{+0.5}\,\mathrm{kpc}$, making this one of the lowest surface brightness galaxies ever found with $\mu =32.3\,\mathrm{mag}\,{\mathrm{arcsec}}^{-2}$. This is the largest, most diffuse galaxy known at this luminosity, suggesting possible tidal interactions with its host.

Evidence from different probes of the stellar initial mass function (IMF) of massive early-type galaxies (ETGs) has repeatedly converged on IMFs more bottom heavy than in the Milky Way (MW). This consensus has come under scrutiny due to often contradictory results from different methods on the level of individual galaxies. In particular, a number of strong lensing probes are ostensibly incompatible with a non-MW IMF. Radial gradients of the IMF—related to gradients of the stellar mass-to-light ratio ϒ—can potentially resolve this issue. We construct Schwarzschild models allowing for ϒ-gradients in seven massive ETGs with MUSE and SINFONI observations. We find dynamical evidence that ϒ increases toward the center for all ETGs. The gradients are confined to subkiloparsec scales. Our results suggest that constant-ϒ models may overestimate the stellar mass of galaxies by up to a factor of 1.5. For all except one galaxy, we find a radius where the total dynamical mass has a minimum. This minimum places the strongest constraints on the IMF outside the center and appears at roughly 1 kpc. We consider the IMF at this radius characteristic for the main body of each ETG. In terms of the IMF mass-normalization α relative to a Kroupa IMF, we find on average an MW-like IMF 〈α_main〉 = 1.03 ± 0.19. In the centers, we find concentrated regions with increased mass normalizations that are less extreme than previous studies suggested, but still point to a Salpeter-like IMF, 〈α_cen〉 = 1.54 ± 0.15.

We study properties of intensity disturbances along polar coronal rays that are associated with plumes below. For this, we draw azimuth–time images of extreme ultraviolet (EUV) emission of 171 Å band observed by the SDO/AIA and white light (WL) observed by the SOHO/LASCO C2 in 2020 July. From the azimuth–time image, we define two tracks in which the EUV intensities were recurrently enhanced during two weeks. The two EUV tracks are rooted at 788 and 814 latitudes, but their projected azimuth angles are changed with time as the Sun rotates. Coherent WL tracks at different altitudes are determined by scaling the azimuth angles of the EUV tracks, accounting for the effect of inclination of coronal rays. From this, we construct time–distance images of WL intensities along WL tracks, whose projected azimuth angle changes along time and altitude, but the intensities are correlated with the EUV intensities measured below. The time–distance images of WL show repeated and inclined intensity features. The propagation speeds in the altitude range 2.3–6 solar radii are calculated to be 159 ± 8 km s⁻¹ and 300 ± 24 km s⁻¹. The EUV and WL intensities are found to be coherent at 1–2 day periods. It is also found that dynamic burst events along the EUV track are responsible for the enhanced emission. We conclude that the variation of the WL intensity along the polar coronal rays is related with the evolution of the EUV intensity below.

Reproduction experiments of radial pyroxene (RP) chondrules were carried out using an Ar–H₂ or Ar gas-jet levitation system in a reducing atmosphere in order to simulate chondrule formation in the protoplanetary disk. The experiments reproduced RP-chondrule textures, consisting of sets of thin pyroxene crystals and mesostasis glass between crystals. However, iron partition coefficients between pyroxene and glassy mesostasis (D_Fe =Fe mol%_pyroxene/Fe mol%_mesostasis) in natural RP chondrules were much higher than that in the experimentally reproduced RP chondrules. The high D_Fe in natural RP chondrules suggests that iron was removed from the mesostasis melt at high temperatures after the growth of pyroxene crystals. We found that many small iron-metal inclusions had formed in the mesostasis glass, indicating that FeO in the high-temperature melt of mesostasis was reduced to metallic iron, and iron in the mesostasis diffused into the newly formed metal inclusions. The formation of the iron-metal inclusions in the mesostasis was reproduced by our experiments in a reducing atmosphere, confirming that D_Fe in natural RP chondrules increased after the growth of RP crystals. Therefore, the D_Fe of RP chondrules can be an indicator to constrain cooling rates and redox states during chondrule formation.

The statistical characteristics of stellar flares at optical bands have received extensive study, but they remain to be studied at soft X-ray bands, in particular for solar-type stars. Here, we present a statistical study of soft X-ray flares on solar-type stars, which can help in understanding the multiwavelength behaviors of stellar flares. We mainly use Chandra Source Catalog Release 2.0, which includes a number of flaring stars with denoted variability, and Gaia Data Release 3, which includes the necessary information for classifying stars. We also develop a set of methods for identifying and classifying stellar soft X-ray flares and estimating their properties. A detailed statistical investigation of 129 flare samples on 103 nearby solar-type stars as selected yields the following main results. (1) The flare energy emitted at the soft X-ray band in our sample ranges from ∼10³³ to ∼10³⁷ erg, and the majority of them are superflares, with the most energetic one having energy of ${6.0}_{-4.7}^{+3.2}\times {10}^{37}\ \mathrm{erg}$. (2) The flare duration is related to its energy as formulated by ${T}_{\mathrm{duration},\mathrm{SXR}}\propto {E}_{\mathrm{flare},\mathrm{SXR}}^{\ 0.201\pm 0.024}$, which is different from those derived at optical and near-IR bands, indicating distinct radiation mechanisms at different bands. (3) The frequency distribution of stellar flares as a function of energy is formulated as ${{dN}}_{\mathrm{flare}}/{{dE}}_{\mathrm{flare},\mathrm{SXR}}\propto {E}_{\mathrm{flare},\mathrm{SXR}}^{\ -1.77}$, which is similar to the results found at other bands and on other types of stars, indicating that the energy emitted at the soft X-ray band could be a constant fraction of the full-band bolometric energy.

Ultralight scalar fields and their noninteracting class, i.e., the so-called fuzzy dark matter (FDM), are dark matter candidates introduced to solve the small-scale problems of the standard cold dark matter. In this paper, we investigate whether the physics of FDM, particularly the quantum pressure that leads to the suppression of structure formation on small scales, could leave significant imprints on the large-scale statistics of matter fluctuations. For this purpose, we utilize the Effective Field Theory of Large Scale Structures (EFT of LSS), wherein small-scale physics is integrated and represented on large scales by only a set of free parameters. These parameters can be determined by fitting them into the cosmological simulations. By fitting the EFT predictions to the simulation data, we determine the value of the speed of sound as a quantitative measure of how UV physics affects large-scale perturbation. We use the Gadget-2 code to study the evolution of 512³ particles in a box with a side length 250 h⁻¹ Mpc. We exploit the suppressed FDM initial power for the FDM universe and perform N-body simulation sufficient to produce accurate—enough for our purpose—results on large scales. In particular, we perform three FDM simulations with different masses and compare their sound speed with the standard cold dark matter (CDM) simulation. We found no difference between the FDM and CDM sound speeds beyond the confidence intervals. However, a consistently increasing trend can be seen in the sound speed for lower masses. This result suggests further investigations using higher-resolution simulations.

The J-region Asymptotic Giant Branch (JAGB) method is a standard candle that leverages the constant luminosities of color-selected, carbon-rich AGB stars, measured in the near-infrared at 1.2 μm. The Chicago-Carnegie Hubble Program has obtained JWST imaging of the SN Ia host galaxies NGC 7250, NGC 4536, and NGC 3972. With these observations, the JAGB method can be studied for the first time using JWST. Lee et al. demonstrated the JAGB magnitude is optimally measured in the outer disks of galaxies, because in the inner regions the JAGB magnitude can vary significantly due to a confluence of reddening, blending, and crowding effects. However, determining where the "outer disk" lies can be subjective. Therefore, we introduce a novel method for systematically selecting the outer disk. In a given galaxy, the JAGB magnitude is first separately measured in concentric regions, and the "outer disk" is then defined as the first radial bin where the JAGB magnitude stabilizes to a few hundredths of a magnitude. After successfully employing this method in our JWST galaxy sample, we find the JAGB stars are well segregated from other stellar populations in color–magnitude space, and have observed dispersions about their individual F115W modes of σ_N7250 = 0.32 mag, σ_N4536 = 0.34 mag, and σ_N3972 = 0.35 mag. These measured dispersions are similar to the scatter measured for the JAGB stars in the LMC using 2MASS data (σ = 0.33 mag). In conclusion, the JAGB stars as observed with JWST clearly demonstrate their considerable power both as high-precision extragalactic distance indicators and as SN Ia supernova calibrators.

We present radio and X-ray studies of A3444 and MS1455.0+2232, two galaxy clusters with radio minihalos in their cool cores. A3444 is imaged using the Giant Metrewave Radio Telescope (GMRT) at 333, 607, and 1300 MHz and the Very Large Array at 1435 MHz. Most of the minihalo is contained within r < 120 kpc, but a fainter extension, stretching out to 380 kpc southwest of the center, is detected at 607 MHz. Using Chandra, we detect four X-ray sloshing cold fronts: three in the cool core at r = 60, 120, and 230 kpc, and a fourth one at r = 400 kpc—in the region of the southwestern radio extension—suggesting that the intracluster medium (ICM) is sloshing on a cluster-wide scale. The radio emission is contained within the envelope defined by these fronts. We also analyzed archival 383 MHz GMRT and Chandra observations of MS 1455.0+2232, which exhibits a known minihalo with its bright part delineated by cold fronts inside the cool core, but with a faint extension beyond the core. Similarly to A3444, we find a cold front at r ∼ 425 kpc, containing the radio emission. Thus the entire diffuse radio emission seen in these clusters appears to be related to large-scale sloshing of the ICM. The radio spectrum of the A3444 minihalo is a power law with a steep index α = 1.0 ± 0.1. The spectrum steepens with increasing distance from the center, as expected if the minihalo originates from reacceleration of relativistic particles by the sloshing-induced turbulence in the ICM.

Recent observations of galaxy clusters and groups with misalignments between their central active galactic nucleus jets and X-ray cavities, or with multiple misaligned cavities, have raised concerns about the jet–bubble connection in cooling cores, and the processes responsible for jet realignment. To investigate the frequency and causes of such misalignments, we construct a sample of 16 cool core galaxy clusters and groups. Using Very Long Baseline Array radio data, we measure the parsec-scale position angle of the jets, and compare it with the position angle of the X-ray cavities detected in Chandra data. Using the overall sample and selected subsets, we consistently find that there is a 30%–38% chance to find a misalignment larger than ΔΨ = 45° when observing a cluster/group with a detected jet and at least one cavity. We determine that projection may account for an apparently large ΔΨ only in a fraction of objects (∼35%), and given that gas dynamical disturbances (such as sloshing) are found in both aligned and misaligned systems, we exclude environmental perturbation as the main driver of cavity–jet misalignment. Moreover, we find that large misalignments (up to ∼90°) are favored over smaller ones (45° ≤ ΔΨ ≤ 70°), and that the change in jet direction can occur on timescales between one and a few tens of Myr. We conclude that misalignments are more likely related to actual reorientation of the jet axis, and we discuss several engine-based mechanisms that may cause these dramatic changes.

Interplanetary coronal mass ejections (ICMEs) are defined as "coherent" if they are capable of responding to external perturbations in a collective manner. This implies that information must be able to propagate across ICME structures, and if this is not the case, single-point in situ measurements cannot be considered as indicative of global ICME properties. Here, we investigate the role of Alfvénic fluctuations (AFs) as mediators of ICME coherence. We consider multipoint magnetic field and plasma measurements of 10 ICMEs observed by the ACE and Wind spacecraft at 1 au at longitudinal separations of 0.5°–0.7°. For each event, we analyze the Alfvénicity in terms of the residual energy and cross helicity of fluctuations, and the coherence in terms of the magnetic correlation between Wind and ACE. We find that ∼65% and 90% of ICME sheaths and magnetic ejecta (MEs), respectively, present extended AFs covering at least 20% of the structure. Cross helicity suggests AFs of solar and interplanetary origin may coexist in the ICME population at 1 au. AFs are mainly concentrated downstream of shocks and in the back of MEs. The magnetic field is poorly correlated within sheaths, while the correlation decreases from the front to the back of the MEs for most magnetic field components. AFs are also associated with lower magnetic field correlations. This suggests either that ICME coherence is not mediated by Alfvén waves, implying that the coherence scale may be smaller than previously predicted, or that the magnetic field correlation is not a measure of coherence.

The Parker Solar Probe (PSP) provides us with an unprecedentedly close approach to the observation of the Sun and hence the possibility of directly understanding the elementary process that occurs on the kinetic scale of particles' collective interaction in solar coronal plasmas. We report a type of weak solar radio burst (SRB) that was detected by PSP when it passed a low-density magnetic channel during its second encounter phase. These weak SRBs have a low starting frequency of ∼20 MHz and a narrow frequency range from a few tens of MHz to a few hundred kHz. Their dynamic spectra display a strongly evolving feature of the intermediate relative drift rate decreasing rapidly from above 0.01 s⁻¹ to below 0.01 s⁻¹. Analyses based on common empirical models of solar coronal plasmas indicate that these weak SRBs originate from a heliocentric distance of ∼1.1–6.1 R_S (the solar radius), a typical solar wind acceleration region with a low-β plasma, and that their sources have a typical motion velocity of ∼v_A (Alfvén velocity) obviously lower than that of the fast electrons required to effectively excite SRBs. We propose that solitary kinetic Alfvén waves with kinetic scales could be responsible for the generation of these small-scale weak SRBs, called solitary wave radiation.

Some quasi-thermal (QT)-dominated gamma-ray bursts (GRBs) could be well described by a multicolor blackbody (mBB) function or a combined model of a BB plus a nonthermal (NT) component. In this analysis, two QT radiation-dominated bursts with known emission properties (GRB 210610B, likely from a hybrid jet, and GRB 210121A, with a spectrum consistent with nondissipative photospheric emission from a pure hot fireball) are used to make a comparison between these two models. To diagnose the magnetization properties of the central engine, the "top-down" approach proposed by Gao and Zhang is adopted. It is found that diagnoses based on these two models provide similar conclusions qualitatively; however, the model with mBB (or mBB+NT) may give more reasonable physical explanations. This implies that impacts from the GRB jet structure and geometrical broadening on the observed spectrum should be considered. However, conservatively, these methods may be not sensitive enough to distinguish between a pure hot fireball and a mildly magnetized hybrid jet. Some other information is necessary to provide more evidence when determining the jet properties for similar GRBs. Based on these considerations, we suggest that the photospheric emission of GRB 221022B is from a hot jet, where dissipation is caused by an internal shock mechanism due to the increasing Lorentz factor with time, which makes its prompt emission display typical evolution from thermal to NT.

We use time-domain simulations of Jupiter observations to test and develop a beam reconstruction pipeline for the Simons Observatory Small Aperture Telescopes. The method relies on a mapmaker that estimates and subtracts correlated atmospheric noise and a beam fitting code designed to compensate for the bias caused by the mapmaker. We test our reconstruction performance for four different frequency bands against various algorithmic parameters, atmospheric conditions, and input beams. We additionally show the reconstruction quality as a function of the number of available observations and investigate how different calibration strategies affect the beam uncertainty. For all of the cases considered, we find good agreement between the fitted results and the input beam model within an ∼1.5% error for a multipole range ℓ = 30–700 and an ∼0.5% error for a multipole range ℓ = 50–200. We conclude by using a harmonic-domain component separation algorithm to verify that the beam reconstruction errors and biases observed in our analysis do not significantly bias the Simons Observatory r-measurement

Despite their shared origin, members of globular clusters display star-to-star variations in composition. The observed pattern of element abundances is unique to these stellar environments and cannot be fully explained by any proposed mechanism. It remains unclear whether stars form with chemical heterogeneity or inherit it from interactions with other members. These scenarios may be differentiated by the dependence of chemical spread on stellar mass; however, obtaining a sufficiently large mass baseline requires abundance measurements on the lower main sequence, which is too faint for spectroscopy even in the nearest globular clusters. We developed a stellar modeling method to obtain precise chemical abundances for stars near the end of the main sequence from multiband photometry, and we applied it to the globular cluster 47 Tucanae. The computational efficiency is attained by matching chemical elements to the model components that are most sensitive to their abundance. We determined [O/Fe] for ∼5000 members below the main-sequence knee at the level of accuracy, comparable to the spectroscopic measurements of evolved members in the literature. The inferred distribution disfavors stellar interactions as the origin of chemical spread; however, an accurate theory of accretion is required to draw a more definitive conclusion. We anticipate that future observations of 47 Tucanae with the James Webb Space Telescope will extend the mass baseline of our analysis into the substellar regime. Therefore, we present predicted color–magnitude diagrams and mass–magnitude relations for the brown dwarf members of 47 Tucanae.

Magnetic flux ropes are a bundle of twisted magnetic field lines produced by internal electric currents, which are responsible for solar eruptions and are the major drivers of geomagnetic storms. As such, it is crucial to develop a numerical model that can capture the entire evolution of a flux rope, from its birth to death, in order to predict whether adverse space weather events might occur or not. In this paper, we develop a data-driven modeling that combines a time-dependent magnetofrictional approach with a thermodynamic magnetohydrodynamic model. Our numerical modeling successfully reproduces the formation and confined eruption of an observed flux rope, and unveils the physical details behind the observations. Regarding the long-term evolution of the active region, our simulation results indicate that the flux cancellation due to collisional shearing plays a critical role in the formation of the flux rope, corresponding to a substantial increase in magnetic free energy and helicity. Regarding the eruption stage, the deformation of the flux rope during its eruption can cause an increase in the downward tension force, which suppresses it from further rising. This finding may shed light on why some torus-unstable flux ropes lead to failed eruptions after large-angle rotations. Moreover, we find that twisted fluxes can accumulate during confined eruptions, which would breed the subsequent eruptive flares.

We present a comprehensive study of the kinematic properties of the different Galactic disk populations, as defined by the chemical abundance ratios and stellar ages, across a large disk volume (4.5 ≤ R ≤ 15.0 kpc and ∣Z∣ ≤ 3.0 kpc), by using the LAMOST-Gaia red clump sample stars. We determine the median velocities for various spatial and population bins, finding large-scale bulk motions; for example, the wave-like behavior in radial velocity, the north–south discrepancy in azimuthal velocity and the warp signal in vertical velocity, and the amplitudes and spatial dependences of these bulk motions show significant variations for different mono-age and mono-abundance populations. The global spatial behaviors of the velocity dispersions clearly show a signal of spiral arms and a signal of the disk perturbation event within 4 Gyr, as well as disk flaring in the outer region (i.e., R ≥ 12 kpc), mostly for young or alpha-poor stellar populations. Our detailed measurements of age/[α/Fe]-velocity dispersion relations for different disk volumes indicate that young/α-poor populations are likely to originate from dynamic heating by both giant molecular clouds and spiral arms, while old/α-enhanced populations require an obvious contribution from other heating mechanisms, such as merger and accretion, or are born in the chaotic mergers of gas-rich systems and/or turbulent interstellar medium.

Ion-scale wave events or wave storms in the solar wind are characterized by enhancements in magnetic field fluctuations as well as coherent magnetic field polarization signatures at or around the local ion cyclotron frequencies. In this paper, we study in detail one such wave event from Parker Solar Probe's (PSP) fourth encounter, consisting of an initial period of left-handed (LH) polarization abruptly transitioning to a strong period of right-handed (RH) polarization, accompanied by a clear core beam structure in both the alpha and proton velocity distribution functions. A linear stability analysis shows that the LH-polarized waves are anti-sunward propagating Alfvén/ion cyclotron waves primarily driven by a proton cyclotron instability in the proton core population, and the RH polarized waves are anti-sunward propagating fast magnetosonic/whistler waves driven by a firehose-like instability in the secondary alpha beam population. The abrupt transition from LH to RH is caused by a drop in the proton core temperature anisotropy. We find very good agreement between the frequencies and polarizations of the unstable wave modes as predicted by linear theory and those observed in the magnetic field spectra. Given the ubiquity of ion-scale wave signatures observed by PSP, this work gives insight into which exact instabilities may be active and mediating energy transfer in wave–particle interactions in the inner heliosphere, as well as highlighting the role a secondary alpha population may play as a rarely considered source of free energy available for producing wave activity.