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

Emerging dimming occurs in isolated solar active regions (ARs) during the early stages of magnetic flux emergence. Observed by the Atmospheric Imaging Assembly, it features a rapid decrease in extreme-ultraviolet (EUV) emission in the 171 Å channel images, and a simultaneous increase in the 211 Å images. Here, we analyze the coronal thermodynamic and magnetic properties to probe its physical origin. We calculate the time-dependent differential emission measures for a sample of 18 events between 2010 and 2012. The emission measure (EM) decrease in the temperature range is well correlated with the EM increase in over eight orders of magnitude. This suggests that the coronal plasma is being heated from the quiet-Sun, sub-MK temperature to 1–2 MK, more typical for ARs. Potential field extrapolation indicates significant change in the local magnetic connectivity: the dimming region is now linked to the newly emerged flux via longer loops. We conclude that emerging dimming is likely caused by coronal heating episodes, powered by reconnection between the emerging and the ambient magnetic fields.

The relative role of turbulence, magnetic fields, and self-gravity in star formation is a subject of intensive debate. We present IRAM 30 m telescope observations of the ¹³CO (1–0) emission in the Serpens G3–G6 molecular cloud and apply to the data a set of statistical methods. These include the probability density functions of column density and the velocity gradients technique. We combine our data with the Planck 353 GHz polarized dust emission observations and Hershel H₂ column density. We suggest that the Serpens G3–G6 south clump is undergoing a gravitational collapse. Our analysis reveals that the gravitational collapse happens at volume density n ≥ 10³ cm⁻³. We estimate the plane-of-the-sky magnetic field strength to be approximately 120 μG using the traditional Davis–Chandrasekhar–Fermi method and 100 μG using a new technique proposed in Lazarian et al (2020). We find that the Serpens G3–G6 south clump’s total magnetic field energy significantly surpasses kinetic energy and gravitational energy. We conclude that the gravitational collapse could be successfully triggered in a supersonic and sub-Alfvénic cloud.

The properties of the first-discovered interstellar object (ISO), 1I/2017 (‘Oumuamua), differ from both solar system asteroids and comets, casting doubt on a protoplanetary disk origin. In this study, we investigate the possibility that it formed with a substantial H₂ ice component in the starless core of a giant molecular cloud. While interstellar solid hydrogen has yet to be detected, this constituent would explain a number of the ISO’s properties. We consider the relevant processes required to build decameter-sized, solid hydrogen bodies and assess the plausibility of growth in various size regimes. Via an energy balance argument, we find the most severe barrier to formation is the extremely low temperature required for the favorability of molecular hydrogen ice. However, if deposition occurs, we find that the turbulence within starless cores is conducive for growth into kilometer-sized bodies on sufficiently short timescales. Then, we analyze mass loss in the interstellar medium and determine the necessary size for a hydrogen object to survive a journey to the solar system as a function of ISO age. Finally, we discuss the implications if the H₂ explanation is correct, and we assess the future prospects of ISO science. If hydrogen ice ISOs do exist, our hypothesized formation pathway would require a small population of porous, 100 μm dust in a starless core region that has cooled to 2.8 K via adiabatic expansion of the surrounding gas and excellent shielding from electromagnetic radiation and cosmic rays.

Young massive clusters play an important role in the evolution of their host galaxies, and feedback from the high-mass stars in these clusters can have profound effects on the surrounding interstellar medium. The nuclear starburst in the nearby galaxy NGC 253 at a distance of 3.5 Mpc is a key laboratory in which to study star formation in an extreme environment. Previous high-resolution (1.9 pc) dust continuum observations from the Atacama Large Millimeter/submillimeter Array (ALMA) discovered 14 compact, massive super star clusters (SSCs) still in formation. We present here ALMA data at 350 GHz with 28 mas (0.5 pc) resolution. We detect blueshifted absorption and redshifted emission (P-Cygni profiles) toward three of these SSCs in multiple lines, including CS 7−6 and H¹³CN 4−3, which represent direct evidence for previously unobserved outflows. The mass contained in these outflows is a significant fraction of the cluster gas masses, which suggests we are witnessing a short but important phase. Further evidence of this is the finding of a molecular shell around the only SSC visible at near-IR wavelengths. We model the P-Cygni line profiles to constrain the outflow geometry, finding that the outflows must be nearly spherical. Through a comparison of the outflow properties with predictions from simulations, we find that none of the available mechanisms completely explains the observations, although dust-reprocessed radiation pressure and O star stellar winds are the most likely candidates. The observed outflows will have a very substantial effect on the clusters’ evolution and star formation efficiency.

We analyze the morphological evolution of open clusters and provide shape parameters for 265 open clusters. The results show that the overall shape of sample clusters becomes more elliptical as they grow older, while their core remains circular or slightly trending to circular. There is a negative correlation of the ellipticities with the number of members of the sample clusters. A significant negative correlation between the overall ellipticities and masses is also detected for the sample clusters with log(age/year) ≥8, suggesting that the overall shapes of the clusters are possibly influenced by the number of members and masses, in addition to the external forces and the surrounding environment. For most young sample clusters, the radial stratification degree of the short-axis direction is greater than that of the long, implying that the radial stratification degree in the two directions within the young sample cluster may be unevenly affected by an internal evolutionary process. Older sample clusters exhibit lower stratification in the tangential direction, which possibly means those clusters may continue to survive for a long time at a low level of stratification. Our analysis shows that the overall shape of the sample clusters may be more susceptible to the influence of Galactic tides toward the Galactic center than the shear forces embedded in Galactic differential rotation. By analyzing the distribution of the ages and number of members of star clusters, we suggest that NGC 6791 may originate from superclusters.

We present sensitive ALMA observations of TWA 3, a nearby, young (∼10 Myr) hierarchical system composed of three pre-main-sequence M3–M4.5 stars. For the first time, we detected ¹²CO and ¹³CO J = 2–1 emissions from the circumbinary protoplanetary disk around TWA 3A. We jointly fit the protoplanetary disk velocity field, stellar astrometric positions, and stellar radial velocities to infer the architecture of the system. The Aa and Ab stars (0.29 ± 0.01 M_⊙ and 0.24 ± 0.01 M_⊙, respectively) comprising the tight (P = 35 days) eccentric (e = 0.63 ± 0.01) spectroscopic binary are coplanar with their circumbinary disk (misalignment <6° with 68% confidence), similar to other short-period binary systems. From models of the spectral energy distribution, we found the inner radius of the circumbinary disk (r_inner = 0.50–0.75 au) to be consistent with theoretical predictions of dynamical truncation r_cav/a_inner ≈ 3. The outer orbit of the tertiary star B (0.40 ± 0.28 M_⊙, a ∼ 65 ± 18 au, e = 0.3 ± 0.2) is not as well constrained as the inner orbit; however, orbits coplanar with the A system are still preferred (misalignment < 20°). To better understand the influence of the B orbit on the TWA 3A circumbinary disk, we performed SPH simulations of the system and found that the outer edge of the gas disk (r_outer = 8.5 ± 0.2 au) is most consistent with truncation from a coplanar, circular, or moderately eccentric orbit, supporting the preference from the joint orbital fit.

Dense cores are the final place where turbulence is dissipated. It has been proposed from theoretical arguments that the nonthermal velocity dispersion should be narrower both for molecular ions (compared to neutrals) and for transitions with higher critical densities. To test these hypotheses, we compare the velocity dispersion of (1–0) (  = 6 × 10⁴) and (  = 2 × 10³), in the dense core Barnard 5. We analyze well-resolved and high signal-to-noise observations of (1,1) and (2,2) obtained with combining Robert C. Byrd Green Bank Telescope (GBT) and Very Large Array (VLA) data, and (1–0) obtained with GBT Argus, which present a similar morphology. Surprisingly, the nonthermal velocity dispersion of the ion is systematically higher than that of the neutral by 20%. The derived sonic Mach number, , has peak values and for and , respectively. This observed difference may indicate that the magnetic field even deep within the dense core is still oscillating, as it is in the turbulent region outside the core. The ions should be more strongly dynamically coupled to this oscillating field than the neutrals, thus accounting for their broader line width. If corroborated by further observations, this finding would shed additional light on the transition to quiescence in dense cores.

From our position embedded within the Milky Way’s interstellar medium, we have limited ability to detect gas at low relative velocities in the extended Galactic halo because those spectral lines are blended with much stronger signals from dense foreground gas. As a result, the content of the Milky Way’s circumgalactic medium (CGM) is poorly constrained at ∣v_LSR∣ ≲150 km s⁻¹. To overcome this complication, the QuaStar survey applies a spectral differencing technique using paired quasar−star sight lines to measure the obscured content of the Milky Way’s CGM for the first time. We present measurements of the C iv doublet (λλ1548, 1550), a rest-frame UV metal line transition, detected in Hubble Space Telescope/Cosmic Origins Spectrograph spectra of 30 halo-star/quasar pairs evenly distributed across the sky at Galactic latitudes ∣b∣ > 30°. The 30 halo stars have well-constrained distances (d ≈ 5–14 kpc) and are paired with quasars separated by <28. We argue that the difference in absorption between the quasar and stellar sight lines originates primarily in the Milky Way’s extended CGM beyond ∼10 kpc. For the Milky Way’s extended, low-velocity CGM (∣v∣ <150 km s⁻¹), we place an upper limit on the mean C iv column density of < 13.39 and find a covering fraction of 20% [6/30], a value significantly lower than the covering fraction for star-forming galaxies at low redshift. Our results suggest either that the bulk of Milky Way’s C iv -traced CGM lies at low Galactic latitudes or that the Milky Way’s CGM is lacking in warm, ionized material compared to low-redshift (z < 0.1) star-forming galaxy halos.

We present initial results from the Cosmic Origins Spectrograph (COS) and Gemini Mapping the Circumgalactic Medium (CGMCGM ≡ CGM²) survey. The CGM² survey consists of 1689 galaxies, all with high-quality Gemini-GMOS spectra, within 1 Mpc of 22 z ≲ 1 quasars, all with a signal-to-noise ratio of ∼10 Hubble Space Telescope/COS G130M+G160M spectra. For 572 of these galaxies with stellar masses 10⁷ M_⊙ < M_⋆ < 10¹¹ M_⊙ and z ≲ 0.5, we show that the H i covering fraction above a threshold of N_HI > 10¹⁴cm⁻² is ≳0.5 within 1.5 virial radii (R_vir ∼ R_200m). We examine the H i kinematics and find that the majority of absorption lies within ±250 km s⁻¹ of the galaxy systemic velocity. We examine H i covering fractions over a range of impact parameters to infer a characteristic size of the CGM, , as a function of galaxy mass. is the impact parameter at which the probability of observing an absorber with N_HI >10¹⁴ cm⁻² is >50%. In this framework, the radial extent of the CGM of M_⋆ > 10^9.9 M_⊙ galaxies is kpc or . Intermediate-mass galaxies with 10^9.2 < M_⋆/M_⊙ < 10^9.9 have an extent of kpc or . Low-mass galaxies, M_⋆ < 10^9.2 M_⊙, show a smaller physical scale of kpc and extend to . Our analysis suggests that using R_vir as a proxy for the characteristic radius of the CGM likely underestimates its extent.

Using six years of spectroscopic monitoring of the luminous quasar HE 0435-4312 (z = 1.2231) with the Southern African Large Telescope, in combination with photometric data (CATALINA, OGLE, SALTICAM, and BMT), we determined a rest-frame time delay of days between the Mg ii broad-line emission and the ionizing continuum using seven different time-delay inference methods. Time-delay artifact peaks and aliases were mitigated using the bootstrap method and prior weighting probability function, as well as by analyzing unevenly sampled mock light curves. The Mg ii emission is considerably variable with a fractional variability of ∼5.4%, which is comparable to the continuum variability (∼4.8%). Because of its high luminosity (L₃₀₀₀ = 10^46.4 erg s⁻¹), the source is beneficial for a further reduction of the scatter along the Mg ii-based radius–luminosity relation and its extended versions, especially when the highly accreting subsample that has an rms scatter of ∼0.2 dex is considered. This opens up the possibility of using the high-accretor Mg ii-based radius–luminosity relation for constraining cosmological parameters. With the current sample of 27 reverberation-mapped sources, the best-fit cosmological parameters (Ω_m, Ω_Λ) = (0.19; 0.62) are consistent with the standard cosmological model within the 1σ confidence level.

