Abstract: Abstract We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems colocated with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to ≤11 known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with six undetected in radio. Overall, ≥236 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power E ̇ decreases to its observed minimum near 10 33 erg s −1 , approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties, and fit results to guide and be compared with modeling results.

Abstract: We present results on the morphological and structural evolution of a total of 3956 galaxies observed with JWST at 1.5 < z < 6.5 in the JWST CEERS observations that overlap with the CANDELS EGS field. This is the biggest visually classified sample observed with JWST yet, ∼20 times larger than previous studies, and allows us to examine in detail how galaxy structure has changed over this critical epoch. All sources were classified by six individual classifiers using a simple classification scheme aimed at producing disk/spheroid/peculiar classifications, whereby we determine how the relative number of these morphologies has evolved since the Universe’s first billion years. Additionally, we explore structural and quantitative morphology measurements using Morfometryka, and show that galaxies with M * > 109 M ⊙ at z > 3 are not dominated by irregular and peculiar structures, either visually or quantitatively, as previously thought. We find a strong dominance of morphologically selected disk galaxies up to z = 6 in this mass range. We also find that the stellar mass and star formation rate densities are dominated by disk galaxies up to z ∼ 6, demonstrating that most stars in the Universe were likely formed in a disk galaxy. We compare our results to theory to show that the fraction of types we find is predicted by cosmological simulations, and that the Hubble Sequence was already in place as early as one billion years after the Big Bang. Additionally, we make our visual classifications public for the community.

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

Abstract: We report the discovery of four galaxy candidates observed 450–600 Myr after the Big Bang with photometric redshifts between z ∼ 8.3 and 10.2 measured using James Webb Space Telescope (JWST) NIRCam imaging of the galaxy cluster WHL0137−08 observed in eight filters spanning 0.8–5.0 μm, plus nine Hubble Space Telescope filters spanning 0.4–1.7 μm. One candidate is gravitationally lensed with a magnification of μ ∼ 8, while the other three are located in a nearby NIRCam module with expected magnifications of μ ≲ 1.1. Using SED fitting, we estimate the stellar masses of these galaxies are typically in the range logM⋆/M⊙ = 8.3–8.7. All appear young, with mass-weighted ages <240 Myr, low dust content A V < 0.15 mag, and specific star formation rates sSFR ∼0.25–10 Gyr−1 for most. One z ∼ 9 candidate is consistent with an age <5 Myr and an sSFR ∼10 Gyr−1, as inferred from a strong F444W excess, implying [O iii ]+H β rest-frame equivalent width ∼2000 Å, although an older z ∼ 10 object is also allowed. Another z ∼ 9 candidate is lensed into an arc 2.″4 long with a magnification of μ ∼ 8. This arc is the most spatially resolved galaxy at z ∼ 9 known to date, revealing structures ∼30 pc across. Follow-up spectroscopy of WHL0137−08 with JWST/NIRSpec will be useful to spectroscopically confirm these high-redshift galaxy candidates and to study their physical properties in more detail.

Abstract: Abstract A gravitational anomaly is found at weak gravitational acceleration g N ≲ 10 −9 m s −2 from analyses of the dynamics of wide binary stars selected from the Gaia DR3 database that have accurate distances, proper motions, and reliably inferred stellar masses. Implicit high-order multiplicities are required and the multiplicity fraction is calibrated so that binary internal motions agree statistically with Newtonian dynamics at a high enough acceleration of ≈10 −8 m s −2 . The observed sky-projected motions and separation are deprojected to the 3D relative velocity v and separation r through a Monte Carlo method, and a statistical relation between the Newtonian acceleration g N ≡ GM / r 2 (where M is the total mass of the binary system) and a kinematic acceleration g ≡ v 2 / r is compared with the corresponding relation predicted by Newtonian dynamics. The empirical acceleration relation at ≲10 −9 m s −2 systematically deviates from the Newtonian expectation. A gravitational anomaly parameter δ obs−newt between the observed acceleration at g N and the Newtonian prediction is measured to be: δ obs−newt = 0.034 ± 0.007 and 0.109 ± 0.013 at g N ≈ 10 −8.91 and 10 −10.15 m s −2 , from the main sample of 26,615 wide binaries within 200 pc. These two deviations in the same direction represent a 10 σ significance. The deviation represents a direct evidence for the breakdown of standard gravity at weak acceleration. At g N = 10 −10.15 m s −2 , the observed to Newton-predicted acceleration ratio is g obs / g pred = 10 2 δ obs − newt = 1.43 ± 0.06 . This systematic deviation agrees with the boost factor that the AQUAL theory predicts for kinematic accelerations in circular orbits under the Galactic external field.

