The 100 most-cited articles in The Astrophysical Journal, 2018

Abstract We present optical light curves, redshifts, and classifications for spectroscopically confirmed Type Ia supernovae (SNe Ia) discovered by the Pan-STARRS1 (PS1) Medium Deep Survey. We detail improvements to the PS1 SN photometry, astrometry, and calibration that reduce the systematic uncertainties in the PS1 SN Ia distances. We combine the subset of PS1 SNe Ia (0.03 < z < 0.68) with useful distance estimates of SNe Ia from the Sloan Digital Sky Survey (SDSS), SNLS, and various low- z and Hubble Space Telescope samples to form the largest combined sample of SNe Ia, consisting of a total of SNe Ia in the range of 0.01 < z < 2.3, which we call the “Pantheon Sample.” When combining Planck 2015 cosmic microwave background (CMB) measurements with the Pantheon SN sample, we find and for the w CDM model. When the SN and CMB constraints are combined with constraints from BAO and local H 0 measurements, the analysis yields the most precise measurement of dark energy to date: and for the CDM model. Tension with a cosmological constant previously seen in an analysis of PS1 and low- z…

Abstract This paper provides an update of our previous scaling relations between galaxy-integrated molecular gas masses, stellar masses, and star formation rates (SFRs), in the framework of the star formation main sequence (MS), with the main goal of testing for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ∼1 mm dust photometry, in a large sample of 1444 star-forming galaxies between z = 0 and 4. The sample covers the stellar mass range log( M * / M ⊙ ) = 9.0–11.8, and SFRs relative to that on the MS, δ MS = SFR/SFR(MS), from 10 −1.3 to 10 2.2 . Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero-point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t depl , defined as the ratio of molecular gas mass to SFR, scales as (1 + z ) −0.6 × ( δ MS) −0.44 and is…

Abstract We present Hubble Space Telescope (HST) photometry of a selected sample of 50 long-period, low-extinction Milky Way Cepheids measured on the same WFC3 F555W-, F814W-, and F160W-band photometric system as extragalactic Cepheids in Type Ia supernova host galaxies. These bright Cepheids were observed with the WFC3 spatial scanning mode in the optical and near-infrared to mitigate saturation and reduce pixel-to-pixel calibration errors to reach a mean photometric error of 5 mmag per observation. We use the new Gaia DR2 parallaxes and HST photometry to simultaneously constrain the cosmic distance scale and to measure the DR2 parallax zeropoint offset appropriate for Cepheids. We find the latter to be −46 ± 13 μas or ±6 μas for a fixed distance scale, higher than found from quasars, as expected for these brighter and redder sources. The precision of the distance scale from DR2 has been reduced by a factor of 2.5 because of the need to independently determine the parallax offset. The best-fit distance scale is 1.006 ± 0.033, relative to the scale from Riess et al. with H 0 = 73.24…

Abstract We study the dust attenuation curves of 230,000 individual galaxies in the local universe, ranging from quiescent to intensely star-forming systems, using GALEX , SDSS, and WISE photometry calibrated on the Herschel ATLAS. We use a new method of constraining SED fits with infrared luminosity (SED+LIR fitting), and parameterized attenuation curves determined with the CIGALE SED-fitting code. Attenuation curve slopes and UV bump strengths are reasonably well constrained independently from one another. We find that attenuation curves exhibit a very wide range of slopes that are on average as steep as the curve slope of the Small Magellanic Cloud (SMC). The slope is a strong function of optical opacity. Opaque galaxies have shallower curves—in agreement with recent radiative transfer models. The dependence of slopes on the opacity produces an apparent dependence on stellar mass: more massive galaxies have shallower slopes. Attenuation curves exhibit a wide range of UV bump amplitudes, from none to Milky Way (MW)-like, with an average strength one-third that of the MW bump.…

Abstract We report a reanalysis of a near-pristine absorption system, located at a redshift toward the quasar Q1243+307, based on the combination of archival and new data obtained with the HIRES echelle spectrograph on the Keck telescope. This absorption system, which has an oxygen abundance [O/H] = −2.769 ± 0.028 (≃1/600 of the solar abundance), is among the lowest metallicity systems currently known where a precise measurement of the deuterium abundance is afforded. Our detailed analysis of this system concludes, on the basis of eight D i absorption lines, that the deuterium abundance of this gas cloud is , which is in very good agreement with the results previously reported by Kirkman et al., but with an improvement on the precision of this single measurement by a factor of ∼3.5. Combining this new estimate with our previous sample of six high precision and homogeneously analyzed D/H measurements, we deduce that the primordial deuterium abundance is or, expressed as a linear quantity, this value corresponds to a one percent determination of the primordial deuterium abundance.…

Abstract BANYAN Σ is a new Bayesian algorithm to identify members of young stellar associations within 150 pc of the Sun. It includes 27 young associations with ages in the range ∼1–800 Myr, modeled with multivariate Gaussians in six-dimensional (6D) XYZUVW space. It is the first such multi-association classification tool to include the nearest sub-groups of the Sco-Cen OB star-forming region, the IC 2602, IC 2391, Pleiades and Platais 8 clusters, and the ρ Ophiuchi, Corona Australis, and Taurus star formation regions. A model of field stars is built from a mixture of multivariate Gaussians based on the Besançon Galactic model. The algorithm can derive membership probabilities for objects with only sky coordinates and proper motion, but can also include parallax and radial velocity measurements, as well as spectrophotometric distance constraints from sequences in color–magnitude or spectral type–magnitude diagrams. BANYAN Σ benefits from an analytical solution to the Bayesian marginalization integrals over unknown radial velocities and distances that makes it more accurate and…

Abstract We present new measurements of the parallax of seven long-period (≥10 days) Milky Way (MW) Cepheid variables (SS CMa, XY Car, VY Car, VX Per, WZ Sgr, X Pup, and S Vul) using one-dimensional astrometric measurements from spatial scanning of Wide-Field Camera 3 on the Hubble Space Telescope (HST). The observations were obtained at ∼6 month intervals over 4 years. The distances are 1.7–3.6 kpc, with a mean precision of 45 μas (signal-to-noise ratio (S/N) ≈ 10) and a best precision of 29 μas (S/N = 14). The accuracy of the parallaxes is demonstrated through independent analyses of >100 reference stars. This raises to 10 the number of long-period Cepheids with significant parallax measurements, 8 obtained from this program. We also present high-precision mean F555W, F814W, and F160W magnitudes of these Cepheids, allowing a direct, zeropoint-independent comparison to >1800 extragalactic Cepheids in the hosts of 19 SNe Ia. This sample addresses two outstanding systematic uncertainties affecting prior comparisons of MW and extragalactic Cepheids used to calibrate the Hubble…

Abstract Rings are the most frequently revealed substructure in Atacama Large Millimeter/submillimeter Array (ALMA) dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution (∼0.″12) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high-contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between…

Abstract We present the catalog of ∼31,500 extragalactic H i line sources detected by the completed Arecibo Legacy Fast ALFA (ALFALFA) survey out to z < 0.06, including both high signal-to-noise ratio (>6.5) detections and ones of lower quality that coincide in both position and recessional velocity with galaxies of known redshift. We review the observing technique, data reduction pipeline, and catalog construction process, focusing on details of particular relevance to understanding the catalog’s compiled parameters. We further describe and make available the digital H i line spectra associated with the cataloged sources. In addition to the extragalactic H i line detections, we report nine confirmed OH megamasers (OHMs) and 10 OHM candidates at 0.16 < z < 0.22 whose OH line signals are redshifted into the ALFALFA frequency band. Because of complexities in data collection and processing associated with the use of a feed-horn array on a complex single-dish antenna in the terrestrial radio frequency interference environment, we also present a list of suggestions and caveats for…

Abstract We present a systematic numerical relativity study of the mass ejection and the associated electromagnetic transients and nucleosynthesis from binary neutron star (NS) mergers. We find that a few 10 −3 M ⊙ of material is ejected dynamically during the mergers. The amount and the properties of these outflows depend on binary parameters and on the NS equation of state (EOS). A small fraction of these ejecta, typically ∼10 −6 M ⊙ , is accelerated by shocks formed shortly after merger to velocities larger than 0.6 c and produces bright radio flares on timescales of weeks, months, or years after merger. Their observation could constrain the strength with which the NSs bounce after merger and, consequently, the EOS of matter at extreme densities. The dynamical ejecta robustly produce second and third r -process peak nuclei with relative isotopic abundances close to solar. The production of light r -process elements is instead sensitive to the binary mass ratio and the neutrino radiation treatment. Accretion disks of up to ∼0.2 M ⊙ are formed after merger, depending on the…

Abstract Transmission spectra are differential measurements that utilize stellar illumination to probe transiting exoplanet atmospheres. Any spectral difference between the illuminating light source and the disk-integrated stellar spectrum due to starspots and faculae will be imprinted in the observed transmission spectrum. However, few constraints exist for the extent of photospheric heterogeneities in M dwarfs. Here we model spot and faculae covering fractions consistent with observed photometric variabilities for M dwarfs and the associated 0.3–5.5 μ m stellar contamination spectra. We find that large ranges of spot and faculae covering fractions are consistent with observations and corrections assuming a linear relation between variability amplitude, and covering fractions generally underestimate the stellar contamination. Using realistic estimates for spot and faculae covering fractions, we find that stellar contamination can be more than 10× larger than the transit depth changes expected for atmospheric features in rocky exoplanets. We also find that stellar spectral…