We present multiband Hubble Space Telescope imaging that spans rest-frame near-ultraviolet through near-infrared wavelengths (–1.1 μm) for 12 compact starburst galaxies at z = 0.4–0.8. These massive galaxies () are driving very fast outflows (–3000 km s⁻¹), and their light profiles are dominated by an extremely compact starburst component (half-light radius ∼ 100 pc). Our goal is to constrain the physical mechanisms responsible for launching these fast outflows by measuring the physical conditions within the central kiloparsec. Based on our stellar population analysis, the central component typically contributes ≈25% of the total stellar mass, and the central escape velocities km s⁻¹ are a factor of two smaller than the observed outflow velocities. This Requires physical mechanisms that can accelerate gas to speeds significantly beyond the central escape velocities, and it makes clear that these fast outflows are capable of traveling into the circumgalactic medium, and potentially beyond. We find central stellar densities kpc⁻² comparable to theoretical estimates of the Eddington limit, and we estimate surface densities within the central kiloparsec comparable to those of compact massive galaxies at . Relative to “red nuggets” and “blue nuggets” at , we find significantly smaller r_e values at a given stellar mass, which we attribute to the dominance of a young stellar component in our sample and the better physical resolution for rest-frame optical observations at versus . We compare to theoretical scenarios involving major mergers and violent disk instability, and we speculate that our galaxies are progenitors of power-law ellipticals in the local universe with prominent stellar cusps.

Pox 186 is an exceptionally small dwarf starburst galaxy hosting a stellar mass of ∼10⁵M_⊙. Undetected in H i (M < 10⁶M_⊙) from deep 21 cm observations and with an [O iii]/[O ii] (5007/3727) ratio of 18.3 ± 0.11, Pox 186 is a promising candidate Lyman continuum emitter. It may be a possible analog of low-mass reionization-era galaxies. We present a spatially resolved kinematic study of Pox 186 and identify two distinct ionized gas components: a broad one with σ > 400 km s⁻¹ and a narrow one with σ < 30 km s⁻¹. We find strikingly different morphologies between the two components and direct evidence of outflows as seen in the high-velocity gas. Possible physical mechanisms driving the creation of high-velocity gas seen in [O iii] are discussed, from outflow geometry to turbulent mixing between a hot (10⁶ K) star-cluster wind and cooler (10⁴ K) gas clouds. We find a modest mass-outflow rate of 0.022 M_⊙ yr⁻¹ with a small mass-loading factor of 0.5, consistent with other low-mass galaxies. Finally, we compare the mass-loading factor of Pox 186 with extrapolations from numerical simulations and discuss possible reasons for the apparent discrepancy between them.

We study the effects of additional cooling due to the emission of a dark matter candidate particle, the dark photon, on the final phases of the evolution of a 15 M_⊙ star and resulting modifications of the pre-supernova neutrino signal. For a substantial portion of the dark photon parameter space the extra cooling speeds up Si burning, which results in a reduced number of neutrinos emitted during the last day before core collapse. This reduction can be described by a systematic acceleration of the relevant timescales and the results can be estimated semi-analytically in good agreement with the numerical simulations. Outside the semi-analytic regime we find more complicated effects. In a narrow parameter range, low-mass dark photons lead to an increase in the number of emitted neutrinos because of additional shell-burning episodes that delay core collapse. Furthermore, relatively strong couplings produce a thermonuclear runaway during O burning, which could result in a complete disruption of the star but requires more detailed simulations to determine the outcome. Our results show that pre-supernova neutrino signals are a potential probe of the dark photon parameter space.

One possible progenitor of short gamma-ray bursts (GRBs) is thought to be from a double neutron star (NS) merger, and the remnant of such a merger may be a supramassive NS, which is supported by rigid rotation and through its survival of hundreds of seconds before collapsing into a black hole (BH). If this is the case, an optical/infrared transient (namely merger-nova) is generated from the ejected materials and it is powered by radioactive decay from r-process, spin-down energy from a supramassive NS, as well as the magnetic wind from a newborn BH. In this paper, we systematically search for the signature of a supramassive NS central engine by analyzing the X-ray emission of short GRBs with internal plateau observed by Swift, and we find that five candidates of short GRBs have such a feature with redshift measurement. Then, we calculate the possible merger-nova emission from those candidates given the typical model parameters by considering the above three energy sources, and compare its brightness with the sensitivity of some optical telescopes. We find that the merger-nova emission of GRB 060801 in K-, r-, and U-bands with variations of M_ej (10⁻⁴–10⁻²M_⊙), κ (0.1–10 cm² g⁻¹), and β (0.1–0.3) is very difficult to detect using the Vera C. Rubin, Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Zwicky Transient Facility, and the Roman Space Telescope (Roman), except for the case of large ejecta mass M_ej = 10⁻²M_⊙. However, we are very hopeful to detect the merger-nova emission of GRBs 090515, 100625A, and 101219A using more sensitive instruments, such as Vera C. Rubin, Pan-STARRS, and Roman. Moreover, the merger-nova emission of GRB 160821B is bright enough to detect in our calculations, and it is also consistent with current real observations of merger-nova emission.

The Habitable-zone Planet Finder (HPF) is a fiber-fed precise radial velocity (RV) spectrograph at the 10 m Hobby–Eberly Telescope (HET). Due to its fixed-altitude design, the HET pupil changes appreciably across a track, leading to significant changes of the fiber far-field illumination. HPF’s fiber scrambler is designed to suppress the impact of these illumination changes on the RVs—but the residual impact on the RV measurements has yet to be probed on-sky. We use GJ 411, a bright early type (M2) M dwarf to probe the effects of far-field input trends due to these pupil variations on HPF RVs. These large changes (∼2x) in the pupil area and centroid present a harsh test of HPF’s far-field scrambling. Our results show that the RVs are effectively decoupled from these extreme far-field input changes due to pupil centroid offsets, attesting to the effectiveness of the scrambler design. This experiment allows us to test the impact of these changes with large pupil variation on-sky, something we would not easily be able to do at a conventional optical telescope. While the pupil and illumination changes expected at these other telescopes are small, scaling from our results enables us to estimate and bound these effects, and show that they are controllable even for the new and next generation of RV instruments in their quest to beat down instrumental noise sources toward the goal of a few .

The cluster Westerlund 1 (Wd1) is host to a large variety of post-main-sequence (MS) massive stars. The simultaneous presence of these stars can only be explained by stellar models if the cluster has a finely tuned age of 4–5 Myr, with several published studies independently claiming ages within this range. At this age, stellar models predict that the cool supergiants (CSGs) should have luminosities of , close to the empirical luminosity limit. Here, we test that prediction using archival data and new photometry from Stratospheric Observatory for Infrared Astronomy to estimate bolometric luminosities for the CSGs. We find that these stars are on average 0.4 dex too faint to be 5 Myr old, regardless of which stellar evolutionary model is used, and instead are indicative of a much older age of Myr. We argue that neither systematic uncertainties in the extinction law nor stellar variability can explain this discrepancy. In reviewing various independent age estimates of Wd1 in the literature, we first show that those based on stellar diversity are unreliable. Second, we reanalyze Wd1's pre-MS stars employing the Damineli extinction law, finding an age of Myr; older than that of previous studies, but which is vulnerable to systematic errors that could push the age close to 10 Myr. However, there remains significant tension between the CSG age and that inferred from the eclipsing binary W13. We conclude that stellar evolutionary models cannot explain Wd1 under the single-age paradigm. Instead, we propose that the stars in the Wd1 region formed over a period of several megayears.

We report results of the first broadband observation of the transient X-ray pulsar GRO J1008−57 performed in the quiescent state. Observations were conducted quasi-simultaneously with NuSTAR and the Mikhail Pavlinsky ART-XC telescope on board SRG right before the beginning of a Type I outburst. GRO J1008−57 was detected in the state with the lowest observed luminosity around several ×10³⁴ erg s⁻¹ and consequently accreting from the cold disk. Timing analysis allowed us to detect pulsations during this state for the first time. The observed pulsed fraction of about 20% is, however, almost three times lower than in brighter states when the accretion proceeds through the standard disk. We traced the evolution of the broadband spectrum of the source on a scale of three orders of magnitude in luminosity and found that at the lowest luminosities the spectrum transforms into the double-hump structure similarly to other X-ray pulsars accreting at low luminosities (X Persei, GX 304–1, and A 0535+262) reinforcing the conclusion that this spectral shape is typical for these objects.

Prompted by a recent paper by Dima and Schad, we reconsider the problem of inferring magnetic properties of the corona using polarimetric observations of magnetic dipole (M1) lines. Dima and Schad point to a potential source of degeneracy in a formalism developed by Plowman, which under some circumstances can lead to the solution being under-determined. Here we clarify the nature of the problem. Its resolution lies in solving for the scattering geometry using the elongation of the observed region of the corona. We discuss some conceptual problems that arise when casting the problem for inversion in the observer’s reference frame, and satisfactorily resolve difficulties identified by Plowman, Dima, and Schad.

The nova rate in the Milky Way remains largely uncertain, despite its vital importance in constraining models of Galactic chemical evolution as well as understanding progenitor channels for Type Ia supernovae. The rate has been previously estimated to be in the range of ≈10–300 yr⁻¹, either based on extrapolations from a handful of very bright optical novae or the nova rates in nearby galaxies; both methods are subject to debatable assumptions. The total discovery rate of optical novae remains much smaller (≈5–10 yr⁻¹) than these estimates, even with the advent of all-sky optical time-domain surveys. Here, we present a systematic sample of 12 spectroscopically confirmed Galactic novae detected in the first 17 months of Palomar Gattini-IR (PGIR), a wide-field near-infrared time-domain survey. Operating in the J band (≈1.2 μm), which is significantly less affected by dust extinction compared to optical bands, the extinction distribution of the PGIR sample is highly skewed to a large extinction values (>50% of events obscured by A_V ≳ 5 mag). Using recent estimates for the distribution of Galactic mass and dust, we show that the extinction distribution of the PGIR sample is commensurate with dust models. The PGIR extinction distribution is inconsistent with that reported in previous optical searches (null-hypothesis probability <0.01%), suggesting that a large population of highly obscured novae have been systematically missed in previous optical searches. We perform the first quantitative simulation of a 3π time-domain survey to estimate the Galactic nova rate using PGIR, and derive a rate of yr⁻¹. Our results suggest that all-sky near-infrared time-domain surveys are well poised to uncover the Galactic nova population.

“Changing-look” quasars are a new class of highly variable active galactic nuclei that have changed their spectral type over surprisingly short timescales of just a few years. The origin of this phenomenon is debated, but is likely to reflect some change in the accretion flow. To investigate the disk–corona systems in these objects, we measure optical/UV–X-ray spectral indices () and Eddington ratios () of 10 previously discovered changing-look quasars at two or more epochs. By comparing these data with simulated results based on the behavior of X-ray binaries, we find possible similarities in spectral indices below the 1% Eddington ratio. We further investigate the Eddington ratios of changing-look quasars before and after their spectral type changes, and find that changing-look quasars cross the 1% Eddington ratio boundary when their broad emission lines disappear/emerge. This is consistent with the disk-wind model as the origin of broad emission lines.

