Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Findings”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

The Role of Outflows, Radiation Pressure, and Magnetic Fields in Massive Star Formation

Stellar feedback in the form of radiation pressure and magnetically driven collimated outflows may limit the maximum mass that a star can achieve and affect the star formation efficiency of massive prestellar cores. Here we present a series of 3D adaptive mesh refinement radiation–magnetohydrodynamic simulations of the collapse of initially turbulent, massive prestellar cores. Our simulations include radiative feedback from both the direct stellar and dust-reprocessed radiation fields, and collimated outflow feedback from the accreting stars. We find that protostellar outflows punch holes in the dusty circumstellar gas along the star’s polar directions, thereby increasing the size of optically thin regions through which radiation can escape. Precession of the outflows as the star’s spin axis changes due to the turbulent accretion flow further broadens the outflow, and causes more material to be entrained. Additionally, the presence of magnetic fields in the entrained material leads to broader entrained outflows that escape the core. We compare the injected and entrained outflow properties and find that the entrained outflow mass is a factor of ∼3 larger than the injected mass and the momentum and energy contained in the entrained material are ∼25% and ∼5% of the injected momentum and energy, respectively. As a result, we find that, when one includes both outflows and radiation pressure, the former are a much more effective and important feedback mechanism, even for massive stars with significant radiative outputs.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectral Variability of VHS J1256–1257b from 1 to 5 μm

Multiwavelength time-resolved observations of rotationally modulated variability from brown dwarfs and giant exoplanets are the most effective method for constraining their heterogeneous atmospheric structures. In a companion paper, we reported the discovery of strong near-infrared variability in HST/WFC3/G141 light curves of the very red L-dwarf companion VHS J1256–1257b. In this paper, we present a follow-up 36 hr Spitzer/IRAC Channel 2 light curve together with an in-depth analysis of the Hubble space telescope (HST) and the Spitzer data. The combined data set provides time-resolved light curves of VHS1256b sampling 1.1 to 4.5 μm. The Spitzer light curve is best fit with a single sine wave with a period of 22.04 ± 0.05 hr and a peak-to-peak amplitude of 5.76 ± 0.04%. Combining the period with a previously measured projected rotational velocity (v sin i), we find that VHS1256b is most consistent with equatorial viewing geometry. The HST/G141+Spitzer spectral energy distribution favors a model with a T{sub eff} of 1000 K and low surface gravity with disequilibrium chemistry. The spectral variability of VHS1256b is consistent with predictions from partly cloudy models, suggesting that heterogeneous clouds are the dominant source of the observed modulations. We find evidence at the 3.3σ level for amplitude variations within the 1.67 μm CH{sub 4} band, which is the first such detection for a variable L-dwarf. We compare the HST/G141 time-resolved spectra of three red L-dwarfs with high-amplitude near-infrared rotational modulations and find that although their time-averaged spectra are similar, their spectroscopic variabilities exhibit notable differences. This diversity reinforces the advantage of time-resolved spectroscopic observations for understanding the atmospheres of brown dwarfs and directly imaged exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗

The Gaia–Kepler Stellar Properties Catalog. II. Planet Radius Demographics as a Function of Stellar Mass and Age

Studies of exoplanet demographics require large samples and precise constraints on exoplanet host stars. Using the homogeneous Kepler stellar properties derived using the Gaia Data Release 2 by Berger et al., we recompute Kepler planet radii and incident fluxes and investigate their distributions with stellar mass and age. We measure the stellar mass dependence of the planet radius valley to be dlogR{sub p}/dlogM{sub ⋆} = 0.26{sub −0.16}{sup +0.21}, consistent with the slope predicted by a planet mass dependence on stellar mass (0.24–0.35) and core-powered mass loss (0.33). We also find the first evidence of a stellar age dependence of the planet populations straddling the radius valley. Specifically, we determine that the fraction of super-Earths (1–1.8 R{sub ⊕}) to sub-Neptunes (1.8–3.5 R{sub ⊕}) increases from 0.61 ± 0.09 at young ages (<1 Gyr) to 1.00 ± 0.10 at old ages (>1 Gyr), consistent with the prediction by core-powered mass loss that the mechanism shaping the radius valley operates over Gyr timescales. Additionally, we find a tentative decrease in the radii of relatively cool (F {sub p} < 150 F{sub ⊕}) sub-Neptunes over Gyr timescales, which suggests that these planets may possess H/He envelopes instead of higher mean molecular weight atmospheres. We confirm the existence of planets within the hot sub-Neptunian “desert” (2.2 R {sub ⊕} < R {sub p} < 3.8 R{sub ⊕}, F {sub p} > 650 F{sub ⊕}) and show that these planets are preferentially orbiting more evolved stars compared to other planets at similar incident fluxes. In addition, we identify candidates for cool (F {sub p} < 20 F{sub ⊕}) inflated Jupiters, present a revised list of habitable zone candidates, and find that the ages of single and multiple transiting planet systems are statistically indistinguishable.