Abstract: Abstract Roughly half of the quasiperiodic eruption (QPE) sources in galactic nuclei exhibit a remarkably regular alternating “long-short” pattern of recurrence times between consecutive flares. We show that a main-sequence star (brought into the nucleus as an extreme mass-ratio inspiral; EMRI) that passes twice per orbit through the accretion disk of the supermassive black hole (SMBH) on a mildly eccentric inclined orbit, each time shocking and ejecting optically thick gas clouds above and below the midplane, naturally reproduces observed properties of QPE flares. Inefficient photon production in the ejecta renders the QPE emission much harder than the blackbody temperature, enabling the flares to stick out from the softer quiescent disk spectrum. Destruction of the star via mass ablation limits the QPE lifetime to decades, precluding a long-lived AGN as the gaseous disk. By contrast, a tidal disruption event (TDE) naturally provides a transient gaseous disk on the requisite radial scale, with a rate exceeding the EMRI inward migration rate, suggesting that many TDEs should host a QPE. This picture is consistent with the X-ray TDE observed several years prior to the QPE appearance from GSN 069. Remarkably, a second TDE-like flare was observed from this event, starting immediately after detectable QPE activity ceased; this event could plausibly result from the (partial or complete) destruction of the QPE-generating star triggered by runaway mass loss, though other explanations cannot be excluded. Our model can also be applied to black hole–disk collisions, such as those invoked in the context of the candidate SMBH binary OJ 287.

Abstract: We report the discovery of 15 exceptionally luminous 10 ≲ z ≲ 14 candidate galaxies discovered in the first 0.28 deg2 of JWST/NIRCam imaging from the COSMOS-Web survey. These sources span rest-frame UV magnitudes of −20.5 > M UV > −22, and thus constitute the most intrinsically luminous z ≳ 10 candidates identified by JWST to date. Selected via NIRCam imaging, deep ground-based observations corroborate their detection and help significantly constrain their photometric redshifts. We analyze their spectral energy distributions using multiple open-source codes and evaluate the probability of low-redshift solutions; we conclude that 12/15 (80%) are likely genuine z ≳ 10 sources and 3/15 (20%) likely low-redshift contaminants. Three of our z ∼ 12 candidates push the limits of early stellar mass assembly: they have estimated stellar masses ∼ 5 × 109 M ⊙, implying an effective stellar baryon fraction of ϵ ⋆ ∼ 0.2−0.5, where ϵ ⋆ ≡ M ⋆/(f b M halo). The assembly of such stellar reservoirs is made possible due to rapid, burst-driven star formation on timescales < 100 Myr where the star formation rate may far outpace the growth of the underlying dark matter halos. This is supported by the similar volume densities inferred for M ⋆ ∼ 1010 M ⊙ galaxies relative to M ⋆ ∼ 109 M ⊙—both about 10−6 Mpc−3—implying they live in halos of comparable mass. At such high redshifts, the duty cycle for starbursts would be of order unity, which could cause the observed change in the shape of the UV luminosity function from a double power law to a Schechter function at z ≈ 8. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understand how, and if, such early massive galaxies push the limits of galaxy formation in the Lambda cold dark matter paradigm.