Abstract We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released 11 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no evidence for a GWB, we place constraints on a population of inspiraling supermassive black hole (SMBH) binaries, a network of decaying cosmic strings, and a primordial GWB. For the first time, we find that the GWB constraints are sensitive to the solar system ephemeris (SSE) model used and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first pulsar-timing array (PTA) constraints that are robust against SSE errors. We thus place a 95% upper limit on the GW-strain amplitude of A GWB < 1.45 × 10−15 at a frequency of f = 1 yr−1 for a fiducial f −2/3 power-law spectrum and with interpulsar correlations modeled. This is a factor of ∼2 improvement over the NANOGrav nine-year limit calculated using the same procedure. Previous PTA upper limits on the GWB (as well as their astrophysical and cosmological…

Abstract The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a novel transit radio telescope operating across the 400–800 MHz band. CHIME is composed of four 20 m × 100 m semicylindrical paraboloid reflectors, each of which has 256 dual-polarization feeds suspended along its axis, giving it a ≳200 deg 2 field of view. This, combined with wide bandwidth, high sensitivity, and a powerful correlator, makes CHIME an excellent instrument for the detection of fast radio bursts (FRBs). The CHIME Fast Radio Burst Project (CHIME/FRB) will search beam-formed, high time and frequency resolution data in real time for FRBs in the CHIME field of view. Here we describe the CHIME/FRB back end, including the real-time FRB search and detection software pipeline, as well as the planned offline analyses. We estimate a CHIME/FRB detection rate of 2–42 FRBs sky –1 day –1 normalizing to the rate estimated at 1.4 GHz by Vander Wiel et al. Likely science outcomes of CHIME/FRB are also discussed. CHIME/FRB is currently operational in a commissioning phase, with science operations expected to…

Abstract We present an analysis of all prime HST legacy fields spanning >800 arcmin2 in the search for z ∼ 10 galaxy candidates and the study of their UV luminosity function (LF). In particular, we present new z ∼ 10 candidates selected from the full Hubble Frontier Field (HFF) data set. Despite the addition of these new fields, we find a low abundance of z ∼ 10 candidates with only nine reliable sources identified in all prime HST data sets that include the HUDF09/12, the HUDF/XDF, all of the CANDELS fields, and now the HFF survey. Based on this comprehensive search, we find that the UV luminosity function decreases by one order of magnitude from z ∼ 8 to z ∼ 10 over a four-magnitude range. This also implies a decrease of the cosmic star formation rate density by an order of magnitude within 170 Myr from z ∼ 8 to z ∼ 10. We show that this accelerated evolution compared to lower redshift can entirely be explained by the fast build up of the dark matter halo mass function at z > 8. Consequently, the predicted UV LFs from several models of galaxy formation are in good agreement with…

Abstract The dense matter equation of state (EOS) determines neutron star (NS) structure but can be calculated reliably only up to one to two times the nuclear saturation density, using accurate many-body methods that employ nuclear interactions from chiral effective field theory constrained by scattering data. In this work, we use physically motivated ansatzes for the speed of sound c S at high density to extend microscopic calculations of neutron-rich matter to the highest densities encountered in stable NS cores. We show how existing and expected astrophysical constraints on NS masses and radii from X-ray observations can constrain the speed of sound in the NS core. We confirm earlier expectations that c S is likely to violate the conformal limit of , possibly reaching values closer to the speed of light c at a few times the nuclear saturation density, independent of the nuclear Hamiltonian. If QCD obeys the conformal limit, we conclude that the rapid increase of c S required to accommodate a 2 M ⊙ NS suggests a form of strongly interacting matter where a description in terms of…

Abstract The central radio source in M87 provides the best opportunity to study jet formation because it has a large angular size for the gravitational radius of the black hole and has a bright jet that is well resolved by very long baseline interferometry observations. We present intensive monitoring observations from 2007 and 2008, plus roughly annual observations that span 17 years, all made with the the Very Long Baseline Array at 43 GHz with a resolution of about 30 by 60 R S . Our high dynamic range images clearly show the wide opening angle structure and the counterjet. The jet and counterjet are nearly symmetric in the inner 1.5 mas (0.12 pc in projection), with both being edge brightened. Both show deviations from parabolic shape in the form of an initial rapid expansion and subsequent contraction followed by further rapid expansion and, beyond the visible counterjet, subsequent collimation. Proper motions and counterjet/jet intensity ratios both indicate acceleration from apparent speeds of ≲0.5 c to ≳2 c in the inner ∼2 mas and suggest a helical flow. The jet displays a…

Abstract We present Atacama Large Millimeter/Sub-Millimeter Array (ALMA) Band 6 observations of a complete sample of protoplanetary disks in the young (∼1–3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12 CO, 13 CO, and C 18 O J = 2–1 lines. The spatial resolution is ∼0.″25 with a medium 3 σ continuum sensitivity of 0.30 mJy, corresponding to M dust ∼ 0.2 M ⊕ . We apply Keplerian masking to enhance the signal-to-noise ratios of our 12 CO zero-moment maps, enabling measurements of gas disk radii for 22 Lupus disks; we find that gas disks are universally larger than millimeter dust disks by a factor of two on average, likely due to a combination of the optically thick gas emission and the growth and inward drift of the dust. Using the gas disk radii, we calculate the dimensionless viscosity parameter, α visc , finding a broad distribution and no correlations with other disk or stellar parameters, suggesting that viscous processes have not yet established quasi-steady states in Lupus disks. By combining our 1.33 mm continuum fluxes with our previous 890 μ m…

Abstract We present a new flexible Bayesian framework for directly inferring the fraction of neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR, z ∼ 6–10) from detections and non-detections of Lyman Alpha (Ly α ) emission from Lyman Break galaxies (LBGs). Our framework combines sophisticated reionization simulations with empirical models of the interstellar medium (ISM) radiative transfer effects on Ly α . We assert that the Ly α line profile emerging from the ISM has an important impact on the resulting transmission of photons through the IGM, and that these line profiles depend on galaxy properties. We model this effect by considering the peak velocity offset of Ly α lines from host galaxies’ systemic redshifts, which are empirically correlated with UV luminosity and redshift (or halo mass at fixed redshift). We use our framework on the sample of LBGs presented in Pentericci et al. and infer a global neutral fraction at z ∼ 7 of , consistent with other robust probes of the EoR and confirming that reionization is ongoing ∼700 Myr after the Big…

Abstract We present the first part of our Disks ARound T Tauri Stars with SPHERE (DARTTS-S) survey: observations of eight T Tauri stars that were selected based on their strong (sub)millimeter excesses using SPHERE/IRDIS polarimetric differential imaging in the J and H bands. All observations successfully detect the disks, which appear vastly different in size, from ≈80 au in scattered light to >400 au, and display total polarized disk fluxes between 0.06% and 0.89% of the stellar flux. For five of these disks, we are able to determine the three-dimensional structure and the flaring of the disk surface, which appears to be relatively consistent across the different disks, with flaring exponents α between ≈1.1 and ≈1.6. We also confirm literature results with regard to the inclination and position angle of several of our disks and are able to determine which side is the near side of the disk in most cases. While there is a clear trend of disk mass with stellar ages (≈1 to >10 Myr), no correlations of disk structures with age were found. There are also no correlations with either…

Abstract We investigate the nature of the double color–magnitude sequence observed in the Gaia DR2 HR diagram of stars with high transverse velocities. The stars in the reddest-color sequence are likely dominated by the dynamically hot tail of the thick disk population. Information from Nissen & Schuster and from the APOGEE survey suggests that stars in the blue-color sequence have elemental abundance patterns that can be explained by this population having a relatively low star formation efficiency during its formation. In dynamical and orbital spaces, such as the “Toomre diagram,” the two sequences show a significant overlap, but with a tendency for stars on the blue-color sequence to dominate regions with no or retrograde rotation and high total orbital energy. In the plane defined by the maximal vertical excursion of the orbits versus their apocenters, stars of both sequences redistribute into discrete wedges. We conclude that stars that are typically assigned to the halo in the solar vicinity are actually both accreted stars lying along the blue sequence in the HR diagram, and…

Abstract We present an updated release of the BaSTI (a Bag of Stellar Tracks and Isochrones) stellar model and isochrone library for a solar-scaled heavy element distribution. The main input physics that have been changed from the previous BaSTI release include the solar metal mixture, electron conduction opacities, a few nuclear reaction rates, bolometric corrections, and the treatment of the overshooting efficiency for shrinking convective cores. The new model calculations cover a mass range between 0.1 and 15 M ⊙ , 22 initial chemical compositions between [Fe/H] = −3.20 and +0.45, with helium to metal enrichment ratio dY / dZ = 1.31. The isochrones cover an age range between 20 Myr and 14.5 Gyr, consistently take into account the pre-main-sequence phase, and have been translated to a large number of popular photometric systems. Asteroseismic properties of the theoretical models have also been calculated. We compare our isochrones with results from independent databases and with several sets of observations to test the accuracy of the calculations. All stellar evolution tracks,…

We present the initial–final mass relation (IFMR) based on the self-consistent analysis of Sirius B and 79 white dwarfs from 13 star clusters. We have also acquired additional signal on eight white dwarfs previously analyzed in the NGC 2099 cluster field, four of which are consistent with membership. These re-observed white dwarfs have masses ranging from 0.72 to 0.97 M ⊙ , with initial masses from 3.0 to 3.65 M ⊙ , where the IFMR has an important change in slope that these new data help to observationally confirm. In total, this directly measured IFMR has small scatter ( σ = 0.06 M ⊙ ) and spans from progenitors of 0.85 to 7.5 M ⊙ . Applying two different stellar evolutionary models to infer two different sets of white dwarf progenitor masses shows that, when the same model is also used to derive the cluster ages, the resulting IFMR has weak sensitivity to the adopted model at all but the highest initial masses (>5.5 M ⊙ ). The nonlinearity of the IFMR is also clearly observed with moderate slopes at lower masses (0.08 M final / M initial ) and higher masses (0.11 M final / M…

Abstract We compile a complete collection of reliable Hubble parameter H ( z ) data to redshift z ≤ 2.36 and use them with the Gaussian Process method to determine continuous H ( z ) functions for various data subsets. From these continuous H ( z )'s, summarizing across the data subsets considered, we find H 0 ∼ 67 ± 4 km s −1 Mpc −1 , more consistent with the recent lower values determined using a variety of techniques. In most data subsets, we see a cosmological deceleration–acceleration transition at 2 σ significance, with the data subsets transition redshifts varying over at 1 σ significance. We find that the flat-ΛCDM model is consistent with the H ( z ) data to a z of 1.5 to 2.0, depending on data subset considered, with 2 σ deviations from flat-ΛCDM above this redshift range. Using the continuous H ( z ) with baryon acoustic oscillation distance-redshift observations, we constrain the current spatial curvature density parameter to be , consistent with a flat universe, but the large error bar does not rule out small values of spatial curvature that are now under debate.

Abstract We report the first detections of the repeating fast radio burst source FRB 121102 above 5.2 GHz. Observations were performed using the 4–8 GHz receiver of the Robert C. Byrd Green Bank Telescope with the Breakthrough Listen digital backend. We present the spectral, temporal, and polarization properties of 21 bursts detected within the first 60 minutes of a total of 6 hr of observations. These observations comprise the highest burst density yet reported in the literature, with 18 bursts being detected in the first 30 minutes. A few bursts clearly show temporal sub-structure with distinct spectral properties. These sub-structures superimpose to provide an enhanced peak signal-to-noise ratio at higher trial dispersion measures. Broad features occur in ∼1 GHz wide subbands that typically differ in peak frequency between bursts within the band. Finer-scale structures (∼10–50 MHz) within these bursts are consistent with the structure expected from Galactic diffractive interstellar scintillation. The bursts exhibit nearly 100% linear polarization, and a large average rotation…

Abstract New millisecond pulsars (MSPs) in compact binaries provide a good opportunity to search for the most massive neutron stars. Their main-sequence companion stars are often strongly irradiated by the pulsar, displacing the effective center of light from their barycenter and making mass measurements uncertain. We present a series of optical spectroscopic and photometric observations of PSR J2215+5135, a “redback” binary MSP in a 4.14 hr orbit, and measure a drastic temperature contrast between the dark/cold (T N = 5660 K) and bright/hot (T D = 8080 K) sides of the companion star. We find that the radial velocities depend systematically on the atmospheric absorption lines used to measure them. Namely, the semi-amplitude of the radial velocity curve (RVC) of J2215 measured with magnesium triplet lines is systematically higher than that measured with hydrogen Balmer lines, by 10%. We interpret this as a consequence of strong irradiation, whereby metallic lines dominate the dark side of the companion (which moves faster) and Balmer lines trace its bright (slower) side. Further,…

Abstract We present results of a deep spectroscopic survey quantifying the statistics of the escape of ionizing radiation from star-forming galaxies at z ∼ 3. We measure the ratio of ionizing to non-ionizing UV flux density , where f 900 is the mean flux density evaluated over the range [880, 910] Å. We quantify the emergent ratio of ionizing to non-ionizing UV flux density by analyzing high signal-to-noise ratio composite spectra formed from subsamples with common observed properties and numbers sufficient to reduce the statistical uncertainty in the modeled IGM+CGM correction to obtain precise values of , including a full-sample average = 0.057 ± 0.006. We show that increases monotonically with , inducing an inverse correlation with UV luminosity as a by-product. We fit the composite spectra using stellar spectral synthesis together with models of the ISM in which a fraction f c of the stellar continuum is covered by gas with column density . We show that the composite spectra simultaneously constrain the intrinsic properties of the stars (L 900/L 1500)int along with f c , , , and…

Abstract We reprise the analysis of Stassun & Torres, comparing the parallaxes of the eclipsing binaries reported in that paper to the parallaxes newly reported in the Gaia second data release (DR2). We find evidence for a systematic offset of −82 ± 33 μ as, in the sense of the Gaia parallaxes being too small, for brightnesses ( G ≲ 12) and for distances (0.03–3 kpc) in the ranges spanned by the eclipsing binary sample. The offset does not appear to depend strongly on distance within this range, though there is marginal evidence that the offset increases (becomes slightly more negative) for distances ≳1 kpc, up to the 3 kpc distances probed by the test sample. The offset reported here is consistent with the expectation that global systematics in the Gaia DR2 parallaxes are below 100 μ as.

Abstract One bottleneck for the exploitation of data from the Kepler mission for stellar astrophysics and exoplanet research has been the lack of precise radii and evolutionary states for most of the observed stars. We report revised radii of 177,911 Kepler stars derived by combining parallaxes from the Gaia Data Release 2 with the DR25 Kepler Stellar Properties Catalog. The median radius precision is ≈8%, a typical improvement by a factor of 4–5 over previous estimates for typical Kepler stars. We find that ≈67% (≈120,000) of all Kepler targets are main-sequence stars, ≈21% (≈37,000) are subgiants, and ≈12% (≈21,000) are red giants, demonstrating that subgiant contamination is less severe than some previous estimates and that Kepler targets are mostly main-sequence stars. Using the revised stellar radii, we recalculate the radii for 2123 confirmed and 1922 candidate exoplanets. We confirm the presence of a gap in the radius distribution of small, close-in planets, but find that the gap is mostly limited to incident fluxes >200 , and its location may be at a slightly larger radius…

Abstract During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral universe, and enables measurements of the hydrogen neutral fraction x H i ( z ) at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Ly α . We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its…

Abstract Detection of the IceCube-170922A neutrino coincident with the flaring blazar TXS 0506+056, the first and only ∼3σ high-energy neutrino source association to date, offers a potential breakthrough in our understanding of high-energy cosmic particles and blazar physics. We present a comprehensive analysis of TXS 0506+056 during its flaring state, using newly collected Swift, NuSTAR, and X-shooter data with Fermi observations and numerical models to constrain the blazar’s particle acceleration processes and multimessenger (electromagnetic (EM) and high-energy neutrino) emissions. Accounting properly for EM cascades in the emission region, we find a physically consistent picture only within a hybrid leptonic scenario, with γ-rays produced by external inverse-Compton processes and high-energy neutrinos via a radiatively subdominant hadronic component. We derive robust constraints on the blazar’s neutrino and cosmic-ray emissions and demonstrate that, because of cascade effects, the 0.1–100 keV emissions of TXS 0506+056 serve as a better probe of its hadronic acceleration and…

Abstract We present a survey of the [C ii ] 158 μ m line and underlying far-infrared (FIR) dust continuum emission in a sample of 27 quasars using the Atacama Large Millimeter Array (ALMA) at resolution. The [C ii ] line was significantly detected (at - σ ) in 23 sources (85%). We find typical line luminosities of , and an average line width of ∼385 . The [C ii ]-to-far-infrared luminosity ratios ([C ii ]/FIR) in our sources span one order of magnitude, highlighting a variety of conditions in the star-forming medium. Four quasar host galaxies are clearly resolved in their [C ii ] emission on a few kpc scales. Basic estimates of the dynamical masses of the host galaxies give masses between 2 × 10 10 and 2 × 10 11 , i.e., more than an order of magnitude below what is expected from local scaling relations, given the available limits on the masses of the central black holes ( , assuming Eddington-limited accretion). In stacked ALMA [C ii ] spectra of individual sources in our sample, we find no evidence of a deviation from a single Gaussian profile. The quasar luminosity does not…

Abstract Black hole (BH) mergers detectable with the Laser Interferometer Gravitational-wave Observatory (LIGO) can occur in active galactic nucleus (AGN) disks. Here we parameterize the merger rates, the mass spectrum, and the spin spectrum of BHs in AGN disks. The predicted merger rate spans ∼10 −3 –10 4 Gpc −1 yr −1 , so upper limits from LIGO (<212 Gpc −1 yr −1 ) already constrain it. The predicted mass spectrum has the form of a broken power law, consisting of a pre-existing BH power-law mass spectrum and a harder power-law mass spectrum resulting from mergers. The predicted spin spectrum is multipeaked with the evolution of retrograde spin BHs in the gas disk playing a key role. We outline the large uncertainties in each of these LIGO observables for this channel and we discuss ways in which they can be constrained in the future.

Abstract We present UV luminosity functions of dropout galaxies at with the complete Hubble Frontier Fields data. We obtain a catalog of ∼450 dropout-galaxy candidates (350, 66, and 40 at , 8, and 9, respectively), with UV absolute magnitudes that reach mag, ∼2 mag deeper than the Hubble Ultra Deep Field detection limits. We carefully evaluate number densities of the dropout galaxies by Monte Carlo simulations, including all lensing effects such as magnification, distortion, and multiplication of images as well as detection completeness and contamination effects in a self-consistent manner. We find that UV luminosity functions at have steep faint-end slopes, , and likely steeper slopes, at . We also find that the evolution of UV luminosity densities shows a non-accelerated decline beyond in the case of , but an accelerated one in the case of . We examine whether our results are consistent with the Thomson scattering optical depth from the Planck satellite and the ionized hydrogen fraction Q H ii at based on the standard analytic reionization model. We find that reionization…

Abstract We perform a search for stellar streams around the Milky Way using the first 3 yr of multiband optical imaging data from the Dark Energy Survey (DES). We use DES data covering ∼5000 deg2 to a depth of g > 23.5 with a relative photometric calibration uncertainty of <1%. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of ∼50 kpc. We search for stellar streams using a matched filter in color–magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of 11 new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extratidal…

Abstract The Asteroid Terrestrial-impact Last Alert System (ATLAS) observes most of the sky every night in search of dangerous asteroids. Its data are also used to search for photometric variability, where sensitivity to variability is limited by photometric accuracy. Since each exposure spans 7.°6 corner to corner, variations in atmospheric transparency in excess of 0.01 mag are common, and 0.01 mag photometry cannot be achieved by using a constant flat-field calibration image. We therefore have assembled an all-sky reference catalog of approximately one billion stars to m ∼ 19 from a variety of sources to calibrate each exposure’s astrometry and photometry. Gaia DR2 is the source of astrometry for this ATLAS Refcat2. The sources of g, r, i, and z photometry include Pan-STARRS DR1, the ATLAS Pathfinder photometry project, ATLAS reflattened APASS data, SkyMapper DR1, APASS DR9, the Tycho-2 catalog, and the Yale Bright Star Catalog. We have attempted to make this catalog at least 99% complete to m < 19, including the brightest stars in the sky. We believe that the systematic errors…

Abstract Nuclear star clusters around a central massive black hole (MBH) are expected to be abundant in stellar black hole (BH) remnants and BH–BH binaries. These binaries form a hierarchical triple system with the central MBH, and gravitational perturbations from the MBH can cause high-eccentricity excitation in the BH–BH binary orbit. During this process, the eccentricity may approach unity, and the pericenter distance may become sufficiently small so that gravitational-wave emission drives the BH–BH binary to merge. In this work, we construct a simple proof-of-concept model for this process, and specifically, we study the eccentric Kozai–Lidov mechanism in unequal-mass, soft BH–BH binaries. Our model is based on a set of Monte Carlo simulations for BH–BH binaries in galactic nuclei, taking into account quadrupole- and octupole-level secular perturbations, general relativistic precession, and gravitational-wave emission. For a typical steady-state number of BH–BH binaries, our model predicts a total merger rate of ∼1–3 −3 −1, depending on the assumed density profile in the…

Abstract We use an end-to-end model of planet formation, thermodynamic evolution, and atmospheric escape to investigate how the statistical imprints of evaporation depend on the bulk composition of planetary cores (rocky versus icy). We find that the population-wide imprints like the location of the “evaporation valley” in the distance–radius plane and the corresponding bimodal radius distribution clearly differ depending on the bulk composition of the cores. Comparison with the observed position of the valley suggests that close-in low-mass Kepler planets have a predominantly Earth-like rocky composition. Combined with the excess of period ratios outside of MMR, this suggests that low-mass Kepler planets formed inside of the water iceline but were still undergoing orbital migration. The core radius becomes visible for planets losing all primordial H/He. For planets in this “triangle of evaporation” in the distance–radius plane, the degeneracy in composition is reduced. In the observed planetary mass–mean density diagram, we identify a trend to more volatile-rich compositions with…

Abstract Extremely irradiated hot Jupiters, exoplanets reaching dayside temperatures >2000 K, stretch our understanding of planetary atmospheres and the models we use to interpret observations. While these objects are planets in every other sense, their atmospheres reach temperatures at low pressures comparable only to stellar atmospheres. In order to understand our a priori theoretical expectations for the nature of these objects, we self-consistently model a number of extreme hot Jupiter scenarios with the PHOENIX model atmosphere code. PHOENIX is well-tested on objects from cool brown dwarfs to expanding supernovae shells, and its expansive opacity database from the UV to far-IR make PHOENIX well-suited to understanding extremely irradiated hot Jupiters. We find several fundamental differences between hot Jupiters at temperatures >2500 K and their cooler counterparts. First, absorption by atomic metals like Fe and Mg, molecules including SiO and metal hydrides, and continuous opacity sources like H−, all combined with the short-wavelength output of early-type host stars, result…

Abstract Gravitational-wave detections have revealed a previously unknown population of stellar mass black holes with masses above 20 M ⊙ . These observations provide a new way to test models of stellar evolution for massive stars. By considering the astrophysical processes likely to determine the shape of the binary black hole mass spectrum, we construct a parameterized model to capture key spectral features that relate gravitational-wave data to theoretical stellar astrophysics. In particular, we model the signature of pulsational pair-instability supernovae, which are expected to cause all stars with initial mass 100 M ⊙ ≲ M ≲ 150 M ⊙ to form ∼40 M ⊙ black holes. This would cause a cutoff in the black hole mass spectrum along with an excess of black holes near 40 M ⊙ . We carry out a simulated data study to illustrate some of the stellar physics that can be inferred using gravitational-wave measurements of binary black holes and demonstrate several such inferences that might be made in the near future. First, we measure the minimum and maximum stellar black hole mass. Second, we…

Abstract The detonation of a sub-Chandrasekhar-mass white dwarf (WD) has emerged as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. Recent studies have suggested that the rapid transfer of a very small amount of helium from one WD to another is sufficient to ignite a helium shell detonation that subsequently triggers a carbon core detonation, yielding a “dynamically driven double-degenerate double-detonation” SN Ia. Because the helium shell that surrounds the core explosion is so minimal, this scenario approaches the limiting case of a bare C/O WD detonation. Motivated by discrepancies in previous literature and by a recent need for detailed nucleosynthetic data, we revisit simulations of naked C/O WD detonations in this paper. We disagree to some extent with the nucleosynthetic results of previous work on sub-Chandrasekhar-mass bare C/O WD detonations; for example, we find that a median-brightness SN Ia is produced by the detonation of a 1.0 WD instead of a more massive and rarer 1.1 WD. The neutron-rich nucleosynthesis in our simulations agrees broadly…

Abstract We constrain the intrinsic architecture of Kepler planetary systems by modeling the observed multiplicities of the transiting planets (tranets) and their transit timing variations (TTVs). We robustly determine that the fraction of Sun-like stars with Kepler-like planets, η Kepler, is 30 ± 3%. Here, Kepler-like planets are planets that have radii R p ≳ R ⊕ and orbital periods P < 400 days. Our result thus significantly revises previous claims that more than 50% of Sun-like stars have such planets. Combined with the average number of Kepler planets per star (∼0.9), we obtain that on average each planetary system has 3.0 ± 0.3 planets within 400 days. We also find that the dispersion in orbital inclinations of planets within a given planetary system, σ i,k , is a steep function of its number of planets, k. This can be parameterized as and we find that −4 < α < −2 at the 2σ level. Such a distribution well describes the observed multiplicities of both transits and TTVs with no excess of single-tranet systems. Therefore we do not find evidence supporting the so-called “Kepler…

Abstract We examine the impact of baryon acoustic oscillation (BAO) scale measurements on the discrepancy between the value of the Hubble constant ( H 0 ) inferred from the local distance ladder and that from Planck cosmic microwave background (CMB) data. While the BAO data alone cannot constrain H 0 , we show that combining the latest BAO results with WMAP , Atacama Cosmology Telescope (ACT), or South Pole Telescope (SPT) CMB data produces values of H 0 that are lower than the distance ladder, independent of Planck , and that this downward pull was less apparent in some earlier analyses that used only angle-averaged BAO scale constraints rather than full anisotropic information. At the same time, the combination of BAO and CMB data also disfavors the lower values of H 0 preferred by the Planck high-multipole temperature power spectrum. Combining galaxy and Ly α forest BAO with a precise estimate of the primordial deuterium abundance produces km s −1 Mpc −1 for the flat model. This value is completely independent of CMB anisotropy constraints and is lower than the latest distance…

Abstract Interferometric imaging now achieves angular resolutions as fine as ∼10 μ as, probing scales that are inaccessible to single telescopes. Traditional synthesis imaging methods require calibrated visibilities; however, interferometric calibration is challenging, especially at high frequencies. Nevertheless, most studies present only a single image of their data after a process of “self-calibration,” an iterative procedure where the initial image and calibration assumptions can significantly influence the final image. We present a method for efficient interferometric imaging directly using only closure amplitudes and closure phases, which are immune to station-based calibration errors. Closure-only imaging provides results that are as noncommittal as possible and allows for reconstructing an image independently from separate amplitude and phase self-calibration. While closure-only imaging eliminates some image information (e.g., the total image flux density and the image centroid), this information can be recovered through a small number of additional constraints. We…

Abstract A large planet orbiting a star in a protoplanetary disk opens a density gap along its orbit due to the strong disk–planet interaction and migrates with the gap in the disk. It is expected that in the ideal case, a gap-opening planet migrates at the viscous drift speed, which is referred to as type II migration. However, recent hydrodynamic simulations have shown that, in general, the gap-opening planet is not locked to the viscous disk evolution. A new physical model is required to explain the migration speed of such a planet. For this reason, we re-examined the migration of a planet in the disk, by carrying out the two-dimensional hydrodynamic simulations in a wide parameter range. We have found that the torque exerted on the gap-opening planet depends on the surface density at the bottom of the gap. The planet migration slows down as the surface density of the bottom of the gap decreases. Using the gap model developed in our previous studies, we have constructed an empirical formula of the migration speed of the gap-opening planets, which is consistent with the results…

Abstract We explore the connection between the UV luminosity functions (LFs) of high-z galaxies and the distribution of stellar masses and star formation histories (SFHs) in their host dark matter halos. We provide a baseline for a redshift-independent star formation efficiency model to which observations and models can be compared. Our model assigns a star formation rate (SFR) to each dark matter halo based on the growth rate of the halo and a redshift-independent star formation efficiency. The dark matter halo accretion rate is obtained from a high-resolution N-body simulation in order to capture the stochasticity in accretion histories and to obtain spatial information for the distribution of galaxies. The halo mass dependence of the star formation efficiency is calibrated at z = 4 by requiring a match to the observed UV LF at this redshift. The model then correctly predicts the observed UV LF at z = 5–10. We present predictions for the UV luminosity and stellar mass functions, JWST number counts, and SFHs. In particular, we find a stellar-to-halo mass relation at z = 4–10 that…

Abstract The density structure surrounding the iron core of a massive star when it dies is known to have a major effect on whether or not the star explodes. Here we repeat previous surveys of presupernova evolution with some important corrections to code physics and four to 10 times better mass resolution in each star. The number of presupernova masses considered is also much larger. Over 4000 models are calculated in the range from 12 to 60 M ⊙ with varying mass loss rates. The core structure is not greatly affected by the increased spatial resolution. The qualitative patterns of compactness measures and their extrema are the same, but with the increased number of models, the scatter seen in previous studies is replaced by several localized branches. More physics-based analyses by Ertl et al. and Müller et al. show these branches with less scatter than the single-parameter characterization of O’Connor & Ott. These branches are particularly apparent for stars in the mass ranges 14–19 and 22–24 M ⊙ . The multivalued solutions are a consequence of interference between several carbon-…

Abstract We perform long-term general relativistic neutrino radiation hydrodynamics simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a torus, which is a canonical remnant formed after the binary neutron star merger. We take into account the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence. The viscous effect plays key roles for the mass ejection from the remnant in two phases of the evolution. In the first t ≲ 10 ms, a differential rotation state of the MNS is changed to a rigidly rotating state. A shock wave caused by the variation of its quasi-equilibrium state induces significant mass ejection of mass ∼(0.5–2.0) for the α-viscosity parameter of 0.01–0.04. For the longer-term evolution with ∼0.1–10 s, a significant fraction of the torus material is ejected. We find that the total mass of the viscosity-driven ejecta ( ) could dominate over that of the dynamical ejecta ( ). The electron fraction, Y e , of the ejecta is always high enough (Y e ≳ 0.25) that this post-merger ejecta is…

Abstract We present an analysis of 15 Type Ia supernovae (SNe Ia) at redshift (9 at ) recently discovered in the CANDELS and CLASH Multi-Cycle Treasury programs using WFC3 on the Hubble Space Telescope. We combine these SNe Ia with a new compilation of ∼1050 SNe Ia, jointly calibrated and corrected for simulated survey biases to produce accurate distance measurements. We present unbiased constraints on the expansion rate at six redshifts in the range based only on this combined SN Ia sample. The added leverage of our new sample at leads to a factor of ∼3 improvement in the determination of the expansion rate at z = 1.5, reducing its uncertainty to ∼20%, a measurement of . We then demonstrate that these six derived expansion rate measurements alone provide a nearly identical characterization of dark energy as the full SN sample, making them an efficient compression of the SN Ia data. The new sample of SNe Ia at usefully distinguishes between alternative cosmological models and unmodeled evolution of the SN Ia distance indicators, placing empirical limits on the latter. Finally,…

Abstract The details of the physical mechanism that drives core-collapse supernovae (CCSNe) remain uncertain. While there is an emerging consensus on the qualitative outcome of detailed CCSN mechanism simulations in 2D, only recently have high-fidelity 3D simulations become possible. Here we present the results of an extensive set of 3D CCSN simulations using high-fidelity multidimensional neutrino transport, high-resolution hydrodynamics, and approximate general relativistic gravity. We employ a state-of-the-art 20 M ☉ progenitor generated using Modules for Experiments in Stellar Astrophysics, and the SFHo equation of state. While none of our 3D CCSN simulations explode within ∼500 ms after core bounce, we find that the presence of large-scale aspherical motion in the Si and O shells aid shock expansion and bring the models closer to the threshold of explosion. We also find some dependence on resolution and geometry (octant versus full 4π). As has been noted in other recent works, we find that the post-shock turbulence plays an important role in determining the overall dynamical…

Abstract Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this “dynamically driven double-degenerate double-detonation” (D 6 ) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its >1000 km s −1 pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia 's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km s −1 , and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is >1000 km s −1 , the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from…

Abstract UV radiation feedback from young massive stars plays a key role in the evolution of giant molecular clouds (GMCs) by photoevaporating and ejecting the surrounding gas. We conduct a suite of radiation hydrodynamic simulations of star cluster formation in marginally bound, turbulent GMCs, focusing on the effects of photoionization and radiation pressure on regulating the net star formation efficiency (SFE) and cloud lifetime. We find that the net SFE depends primarily on the initial gas surface density, Σ0, such that the SFE increases from 4% to 51% as Σ0 increases from 13 to . Cloud destruction occurs within 2–10 Myr after the onset of radiation feedback, or within 0.6–4.1 freefall times (increasing with Σ0). Photoevaporation dominates the mass loss in massive, low surface density clouds, but because most photons are absorbed in an ionization-bounded Strömgren volume, the photoevaporated gas fraction is proportional to the square root of the SFE. The measured momentum injection due to thermal and radiation pressure forces is proportional to , and the ejection of neutrals…

Abstract We report on the design and performance of the Keck Cosmic Web Imager (KCWI), a general purpose optical integral field spectrograph that has been installed at the Nasmyth port of the 10 m Keck II telescope on Maunakea, Hawaii. The novel design provides blue-optimized seeing-limited imaging from 350–560 nm with configurable spectral resolution from 1000–20,000 in a field of view up to 20″ × 33″. Selectable volume phase holographic (VPH) gratings and high-performance dielectric, multilayer silver, and enhanced-aluminum coatings provide end-to-end peak efficiency in excess of 45% while accommodating the future addition of a red channel that will extend wavelength coverage to 1 micron. KCWI takes full advantage of the excellent seeing and dark sky above Maunakea with an available nod-and-shuffle observing mode. The instrument is optimized for observations of faint, diffuse objects such as the intergalactic medium or cosmic web. In this paper, a detailed description of the instrument design is provided with measured performance results from the laboratory test program and 10…

Abstract We measure the velocity dispersion, σ, and surface density, Σ, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution 45–120 pc, matched to the size of individual giant molecular clouds. Combining 11 galaxies from the PHANGS-ALMA survey with four targets from the literature, we characterize ∼30,000 independent sightlines where CO is detected at good significance. Σ and σ show a strong positive correlation, with the best-fit power-law slope close to the expected value for resolved, self-gravitating clouds. This indicates only a weak variation in the virial parameter α vir ∝ σ 2/Σ, which is ∼1.5–3.0 for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure P turb ∝ Σσ 2, which spans ∼5 dex across our sample. We find Σ, σ, and P turb to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kiloparsec of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the σ–Σ relation flattens at Σ ≲ 10 M ⊙ pc−2,…

Abstract Solar flares originate from the release of the energy stored in the magnetic field of solar active regions, the triggering mechanism for these flares, however, remains unknown. For this reason, the conventional solar flare forecast is essentially based on the statistic relationship between solar flares and measures extracted from observational data. In the current work, the deep learning method is applied to set up the solar flare forecasting model, in which forecasting patterns can be learned from line-of-sight magnetograms of solar active regions. In order to obtain a large amount of observational data to train the forecasting model and test its performance, a data set is created from line-of-sight magnetogarms of active regions observed by SOHO /MDI and SDO /HMI from 1996 April to 2015 October and corresponding soft X-ray solar flares observed by GOES . The testing results of the forecasting model indicate that (1) the forecasting patterns can be automatically reached with the MDI data and they can also be applied to the HMI data; furthermore, these forecasting patterns…

Abstract We present a catalog of 10,718 objects in the COSMOS field, observed through multi-slit spectroscopy with the Deep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck II telescope in the wavelength range ∼5500–9800 Å. The catalog contains 6617 objects with high-quality spectra (two or more spectral features), and 1798 objects with a single spectroscopic feature confirmed by the photometric redshift. For 2024 typically faint objects, we could not obtain reliable redshifts. The objects have been selected from a variety of input catalogs based on multi-wavelength observations in the field, and thus have a diverse selection function, which enables the study of the diversity in the galaxy population. The magnitude distribution of our objects is peaked at I AB ∼ 23 and K AB ∼ 21, with a secondary peak at K AB ∼ 24. We sample a broad redshift distribution in the range 0 < z < 6, with one peak at z ∼ 1, and another one around z ∼ 4. We have identified 13 redshift spikes at z > 0.65 with chance probabilities < 4 × 10 −4 , some of which are clearly related to protocluster…

Abstract Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site between core-collapse supernovae and neutron star–neutron star (NS–NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by quantifying the frequency and yields required to reproduce the amount of europium (Eu) observed in galaxies. Although NS–NS mergers have become popular candidates, their required frequency (or rate) needs to be consistent with that obtained from gravitational wave measurements. Here, we address the first NS–NS merger detected by LIGO/Virgo (GW170817) and its associated gamma-ray burst and analyze their implication for the origin of r-process elements. The range of NS–NS merger rate densities of 320–4740 Gpc−3 yr−1 provided by LIGO/Virgo is remarkably consistent with the range required by GCE to explain the Eu abundances in the Milky Way with NS–NS mergers, assuming the solar…

Abstract We use a newly assembled sample of 3545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z = 1.5–2.5 to constrain the relationship between UV slope (β) and dust attenuation (L IR/L UV ≡ IRX). Our sample significantly extends the range of L UV and β probed in previous UV-selected samples, including those as faint as M 1600 = −17.4 ( ) and −2.6 ≲ β ≲ 0.0. IRX is measured using stacks of deep Herschel data, and the results are compared with predictions of the IRX−β relation for different assumptions of the stellar population model and obscuration curve. We find that z = 1.5–2.5 galaxies have an IRX−β relation that is consistent with the predictions for an SMC curve if we invoke subsolar-metallicity models currently favored for high-redshift galaxies, while the commonly assumed starburst curve overpredicts the IRX at a given β by a factor of ≳3. IRX is roughly constant with L UV for L UV ≳ 3 × 109 L ⊙. Thus, the commonly observed trend of fainter galaxies having bluer β may simply reflect bluer intrinsic slopes for such galaxies, rather than…

Abstract We present new measurements of the quasar luminosity function (LF) at z ∼ 6 over an unprecedentedly wide range of the rest-frame ultraviolet luminosity M 1450 from −30 to −22 mag. This is the fifth in a series of publications from the Subaru High- z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multiband imaging data produced by the Hyper Suprime-Cam Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at 5.7 ≤ z ≤ 6.5, which includes 48 SHELLQs quasars discovered over 650 deg 2 and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada–France–Hawaii Quasar Survey (including one overlapping object). This is the largest sample of z ∼ 6 quasars with a well-defined selection function constructed to date, which has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope , a bright-end slope , a break magnitude , and a characteristic space density Gpc −3 mag −1 . Integrating this best-fit model…

Abstract We consider implications of high-energy neutrino emission from blazar flares, including the recent event IceCube-170922A and the 2014–2015 neutrino flare that could originate from TXS 0506+056. First, we discuss their contribution to the diffuse neutrino intensity taking into account various observational constraints. Blazars are likely to be subdominant in the diffuse neutrino intensity at sub-PeV energies, and we show that blazar flares like those of TXS 0506+056 could make ≲1%–10% of the total neutrino intensity. We also argue that the neutrino output of blazars can be dominated by the flares in the standard leptonic scenario for their γ-ray emission, and energetic flares may still be detected with a rate of . Second, we consider multi-messenger constraints on the source modeling. We show that luminous neutrino flares should be accompanied by luminous broadband cascade emission, emerging also in X-rays and γ-rays. This implies that not only γ-ray telescopes like Fermi but also X-ray sky monitors such as Swift and MAXI are critical to test the canonical picture based on…

Abstract We present structural parameters from a wide-field homogeneous imaging survey of Milky Way satellites carried out with the MegaCam imagers on the 3.6 m Canada–France–Hawaii Telescope and 6.5 m Magellan-Clay telescope. Our survey targets an unbiased sample of “outer halo” satellites (i.e., substructures having galactocentric distances greater than 25 kpc) and includes classical dSph galaxies, ultra-faint dwarfs, and remote globular clusters. We combine deep, panoramic gr imaging for 44 satellites and archival gr imaging for 14 additional objects (primarily obtained with the DECam instrument as part of the Dark Energy Survey) to measure photometric and structural parameters for 58 outer halo satellites. This is the largest and most uniform analysis of Milky Way satellites undertaken to date and represents roughly three-quarters (58/81 ≃ 72%) of all known outer halo satellites. We use a maximum-likelihood method to fit four density laws to each object in our survey: exponential, Plummer, King, and Sérsic models. We systematically examine the isodensity contour maps and…

Abstract The rare case of changing-look (CL) active galactic nuclei (AGNs), with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at 0.08 < z < 0.58 , which doubles the number of such objects known to date. These new CL AGNs were discovered in various ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer ( WISE ), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of the transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest-frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of…

Abstract Turbulence is a fundamental parameter in models of grain growth during the early stages of planet formation. As such, observational constraints on its magnitude are crucial. Here we self-consistently analyze ALMA CO(2–1), SMA CO(3–2), and SMA CO(6–5) observations of the disk around TW Hya and find an upper limit on the turbulent broadening of <0.08 c s ( α < 0.007 for α defined only within 2–3 pressure scale heights above the midplane), lower than the tentative detection previously found from an analysis of the CO(2–1) data. We examine in detail the challenges of image plane fitting versus directly fitting the visibilities, while also considering the role of the vertical temperature gradient, systematic uncertainty in the amplitude calibration, and assumptions about the CO abundance, as potential sources of the discrepancy in the turbulence measurements. These tests result in variations of the turbulence limit between <0.04 c s and <0.13 c s , consistently lower than the 0.2–0.4 c s found previously. Having ruled out numerous factors, we restrict the source of the…

Abstract The basic unified model of active galactic nuclei (AGNs) invokes an anisotropic obscuring structure, usually referred to as a torus, to explain AGN obscuration as an angle-dependent effect. We present a new grid of X-ray spectral templates based on radiative transfer calculations in neutral gas in an approximately toroidal geometry, appropriate for CCD-resolution X-ray spectra (FWHM ≥ 130 eV). Fitting the templates to broadband X-ray spectra of AGNs provides constraints on two important geometrical parameters of the gas distribution around the supermassive black hole: the average column density and the covering factor. Compared to the currently available spectral templates, our model is more flexible, and capable of providing constraints on the main torus parameters in a wider range of AGNs. We demonstrate the application of this model using hard X-ray spectra from NuSTAR (3–79 keV) for four AGNs covering a variety of classifications: 3C 390.3, NGC 2110, IC 5063, and NGC 7582. This small set of examples was chosen to illustrate the range of possible torus configurations,…

Abstract Merging binaries consisting of two neutron stars (NSs) or an NS and a stellar-mass black hole typically form a massive accretion torus around the remnant black hole or long-lived NS. Outflows from these neutrino-cooled accretion disks represent an important site for r-process nucleosynthesis and the generation of kilonovae. We present the first three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulations including weak interactions and a realistic equation of state of such accretion disks over viscous timescales (380 ms). We witness the emergence of steady-state MHD turbulence, a magnetic dynamo with an ∼20 ms cycle, and the generation of a “hot” disk corona that launches powerful thermal outflows aided by the energy released as free nucleons recombine into α-particles. We identify a self-regulation mechanism that keeps the midplane electron fraction low (Y e ∼ 0.1) over viscous timescales. This neutron-rich reservoir, in turn, feeds outflows that retain a sufficiently low value of Y e ≈ 0.2 to robustly synthesize third-peak r-process elements. The…

Abstract As one paper in a series reporting on a large reverberation mapping campaign of super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs), we present the results of 10 SEAMBHs monitored spectroscopically during 2015–2017. Six of them are observed for the first time, and have generally higher 5100 Å luminosities than the SEAMBHs monitored in our campaign from 2012 to 2015; the remaining four are repeat observations to check if their previous lags change. Similar to the previous SEAMBHs, the Hβ time lags of the newly observed objects are shorter than the values predicted by the canonical R Hβ –L 5100 relation of sub-Eddington AGNs, by factors of ∼2–6, depending on the accretion rate. The four previously observed objects have lags consistent with previous measurements. We provide linear regressions for the R Hβ –L 5100 relation, solely for the SEAMBH sample and for low-accretion AGNs. We find that the relative strength of Fe ii and the profile of the Hβ emission line can be used as proxies of accretion rate, showing that the shortening of Hβ lags…

Abstract We report the first science results from the newly completed Expanded Owens Valley Solar Array (EOVSA), which obtained excellent microwave (MW) imaging spectroscopy observations of SOL2017-09-10, a classic partially occulted solar limb flare associated with an erupting flux rope. This event is also well-covered by the Reuven Ramaty High Energy Solar Spectroscopic Imager ( RHESSI ) in hard X-rays (HXRs). We present an overview of this event focusing on MW and HXR data, both associated with high-energy nonthermal electrons, and we discuss them within the context of the flare geometry and evolution revealed by extreme ultraviolet observations from the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory . The EOVSA and RHESSI data reveal the evolving spatial and energy distribution of high-energy electrons throughout the entire flaring region. The results suggest that the MW and HXR sources largely arise from a common nonthermal electron population, although the MW imaging spectroscopy provides information over a much larger volume of the corona.

Abstract We estimate the 21 cm radio background from accretion onto the first intermediate-mass black holes between z ≈ 30 and z ≈ 16. Combining potentially optimistic, but plausible, scenarios for black hole formation and growth with empirical correlations between luminosity and radio emission observed in low-redshift active galactic nuclei, we find that a model of black holes forming in molecular cooling halos is able to produce a 21 cm background that exceeds the cosmic microwave background (CMB) at z ≈ 17, though models involving larger halo masses are not entirely excluded. Such a background could explain the surprisingly large amplitude of the 21 cm absorption feature recently reported by the EDGES collaboration. Such black holes would also produce significant X-ray emission and contribute to the 0.5–2 keV soft X-ray background at the level of ≈10−13–10−12 erg s−1 cm−2 deg−2, consistent with existing constraints. In order to avoid heating the intergalactic medium (IGM) over the EDGES trough, these black holes would need to be obscured by hydrogen column depths of N H ∼ 5 ×…

Abstract The neutron star (NS) merger GW170817 was followed over several days by optical-wavelength (“blue”) kilonova (KN) emission likely powered by the radioactive decay of light r-process nuclei synthesized by ejecta with a low neutron abundance (electron fraction Y e ≈ 0.25–0.35). While the composition and high velocities of the blue KN ejecta are consistent with shock-heated dynamical material, the large quantity is in tension with the results of numerical simulations. We propose an alternative ejecta source: the neutrino-heated, magnetically accelerated wind from the strongly magnetized hypermassive NS (HMNS) remnant. A rapidly spinning HMNS with an ordered surface magnetic field of strength B ≈ (1–3) × 1014 G and lifetime t rem ∼ 0.1–1 s can simultaneously explain the velocity, total mass, and electron fraction of the blue KN ejecta. The inferred HMNS lifetime is close to its Alfvén crossing time, suggesting that global magnetic torques could be responsible for bringing the HMNS into solid-body rotation and instigating its gravitational collapse. Different origins for the KN…

Abstract Polaris, the nearest and brightest Cepheid, is a potential anchor point for the Leavitt period–luminosity relation. However, its distance is a matter of contention, with recent advocacy for a parallax of ∼10 mas, in contrast with the Hipparcos measurement of 7.54 ± 0.11 mas. We report an independent trigonometric parallax determination, using the Fine Guidance Sensors (FGS) on the Hubble Space Telescope. Polaris itself is too bright for FGS, so we measured its eighth-magnitude companion Polaris B, relative to a network of background reference stars. We converted the FGS relative parallax to absolute, using estimated distances to the reference stars from ground-based photometry and spectral classification. Our result, 6.26 ± 0.24 mas, is even smaller than that found by Hipparcos. We note other objects for which Hipparcos appears to have overestimated parallaxes, including the well-established case of the Pleiades. We consider possible sources of systematic error in the FGS parallax, but find no evidence they are significant. If our “long” distance is correct, the high…

Abstract Current and upcoming radio telescopes will map the spatial distribution of cosmic neutral hydrogen (H i) through its 21 cm emission. In order to extract the maximum information from these surveys, accurate theoretical predictions are needed. We study the abundance and clustering properties of H i at redshifts z ≤ 5 using TNG100, a large state-of-the-art magnetohydrodynamic simulation of a 75 h −1 Mpc box size, which is part of the IllustrisTNG Project. We show that most of the H i lies within dark matter halos, and we provide fits for the halo H i mass function, i.e., the mean H i mass hosted by a halo of mass M at redshift z. We find that only halos with circular velocities larger than ≃30 km s−1 contain H i. While the density profiles of H i exhibit a large halo-to-halo scatter, the mean profiles are universal across mass and redshift. The H i in low-mass halos is mostly located in the central galaxy, while in massive halos the H i is concentrated in the satellites. Our simulation reproduces the bias value of damped Lyα systems from observations. We show that the H i and…

Abstract We report the identification of an extreme protocluster of galaxies in the early universe whose core (nicknamed Distant Red Core, DRC, because of its very red color in Herschel SPIRE bands) is formed by at least 10 dusty star-forming galaxies (DSFGs), spectroscopically confirmed to lie at z spec = 4.002 via detection of [C i ](1–0), 12 CO(6–5), 12 CO(4–3), 12 CO(2–1), and H 2 O ( 2 11 – 2 02 ) emission lines with ALMA and ATCA. These DSFGs are distributed over a 260 kpc × 310 kpc region and have a collective obscured star formation rate (SFR) of ∼ 6500 M ⊙ yr − 1 , considerably higher than those seen before in any protocluster at z ≳ 4 . Most of the star formation is taking place in luminous DSFGs since no Ly α emitters are detected in the protocluster core, apart from a Ly α blob located next to one of the DRC components, extending over 60 kpc . The total obscured SFR of the protocluster could rise to SFR ∼ 14,400 M ⊙ yr − 1 if all the members of an overdensity of bright DSFGs discovered around DRC in a wide-field Large APEX BOlometer CAmera 870 μ m image are part of the…

Abstract Gas blown away from galactic disks by supernova (SN) feedback plays a key role in galaxy evolution. We investigate outflows utilizing the solar neighborhood model of our high-resolution, local galactic disk simulation suite, TIGRESS. In our numerical implementation, star formation and SN feedback are self-consistently treated and well resolved in the multiphase, turbulent, magnetized interstellar medium. Bursts of star formation produce spatially and temporally correlated SNe that drive strong outflows, consisting of hot ( ) winds and warm ( ) fountains. The hot gas at distance from the midplane has mass and energy fluxes nearly constant with d. The hot flow escapes our local Cartesian box barely affected by gravity, and is expected to accelerate up to terminal velocity of . The mean mass and energy loading factors of the hot wind are 0.1 and 0.02, respectively. For warm gas, the mean outward mass flux through is comparable to the mean star formation rate, but only a small fraction of this gas is at velocity . Thus, the warm outflows eventually fall back as inflows. The…

Abstract We present measurements of the E-mode polarization angular auto-power spectrum (EE) and temperature–E-mode cross-power spectrum (TE) of the cosmic microwave background (CMB) using 150 GHz data from three seasons of SPTpol observations. We report the power spectra over the spherical harmonic multipole range and detect nine acoustic peaks in the EE spectrum with high signal-to-noise ratio. These measurements are the most sensitive to date of the EE and TE power spectra at and , respectively. The observations cover 500 , a fivefold increase in area compared to previous SPTpol analyses, which increases our sensitivity to the photon diffusion damping tail of the CMB power spectra enabling tighter constraints on ΛCDM model extensions. After masking all sources with unpolarized flux mJy, we place a 95% confidence upper limit on residual polarized point-source power of at , suggesting that the EE damping tail dominates foregrounds to at least with modest source masking. We find that the SPTpol data set is in mild tension with the ΛCDM model ( ), and different data splits prefer…

Abstract We report the detection of 72 new pulses from the repeating fast radio burst FRB 121102 in Breakthrough Listen C -band (4–8 GHz) observations at the Green Bank Telescope. The new pulses were found with a convolutional neural network in data taken on 2017 August 26, where 21 bursts have been previously detected. Our technique combines neural network detection with dedispersion verification. For the current application, we demonstrate its advantage over a traditional brute-force dedispersion algorithm in terms of higher sensitivity, lower false-positive rates, and faster computational speed. Together with the 21 previously reported pulses, this observation marks the highest number of FRB 121102 pulses from a single observation, totaling 93 pulses in five hours, including 45 pulses within the first 30 minutes. The number of data points reveals trends in pulse fluence, pulse detection rate, and pulse frequency structure. We introduce a new periodicity search technique, based on the Rayleigh test, to analyze the time of arrivals (TOAs), with which we exclude with 99% confidence…

Abstract We present a catalog of extended low surface brightness galaxies (LSBGs) identified in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first ∼200 deg 2 of the survey, we have uncovered 781 LSBGs, spanning red ( g − i ≥ 0.64) and blue ( g − i < 0.64) colors and a wide range of morphologies. Since we focus on extended galaxies ( r eff = 2.″5–14″), our sample is likely dominated by low-redshift objects. We define LSBGs to have mean surface brightnesses mag arcsec −2 , which allows nucleated galaxies into our sample. As a result, the central surface brightness distribution spans a wide range of μ 0 ( g ) = 18–27.4 mag arcsec −2 , with 50% and 95% of galaxies fainter than 24.3 and 22 mag arcsec −2 , respectively. Furthermore, the surface brightness distribution is a strong function of color, with the red distribution being much broader and generally fainter than that of the blue LSBGs, and this trend shows a clear correlation with galaxy morphology. Red LSBGs typically have smooth light profiles that are well characterized by…

Abstract The Wide Field InfraRed Survey Telescope ( WFIRST ) was the highest-ranked large space-based mission of the 2010 New Worlds , New Horizons decadal survey. It is now a NASA mission in formulation with a planned launch in the mid 2020s. A primary mission objective is to precisely constrain the nature of dark energy through multiple probes, including Type Ia supernovae (SN Ia). Here, we present the first realistic simulations of the WFIRST SN survey based on current hardware specifications and using open-source tools. We simulate SN light curves and spectra as viewed by the WFIRST wide-field channel (WFC) imager and integral field channel (IFC) spectrometer, respectively. We examine 11 survey strategies with different time allocations between the WFC and IFC, two of which are based upon the strategy described by the WFIRST Science Definition Team, which measures SN distances exclusively from IFC data. We propagate statistical and, crucially, systematic uncertainties to predict the Dark Energy Task Force figure of merit (FoM) for each strategy. Of the strategies investigated,…

Abstract Recent modeling of hydrogen-rich Type II supernova (SN II) light curves suggests the presence of dense circumstellar material (CSM) surrounding the exploding progenitor stars. This has important implications for the activity and structure of massive stars near the end of their lives. Since previous work focused on just a few events, here we expand to a larger sample of 20 well-observed SNe II. For each event we are able to constrain the progenitor zero-age main-sequence (ZAMS) mass, explosion energy, and the mass and radial extent of the dense CSM. We then study the distribution of each of these properties across the full sample of SNe. The inferred ZAMS masses are found to be largely consistent with a Salpeter distribution with minimum and maximum masses of 10.4 and 22.9 M ⊙ , respectively. We also compare the individual ZAMS masses we measure with specific SNe II that have pre-explosion imaging to check their consistency. Our masses are generally comparable to or higher than the pre-explosion imaging masses, potentially helping ease the red supergiant problem. The…

Abstract We investigate r-process nucleosynthesis in 3D general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. The simulations include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission and absorption. We track the composition of the ejected material using the nuclear reaction network SkyNet. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae (CCSN) are important for their nucleosynthetic signature. We find that production of r-process material beyond the second peak is reduced by a factor of 100 when the magnetorotational jets produced by the rapidly rotating core undergo a kink instability. Our results indicate that 3D magnetorotationally powered CCSNe are robust r-process sources only if they are obtained by the collapse of cores with unrealistically large precollapse magnetic fields of the order of 1013 G. Additionally, a comparison simulation that we restrict to…

Abstract We present an analysis of the final data release of the Carnegie Supernova Project I, focusing on the absolute calibration of the luminosity–decline rate relation for Type Ia supernovae (SNe Ia) using new intrinsic color relations with respect to the color-stretch parameter, s BV , enabling improved dust extinction corrections. We investigate to what degree the so-called fast-declining SNe Ia can be used to determine accurate extragalactic distances. We estimate the intrinsic scatter in the luminosity–decline rate relation and find it ranges from ±0.13 mag to ±0.18 mag with no obvious dependence on wavelength. Using the Cepheid variable star data from the SH0ES project, the SN Ia distance scale is calibrated and the Hubble constant is estimated using our optical and near-infrared sample, and these results are compared to those determined exclusively from a near-infrared subsample. The systematic effect of the supernova’s host galaxy mass is investigated as a function of wavelength and is found to decrease toward redder wavelengths, suggesting this effect may be due to dust…

Abstract We present a “super-deblended” far-infrared (FIR) to (sub)millimeter photometric catalog in the Cosmic Evolution Survey (COSMOS), prepared with the method recently developed by Liu et al., with key adaptations. We obtain point-spread function fitting photometry at fixed prior positions including 88,008 galaxies detected in VLA 1.4, 3 GHz, and/or MIPS 24 μ m images. By adding a specifically carved mass-selected sample (with an evolving stellar mass limit), a highly complete prior sample of 194,428 galaxies is achieved for deblending FIR/(sub)mm images. We performed “active” removal of nonrelevant priors at FIR/(sub)mm bands using spectral energy distribution fitting and redshift information. In order to cope with the shallower COSMOS data, we subtract from the maps the flux of faint nonfitted priors and explicitly account for the uncertainty of this step. The resulting photometry (including data from Spitzer , Herschel , SCUBA2, AzTEC, MAMBO, and NSF’s Karl G. Jansky Very Large Array at 3 and 1.4 GHz) displays well-behaved quasi-Gaussian uncertainties calibrated from Monte…

Abstract The extremely high brightness temperatures of pulsars and fast radio bursts (FRBs) require their radiation mechanisms to be coherent. Coherent curvature radiation from bunches has been long discussed as the mechanism for radio pulsars and recently for FRBs. Assuming that bunches are already generated in pulsar magnetospheres, we calculate the spectrum of coherent curvature radiation under a three-dimensional magnetic field geometry. Different from previous works assuming parallel trajectories and a monoenergetic energy distribution of electrons, we consider a bunch characterized by its length, curvature radius of the trajectory family, bunch opening angle, and electron energy distribution. We find that the curvature radiation spectra of the bunches are characterized by a multisegment broken power law, with the break frequencies depending on bunch properties and trajectory configuration. We also emphasize that in a pulsar magnetosphere, only the fluctuation of net charges with respect to the background (Goldreich–Julian) outflow can make a contribution to coherent radiation.…

Abstract We study gravitational waves (GWs) from a set of 2D multigroup neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Our goal is to systematize the current knowledge about the post-bounce CCSN GW signal and recognize the templatable features that could be used by the ground-based laser interferometers. We demonstrate that, starting from ∼400 ms after core bounce, the dominant GW signal represents the fundamental quadrupole ( l = 2) oscillation mode ( f -mode) of the proto–neutron star (PNS), which can be accurately reproduced by a linear perturbation analysis of the angle-averaged PNS profile. Before that, in the time interval between ∼200 and ∼400 ms after bounce, the dominant mode has two radial nodes and represents a g -mode. We associate the high-frequency noise in the GW spectrograms above the main signal with p -modes, while below the dominant frequency there is a region with very little power. The collection of models presented here summarizes the dependence of the CCSN GW signal on the progenitor mass, equation of state, many-body…

Abstract The recent detection of a binary neutron star merger and the clear evidence of the decay of radioactive material observed in this event have, after 60 years of effort, provided an astrophysical site for the rapid neutron-capture (r-) process which is responsible for the production of the heaviest elements in our universe. However, observations of metal-poor stars with highly enhanced r-process elements have revealed abundance patterns suggesting that multiple sites may be involved. To address this issue, and to advance our understanding of the r-process, we have initiated an extensive search for bright (V < 13.5), very metal-poor ([Fe/H] < −2) stars in the Milky Way halo exhibiting strongly enhanced r-process signatures. This paper presents the first sample collected in the southern hemisphere using the echelle spectrograph on du Pont 2.5 m telescope at Las Campanas Observatory. We have observed and analyzed 107 stars with −3.13 < [Fe/H] < −0.79. Of those, 12 stars are strongly enhanced in heavy r-process elements (r-II), 42 stars show moderate enhancements of heavy…

Abstract Solar eruptions are the main driver of space-weather disturbances at Earth. Extreme events are of particular interest, not only because of the scientific challenges they pose, but also because of their possible societal consequences. Here we present a magnetohydrodynamic (MHD) simulation of the 2000 July 14 “Bastille Day” eruption, which produced a very strong geomagnetic storm. After constructing a “thermodynamic” MHD model of the corona and solar wind, we insert a magnetically stable flux rope along the polarity inversion line of the eruption’s source region and initiate the eruption by boundary flows. More than 1033 erg of magnetic energy is released in the eruption within a few minutes, driving a flare, an extreme-ultraviolet wave, and a coronal mass ejection (CME) that travels in the outer corona at ≈1500 km s−1, close to the observed speed. We then propagate the CME to Earth, using a heliospheric MHD code. Our simulation thus provides the opportunity to test how well in situ observations of extreme events are matched if the eruption is initiated from a stable magnetic…

Abstract We developed a simple, physical, and self-consistent cloud model for brown dwarfs and young giant exoplanets. We compared different parametrizations for the cloud particle size, by fixing either particle radii or the mixing efficiency (parameter f sed ), or by estimating particle radii from simple microphysics. The cloud scheme with simple microphysics appears to be the best parametrization by successfully reproducing the observed photometry and spectra of brown dwarfs and young giant exoplanets. In particular, it reproduces the L–T transition, due to the condensation of silicate and iron clouds below the visible/near-IR photosphere. It also reproduces the reddening observed for low-gravity objects, due to an increase of cloud optical depth for low gravity. In addition, we found that the cloud greenhouse effect shifts chemical equilibrium, increasing the abundances of species stable at high temperature. This effect should significantly contribute to the strong variation of methane abundance at the L–T transition and to the methane depletion observed on young exoplanets.…

Abstract We examine nucleosynthesis in the innermost neutrino-processed ejecta (a few ) of self-consistent two-dimensional explosion models of core-collapse supernovae (CCSNe) for six progenitor stars with different initial masses. Three models have initial masses near the low-mass end of the SN range of (e8.8; electron-capture SN), (z9.6), and (u8.1), with initial metallicities of 1, 0, and 10−4 times the solar metallicity, respectively. The other three are solar-metallicity models with initial masses of (s11), (s15), and (s27). The low-mass models e8.8, z9.6, and u8.1 exhibit high production factors (nucleosynthetic abundances relative to the solar abundances) of 100–200 for light trans-Fe elements from Zn to Zr. This is associated with an appreciable ejection of neutron-rich matter in these models. Remarkably, the nucleosynthetic outcomes for the progenitors e8.8 and z9.6 are almost identical, including interesting productions of 48Ca and 60Fe, irrespective of their quite different (O–Ne–Mg and Fe) cores prior to collapse. In the more massive models s11, s15, and s27, several…

Abstract The compositions of stars are a critical diagnostic tool for many topics in astronomy such as the evolution of our Galaxy, the formation of planets, and the uniqueness of the Sun. Previous spectroscopic measurements indicate a large intrinsic variation in the elemental abundance patterns of stars with similar overall metal content. However, systematic errors arising from inaccuracies in stellar models are known to be a limiting factor in such studies, and thus it is uncertain to what extent the observed diversity of stellar abundance patterns is real. Here we report the abundances of 30 elements with precisions of 2% for 79 Sun-like stars within 100 pc. Systematic errors are minimized in this study by focusing on solar twin stars and performing a line-by-line differential analysis using high-resolution, high-signal-to-noise spectra. We resolve [X/Fe] abundance trends in galactic chemical evolution at precisions of 10 −3 dex Gyr −1 and reveal that stars with similar ages and metallicities have nearly identical abundance patterns. Contrary to previous results, we find that…

Abstract The Pan-Andromeda Archaeological Survey is a survey of >400 square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of an L ⋆ galaxy group to large projected galactocentric radii. Here, we collate and summarize the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyr. We suggest that a few of the substructures farthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured data set; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo and identify three new faint structures. M31's halo, in…

Abstract We analyze high angular resolution ALMA observations of the TW Hya disk to place constraints on the CO and dust properties. We present new, sensitive observations of the 12CO J = 3 − 2 line at a spatial resolution of 8 au (0.″14). The CO emission exhibits a bright inner core, a shoulder at r ≈ 70 au, and a prominent break in slope at r ≈ 90 au. Radiative transfer modeling is used to demonstrate that the emission morphology can be reasonably reproduced with a 12CO column density profile featuring a steep decrease at r ≈ 15 au and a secondary bump peaking at r ≈ 70 au. Similar features have been identified in observations of rarer CO isotopologues, which trace heights closer to the midplane. Substructure in the underlying gas distribution or radially varying CO depletion that affects much of the disk’s vertical extent may explain the shared emission features of the main CO isotopologues. We also combine archival 1.3 mm and 870 μm continuum observations to produce a spectral index map at a spatial resolution of 2 au. The spectral index rises sharply at the continuum emission…

Abstract We developed a solar flare prediction model using a deep neural network (DNN) named Deep Flare Net (DeFN). This model can calculate the probability of flares occurring in the following 24 hr in each active region, which is used to determine the most likely maximum classes of flares via a binary classification (e.g., ≥M class versus 5 observation images taken during 2010–2015 by the Solar Dynamic Observatory , we automatically detected sunspots and calculated 79 features for each region, to which flare occurrence labels of X-, M-, and C-class were attached. We adopted the features used in Nishizuka et al. (2017) and added some features for operational prediction: coronal hot brightening at 131 Å ( T ≥ 10 7 K) and the X-ray and 131 Å intensity data 1 and 2 hr before an image. For operational evaluation, we divided the database into two for training and testing: the data set in 2010–2014 for training, and the one in 2015 for testing. The DeFN model consists of deep multilayer neural networks formed by adapting skip connections and batch normalizations. To statistically predict…

Abstract We report a new changing-look quasar, WISE J105203.55+151929.5 at z = 0.303, found by identifying highly mid-IR-variable quasars in the Wide-field Infrared Survey Explorer ( WISE )/Near-Earth Object WISE Reactivation (NEOWISE) data stream. Compared to multiepoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at 3.4 and 4.6 μ m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of 15 decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 that, when compared with the 2006 archival SDSS spectrum, reveals that the broad H β emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes…

Abstract Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 mm wavelength with 43 × 39 mas angular resolution (6.9 × 6.2 au) and 20 μJy beam−1 rms. The central submillimeter emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to ∼4× the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed…

Abstract We present updated stellar population models appropriate for old ages (>1 Gyr) and covering a wide range in metallicities (−1.5 ≲ [Fe/H] ≲ 0.3). These models predict the full spectral variation associated with individual element abundance variation as a function of metallicity and age. The models span the optical–NIR wavelength range (0.37–2.4 μ m), include a range of initial mass functions, and contain the flexibility to vary 18 individual elements including C, N, O, Mg, Si, Ca, Ti, and Fe. To test the fidelity of the models, we fit them to integrated light optical spectra of 41 Galactic globular clusters (GCs). The value of testing models against GCs is that their ages, metallicities, and detailed abundance patterns have been derived from the Hertzsprung–Russell diagram in combination with high-resolution spectroscopy of individual stars. We determine stellar population parameters from fits to all wavelengths simultaneously (“full spectrum fitting”), and demonstrate explicitly with mock tests that this approach produces smaller uncertainties at fixed signal-to-noise ratio…

Abstract We develop and apply a model to quantify the global efficiency of radial orbit migration among stars in the Milky Way disk. This model parameterizes the possible star formation and enrichment histories and radial birth profiles, and combines them with a migration model that relates present-day orbital radii to birth radii through a Gaussian probability, broadening with age τ as . Guided by observations, we assume that stars are born with an initially tight age–metallicity relation at given radius, which becomes subsequently scrambled by radial orbit migration, thereby providing a direct observational constraint on radial orbit migration strength . We fit this model with Markov Chain Monte Carlo sampling of the observed age–metallicity distribution of low- α red clump stars with Galactocentric radii between 5 and 14 kpc from APOGEE DR12, sidestepping the complex spatial selection function and accounting for the considerable age uncertainties. This simple model reproduces the observed data well, and we find a global (in radius and time) radial orbit migration efficiency in…

Abstract Nearly every massive galaxy harbors a supermassive black hole (SMBH) in its nucleus. SMBH masses are millions to billions of solar mass, and they correlate with properties of spheroids of their host galaxies. While the SMBH growth channels, mergers, and gas accretion are well established, their origin remains uncertain: they could have emerged either from massive “seeds” (105–106 M ⊙) formed by direct collapse of gas clouds in the early universe or from smaller (100 M ⊙) BHs, end products of first stars. The latter channel would leave behind numerous intermediate-mass BHs (IMBHs, 102–105 M ⊙). Although many IMBH candidates have been identified, none are accepted as definitive; thus, their very existence is still debated. Using data mining in wide-field sky surveys and applying dedicated analysis to archival and follow-up optical spectra, we identified a sample of 305 IMBH candidates having masses 3 × 10 4 M ⊙ < M BH < 2 × 10 5 M ⊙ , which reside in galaxy centers and are accreting gas that creates characteristic signatures of a type I active galactic nucleus (AGN). We…

Abstract The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations, as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of ∼100 gravitational radii ( r g ) for various black hole (BH) spins. As is known, the funnel edge is approximately determined by the following equipartitions: (i) the magnetic and rest-mass energy densities and (ii) the gas and magnetic pressures. Our numerical results give an additional factor that they follow the outermost parabolic streamline of the FFE solution, which is anchored to the event horizon on the equatorial plane. We also show that the matter-dominated, nonrelativistic corona/wind plays a dynamical role in shaping the funnel jet into the parabolic geometry. We confirm a quantitative overlap between the outermost parabolic streamline of the FFE jet and the edge of the jet sheath in very long baseline interferometry (VLBI) observations at ∼(10 1 –10 5 ) r g , suggesting that the M87 jet is likely…

Abstract Within the parameter space of the equation of state (EOS) of dense neutron-rich matter limited by existing constraints mainly from terrestrial nuclear experiments, we investigate how the neutron star maximum mass M max > 2.01 ± 0.04 M ⊙ , radius 10.62 km < R 1.4 < 12.83 km and tidal deformability Λ 1.4 ≤ 800 of canonical neutron stars together constrain the EOS of dense neutron-rich nucleonic matter. While the 3D parameter space of K sym (curvature of nuclear symmetry energy), J sym , and J 0 (skewness of the symmetry energy and EOS of symmetric nuclear matter, respectively) is narrowed down significantly by the observational constraints, more data are needed to pin down the individual values of K sym , J sym , and J 0 . The J 0 largely controls the maximum mass of neutron stars. While the EOS with J 0 = 0 is sufficiently stiff to support neutron stars as massive as 2.37 M ⊙ , supporting the hypothetical ones as massive as 2.74 M ⊙ (composite mass of GW170817) requires J 0 to be larger than its currently known maximum value of about 400 MeV and beyond the causality limit.…

Abstract We perform global particle-in-cell simulations of pulsar magnetospheres, including pair production, ion extraction from the surface, frame-dragging corrections, and high-energy photon emission and propagation. In the case of oblique rotators, the effects of general relativity increase the fraction of the open field lines that support active pair discharge. We find that the plasma density and particle energy flux in the pulsar wind are highly non-uniform with latitude. A significant fraction of the outgoing particle energy flux is carried by energetic ions, which are extracted from the stellar surface. Their energies may extend up to a large fraction of the open field line voltage, making them interesting candidates for ultra-high-energy cosmic rays. We show that pulsar gamma-ray radiation is dominated by synchrotron emission, produced by particles that are energized by relativistic magnetic reconnection close to the Y-point and in the equatorial current sheet. In most cases, the calculated light curves contain two strong peaks, which is in general agreement with Fermi…

Abstract We construct z ∼ 6–7, 8, and 9 faint Lyman break galaxy samples (334, 61, and 37 galaxies, respectively) with accurate size measurements with the software glafic from the complete Hubble Frontier Fields (HFF) cluster and parallel fields data. These are the largest samples hitherto and reach down to the faint ends of recently obtained deep luminosity functions. At faint magnitudes, however, these samples are highly incomplete for galaxies with large sizes, implying that derivation of the luminosity function sensitively depends on the intrinsic size–luminosity relation. We thus conduct simultaneous maximum-likelihood estimation of luminosity function and size–luminosity relation parameters from the observed distribution of galaxies on the size–luminosity plane with the help of a completeness map as a function of size and luminosity. At z ∼ 6–7, we find that the intrinsic size–luminosity relation expressed as r e ∝ L β has a notably steeper slope of than those at lower redshifts, which in turn implies that the luminosity function has a relatively shallow faint-end slope of .…

Abstract The TRAPPIST-1 planetary system provides an unprecedented opportunity to study terrestrial exoplanet evolution with the James Webb Space Telescope (JWST) and ground-based observatories. Since M dwarf planets likely experience extreme volatile loss, the TRAPPIST-1 planets may have highly evolved, possibly uninhabitable atmospheres. We used a versatile, 1D terrestrial planet climate model with line-by-line radiative transfer and mixing length convection (VPL Climate) coupled to a terrestrial photochemistry model to simulate environmental states for the TRAPPIST-1 planets. We present equilibrium climates with self-consistent atmospheric compositions and observational discriminants of postrunaway, desiccated, 10–100 bar O2- and CO2-dominated atmospheres, including interior outgassing, as well as for water-rich compositions. Our simulations show a range of surface temperatures, most of which are not habitable, although an aqua planet TRAPPIST-1 e could maintain a temperate surface given Earth-like geological outgassing and CO2. We find that a desiccated TRAPPIST-1 h may produce…