Table of contents

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the l = m = 2 mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the l = 2, m = 1, 2 modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found, so we present 95% credible upper limits on the strain amplitudes h₀ for the single-harmonic search along with limits on the pulsars' mass quadrupole moments Q₂₂ and ellipticities ε. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437−4715 and J0711−6830, which have spin-down ratios of 0.87 and 0.57, respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars, our limits are factors of ∼100 and ∼20 more constraining than their spin-down limits, respectively. For the dual-harmonic searches, new limits are placed on the strain amplitudes C₂₁ and C₂₂. For 23 pulsars, we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

Observations with the Fermi Large Area Telescope (LAT) of the gamma-ray source 4FGL J1702.7−5655, previously classified as a candidate millisecond pulsar, show highly significant modulation at a period of 0.2438033 days (∼5.85 hr). Further examination of the folded light curve indicates the presence of narrow eclipses, suggesting that this is a redback binary system. An examination of the long-term properties of the modulation over 13 years of LAT observations indicates that the orbital modulation of the gamma rays changed from a simple eclipse before early 2013 to a broader, more easily detected quasi-sinusoidal modulation. In addition, the time of the eclipse shifts to ∼0.05 later in phase. This change in the orbital modulation properties, however, is not accompanied by a significant overall change in gamma-ray flux or spectrum. The quasi-sinusoidal component peaks ∼0.5 out of phase with the eclipse, which would indicate inferior conjunction of the compact object in the system. Swift X-ray Telescope observations reveal a possible X-ray counterpart within the LAT error ellipse. However, radio observations obtained with the Australia Telescope Compact Array do not detect a source in the region. 4FGL J1702.7−5655 appears to have changed its state in 2013, perhaps related to changes in the intrabinary shock in the system. We compare the properties of 4FGL J1702.7−5655 to those of other binary millisecond pulsars that have exhibited orbital modulation in gamma-rays.

We introduce a new method for detecting ultra-diffuse galaxies by searching for over-densities in intergalactic globular cluster populations. Our approach is based on an application of the log-Gaussian Cox process, which is a commonly used model in the spatial statistics literature but rarely used in astronomy. This method is applied to the globular cluster data obtained from the PIPER survey, a Hubble Space Telescope imaging program targeting the Perseus cluster. We successfully detect all confirmed ultra-diffuse galaxies with known globular cluster populations in the survey. We also identify a potential galaxy that has no detected diffuse stellar content. Preliminary analysis shows that it is unlikely to be merely an accidental clump of globular clusters or other objects. If confirmed, this system would be the first of its kind. Simulations are used to assess how the physical parameters of the globular cluster systems within ultra-diffuse galaxies affect their detectability using our method. We quantify the correlation of the detection probability with the total number of globular clusters in the galaxy and the anticorrelation with increasing half-number radius of the globular cluster system. The Sérsic index of the globular cluster distribution has little impact on detectability.

In this work, we compile a sample of 449 Fermi blazars with the luminosity of the broadline region, the black hole mass, the beam radio luminosity, and the jet power; obtain the beam power and the black hole spin; investigate the dividing line between BL Lacertae objects (BL Lacs) and flat-spectrum radio quasars (FSRQs) to identify the discrepancies in their dominant mechanisms; and revisit the dependence of the jet power on the disk accretion luminosity, the black hole mass, and the black hole spin. We come to the following conclusions. (1) A boundary of log (L_BLR/L_Edd) = −3.14, separating the BL Lacs and the FSRQs, is obtained from the Bayesian analysis, which is consistent with the results from the literature. We employ the boundary to divide the blazar candidates of uncertain types into candidates for BL Lacs or FSRQs, and we find five changing-look blazars at the same time. (2) A strong correlation is found between black hole mass and intrinsic γ-ray luminosity, but a weaker correlation is found between black hole mass and observed γ-ray luminosity. The latter is weakened by jet effects: it is apparently weak for BL Lacs that have disordered amplification of the Doppler factor, since their mechanism is dominated by jets, while it is moderate for FSRQs, since their mechanism is dominated by accretion processes. (3) The jets of both FSRQs and BL Lacs are likely governed by the Blandford–Znajek mechanism.

We construct a physically parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of Type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an autoencoder that is interpreted probabilistically after training using a normalizing flow. We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multistage training setup alongside our physically parameterized network, we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including the automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an rms of 0.091 ± 0.010 mag, which corresponds to 0.074 ± 0.010 mag if peculiar velocity contributions are removed. Trained models and codes are released at https://github.com/georgestein/suPAErnova.

A major obstacle to detecting and characterizing long-period, low-mass exoplanets is the intrinsic radial-velocity (RV) variability of host stars. To better understand RV variability, we estimate disk-averaged RV variations of the Sun over its magnetic cycle, from the Fe i line observed by SDO/HMI, using a physical model for rotationally modulated magnetic activity that was previously validated against HARPS-N solar observations. We estimate the unsigned magnetic flux and show that a linear fit to it reduces the rms of RV variations by 62%, i.e., a factor of 2.6. We additionally apply the FF' method, which predicts RV variations based on a star's photometric variations. At cycle maximum, we find that additional processes must be at play beyond suppression of convective blueshift and velocity imbalances resulting from brightness inhomogeneities, in agreement with recent studies of RV variations. By modeling RV variations over the magnetic cycle using a linear fit to the unsigned magnetic flux, we recover injected planets at a period of ≈300 days with RV semi-amplitudes down to 0.3 m s⁻¹. To reach 0.1 m s⁻¹, we will need to identify and model additional phenomena that are not well traced by $| {\hat{B}}_{\mathrm{obs}}|$ or FF'. This study motivates ongoing and future efforts to develop observation and analysis techniques to measure the unsigned magnetic flux at high precision in slowly rotating, relatively inactive stars like the Sun. We conclude that the unsigned magnetic flux is an excellent proxy for rotationally modulated, activity-induced RV variations, and could become key to confirming and characterizing Earth analogs.

An improved Amati correlation was constructed in ApJ 931 (2022) 50 by us recently. In this paper, we further study constraints on the ΛCDM and wCDM models from the gamma-ray bursts (GRBs) standardized with the standard and improved Amati correlations, respectively. By using the Pantheon Type Ia supernova sample to calibrate the latest A220 GRB data set, the GRB Hubble diagram is obtained model-independently. We find that at the high-redshift region (z > 1.4) the GRB distance modulus from the improved Amati correlation is larger apparently than that from the standard Amati one. The GRB data from the standard Amati correlation only give a lower bound limit on the present matter density parameter Ω_m0, while the GRBs from the improved Amati correlation constrain the Ω_m0 with the 68% confidence level to be ${0.308}_{-0.230}^{+0.066}$ and ${0.307}_{-0.290}^{+0.057}$ in the ΛCDM and wCDM models, respectively, which are very consistent with those given by other current popular observational data including baryon acoustic oscillation, cosmic microwave background (CMB) radiation, and so on. Once the H(z) data are added in our analysis, the constraint on the Hubble constant H₀ can be achieved. We find that two different correlations provide slightly different H₀ results but the marginalized mean values seem to be close to that from the Planck 2018 CMB radiation observations.

We investigate various dynamic processes including magnetic reconnection, chromospheric evaporation, and coronal rain draining in two limb solar flares through imaging and spectroscopic observations from the Interface Region Imaging Spectrograph (IRIS) and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory. In the early phase of the flares, a bright and dense loop-top structure with a cusp-like shape can be seen in multiwavelength images, which is cospatial with the hard X-ray 25–50 keV emission. In particular, intermittent magnetic reconnection downflows are detected in the time–space maps of AIA 304 Å. The reconnection downflows are manifested as redshifts on one half of the loops and blueshifts on the other half in the IRIS Si iv 1393.76 Å line due to a projection effect. The Si iv profiles exhibit complex features (say, multipeak) with a relatively larger width at the loop-top region. During the impulsive phase, chromospheric evaporation is observed in both AIA images and the IRIS Fe xxi 1354.08 Å line. Upward motions can be seen from AIA 131 Å images. The Fe xxi line is significantly enhanced and shows a good Gaussian shape. In the gradual phase, warm rains are observed as downward moving plasmas in AIA 304 Å images. Both the Si iv and Fe xxi lines show a relatively symmetric shape with a larger width around the loop top. These results provide observational evidence for various dynamic processes involved in the energy release process of solar flares and are crucial to the understanding of this process.