Abstract: We present the results of an all-sky search for continuous gravitational waves in the public LIGO O3 data. The search covers signal frequencies 20.0 Hz ≤ f ≤ 800.0 Hz and a spin-down range down to −2.6 × 10−9 Hz s−1, motivated by detectability studies on synthetic populations of Galactic neutron stars. This search is the most sensitive all-sky search to date in this frequency/spin-down region. The initial search was performed using the first half of the public LIGO O3 data (O3a), utilizing graphical processing units provided in equal parts by the volunteers of the Einstein@Home computing project and by the ATLAS cluster. After a hierarchical follow-up in seven stages, 12 candidates remain. Six are discarded at the eighth stage, by using the remaining O3 LIGO data (O3b). The surviving six can be ascribed to continuous-wave fake signals present in the LIGO data for validation purposes. We recover these fake signals with very high accuracy with our last stage search, which coherently combines all O3 data. Based on our results, we set upper limits on the gravitational-wave amplitude h 0 and translate these into upper limits on the neutron star ellipticity and on the r-mode amplitude. The most stringent upper limits are at 203 Hz, with h 0 = 8.1 × 10−26 at the 90% confidence level. Our results exclude isolated neutron stars rotating faster than 5 ms with ellipticities greater than 5×10−8d100pc within a distance d from Earth and r-mode amplitudes α≥10−5d100pc for neutron stars spinning faster than 150 Hz.

Abstract: A Type Ia supernova (SN) at z = 1.78 was discovered in James Webb Space Telescope Near Infrared Camera imaging of the galaxy cluster PLCK G165.7+67.0 (G165; z = 0.35). The SN is situated 1.5–2 kpc from the host-galaxy nucleus and appears in three different locations as a result of gravitational lensing by G165. These data can yield a value for Hubble’s constant using time delays from this multiply imaged SN Ia that we call “SN H0pe.” Over the cluster, we identified 21 image multiplicities, confirmed five of them using the Near-Infrared Spectrograph, and constructed a new lens model that gives a total mass within 600 kpc of (2.6 ± 0.3) × 1014 M ⊙. The photometry uncovered a galaxy overdensity coincident with the SN host galaxy. NIRSpec confirmed six member galaxies, four of which surround the SN host galaxy with relative velocity ≲900 km s−1 and projected physical extent ≲33 kpc. This compact galaxy group is dominated by the SN host galaxy, which has a stellar mass of (5.0 ± 0.1) × 1011 M ⊙. The group members have specific star formation rates of 2–260 Gyr−1 derived from the Hα-line fluxes corrected for stellar absorption, dust extinction, and slit losses. Another group centered on a strongly lensed dusty star-forming galaxy is at z = 2.24. The total (unobscured and obscured) SFR of this second galaxy group is estimated to be (≳ 100 M ⊙ yr−1), which translates to a supernova rate of ∼1 SNe yr−1, suggesting that regular monitoring of this cluster may yield additional SNe.

Abstract: It is found that Gaia Data Release 3 binary stars selected with stringent requirements on astrometric measurements and radial velocities naturally satisfy Newtonian dynamics without hidden close companions when the projected separation s ≲ 2 kau, showing that pure binaries can be selected. It is then found that pure binaries selected with the same criteria show a systematic deviation from the Newtonian expectation when s ≳ 2 kau. When both proper motions and parallaxes are required to have precision better than 0.005 and radial velocities better than 0.2, I obtain 2463 statistically pure binaries within a “clean” G-band absolute magnitude range. From this sample, I obtain an observed-to-Newtonian predicted kinematic acceleration ratio of γg=gobs/gpred=1.49−0.19+0.21 for acceleration ≲10−10 m s−2, in excellent agreement with 1.49 ± 0.07 for a much larger general sample with the amount of hidden close companions self-calibrated. I also investigate the radial profile of stacked sky-projected relative velocities without a deprojection to the 3D space. The observed profile matches the Newtonian predicted profile for s ≲ 2 kau without any free parameters, but shows a clear deviation at a larger separation with a significance of ≈5.0σ. The projected velocity boost factor for s ≳ 5 kau is measured to be γvp=1.20±0.06 (stat) ±0.05 (sys) matching γg . Finally, for a small sample of 40 binaries with exceptionally precise radial velocities (fractional errors <0.005), the directly measured relative velocities in the 3D space also show a boost at larger separations. These results robustly confirm the recently reported gravitational anomaly at low acceleration for a general sample.

Abstract: Abstract The IceCube Neutrino Observatory sends realtime neutrino alerts with a high probability of being astrophysical in origin. We present a new method to correlate these events and possible candidate sources using 2089 blazars from the Fermi-LAT 4LAC-DR2 catalog and with 3413 active galactic nuclei (AGNs) from the Radio Fundamental Catalog. No statistically significant neutrino emission was found in any of the catalog searches. The result suggests that a small fraction, <1%, of the studied AGNs emit neutrinos that pass the alert criteria, and is compatible with prior evidence for neutrino emission presented by IceCube and other authors from sources such as TXS 0506 + 056 and PKS 1502 + 106. We also present cross-checks to other analyses that claim a significant correlation using similar data samples.