79 ASTRONOMY AND ASTROPHYSICS↗

Untangling the Galaxy. II. Structure within 3 kpc

We present the results of the hierarchical clustering analysis of the Gaia DR2 data to search for clusters, comoving groups, and other stellar structures. The current paper builds on the sample from the previous work, extending it in distance from 1 to 3 kpc and increasing the number of identified structures up to 8292. To aid in the analysis of the population properties, we developed a neural network called Auriga to robustly estimate the age, extinction, and distance of a stellar group based on the input photometry and parallaxes of the individual members. We apply Auriga to derive the properties of not only the structures found in this paper, but also previously identified open clusters. Through this work, we examine the temporal structure of the spiral arms. Specifically, we find that the Sagittarius Arm has moved by >500 pc in the last 100 Myr and the Perseus Arm has been experiencing a relative lull in star formation activity over the last 25 Myr. We confirm the findings of the previous paper on the transient nature of the spiral arms, with the timescale of transition of a few 100 Myr. Finally, we find a peculiar ~1 Gyr old stream of stars that appears to be heliocentric. Its origin is unclear.

79 ASTRONOMY AND ASTROPHYSICS↗

A Collage of Small Planets from the Lick–Carnegie Exoplanet Survey: Exploring the Super-Earth and Sub-Neptune Mass Regime

Analysis of new precision radial velocity (RV) measurements from the Lick Automated Planet Finder and Keck HIRES has yielded the discovery of three new exoplanet candidates orbiting the nearby stars HD 190007 and HD 216520. We also report new velocities from the APF and the Planet Finder Spectrograph and updated orbital fits for the known exoplanet host stars GJ 686 and HD 180617. Of the newly discovered planets, HD 190007 b has a period of P = 11.72 days, an RV semiamplitude of K = 5.64 ± 0.55 m s{sup −1}, a minimum mass of M {sub pl} = 16.46 ± 1.66 M {sub ⊕}, and orbits the slightly metal-rich, active K4V star HD 190007. For HD 216520 b, we find P = 35.45 days, K = 2.28 ± 0.20 m s{sup −1}, and M {sub pl} = 10.26 ± 0.99 M {sub ⊕}, while for HD 216520 c, P = 154.43 days, K = 1.29 ± 0.22 m s{sup −1}, and M {sub pl} = 9.44 ± 1.63 M {sub ⊕}. Both planets orbit the slightly metal-poor, inactive K0V star HD 216520. Our updated best-fit models for HD 180617 b and GJ 686 b are in good agreement with the published results. For HD 180617 b, we obtain P = 105.91 days and M {sub pl} = 12.214 ± 1.05 M {sub ⊕}. For GJ 686 b, we find P = 15.53 days and M {sub pl} = 6.624 ± 0.432 M {sub ⊕}. Using an injection-recovery exercise, we find that HD 190007 b and HD 216520 b are unlikely to have additional planets with masses and orbital periods within a factor of 2, in marked contrast to ∼85% of planets in this mass and period range discovered by Kepler.

79 ASTRONOMY AND ASTROPHYSICS↗

An Intermediate-age Alpha-rich Galactic Population in K2

We explore the relationships between the chemistry, ages, and locations of stars in the Galaxy using asteroseismic data from the K2 mission and spectroscopic data from the Apache Point Galactic Evolution Experiment survey. Previous studies have used giant stars in the Kepler field to map the relationship between the chemical composition and the ages of stars at the solar circle. Consistent with prior work, we find that stars with high [α/Fe] have distinct, older ages in comparison to stars with low [α/Fe]. We provide age estimates for red giant branch (RGB) stars in the Kepler field, which support and build upon previous age estimates by taking into account the effect of α-enrichment on opacity. Including this effect for [α/Fe]-rich stars results in up to 10% older ages for low-mass stars relative to corrected solar mixture calculations. This is a significant effect that Galactic archeology studies should take into account. Looking beyond the Kepler field, we estimate ages for 735 RGB stars from the K2 mission, mapping age trends as a function of the line of sight. We find that the age distributions for low- and high-[α/Fe] stars converge with increasing distance from the Galactic plane, in agreement with suggestions from earlier work. We find that K2 stars with high [α/Fe] appear to be younger than their counterparts in the Kepler field, overlapping more significantly with a similarly aged low-[α/Fe] population. This observation may suggest that star formation or radial migration proceeds unevenly in the Galaxy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fundamental Parameters of the Eclipsing Binary DD CMa and Evidence for Mass Exchange

