Researcher profile

Robert J. De Rosa

Robert J. De Rosa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
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

21 published item(s)

preprint2023arXiv

First on-sky results of ERIS at VLT

ERIS (Enhanced Resolution Imager and Spectrograph) is a new adaptive optics instrument installed at the Cassegrain focus of the VLT-UT4 telescope at the Paranal Observatory in Chile. ERIS consists of two near-infrared instruments: SPIFFIER, an integral field unit (IFU) spectrograph covering J to K bands, and NIX, an imager covering J to M bands. ERIS has an adaptive optics system able to work with both LGS and NGS. The Assembly Integration Verification (AIV) phase of ERIS at the Paranal Observatory was carried out starting in December 2021, followed by several commissioning runs in 2022. This contribution will describe the first preliminary results of the on-sky performance of ERIS during its commissioning and the future perspectives based on the preliminary scientific results.

preprint2022arXiv

Deep orbital search for additional planets in the HR 8799 system

The HR 8799 system hosts four massive planets orbiting 15 and 80 AU. Studies of the system's orbital stability and its outer debris disk open the possibility of additional planets, both interior to and exterior to the known system. Reaching a sufficient sensitivity to search for interior planets is very challenging due to the combination of bright quasi static speckle noise close to the stellar diffraction core and relatively fast orbital motion. In this work, we present a deep L-band imaging campaign using NIRC2 at Keck comprising 14 observing sequences. We further re-reduce archival data for a total of 16.75 hours, one of the largest uniform datasets of a single direct imaging target. Using a Bayesian modeling technique for detecting planets in images while compensating for plausible orbital motion, we then present deep limits on the existence of additional planets in the HR 8799 system. The final combination shows a tentative candidate, consistent with 4-7 $M_{jup}$ at 4-5 AU, detected with an equivalent false alarm probability better than $3σ$. This analysis technique is widely applicable to archival data and to new observations from upcoming missions that revisit targets at multiple epochs.

preprint2022arXiv

Nancy Grace Roman Space Telescope Coronagraph Instrument Observation Calibration Plan

NASA's next flagship mission, the Nancy Grace Roman Space Telescope, is a 2.4-meter observatory set to launch no later than May 2027. Roman features two instruments: the Wide Field Imager and the Coronagraph Instrument. Roman's Coronagraph is a Technology Demonstration that will push the current capabilities of direct imaging to smaller contrast ratios ($\sim$10$^{-9}$) and inner-working angles (3~$λ$/D). In order to achieve this high precision, Roman Coronagraph data must be calibrated to remove as many potential sources of error as possible. Here we present a detailed overview of the Nancy Grace Roman Space Telescope Coronagraph Instrument Observation Calibration Plan including identifying potential sources of error and how they will be mitigated via on-sky calibrations.

preprint2022arXiv

Search for Stellar Flybys in the Sco-Cen OB Association with the Gaia DR2

High-contrast imaging studies of debris disks have revealed a significant diversity in their morphologies, including large-scale asymmetries. Theories involving stellar flybys, an external source of gravitational disturbance, have offered a plausible explanation for the origin of these morphological variations. Our study is an experiment to gain empirical evidence that has been lacking from such theories. We explore this paradigm by using astrometric and radial velocity measurements from the Gaia DR2 and ground-based observations to trace the trajectories of 625 stars in the Sco-Cen OB Association from 5 Myr in the past to 2 Myr in the future. We identified 119 stars that had at least one past flyby event occurring within one Hill radius, and 23 of these experienced flybys within 0.5 Hill radius. We found no evidence of a significant correlation between the presence of flyby events and infrared excess detections, although the sample is not uniformly sensitive to infrared excess emission. Ten stars that had past flyby events host resolved circumstellar disks that appear relatively symmetric in the existing data except for the circumbinary disk surrounding HD 106906. We determined the trajectory and relative velocity of each of these flyby events, and compared these to the geometry of the spatially-resolved disks. Future work is needed to measure the kinematics of lower mass stars and to improve sensitivity to circumstellar disks for the entire sample.

preprint2020arXiv

A significant mutual inclination between the planets within the $π$ Mensae system

Measuring the geometry of multi-planet extrasolar systems can provide insight into their dynamical history and the processes of planetary formation. Such measurements are challenging for systems detected through indirect techniques such as radial velocity and transit, having only been measured for a handful of systems to-date. We aimed to place constraints on the orbital geometry of the outer planet in the $π$ Mensae system, a G0V star at 18.3 pc host to a wide-orbit super-jovian ($M\sin i = 10.02\pm0.15$ $M_{\rm Jup}$) with a 5.7-year period and an inner transiting super-earth ($M=4.82\pm0.85$ $M_\oplus$) with a 6.3-d period. We combined astrometric measurements from the Hipparcos and Gaia satellites with a precisely determined spectroscopic orbit in an attempt to constrain the inclination of the orbital plane of the outer planet. We measured an inclination of $i_b=49.9_{-4.5}^{+5.3}$ deg for the orbital plane of $π$ Mensae b, leading to a direct measurement of its mass of $13.01_{-0.95}^{+1.03}$ $M_{\rm Jup}$. We found a significant mutual inclination between the orbital planes of the two planets; a 95% credible interval for $i_{\rm mut}$ of between $34.5^\circ$ and $140.6^\circ$ after accounting for the unknown position angle of the orbit of $π$ Mensae c, strongly excluding a co-planar scenario for the two planets within this system. All orbits are stable in the present-day configuration, and secular oscillations of planet c's eccentricity are quenched by general relativistic precession. Planet c may have undergone high eccentricity tidal migration triggered by Kozai-Lidov cycles, but dynamical histories involving disk migration or in situ formation are not ruled out. Nonetheless, this system provides the first direct evidence that giant planets with large mutual inclinations have a role to play in the origins and evolution of some super-Earth systems.

