Source author record

David O. Cook

David O. Cook appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2020arXiv

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

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. 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 of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

preprint2019arXiv

GROWTH on S190425z: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

preprint2019arXiv

GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-Infrared Counterpart to a Neutron Star-Black Hole Merger

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $κ< 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ χ< 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

preprint2016arXiv

The Connection Between Galaxy Environment and the Luminosity Function Slopes of Star-Forming Regions

We present the first study of GALEX far ultra-violet (FUV) luminosity functions of individual star-forming regions within a sample of 258 nearby galaxies spanning a large range in total stellar mass and star formation properties. We identify ~65,000 star-forming regions (i.e., FUV sources), measure each galaxy's luminosity function, and characterize the relationships between the luminosity function slope (alpha) and several global galaxy properties. A final sample of 82 galaxies with reliable luminosity functions are used to define these relationships and represent the largest sample of galaxies with the largest range of galaxy properties used to study the connection between luminosity function properties and galaxy environment. We find that alpha correlates with global star formation properties, where galaxies with higher star formation rates and star formation rate densities (Sigma_SFR) tend to have flatter luminosity function slopes. In addition, we find that neither stochastic sampling of the luminosity function in galaxies with low-number statistics nor the effects of blending due to distance can fully account for these trends. We hypothesize that the flatter slopes in high Sigma_SFR galaxies is due to higher gas densities and higher star formation efficiencies which result in proportionally greater numbers of bright star-forming regions. Finally, we create a composite luminosity function composed of star-forming regions from many galaxies and find a break in the luminosity function at brighter luminosities. However, we find that this break is an artifact of varying detection limits for galaxies at different distances.

preprint2015arXiv

Hierarchical Star Formation across the ring galaxy NGC 6503

We present a detailed clustering analysis of the young stellar population across the star-forming ring galaxy NGC 6503, based on the deep HST photometry obtained with the Legacy ExtraGalactic UV Survey (LEGUS). We apply a contour-based map analysis technique and identify in the stellar surface density map 244 distinct star-forming structures at various levels of significance. These stellar complexes are found to be organized in a hierarchical fashion with 95% being members of three dominant super-structures located along the star-forming ring. The size distribution of the identified structures and the correlation between their radii and numbers of stellar members show power-law behaviors, as expected from scale-free processes. The self-similar distribution of young stars is further quantified from their autocorrelation function, with a fractal dimension of ~1.7 for length-scales between ~20 pc and 2.5 kpc. The young stellar radial distribution sets the extent of the star-forming ring at radial distances between 1 and 2.5 kpc. About 60% of the young stars belong to the detected stellar structures, while the remaining stars are distributed among the complexes, still inside the ring of the galaxy. The analysis of the time-dependent clustering of young populations shows a significant change from a more clustered to a more distributed behavior in a time-scale of ~60 Myr. The observed hierarchy in stellar clustering is consistent with star formation being regulated by turbulence across the ring. The rotational velocity difference between the edges of the ring suggests shear as the driving mechanism for this process. Our findings reveal the interesting case of an inner ring forming stars in a hierarchical fashion.

preprint2015arXiv

The Stellar Halo and Tidal Streams of Messier 63

We present new near-infrared (NIR) observations of M63 from the Extended Disk Galaxy Exploration Science (EDGES) Survey. The extremely deep 3.6 $μ$m mosaic reaches 29 AB mag arcsec$^{-2}$ at the outer reaches of the azimuthally-averaged surface brightness profile. At this depth the consequences of galactic accretion are found within a nearby tidal stream and an up-bending break in the slope of the surface brightness profile. This break occurs at a semi-major axis length of $\sim$8', and is evidence of either an enhanced outer disc or an inner stellar halo. Simulations of galaxy evolution, along with our observations, support an inner halo as the explanation for the up-bending break. The mass of this halo component is the largest found in an individual galaxy thus far. Additionally, our observations detect a nearby tidal stream. The mass of the stream suggests that a handful of such accretion events are necessary to populate the inner stellar halo. We also find that the accretion rate of the galaxy from the stream alone underestimates the accretion rate required to build M63's inner stellar halo.

preprint2014arXiv

Empirical ugri-UBVRc Transformations for Galaxies

We present empirical color transformations between Sloan Digital Sky Survey ugri and Johnson-Cousins UBVRc photometry for nearby galaxies (D < 11 Mpc). We use the Local Volume Legacy (LVL) galaxy sample where there are 90 galaxies with overlapping observational coverage for these two filter sets. The LVL galaxy sample consists of normal, non-starbursting galaxies. We also examine how well the LVL galaxy colors are described by previous transformations derived from standard calibration stars and model-based galaxy templates. We find significant galaxy color scatter around most of the previous transformation relationships. In addition, the previous transformations show systematic offsets between transformed and observed galaxy colors which are visible in observed color-color trends. The LVL-based $galaxy$ transformations show no systematic color offsets and reproduce the observed color-color galaxy trends.

preprint2014arXiv

Spitzer Local Volume Legacy (LVL) SEDs and Physical Properties