We present a detailed photometric and spectroscopic analysis of DD CMa, based on published survey photometry and new spectroscopic data. We find an improved orbital period of P {sub o} = 2.0084530(6) days. Our spectra reveal Hβ and Hα absorptions with weak emission shoulders, and we also find a color excess in the Wide-field Infrared Survey Explorer multiband photometry, interpreted as signatures of circumstellar matter. We model the V-band orbital light curve derived from the ASAS and ASAS-SN surveys, assuming a semidetached configuration and using the mass ratio and temperature of the hotter star derived from our spectroscopic analysis. Our model indicates that the system consists of a B2.5 dwarf and a B9 giant of radii 3.2 and 3.7 R {sub ⊙}, respectively, orbiting in a circular orbit of radius 6.75 R {sub ⊙}. We also found M {sub c} = 1.7 ± 0.1 M {sub ⊙}, T {sub c} = 11,350 ± 100 K, and M {sub h} = 6.4 ± 0.1 M {sub ⊙}, T {sub h} = 20,000 ± 500 K, for the cooler and hotter star, respectively. We find broad single emission peaks in Hα and Hβ after subtracting the synthetic stellar spectra. Our results are consistent with mass exchange between the stars and suggest the existence of a stream of gas being accreted onto the early B-type star.

47 OTHER INSTRUMENTATION↗

The Rarity of Very Red Trans-Neptunian Objects in the Scattered Disk

We investigate the origins of the photometrically very red (VR) and less red (LR) trans-Neptunian objects (TNOs). We first reanalyze the data set of Marsset et al. and find that in addition to the known color–inclination correlation in hot TNOs, a similar trend exists for color–eccentricity. We show that VR TNOs are sharply constrained to eccentricities <0.42 and inclinations <21°, leading to a paucity of VR scattered disk and distant mean motion resonance objects. We then interpret these findings using N-body simulations accounting for Neptune’s outward migration into a massless particles disk and find that these observations are best reproduced with an LR-to-VR color transition line between ∼38 and 42 au in the primordial disk, separating the objects’ formation locations. For an initial surface density profile (Σ ∝ 1/r {sup 2}), a color transition around 38 au is needed to explain the high abundance of VR plutinos, but it creates too many VR scattered disk objects, while a transition line around 42 au seems to better reproduce the scattered disk colors but creates virtually no VR plutinos. Our simulations furthermore show that the rarity of VR particles at high eccentricity is possibly due to the absence of sweeping higher-order MMRs, and secular resonances, beyond 42 au. Inspecting individual populations, we show that the majority of VR SDOs originate as objects trapped in Neptune’s second- and third-order MMRs. These then evolve due to diffusion, scattering, Kozai–Lidov cycles, and secular resonances into their current orbits. Future unbiased color surveys are crucial to better constrain the TNOs dynamical origins.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Far-ultraviolet to FIR Spectral-energy Distribution Modeling of the Stellar Formation History of the M31 Bulge

M31 is being surveyed at far- and near-ultraviolet with the UVIT telescope on AstroSat. The central bulge of M31 was observed in the N279N (275–280 nm), N219M (200–240 nm), F172M (160–185 nm), F169M (145–175 nm), and F148W (120–180 nm) filters. These images are made publicly available here. The UVIT data are supplemented with Sloan Digital Sky Survey data in optical, Spitzer data in near-infrared and Herschel data in mid- and far-infrared. The resulting far-ultraviolet to IR spectral-energy distributions for the bulge and for 10 subregions, are modeled using combinations of simple stellar populations and with CIGALE models. We find a dominant old (10–12 Gyr) metal-rich ([Z/H] ~ 0.3) population and a younger (600 Myr) solar abundance ([Z/H] ~ 0) population throughout the bulge. For the innermost 120" we find an additional very young (25 Myr) metal-poor ([Z/H] ~ –0.7) population. The results are consistent with the most recent stellar population studies of the bulge, which find the two populations for the whole bulge and a third young population in the innermost bulge.