The evolution and the stability of mass transfer of CO+He white dwarf (WD) binaries are not well understood. Observationally they may emerge as AM CVn binaries and are important gravitational wave (GW) emitters. In this work, we have modeled the evolution of double WD binaries with accretor masses of 0.50–1.30 M_⊙ and donor masses of 0.17–0.45 M_⊙ using the detailed stellar evolution code mesa. We find that the evolution of binaries with same donor masses but different accretor masses is very similar and binaries with same accretor masses but larger He donor masses have larger maximum mass transfer rates and smaller minimum orbital periods. We also demonstrate that the GW signal from AM CVn binaries can be detected by spaceborne GW observatories, such as LISA and TianQin. There is a linear relation between the donor mass and gravitational wave frequency during the mass transfer phase. In our calculation, all binaries can have dynamically stable mass transfer, which is very different from previous studies. The threshold donor mass of Eddington-limited mass transfer for a given accretor WD mass is lower than previous studies. Assuming that a binary may enter a common envelope if the mass transfer rate exceeds the maximum stable burning rate of He, we provide a new criterion for double WDs surviving mass transfer, which is below the threshold of the Eddington limit. Finally, we find that some systems with oxygen–neon (ONe) WDs in our calculation may evolve into detached binaries consisting of neutron stars and extremely low-mass He WDs, and further ultracompact X-ray binaries.

One month after launching the Gravitational Wave High-energy Electromagnetic Counterpart All-sky Monitor, a bright thermonuclear X-ray burst from 4U 0614+09 was observed on 2021 January 24. We report the time-resolved spectroscopy of the burst and a hint of burst oscillation at 413 Hz with a fractional amplitude ∼2.0% (rms). This coincides with the burst oscillation previously discovered with Swift/Burst Alert Telescope (Strohmayer et al. 2008), and therefore supports the spin frequency of this source. This burst is a bright one in the normal bursts detected from 4U 0614+09, which leads to an upper limit of distance estimation of 3.1 kpc. The folded light curve during the burst oscillation shows a sinusoidal structure, which is consistent with previous observations.

We investigate the effect of dark stars (DSs) on the reionization history of the universe, and the interplay between them and feedback due to Lyman–Werner (LW) radiation in reducing the cosmic microwave background (CMB) optical depth to a value within the τ = 0.054 ± 0.007 range measured by Planck. We use a semianalytic approach to evaluate reionization histories and CMB optical depths, which includes Population II stars in atomic cooling halos and Population III stars in minihalos with LW feedback, preceded by a DS phase. We show that while LW feedback by itself can reduce the integrated optical depth to the last scattering surface to ∼0.05 only if the Population III star formation efficiency is less than ∼0.2%, the inclusion of a population of DSs can naturally lead to the measured CMB optical depth for much larger Population III star formation efficiencies ≳1%.

NGC 5033 is an intriguing Seyfert galaxy because its subclassification may change with time, and because optical and submillimeter observations find that the massive black hole does not sit at the dynamical center of the galaxy, pointing to a past merger. We obtained a new optical spectrum of NGC 5033 using the 200'' Hale telescope at Palomar that clearly reveals a broad Hβ line (FWHM = 5400 ± 300 km s⁻¹). This signals a clear view of the optical broad line region and requires Seyfert-1.5 designation. Some spectra obtained in the past suggest a Seyfert-1.9 classification, potentially signaling a variable or "changing-look" geometry. Our analysis of a 2019 Chandra spectrum of the massive black hole reveals very little obscuration, also suggesting a clean view of the central engine. However, the narrow Fe Kα emission line is measured to have an equivalent width of EW$\,=\,{460}_{-90}^{+100}$ eV. This value is extremely high compared to typical values in unobscured active galactic nuclei . Indeed, the line is persistently strong in NGC 5033: the line equivalent width in a 2002 XMM-Newton snapshot is EW$\,=\,{250}_{-40}^{+40}$ eV, similar to the EW$\,=\,{290}_{-100}^{+100}$ eV equivalent width measured using ASCA in 1999. These results can likely be explained through a combination of an unusually high covering factor for reflection, and fluxes that are seen out of phase owing to light travel times. We examine the possibility that NGC 5033 may strengthen evidence for the X-ray Baldwin effect.

We describe a new automated technique for active region emergence in coronal magnetic field models, based on the inversion of the electric field locally from a single line-of-sight magnetogram for each region. The technique preserves the arbitrary shapes of magnetic field distribution associated with individual active regions and incorporates emerging magnetic helicity (twist) in a parametrized manner through a noninductive electric field component. We test the technique with global magnetofrictional simulations of the coronal magnetic field during Solar Cycle 24 Maximum from 2011 June 1 to 2011 December 31. The active regions are determined in a fully automated and objective way using Spaceweather HMI Active Region Patch (SHARP) data. Our primary aim is to constrain two free parameters in the emergence algorithm: the duration of emergence and the twist parameter for each individual active region. While the duration has a limited effect on the resulting coronal magnetic field, changing the sign and amplitude of the twist parameters profoundly influences the amount of nonpotentiality generated in the global coronal magnetic field. We explore the possibility of constraining both the magnitude and sign of the twist parameter using estimates of the current helicity derived from vector magnetograms and supplied in the SHARP metadata for each region. Using the observed sign of twist for each region reduces the overall nonpotentiality in the corona, highlighting the importance of scatter in the emerging active region helicities.

M54 is a prototype for a globular cluster embedded in a dark matter halo. Gaia Early Data Release 3 photometry and proper motions separate the old, metal-poor stars from the more metal-rich and younger dwarf galaxy stars. The metal-poor stars dominate the inner 50 pc, with a velocity dispersion profile that declines to a minimum around 30 pc then rises back to nearly the central velocity dispersion, as expected for a globular cluster at the center of a cold dark matter (CDM) cosmology dark matter halo. The Jeans equation mass analysis of the three separate stellar populations gives consistent masses that rise approximately linearly with radius to 1 kpc. These data are compatible with an infalling CDM dark matter halo reduced to ≃3 × 10⁸M_⊙ at the 50 kpc apocenter 2.3 Gyr ago, with a central globular cluster surrounded by the remnant of a dwarf galaxy. Tides gradually remove material beyond 1 kpc but have little effect on the stars and dark matter within 300 pc of the center. M54 appears to be a "transitional" system between globular clusters with and without local dark halos whose evolution within the galaxy depends on the time of accretion and orbital pericenter.

We present the third discovery from the COol Companions ON Ultrawide orbiTS (COCONUTS) program, the COCONUTS-3 system, composed of the young M5 primary star UCAC4 374−046899 and the very red L6 dwarf WISEA J081322.19−152203.2. These two objects have a projected separation of $61^{\prime\prime}$ (1891 au) and are physically associated given their common proper motions and estimated distances. The primary star, COCONUTS-3A, has a mass of 0.123 ± 0.006 M_⊙, and we estimate its age as 100 Myr to 1 Gyr based on its stellar activity (via Hα and X-ray emission), kinematics, and spectrophotometric properties. We derive its bulk metallicity as 0.21 ± 0.07 dex using empirical calibrations established by older and higher-gravity M dwarfs and find that this [Fe/H] could be slightly underestimated according to PHOENIX models given COCONUTS-3A's younger age. The companion, COCONUTS-3B, has a near-infrared spectral type of L6 ± 1 int-g, and we infer physical properties of T_eff = ${1362}_{-73}^{+48}$ K, $\mathrm{log}(g)$ = ${4.96}_{-0.34}^{+0.15}$ dex, $R\,=\,{1.03}_{-0.06}^{+0.12}$R_Jup, and $M\,=\,{39}_{-18}^{+11}$M_Jup using its bolometric luminosity, its host star's age, and hot-start evolution models. We construct cloudy atmospheric model spectra at the evolution-based physical parameters and compare them to COCONUTS-3B's spectrophotometry. We find that this companion possesses ample condensate clouds in its photosphere (f_sed = 1) with the data–model discrepancies likely due to the models using an older version of the opacity database. Compared to field-age L6 dwarfs, COCONUTS-3B has fainter absolute magnitudes and a 120 K cooler T_eff. Also, the J − K color of this companion is among the reddest for ultracool benchmarks with ages older than a few hundred megayears. COCONUTS-3 likely formed in the same fashion as stellar binaries given the companion-to-host mass ratio of 0.3 and represents a valuable benchmark to quantify the systematics of substellar model atmospheres.

The next generation of wide-field cosmic microwave background (CMB) surveys are uniquely poised to open a new window into time-domain astronomy in the millimeter band. Here, we explore the discovery phase space for extragalactic transients with near-term and future CMB experiments to characterize the expected population. We use existing millimeter-band light curves of known transients (gamma-ray bursts, tidal disruption events, fast blue optical transients (FBOTs), neutron star mergers) and theoretical models, in conjunction with known and estimated volumetric rates. Using Monte Carlo simulations of various CMB survey designs (area, cadence, depth, duration) we estimate the detection rates and the resulting light-curve characteristics. We find that existing and near-term surveys will find tens to hundreds of long-duration gamma-ray bursts (LGRBs), driven primarily by detections of the reverse shock emission, and including off-axis LGRBs. Next-generation experiments (CMB-S4, CMB-HD) will find tens of FBOTs in the nearby universe and will detect a few tidal disruption events. CMB-HD will additionally detect a small number of short gamma-ray bursts, where these will be discovered within the detection volume of next-generation gravitational wave experiments like the Cosmic Explorer.