We present the largest and deepest late-time radio and millimeter survey to date of superluminous supernovae (SLSNe) and long-duration gamma-ray bursts (LGRBs) to search for associated nonthermal synchrotron emission. Using the Karl G. Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA), we observed 43 sources at 6 and 100 GHz on a timescale of ∼ 1–19 yr post-explosion. We do not detect radio/millimeter emission from any of the sources, with the exception of a 6 GHz detection of PTF10hgi, as well as the detection of 6 GHz emission near the location of the SLSN PTF12dam, which we associate with its host galaxy. We use our data to place constraints on central engine emission due to magnetar wind nebulae and off-axis relativistic jets. We also explore nonrelativistic emission from the SN ejecta, and place constraints on obscured star formation in the host galaxies. In addition, we conduct a search for fast radio bursts (FRBs) from some of the sources using VLA phased-array observations; no FRBs are detected to a limit of 16 mJy (7σ; 10 ms duration) in about 40 minutes on source per event. A comparison to theoretical models suggests that continued radio monitoring may lead to detections of persistent radio emission on timescales of ≳ a decade.

Using neural networks, we integrate the ability to account for Doppler smearing due to a pulsar’s orbital motion with the pulsar population synthesis package psrpoppy to develop accurate modeling of the observed binary pulsar population. As a first application, we show that binary neutron star systems where the two components have highly unequal mass are, on average, easier to detect than systems that are symmetric in mass. We then investigate the population of ultracompact (1.5 minutes ≤ P_b ≤ 15 minutes) neutron star–white dwarf (NS–WD) and double neutron star (DNS) systems, which are promising sources for the Laser Interferometer Space Antenna gravitational-wave detector. Given the nondetection of these systems in radio surveys thus far, we estimate a 95% confidence upper limit of ∼1450 and ∼1100 ultracompact NS–WD and DNS systems in the Milky Way that are beaming toward the Earth, respectively. We also show that using survey integration times in the range 20 s–200 s with time-domain resampling will maximize the signal-to-noise ratio as well as the probability of detection of these ultracompact binary systems. Among all the large-scale radio pulsar surveys, those that are currently being carried out using archival data collected with the Arecibo radio telescope have a ∼50%–80% chance of detecting at least one of these systems using current integration integration times and ∼80%–95% using optimal integration times in the next several years.

The progenitors of Type Ia supernovae (SNe Ia) are debated, particularly the evolutionary state of the binary companion that donates mass to the exploding carbon–oxygen white dwarf. In our previous work, we presented hydrodynamic models and optically thin radio synchrotron light curves of SNe Ia interacting with detached, confined shells of CSM, representing CSM shaped by novae. In this work, we extend these light curves to the optically thick regime, considering both synchrotron self-absorption and free–free absorption. We obtain simple formulae to describe the evolution of optical depth seen in the simulations, allowing optically thick light curves to be approximated for arbitrary shell properties. We then demonstrate the use of this tool by interpreting published radio data. First, we consider the nondetection of PTF11kx—an SN Ia known to have a detached, confined shell—and we find that the nondetection is consistent with current models for its CSM, and that observations at a later time would have been useful for this event. Second, we statistically analyze an ensemble of radio nondetections for SNe Ia with no signatures of interaction. We find that shells with masses (10⁻⁴–0.3) M_⊙ located (10¹⁵–10¹⁶) cm from the progenitor are currently not well constrained by radio datasets, due to their dim, rapidly evolving light curves.

Recent studies of gamma-ray, cosmic-ray, and radio data put stringent constraints on the fraction of primordial black holes (PBHs) in our universe. In this article, we propose a new indirect method in using the X-ray luminosity data of cool-core clusters to constrain the evaporating PBH fraction for the monochromatic, log-normal and power-law mass distributions. The present results show that the amount of evaporating PBHs only constitutes a minor component of dark matter for a large parameter space. The constraints are consistent with and close to that obtained from other cosmic-ray and multiwavelength observations.

Chromospheric condensation is a brief episode of downflow often accompanying energy release and evaporation in a solar flare. While this component of a flare reflects the energy release process only indirectly, it can be observed at high spatial and temporal resolution, even from the ground. It appears in spectroscopic observations of cooler lines, formed below ∼10⁵ K, as a redshift that peaks and decays after less than 1 minute. In order to use this signature to infer characteristics of solar flare energy release, it is important to establish quantitative relationships with properties of the condensation. The initial investigation reported here does so after restricting consideration to energy transport via thermal conduction into a simplified, stratified chromosphere. We develop an analytical expression for the decay of a condensation propagating into a stratified atmosphere. This model provides a relationship between shock velocity and preshock density structure. We also use one-dimensional gasdynamic simulations to explore the dynamics of these shocks as they penetrate into the stratified chromosphere. We find that the peak downflow speed primarily reflects the energy flux into the chromosphere, while the product of this velocity and the redshift duration is proportional to the preshock density scale height as H ≃ 0.6u₀τ.

Future microcalorimeter X-ray observations will resolve spectral features in unmatched detail. Understanding of line formation processes in X-rays deserves much attention. The purpose of this paper is to discuss such processes in the presence of a photoionizing source. Line formation processes in one- and two-electron species are broadly categorized into four cases. Case A occurs when the Lyman line optical depths are very small and photoexcitation does not occur. Line photons escape the cloud without any scattering. Case B occurs when the Lyman line optical depths are large enough for photons to undergo multiple scatterings. Case C occurs when a broadband continuum source strikes an optically thin cloud. The Lyman lines are enhanced by induced radiative excitation of the atoms/ions by continuum photons, also known as continuum pumping. A fourth, less studied scenario, where the Case B spectrum is enhanced by continuum pumping, is called Case D. Here, we establish the mathematical foundation of Cases A, B, C, and D in an irradiated cloud with Cloudy. We also show the total X-ray emission spectrum for all four cases within the energy range 0.1–10 keV at the resolving power of XRISM around 6 keV. Additionally, we show that the combined effect of electron scattering and partial blockage of continuum pumping reduces the resonance line intensities. Such reduction increases with column density and can serve as an important tool to measure the column density/optical depth of the cloud.

The color–magnitude diagrams of young star clusters show that, particularly at ultraviolet wavelengths, their upper main sequences (MSs) bifurcate into a sequence comprising the bulk population and a blue periphery. The spatial distribution of stars is crucial to understanding the reasons for these distinct stellar populations. This study uses high-resolution photometric data obtained with the Hubble Space Telescope to study the spatial distributions of the stellar populations in seven Magellanic Cloud star clusters. The cumulative radial number fractions of blue stars within four clusters are strongly anticorrelated with those of the high mass ratio binaries in the bifurcated region, with negative Pearson coefficients <−0.7. Those clusters are generally young or in an early dynamical evolutionary stage. In addition, a supporting N-body simulation suggests that the increasing percentage of blue MS stars from the cluster centers to their outskirts may be associated with the dissolution of soft binaries. This study provides a different perspective to explore the MS bimodalities in young clusters and adds extra puzzles. A more comprehensive study combined with detailed simulations is needed in the future.

The Parker Solar Probe (PSP) routinely observes magnetic field deflections in the solar wind at distances less than 0.3 au from the Sun. These deflections are related to structures commonly called “switchbacks” (SBs), whose origins and characteristic properties are currently debated. Here, we use a database of visually selected SB intervals—and regions of solar wind plasma measured just before and after each SB—to examine plasma parameters, turbulent spectra from inertial to dissipation scales, and intermittency effects in these intervals. We find that many features, such as perpendicular stochastic heating rates and turbulence spectral slopes are fairly similar inside and outside of SBs. However, important kinetic properties, such as the characteristic break scale between the inertial to dissipation ranges differ inside and outside these intervals, as does the level of intermittency, which is notably enhanced inside SBs and in their close proximity, most likely due to magnetic field and velocity shears observed at the edges. We conclude that the plasma inside and outside of an SB, in most of the observed cases, belongs to the same stream, and that the evolution of these structures is most likely regulated by kinetic processes, which dominate small-scale structures at the SB edges.

Convection and turbulence in core-collapse supernovae (CCSNe) are inherently three-dimensional (3D) in nature. However, 3D simulations of CCSNe are computationally demanding. Thus, it is valuable to modify simulations in spherical symmetry to incorporate 3D effects using some parametric model. In this paper, we report on the formulation and implementation of general relativistic neutrino-driven turbulent convection in the spherically symmetric core-collapse supernova code GR1D. This is based upon the recently proposed method of Simulated Turbulence in Reduced Dimensionality (STIR) in Newtonian simulations from Couch et al. (2020). When the parameters of this model are calibrated to 3D simulations, we find that our general relativistic formulation of STIR requires larger turbulent eddies to achieve a shock evolution similar to the original STIR model. We also find that general relativity may alter the correspondence between progenitor mass and successful versus failed explosions.

We present optical and near-infrared observations of SN 2019ehk, which was initially reported as a Type Ib supernova (SN). We show that it evolved to a Ca-rich transient according to its spectral properties and evolution in late phases. However, it shows a few properties distinct from those of the canonical Ca-rich transients: a short-duration first peak in the light curve, high peak luminosity, and association with a star-forming environment. Indeed, some of these features are shared with iPTF14gqr and iPTF16hgs, which are candidates for a special class of core-collapse SNe: the so-called ultra-stripped envelope SNe, i.e., a relatively low-mass He (or C+O) star explosion in a binary as a precursor of short-period double neutron star (NS) binaries. The estimated ejecta mass (0.4M_⊙) and explosion energy (1.7 × 10⁵⁰ erg) are consistent with this scenario. The analysis of the first peak suggests the existence of dense circumstellar material in the vicinity of the progenitor, implying a CCSN origin. Based on this analysis, we suggest SN 2019ehk is another candidate for a low-mass He star explosion. It might create a double NS binary, but with a wide separation. These candidates for low-mass stripped envelope SNe, including ultra-stripped envelope SN candidates, seem to form a subpopulation among Ca-rich transients, associated with young population. We propose that the key to distinguishing this population is the early first peak in their light curves.

Building on recent work by Chandar et al., we construct X-ray luminosity functions (XLFs) for different classes of X-ray binary (XRB) donors in the nearby star-forming galaxy M83 through a novel methodology. Rather than classifying low- versus high-mass XRBs based on the scaling of the number of X-ray sources with stellar mass and star formation rate, respectively, we utilize multiband Hubble Space Telescope imaging data to classify each Chandra-detected compact X-ray source as a low-mass (i.e., donor mass ≲3 M_⊙), high-mass (donor mass ≳8M_⊙), or intermediate-mass XRB based on either the location of its candidate counterpart on optical color–magnitude diagrams or the age of its host star cluster. In addition to the standard (single and/or truncated) power-law functional shape, we approximate the resulting XLFs with a Schechter function. We identify a marginally significant (at the 1σ-to-2σ level) exponential downturn for the high-mass XRB XLF, at (in log CGS units). In contrast, the low- and intermediate-mass XRB XLFs, as well as the total XLF of M83, are formally consistent with sampling statistics from a single power law. Our method suggests a non-negligible contribution from low- and possibly intermediate-mass XRBs to the total XRB XLF of M83, i.e., between 20% and 50%, in broad agreement with X-ray-based XLFs. More generally, we caution against considerable contamination from X-ray emitting supernova remnants to the published, X-ray-based XLFs of M83, and possibly all actively star-forming galaxies.

In this paper, we collect a sample of stars observed both in LAMOST and Gaia, which have colors implying a temperature hotter than 7000 K. We train a machine-learning algorithm on LAMOST spectroscopic data which has been tagged with stellar classifications and metallicities, and use this machine to construct a catalog of blue horizontal branch stars (BHBs), together with metallicity information. Another machine is trained using Gaia parallaxes to predict absolute magnitudes for these stars. The final catalog of 13,693 BHBs is thought to be about 86% pure, with σ_[Fe/H] ∼ 0.35 dex, and σ_G ∼ 0.31 mag. These values are confirmed via comparison to globular clusters, although a covariance error seems to affect our magnitude and abundance estimates. We analyze a subset of this catalog in the Galactic Halo. We find that BHB populations in the outer halo appear redder, which could imply a younger population, and that the metallicity gradient is relatively flat around [Fe/H] = −1.9 dex over our sample footprint. We find that our metal-rich BHB stars are on more radial velocity dispersion-dominated orbits (β ∼ 0.70) at all radii than our metal-poor BHB stars (β ∼ 0.62).