Abstract: Abstract We examine the factors responsible for the variation in the ionization parameter ( U ) of high-redshift star-forming galaxies based on medium-resolution JWST/NIRSpec observations obtained by the Cosmic Evolution Early Release Science survey. The sample consists of 48 galaxies with redshifts z spec = 2.7−6.3, which are largely representative of typical galaxies at these redshifts. The [S ii ] λ λ 6718, 6733 doublet is used to estimate electron densities ( n e ), and dust-corrected H α luminosities are used to compute ionizing photon rates ( Q ). Using composite spectra of galaxies in bins of [O iii ] λ λ 4960, 5008/[O ii ] λ λ 3727, 3730 (O32) as a proxy for U , we determine that galaxies with higher O32 have 〈 n e 〉 ≃ 500 cm −3 that are ≳5 × larger than that of lower-O32 galaxies. We do not find a significant difference in 〈 Q 〉 between low- and high-O32 galaxies. Photoionization modeling indicates a large spread in log U of ≈1.5 dex at a fixed Z neb . On the other hand, the data indicate a highly significant correlation between U and star-formation-rate surface density (Σ SFR ), which appears to be redshift invariant at z ∼ 1.6−6.3, and possibly up to z ∼ 9.5. We consider several avenues through which metallicity and Σ SFR (or gas density) may influence U , including variations in n e and Q , internal dust extinction of ionizing photons, and the effects of gas density on the volume filling fraction. Based on these considerations, we conclude that gas density may play a more central role than metallicity in modulating U at these redshifts.

Abstract: We present a multiwavelength analysis using the Submillimeter Array (SMA), James Clerk Maxwell Telescope, NOEMA, JWST, the Hubble Space Telescope (HST), and the Spitzer Space Telescope of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A 1489. These SMA-located sources both lie at z = 4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at ≲ 2 μm, even with deep JWST/NIRCam observations. We investigate their stellar, interstellar medium, and dynamical properties, including a pixel-level spectral energy distribution analysis to derive subkiloparsec-resolution stellar-mass and A V maps. We find that 850.1 is one of the most massive and highly obscured, A V ∼ 5, galaxies known at z > 4 with M * ∼1011.8 M ⊙ (likely forming at z > 6), and 850.2 is one of the least massive and least obscured, A V ∼ 1, members of the z > 4 dusty star-forming population. The diversity of these two dust-mass-selected galaxies illustrates the incompleteness of galaxy surveys at z ≳ 3–4 based on imaging at ≲ 2 μm, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar-mass distribution, R e mass ∼1 kpc, but its expected evolution means that it matches both the properties of massive, quiescent galaxies at z ∼ 1.5 and ultramassive early-type galaxies at z ∼ 0. We suggest that 850.1 is the central galaxy of a group in which 850.2 is a satellite that will likely merge in the near future. The stellar morphology of 850.1 shows arms and a linear bar feature that we link to the active dynamical environment it resides within.

Abstract: Abstract The mean free path of ionizing photons, λ mfp , is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of λ mfp from QSO spectra over the redshift range 5 < z < 6, including the first measurements at z ≃ 5.3 and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI observations of QSOs near z ∼ 5.5, for which we also acquire new [C ii ] 158 μ m redshifts with ALMA. By measuring the Lyman continuum transmission profile in stacked QSO spectra, we find λ mfp = 9.33 − 1.80 + 2.06 , 5.40 − 1.40 + 1.47 , 3.31 − 1.34 + 2.74 , and 0.81 − 0.48 + 0.73 pMpc at z = 5.08, 5.31, 5.65, and 5.93, respectively. Our results demonstrate that λ mfp increases steadily and rapidly with time over 5 < z < 6. Notably, we find that λ mfp deviates significantly from predictions based on a fully ionized and relaxed IGM as late as z = 5.3. By comparing our results to model predictions and indirect λ mfp constraints based on IGM Ly α opacity, we find that the evolution of λ mfp is consistent with scenarios wherein the IGM is still undergoing reionization and/or retains large fluctuations in the ionizing UV background well below redshift 6.

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