Spatially extended halos of H i Lyα emission are now ubiquitously found around high-redshift star-forming galaxies. But our understanding of the nature and powering mechanisms of these halos is still hampered by the complex radiative transfer effects of the Lyα line and limited angular resolution. In this paper, we present resolved Multi Unit Spectroscopic Explorer (MUSE) observations of SGAS J122651.3+215220, a strongly lensed pair of L* galaxies at z = 2.92 embedded in a Lyα halo of L_Lyα = (6.2 ± 1.3) × 10⁴² erg s⁻¹. Globally, the system shows a line profile that is markedly asymmetric and redshifted, but its width and peak shift vary significantly across the halo. By fitting the spatially binned Lyα spectra with a collection of radiative transfer galactic wind models, we infer a mean outflow expansion velocity of ≈211 km s⁻¹, with higher values preferentially found on both sides of the system's major axis. The velocity of the outflow is validated with the blueshift of low-ionization metal absorption lines in the spectra of the central galaxies. We also identify a faint (M₁₅₀₀ ≈ −16.7) companion detected in both Lyα and the continuum, whose properties are in agreement with a predicted population of satellite galaxies that contribute to the extended Lyα emission. Finally, we briefly discuss the impact of the interaction between the central galaxies on the properties of the halo and the possibility of in situ fluorescent Lyα production.

We present a joint analysis of the power spectra of the Planck Compton y parameter map and the projected galaxy density field using the Wide Field Infrared Survey Explorer (WISE) all-sky survey. We detect the statistical correlation between WISE and Planck data (gy) with a significance of 21.8σ. We also measure the autocorrelation spectrum for the thermal Sunyaev–Zel'dovich (tSZ) (yy) and the galaxy density field maps (gg) with a significance of 150σ and 88σ, respectively. We then construct a halo model and use the measured correlations ${C}_{{\ell }}^{\mathrm{gg}}$, ${C}_{{\ell }}^{{yy}}$, and ${C}_{{\ell }}^{{\rm{g}}y}$ to constrain the tSZ mass bias $B\equiv {M}_{500}/{M}_{500}^{\mathrm{tSZ}}$. We also fit for the galaxy bias, which is included with explicit redshift and multipole dependencies as ${b}_{{\rm{g}}}{(z,{\ell })={b}_{{\rm{g}}}^{0}(1+z)}^{\alpha }{({\ell }/{{\ell }}_{0})}^{\beta }$, with ℓ₀ = 117. We obtain the constraints to be B = 1.50 ± 0.07(stat) ± 0.34(sys), i.e., 1 − b_H = 0.67 ± 0.03(stat) ± 0.16(sys) (68% confidence level) for the hydrostatic mass bias, and ${b}_{{\rm{g}}}^{0}={1.28}_{-0.04}^{+0.03}(\mathrm{stat})\pm 0.11(\mathrm{sys})$, with $\alpha ={0.20}_{-0.07}^{+0.11}(\mathrm{stat})\pm 0.10(\mathrm{sys})$ and β = 0.45 ±0.01(stat) ± 0.02(sys) for the galaxy bias. Incoming data sets from future CMB and galaxy surveys (e.g., Rubin Observatory) will allow probing the large-scale gas distribution in more detail.

The center of the nearby galaxy NGC 253 hosts a population of more than a dozen super star clusters (SSCs) that are still in the process of forming. The majority of the star formation of the burst is concentrated in these SSCs, and the starburst is powering a multiphase outflow from the galaxy. In this work, we measure the 350 GHz dust continuum emission toward the center of NGC 253 at 47 mas (0.8 pc) resolution using data from the Atacama Large Millimeter/submillimeter Array. We report the detection of 350 GHz (dust) continuum emission in the outflow for the first time, associated with the prominent South-West streamer. In this feature, the dust emission has a width of ≈8 pc, is located at the outer edge of the CO emission, and corresponds to a molecular gas mass of ∼(8–17)×10⁶M_⊙. In the starburst nucleus, we measure the resolved radial profiles, sizes, and molecular gas masses of the SSCs. Compared to previous work at the somewhat lower spatial resolution, the SSCs here break apart into smaller substructures with radii 0.4–0.7 pc. In projection, the SSCs, dust, and dense molecular gas appear to be arranged as a thin, almost linear, structure roughly 155 pc in length. The morphology and kinematics of this structure can be well explained as gas following x₂ orbits at the center of a barred potential. We constrain the morpho-kinematic arrangement of the SSCs themselves, finding that an elliptical, angular-momentum-conserving ring is a good description of both the morphology and kinematics of the SSCs.

We present the Waves in Nearby Disk galaxies Survey (WiNDS) consisting of 40 nearby low-inclination disk galaxies observed through H_α high-resolution Fabry–Perot interferometry. WiNDS consists of 12 new galaxy observations and 28 data archived observations obtained from different galaxy surveys. We derive two-dimensional line-of-sight velocity fields that are analyzed to identify the possible presence of vertical velocity flows in the galactic disks of these low-inclination late-type galaxies using velocity residual maps, derived from the subtraction of an axisymmetric rotation model from a rotational velocity map. Large and globally coherent flows in the line-of-sight velocity of nearly face-on galaxies can be associated with large vertical displacement of the disk with respect to its midplane. Our goal is to characterize how frequent vertical perturbations, such as those observed in the Milky Way, arise in the local universe. Our currently available data have allowed us to identify 20% of WiNDS galaxies with strong velocity perturbations that are consistent with vertically perturbed galactic disks.

We seek signatures of the current experimental ¹²C ${\left(\alpha ,\gamma \right)}^{16}{\rm{O}}$ reaction rate probability distribution function in the pulsation periods of carbon–oxygen white dwarf (WD) models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it forms during the evolution of low-mass stars under radiative helium-burning conditions, mitigating the impact of convective mixing processes. We make direct quantitative connections between the pulsation periods of the identified trapped g-modes in variable WD models and the current experimental ¹²C ${\left(\alpha ,\gamma \right)}^{16}{\rm{O}}$ reaction rate probability distribution function. We find an average spread in relative period shifts of ΔP/P ≃ ±2% for the identified trapped g-modes over the ±3σ uncertainty in the ¹²C ${\left(\alpha ,\gamma \right)}^{16}{\rm{O}}$ reaction rate probability distribution function—across the effective temperature range of observed DAV and DBV WDs and for different WD masses, helium shell masses, and hydrogen shell masses. The g-mode pulsation periods of observed WDs are typically given to six to seven significant figures of precision. This suggests that an astrophysical constraint on the ¹²C ${\left(\alpha ,\gamma \right)}^{16}{\rm{O}}$ reaction rate could, in principle, be extractable from the period spectrum of observed variable WDs.

Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the magnetic field, which is a key component of the launching process, is still an open question. Here we address the question of how the magnetic field required for jet launching is generated and maintained by a dynamo process. By carrying out nonideal MHD simulations (PLUTO code), we investigate how the feedback of the generated magnetic field on the mean-field dynamo affects the disk and jet properties. We find that a stronger quenching of the dynamo leads to a saturation of the magnetic field at a lower disk magnetization. Nevertheless, we find that, while applying different dynamo feedback models, the overall jet properties remain unaffected. We then investigate a feedback model that encompasses a quenching of the magnetic diffusivity. Our modeling considers a more consistent approach for mean-field dynamo modeling simulations, as the magnetic quenching of turbulence should be considered for both a turbulent dynamo and turbulent magnetic diffusivity. We find that, after the magnetic field is saturated, the Blandford–Payne mechanism can work efficiently, leading to more collimated jets, which move, however, with slower speed. We find strong intermittent periods of flaring and knot ejection for low Coriolis numbers. In particular, flux ropes are built up and advected toward the inner disk thereby cutting off the inner disk wind, leading to magnetic field reversals, reconnection and, the emergence of intermittent flares.

Due to the latest advances in technology, telescopes with significant sky coverage will produce millions of astronomical alerts per night that must be classified both rapidly and automatically. Currently, classification consists of supervised machine-learning algorithms whose performance is limited by the number of existing annotations of astronomical objects and their highly imbalanced class distributions. In this work, we propose a data augmentation methodology based on generative adversarial networks (GANs) to generate a variety of synthetic light curves from variable stars. Our novel contributions, consisting of a resampling technique and an evaluation metric, can assess the quality of generative models in unbalanced data sets and identify GAN-overfitting cases that the Fréchet inception distance does not reveal. We applied our proposed model to two data sets taken from the Catalina and Zwicky Transient Facility surveys. The classification accuracy of variable stars is improved significantly when training with synthetic data and testing with real data with respect to the case of using only real data.

