Engineering topics

Czekala, Ian

Publications and source records attributed to Czekala, Ian.

11 records

DOE OSTI2020

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

We report the results of a ∼4 yr direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager (GPI) Exoplanet Survey. We targeted nearby (≲150 pc), young (≲500 Myr) stars with high infrared (IR) excesses (L {sub IR}/L {sub ⋆} > 10{sup −5}), including 38 with previously resolved disks. Observations were made using the GPI high-contrast integral field spectrograph in H-band (1.6 μm) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and 3 protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks except HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of nondetections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low-inclination disks (≲70°). Based on postsurvey statistics, we improved upon our presurvey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing H-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

The Gemini Planet Imager View of the HD 32297 Debris Disk

We present new H-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager. The disk is detected in total and polarized intensity down to a projected angular separation of 0.″15, or 20 au. On the other hand, the large-scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of ≈100 au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high-fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics from 10 to 100 au

Here, we introduce a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey. This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semimajor axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M * > 1.5 M ⊙ more likely to host planets with masses between 2 and 13$M$ Jup and semimajor axes of 3–100 au at 99.92% confidence. We fit a double power-law model in planet mass ($m$) and semimajor axis ($a$) for planet populations around high-mass stars (M * > 1.5 M ⊙ ) of the form ${d}^{2}N/({dm}\,{da})\propto {m}^{\alpha }\,{a}^{\beta }$, finding α = –2.4 ± 0.8 and β = –2.0 ± 0.5, and an integrated occurrence rate of ${9}_{-4}^{+5}$% between 5–13M Jup and 10–100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with ${0.8}_{-0.5}^{+0.8}$% of stars hosting a brown dwarf companion between 13–80$M$ Jup and 10–100 au. Brown dwarfs also appear to be distributed differently in mass and semimajor axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semimajor axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the radial velocity method, our results are consistent with a peak in occurrence of giant planets between ~1 and 10 au. These trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability, according to this work.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Detecting Weak Spectral Lines in Interferometric Data through Matched Filtering

Modern radio interferometers enable observations of spectral lines with unprecedented spatial resolution and sensitivity. In spite of these technical advances, many lines of interest are still at best weakly detected and therefore necessitate detection and analysis techniques specialized for the low signal-to-noise ratio (S/N) regime. Matched filters can leverage knowledge of the source structure and kinematics to increase sensitivity of spectral line observations. Application of the filter in the native Fourier domain improves S/N while simultaneously avoiding the computational cost and ambiguities associated with imaging, making matched filtering a fast and robust method for weak spectral line detection. We demonstrate how an approximate matched filter can be constructed from a previously observed line or from a model of the source, and we show how this filter can be used to robustly infer a detection significance for weak spectral lines. When applied to ALMA Cycle 2 observations of CH{sub 3}OH in the protoplanetary disk around TW Hya, the technique yields a ≈53% S/N boost over aperture-based spectral extraction methods, and we show that an even higher boost will be achieved for observations at higher spatial resolution. A Python-based open-source implementation of this technique is available under the MIT license at http://github.com/AstroChem/VISIBLE.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Characterizing 51 Eri b from 1 to 5 μm: A Partly Cloudy Exoplanet

We present spectrophotometry spanning 1–5 μm of 51 Eridani b, a 2–10 M{sub Jup} planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new K1 (1.90–2.19 μm) and K2 (2.10–2.40 μm) spectra taken with the Gemini Planet Imager as well as an updated L {sub P} (3.76 μm) and new M {sub S} (4.67 μm) photometry from the NIRC2 Narrow camera. The new data were combined with J (1.13–1.35 μm) and H (1.50–1.80 μm) spectra from the discovery epoch with the goal of better characterizing the planet properties. The 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4 and T8), and we propose that 51 Eri b might be in the process of undergoing the transition from L-type to T-type. We used two complementary atmosphere model grids including either deep iron/silicate clouds or sulfide/salt clouds in the photosphere, spanning a range of cloud properties, including fully cloudy, cloud-free, and patchy/intermediate-opacity clouds. The model fits suggest that 51 Eri b has an effective temperature ranging between 605 and 737 K, a solar metallicity, and a surface gravity of log(g) = 3.5–4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the spectral energy distribution (SED). From the model atmospheres, we infer a luminosity for the planet of −5.83 to −5.93 (logL/L{sub ⊙}), leaving 51 Eri b in the unique position of being one of the only directly imaged planets consistent with having formed via a cold-start scenario. Comparisons of the planet SED against warm-start models indicate that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 M{sub ⊕}.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Improving and Assessing Planet Sensitivity of the GPI Exoplanet Survey with a Forward Model Matched Filter