preprint2020arXiv

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

We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager 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 three 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 but 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 non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 °$). Based on post-survey statistics, we improved upon our pre-survey 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.

preprint2020arXiv

HD 165054: an astrometric calibration field for high-contrast imagers in Baade's Window

We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade&#39;s Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers.

preprint2020arXiv

Keck/NIRC2 $L$&#39;-Band Imaging of Jovian-Mass Accreting Protoplanets around PDS 70

We present $L$&#39;-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets&#39; SEDs.

preprint2020arXiv

Multiband GPI Imaging of the HR 4796A Debris Disk

We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.

preprint2020arXiv

Multiband Polarimetric Imaging of HR 4796A with the Gemini Planet Imager

HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models.

preprint2020arXiv

Revised Astrometric Calibration of the Gemini Planet Imager

We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry.

preprint2020arXiv

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 (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15&#34;, or 20au. 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 $\approx$100au 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.

preprint2020arXiv

Understanding the Impacts of Stellar Companions on Planet Formation and Evolution: A Survey of Stellar and Planetary Companions within 25 pc

We explore the impact of outer stellar companions on the occurrence rate of giant planets detected with radial velocities. We searched for stellar and planetary companions to a volume-limited sample of solar-type stars within 25 pc. Using adaptive optics imaging from the Lick 3m and Palomar 200&#34; Telescopes, we characterized the multiplicity of our sample stars, down to the bottom of the main sequence. With these data, we confirm field star multiplicity statistics from previous surveys. We combined three decades of radial velocity data from the California Planet Search with new RV data from Keck/HIRES and APF/Levy to search for planets in the same systems. Using an updated catalog of both stellar and planetary companions and injection/recovery tests to determine our sensitivity, we measured the occurrence rate of planets among the single and multiple star systems. We found that planets with masses of 0.1-10 $M_{Jup}$ and semi-major axes of 0.1-10 AU have an occurrence rate of $0.18^{+0.04}_{-0.03}$ planets per single star, and $0.12\pm0.04$ planets per binary primary. Only one planet-hosting binary system in our sample had a binary separation $<100$ AU, and none had a separation $<50$ AU. We found planet occurrence rates of $0.20^{+0.07}_{-0.06}$ planets per star for binaries with separation $a_B > 100$ AU, and $0.04^{+0.04}_{-0.02}$ planets per star for binaries with separation $a_B<100$ AU. The similarity in the planet occurrence rate around single stars and wide primaries implies that wide binary systems should host more planets than single star systems, since they have more potential host stars. We estimated a system-wide planet occurrence rate of 0.3 planets per wide binary system for binaries with separations $a_B > 100$ AU. Finally, we found evidence that giant planets in binary systems have a different semi-major axis distribution than their counterparts in single star systems.

preprint2019arXiv

An updated visual orbit of the directly-imaged exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia

We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$σ$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$).

preprint2019arXiv

Atmospheric Characterization and Further Orbital Modeling of $κ$ And b

We present $κ$ Andromeda b&#39;s photometry and astrometry taken with Subaru/SCExAO+HiCIAO and Keck/NIRC2, combined with recently published SCExAO/CHARIS low-resolution spectroscopy and published thermal infrared photometry to further constrain the companion&#39;s atmospheric properties and orbit. $κ$ And b&#39;s Y/Y-K colors are redder than field dwarfs, consistent with its youth and lower gravity. Empirical comparisons of its Y-band photometry and CHARIS spectrum to a large spectral library of isolated field dwarfs reaffirm the conclusion from Currie et al. (2018) that it likely has a low gravity but admit a wider range of most plausible spectral types (L0-L2). Our gravitational classification also suggests that the best-fit objects for $κ$ And b may have lower gravity than those previously reported. Atmospheric models lacking dust/clouds fail to reproduce its entire 1--4.7 $μm$ spectral energy distribution, cloudy atmosphere models with temperatures of $\sim$ 1700--2000 $K$ better match $κ$ And b data. Most well-fitting model comparisons favor 1700--1900 $K$, a surface gravity of log(g) $\sim$ 4--4.5, and a radius of 1.3--1.6\,$R_{\rm Jup}$; the best-fit model (DRIFT-Phoenix) yields the coolest and lowest-gravity values: $T_{\rm eff}$=1700 K and $\log g$=4.0. An update to $κ$ And b&#39;s orbit with ExoSOFT using new astrometry spanning seven years reaffirms its high eccentricity ($0.77\pm0.08$). We consider a scenario where unseen companions are responsible for scattering $κ$ And b to a wide separation and high eccentricity. If three planets, including $κ$ And b, were born with coplanar orbits and one of them was ejected by gravitational scattering, a potential inner companion with mass $\gtrsim10M_{\rm Jup}$ could be located at $\lesssim$ 25 au.

