BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Daniel Foreman-Mackey

Daniel Foreman-Mackey contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-ph.IMastro-ph.SRastro-ph.GA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

18 published item(s)

preprint2022arXiv

A Bayesian Population Model for the Observed Dust Attenuation in Galaxies

Dust plays a pivotal role in determining the observed spectral energy distribution (SED) of galaxies. Yet our understanding of dust attenuation is limited and our observations suffer from the dust-metallicity-age degeneracy in SED fitting (single galaxies), large individual variances (ensemble measurements), and the difficulty in properly dealing with uncertainties (statistical considerations). In this study, we create a population Bayesian model to rigorously account for correlated variables and non-Gaussian error distributions and demonstrate the improvement over a simple Bayesian model. We employ a flexible 5-D linear interpolation model for the parameters that control dust attenuation curves as a function of stellar mass, star formation rate (SFR), metallicity, redshift, and inclination. Our setup allows us to determine the complex relationships between dust attenuation and these galaxy properties simultaneously. Using Prospector fits of nearly 30,000 3D-HST galaxies, we find that the attenuation slope ($n$) flattens with increasing optical depth ($τ$), though less so than in previous studies. $τ$ increases strongly with SFR, though when $\log~{\rm SFR}\lesssim 0$, $τ$ remains roughly constant over a wide range of stellar masses. Edge-on galaxies tend to have larger $τ$ than face-on galaxies, but only for $\log~M_*\gtrsim 10$, reflecting the lack of triaxiality for low-mass galaxies. Redshift evolution of dust attenuation is strongest for low-mass, low-SFR galaxies, with higher optical depths but flatter curves at high redshift. Finally, $n$ has a complex relationship with stellar mass, highlighting the intricacies of the star-dust geometry. We have publicly released software (https://github.com/Astropianist/DustE) for users to access our population model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Aligned Orbit for the Young Planet V1298 Tau b

The alignment of planetary orbits with respect to the stellar rotation preserves information on their dynamical histories. Measuring this angle for young planets help illuminate the mechanisms that create misaligned orbits for older planets, as different processes could operate over timescales ranging from a few Myr to a Gyr. We present spectroscopic transit observations of the young exoplanet V1298 Tau b; we update the age of V1298 Tau to be $28\pm4$ Myr based on Gaia EDR3 measurements. We observed a partial transit with Keck/HIRES and LBT/PEPSI, and detected the radial velocity anomaly due to the Rossiter-McLaughlin effect. V1298 Tau~b has a prograde, well-aligned orbit, with $λ= 4_{-10}^{+7 \circ}$. By combining the spectroscopically-measured $v\sin i_{\star}$ and the phtometrically-measured rotation period of the host star we also find that the orbit is aligned in 3D, $ψ= 8_{-7}^{+4 \circ}$ deg. Finally, we combine our obliquity constraints with a previous measurement for the interior planet V1298 Tau c to constrain the mutual inclination between the two planets to be $i_{\mathrm{mut}}=0^{\circ} \pm 19^{\circ}$. This measurements adds to the growing number of well-aligned planets at young ages, hinting that misalignments may be generated over timescales of longer than tens of Myr. The number of measurements, however, is still small, and this population may not be representative of the older planets that have been observed to date. We also present the derivation of the relationship between $i_{\mathrm{mut}}$, $λ$, and $i$ for two planets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The 3D Galactocentric velocities of Kepler stars: marginalizing over missing RVs

Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e. 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to directly calculate 3D stellar velocities, however one can infer 3D stellar velocities by marginalizing over the missing RV dimension. In this paper, we infer the 3D velocities of stars in the Kepler field in Cartesian Galactocentric coordinates (vx, vy, vz). We directly calculate velocities for around a quarter of all Kepler targets, using RV measurements available from the Gaia, LAMOST and APOGEE spectroscopic surveys. Using the velocity distributions of these stars as our prior, we infer velocities for the remaining three-quarters of the sample by marginalizing over the RV dimension. The median uncertainties on our inferred vx, vy, and vz velocities are around 4, 18, and 4 km/s, respectively. We provide 3D velocities for a total of 148,590 stars in the Kepler field. These 3D velocities could enable kinematic age-dating, Milky Way stellar population studies, and other scientific studies using the benchmark sample of well-studied Kepler stars. Although the methodology used here is broadly applicable to targets across the sky, our prior is specifically constructed from and for the Kepler field. Care should be taken to use a suitable prior when extending this method to other parts of the Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Thresher: Lucky Imaging without the Waste

In traditional lucky imaging (TLI), many consecutive images of the same scene are taken with a high frame-rate camera, and all but the sharpest images are discarded before constructing the final shift-and-add image. Here we present an alternative image analysis pipeline -- The Thresher -- for these kinds of data, based on online multi-frame blind deconvolution. It makes use of all available data to obtain a best estimate of the astronomical scene in the context of reasonable computational limits; it does not require prior estimates of the point-spread functions in the images, or knowledge of point sources in the scene that could provide such estimates. Most importantly, the scene it aims to return is the optimum of a justified scalar objective based on the likelihood function. Because it uses the full set of images in the stack, The Thresher outperforms TLI in signal-to-noise; as it accounts for the individual-frame PSFs, it does this without loss of angular resolution. We demonstrate the effectiveness of our algorithm on both simulated data and real Electron-Multiplying CCD images obtained at the Danish 1.54m telescope (hosted by ESO, La Silla). We also explore the current limitations of the algorithm, and find that for the choice of image model presented here, non-linearities in flux are introduced into the returned scene. Ongoing development of the software can be viewed at https://github.com/jah1994/TheThresher.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The unpopular Package: a Data-driven Approach to De-trend TESS Full Frame Image Light Curves

The majority of observed pixels on the Transiting Exoplanet Survey Satellite (TESS) are delivered in the form of full frame images (FFI). However, the FFIs contain systematic effects such as pointing jitter and scattered light from the Earth and Moon that must be removed before downstream analysis. We present unpopular, an open-source Python package to de-trend TESS FFI light curves based on the causal pixel model method. Under the assumption that shared flux variations across multiple distant pixels are likely to be systematics, unpopular removes these common (i.e., popular) trends by modeling the systematics in a given pixel's light curve as a linear combination of light curves from many other distant pixels. To prevent overfitting we employ ridge regression and a train-and-test framework where the data points being de-trended are separated from those used to obtain the model coefficients. We also allow for simultaneous fitting with a polynomial model to capture any long-term astrophysical trends. We validate our method by de-trending different sources (e.g., supernova, tidal disruption event, exoplanet-hosting star, fast rotating star) and comparing our light curves to those obtained by other pipelines when appropriate. We also show that unpopular is able to preserve sector-length astrophysical signals, allowing for the extraction of multi-sector light curves from the FFI data. The unpopular source code and tutorials are freely available online.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

V1298 Tau with TESS: Updated Ephemerides, Radii, and Period Constraints from a Second Transit of V1298 Tau e

V1298 Tau is a young (20-30 Myr) solar-mass K star hosting four transiting exoplanets with sizes between 0.5 - 0.9 $R_J$. Given the system's youth, it provides a unique opportunity to understand the evolution of planetary radii at different separations in the same stellar environment. V1298 Tau was originally observed 6 years ago during K2 Campaign 4. Now, V1298 Tau has been re-observed during the extended mission of NASA's Transiting Exoplanet Survey Satellite (TESS). Here, we present new photometric observations of V1298 Tau from the 10-minute TESS Full-Frame Images. We use the TESS data to update the ephemerides for V1298 Tau bcde as well as compare newly observed radii to those measured from the K2 light curve, finding shallower transits for V1298 Tau bcd in the redder TESS bandpass at the $1-2σ$ level. We suspect the difference in radii is due to starspot-crossing events or contamination from nearby faint stars on the same pixels as V1298 Tau. Additionally, we catch a second transit of V1298 Tau e and present a new method for deriving the marginalized posterior probability of a planet's period from two transits observed years apart. We find the highest probability period for V1298 Tau e to be in a near 2:1 mean motion resonance with V1298 Tau b which, if confirmed, could make V1298 Tau bcde a 4 planet resonant chain. V1298 Tau is the target of several ongoing and future observations. These updated ephemerides will be crucial for scheduling future transit observations and interpreting future Doppler tomographic or transmission spectroscopy signals.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Analytic Light Curves in Reflected Light: Phase Curves, Occultations, and Non-Lambertian Scattering for Spherical Planets and Moons

We derive efficient, closed form, differentiable, and numerically stable solutions for the flux measured from a spherical planet or moon seen in reflected light, either in or out of occultation. Our expressions apply to the computation of scattered light phase curves of exoplanets, secondary eclipse light curves in the optical, or future measurements of planet-moon and planet-planet occultations, as well as to photometry of solar system bodies. We derive our solutions for Lambertian bodies illuminated by a point source, but extend them to model illumination sources of finite angular size and rough surfaces with phase-dependent scattering. Our algorithm is implemented in Python within the open-source starry mapping framework and is designed with efficient gradient-based inference in mind. The algorithm is 4-5 orders of magnitude faster than direct numerical evaluation methods and about 10 orders of magnitude more precise. We show how the techniques developed here may one day lead to the construction of two-dimensional maps of terrestrial planet surfaces, potentially enabling the detection of continents and oceans on exoplanets in the habitable zone.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multi-Wavelength Photometry Derived from Monochromatic Kepler Data

The Kepler mission has provided a wealth of data, revealing new insights in time-domain astronomy. However, Kepler's single band-pass has limited studies to a single wavelength. In this work we build a data-driven, pixel-level model for the Pixel Response Function (PRF) of Kepler targets, modeling the image data from the spacecraft. Our model is sufficiently flexible to capture known detector effects, such as non-linearity, intra-pixel sensitivity variations, and focus change. In theory, the shape of the Kepler PRF should also be weakly wavelength dependent, due to optical chromatic aberration and wavelength dependent detector response functions. We are able to identify these predicted shape changes to the PRF using the residuals between Kepler data and our model. In this work, we show that these PRF changes correspond to wavelength variability in Kepler targets using a small sample of eclipsing binaries. Using our model, we demonstrate that pixel-level light curves of eclipsing binaries show variable eclipse depths, ellipsoidal modulation and limb darkening. These changes at the pixel level are consistent with multi-wavelength photometry. Our work suggests each pixel in the Kepler data of a single target has a different effective wavelength, ranging from $\approx$ 550-750 $nm$. In this proof of concept, we demonstrate our model, and discuss possible use cases for the wavelength dependent Pixel Response Function of Kepler. These use cases include characterizing variable systems, and vetting exoplanet discoveries at the pixel level. The chromatic PRF of Kepler is due to weak wavelength dependence in the optical systems and detector of the telescope, and similar chromatic PRFs are expected in other similar telescopes, notably the NASA TESS telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Occultation mapping of Io's surface in the near-infrared I: Inferring static maps

Jupiter's moon Io is the most volcanically active body in the Solar System with hundreds of active volcanoes varying in intensity on different timescales. Io has been observed during occultations by other Galilean moons and Jupiter since the 1980s, using high-cadence near infrared photometry. These observations encode a wealth of information about the volcanic features on its surface. We built a generative model for the observed occultations using the code starry which enables fast, analytic, and differentiable computation of occultation light curves in emitted and reflected light. Our probabilistic Bayesian model is able to recover known hotspots on the surface of Io using only two light curves and without any assumptions on the locations, shapes or the number of spots. The methods we have developed are also directly applicable to the problem of mapping the surfaces of stars and exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Refining the transit timing and photometric analysis of TRAPPIST-1: Masses, radii, densities, dynamics, and ephemerides

We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST and K2 transit time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer light curves to derive the density of the host star and the planet densities. We find that all seven planets' densities may be described with a single rocky mass-radius relation which is depleted in iron relative to Earth, with Fe 21 wt% versus 32 wt% for Earth, and otherwise Earth-like in composition. Alternatively, the planets may have an Earth-like composition, but enhanced in light elements, such as a surface water layer or a core-free structure with oxidized iron in the mantle. We measure planet masses to a precision of 3-5%, equivalent to a radial-velocity (RV) precision of 2.5 cm/sec, or two orders of magnitude more precise than current RV capabilities. We find the eccentricities of the planets are very small; the orbits are extremely coplanar; and the system is stable on 10 Myr timescales. We find evidence of infrequent timing outliers which we cannot explain with an eighth planet; we instead account for the outliers using a robust likelihood function. We forecast JWST timing observations, and speculate on possible implications of the planet densities for the formation, migration and evolution of the planet system.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

starry_process: Interpretable Gaussian processes for stellar light curves

In this note we present the starry_process code, which implements an interpretable Gaussian process (GP) for modeling variability in stellar light curves. As dark starspots rotate in and out of view, the total flux received from a distant star will change over time. Unresolved flux time series therefore encode information about the spatial structure of features on the stellar surface. The starry_process software package allows one to easily model the flux variability due to starspots, whether one is interested in understanding the properties of these spots or marginalizing over the stellar variability when it is treated as a nuisance signal. The main difference between the GP implemented here and typical GPs used to model stellar variability is the explicit dependence of our GP on physical properties of the star, such as its period, inclination, and limb darkening coefficients, and on properties of the spots, such as their radius and latitude distributions. This code is the Python implementation of the interpretable GP algorithm developed in Luger, Foreman-Mackey, and Hedges (2021).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A planet within the debris disk around the pre-main-sequence star AU Microscopii

AU Microscopii (AU Mic) is the second closest pre main sequence star, at a distance of 9.79 parsecs and with an age of 22 million years. AU Mic possesses a relatively rare and spatially resolved3 edge-on debris disk extending from about 35 to 210 astronomical units from the star, and with clumps exhibiting non-Keplerian motion. Detection of newly formed planets around such a star is challenged by the presence of spots, plage, flares and other manifestations of magnetic activity on the star. Here we report observations of a planet transiting AU Mic. The transiting planet, AU Mic b, has an orbital period of 8.46 days, an orbital distance of 0.07 astronomical units, a radius of 0.4 Jupiter radii, and a mass of less than 0.18 Jupiter masses at 3 sigma confidence. Our observations of a planet co-existing with a debris disk offer the opportunity to test the predictions of current models of planet formation and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Exploring the evolution of stellar rotation using Galactic kinematics

The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Forward modeling the orbits of companions to pulsating stars from their light travel time variations

Mutual gravitation between a pulsating star and an orbital companion leads to a time-dependent variation in path length for starlight traveling to Earth. These variations can be used for coherently pulsating stars, such as the δ Scuti variables, to constrain the masses and orbits of their companions. Observing these variations for δ Scuti stars has previously relied on subdividing the light curve and measuring the average pulsation phase in equally sized subdivisions, which leads to under-sampling near periapsis. We introduce a new approach that simultaneously forward-models each sample in the light curve and show that this method improves upon current sensitivity limits - especially in the case of highly eccentric and short-period binaries. We find that this approach is sensitive enough to observe Jupiter mass planets around δ Scuti stars under ideal conditions, and use gravity-mode pulsations in the subdwarf B star KIC 7668647 to detect its companion without radial velocity data. We further provide robust detection limits as a function of the SNR of the pulsation mode and determine that the minimum detectable light travel time amplitude for a typical Kepler δ Scuti is around 2 s. This new method significantly enhances the application of light travel time variations to detecting short period binaries with pulsating components, and pulsating A-type exoplanet host stars, especially as a tool for eliminating false positives.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Assembly History of M87 Through Radial Variations in Chemical Abundances of its Field Star And Globular Cluster Populations

We present an extensive study of spectroscopically-derived chemical abundances for M87 and its globular cluster (GC) system. Using observations from the Mitchell spectrograph at McDonald, LRIS at Keck, and Hectospec on the MMT, we derive new metallicity gradients from $\sim 2$ to $140$ kpc. We use a novel hierarchical statistical framework to simultaneously separate the GC system into subpopulations while measuring the metallicity gradients of those subpopulations. We create physically-motivated spectral stacks of the GC subpopulations by leveraging the output of this statistical framework to perform the first application of abundance tagging in a massive ETG to better constrain the origins of the GC subpopulations and, thus, the assembly history of M87. We find a metal-poor, $α$-enhanced population of GCs in both in the inner and outer halo unanticipated by current cosmological simulations of galaxy evolution. We use the remarkably flat metallicity gradients we find for both the metal-rich and metal-poor GC subpopulations in the inner halo as tentative evidence that some amount of the metal-poor GCs formed directly in the halo of M87 at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Analytic Planetary Transit Light Curves and Derivatives for Stars with Polynomial Limb Darkening

We derive analytic, closed-form solutions for the light curve of a planet transiting a star with a limb darkening profile which is a polynomial function of the stellar elevation, up to arbitrary integer order. We provide improved analytic expressions for the uniform, linear, and quadratic limb-darkened cases, as well as novel expressions for higher order integer powers of limb darkening. The formulae are crafted to be numerically stable over the expected range of usage. We additionally present analytic formulae for the partial derivatives of instantaneous flux with respect to the radius ratio, impact parameter, and limb darkening coefficients. These expressions are rapid to evaluate, and compare quite favorably in speed and accuracy to existing transit light curve codes. We also use these expressions to numerically compute the first partial derivatives of exposure-time averaged transit light curves with respect to all model parameters. An additional application is modeling eclipsing binary or eclipsing multiple star systems in cases where the stars may be treated as spherically symmetric. We provide code which implements these formulae in C++, Python, IDL, and Julia, with tests and examples of usage.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

NGTS clusters survey. I. Rotation in the young benchmark open cluster Blanco 1

We determine rotation periods for 127 stars in the ~115 Myr old Blanco 1 open cluster using ~200 days of photometric monitoring with the Next Generation Transit Survey (NGTS). These stars span F5-M3 spectral types (1.2 $\gtrsim M \gtrsim$ 0.3 M$_{\odot}$) and increase the number of known rotation periods in Blanco 1 by a factor of four. We determine rotation periods using three methods: Gaussian process (GP) regression, generalised autocorrelation (G-ACF) and Lomb-Scargle (LS) periodograms, and find that GPs and G-ACF are more applicable to evolving spot modulation patterns. Between mid-F and mid-K spectral types, single stars follow a well-defined rotation sequence from ~2 to 10 days, whereas stars in photometric multiple systems typically rotate faster. This may suggest that the presence of a moderate-to-high mass ratio companion inhibits angular momentum loss mechanisms during the early pre-main sequence, and this signature has not been erased at ~100 Myr. The majority of mid-F to mid-K stars display evolving modulation patterns, whereas most M stars show stable modulation signals. This morphological change coincides with the shift from a well-defined rotation sequence (mid-F to mid-K stars) to a broad rotation period distribution (late-K and M stars). Finally, we compare our rotation results for Blanco 1 to the similarly-aged Pleiades: the single star populations in both clusters possess consistent rotation period distributions, which suggests that the angular momentum evolution of stars follows a well-defined pathway that is, at least for mid-F to mid-K stars, strongly imprinted by ~100 Myr.

0 citations0 reviewsNo code linkedOpen access