We present a new matched-filter algorithm for direct detection of point sources in the immediate vicinity of bright stars. The stellar point-spread function (PSF) is first subtracted using a Karhunen-Loéve image processing (KLIP) algorithm with angular and spectral differential imaging (ADI and SDI). The KLIP-induced distortion of the astrophysical signal is included in the matched-filter template by computing a forward model of the PSF at every position in the image. To optimize the performance of the algorithm, we conduct extensive planet injection and recovery tests and tune the exoplanet spectra template and KLIP reduction aggressiveness to maximize the signal-to-noise ratio (S/N) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal S/N loss. We also developed a complete pipeline for the automated detection of point-source candidates, the calculation of receiver operating characteristics (ROC), contrast curves based on false positives, and completeness contours. We process in a uniform manner more than 330 data sets from the Gemini Planet Imager Exoplanet Survey and assess GPI typical sensitivity as a function of the star and the hypothetical companion spectral type. This work allows for the first time a comparison of different detection algorithms at a survey scale accounting for both planet completeness and false-positive rate. We show that the new forward model matched filter allows the detection of 50% fainter objects than a conventional cross-correlation technique with a Gaussian PSF template for the same false-positive rate.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Disentangling Time-series Spectra with Gaussian Processes: Applications to Radial Velocity Analysis

Measurements of radial velocity variations from the spectroscopic monitoring of stars and their companions are essential for a broad swath of astrophysics; these measurements provide access to the fundamental physical properties that dictate all phases of stellar evolution and facilitate the quantitative study of planetary systems. The conversion of those measurements into both constraints on the orbital architecture and individual component spectra can be a serious challenge, however, especially for extreme flux ratio systems and observations with relatively low sensitivity. Gaussian processes define sampling distributions of flexible, continuous functions that are well-motivated for modeling stellar spectra, enabling proficient searches for companion lines in time-series spectra. We introduce a new technique for spectral disentangling, where the posterior distributions of the orbital parameters and intrinsic, rest-frame stellar spectra are explored simultaneously without needing to invoke cross-correlation templates. To demonstrate its potential, this technique is deployed on red-optical time-series spectra of the mid-M-dwarf binary LP661-13. We report orbital parameters with improved precision compared to traditional radial velocity analysis and successfully reconstruct the primary and secondary spectra. We discuss potential applications for other stellar and exoplanet radial velocity techniques and extensions to time-variable spectra. The code used in this analysis is freely available as an open-source Python package.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Placing the Spotted T Tauri Star LkCa 4 on an HR Diagram

Ages and masses of young stars are often estimated by comparing their luminosities and effective temperatures to pre-main-sequence stellar evolution tracks, but magnetic fields and starspots complicate both the observations and evolution. To understand their influence, we study the heavily spotted weak-lined T-Tauri star LkCa 4 by searching for spectral signatures of radiation originating from the starspot or starspot groups. We introduce a new methodology for constraining both the starspot filling factor and the spot temperature by fitting two-temperature stellar atmosphere models constructed from Phoenix synthetic spectra to a high-resolution near-IR IGRINS spectrum. Clearly discernable spectral features arise from both a hot photospheric component T{sub hot} ∼ 4100 K and a cool component T{sub cool} ∼ 2700–3000 K, which covers ∼80% of the visible surface. This mix of hot and cool emission is supported by analyses of the spectral energy distribution, rotational modulation of colors and of TiO band strengths, and features in low-resolution optical/near-IR spectroscopy. Although the revised effective temperature and luminosity make LkCa 4 appear to be much younger and of much lower mass than previous estimates from unspotted stellar evolution models, appropriate estimates will require the production and adoption of spotted evolutionary models. Biases from starspots likely afflict most fully convective young stars and contribute to uncertainties in ages and age spreads of open clusters. In some spectral regions, starspots act as a featureless “veiling” continuum owing to high rotational broadening and heavy line blanketing in cool star spectra. Some evidence is also found for an anticorrelation between the velocities of the warm and cool components.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

ALMA MEASUREMENTS OF CIRCUMSTELLAR MATERIAL IN THE GQ LUP SYSTEM