We present new proper-motion measurements of optically emitting oxygen-rich knots of supernova remnant 1E 0102.2-7219 (E0102), which are used to estimate the remnant’s center of expansion and age. Four epochs of high-resolution Hubble Space Telescope images spanning 19 yr were retrieved and analyzed. We found a robust center of expansion of α = 1^h04^m0248 and δ = −72°01′5392 (J2000) with 1σ uncertainty of 177 using 45 knots from images obtained with the Advanced Camera for Surveys using the F475W filter in 2003 and 2013 having the highest signal-to-noise ratio. We also estimate an upper limit explosion age of 1738 ± 175 yr by selecting the knots with the highest proper motions and these knots are assumed to be the least decelerated. We find evidence of an asymmetry in the proper motions of the knots as a function of position angle. We conclude that these asymmetries were most likely caused by interaction between E0102's original supernova blast wave and an inhomogeneous surrounding environment, as opposed to intrinsic explosion asymmetry. The observed nonhomologous expansion suggests that the use of a free expansion model inaccurately offsets the center of expansion and leads to an overestimated explosion age. We discuss our findings as they compare to previous age and center of expansion estimates of E0102, as well as their relevance to a recently identified candidate central compact object.

We present the results of ALMA observations toward the low-mass Class 0 binary system VLA 1623Aab in the Ophiuchus molecular cloud in ¹²CO, ¹³CO, and C¹⁸O(2–1) lines. Our ¹²CO (J = 2–1) data reveal that the VLA 1623 outflow consists of twin spatially overlapped outflows/jets. The redshifted northwestern jet exhibits three cycles of wiggle with a spatial period of 1360 ± 10 au, corresponding to a time period of 180 yr. The wiggle-like structure is also found in the position–velocity (PV) diagram, showing an amplitude in the velocity of about 0.9 km s⁻¹. Both the period and velocity amplitude of the wiggle are roughly consistent with those expected from the binary parameters, i.e., the orbital period (460 ± 20 yr) and the Keplerian velocity (2.2 km s⁻¹). Our ¹³CO and C¹⁸O images show a dense gas nature in the two centimeter/millimeter sources, VLA 1623B and W, and its relation to the outflows, and strongly support the previous interpretation that both are shocked cloudlets. The driving sources of the twin molecular outflows are, therefore, likely to be the VLA 1623Aab binary. The outflow axes of the two molecular outflows are estimated to be inclined by 70° to each other across the plane of sky, implying that protostellar disks are also misaligned by . Such nature together with a small binary separation of 34 au in one of the youngest protobinary systems seems difficult to explain by disk fragmentation in quiescent environments. Other effects such as turbulence probably play roles.

Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov and conservation of the Walker–Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images.

We continue to investigate two-dimensional laterally propagating flames in type I X-ray bursts using fully compressible hydrodynamics simulations. In the current study we relax previous approximations where we artificially boosted the flames. We now use more physically realistic reaction rates, thermal conductivities, and rotation rates, exploring the effects of neutron star rotation rate and thermal structure on the flame. We find that at lower rotation rates the flame becomes harder to ignite, whereas at higher rotation rates the nuclear burning is enhanced by increased confinement from the Coriolis force and the flame propagates steadily. At higher crustal temperatures, the flame moves more quickly and accelerates as it propagates through the atmosphere. If the temperature is too high, instead of a flame propagating across the surface the entire atmosphere burns uniformly. Our findings could have implications for the relationship between observed burst rise times and neutron star rotation and accretion rates. All of the software used for these simulations is freely available.

The largest moon of Neptune, Triton, possesses a cold and hazy atmosphere. Since the discovery of the near-surface haze layer during the Voyager fly in 1989, the haze formation mechanism has not been investigated in detail. Here we provide the first haze microphysical model on Triton. Our model solves the evolution of both size and porosity distributions of haze particles in a self-consistent manner. We simulated the formation of sphere and aggregate hazes with and without condensation of the C₂H₄ ice. The haze particles can grow into fractal aggregates with mass-equivalent sphere sizes of ∼0.1–1 μm and fractal dimensions of D_f = 1.8–2.2. The ice-free hazes cannot simultaneously explain both UV and visible observations of Voyager 2, while including the condensation of C₂H₄ ices provides two better solutions. For ice aggregates, the required total haze mass flux is ∼2 × 10⁻¹⁵ g cm⁻² s⁻¹. For the icy sphere scenario, the column-integrated C₂H₄ production rate is ∼8 × 10⁻¹⁵ g cm⁻² s⁻¹, and the ice-free mass flux is ∼6 × 10⁻¹⁷ g cm⁻² s⁻¹. The UV occultation observations at short wavelengths, <0.15 μm, may slightly favor the icy aggregates. Observations of the haze optical depth and the degree of forward scattering in UV and visible should be able to distinguish whether Triton’s hazes are icy spheres or ice aggregates in future Triton missions.

We present new high-resolution (R = λ/Δλ ∼ 14,000) spectra of the two brightest He ii-transparent quasars in the far-UV at z_em > 3.5, HE2QS J2311−1417 (z_em = 3.70) and HE2QS J1630+0435 (z_em = 3.81), obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. In the predominantly saturated He ii absorption spectra, both sight lines show several isolated resolved (FWHM ≳ 50 km s⁻¹) transmission spikes in He ii Lyα and He ii Lyβ. The incidence of such spikes decreases with increasing redshift, but both sight lines show significant spikes at z ≳ 3.5, signaling the presence of fully ionized regions in the z ≳ 3.5 intergalactic medium (IGM). We employ an automated algorithm to compare the number of detected He ii transmission spikes to predictions from mock spectra created from the outputs of a cubic (146 cMpc)³ optically thin Nyx hydrodynamical simulation, assuming a range of UV background photoionization rates Γ_{He II}. From the incidence of Lyα and Lyβ transmission spikes we infer similar photoionization rates of at 3.51 < z < 3.66 and Γ_{He II} = (0.9 ± 0.3) × 10⁻¹⁵ s⁻¹ at 3.460 < z < 3.685, respectively. Because the transmission spikes indicate fully ionized regions at z ≳ 3.5 along both lines of sight, our observations provide further evidence that He ii reionization had substantially progressed at these redshifts.

The recent advent of the Event Horizon Telescope (EHT) has made direct imaging of supermassive black holes a reality. Simulated images of black holes produced via general relativistic ray tracing and radiative transfer provide a key counterpart to these observational efforts. Black hole images have a wide range of physically interesting image structures, ranging from extremely fine scales in their lensed “photon rings” to the very large scales in their relativistic jets. The multiscale nature of the black hole system is therefore suitable for a multiscale approach to generate simulated images that capture all key elements of the system. Here, we present a prescription for adaptive ray tracing, which enables efficient computation of extremely high-resolution images of black holes. Using the polarized ray-tracing code ipole, we image a combination of semianalytic and general relativistic magnetohydrodynamic (GRMHD) models, and we show that images can be reproduced with a mean squared error of less than 0.1% even after tracing 12× fewer rays. We then use adaptive ray tracing to explore the properties of the photon ring. We illustrate the behavior of individual subrings in GRMHD simulations, and we explore their signatures in interferometric visibilities.

The exact mechanism for the production of fast γ-ray variability in blazars remains debated. Magnetic reconnection, in which plasmoids filled with relativistic particles and magnetic fields are formed, is a viable candidate to explain the broadband electromagnetic spectrum and variability of these objects. Using state-of-the-art magnetic reconnection simulations, we generate realistic γ-ray light curves that would be observed with the Fermi Large Area Telescope. A comparison with observed γ-ray flares from flat spectrum radio quasars (FSRQs) reveals that magnetic reconnection events lead to comparable flux levels and variability patterns, in particular, when the reconnection layer is slightly misaligned with the line of sight. Emission from fast plasmoids moving close to the line of sight could explain the fast variability on the timescales of minutes for which evidence has been found in observations of FSRQs. Our results motivate improvements in existing radiative transfer simulations as well as dedicated searches for fast variability as evidence for magnetic reconnection events.

We present a compendium of disk galaxy scaling relations and a detailed characterization of their intrinsic scatter. Observed scaling relations are typically characterized by their slope, intercept, and scatter; however, these parameters are a mixture of observational errors and astrophysical processes. We introduce a novel Bayesian framework for computing the intrinsic scatter of scaling relations that accounts for nonlinear error propagation and covariant uncertainties. Bayesian intrinsic scatters are ∼25% more accurate than those obtained with a first-order classical method, which systematically underestimates the true intrinsic scatter. Structural galaxy scaling relations based on velocity (V_23.5), size (R_23.5), luminosity (L_23.5), color (g − z), central stellar surface density (Σ₁), stellar mass (M_*), dynamical mass (M_dyn), stellar angular momentum (j_*), and dynamical angular momentum (j_dyn) are examined to demonstrate the power and importance of the Bayesian formalism. Our analysis is based on a diverse selection of over 1000 late-type galaxies from the Photometry and Rotation Curve Observations from Extragalactic Surveys compilation with deep optical photometry and extended rotation curves. We determine the tightest relation for each parameter by intrinsic orthogonal scatter, finding M_* − V_23.5, R_23.5 − j_*, and L_23.5 − j_dyn to be especially tight. The scatter of the R_23.5 − L_23.5, V_23.5 − (g − z), and R_23.5 − j_dyn relations is mostly intrinsic, making them ideal for galaxy formation and evolutionary studies. Our code to compute the Bayesian intrinsic scatter of any scaling relation is also presented. We quantify the correlated nature of many uncertainties in galaxy scaling relations and scrutinize the uncertain nature of disk inclination corrections and their effect on scatter estimates.

X-ray reverberation mapping has emerged as a new tool to probe accretion in active galactic nuclei (AGN), providing a potentially powerful probe of accretion at the black hole scale. The lags, along with relativistic spectral signatures are often interpreted in light of the lamppost model. Focusing specifically on testing the prediction of the relativistic reverberation model, we have targeted several of the brightest Seyfert Galaxies in X-rays with different observing programs. Here, we report the results from two large campaigns with NuSTAR targeting MCG-5-23-16 and SWIFT J2127.4+5654 to test the model predictions in the 3–50 keV band. These are two of three sources that showed indications of a delayed Compton hump in early data. With triple the previously analyzed exposures, we find no evidence for relativistic reverberation in MCG-5-23-16, and the energy-dependent lags are consistent with a log-linear continuum. In SWIFT J2127.4+5654, although a continuum-only model explains the data, the relativistic reverberation model provides a significant improvement to the energy and frequency-dependent lags, but with parameters that are not consistent with the time-averaged spectrum. This adds to mounting evidence showing that the lag data is not consistent with a static lamppost model.

Dynamical friction is typically regarded as a secular process, in which the subject (“perturber”) evolves very slowly (secular approximation) and has been introduced to the host over a long time (adiabatic approximation). These assumptions imply that dynamical friction arises from the LBK torque with nonzero contribution only from purely resonant orbits. However, dynamical friction is only of astrophysical interest if its timescale is shorter than the age of the universe. In this paper we therefore relax the adiabatic and secular approximations. We first derive a generalized LBK torque, which reduces to the LBK torque in the adiabatic limit, and show that it gives rise to transient oscillations due to nonresonant orbits that slowly damp out, giving way to the LBK torque. This is analogous to how a forced, damped oscillator undergoes transients before settling to a steady state, except that here the damping is due to phase mixing rather than dissipation. Next, we present a self-consistent treatment that properly accounts for time dependence of the perturber potential and circular frequency (memory effect), which we use to examine orbital decay in a cored galaxy. We find that the memory effect results in a phase of accelerated, super-Chandrasekhar friction before the perturber stalls at a critical radius, R_crit, in the core (core stalling). Inside R_crit the torque flips sign, giving rise to dynamical buoyancy, which counteracts friction and causes the perturber to stall. This phenomenology is consistent with N-body simulations, but has thus far eluded proper explanation.