Parker Solar Probe's (PSP's) unique orbital path allows us to observe the solar wind closer to the Sun than ever before. Essential to advancing our knowledge of solar wind and energetic particle formation is identifying the sources of PSP observations. We report on results for the first two PSP solar encounters derived using the Wang–Sheeley–Arge (WSA) model driven by Air Force Data Assimilative Photospheric Flux Transport (ADAPT) model maps. We derive the coronal magnetic field and the 1 R_⊙ source regions of the PSP-observed solar wind. We validate our results with the solar wind speed and magnetic polarity observed at PSP. When modeling results are very reliable, we derive time series of model-derived spacecraft separation from the heliospheric current sheet, magnetic expansion factor, coronal hole boundary distance, and photospheric field strength along the field lines estimated to be connected to the spacecraft. We present new results for Encounter 1, which show time evolution of the far-side mid-latitude coronal hole that PSP corotates with. We discuss how this evolution coincides with solar wind speed, density, and temperature observed at the spacecraft. During Encounter 2, a new active region emerges on the solar far side, making it difficult to model. We show that ADAPT-WSA output agrees well with PSP observations once this active region rotates onto the near side, allowing us to reliably estimate the solar wind sources retrospectively for most of the encounter. We close with ways in which coronal modeling enables scientific interpretation of these encounters that would otherwise not have been possible.

On 2017 September 2 MAXI J1535–571 went into outburst and peaked at ∼5 Crab in the 2–20 keV energy range. Early in the flare, the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) performed target of opportunity pointings and monitored the source as it transitioned from the hard state to the soft state. Using quasi-simultaneous observations from MAXI/GSC and INTEGRAL/SPI, we studied the temporal and spectral evolution of MAXI J1535–571 in the 2–500 keV range. Early spectra show a Comptonized spectrum and a high-energy component dominant above ∼150 keV. CompTT fits to the SPectrometer on INTEGRAL (SPI) data found electron temperatures (kT_e ) evolve from ∼31 keV to 18 keV with a tied optical depth (τ ∼ 0.85) or τ evolving from ∼1.2–0.65 with a tied kT_e (∼24 keV). To investigate the nature of the high-energy component, we performed a spectral decomposition of the 100–400 keV energy band. The CompTT flux varies significantly during the hard state while the high-energy component flux is consistent with a constant flux. This result suggests that the two components originate from different locations, which favors a jet origin interpretation for the high-energy component over a hybrid corona interpretation. Lastly, two short rebrightenings during the hard-to-soft transition are compared to similar events reported in MAXI J1820+070.

By performing two-dimensional axisymmetric general relativistic radiation magnetohydrodynamics simulations with spin parameter a* varying from −0.9 to 0.9, we investigate the dependence on the black hole spin of the energy flow from a supercritical accretion disk around a stellar mass black hole. It is found that optically and geometrically thick disks form near the equatorial plane, and a part of the disk matter is launched from the disk surface in all models. The gas ejection is mainly driven by the radiative force, but magnetic force cannot be neglected when ∣a*∣ is large. The energy outflow efficiency (total luminosity normalized by ${\dot{M}}_{\mathrm{in}}{c}^{2};$${\dot{M}}_{\mathrm{in}}$ and c are the mass-accretion rate at the event horizon and the light speed) is higher for rotating black holes than for nonrotating black holes. This is 0.7% for a* = −0.7, 0.3% for a* = 0, and 5% for a* = 0.7 for ${\dot{M}}_{\mathrm{in}}\sim 100{L}_{\mathrm{Edd}}/{c}^{2}$ (L_Edd is the Eddington luminosity). Furthermore, although the energy is mainly released by radiation when a* ∼ 0, the Poynting power increases with ∣a*∣ and exceeds the radiative luminosity for models with a* ≥ 0.5 and a* ≤ −0.7. The faster the black hole rotates, the higher the power ratio of the kinetic luminosity to the isotropic luminosity tends to be. This implies that objects with a high (low) power ratio may have rapidly (slowly) rotating black holes. Among ultraluminous X-ray sources, IC342 X-1, is a candidate with a rapidly rotating black hole.

A promising astrophysical site to produce the lighter heavy elements of the first r-process peak (Z = 38 − 47) is the moderately neutron-rich (0.4 < Y_e < 0.5) neutrino-driven ejecta of explosive environments, such as core-collapse supernovae and neutron star mergers, where the weak r-process operates. This nucleosynthesis exhibits uncertainties from the absence of experimental data from (α, xn) reactions on neutron-rich nuclei, which are currently based on statistical model estimates. In this work, we report on a new study of the nuclear reaction impact using a Monte Carlo approach and improved (α, xn) rates based on the Atomki-V2 α optical model potential. We compare our results with observations from an up-to-date list of metal-poor stars with [Fe/H] < −1.5 to find conditions of the neutrino-driven wind where the lighter heavy elements can be synthesized. We identified a list of (α, xn) reaction rates that affect key elemental ratios in different astrophysical conditions. Our study aims to motivate more nuclear physics experiments on (α, xn) reactions using the current and new generation of radioactive beam facilities and also more observational studies of metal-poor stars.

The Galactic disk exhibits complex chemical and dynamical substructure thought to be induced by the bar, spiral arms, and satellites. Here, we explore the chemical signatures of bar resonances in action and velocity space, and characterize the differences between the signatures of corotation (CR) and higher-order resonances using test particle simulations. Thanks to recent surveys, we now have large data sets containing metallicities and kinematics of stars outside the solar neighborhood. We compare the simulations to the observational data from Gaia EDR3 and LAMOST DR5 and find weak evidence for a slow bar with the "hat" moving group (250 km s⁻¹ ≲ v_ϕ ≲ 270 km s⁻¹) associated with its outer Lindblad resonance and "Hercules" (170 km s⁻¹ ≲ v_ϕ ≲ 195 km s⁻¹) with CR. While constraints from current data are limited by their spatial footprint, stars closer in azimuth than the Sun to the bar's minor axis show much stronger signatures of the bar's outer Lindblad and CR resonances in test particle simulations. Future data sets with greater azimuthal coverage, including the final Gaia data release, will allow reliable chemodynamical identification of bar resonances.

We present an analysis of 10 ks snapshot Chandra observations of 12 shocked post-starburst galaxies, which provide a window into the unresolved question of active galactic nuclei (AGN) activity in post-starburst galaxies and its role in the transition of galaxies from active star formation to quiescence. While seven of the 12 galaxies have statistically significant detections (with two more marginal detections), the brightest only obtained 10 photons. Given the wide variety of hardness ratios in this sample, we chose to pursue a forward-modeling approach to constrain the intrinsic luminosity and obscuration of these galaxies, rather than stacking. We constrain the intrinsic luminosity of obscured power laws based on the total number of counts and spectral shape, itself mostly set by the obscuration, with hardness ratios consistent with the data. We also tested thermal models. While all the galaxies have power-law models consistent with their observations, a third of the galaxies are better fit as an obscured power law and another third are better fit as thermal emission. If these post-starburst galaxies, early in their transition, contain AGNs, then these are mostly confined to lower obscuration (N_H ≤ 10²³ cm⁻²) and lower luminosity (L_{2−10 keV} ≤ 10⁴² erg s⁻¹). Two galaxies, however, are clearly best fit as significantly obscured AGNs. At least half of this sample shows evidence of at least low-luminosity AGN activity, though none could radiatively drive out the remaining molecular gas reservoirs. Therefore, these AGNs are more likely along for the ride, having been fed gas by the same processes driving the transition.