79 ASTRONOMY AND ASTROPHYSICS↗

The Galactic Distribution of Phosphorus: A Survey of 163 Disk and Halo Stars*

Phosphorus (P) is a critical element for life on Earth, yet the cosmic production sites of P are relatively uncertain. To understand how P has evolved in the solar neighborhood, we measured abundances for 163 FGK stars over a range of –1.09 < [Fe/H] < 0.47 using observations from the Habitable-zone Planet Finder instrument on the Hobby–Eberly Telescope. Atmospheric parameters were calculated by fitting a combination of astrometry, photometry, and Fe I line equivalent widths. Phosphorus abundances were measured by matching synthetic spectra to a P I feature at 10529.52 Å. Our [P/Fe] ratios show that chemical evolution models generally underpredict P over the observed metallicity range. Additionally, we find that the [P/Fe] differs by ~0.1 dex between thin disk and thick disk stars that were identified with kinematics. The P abundances were compared with α-elements, iron-peak, odd-Z, and s-process elements, and we found that the evolution of P in the disk most strongly resembles that of the α-elements. We also find that molar P/C and N/C ratios for our sample match the scatter seen from other abundance studies. Finally, we measure a [P/Fe] = 0.09 ± 0.1 ratio in one low-α halo star and probable Gaia–Sausage–Enceladus member, an abundance ratio ~0.3–0.5 dex lower than the other Milky Way disk and halo stars at similar metallicities. Overall, we find that P is likely most significantly produced by massive stars in core-collapse supernovae, based on the largest P abundance survey to date.

79 ASTRONOMY AND ASTROPHYSICS↗

A Search for Predicted Astrometric Microlensing Events by Nearby Brown Dwarfs*

Gravitational microlensing has the potential to provide direct gravitational masses of single, free-floating brown dwarfs, independent of evolutionary and atmospheric models. The proper motions and parallaxes of nearby brown dwarfs can be used to predict close future alignments with distant background stars that cause a microlensing event. Targeted astrometric follow up of the predicted microlensing events permits the brown dwarf’s mass to be measured. Predicted microlensing events are typically found via searching for a peak threshold signal using an estimate of the lens mass. We develop a novel method that finds predicted events that instead will lead to a target lens-mass precision. The main advantage of our method is that it does not require a lens-mass estimate. We use this method to search for predicted astrometric microlensing events occurring between 2014 and 2032 using a catalog of 1225 low-mass star and brown-dwarf lenses in the Solar Neighborhood of spectral type M6 or later and a background source catalog from DECaLS Data Release 9. The background source catalog extends to g = 23.95, providing a more dense catalog compared to Gaia. Our search did not reveal any upcoming microlensing events. We estimate the rate of astrometric microlensing event for brown dwarfs in the Legacy Survey and find it to be low ~10 -5 yr -1 . We recommend carrying out targeted searches for brown dwarfs in front of the Galactic Bulge and Plane to find astrometric microlensing events that will allow the masses of single, free-floating brown dwarfs to be measured.

79 ASTRONOMY AND ASTROPHYSICS↗

The Dark Energy Survey Supernova Program: Corrections on Photometry Due to Wavelength-dependent Atmospheric Effects

Abstract Wavelength-dependent atmospheric effects impact photometric supernova flux measurements for ground-based observations. We present corrections on supernova flux measurements from the Dark Energy Survey Supernova Program’s 5YR sample (DES-SN5YR) for differential chromatic refraction (DCR) and wavelength-dependent seeing, and we show their impact on the cosmological parameters w and Ω m . We use g − i colors of Type Ia supernovae to quantify astrometric offsets caused by DCR and simulate point-spread functions (PSFs) using the GalSIM package to predict the shapes of the PSFs with DCR and wavelength-dependent seeing. We calculate the magnitude corrections and apply them to the magnitudes computed by the DES-SN5YR photometric pipeline. We find that for the DES-SN5YR analysis, not accounting for the astrometric offsets and changes in the PSF shape cause an average bias of +0.2 mmag and −0.3 mmag, respectively, with standard deviations of 0.7 mmag and 2.7 mmag across all DES observing bands ( griz ) throughout all redshifts. When the DCR and seeing effects are not accounted for, we find that w and Ω m are lower by less than 0.004 ± 0.02 and 0.001 ± 0.01, respectively, with 0.02 and 0.01 being the 1 σ statistical uncertainties. Although we find that these biases do not limit the constraints of the DES-SN5YR sample, future surveys with much higher statistics, lower systematics, and especially those that observe in the u band will require these corrections as wavelength-dependent atmospheric effects are larger at shorter wavelengths. We also discuss limitations of our method and how they can be better accounted for in future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