Total solar eclipses (TSEs) provide a unique opportunity to quantify the properties of the K-corona (electrons), F-corona (dust), and E-corona (ions) continuously from the solar surface out to a few solar radii. We apply a novel inversion method to separate emission from the K- and F-corona continua using unpolarized total brightness (tB) observations from five 0.5 nm bandpasses acquired during the 2019 July 2 TSE between 529.5 and 788.4 nm. The wavelength dependence relative to the photosphere (i.e., color) of the F-corona itself is used to infer the tB of the K- and F-corona for each line of sight. We compare our K-corona emission results with the Mauna Loa Solar Observatory (MLSO) K-Cor polarized brightness (pB) observations from the day of the eclipse, and the forward modeled K-corona intensity from the Predictive Science Inc. (PSI) magnetohydrodynamic (MHD) model prediction. Our results are generally consistent with previous work and match both the MLSO data and PSI-MHD predictions quite well, supporting the validity of our approach and of the PSI-MHD model. However, we find that the tB of the F-corona is higher than expected in the low corona, perhaps indicating that the F-corona is slightly polarized—challenging the common assumption that the F-corona is entirely unpolarized.

To determine the importance of merging galaxies to galaxy evolution, it is necessary to design classification tools that can identify the different types and stages of merging galaxies. Previously, using GADGET-3/SUNRISE simulations of merging galaxies and linear discriminant analysis (LDA), we created an accurate merging galaxy classifier based on imaging predictors. Here, we develop a complementary tool, based on stellar kinematic predictors, derived from the same simulation suite. We design mock stellar velocity and velocity dispersion maps to mimic the specifications of the Mapping Nearby Galaxies at Apache Point (MaNGA) integral field spectroscopy (IFS) survey, and utilize an LDA to create a classification, based on a linear combination of 11 kinematic predictors. The classification varies significantly with mass ratio; the major (minor) merger classifications have a mean statistical accuracy of 80% (70%), a precision of 90% (85%), and a recall of 75% (60%). The major mergers are best identified by predictors that trace global kinematic features, while the minor mergers rely on local features that trace a secondary stellar component. While the kinematic classification is less accurate than the imaging classification, the kinematic predictors are better at identifying post-coalescence mergers. A combined imaging + kinematic classification has the potential to reveal more complete merger samples from imaging and IFS surveys such as MaNGA. We note that since the suite of simulations used to train the classifier covers a limited range of galaxy properties (i.e., the galaxies are of intermediate mass, and disk-dominated), the results may not be applicable to all MaNGA galaxies.

Spectroscopic detection of narrow emission lines traces the presence of circumstellar mass distributions around massive stars exploding as core-collapse supernovae. Transient emission lines disappearing shortly after the supernova explosion suggest that the material spatial extent is compact and implies an increased mass loss shortly prior to explosion. Here, we present a systematic survey for such transient emission lines (Flash Spectroscopy) among Type II supernovae detected in the first year of the Zwicky Transient Facility survey. We find that at least six out of ten events for which a spectrum was obtained within two days of the estimated explosion time show evidence for such transient flash lines. Our measured flash event fraction (>30% at 95% confidence level) indicates that elevated mass loss is a common process occurring in massive stars that are about to explode as supernovae.

We describe and discuss remarkable infrared spectra, covering key portions of the 2–5 μm wavelength interval, of the probable OH/IR supergiant 2MASS J17470898−2829561 (2M1747), located in direction of the Sgr B molecular cloud complex within the Central Molecular Zone (CMZ) of the Galaxy. This star was originally singled out for examination based on its suitability for spectroscopy of lines of in the CMZ. Analysis of the spectra shows that 2M1747 is deeply embedded within Sgr B1, with A_V ≳ 100 mag, making it the only star within Sgr B for which infrared spectra have been obtained at present, and thereby a unique infrared probe of the dense interstellar medium within the CMZ. Despite the high extinction, spectra of 2M1747 reveal a veiled photosphere in the K band and circumstellar gas in the M band, giving clues as to its nature. Its 3.5–4.0 μm spectrum contains the strongest absorption lines of observed toward any object to date. The 4.5–4.8 μm spectrum has impressively deep and wide absorption lines of interstellar CO, most of which arise in dense gas within Sgr B1. The 3–5 μm spectrum also contains several solid-state absorption features, which are characteristic of both dense and diffuse clouds, and which raise questions about the identifications of some of these features. We discuss the nature of the star, the extinction to it, the extinction law for dust in the CMZ, and the identifications of the various solid-state features and where they are produced along this complex line of sight.

Plasmoids—magnetized quasi-circular structures formed self-consistently in reconnecting current sheets—were previously considered to be the graveyards of energetic particles. In this paper, we demonstrate the important role of plasmoids in shaping the particle energy spectrum in relativistic reconnection (i.e., with upstream magnetization σ_up ≫ 1). Using 2D particle-in-cell simulations in pair plasmas with σ_up = 10 and 100, we study a secondary particle energization process that takes place inside compressing plasmoids. We demonstrate that plasmoids grow in time, while their interiors compress, amplifying the internal magnetic field. The magnetic field felt by particles injected in an isolated plasmoid increases linearly with time, which leads to particle energization as a result of magnetic moment conservation. For particles injected with a power-law distribution function, this energization process acts in such a way that the shape of the injected power law is conserved, while producing an additional nonthermal tail f(E) ∝ E⁻³ at higher energies, followed by an exponential cutoff. The cutoff energy, which increases with time as , can greatly exceed σ_upm_ec². We analytically predict the secondary acceleration timescale and the shape of the emerging particle energy spectrum, which can be of major importance in certain astrophysical systems, such as blazar jets.

Giant, star-forming clumps are a common feature prevalent among high-redshift star-forming galaxies and play a critical role in shaping their chaotic morphologies and yet, their nature and role in galaxy evolution remains to be fully understood. A majority of the effort to study clumps has been focused at high redshifts, and local clump studies have often suffered from small sample sizes. In this work, we present an analysis of clump properties in the local universe, and for the first time, performed with a statistically significant sample. With the help of the citizen science-powered Galaxy Zoo: Hubble project, we select a sample of 92 z < 0.06 clumpy galaxies in Sloan Digital Sky Survey Stripe 82 galaxies. Within this sample, we identify 543 clumps using a contrast-based image analysis algorithm and perform photometry as well as estimate their stellar population properties. The overall properties of our z < 0.06 clump sample are comparable to the high-redshift clumps. However, contrary to the high-redshift studies, we find no evidence of a gradient in clump ages or masses as a function of their galactocentric distances. Our results challenge the inward migration scenario for clump evolution for the local universe, potentially suggesting a larger contribution of ex situ clumps and/or longer clump migration timescales.

The purpose of this paper is to study the behavior of magnetohydrodynamic (MHD) wave modes that propagate in compressible magnetic flux tubes with an elliptical cross section embedded in a magnetic environment. The dispersion relation that describes the behavior of MHD wave modes permitted in an elliptical magnetic flux tube is solved numerically. Distortion of the spatial structure of the purely real eigenmodes from the well-known circular flux tube model has been considered. It has been studied under both photospheric and coronal conditions. It has been shown that (i) solutions in the form of even Mathieu functions are more sensitive to the value of eccentricity than solutions with the form of odd Mathieu functions; (ii) if the ellipticity of the cross section of the magnetic flux tube increases, a sausage mode (m = 0) cannot be easily identified; (iii) even solutions that correspond to the fluting mode (m = 3) can be misinterpreted as a kink mode (m = 1) due to their similarities. In contrast to the fluting modes that are polarized along the major axis and strongly depend on the ellipticity of the magnetic flux tube, the kink and sausage surface modes are practically unaffected by ellipticity. Several examples of the spatial structure of the eigenmodes permitted in the pores and sunspots have been visualized. The solutions obtained in the approximation of cylindrical symmetry are in agreement with previous studies.

A small number of interplanetary coronal mass ejections (ICMEs) have been identified that contain measurable contributions from prominence plasma. In situ measurements from during these events are marked by the presence of unusually low-charge states of C, O, and Fe, representing ionization equilibrium formation temperatures of ∼10⁴–10⁵ K, consistent with prominence material observed at the Sun. We present a thorough analysis of the elemental abundances of a wide variety of heavy ions, measured by Advanced Composition Explorer/SWICS, in prominence material observed in the solar wind. We find that prominence material observed in situ tends to be more enriched in heavy ions than the surrounding ICME plasma and the fast and slow solar wind. We also find that the material is on average moderately enhanced in low first ionization potential elements compared to photospheric abundances, with values that lie between fast and slow solar wind. In rare instances, where in situ prominence material is observed to have clear, persistent, low-charge states over longer periods of time, it exhibits elemental abundances that are photospheric in nature. However, in most prominence events we see indications that the associated material contains a mixture of prominence and adjacent ICME plasma. The anomalous behavior of the elemental and ionic composition in ICMEs with and without prominence material can be used to study physical processes that occur during CME initiation and release.

We present high-resolution Magellan/MIKE spectra of 22 bright (9 < V < 13.5) metal-poor stars (−3.18 < [Fe/H] < −1.37) in three different stellar streams, the Helmi debris stream, the Helmi trail stream, and the ω Centauri progenitor stream. We augment our Helmi debris sample with results for 10 stars by Roederer et al. for a total of 32 stars. Detailed chemical abundances of light elements as well as heavy neutron-capture elements have been determined for our 22 stars. All three streams contain carbon-enhanced stars. For 13 stars, neutron-capture element lines were detectable, and they all show signatures in agreement with the scaled solar r-process pattern, albeit with a large spread of −0.5 < [Eu/Fe] < +1.3. Eight of these stars show an additional small s-process contribution superposed onto their r-process pattern. This could be discerned because of the relatively high signal-to-noise ratio of the spectra given that the stars are close by in the halo. Our results suggest that the progenitors of these streams experienced one or more r-process events early on, such as a neutron star merger or another prolific r-process source. This widely enriched these host systems before their accretion by the Milky Way. The small s-process contribution suggests the presence of asymptotic giant branch stars and associated local (inhomogeneous) enrichment as part of the ongoing chemical evolution by low-mass stars. Stars in stellar streams may thus be a promising avenue for studying the detailed history of large dwarf galaxies and their role in halo assembly with easily accessible targets for high-quality spectra of many stars.

Due to their small mass, subsolar mass black hole binaries would have to be primordial in origin instead of the result of stellar evolution. Soon after formation in the early universe, primordial black holes can form binaries after decoupling from the cosmic expansion. Alternatively, primordial black holes as dark matter could also form binaries in the late universe due to dynamical encounters and gravitational-wave braking. A significant feature for this channel is the possibility that some sources retain nonzero eccentricity in the LIGO/Virgo band. Assuming all dark matter is primordial black holes with a delta function mass distribution, 1M_⊙–1M_⊙ binaries formed in this late-universe channel can be detected by Advanced LIGO and Virgo with their design sensitivities at a rate of yr⁻¹, where 12%(3%) of events have eccentricity at a gravitational-wave frequency of 10 Hz, e^{10 Hz} ≥ 0.01(0.1), and nondetection can constrain the binary formation rate within this model. Third generation detectors would be expected to detect subsolar mass eccentric binaries as light as 0.01M_⊙ within this channel, if they accounted for the majority of the dark matter. Furthermore, we use simulated gravitational-wave data to study the ability to search for eccentric gravitational-wave signals using a quasi-circular waveform template bank with Advanced LIGO design sensitivity. For a match-filtering targeted search, assuming binaries with a delta function mass of 0.1(1)M_⊙ and the eccentricity distribution derived from this late-universe formation channel, 41%(6%) of the signals would be missed compared to the ideal detection rate due to the mismatch in the gravitational-wave signal from eccentricity.