Abstract: Dust particle sizes constrained from dust continuum and polarization observations by radio interferometry are inconsistent by at least an order of magnitude. Motivated by porous dust observed in small Solar System bodies (e.g., from the Rosetta mission), we explore how the dust particle's porosity affects the estimated particle sizes from these two methods. Porous particles have lower refractive indices, which affect both opacity and polarization fraction. With weaker Mie interference patterns, the porous particles have lower opacity at mm wavelengths than the compact particles if the particle size exceeds several hundred microns. Consequently, the inferred dust mass using porous particles can be up to a factor of six higher. The most significant difference between compact and porous particles is their scattering properties. The porous particles have a wider range of particle sizes with high linear polarization from dust self-scattering, allowing mm-cm-sized particles to explain polarization observations. With a Bayesian approach, we use porous particles to fit HL Tau disk's multi-wavelength continuum and mm-polarization observations from ALMA and VLA. The moderately porous particles with sizes from 1 mm-1 m can explain both continuum and polarization observations, especially in the region between 20-60 au. If the particles in HL Tau are porous, the porosity should be from 70% to 97% from current polarization observations. We also predict that future observations of the self-scattering linear polarization at longer wavelengths (e.g., ALMA B1 and ngVLA) have the potential to further constrain the particle's porosity and size.

Abstract: Abstract We present a comprehensive study of 29 short gamma-ray bursts (SGRBs) observed ≈0.8−60 days postburst using Chandra and XMM-Newton. We provide the inferred distributions of the SGRB jet opening angles and true event rates to compare against neutron star merger rates. We perform a uniform analysis and modeling of their afterglows, obtaining 10 opening angle measurements and 19 lower limits. We report on two new opening angle measurements (SGRBs 050724A and 200411A) and eight updated values, obtaining a median value of 〈 θ j 〉 ≈ 6.°1 [−3.°2, +9.°3] (68% confidence on the full distribution) from jet measurements alone. For the remaining events, we infer θ j ≳ 0.°5–26°. We uncover a population of SGRBs with wider jets of θ j ≳ 10° (including two measurements of θ j ≳ 15°), representing ∼28% of our sample. Coupled with multiwavelength afterglow information, we derive a total true energy of 〈 E true,tot 〉 ≈ 10 49 –10 50 erg, which is consistent with magnetohydrodynamic jet launching mechanisms. Furthermore, we determine a range for the beaming-corrected event rate of R true ≈ 360 − 1800 Gpc −3 yr −1 , set by the inclusion of a population of wide jets on the low end, and the jet measurements alone on the high end. From a comparison with the latest merger rates, our results are consistent with the majority of SGRBs originating from binary neutron star mergers. However, our inferred rates are well above the latest neutron star–black hole merger rates, consistent with at most a small fraction of SGRBs originating from such mergers.

Abstract: Abstract We characterize the stellar mass of J2239+0207, a z ∼ 6.25 sub-Eddington quasar ( M 1450 = −24.6), using dedicated JWST/NIRCam medium-band observations of a nearby point-spread function star to remove the central point source and reveal the underlying galaxy emission. We detect the host galaxy in two bands longward of the Balmer break, obtaining a stellar mass of ∼10 10 M ⊙ , more than an order of magnitude less than this quasar’s existing measured [C ii ] dynamical mass. We additionally calculate the mass of J2239+0207's central supermassive black hole using JWST/NIRSpec integral field unit observations, and determine that the black hole is ∼15 times more massive than predicted by the local M BH – M * relation, similar to many high-redshift quasars with dynamical masses determined via millimeter-wave line widths. We carefully consider potential selection effects at play, and find that even when z ∼ 6 quasars are compared to a local sample with similarly determined dynamical masses, many of the high-redshift quasars appear to possess overmassive black holes. We conclude z ∼ 6 quasars are likely to have a larger spread about the M BH – M * relation than observed in the local Universe.