X-Ray Emission of Ultraviolet Variable Active Galactic Nucleus Candidates

We analyze the X-ray properties of a sample of 23 high-probability active galactic nucleus (AGN) candidates with ultraviolet variability identified in Wasleske et al. Using data from the Chandra X-ray Observatory and XMM-Newton Observatory, we that find 11/23 nuclei are X-ray detected. We use spectral energy distribution modeling to compute star formation rates and show that the X-ray luminosities are typically in excess of the X-ray emission expected from star formation by at least an order of magnitude. Interestingly, this sample shows a diversity of optical spectroscopic properties. We explore possible reasons for why some objects lack optical spectroscopic signatures of black hole activity while still being UV-variable and X-ray bright. We find that host galaxy stellar emission and obscuration from gas and dust are potential factors. We study where this sample falls on relationships such as α OX - L 2500 and L X - L IR and find that some of the sample falls outside the typical scatter for these relations, indicating that they differ from the standard quasar population. With the diversity of optical spectroscopic signatures and varying impacts of dust and stellar emissions on our sample, these results emphasize the strength of variability in selecting the most complete set of AGN, regardless of other host galaxy properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Signatures of a Tidally Induced Spiral Arm at the Anticenter of the Milky Way and a Kinematically Extended Anticenter Stream Using DESI Data Release 2

Using the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey, we examine the six-dimensional space of the anticenter region of the Milky Way stellar disk (150° < Galactic longitude < 220°) using 61,883 main-sequence turnoff stars. We focus on two well-known stellar overdensities in the anticenter: the Monoceros Ring (MRi) and Anticenter Stream (ACS). We find that the MRi overdensity has kinematic signatures consistent with a tidally induced spiral arm, a type of dynamic spiral arm created by an interaction with a satellite galaxy, most likely the Sagittarius dwarf spheroidal galaxy (Sgr). We use the kinematics of the MRi to calculate the two most recent passage times of Sgr, finding 0.25 ± 0.09 Gyr and 1.10 ± 0.23 Gyr from the present day. We validate that the ACS is kinematically decoupled from the MRi because they are moving in opposite radial and vertical directions. We find that the kinematics associated with the ACS extends beyond our defined overdensity. The features we see in the ACS region are likely part of a broader distribution of stars with the same kinematic signature as detected in other places, like the vertical wave in the outer disk and phase spiral.