In 2019 July, the IceCube experiment detected a high-energy neutrino from the direction of the powerful blazar PKS 1502+106. We perform multiwavelength and multimessenger modeling of this source, using a fully selfconsistent one-zone model that includes the contribution of external radiation fields typical of flat-spectrum radio quasars. We identify three different activity states of the blazar: one quiescent state and two flaring states with hard and soft gamma-ray spectra. We find two hadronic models that can describe the multiwavelength emission during all three states: a leptohadronic model with a contribution from photohadronic processes to X-rays and gamma-rays, and a proton synchrotron model, where the emission from keV to 10 GeV comes from proton synchrotron radiation. Both models predict a substantial neutrino flux that is correlated with the gamma-ray and soft X-ray fluxes. Our results are compatible with the detection of a neutrino during the quiescent state, based on event rate statistics. We conclude that the soft X-ray spectra observed during bright flares strongly suggest a hadronic contribution, which can be interpreted as additional evidence for cosmic-ray acceleration in the source independently of neutrino observations. We find that more arguments can be made in favor of the leptohadronic model vis-a-vis the proton synchrotron scenario, such as a lower energetic demand during the quiescent state. However, the same leptohadronic model would be disfavored for flaring states of PKS 1502+106 if no IceCube events were found from the direction of the source before 2010, which would require an archival search.

The growth of galaxy clusters is energetic and may trigger and/or quench star formation and black hole activity. The ENISALA⁴project is a collection of multiwavelength observations aimed at understanding how large-scale structure drives galaxy and black hole evolution. Here, we introduce optical spectroscopy of over 800 Hα emission-line galaxies, selected in 14 z ∼ 0.15–0.31 galaxy clusters, spanning a range of masses and dynamical states. We investigate the nature of the emission lines in relation to the host galaxy properties, its location within the cluster, and the properties of the parent cluster. We uncover remarkable differences between mergers and relaxed clusters. The majority of Hα emission-line galaxies in merging cluster fields are located within 3 Mpc of their center. A large fraction of these line emitters in merging clusters are powered by star formation irrespective of cluster-centric radius, while the rest are powered by active galactic nuclei (AGNs). Star-forming galaxies are rare within 3 Mpc of relaxed clusters and AGNs are most abundant at their outskirts (∼1.5–3 Mpc). We discover a population of star-forming galaxies with large equivalent widths and blue UV–optical colors found exclusively in the merging clusters in our sample. The widespread emission-line activity in merging clusters is likely supported by triggered activity in recently accreted, gas-rich galaxies. By contrast, our observations for relaxed clusters match established models in which black hole activity is enhanced at the virial radius and star formation is quenched within the infall region. We conclude that emission-line galaxies experience distinct evolutionary paths in merging and relaxed clusters.

Besides the spirals induced by the Lindblad resonances, planets can generate a family of tightly wound spirals through buoyancy resonances. The excitation of buoyancy resonances depends on the thermal relaxation timescale of the gas. By computing timescales of various processes associated with thermal relaxation, namely, radiation, diffusion, and gas–dust collision, we show that the thermal relaxation in protoplanetary disks’ surface layers (Z/R ≳ 0.1) and outer disks (R ≳ 100 au) is limited by infrequent gas–dust collisions. The use of the isothermal equation of state or rapid cooling, common in protoplanetary disk simulations, is therefore not justified. Using three-dimensional hydrodynamic simulations, we show that the collision-limited slow thermal relaxation provides favorable conditions for buoyancy resonances to develop. Buoyancy resonances produce predominantly vertical motions, whose magnitude at the ¹²CO emission surface is of the order of 100 m s⁻¹ for Jovian-mass planets, sufficiently large to detect using molecular line observations with ALMA. We generate synthetic observations and describe characteristic features of buoyancy resonances in Keplerian-subtracted moment maps and velocity channel maps. Based on the morphology and magnitude of the perturbation, we propose that the tightly wound spirals observed in TW Hya could be driven by a (sub-)Jovian-mass planet at 90 au. We discuss how non-Keplerian motions driven by buoyancy resonances can be distinguished from those driven by other origins. We argue that observations of multiple lines tracing different heights, with sufficiently high spatial/spectral resolution and sensitivity to separate the emission arising from the near and far sides of the disk, will help constrain the origin of non-Keplerian motions.

We revisit the question of what mechanism is responsible for the spins of halos of dark matter. The answer to this question is of high importance for modeling galaxy intrinsic alignment, which can potentially contaminate current and future lensing data. In particular, we show that when the dark matter halos pass nearly each other in dense environments—namely halo assemblies—they swing and spin each other via exerting mutual tidal torques. We show that this has a significant contribution to the spin of dark matter halos comparable to that calculated by the so-called tidal torque theory. We use the results of the state-of-the-art simulation of Illutris to check the prediction of this theory against the simulation data.

We reported the first glitch of PSR B2021+51 detected with the Shanghai Tian Ma Radio Telescope that occurred around MJD 58289.1 (2018 June 20). The spin frequency and its derivative suffered a jump of about 7.04 × 10⁻¹⁰ Hz and 2.6 × 10⁻¹⁸ s⁻², respectively. The pulse width at the 10% intensity level (W₁₀) of the mean pulse profile significantly decreased right before the glitch and then increased after the glitch, accompanied by an enhancement in the trailing peak of the normalized profile. The mean pulse profiles were decomposed into three Gaussian components, i.e., C1, C2, and C3 from left to right by peak phase. The width of C3 (w₃) had an apparent maximum right before the glitch and its weighted mean value became slightly smaller than the pre-glitch value after the glitch, while the weighted mean values of both w₁ and w₂ had no apparent change. The significant decrease of the mean pulse width right before the glitch is mainly due to the movement of C3 toward C1 and C2, while the increase of the mean pulse width after the glitch results from the movement of C2 toward C3. Changes in the widths of mean pulse profiles and in the widths and phase separations of Gaussian components associated with the glitch convinced us of a connection between the changes in characteristics of the pulsar emission zone and the glitch activity. We discuss a possible interpretation invoking the movements of flux tubes in the emission zone.

The ²³Na(α, p)²⁶Mg reaction has been identified as having a significant impact on the nucleosynthesis of several nuclei between Ne and Ti in Type Ia supernovae, and of ²³Na and ²⁶Al in massive stars. The reaction has been subjected to renewed experimental interest recently, motivated by high uncertainties in early experimental data and in the statistical Hauser-Feshbach models used in reaction rate compilations. Early experiments were affected by target deterioration issues and unquantifiable uncertainties. Three new independent measurements instead are utilizing inverse kinematics and Rutherford scattering monitoring to resolve this. In this work we present directly measured angular distributions of the emitted protons to eliminate a discrepancy in the assumptions made in the recent reaction rate measurements, which results in cross sections differing by a factor of 3. We derive a new combined experimental reaction rate for the ²³Na(α, p)²⁶Mg reaction with a total uncertainty of 30% at relevant temperatures. Using our new ²³Na(α, p)²⁶Mg rate, the ²⁶Al and ²³Na production uncertainty is reduced to within 8%. In comparison, using the factor of 10 uncertainty previously recommended by the rate compilation STARLIB, ²⁶Al and ²³Na production was changing by more than a factor of 2. In Type Ia supernova conditions, the impact on production of ²³Na is constrained to within 15%.

Radio active galactic nuclei (RAGNs) are mainly found in dense structures (i.e., clusters/groups) at redshifts of z < 2 and are commonly used to detect protoclusters at higher redshift. Here, we attempt to study the host and environmental properties of two relatively faint (L_{1.4 GHz} ∼ 10²⁵ W Hz⁻¹) RAGNs in a known protocluster at z = 3.3 in the PCl J0227-0421 field, detected using the latest radio observation obtained as part of the Observations of Redshift Evolution in Large-Scale Environments (ORELSE) survey. Using new spectroscopic observations obtained from the Keck/Multi-Object Spectrometer for Infra-Red Exploration as part of the Charting Cluster Construction with the VIMOS Ultra-Deep Survey (VUDS) and ORELSE (C3VO) survey and previous spectroscopic data obtained as part of the VIMOS-Very Large Telescope Deep Survey and VUDS, we revise the three-dimensional overdensity field around this protocluster. The protocluster is embedded in a large-scale overdensity protostructure. This protostructure has an estimated total mass of ∼2.6 × 10¹⁵M_⊙ and contains several overdensity peaks. Both RAGNs are hosted by very bright and massive galaxies, while their hosts show extreme differences in color, indicating that they are of different ages and are in different evolutionary stages. Furthermore, we find that they are not in the most locally dense parts of the protostructure, but are fairly close to the centers of their parent overdensity peaks. We propose a scenario where merging might have already happened in both cases, which lowered the local density of their surrounding area and boosted their stellar mass. This work is the first time that two RAGNs at low luminosity have been found and studied within a high-redshift protostructure.

We present new observational evidence for one of the most important three wave interactions, called the electrostatic decay instability (ESD) in the source regions of complex solar type III radio bursts (L is the electron beam-excited Langmuir wave, and and S are the ESD excited daughter Langmuir and ion sound waves, respectively). The STEREO in situ wave observations in the source regions of complex type III bursts show that Langmuir waves often occur as one-dimensional magnetic field aligned beat-type wave packets, with peak intensities well in excess of the threshold for excitation of ESD, and with spectra containing (a) two closely spaced narrow peaks (L and ) corresponding probably to the beating modes responsible for the beat patterns at frequencies very close to the local electron plasma frequency, f_pe, and (b) narrow peaks at ion sound frequencies, f_S, which are very close to beat frequencies. Using the FFT and higher order spectral techniques, we show that the frequency, wavevector and phase resonance conditions required for excitation of ESD are well satisfied for these wave packets, and the speeds of electron beams derived from the resonance conditions agree reasonably well with those derived from the drift rates of the associated type III events. We also show that the merging of (L) and () most probably is the excitation mechanism of the second harmonic radio emission of these type III bursts.

We present new empirical constraints on the evolution of , the cosmological mass density of molecular hydrogen, back to z ≈ 2.5. We employ a statistical approach measuring the average observed 850 μm flux density of near-infrared selected galaxies as a function of redshift. The redshift range considered corresponds to a span where the 850 μm band probes the Rayleigh–Jeans tail of thermal dust emission in the rest frame, and can therefore be used as an estimate of the mass of the interstellar medium. Our sample comprises of ≈150,000 galaxies in the UK InfraRed Telescope Infrared Deep Sky Survey Ultra-Deep Survey field with near-infrared magnitudes K_AB ≤ 25 mag and photometric redshifts with corresponding probability distribution functions derived from deep 12-band photometry. With a sample approximately 2 orders of magnitude larger than in previous works we significantly reduce statistical uncertainties on to z ≈ 2.5. Our measurements are in broad agreement with recent direct estimates from blank field molecular gas surveys, finding that the epoch of molecular gas coincides with the peak epoch of star formation with at z ≈ 2. We demonstrate that can be broadly modeled by inverting the star formation rate (SFR) density with a fixed or weakly evolving star formation efficiency. This “constant efficiency” model shows a similar evolution to our statistically derived , indicating that the dominant factor driving the peak star formation history at z ≈ 2 is a larger supply of molecular gas in galaxies rather than a significant evolution of the SFR efficiency within individual galaxies.

We compare maps of scattering polarization signals obtained from three-dimensional (3D) radiation transfer calculations in a magnetoconvection model of the solar atmosphere using formal solvers based on the “short characteristics” (SC) and the “long characteristics” (LC) methods. The SC method requires less computational work, but it is known to introduce spatial blurring in the emergent radiation for inclined lines of sight. For polarized radiation this effect is generally more severe due to it being a signed quantity and to the sensitivity of the scattering polarization to the model’s inhomogeneities. We study the differences in the polarization signals of the emergent spectral line radiation calculated with such formal solvers. We take as a case study already published results of the scattering polarization in the Sr i 4607 Å line obtained with the SC method, demonstrating that in high-resolution grids it is accurate enough for that type of study. In general, the LC method is the preferred one for accurate calculations of the emergent radiation, which is the reason why it is now one of the options in the public version of the 3D radiative transfer code PORTA.