Abstract: We derive the bolometric luminosities (L bol) of 865 field-age and 189 young ultracool dwarfs (spectral types M6–T9, including 40 new discoveries presented here) by directly integrating flux-calibrated optical to mid-infrared (MIR) spectral energy distributions (SEDs). The SEDs consist of low-resolution (R ∼ 150) near-infrared (NIR; 0.8–2.5μm) spectra (including new spectra for 97 objects), optical photometry from the Pan-STARRS1 survey, and MIR photometry from the CatWISE2020 survey and Spitzer/IRAC. Our L bol calculations benefit from recent advances in parallaxes from Gaia, Spitzer, and UKIRT, as well as new parallaxes for 19 objects from CFHT and Pan-STARRS1 presented here. Coupling our L bol measurements with a new uniform age analysis for all objects, we estimate substellar masses, radii, surface gravities, and effective temperatures (T eff) using evolutionary models. We construct empirical relationships for L bol and T eff as functions of spectral type and absolute magnitude, determine bolometric corrections in optical and infrared bandpasses, and study the correlation between evolutionary model-derived surface gravities and NIR gravity classes. Our sample enables a detailed characterization of BT-Settl and ATMO 2020 atmospheric model systematics as a function of spectral type and position in the NIR color–magnitude diagram. We find the greatest discrepancies between atmospheric and evolutionary model-derived T eff (up to 800 K) and radii (up to 2.0 R Jup) at the M/L spectral type transition boundary. With 1054 objects, this work constitutes the largest sample to date of ultracool dwarfs with determinations of their fundamental parameters.

Abstract: Abstract We study current bounds on strong first-order phase transitions (PTs) along the equation of state (EOS) of dense strongly interacting matter in neutron stars, under the simplifying assumption that on either side of the PT, the EOS can be approximated by a simple polytropic form. We construct a large ensemble of possible EOSs of this form, anchor them to chiral effective field theory calculations at nuclear density and perturbative Quantum Chromodynamics at high densities, and subject them to astrophysical constraints from high-mass pulsars and gravitational-wave observations. Within this setup, we find that a PT permits neutron-star solutions with larger radii, but only if the transition begins below twice nuclear saturation density. We also identify a large parameter space of allowed PTs currently unexplored by numerical-relativity studies. Additionally, we locate a small region of parameter space allowing twin-star solutions, though we find them to only marginally pass the current astrophysical constraints. Finally, we find that sizeable cores of high-density matter beyond the PT may be located in the centers of some stable neutron stars, primarily those with larger masses.

Abstract: The carbon (3P) atom is a reactive species that, according to laboratory experiments and theoretical calculations, condensates with interstellar ice components. This fact is of uttermost importance for the chemistry in the interstellar medium (ISM) because the condensation reaction is barrierless, and the subsequent species formed are still reactive given their open-shell character. Carbon condensation on CO-rich ices forms the C=C=O (3Σ−) species, which can be easily hydrogenated twice to form ketene (H2CCO). Ketene is very reactive in terrestrial conditions, usually found as an intermediate that is difficult to isolate in chemical synthesis laboratories. These characteristics suggest that ketene can be a good candidate to form interstellar complex organic molecules via a two-step process, i.e., its activation followed by a radical–radical coupling. In this work, reactions between ketene and atomic H and the OH and NH2 radicals on a CO-rich ice model have been explored by means of quantum chemical calculations complemented by kinetic calculations to evaluate if they are favorable in the ISM. Results indicate that the addition of H to ketene (helped by tunneling) to form the acetyl radical (CH3CO) is the most preferred path as the reactions with OH and NH2 possess activation energies (≥9 kJ mol−1) hard to surmount in the ISM conditions unless external processes provide energy to the system. Thus, acetaldehyde (CH3CHO) and, probably, ethanol (CH3CH2OH) formation via further hydrogenations, are the possible unique operating synthetic routes. Moreover, from the computed, relatively large binding energies of OH and NH2 on CO ice, slow diffusion is expected, hampering possible radical–radical couplings with CH3CO. The astrophysical implications of these findings are discussed considering the incoming James Webb Space Telescope observations.