Lambert, Mika [University of California, Santa Cru↗

The X-Ray Fundamental Plane of the Platinum Sample, the Kilonovae, and the SNe Ib/c Associated with GRBs

A large fraction of gamma-ray burst (GRB) lightcurves (LCs) show X-ray plateaus. We analyze all GRBs with known redshifts presenting plateaus observed by The Neil Gehrels Swift Observatory from its launch until 2019 August. The fundamental plane relation between the rest-frame time and X-ray luminosity at the end of the plateau emission and the peak prompt luminosity holds for all the GRB classes when selection biases and cosmological evolutions are applied. We have discovered two important findings: (1) a new class of long GRBs with good data coverage: the platinum sample; and (2) the platinum, the SNe-LGRB, and the KN-SGRB samples, yield the smallest intrinsic scatter with σ platinum , GRB-SNe = 0.22 ± 0.10 and σ KN-SGRB = 0.24 ± 0.12. The SNe-LGRBs are composed of GRBs associated spectroscopically with the SNe Ib,c, the KN-SGRBs are composed by eight GRBs associated with kilonovae or where there could have been such an association. The highest correlation coefficients are yielded for the SN-LGRB-ABC sample, which includes GRBs spectroscopically associated with SNe Ib/c or with a clear optical bump in the LC resembling the SNe Ib/c, ( ), the SN-LGRBs ( ), and the KN-SGRBs ( ) when the redshift evolution is considered. We find that these category planes are reliable candidates to use as cosmological tools. Furthermore, the distance from the gold fundamental plane is a crucial discriminant among classes. In fact, we find that the distributions of the distances of the SNe-LGRB, SNe-LGRB-ABC, KN-SGRB, and SGRB samples from the gold fundamental plane are statistically different from the distribution of the gold GRBs' distances from the gold fundamental plane with and without considering evolution cases. <

79 ASTRONOMY AND ASTROPHYSICS↗

A Template-based Approach to the Photometric Classification of SN 1991bg-like Supernovae in the SDSS-II Supernova Survey

The use of SNe Ia to measure cosmological parameters has grown significantly over the past two decades. However, there exists a significant diversity in the SN Ia population that is not well understood. Overluminous SN 1991T-like and subluminous SN 1991bg-like objects are two characteristic examples of peculiar SNe. The identification and classification of such objects is an important step in studying what makes them unique from the remaining SN population. With the upcoming Vera C. Rubin Observatory promising on the order of a million new SNe over a 10 year survey, spectroscopic classifications will be possible for only a small subset of observed targets. As such, photometric classification has become an increasingly important concern in preparing for the next generation of astronomical surveys. Using observations from the Sloan Digital Sky Survey II (SDSS-II) SN Survey, we apply here an empirically based classification technique targeted at the identification of SN 1991bg-like SNe in photometric data sets. By performing dedicated fits to photometric data in the rest-frame redder and bluer bandpasses, we classify 16 previously unidentified 91bg-like SNe. Using SDSS-II host galaxy measurements, we find that these SNe are preferentially found in host galaxies with an older average stellar age than the hosts of normal SNe Ia. We also find that these SNe are found at a further physical distance from the center of their host galaxies. We find no statistically significant bias in host galaxy mass or specific star formation rate for these targets.

Astronomy & Astrophysics↗

The Distant, Galaxy Cluster Environment of the Short GRB 161104A at $z\sim 0.8$ and a Comparison to the Short GRB Host Population

We present optical observations of the Swift short-duration gamma-ray burst (GRB) GRB 161104A and its host galaxy at z = 0.793 ± 0.003. We model the multiband photometry and spectroscopy with the stellar population inference code Prospector and explore the posterior using nested sampling. We find a mass-weighted age of $t_m = 2.12^{+0.23}_{-0.21}$ Gyr, stellar mass of log(M/M ⊙ ), metallicity of log(Z/Z ⊙ ), dust extinction of $A_V = 0.08^{+0.08}_{-0.05}$, and low star formation rate of 9.9 × 10 -2 M ⊙ yr -1 . These properties, along with a prominent 4000 Å break and optical absorption lines, classify this host as an early-type, quiescent galaxy. Using Dark Energy Survey galaxy catalogs, we demonstrate that the host of GRB 161104A resides on the outskirts of a galaxy cluster at z ≈ 0.8, situated ≈1 Mpc from the likely brightest cluster galaxy. We also present new modeling for 20 additional short GRB hosts (≈33% of which are early-type galaxies), finding population medians of and $t_m = 1.07^{+1.98}_{-0.67}$ Gyr (68% confidence). We further find that the host of GRB 161104A is more distant, less massive, and younger than the four other short GRB hosts known to be associated with galaxy clusters. Cluster short GRBs have faint afterglows, in the lower ≈11% (≈30%) of observed X-ray (optical) luminosities. We place a lower limit on the fraction of short GRBs in galaxy clusters versus those in the field of ≈5%–13%, consistent with the fraction of stellar mass of ≈10%–20% in galaxy clusters at redshifts 0.1 ≤ z ≤ 0.8. Future studies that take advantage of wider-field and deeper cluster surveys are needed to understand the true rate of short GRBs in clusters and their effect on heavy-element enrichment in the intracluster medium.

79 ASTRONOMY AND ASTROPHYSICS↗

Kilonova Luminosity Function Constraints Based on Zwicky Transient Facility Searches for 13 Neutron Star Merger Triggers during O3

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run (O3). We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization area of 4480 deg 2 , a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10 -25 yr -1 . The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (-16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than-16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day -1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than <57% (<89%) of putative kilonovae could be brighter than -16.6 mag assuming flat evolution (fading by 1 mag day -1 ) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than <68% of putative kilonovae could be brighter than -16.6 mag. Comparing to model grids, we find that some kilonovae must have M ej < 0.03 M ⊙ , X lan > 10 -4 , or φ > 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of-16 mag would constrain the maximum fraction of bright kilonovae to <25%.

79 ASTRONOMY AND ASTROPHYSICS↗