-
Bifurcations in the Herd Immunity Threshold for Discrete-Time Models of Epidemic Spread
Authors:
Sinan A. Ozbay,
Bjarke F. Nielsen,
Maximilian M. Nguyen
Abstract:
We performed a thorough sensitivity analysis of the herd immunity threshold for discrete-time SIR compartmental models with a static network structure. We find unexpectedly that these models violate classical intuition which holds that the herd immunity threshold should monotonically increase with the transmission parameter. We find the existence of bifurcations in the herd immunity threshold in t…
▽ More
We performed a thorough sensitivity analysis of the herd immunity threshold for discrete-time SIR compartmental models with a static network structure. We find unexpectedly that these models violate classical intuition which holds that the herd immunity threshold should monotonically increase with the transmission parameter. We find the existence of bifurcations in the herd immunity threshold in the high transmission probability regime. The extent of these bifurcations is modulated by the graph heterogeneity, the recovery parameter, and the network size. In the limit of large, well-mixed networks, the behavior approaches that of difference equation models, suggesting this behavior is a universal feature of all discrete-time SIR models. These results suggest careful attention is needed in both selecting the assumptions on how to model time and heterogeneity in epidemiological models and the subsequent conclusions that can be drawn.
△ Less
Submitted 24 February, 2023; v1 submitted 13 December, 2022;
originally announced December 2022.
-
Parameterizing Network Graph Heterogeneity using a Modified Weibull Distribution
Authors:
Sinan A. Ozbay,
Maximilian M. Nguyen
Abstract:
We present a simple method to quantitatively capture the heterogeneity in the degree distribution of a network graph using a single parameter $σ$. Using an exponential transformation of the shape parameter of the Weibull distribution, this control parameter allows the degree distribution to be easily interpolated between highly symmetric and highly heterogeneous distributions on the unit interval.…
▽ More
We present a simple method to quantitatively capture the heterogeneity in the degree distribution of a network graph using a single parameter $σ$. Using an exponential transformation of the shape parameter of the Weibull distribution, this control parameter allows the degree distribution to be easily interpolated between highly symmetric and highly heterogeneous distributions on the unit interval. This parameterization of heterogeneity also recovers several other canonical distributions as intermediate special cases, including the Gaussian, Rayleigh, and exponential distributions. We then outline a general graph generation algorithm to produce graphs with a desired amount of heterogeneity. The utility of this formulation of a heterogeneity parameter is demonstrated with examples relating to epidemiological modeling and spectral analysis.
△ Less
Submitted 20 December, 2022; v1 submitted 13 December, 2022;
originally announced December 2022.
-
APLC-Optimization: an apodized pupil Lyot coronagraph design survey toolkit
Authors:
Bryony F. Nickson,
Emiel H. Por,
Meiji M. Nguyen,
Remi Soummer,
Iva Laginja,
Ananya Sahoo,
Laurent Pueyo,
Kathryn St. Laurent,
Mamadou N'Diaye,
Neil T. Zimmerman,
James Noss,
Marshall Perrin
Abstract:
We present a publicly available software package developed for exploring apodized pupil Lyot coronagraph (APLC) solutions for various telescope architectures. In particular, the package optimizes the apodizer component of the APLC for a given focal-plane mask and Lyot stop geometry to meet a set of constraints (contrast, bandwidth etc.) on the coronagraph intensity in a given focal-plane region (i…
▽ More
We present a publicly available software package developed for exploring apodized pupil Lyot coronagraph (APLC) solutions for various telescope architectures. In particular, the package optimizes the apodizer component of the APLC for a given focal-plane mask and Lyot stop geometry to meet a set of constraints (contrast, bandwidth etc.) on the coronagraph intensity in a given focal-plane region (i.e. dark zone). The package combines a high-contrast imaging simulation package HCIPy with a third-party mathematical optimizer (Gurobi) to compute the linearly optimized binary mask that maximizes transmission. We provide examples of the application of this toolkit to several different telescope geometries, including the Gemini Planet Imager (GPI) and the High-contrast imager for Complex Aperture Telescopes (HiCAT) testbed. Finally, we summarize the results of a preliminary design survey for the case of a 6~m aperture off-axis space telescope, as recommended by the 2020 NASA Decadal Survey, exploring APLC solutions for different segment sizes. We then use the Pair-based Analytical model for Segmented Telescope Imaging from Space (PASTIS) to perform a segmented wavefront error tolerancing analysis on these solutions.
△ Less
Submitted 6 October, 2022;
originally announced October 2022.
-
GPI 2.0: Optical Designs for the Upgrade of the Gemini Planet Imager Coronagraphic system
Authors:
Meiji M. Nguyen,
Bryony F. Nickson,
Emiel H. Por,
Remi Soummer,
John G. Hagopian,
Bruce Macintosh,
Jeffrey Chilcote,
Laurent Pueyo,
Marshall Perrin,
Quinn Konopacky
Abstract:
The Gemini Planet Imager (GPI) is an integral field spectrograph (IFS) and coronagraph that is one of the few current generation instruments optimized for high-contrast direct imaging of substellar companions. The instrument is in the process of being upgraded and moved from its current mount on the Gemini South Observatory in Cerro Pachon, Chile, to its twin observatory, Gemini North, on Mauna Ke…
▽ More
The Gemini Planet Imager (GPI) is an integral field spectrograph (IFS) and coronagraph that is one of the few current generation instruments optimized for high-contrast direct imaging of substellar companions. The instrument is in the process of being upgraded and moved from its current mount on the Gemini South Observatory in Cerro Pachon, Chile, to its twin observatory, Gemini North, on Mauna Kea (a process colloquially dubbed 'GPI 2.0'). We present the designs that have been developed for the part of GPI 2.0 that pertains to upgrading various optical components of the GPI coronagraphic system. More specifically, we present new designs for the apodizer and Lyot stop (LS) that achieve better raw contrast at the inner working angle of the dark zone as well as improved core throughput while retaining a similar level of robustness to LS misalignment. To generate these upgraded designs, we use our own publicly available software package called APLC-Optimization that combines a commercial linear solver (Gurobi) with a high contrast imaging simulation package (HCIPy) in order to iteratively propagate light through a simulated model of an apodized phase lyot coronagraph (APLC), optimizing for the best coronagraph performance metrics. The designs have recently finished being lithographically printed by a commercial manufacturer and will be ready for use when GPI 2.0 goes on-sky in 2023.
△ Less
Submitted 26 September, 2022;
originally announced September 2022.
-
First detection of orbital motion for HD 106906 b: A wide-separation exoplanet on a Planet Nine-like orbit
Authors:
Meiji M. Nguyen,
Robert J. De Rosa,
Paul Kalas
Abstract:
HD 106906 is a 15 Myr old short-period (49 days) spectroscopic binary that hosts a wide-separation (737 au) planetary-mass ($\sim11\,M_{\rm Jup}$) common proper motion companion, HD 106906 b. Additionally, a circumbinary debris disk is resolved at optical and near-infrared wavelengths that exhibits a significant asymmetry at wide separations that may be driven by gravitational perturbations from t…
▽ More
HD 106906 is a 15 Myr old short-period (49 days) spectroscopic binary that hosts a wide-separation (737 au) planetary-mass ($\sim11\,M_{\rm Jup}$) common proper motion companion, HD 106906 b. Additionally, a circumbinary debris disk is resolved at optical and near-infrared wavelengths that exhibits a significant asymmetry at wide separations that may be driven by gravitational perturbations from the planet. In this study we present the first detection of orbital motion of HD 106906 b using Hubble Space Telescope images spanning a 14 yr period. We achieve high astrometric precision by cross-registering the locations of background stars with the Gaia astrometric catalog, providing the subpixel location of HD 106906 that is either saturated or obscured by coronagraphic optical elements. We measure a statistically significant $31.8\pm7.0$ mas eastward motion of the planet between the two most constraining measurements taken in 2004 and 2017. This motion enables a measurement of the inclination between the orbit of the planet and the inner debris disk of either $36_{-14}^{+27}$ deg or $44_{-14}^{+27}$ deg, depending on the true orientation of the orbit of the planet. There is a strong negative correlation between periastron and mutual inclination; orbits with smaller periastra are more misaligned with the disk plane. With a periastron of $510_{-320}^{+480}$ au, HD 106906 b is likely detached from the planetary region within 100 au radius, showing that a Planet Nine-like architecture can be established very early in the evolution of a planetary system.
△ Less
Submitted 8 December, 2020;
originally announced December 2020.
-
HD 165054: an astrometric calibration field for high-contrast imagers in Baade's Window
Authors:
Meiji M. Nguyen,
Robert J. De Rosa,
Jason J. Wang,
Thomas M. Esposito,
Paul Kalas,
James R. Graham,
Bruce Macintosh,
Vanessa P. Bailey,
Travis Barman,
Joanna Bulger,
Jeffrey Chilcote,
Tara Cotten,
Rene Doyon,
Gaspard Duchêne,
Michael P. Fitzgerald,
Katherine B. Follette,
Benjamin L. Gerard,
Stephen J. Goodsell,
Alexandra Z. Greenbaum,
Pascale Hibon,
Justin Hom,
Li-Wei Hung,
Patrick Ingraham,
Quinn Konopacky,
James E. Larkin
, et al. (29 additional authors not shown)
Abstract:
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'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…
▽ More
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'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.
△ Less
Submitted 6 April, 2020;
originally announced April 2020.
-
The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet beta Pictoris b from Combined Direct Imaging and Astrometry
Authors:
Eric L. Nielsen,
Robert J. De Rosa,
Jason J. Wang,
Johannes Sahlmann,
Paul Kalas,
Gaspard Duchene,
Julien Rameau,
Mark S. Marley,
Didier Saumon,
Bruce Macintosh,
Maxwell A. Millar-Blanchaer,
Meiji M. Nguyen,
S. Mark Ammons,
Vanessa P. Bailey,
Travis Barman,
Joanna Bulger,
Jeffrey Chilcote,
Tara Cotten,
Rene Doyon,
Thomas M. Esposito,
Michael P. Fitzgerald,
Katherine B. Follette,
Benjamin L. Gerard,
Stephen J. Goodsell,
James R. Graham
, et al. (29 additional authors not shown)
Abstract:
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…
▽ More
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'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.
△ Less
Submitted 25 November, 2019;
originally announced November 2019.
-
An updated visual orbit of the directly-imaged exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia
Authors:
Robert J. De Rosa,
Eric L. Nielsen,
Jason J. Wang,
S. Mark Ammons,
Gaspard Duchêne,
Bruce Macintosh,
Meiji M. Nguyen,
Julien Rameau,
Vanessa P. Bailey,
Travis Barman,
Joanna Bulger,
Jeffrey Chilcote,
Tara Cotten,
Rene Doyon,
Thomas M. Esposito,
Michael P. Fitzgerald,
Katherine B. Follette,
Benjamin L. Gerard,
Stephen J. Goodsell,
James R. Graham,
Alexandra Z. Greenbaum,
Pascale Hibon,
Justin Hom,
Li-Wei Hung,
Patrick Ingraham
, et al. (30 additional authors not shown)
Abstract:
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…
▽ More
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}$).
△ Less
Submitted 22 October, 2019;
originally announced October 2019.
-
Revised Astrometric Calibration of the Gemini Planet Imager
Authors:
Robert J. De Rosa,
Meiji M. Nguyen,
Jeffrey Chilcote,
Bruce Macintosh,
Marshall D. Perrin,
Quinn Konopacky,
Jason J. Wang,
Gaspard Duchêne,
Eric L. Nielsen,
Julien Rameau,
S. Mark Ammons,
Vanessa P. Bailey,
Travis Barman,
Joanna Bulger,
Tara Cotten,
Rene Doyon,
Thomas M. Esposito,
Michael P. Fitzgerald,
Katherine B. Follette,
Benjamin L. Gerard,
Stephen J. Goodsell,
James R. Graham,
Alexandra Z. Greenbaum,
Pascale Hibon,
Li-Wei Hung
, et al. (27 additional authors not shown)
Abstract:
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 redu…
▽ More
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.
△ Less
Submitted 8 April, 2020; v1 submitted 18 October, 2019;
originally announced October 2019.