preprint2019arXiv

Detection of a low-mass stellar companion to the accelerating A2IV star HR 1645

The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions.

preprint2019arXiv

First Resolved Scattered-Light Images of Four Debris Disks in Scorpius-Centaurus with the Gemini Planet Imager

We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems.

preprint2019arXiv

Model of the Search For Extraterrestrial Intelligence with Coronagraphic Imaging

We present modeled detection limits of the Gemini Planet Imager (GPI) and the Wide-Field Infrared Space Telescope (WFIRST) to an optical and infrared laser which could be used by an extraterrestrial civilization to signal their presence. GPI and WFIRST could utilize a coronagraph to search for extraterrestrial intelligence (SETI) in the present and future. We use archival data for GPI stars and simulated WFIRST observations to find the detectable flux ratio of a laser signal to residual scattered starlight around the target star. This flux ratio is then converted to detectable power as a function of distance from the parent star. For GPI, we assume a monochromatic laser wavelength of 1.55 $μ$m, and a wavelength of 575 nm for WFIRST. We assume the lasers are projected through a 10-m aperture, and that the intensity of the laser beam follows a Gaussian profile. Our analysis is performed on 6 stars with spectral types later than F within 20 pc (with an emphasis on solar analogs at different distances). The most notable result is the detection limit for $τ$ Ceti, a G5V star with four known exoplanets, two of those within the habitable zone (HZ). The result shows that a 24 kW laser is detectable from $τ$ Ceti from outside of the HZ with GPI and a 7.3 W laser is detectable from within $τ$ Ceti&#39;s HZ by WFIRST.

preprint2019arXiv

orbitize!: A Comprehensive Orbit-fitting Software Package for the High-contrast Imaging Community

orbitize! is an open-source, object-oriented software package for fitting the orbits of directly-imaged objects. It packages the Orbits for the Impatient (OFTI) algorithm and a parallel-tempered Markov Chain Monte Carlo (MCMC) algorithm into a consistent and intuitive Python API. orbitize! makes it easy to run standard astrometric orbit fits; in less than 10 lines of code, users can read in data, perform one fit using OFTI and another using MCMC, and make two publication-ready figures. Extensive pedagogical tutorials, intended to be navigable by both orbit-fitting novices and seasoned experts, are available on our documentation page. We have designed the orbitize! API to be flexible and easy to use/modify for unique cases. orbitize! was designed by members of the exoplanet imaging community to be a central repository for algorithms, techniques, and know-how developed by this community. We intend for it to continue to expand and change as the field progresses and new techniques are developed, and call for community involvement in this process. Complete and up-to-date documentation is available at orbitize.info.

preprint2019arXiv

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet beta Pictoris b from Combined Direct Imaging and Astrometry

We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet&#39;s Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.

preprint2019arXiv

The Possible Astrometric Signature of a Planetary-mass Companion to the Nearby Young Star TW Piscis Austrini (Fomalhaut B): Constraints from Astrometry, Radial Velocities, and Direct Imaging

We present constraints on the presence of substellar companions to the nearby ($d\sim7.6$ pc), young ($440\pm40$ Myr) K4Ve star TW Piscis Austrini, the wide ($\sim$0.3 pc) companion to the A4V star Fomalhaut. We combined absolute astrometry from Hipparcos and Gaia with literature radial velocity measurements and dedicated high-contrast imaging observations obtained with Keck/NIRC2 to achieve sensitivity to brown dwarf and planetary-mass companions ($\gtrsim2$ M$_{\rm Jup}$) over many decades of orbital period ($\lesssim 10^3$ yr). The significant astrometric acceleration measured between the Hipparcos and Gaia catalogues, reported previously in the literature, cannot be explained by the orbital motion of TW PsA around the barycenter of the Fomalhaut triple system. Instead, we find that it is consistent with the reflex motion induced by an orbiting substellar companion. The combination of astrometry, radial velocities, and a deep $L^{\prime}$ imaging dataset leads to a constraint on the companion mass of $1.2_{-0.6}^{+0.7}$ M$_{\rm Jup}$. However, the period of the companion is poorly constrained, with a highly multi-modal period posterior distribution due to aliasing with the 24.25 year baseline between Hipparcos and Gaia. If confirmed through continued astrometric or spectroscopic monitoring, or via direct detection, the companion to TW PsA would represent a choice target for detailed atmospheric characterization with high-contrast instruments on the upcoming James Webb Space Telescope and Wide Field Infrared Survey Telescope.