We present the panchromatic spectral energy distributions (SEDs) of the Local Volume Legacy (LVL) survey which consists of 258 nearby galaxies ($D<$11 Mpc). The wavelength coverage spans the ultraviolet to the infrared (1500 $\textrmÅ$ to 24 $μ$m) which is utilized to derive global physical properties (i.e., star formation rate, stellar mass, internal extinction due to dust.). With these data, we find color-color relationships and correlated trends between observed and physical properties (i.e., optical magnitudes and dust properties, optical color and specific star formation rate, and ultraviolet-infrared color and metallicity). The SEDs are binned by different galaxy properties to reveal how each property affects the observed shape of these SEDs. In addition, due to the volume-limited nature of LVL, we utilize the dwarf-dominated galaxy sample to test star formation relationships established with higher-mass galaxy samples. We find good agreement with the star-forming "main-sequence" relationship, but find a systematic deviation in the infrared "main-sequence" at low luminosities. This deviation is attributed to suppressed polycyclic aromatic hydrocarbon (PAH) formation in low metallicity environments and/or the destruction of PAHs in more intense radiation fields occurring near a suggested threshold in sSFR at a value of log($sSFR$) $\sim$ $-$10.2.

preprint2014arXiv

The Spitzer Local Volume Legacy (LVL) Global Optical Photometry

We present the global optical photometry of 246 galaxies in the Local Volume Legacy (LVL) survey. The full volume-limited sample consists of 258 nearby (D < 11 Mpc) galaxies whose absolute B-band magnitude span a range of -9.6 < M_B < -20.7 mag. A composite optical (UBVR) data set is constructed from observed UBVR and SDSS ugriz imaging, where the ugriz magnitudes are transformed into UBVR. We present photometry within three galaxy apertures defined at UV, optical, and IR wavelengths. Flux comparisons between these apertures reveal that the traditional optical R25 galaxy apertures do not fully encompass extended sources. Using the larger IR apertures we find color-color relationships where later-type spiral and irregular galaxies tend to be bluer than earlier-type galaxies. These data provide the missing optical emission from which future LVL studies can construct the full panchromatic (UV-optical-IR) spectral energy distributions.

preprint2012arXiv

The ACS Nearby Galaxy Survey Treasury. X. Quantifying the Star Cluster Formation Efficiency of Nearby Dwarf Galaxies

We study the relationship between the field star formation and cluster formation properties in a large sample of nearby dwarf galaxies. We use optical data from the Hubble Space Telescope and from ground-based telescopes to derive the ages and masses of the young (t_age < 100Myr) cluster sample. Our data provides the first constraints on two proposed relationships between the star formation rate of galaxies and the properties of their cluster systems in the low star formation rate regime. The data show broad agreement with these relationships, but significant galaxy-to-galaxy scatter exists. In part, this scatter can be accounted for by simulating the small number of clusters detected from stochastically sampling the cluster mass function. However, this stochasticity does not fully account for the observed scatter in our data suggesting there may be true variations in the fraction of stars formed in clusters in dwarf galaxies. Comparison of the cluster formation and the brightest cluster in our sample galaxies also provide constraints on cluster destruction models.

preprint2010arXiv

The Wyoming Survey for H-alpha. III. A Multi-wavelength Look at Attenuation by Dust in Galaxies out to z~0.4

We report results from the Wyoming Survey for H-alpha (WySH), a comprehensive four-square degree survey to probe the evolution of star-forming galaxies over the latter half of the age of the Universe. We have supplemented the H-alpha data from WySH with infrared data from the Spitzer Wide-area Infrared Extragalactic (SWIRE) Survey and ultraviolet data from the Galaxy Evolution Explorer (GALEX) Deep Imaging Survey. This dataset provides a multi-wavelength look at the evolution of the attenuation by dust, and here we compare a traditional measure of dust attenuation (L(TIR)/L(FUV)) to a diagnostic based on a recently-developed robust star formation rate (SFR) indicator, [H-alpha_obs+24-micron]/H-alpha_obs. With such data over multiple epochs, the evolution in the attenuation by dust with redshift can be assessed. We present results from the ELAIS-N1 and Lockman Hole regions at z~0.16, 0.24, 0.32 and 0.40. While the ensemble averages of both diagnostics are relatively constant from epoch to epoch, each epoch individually exhibits a larger attenuation by dust for higher star formation rates. Hence, an epoch to epoch comparison at a fixed star formation rate suggests a mild decrease in dust attenuation with redshift.

preprint2010arXiv

The Wyoming Survey for H-alpha. III. H-alpha Luminosity Functions at z ~ 0.16, 0.24, 0.32, and 0.40

The Wyoming Survey for H-alpha, or WySH, is a large-area, ground-based imaging survey for H-alpha-emitting galaxies at redshifts of z ~ 0.16, 0.24, 0.32, and 0.40. The survey spans up to four square degrees in a set of fields of low Galactic cirrus emission, using twin narrowband filters at each epoch for improved stellar continuum subtraction. H-alpha luminosity functions are presented for each Delta(z) ~ 0.02 epoch based on a total of nearly 1200 galaxies. These data clearly show an evolution with lookback time in the volume-averaged cosmic star formation rate. Integrals of Schechter fits to the incompleteness- and extinction-corrected H-alpha luminosity functions indicate star formation rates per co-moving volume of 0.010, 0.013, 0.020, 0.022 h_70 M_sun yr^{-1} Mpc^{-3} at z ~ 0.16, 0.24, 0.32, and 0.40, respectively. Statistical and systematic measurement uncertainties combined are on the order of 25% while the effects of cosmic variance are at the 20% level. The bulk of this evolution is driven by changes in the characteristic luminosity L_* of the H-alpha luminosity functions, with L_* for the earlier two epochs being a factor of two larger than L_* at the latter two epochs; it is more difficult with this data set to decipher systematic evolutionary differences in the luminosity function amplitude and faint-end slope. Coupling these results with a comprehensive compilation of results from the literature on emission line surveys, the evolution in the cosmic star formation rate density over 0 < z < 1.5 is measured to be rho_dot_SFR(z) = rho_dot_SFR(0) (1+z)^{3.4+/-0.4}.