Abstract: Understanding the coevolution of supermassive black holes and their host systems requires a comprehensive census of active galactic nuclei (AGNs) behavior across a wide range of redshift, luminosity, obscuration level, and galaxy properties. We report significant progress with JWST toward this goal from the Systematic Mid-infrared Instrument Legacy Extragalactic Survey (SMILES). Based on comprehensive spectral energy distribution (SED) analysis of 3273 MIRI-detected sources, we identify 217 AGN candidates over a survey area of ∼34 arcmin2, including a primary sample of 111 AGNs in normal massive galaxies (M * > 109.5 M ☉) at z ∼ 0–4, an extended sample of 86 AGN candidates in low-mass galaxies (M * < 109.5 M ☉), and a high-z sample of 20 AGN candidates at z ∼ 4–8.4. Notably, about 80% of our MIRI-selected AGN candidates are new discoveries despite the extensive pre-JWST AGN searches. Even among the massive galaxies where the previous AGN search is believed to be thorough, 34% of the MIRI AGN identifications are new, highlighting the impact of obscuration on previous selections. By combining our results with the efforts at other wavelengths, we build the most complete AGN sample to date and examine the relative performance of different selection techniques. We find the obscured AGN fraction increases from L AGN,bol ∼ 1010 L ⊙ to 1011 L ⊙ and then drops toward higher luminosity. Additionally, the obscured AGN fraction gradually increases from z ∼ 0 to z ∼ 4 with most high-z AGNs obscured. We discuss how AGN obscuration, intrinsic SED variations, galaxy contamination, survey depth, and selection techniques complicate the construction of a complete AGN sample.

Abstract: Solar radio bursts are strongly affected by radio-wave scattering on density inhomogeneities, changing their observed time characteristics, sizes, and positions. The same turbulence causes angular broadening and scintillation of galactic and extragalactic compact radio sources observed through the solar atmosphere. Using large-scale simulations of radio-wave transport, the characteristics of anisotropic density turbulence from 0.1 R ⊙ to 1 au are explored. For the first time, a profile of heliospheric density fluctuations is deduced that accounts for the properties of extrasolar radio sources, solar radio bursts, and in situ density fluctuation measurements in the solar wind at 1 au. The radial profile of the spectrum-weighted mean wavenumber of density fluctuations (a quantity proportional to the scattering rate of radio waves) is found to have a broad maximum at around (4–7) R ⊙, where the slow solar wind becomes supersonic. The level of density fluctuations at the inner scale (which is consistent with the proton resonance scale) decreases with heliocentric distance as 〈δni2〉(r)≃2×107r/R⊙−1−3.7 cm−6. Due to scattering, the apparent positions of solar burst sources observed at frequencies between 0.1 and 300 MHz are computed to be essentially cospatial and to have comparable sizes, for both fundamental and harmonic emission. Anisotropic scattering is found to account for the shortest solar radio burst decay times observed, and the required wavenumber anisotropy is q ∥/q ⊥ = 0.25–0.4, depending on whether fundamental or harmonic emission is involved. The deduced radio-wave scattering rate paves the way to quantify intrinsic solar radio burst characteristics.