Understanding the halo–galaxy connection is fundamental in order to improve our knowledge on the nature and properties of dark matter. In this work, we build a model that infers the mass of a halo given the positions, velocities, stellar masses, and radii of the galaxies it hosts. In order to capture information from correlations among galaxy properties and their phase space, we use Graph Neural Networks (GNNs), which are designed to work with irregular and sparse data. We train our models on galaxies from more than 2000 state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations project. Our model, which accounts for cosmological and astrophysical uncertainties, is able to constrain the masses of the halos with a ∼0.2 dex accuracy. Furthermore, a GNN trained on a suite of simulations is able to preserve part of its accuracy when tested on simulations run with a different code that utilizes a distinct subgrid physics model, showing the robustness of our method. The PyTorch Geometric implementation of the GNN is publicly available on GitHub (https://github.com/PabloVD/HaloGraphNet).

We present high-cadence optical, ultraviolet (UV), and near-infrared data of the nearby (D ≈ 23 Mpc) Type II supernova (SN) 2021yja. Many Type II SNe show signs of interaction with circumstellar material (CSM) during the first few days after explosion, implying that their red supergiant (RSG) progenitors experience episodic or eruptive mass loss. However, because it is difficult to discover SNe early, the diversity of CSM configurations in RSGs has not been fully mapped. SN 2021yja, first detected within ≈ 5.4 hours of explosion, shows some signatures of CSM interaction (high UV luminosity and radio and x-ray emission) but without the narrow emission lines or early light-curve peak that can accompany CSM. Here we analyze the densely sampled early light curve and spectral series of this nearby SN to infer the properties of its progenitor and CSM. We find that the most likely progenitor was an RSG with an extended envelope, encompassed by low-density CSM. We also present archival Hubble Space Telescope imaging of the host galaxy of SN 2021yja, which allows us to place a stringent upper limit of ≲ 9 M_☉ on the progenitor mass. However, this is in tension with some aspects of the SN evolution, which point to a more massive progenitor. Our analysis highlights the need to consider progenitor structure when making inferences about CSM properties, and that a comprehensive view of CSM tracers should be made to give a fuller view of the last years of RSG evolution.

We examine variations in energetic storm particle (ESP) heavy ion peak intensities and energy spectra at CME-driven interplanetary shocks. We focus on their dependence with heliolongitude relative to the source region of their associated CMEs, and with CME speed, for events observed in Solar Cycle 24 at the STEREO-A, STEREO-B, and/or ACE spacecraft. We find that observations of ESP events at 1 au are organized by longitude relative to their CME solar source. The ESP event longitude distribution also showed organization with CME speed. The near-Sun CME speeds (V_i) for these events ranged from ∼560 to 2650 km s⁻¹ while the average CME transit speeds to 1 au were significantly slower. The angular width of the events had a clear threshold at V_i of ∼1300 km s⁻¹, above which events showed significantly larger angular extension compared to events with speeds below. High-speed events also showed larger heavy ion peak intensities near the nose of the shock compared to the flanks while their spectral index was smaller near the nose and larger near the flanks. This organization for events with V_i < 1300 km s⁻¹ was not as clear. These ESP events were observed over a narrower range of longitudes though the heavy ion peak intensities still appeared largest near the nose of the shock. Their heavy ion spectra showed no clear organization with longitude. These observations highlight the impact of spacecraft position relative to the CME source longitude and V_i on the properties of ESP events at 1 au.

We study the properties of oscillatory double-diffusive convection (ODDC) in the presence of a uniform vertical background magnetic field. ODDC takes place in stellar regions that are unstable according to the Schwarzschild criterion and stable according to the Ledoux criterion (sometimes called semiconvective regions), which are often predicted to reside just outside the core of intermediate-mass main-sequence stars. Previous hydrodynamic studies of ODDC have shown that the basic instability saturates into a state of weak wave-like convection, but that a secondary instability can sometimes transform it into a state of layered convection, where layers then rapidly merge and grow until the entire region is fully convective. We find that magnetized ODDC has very similar properties overall, with some important quantitative differences. A linear stability analysis reveals that the fastest-growing modes are unaffected by the field, but that other modes are. Numerically, the magnetic field is seen to influence the saturation of the basic instability, overall reducing the turbulent fluxes of temperature and composition. This in turn affects layer formation, usually delaying it, and occasionally suppressing it entirely for sufficiently strong fields. Further work will be needed, however, to determine the field strength above which layer formation is actually suppressed in stars. Potential observational implications are briefly discussed.

We determine the distance to the Sculptor Dwarf Spheroidal via three Population II stellar distance indicators: (a) the Tip of the Red Giant Branch (TRGB), (b) RR Lyrae variables (RRLs), and (c) the ridgeline of the blue horizontal branch (HB). High signal-to-noise, wide-field VI imaging that covers an area $48^{\prime} \,\times \,48^{\prime}$ and reaches a photometric depth approximately 2 mag fainter than the HB was acquired with the Magellan-Baade 6.5 m telescope. The true modulus derived from Sculptor's TRGB is found to be ${\mu }_{o}^{\mathrm{TRGB}}=19.59\,\pm \,{0.07}_{\mathrm{stat}}\,\pm \,{0.05}_{\mathrm{sys}}$ mag. Along with periods adopted from the literature, newly acquired RRL phase points are fit with template light curves to determine ${\mu }_{{W}_{I,V-I}}^{\mathrm{RRL}}=19.60\pm {0.01}_{\mathrm{stat}}\pm {0.05}_{\mathrm{sys}}$ mag. Finally, the HB distance is found to be ${\mu }_{o}^{\mathrm{HB}}\,=19.54\pm {0.03}_{\mathrm{stat}}\pm {0.09}_{\mathrm{sys}}$ mag. Absolute calibrations of each method are anchored by independent geometric zero-points, utilize a different class of stars, and are determined from the same photometric calibration.

Many processes during the evolution of protoplanetary disks and during planet formation are highly sensitive to the sizes of dust particles that are present in the disk: the efficiency of dust accretion in the disk and volatile transport on dust particles, gravoturbulent instabilities leading to the formation of planetesimals, or the accretion of pebbles onto large planetary embryos to form giant planets are typical examples of processes that depend on the sizes of the dust particles involved. Furthermore, radiative properties like absorption or scattering opacities depend on the particle sizes. To interpret observations of dust in protoplanetary disks, a proper estimate of the dust particle sizes is needed. We present DustPy: a Python package to simulate dust evolution in protoplanetary disks. DustPy solves gas and dust transport including viscous advection and diffusion as well as collisional growth of dust particles. DustPy is written with a modular concept, such that every aspect of the model can be easily modified or extended to allow for a multitude of research opportunities.

We present the discovery and multiwavelength characterization of SRGA J181414.6-225604, a Galactic hard X-ray transient discovered during the ongoing SRG/ART-XC sky survey. Using data from the Palomar Gattini-IR survey, we identify a spatially and temporally coincident variable infrared (IR) source, IRAS 18111-2257, and classify it as a very-late-type (M7–M8), long-period (1502 ± 24 days), and luminous (M_K ≈ −9.9 ± 0.2) O-rich Mira donor star located at a distance of ≈14.6^+2.9_−2.3 kpc. Combining multicolor photometric data over the last ≈25 yr, we show that the IR counterpart underwent a recent (starting ≈800 days before the X-ray flare) enhanced mass-loss (reaching ≈2.1 × 10⁻⁵M_⊙ yr⁻¹) episode, resulting in an expanding dust shell obscuring the underlying star. Multi-epoch follow-up observations from Swift, NICER, and NuSTAR reveal a ≈200 day long X-ray outburst reaching a peak luminosity of L_X ≈ 2.5 × 10³⁶ erg s⁻¹, characterized by a heavily absorbed (N_H ≈ 6 × 10²² cm⁻²) X-ray spectrum consistent with an optically thick Comptonized plasma. The X-ray spectral and timing behavior suggest the presence of clumpy wind accretion, together with a dense ionized nebula overabundant in silicate material surrounding the compact object. Together, we show that SRGA J181414.6-225604 is a new symbiotic X-ray binary in outburst, triggered by an intense dust-formation episode of a highly evolved donor. Our results offer the first direct confirmation for the speculated connection between enhanced late-stage donor mass loss and the active lifetimes of symbiotic X-ray binaries.

In the two-phase scenario of galaxy formation, a galaxy's stellar mass growth is first dominated by in-situ star formation, and subsequently by accretion. We analyze the radial distribution of the accreted stellar mass in ∼500 galaxies from the (48 Mpc/h)³ box volume of the hydrodynamical cosmological simulation Magneticum, in a stellar-mass range of 10¹⁰ to 10¹²M_⊙. We find that higher-mass galaxies have larger accreted fractions, as found in previous works, but predict generally higher accretion fractions for low-mass galaxies. Based on the 3D radial distribution of the accreted and in-situ components, we define six galaxy classes, from completely accretion to completely in-situ dominated, and measure the transition radii between in-situ and accretion-dominated regions for galaxies that reveal a transition. About 70% of our galaxies have one transition radius. However, about 10% of the galaxies are accretion dominated everywhere, and about 13% have two transition radii, with the center and the outskirts both being accretion dominated. We show that these classes are strongly correlated with the galaxy merger histories, especially with the cold gas fraction at the time of merging. We find high total in-situ (low accretion) fractions to be associated with smaller, lower-mass galaxies, lower central dark-matter fractions, and larger transition radii. Finally, we show that the dips in observed surface brightness profiles seen in many early-type galaxies do not correspond to the transition from in-situ to accretion-dominated regions, and that any inferred mass fractions are not indicative of the true accreted mass but contain information about the galaxies' dry-merger history.

The bubble nebula surrounding NGC 1313 X-2 is believed to be powered by high velocity winds from the central ultraluminous X-ray source (ULX) as a result of supercritical accretion. With the Multi-Unit Spectroscopic Explorer (MUSE) observation of the nebula, we find enhanced [O iii] emission at locations spatially coincident with clusters of stars and the central X-ray source, suggesting that photoionization in addition to shock ionization plays an important role in powering the nebula. The X-ray luminosity of the ULX and the number of massive stars in the nebula region can account for the required ionizing luminosity derived with MAPPINGS V, which also confirms that pure shocks cannot explain the observed emission line ratios.

We present a statistical study of the in situ acceleration of solar wind suprathermal electrons at the Earth's bow shock, utilizing the Wind measurements in ambient solar wind and MMS1 measurements around the bow shock from 2015 September to 2017 December. All the selected 74 cases show significant suprathermal electron acceleration at the bow shock. The observed power-law indexes of accelerated electron energy spectra are significantly larger than the first-order Fermi acceleration prediction and the flux enhancement ratio peaks near a 90° pitch angle, suggesting that the shock drift acceleration (SDA) process plays a crucial role in accelerating suprathermal electrons at the bow shock. According to the observed electron spectral characteristics, the 74 cases can be classified into Types 1, 2, 3 and 4. The electron acceleration efficiency roughly increases from Type 1 to Type 4. For the Type 4 cases with strong electron acceleration, the shocked suprathermal electrons show a double-power-law energy spectrum bending downwards at a break near ∼65 keV with a low-energy spectral index of ∼3.1 and high-energy index of ∼7.6. The observed break energy is comparable to a critical electron energy ${\varepsilon }_{\mathrm{ramp}}^{{dn}}$ with its cross-shock gyrodiameter equal to the shock's ramp thickness D_ramp. At energies below (above) ${\varepsilon }_{\mathrm{ramp}}^{{dn}}$, the accelerated electrons can be effectively confined into (can easily escape from) the SDA process with a roughly energy-independent (energy-decreasing) drift time and a probably conserved (nonadiabatically reduced) magnetic moment, since their gyrodiameter is less (greater) than D_ramp. Therefore, such an SDA process could produce the observed double-power-law spectrum bending downwards at a break energy that is associated with D_ramp.

In 2020, the Super-Kamiokande (SK) experiment moved to a new stage (SK-Gd) in which gadolinium (Gd) sulfate octahydrate was added to the water in the detector, enhancing the efficiency to detect thermal neutrons and consequently improving the sensitivity to low energy electron anti-neutrinos from inverse beta decay (IBD) interactions. SK-Gd has the potential to provide early alerts of incipient core-collapse supernovae through detection of electron anti-neutrinos from thermal and nuclear processes responsible for the cooling of massive stars before the gravitational collapse of their cores. These pre-supernova neutrinos emitted during the silicon burning phase can exceed the energy threshold for IBD reactions. We present the sensitivity of SK-Gd to pre-supernova stars and the techniques used for the development of a pre-supernova alarm based on the detection of these neutrinos in SK, as well as prospects for future SK-Gd phases with higher concentrations of Gd. For the current SK-Gd phase, high-confidence alerts for Betelgeuse could be issued up to 9 hr in advance of the core collapse itself.

Magnetohydrodynamic waves are ubiquitously detected in the finely structured solar atmosphere. At the same time, our Sun is a highly dynamic plasma environment, giving rise to flows of various magnitudes, which can lead to the instability of waveguides. Recent studies have employed the method of introducing waveguide asymmetry to generalize "classical" symmetric descriptions of the fine structuring within the solar atmosphere, with some of them introducing steady flows as well. Building on these recent studies, here we investigate the magnetoacoustic waves guided by a magnetic slab within an asymmetric magnetic environment, in which the slab is under the effect of a steady flow. We provide an analytical investigation of how the phase speeds of the guided waves are changed, and where possible, determine the limiting flow speeds required for the onset of the Kelvin–Helmholtz instability. Furthermore, we complement the study with initial numerical results, which allows us to demonstrate the validity of our approximations and extend the investigation to a wider parameter regime. This configuration is part of a series of studies aimed to generalize, step-by-step, well-known symmetric waveguide models and understand the additional physics stemming from introducing further sources of asymmetry.

Active galactic nuclei (AGN) show a range of morphologies and dynamical properties, which are determined not only by parameters intrinsic to the central engine but also their interaction with the surrounding environment. We investigate the connection of kiloparsec scale AGN jet properties to their intrinsic parameters and surroundings. This is done using a suite of 40 relativistic hydrodynamic simulations spanning a wide range of engine luminosities and opening angles. We explore AGN jet propagation with different ambient density profiles, including r⁻² (self-similar solution) and r⁻¹, which is more relevant for AGN host environments. While confirmation awaits future 3D studies, the Fanaroff–Riley (FR) morphological dichotomy arises naturally in our 2D models. Jets with low energy density compared to the ambient medium produce a center-brightened emissivity distribution, while emissivity from relatively higher energy density jets is dominated by the jet head. We observe recollimation shocks in our simulations that can generate bright spots along the spine of the jet, providing a possible explanation for "knots" observed in AGN jets. We additionally find a scaling relation between the number of knots and the jet-head-to-surroundings energy density ratio. This scaling relation is generally consistent with the observations of the jets in M87 and Cygnus A. Our model also correctly predicts M87 as FRI and Cygnus A as FRII. Our model can be used to relate jet dynamical parameters such as jet head velocity, jet opening angle, and external pressure to jet power, and ambient density estimates.

Measuring the evolution of X-ray emission from pre-main-sequence (PMS) stars gives insight into two issues: the response of magnetic dynamo processes to changes in the interior structure, and the effects of high-energy radiation on protoplanetary disks and primordial planetary atmospheres. We present a sample of 6003 stars with ages 7–25 Myr in 10 nearby open clusters from Chandra X-ray and Gaia-EDR3 surveys. Combined with previous results in large samples of younger (≲5 Myr) stars in MYStIX and SFiNCs star-forming regions, mass-stratified activity-age relations are derived for the early phases of stellar evolution. X-ray luminosity (L_X) is constant during the first few Myr, possibly due to the presence of extended X-ray coronas insensitive to temporal changes in stellar size. L_X then decays during the 7–25 Myr period, more rapidly as stellar mass increases. This decay is interpreted as decreasing efficiency of the α² dynamo as radiative cores grow and a solar-type αΩ dynamo emerges. For more massive 3.5–7 M_⊙ fully radiative stars, the X-ray emission plummets—indicating the lack of an effective magnetic dynamo. The findings provide improved measurements of high-energy radiation effects on circumstellar material, first for the protoplanetary disk and then for the atmospheres of young planets. The observed X-ray luminosities can be so high that an inner Earth-mass rocky, unmagnetized planet around a solar-mass PMS star might lose its primary and secondary atmospheres within a few (several) million years. PMS X-ray emission may thus have a significant impact on the evolution of early-planetary atmospheres and the conditions promoting the rise of habitability.

We report the Fermi LAT γ-ray detection of the 2021 outburst of the symbiotic recurrent nova RS Ophiuchi. In this system, unlike classical novae from cataclysmic binaries, the ejecta from the white dwarf form shocks when interacting with the dense circumstellar wind environment of the red giant companion. We find the LAT spectra from 50 MeV to ∼20–23 GeV, the highest-energy photons detected in some subintervals, are consistent with π⁰-decay emission from shocks in the ejecta as proposed by Tatischeff & Hernanz for its previous 2006 outburst. The LAT light curve displayed a fast rise to its peak >0.1 GeV flux of ≃6 × 10⁻⁶ ph cm⁻² s⁻¹ beginning on day 0.745 after its optically constrained eruption epoch of 2021 August 8.50. The peak lasted for ∼1 day and exhibited a power-law decline up to the final LAT detection on day 45. We analyze the data on shorter timescales at early times and found evidence of an approximate doubling of emission over ∼200 minutes at day 2.2, possibly indicating a localized shock-acceleration event. Comparing the data collected by the American Association of Variable Star Observers, we measured a constant ratio of ∼ 2.8 × 10⁻³ between the γ-ray and optical luminosities except for a ∼5×smaller ratio within the first day of the eruption likely indicating attenuation of γ rays by ejecta material and lower high-energy proton fluxes at the earliest stages of the shock development. The hard X-ray emission due to bremsstrahlung from shock-heated gas traced by the Swift-XRT 2–10 keV light curve peaked at day ∼6, later than at GeV and optical energies. Using X-ray derived temperatures to constrain the velocity profile, we find the hadronic model reproduces the observed >0.1 GeV light curve.

Solar flares create adverse space weather impacting space- and Earth-based technologies. However, the difficulty of forecasting flares, and by extension severe space weather, is accentuated by the lack of any unique flare trigger or a single physical pathway. Studies indicate that multiple physical properties contribute to active region flare potential, compounding the challenge. Recent developments in machine learning (ML) have enabled analysis of higher-dimensional data leading to increasingly better flare forecasting techniques. However, consensus on high-performing flare predictors remains elusive. In the most comprehensive study to date, we conduct a comparative analysis of four popular ML techniques (k nearest neighbors, logistic regression, random forest classifier, and support vector machine) by training these on magnetic parameters obtained from the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory for the entirety of solar cycle 24. We demonstrate that the logistic regression and support vector machine algorithms perform extremely well in forecasting active region flaring potential. The logistic regression algorithm returns the highest true skill score of 0.967 ± 0.018, possibly the highest classification performance achieved with any strictly parametric study. From a comparative assessment, we establish that magnetic properties like total current helicity, total vertical current density, total unsigned flux, R_VALUE, and total absolute twist are the top-performing flare indicators. We also introduce and analyze two new performance metrics, namely, severe and clear space weather indicators. Our analysis constrains the most successful ML algorithms and identifies physical parameters that contribute most to active region flare productivity.

In this paper, we develop a 3D implicit single-fluid magnetohydrodynamic (MHD) model to simulate the steady-state solar corona with a wide range of Mach numbers and low plasma β. We employ a low-dissipation advection upstream splitting method (AUSM) to calculate the convective flux in the regions of low Mach numbers for a high resolution, and hybridize the AUSM with Harten-Lax-van Leer Riemann solver in the regions of high Mach numbers to improve the solver's robustness. The inner boundary condition of no backflow is implemented by numerical flux. A reconstruction method based on the divergence-free radial basis function is adopted to enhance the divergence-free constraint of magnetic field. Also, an anisotropic thermal conduction term is considered; the positivity-preserving reconstruction method is used to prevent the presence of negative thermal pressure and plasma density, and the implicit lower-upper symmetric Gauss Seidel method is implemented for a better convergence rate. After establishing the implicit solar wind MHD model, we employ it to simulate steady-state solar coronal structures in Carrington rotations 2177 and 2212. The simulations demonstrate that the MHD model's computational efficiency is desirable, and the modeled results are basically in agreement with the solar coronal observations and the mapped in situ measurements from the OMNI archive. Consequently, this implicit MHD model is promising to simulate a complex plasma environment with high-intensity magnetic field and wide-ranging Mach numbers.

We present a new methodology—the Keplerian Optical Dynamics Analysis (KODA)—for quantifying the dynamics of erupting magnetic structures in the solar corona. The technique involves adaptive spatiotemporal tracking of propagating intensity gradients and their characterization in terms of time-evolving Keplerian areas swept out by the position vectors of moving plasma blobs. Whereas gravity induces purely ballistic motions consistent with Kepler's second law, noncentral forces such as the Lorentz force introduce nonzero torques resulting in more complex motions. KODA algorithms enable direct evaluation of the line-of-sight component of the net torque density from the image-plane projection of the areal acceleration. The method is applied to the prominence eruption of 2011 June 7, observed by the Solar Dynamics Observatory's Atmospheric Imaging Assembly. Results obtained include quantitative estimates of the magnetic forces, field intensities, and blob masses and energies across a vast region impacted by the postreconnection redistribution of the prominence material. The magnetic pressure and energy are strongly dominant during the early, rising phase of the eruption, while the dynamic pressure and kinetic energy become significant contributors during the subsequent falling phases. Measured intensive properties of the prominence blobs are consistent with those of typical active-region prominences; measured extensive properties are compared with those of the whole pre-eruption prominence and the post-eruption coronal mass ejection of 2011 June 7, all derived by other investigators and techniques. We show that KODA provides valuable information on spatially and temporally dependent characteristics of coronal eruptions that is not readily available via alternative means, thereby shedding new light on the environment and evolution of these solar events.

Warped disk galaxies are classified into two morphologies: S and U types. Conventional theories routinely attribute both types to galactic tidal interaction and/or gas accretion, but reproducing U types in simulations is extremely challenging. Here we investigate whether both types are governed by the same mechanisms using the most extensive sample of ∼8000 nearby (0.02 < z < 0.06) massive (M_*/M_⊙ > 10⁹) edge-on disks from the Sloan Digital Sky Survey. We find that U types show on average bluer optical colors and a higher specific star formation rate (sSFR) than S types, with more strongly warped U types having a higher sSFR. We also find that while the S-type warp properties correlate with the tidal force by the nearest neighbor regardless of the environment, there is no such correlation for U types in groups/clusters, suggesting a nontidal mechanism could be at play for U types, such as ram pressure stripping (RPS). Indeed, U types are more common in groups/clusters than in fields and they have stellar mass, gas fraction, sSFR enhancement, and phase-space distribution closely analogous to RPS-induced jellyfish galaxies in clusters. We furthermore show that the stellar disks of most RPS galaxies in the IllustrisTNG simulation are warped in a U shape and bent in the opposite direction of stripped gas tails, satisfying theoretical expectations for stellar warps embedded in jellyfishes. We therefore suggest that despite the majority of U types that live in fields being still less explained, RPS can be an alternative origin for those in groups/clusters.

We present GIGA-Lens: a gradient-informed, GPU-accelerated Bayesian framework for modeling strong gravitational lensing systems, implemented in TensorFlow and JAX. The three components, optimization using multistart gradient descent, posterior covariance estimation with variational inference, and sampling via Hamiltonian Monte Carlo, all take advantage of gradient information through automatic differentiation and massive parallelization on graphics processing units (GPUs). We test our pipeline on a large set of simulated systems and demonstrate in detail its high level of performance. The average time to model a single system on four Nvidia A100 GPUs is 105 s. The robustness, speed, and scalability offered by this framework make it possible to model the large number of strong lenses found in current surveys and present a very promising prospect for the modeling of ${ \mathcal O }({10}^{5})$ lensing systems expected to be discovered in the era of the Vera C. Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope.

We use panoramic optical spectroscopy obtained with the Very Large Telescope/MUSE to investigate the nature of five candidate extremely isolated low-mass star-forming regions (Blue Candidates; hereafter, BCs) toward the Virgo cluster of galaxies. Four of the five (BC1, BC3, BC4, and BC5) are found to host several H ii regions and to have radial velocities fully compatible with being part of the Virgo cluster. All the confirmed candidates have mean metallicity significantly in excess of that expected from their stellar mass, indicating that they originated from gas stripped from larger galaxies. In summary, these four candidates share the properties of the prototype system SECCO 1, suggesting the possible emergence of a new class of stellar systems, intimately linked to the complex duty cycle of gas within clusters of galaxies. A thorough discussion of the nature and evolution of these objects is presented in a companion paper, where the results obtained here from the MUSE data are complemented with Hubble Space Telescope (optical) and Very Large Array (Hi) observations.

We discuss five blue stellar systems in the direction of the Virgo cluster, analogous to the enigmatic object SECCO 1 (AGC 226067). These objects were identified based on their optical and UV morphology and followed up with H i observations with the Very Large Array (and Green Bank Telescope), Multi Unit Spectroscopic Explorer (on the Very Large Telescope) optical spectroscopy, and Hubble Space Telescope imaging. These new data indicate that one system is a distant group of galaxies. The remaining four are extremely low mass (M_* ∼ 10⁵M_⊙), are dominated by young blue stars, have highly irregular and clumpy morphologies, are only a few kiloparsecs across, yet host an abundance of metal-rich, $12+\mathrm{log}({\rm{O}}/{\rm{H}})\gt 8.2$, H ii regions. These high metallicities indicate that these stellar systems formed from gas stripped from much more massive galaxies. Despite the young age of their stellar populations, only one system is detected in H i, while the remaining three have minimal (if any) gas reservoirs. Furthermore, two systems are surprisingly isolated and have no plausible parent galaxy within ∼30' (∼140 kpc). Although tidal stripping cannot be conclusively excluded as the formation mechanism of these objects, ram pressure stripping more naturally explains their properties, in particular their isolation, owing to the higher velocities, relative to the parent system, that can be achieved. Therefore, we posit that most of these systems formed from ram-pressure-stripped gas removed from new infalling cluster members and survived in the intracluster medium long enough to become separated from their parent galaxies by hundreds of kiloparsecs and that they thus represent a new type of stellar system.

Recent evidence suggests that high-redshift Lyα emitting galaxies (LAEs) with $\mathrm{log}L(\mathrm{Ly}\alpha )\gt 43.5\,\mathrm{erg}\ {{\rm{s}}}^{-1}$, referred to as ultraluminous LAEs (ULLAEs), may show less evolution than lower-luminosity LAEs in the redshift range z = 5.7–6.6. Here we explore the redshift evolution of the velocity widths of the Lyα emission lines in LAEs over this redshift interval. We use new wide-field, narrowband observations from Subaru/Hyper Suprime-Cam to provide a sample of 24 z = 6.6 and 12 z = 5.7 LAEs with $\mathrm{log}L(\mathrm{Ly}\alpha )\gt 43\,\mathrm{erg}\ {{\rm{s}}}^{-1}$, all of which have follow-up spectroscopy from Keck/DEIMOS. Combining with archival lower-luminosity data, we find a significant narrowing of the Lyα lines in LAEs at $\mathrm{log}L(\mathrm{Ly}\alpha )\lt 43.25\,\mathrm{erg}\ {{\rm{s}}}^{-1}$—somewhat lower than the usual ULLAE definition—at z = 6.6 relative to those at z = 5.7, but we do not see this in higher-luminosity LAEs. As we move to higher redshifts, the increasing neutrality of the intergalactic medium should increase the scattering of the Lyα lines, making them narrower. The absence of this effect in the higher-luminosity LAEs suggests they may lie in more highly ionized regions, self-shielding from the scattering effects of the intergalactic medium.

To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass M_GC ∼ 10⁵M_⊙ forms within ≈3 Myr when two GMCs with mass M_GMC = 3.6 × 10⁵M_⊙ collide. The GC candidate does not form during a single collapsing event but emerges due to the mergers of local dense gas clumps and gas accretion. The momentum transfer due to the absorption of the ionizing radiation is the dominant feedback process that suppresses the gas collapse, and photoionization becomes efficient once a sufficient number of stars form. The cluster mass is larger by a factor of ∼2 when the radiation feedback is neglected, and the difference is slightly more pronounced (16%) when extreme Lyα feedback is considered in the fiducial run. In the simulations with radiation feedback, supernovae explode after the star-forming clouds are dispersed, and their metal ejecta are not instantaneously recycled to form stars.

Interpreting the short-timescale variability of the accreting, young, low-mass stars known as Classical T Tauri stars remains an open task. Month-long, continuous light curves from the Transiting Exoplanet Survey Satellite (TESS) have become available for hundreds of T Tauri stars. With this vast data set, identifying connections between the variability observed by TESS and short-timescale accretion variability is valuable for characterizing the accretion process. To this end, we obtained short-cadence TESS observations of 14 T Tauri stars in the Taurus star formation region along with simultaneous ground-based, UBVRI-band photometry to be used as accretion diagnostics. In addition, we combine our data set with previously published simultaneous near-UV–near-IR Hubble Space Telescope spectra for one member of the sample. We find evidence that much of the short-timescale variability observed in the TESS light curves can be attributed to changes in the accretion rate, but note significant scatter between separate nights and objects. We identify hints of time lags within our data set that increase at shorter wavelengths, which we suggest may be evidence of longitudinal density stratification of the accretion column. Our results highlight that contemporaneous, multiwavelength observations remain critical for providing context for the observed variability of these stars.

Solar magnetic fields comprise an 11 yr activity cycle, represented by the number of sunspots. The maintenance of such a solar magnetic field can be attributed to fluid motion in the convection zone, i.e., a dynamo. This study conducts the mean-field analyses of the global solar dynamo simulation presented by Hotta et al. (2016). Although the study succeeds in producing coherent large-scale magnetic fields at high Reynolds numbers, the detailed physics of the maintenance of these fields have not been fully understood. This study extracts the α tensor and the turbulent magnetic diffusivity tensor β through mean-field analyses. The turbulent magnetic diffusivity exhibits a significant decrease toward high Reynolds numbers. The decrease in the turbulent magnetic diffusivity suppresses the energy conversion of large-scale field to small-scale field. This implies that the decrease in the turbulent magnetic diffusivity contributes to the maintenance of a large-scale magnetic field at high Reynolds numbers. A significant downward turbulent pumping is observed; it is enhanced in the weak phase of the large-scale field. This study proposes a cyclic reversal process of a large-scale field, which is dominantly driven by the α effect and is possibly triggered by downward pumping.

Gas accretion onto the circumplanetary disks and the source region of accreting gas are important to reveal dust accretion that leads to satellite formation around giant planets. We performed local three-dimensional high-resolution hydrodynamic simulations of an isothermal and inviscid gas flow around a planet to investigate the planetary-mass dependence of the gas accretion bandwidth and gas accretion rate onto circumplanetary disks. We examined cases with various planetary masses corresponding to M_p = 0.05–1M_Jup at 5.2 au, where M_Jup is the current Jovian mass. We found that the radial width of the gas accretion band is proportional to ${M}_{{\rm{p}}}^{1/6}$ for the low-mass regime with M_p ≲ 0.2M_Jup while it is proportional to M_p for the high-mass regime with M_p ≳ 0.2M_Jup. We found that the ratio of the mass accretion rate onto the circumplanetary disk to that into the Hill sphere is about 0.4 regardless of the planetary mass for the cases we examined. Combining our results with the gap model obtained from global hydrodynamic simulations, we derive a semi-analytical formula of mass accretion rate onto circumplanetary disks. We found that the mass dependence of our three-dimensional accretion rates is the same as the previously obtained two-dimensional case, although the qualitative behavior of accretion flow onto the circumplanetary disk is quite different between the two cases.

High-resolution observations of the Perseus B5 "core" have previously revealed that this subsonic region actually consists of several filaments that are likely in the process of forming a quadruple stellar system. Since subsonic filaments are thought to be produced at the ∼0.1 pc sonic scale by turbulent compression, a detailed kinematic study is crucial to test such a scenario in the context of core and star formation. Here we present a detailed kinematic follow-up study of the B5 filaments at a 0.009 pc resolution using the VLA and GBT combined observations fitted with multicomponent spectral models. Using precisely identified filament spines, we find a remarkable resemblance between the averaged width profiles of each filament and Plummer-like functions, with filaments possessing FWHM widths of ∼0.03 pc. The velocity dispersion profiles of the filaments also show decreasing trends toward the filament spines. Moreover, the velocity gradient field in B5 appears to be locally well ordered (∼0.04 pc) but globally complex, with kinematic behaviors suggestive of inhomogeneous turbulent accretion onto filaments and longitudinal flows toward a local overdensity along one of the filaments.

In a recent paper, we investigated possible systematic uncertainties related to the Cepheid color–luminosity calibration method and their influence on the tension between the Hubble constant as inferred from distances to Type Ia supernovae and the cosmic microwave background as measured with the Planck satellite. Here, we study the impact of other sources of uncertainty in the supernova distance ladder, including Cepheid temperature and metallicity variations, supernova magnitudes, and GAIA parallax distances. Using Cepheid data in 19 Type Ia supernova host galaxies from Riess et al., anchor data from Riess et al., and a set of recalibrated Milky Way Cepheid distances, we obtain H₀ = 71.9 ± 2.2 km s⁻¹ Mpc⁻¹, 2.0σ from the Planck value. Excluding Cepheids with estimated color excesses $\hat{E}(V-I)=0.15$ mag to mitigate the impact of the Cepheid color–luminosity calibration, the inferred Hubble constant is H₀ = 68.1 ± 2.6 km s⁻¹ Mpc⁻¹, removing the tension with the Planck value.

We carried out Chandra, Hubble Space Telescope (HST), and Very Large Array observations of four MOJAVE blazars that have previously been classified as hybrid (FR I/II) blazars in terms of radio morphology but not total radio power. The motivation of this study is to determine the X-ray emission mechanism in jets, these being different in FR I and FR II jets. We detected X-ray jet emission with sufficient signal to noise in two blazars, viz. PKS 0215+015 and TXS 0730+504. We carried out spectral energy distribution modeling of the broadband emission from the jet regions in these sources and found that a single synchrotron emission model is ruled out due to the deep upper limits obtained from HST optical and IR data. The inverse Compton-cosmic microwave background model can reproduce the X-ray jet emission in both sources although the model requires extreme jet parameters. Both our sources possess FR II-like radio powers and our results are consistent with previous studies suggesting that radio power is more important than FR morphology in determining the emission mechanism of X-ray jets.

The conservation of wave action in moving plasmas has been well known for over half a century. However, wave action is not conserved when multiple wave modes propagate and coexist close to the degeneration condition (where the sound speed equals the Alfvén speed, i.e., plasma β ∼ 1). Here, we show that the violation of conservation is due to wave mode conversion, and that the total wave action summed over the interacting modes is still conserved. Though the result is general, we focus on MHD waves and identify three distinctive mode-conversion mechanisms, i.e., degeneracy, linear mode conversion, and resonance, and provide an intuitive physical picture for the mode-conversion processes. We use one-dimensional MHD simulations with the Expanding Box Model to simulate the nonlinear evolution of monochromatic MHD waves in the expanding solar wind. The simulation results validate the theory; total wave action therefore remains an interesting diagnostic for studies of waves and turbulence in the solar wind.

We report measurements of the gravitationally lensed secondary image—the first in an infinite series of so-called "photon rings"—around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford–Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of ${4.20}_{-0.06}^{+0.12}\,\mu \mathrm{as}$ and a corresponding black hole mass of (7.13 ± 0.39) × 10⁹M_⊙, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.