The presence of planets or substellar objects still embedded in their native protoplanetary disks is indirectly suggested by disk substructures like gaps, cavities, and spirals. However, these companions are rarely detected. We present Very Large Telescope/NACO high-contrast images in the J, H, K_S, and bands of the young star [BHB2007]-1 probing the inclined disk in scattered light and revealing the probable presence of a companion. The point source is detected in the band in spatial correspondence with complementary Very Large Array observations. This object is constrained to have a mass in the range of 37–47 M_Jup and is located at 50 au from the central star, inside the 70 au-large disk cavity recently imaged by the Atacama Large Millimeter/submillimeter Array (ALMA), that is absent from our NACO data (down to 20 au). This mass range is compatible with the upper end derived from the size of the ALMA cavity. The NIR disk brightness is highly asymmetric around the minor axis, with the southern side 5.5 times brighter than the northern side. The constant amount of asymmetry across all wavelengths suggests that it is due to a shadow cast by a misaligned inner disk. The massive companion that we detect could, in principle, explain the possible disk misalignment, as well as the different cavity sizes inferred by the NACO and ALMA observations. The confirmation and characterization of the companion is entrusted to future observations.

Two magnetic braking models are implemented in MESA for use in the MIST stellar model grids. Stars less than about 1.3 solar masses are observed to spin down over time through interaction with their magnetized stellar winds (i.e., magnetic braking). This is the basis for gyrochronology and is fundamental to the evolution of lower-mass stars. The detailed physics behind magnetic braking are uncertain, as are 1D stellar evolution models. Thus, we calibrate our models and compare to data from open clusters. Each braking model tested here is capable of reproducing aspects of the data, with important distinctions; neither fully accounts for the observations. The Matt et al. prescription matches the slowly rotating stars observed in open clusters but tends to overestimate the presence of rapidly rotating stars. The Garraffo et al. prescription often produces too much angular momentum loss to accurately match the observed slow sequence for lower-mass stars but reproduces the bimodal nature of slowly and rapidly rotating stars observed in open clusters fairly well. Our models additionally do not reproduce the observed solar lithium depletion, corroborating previous findings that effects other than rotation may be important. We find additional evidence that some level of mass dependency may be missing in these braking models to match the rotation periods observed in clusters older than 1 Gyr better.

We present ¹²CO (J = 2–1), ¹³CO (J = 2–1), and C¹⁸O (J = 2–1) observations toward GMC-8, one of the most massive giant molecular clouds (GMCs) in M33 using ALMA with an angular resolution of 044 × 027 (∼2 pc × 1 pc). The earlier studies revealed that its high-mass star formation is inactive in spite of a sufficient molecular reservoir with a total mass of ∼10⁶M_⊙. The high-angular resolution data enable us to resolve this peculiar source down to a molecular clump scale. One of the GMC’s remarkable features is that a round-shaped gas structure (the “Main cloud”) extends over the ∼50 pc scale, which is quite different from the other two active star-forming GMCs dominated by remarkable filaments/shells obtained by our series of studies in M33. The fraction of the relatively dense gas traced by the ¹³CO data with respect to the total molecular mass is only ∼2%, suggesting that their spatial structure and the density are not well developed to reach an active star formation. The CO velocity analysis shows that the GMC is composed of a single component as a whole, but we found some local velocity fluctuations in the Main cloud and extra blueshifted components at the outer regions. Comparing the CO with previously published large-scale H i data, we suggest that an external atomic gas flow supplied a sufficient amount of material to grow the GMC up to ∼10⁶M_⊙.

The Atacama Large Millimeter/submillimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Band 6 scan (212–272 GHz) covers potential [C ii] emission in galaxies at 6 ≤ z ≤ 8 throughout a 2.9 arcmin² area. By selecting on known Lyα emitters (LAEs) and photometric dropout galaxies in the field, we perform targeted searches down to a 5σ [C ii] luminosity depth L_{[C II]} ∼ 2.0 × 10⁸L_⊙, corresponding roughly to star formation rates (SFRs) of 10–20 M_⊙ yr⁻¹ when applying a locally calibrated conversion for star-forming galaxies, yielding zero detections. While the majority of galaxies in this sample are characterized by lower SFRs, the resulting upper limits on [C ii] luminosity in these sources are consistent with the current literature sample of targeted ALMA observations of z = 6–7 LAEs and Lyman-break galaxies (LBGs), as well as the locally calibrated relations between L_{[C ii]} and SFR—with the exception of a single [C ii]-deficient, UV-luminous LBG. We also perform a blind search for [C ii]-bright galaxies that may have been missed by optical selections, resulting in an upper limit on the cumulative number density of [C ii] sources with L_{[C II]} > 2.0 × 10⁸L_⊙ (5σ) to be less than 1.8 × 10⁻⁴ Mpc⁻³ (90% confidence level). At this luminosity depth and volume coverage, we present an observed evolution of the [C ii] luminosity function from z = 6–8 to z ∼ 0 by comparing the ASPECS measurement to literature results at lower redshift.

We present simultaneous measurements of emission from dust continuum at 230 GHz and the J = 2–1 ¹²CO, ¹³CO, and C¹⁸O isotopologues at ∼15 pc resolution from individual giant molecular clouds (GMCs) in the Andromeda galaxy (M31). These observations were obtained in an ongoing survey of this galaxy being conducted with the Submillimeter Array. Initial results describing the continuum and ¹²CO emission have been published earlier. Here, we primarily analyze the observations of ¹³CO and C¹⁸O emission and compare them to the measurements of dust continuum and ¹²CO emission. We also report additional dust continuum and CO measurements from newly added GMCs to the M31 sample. We detect spatially resolved ¹³CO emission with high signal-to-noise ratios in 31 objects. We find the extent of the ¹³CO emission to be nearly comparable to that of ¹²CO, typically covering 75% of the area of the ¹²CO emission. We derive ¹³CO and C¹⁸O abundances of 2.9 × 10⁻⁶ and 4.4 × 10⁻⁷ relative to H₂, respectively, by comparison with hydrogen column densities of the same regions derived from the dust continuum observations assuming a Milky Way gas-to-dust ratio. We find the isotopic abundance ratio [¹³CO]/[C¹⁸O] = 6.7 ± 2.9 to be consistent with the Milky Way value (8.1). Finally, we derive the mass-to-light conversion factors for all three CO species to be α₁₂ = 8.7 ± 3.9, α₁₃ = 48.9 ± 20.4, and M_⊙ (K km s⁻¹ pc²)⁻¹ for the J = 2–1 transitions of ¹²CO, ¹³CO, and C¹⁸O, respectively.

The prediction of the equation of state of hot, dense nuclear matter is one of the most complicated and interesting problems in nuclear astrophysics. At the same time, knowledge of it is the basic ingredient for some of the most interesting studies. In the present work, we concentrate our study on the construction of the equation of state of hot, dense nuclear matter, related mainly to the interior of the neutron star. We employ a theoretical nuclear model, which includes momentum-dependent interaction among the nucleons, along with state-of-the-art microscopic calculations. Thermal effects are introduced in a self-consistent way, and a set of isothermal and isentropic equations of state are predicted. The predicted equations of state are used in order to acquire and extend the knowledge of the thermal effect on both nonrotating and rapidly rotating with the Kepler frequency neutron stars. The simultaneous study of thermal and rotation effects provides useful information on some of the most important quantities, including the mass (gravitational and baryon) and radius, the Kepler frequency and Kerr parameter, the moment of inertia, etc. These quantities are directly related to studies of protoneutron stars and mainly the hot and rapidly rotating remnant of a binary neutron star merger. Data from the late observations of binary neutron star mergers and the present study may offer useful tools for investigation and help in providing possible constraints on the equation of state of nuclear matter.

We study the spectral diversity of Type Ia supernovae (SNe Ia) at maximum light using high signal-to-noise spectrophotometry of 173 SNe Ia from the Nearby Supernova Factory. We decompose the diversity of these spectra into different extrinsic and intrinsic components, and we construct a nonlinear parameterization of the intrinsic diversity of SNe Ia that preserves pairings of “twin” SNe Ia. We call this parameterization the “Twins Embedding.” Our methodology naturally handles highly nonlinear variability in spectra, such as changes in the photosphere expansion velocity, and uses the full spectrum rather than being limited to specific spectral line strengths, ratios, or velocities. We find that the time evolution of SNe Ia near maximum light is remarkably similar, with 84.6% of the variance in common to all SNe Ia. After correcting for brightness and color, the intrinsic variability of SNe Ia is mostly restricted to specific spectral lines, and we find intrinsic dispersions as low as ∼0.02 mag between 6600 and 7200 Å. With a nonlinear three-dimensional model plus one dimension for color, we can explain 89.2% of the intrinsic diversity in our sample of SNe Ia, which includes several different kinds of “peculiar” SNe Ia. A linear model requires seven dimensions to explain a comparable fraction of the intrinsic diversity. We show how a wide range of previously established indicators of diversity in SNe Ia can be recovered from the Twins Embedding. In a companion article, we discuss how these results can be applied to the standardization of SNe Ia for cosmology.

We show how spectra of Type Ia supernovae (SNe Ia) at maximum light can be used to improve cosmological distance estimates. In a companion article, we used manifold learning to build a three-dimensional parameterization of the intrinsic diversity of SNe Ia at maximum light that we call the “Twins Embedding.” In this article, we discuss how the Twins Embedding can be used to improve the standardization of SNe Ia. With a single spectrophotometrically calibrated spectrum near maximum light, we can standardize our sample of SNe Ia with an rms of 0.101 ± 0.007 mag, which corresponds to 0.084 ± 0.009 mag if peculiar velocity contributions are removed and to 0.073 ± 0.008 mag if a larger reference sample were obtained. Our techniques can standardize the full range of SNe Ia, including those typically labeled as peculiar and often rejected from other analyses. We find that traditional light-curve width + color standardization such as SALT2 is not sufficient. The Twins Embedding identifies a subset of SNe Ia, including, but not limited to, 91T-like SNe Ia whose SALT2 distance estimates are biased by 0.229 ± 0.045 mag. Standardization using the Twins Embedding also significantly decreases host-galaxy correlations. We recover a host mass step of 0.040 ± 0.020 mag compared to 0.092 ± 0.026 mag for SALT2 standardization on the same sample of SNe Ia. These biases in traditional standardization methods could significantly impact future cosmology analyses if not properly taken into account.

We combine asteroseismology, optical high-resolution spectroscopy, and kinematic analysis for 26 halo red giant branch stars in the Kepler field in the range of −2.5 < [Fe/H] < −0.6. After applying theoretically motivated corrections to the seismic scaling relations, we obtain an average mass of 0.97 ± 0.03 M_⊙ for our sample of halo stars. Although this maps into an age of ∼7 Gyr, significantly younger than independent age estimates of the Milky Way stellar halo, we considered this apparently young age to be due to the overestimation of stellar mass in the scaling relations. There is no significant mass dispersion among lower red giant branch stars (log g > 2), which constrains the relative age dispersion to <18%, corresponding to <2 Gyr. The precise chemical abundances allow us to separate the stars with [Fe/H] > −1.7 into two [Mg/Fe] groups. While the [α/Fe] and [Eu/Mg] ratios are different between the two subsamples, [s/Eu], where s stands for Ba, La, Ce, and Nd, does not show a significant difference. These abundance ratios suggest that the chemical evolution of the low-Mg population is contributed by Type Ia supernovae, but not by low- to intermediate-mass asymptotic giant branch stars, providing a constraint on its star formation timescale as 100 Myr < τ < 300 Myr. We also do not detect any significant mass difference between the two [Mg/Fe] groups, thus suggesting that their formation epochs are not separated by more than 1.5 Gyr.