Abstract: Abstract We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), α -element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low- α disk, high- α disk, and total population independently. The low- α disk exhibits a negative radial metallicity gradient of −0.06 ± 0.001 dex kpc −1 , which flattens with distance from the midplane. The high- α disk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]–[Fe/H] plane and in the [Mg/Fe]–age relation persist across the entire disk. The age estimates have typical uncertainties of ∼0.15 in log(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy.

Abstract: The observation of gravitational waves from multiple compact binary coalescences by the LIGO–Virgo–KAGRA detector networks has enabled us to infer the underlying distribution of compact binaries across a wide range of masses, spins, and redshifts. In light of the new features found in the mass spectrum of binary black holes and the uncertainty regarding binary formation models, nonparametric population inference has become increasingly popular. In this work, we develop a data-driven clustering framework that can identify features in the component mass distribution of compact binaries simultaneously with those in the corresponding redshift distribution, from gravitational-wave data in the presence of significant measurement uncertainties, while making very few assumptions about the functional form of these distributions. Our generalized model is capable of inferring correlations among various population properties, such as the redshift evolution of the shape of the mass distribution itself, in contrast to most existing nonparametric inference schemes. We test our model on simulated data and demonstrate the accuracy with which it can reconstruct the underlying distributions of component masses and redshifts. We also reanalyze public LIGO–Virgo–KAGRA data from events in GWTC-3 using our model and compare our results with those from some alternative parametric and nonparametric population inference approaches. Finally, we investigate the potential presence of correlations between mass and redshift in the population of binary black holes in GWTC-3 (those observed by the LIGO–Virgo–KAGRA detector network in their first three observing runs), without making any assumptions about the specific nature of these correlations.

Abstract: The gravitational three-body problem is a fundamental problem in physics and has significant applications to astronomy. Three-body configurations are often considered stable as long the system is hierarchical; that is, the two orbital distances are well-separated. However, instability, which is often associated with significant energy exchange between orbits, takes time to develop. Assuming two massive objects in a circular orbit and a test particle in an eccentric orbit, we develop an analytical formula estimating the time it takes for the test particle’s orbital energy to change by an order of itself. We show its consistency with results from N-body simulations. For eccentric orbits in particular, the instability is primarily driven not by close encounters of the test particle with one of the other bodies, but by the fundamental susceptibility of eccentric orbits to exchange energy at their periapsis. Motivated by recent suggestions that the galactic center may host an intermediate-mass black hole (IMBH) as a companion to the massive black hole Sgr A*, we use our timescale to explore the parameter space that could harbor an IMBH for the lifetime of the S-cluster of stars surrounding Sgr A*. Furthermore, we show that the orbit of an S-star can be stable for long timescales in the presence of other orbital crossing stars, thus suggesting that the S-cluster may be stable for the lifetimes of its member stars.

Abstract: Luminous quasars are powerful targets to investigate the role of feedback from supermassive black holes (BHs) in regulating the growth phases of BHs themselves and of their host galaxies, up to the highest redshifts. Here we investigate the cosmic evolution of the occurrence and kinematics of BH-driven outflows, as traced by broad absorption line (BAL) features, due to the C iv ionic transition. We exploit a sample of 1935 quasars at z = 2.1–6.6 with bolometric luminosity log(L bol/erg s−1) ≳ 46.5, drawn from the Sloan Digital Sky Survey and from the X-Shooter legacy survey of Quasars at the Reionization Epoch (XQR-30). We consider rest-frame optical bright quasars to minimize observational biases due to quasar selection criteria. We apply a homogeneous BAL-identification analysis, based on employing composite template spectra to estimate the quasar intrinsic emission. We find a BAL quasar fraction close to 20% at z ∼ 2–4, while it increases to almost 50% at z ∼ 6. The velocity and width of the BAL features also increase at z ≳ 4.5. We exclude the possibility that the redshift evolution of the BAL properties is due to differences in terms of quasar luminosity and accretion rate. These results suggest significant BH feedback occurring in the 1 Gyr old universe, likely affecting the growth of BHs and, possibly, of their host galaxies, as supported by models of early BH and galaxy evolution.

Abstract: Abstract Brown dwarf spectra offer vital testbeds for our understanding of the chemical and physical processes that sculpt substellar atmospheres. Recently, atmospheric retrieval approaches have been successfully applied to low-resolution ( R ∼ 100) spectra of L, T, and Y dwarfs, yielding constraints on the chemical abundances and temperature structures of these atmospheres. Medium-resolution ( R ∼ 10 3 ) spectra of brown dwarfs offer additional insight, as molecular features are more easily disentangled and the thermal structure of the upper atmosphere is better probed. We present results from a GPU-based retrieval analysis of a high signal-to-noise, medium-resolution ( R ∼ 6000) FIRE spectrum from 0.85 to 2.5 μ m of the T9 dwarf UGPS J072227.51–054031.2. At 60× higher spectral resolution than previous brown dwarf retrievals, a number of novel challenges arise. We examine the effect of different opacity sources, in particular for CH 4 . Furthermore, we find that flaws in the data like errors from order stitching can bias our constraints. We compare these retrieval results to those for an R ∼ 100 spectrum of the same object, revealing how constraints on atmospheric abundances and temperatures improve by an order of magnitude or more with increased spectral resolution. In particular, we can constrain the abundance of H 2 S, which is undetectable at lower spectral resolution. While these medium-resolution retrievals offer the potential of precise, stellar-like constraints on atmospheric abundances (∼0.02 dex), our retrieved radius is unphysically small ( R = 0.50 − 0.01 + 0.01 R Jup ), indicating shortcomings with our modeling framework. This work is an initial investigation into brown dwarf retrievals at medium spectral resolution, offering guidance for future ground-based studies and JWST observations.