We present Atacama Large Millimeter/submillimeter Array observations of the GQ Lup system, a young Sun-like star with a substellar-mass companion in a wide-separation orbit. These observations of 870 μ m continuum and CO J = 3–2 line emission with beam size ∼0.″3 (∼45 au) resolve the disk of dust and gas surrounding the primary star, GQ Lup A, and provide deep limits on any circumplanetary disk surrounding the companion, GQ Lup b. The circumprimary dust disk is compact with an FWHM of 59 ± 12 au, while the gas has a larger extent with a characteristic radius of 46.5 ± 1.8 au. By forward-modeling the velocity field of the circumprimary disk based on the CO emission, we constrain the mass of GQ Lup A to be M {sub *} = (1.03 ± 0.05) ∗ ( d /156 pc) M {sub ⊙}, where d is a known distance, and determine that we view the disk at an inclination angle of 60.°5 ± 0.°5 and a position angle of 346° ± 1°. The 3 σ upper limit on the 870 μ m flux density of any circumplanetary disk associated with GQ Lup b of <0.15 mJy implies an upper limit on the dust disk mass of <0.04 M {sub ⊕} for standard assumptions about optically thin emission. We discuss proposed mechanisms for the formation of wide-separation substellar companions given the non-detection of circumplanetary disks around GQ Lup b and other similar systems.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

THE COUPLED PHYSICAL STRUCTURE OF GAS AND DUST IN THE IM Lup PROTOPLANETARY DISK

The spatial distribution of gas and solids in protoplanetary disks determines the composition and formation efficiency of planetary systems. A number of disks show starkly different distributions for the gas and small grains compared to millimeter–centimeter-sized dust. We present new Atacama Large Millimeter/Submillimeter Array observations of the dust continuum, CO, {sup 13}CO, and C{sup 18}O in the IM Lup protoplanetary disk, one of the first systems where this dust–gas dichotomy was clearly seen. The {sup 12}CO is detected out to a radius of 970 au, while the millimeter continuum emission is truncated at just 313 au. Based upon these data, we have built a comprehensive physical and chemical model for the disk structure, which takes into account the complex, coupled nature of the gas and dust and the interplay between the local and external environment. We constrain the distributions of gas and dust, the gas temperatures, the CO abundances, the CO optical depths, and the incident external radiation field. We find that the reduction/removal of dust from the outer disk exposes this region to higher stellar and external radiation and decreases the rate of freeze-out, allowing CO to remain in the gas out to large radial distances. We estimate a gas-phase CO abundance of 5% of the interstellar medium value and a low external radiation field ( G {sub 0} ≲ 4). The latter is consistent with that expected from the local stellar population. We additionally find tentative evidence for ring-like continuum substructure, suggestions of isotope-selective photodissociation, and a diffuse gas halo.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

CONSTRUCTING A FLEXIBLE LIKELIHOOD FUNCTION FOR SPECTROSCOPIC INFERENCE

We present a modular, extensible likelihood framework for spectroscopic inference based on synthetic model spectra. The subtraction of an imperfect model from a continuously sampled spectrum introduces covariance between adjacent datapoints (pixels) into the residual spectrum. For the high signal-to-noise data with large spectral range that is commonly employed in stellar astrophysics, that covariant structure can lead to dramatically underestimated parameter uncertainties (and, in some cases, biases). We construct a likelihood function that accounts for the structure of the covariance matrix, utilizing the machinery of Gaussian process kernels. This framework specifically addresses the common problem of mismatches in model spectral line strengths (with respect to data) due to intrinsic model imperfections (e.g., in the atomic/molecular databases or opacity prescriptions) by developing a novel local covariance kernel formalism that identifies and self-consistently downweights pathological spectral line “outliers.” By fitting many spectra in a hierarchical manner, these local kernels provide a mechanism to learn about and build data-driven corrections to synthetic spectral libraries. An open-source software implementation of this approach is available at http://iancze.github.io/Starfish, including a sophisticated probabilistic scheme for spectral interpolation when using model libraries that are sparsely sampled in the stellar parameters. We demonstrate some salient features of the framework by fitting the high-resolution V-band spectrum of WASP-14, an F5 dwarf with a transiting exoplanet, and the moderate-resolution K-band spectrum of Gliese 51, an M5 field dwarf.

79 ASTRONOMY AND ASTROPHYSICS↗