The high-frequency radio sky has historically remained largely unexplored due to the typical faintness of sources in this regime, and the modest survey speed compared to observations at lower frequencies. However, high-frequency radio surveys offer an invaluable tracer of high-redshift star formation, as they directly target the faint radio free–free emission. We present deep continuum observations at 34 GHz in the COSMOS and GOODS-North fields from the Karl G. Jansky Very Large Array (VLA), as part of the COLDz survey. The deep COSMOS mosaic spans down to σ = 1.3 μJy beam⁻¹, while the wider GOODS-N observations cover to σ = 5.3 μJy beam⁻¹. We detect a total of 18 galaxies at 34 GHz, of which nine show radio emission consistent with being powered by star formation; although for two sources, this is likely due to thermal emission from dust. Utilizing deep ancillary radio data at 1.4, 3, 5, and 10 GHz, we decompose the spectra of the remaining seven star-forming galaxies into their synchrotron and thermal free–free components, and find typical thermal fractions and synchrotron spectral indices comparable to those observed in local star-forming galaxies. We further determine free–free star formation rates (SFRs), and show that these are in agreement with SFRs from spectral energy distribution-fitting and the far-infrared/radio correlation. Our observations place strong constraints on the high-frequency radio emission in typical galaxies at high redshift, and provide some of the first insights into what is set to become a key area of study with future radio facilities, such as the Square Kilometer Array Phase 1 and next-generation VLA.

We present elemental abundance results for HE 2148−2039 and HE 2155−2043 based on a detailed high-resolution spectroscopic analysis. The high-resolution Subaru/High Dispersion Spectrograph (HDS) spectra used for our analysis have a resolution of R ∼ 60,000. Although limited information based on photometry and low-resolution spectroscopy is available, we present for the first time an abundance analysis based on high-resolution spectra for both objects. Our analysis shows that the two objects are extremely metal-poor with [Fe/H] < −3. Among the neutron-capture elements, abundances of only Sr and Ba could be determined in our program stars. For both the objects [Ba/Fe] is found to be <0. While strontium is underabundant in HE 2148−2039 with [Sr/Fe] ∼ −2.02, Sr is near solar in HE 2155−2043. The locations of the program stars in the absolute carbon abundance, A(C) versus [Fe/H] diagram, show that HE 2148−2039 is a carbon-enhanced metal-poor (CEMP)-no Group II object and HE 2155−2043 is a CEMP-no Group III object. Observed [Sr/Ba] ratios are characteristics of a fast rotating massive star progenitor for HE 2155−2043 and a metal-poor asymptotic giant branch (AGB) star for HE 2148−2039. The estimated [Sc/Mn] as well as [C/Cr] ratios in HE 2155−2043 show that the surface chemical composition of this object is mono-enriched. The surface chemical composition of HE 2148−2039 is also found to be mono-enriched based on the [Mg/C] ratio. With respect to their locations in the [C/N] versus T_eff diagram, HE 2148−2039 shows signatures of mixing, and HE 2155−2043 falls in the unmixed region of the [C/N] versus T_eff plot. Kinematic analysis shows that both objects belong to the Galactic halo population.

Observational evidence shows that coronal jets can hit prominences and set them in motion. The impact leads to large-amplitude oscillations (LAOs) in the prominence. In this paper, we attempt to understand this process via 2.5D MHD numerical experiments. In our model, the jets are generated in a sheared magnetic arcade above a parasitic bipolar region located in one of the footpoints of the filament channel (FC) supporting the prominence. The shear is imposed at velocities not far above the observed photospheric values; this leads to a multiple reconnection process, as obtained in previous jet models. Both a fast Alfvénic perturbation and a slower supersonic front preceding a plasma jet are issued from the reconnection site; in the later phase, a more violent (eruptive) jet is produced. The perturbation and jets run along the FC; they are partially reflected at the prominence, and partially transmitted through it. This results in a pattern of counter-streaming flows along the FC, and oscillations in the prominence. The oscillations are LAOs (i.e., with amplitudes above 10 km s⁻¹) in some areas of the prominence, both in the longitudinal and transverse directions. In some field lines, the impact is so strong that the prominence mass is brought out of the dip and down to the chromosphere along the FC. Two cases are studied, with respect to arcades at different heights above the parasitic bipolar region, leading to different heights for the region of the prominence perturbed by the jets. The obtained oscillation amplitudes and periods are in general agreement with the observations.

We use a time-dependent hydrodynamic code and a non-LTE Monte Carlo code to model disk dissipation for the Be star 66 Ophiuchi. We compiled 63 years of observations from 1957 to 2020 to encompass the complete history of the growth and subsequent dissipation of the star’s disk. Our models are constrained by new and archival photometry, spectroscopy, and polarization observations, allowing us to model the disk dissipation event. Using Markov Chain Monte Carlo methods, we find that the properties of 66 Oph are consistent with those of a standard B2Ve star. We computed a grid of 61,568 Be star disk models to constrain the density profile of the disk before dissipation using observations of the Hα line profile and spectral energy distribution. We find at the onset of dissipation the disk has a base density of 2.5 × 10⁻¹¹ g cm⁻³ with a radial power-law index of n = 2.6. Our models indicate that after 21 yr of disk dissipation 66 Oph’s outer disk remained present and bright in the radio. We find an isothermal disk with constant viscosity with an α = 0.4 and an outer disk radius of ∼115 stellar radii best reproduces the rate of 66 Oph’s disk dissipation. We determined the interstellar polarization in the direction of the star in the V band is p = 0.63 ± 0.02% with a polarization position angle of θ_IS ≈ 857 ± 07. Using the Stokes QU diagram, we find the intrinsic polarization position angle of 66 Oph’s disk is θ_int ≈ 98° ± 3°.

An accurate and systematic quantification of the state variables in the Mars upper atmosphere is important for the study of Mars volatile evolution. In this study, we perform a comprehensive analysis of the Lyα limb scans observed by MAVEN during Martian years 32–35 to quantify the Mars upper atmospheric temperature, the CO₂ density, the H density, and the H thermal escape rate. A previously uncharacterized feature of the proton aurora at Mars is reported, which affects the overall shape of the limb scans without introducing peak intensity structures. Our inversion results excluding such scans demonstrate an unambiguous and systematic consistency of the Mars exobase temperature (∼150–280 K) derived from the dayside Lyα emission with those derived from the CO₂⁺ airglow emission and those from atmospheric drag measurements, resolving a long-standing discrepancy between those different estimates. The temperature is shown to be highly predictable over different solar cycles and seasons. However, the H density and escape rate can be episodically enhanced by a factor of ∼2–3, likely driven by dust storms that coincided with the enhancements. During the three Martian years, the averaged thermal escape rates between Mars solar longitude 0°–180° and 210°–330° are 0.6 × 10⁷ cm⁻² s⁻¹ and 9.0 × 10⁷ cm⁻² s⁻¹, respectively, a factor of 15 seasonal difference, with a factor of ∼7.5 and ∼2 due to temperature and H density variations, respectively. These temperatures and thermal escape rates derived from the thermospheric Lyα emission are generally lower than those derived from the exospheric Lyα emission in the literature, indicating the presence of a hot H population in the Mars exosphere.

Using the first long-term photometry from the Transiting Exoplanet Survey Satellite, we have carried out a detailed time-resolved timing analysis of an intermediate polar TX Col. The power spectra of almost 52 days of continuous time-series data reveal the orbital period of 5.691 ± 0.006 hr, spin period of 1909.5 ± 0.2 s, and beat period of 2105.76 ± 0.25 s, which is consistent with the earlier results. We have also found the presence of quasi-periodic oscillations (QPOs) for a few days with a period of 5850–5950 s, which appears to be due to the beating of the Keplerian period with the spin period of the white dwarf. The continuous data allowed us to look thoroughly at the day-wise evolution of the system’s accretion geometry. We report here that the TX Col changes its accretion mechanism even on a timescale of one day, confirming its variable disk-overflow accretion nature. For the majority of the time, it was found to be disk-overflow system with stream-fed dominance; however, pure disk-fed and pure stream-fed accretions cannot be ruled out during the observations.

Observations proved that a distributed structure named a supra-arcade fan (SAF) exists above post-flare loops in solar eruptions. The locations of the SAF are spatially consistent with various emission sources. Termination shocks (TSs) that are often regarded as an efficient driver for particle acceleration possibly exist in the SAF. We performed the numerical simulations of solar flares based on the standard flare model to study the dynamical and thermal manifestations of the SAF, as well as the possibility of detecting TSs in extreme-ultraviolet (EUV) images. In the simulations, the SAF and TSs can be clearly identified. The motion history and temperature evolution of plasmas inside the SAF indicate that the mass of the SAF comes from the corona and the plasmas are heated in the current sheet. The height of the SAF decreases with the speed of about 64.6 km s⁻¹ when the rate of magnetic reconnection quickly increases, and then increases with a slightly lower velocity of about 50.5 km s⁻¹ after the peak of the rate of magnetic reconnection. The descent−ascent path of the SAF is due to the unbalance of the Lorentz force and the pressure force inside the magnetic loops. In synthetic EUV images, emission intensity variations in the area surrounding TSs are significant, indicating that, depending on the viewing angle, TSs could be identifiable in EUV observations. The results of numerical simulations are generally consistent with observations, helping us to better understand the characteristics of the SAF and the physical natures behind it.

High-mass gravitational-wave events in the neutron-star mass range, such as GW190425, have recently started to be detected by the LIGO/Virgo detectors. If the masses of the two binary components fall in the neutron-star mass range, such a system is typically classified as a binary neutron-star system, although the detected gravitational-wave signal may be too noisy to clearly establish a neutron-star nature of the high-mass component in the binary and rule out a black hole–neutron star system for such an event. We show that high-mass binary neutron-star mergers with a very massive neutron-star primary close to the maximum-mass limit, m₁ ≳ 2.2 M_⊙, produce fast dynamical mass ejecta from the spin-up of the primary star at merger. By simulating the merger of black hole–neutron star systems of exactly the same masses and spins, we show that these fast ejecta are entirely absent if the primary is instead a black hole. In addition, we find that both systems leave almost identical amounts of baryon mass behind, which is not immediately accreted by the black hole. This implies that both systems will likely have comparable electromagnetic afterglow emission stemming from the remnant disk. Hence, fast ejecta at merger have the potential to distinguish neutron stars from black holes in high-mass mergers, although these ejecta may be challenging to detect observationally.

Solar and stellar flares are powerful events that produce intense radiation across the electromagnetic spectrum. Multiwavelength observations are highly important for understanding the nature of flares, because different flare-related processes reveal themselves in different spectral ranges. To study the correlation between thermal and nonthermal processes in stellar flares, we have searched the databases of Kepler (optical observations) and XMM-Newton (soft X-rays) for the flares observed simultaneously with both instruments; nine distinctive flares (with energies exceeding 10³³ erg) on three stars (of K-M spectral classes) have been found. We have analyzed and compared the flare parameters in the optical and X-ray spectral ranges; we have also compared the obtained results with similar observations of solar flares. Most of the studied stellar flares released more energy in the optical range than in X-rays. In one flare, X-ray emission strongly dominated, which could be caused either by a soft spectrum of energetic electrons or by a near-limb position of this flare. The X-ray flares were typically delayed with respect to and shorter than their optical counterparts, which is partially consistent with the Neupert effect. Using the scaling laws based on the magnetic reconnection theory, we have estimated the characteristic magnetic field strengths in the stellar active regions and the sizes of these active regions as about 25–70 G and 250,000–500,000 km, respectively. The observed stellar superflares appear to be scaled-up versions of solar flares, with a similar underlying mechanism and nearly the same characteristic magnetic field values, but with much larger active region sizes.