Showing 1–50 of 149 results for author: Brewer, J

The Lowell Observatory Solar Telescope: A fiber feed into the EXtreme PREcision Spectrometer

Authors: Joe Llama, Lily L. Zhao, John M. Brewer, Andrew Szymkowiak, Debra A. Fischer, Michael Collins, Jake Tiegs, Frank Cornelius

Abstract: The signal induced by a temperate, terrestrial planet orbiting a Sun-like star is an order of magnitude smaller than the host stars' intrinsic variability. Understanding stellar activity is, therefore, a fundamental obstacle in confirming the smallest exoplanets. We present the Lowell Observatory Solar Telescope (LOST), a solar feed for the EXtreme PREcision Spectrometer (EXPRES) at the 4.3-m Lowe… ▽ More The signal induced by a temperate, terrestrial planet orbiting a Sun-like star is an order of magnitude smaller than the host stars' intrinsic variability. Understanding stellar activity is, therefore, a fundamental obstacle in confirming the smallest exoplanets. We present the Lowell Observatory Solar Telescope (LOST), a solar feed for the EXtreme PREcision Spectrometer (EXPRES) at the 4.3-m Lowell Discovery Telescope (LDT). EXPRES is one of the newest high-resolution spectrographs that accurately measure extreme radial velocity. With LOST/EXPRES, we observe disk-integrated sunlight autonomously throughout the day. In clear conditions, we achieve a ~137,500 optical spectrum of the Sun with a signal-to-noise of 500 in ~150s. Data is reduced using the standard EXPRES pipeline with minimal modification to ensure the data are comparable to the observations of other stars with the LDT. During the first three years of operation, we find a daily RMS of 71 cm/s. Additionally, having two EPRV spectrometers located in Arizona gives us an unprecedented opportunity to benchmark the performance of these planet-finders. We find a RMS of just 55 cm/s when comparing data taken simultaneously with EXPRES and NEID. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: SPIE Astronomical Telescopes & Instrumentation proceedings paper

Rotation of the Globular Cluster Population of the Dark Matter Deficient Galaxy NGC 1052-DF4: Implication for the Total Mass

Authors: Yuan, Li, Brendon J. Brewer, Geraint F. Lewis

Abstract: We explore the globular cluster population of NGC 1052-DF4, a dark matter deficient galaxy, using Bayesian inference to search for the presence of rotation. The existence of such a rotating component is relevant to the estimation of the mass of the galaxy, and therefore the question of whether NGC 1052-DF4 is truly deficient of dark matter, similar to NGC 1052-DF2 another galaxy in the same group.… ▽ More We explore the globular cluster population of NGC 1052-DF4, a dark matter deficient galaxy, using Bayesian inference to search for the presence of rotation. The existence of such a rotating component is relevant to the estimation of the mass of the galaxy, and therefore the question of whether NGC 1052-DF4 is truly deficient of dark matter, similar to NGC 1052-DF2 another galaxy in the same group. The rotational characteristics of seven globular clusters in NGC 1052-DF4 were investigated, finding that a non-rotating kinematic model has a higher Bayesian evidence than a rotating model, by a factor of approximately 2.5. In addition, we find that under the assumption of rotation, its amplitude must be small. This distinct lack of rotation strengthens the case that, based on its intrinsic velocity dispersion, NGC 1052-DF4 is a truly dark matter deficient galaxy. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 9 pages 6 figures. Accepted for publication in PASA

The Extreme Stellar-Signals Project III. Combining Solar Data from HARPS, HARPS-N, EXPRES, and NEID

Authors: Lily L. Zhao, Xavier Dumusque, Eric B. Ford, Joe Llama, Annelies Mortier, Megan Bedell, Khaled Al Moulla, Chad F. Bender, Cullen H. Blake, John M. Brewer, Andrew Collier Cameron, Rosario Cosentino, Pedro Figueira, Debra A. Fischer, Adriano Ghedina, Manuel Gonzalez, Samuel Halverson, Shubham Kanodia, David W. Latham, Andrea S. J. Lin, Gaspare Lo Curto, Marcello Lodi, Sarah E. Logsdon, Christophe Lovis, Suvrath Mahadevan , et al. (15 additional authors not shown)

Abstract: We present an analysis of Sun-as-a-star observations from four different high-resolution, stabilized spectrographs -- HARPS, HARPS-N, EXPRES, and NEID. With simultaneous observations of the Sun from four different instruments, we are able to gain insight into the radial velocity precision and accuracy delivered by each of these instruments and isolate instrumental systematics that differ from true… ▽ More We present an analysis of Sun-as-a-star observations from four different high-resolution, stabilized spectrographs -- HARPS, HARPS-N, EXPRES, and NEID. With simultaneous observations of the Sun from four different instruments, we are able to gain insight into the radial velocity precision and accuracy delivered by each of these instruments and isolate instrumental systematics that differ from true astrophysical signals. With solar observations, we can completely characterize the expected Doppler shift contributed by orbiting Solar System bodies and remove them. This results in a data set with measured velocity variations that purely trace flows on the solar surface. Direct comparisons of the radial velocities measured by each instrument show remarkable agreement with residual intra-day scatter of only 15-30 cm/s. This shows that current ultra-stabilized instruments have broken through to a new level of measurement precision that reveals stellar variability with high fidelity and detail. We end by discussing how radial velocities from different instruments can be combined to provide powerful leverage for testing techniques to mitigate stellar signals. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 17 pages, 9 figures, accepted for publication

Refining the Stellar Parameters of $τ$ Ceti: a Pole-on Solar Analog

Authors: Maria Korolik, Rachael M. Roettenbacher, Debra A. Fischer, Stephen R. Kane, Jean M. Perkins, John D. Monnier, Claire L. Davies, Stefan Kraus, Jean-Baptiste Le Bouquin, Narsireddy Anugu, Tyler Gardner, Cyprien Lanthermann, Gail H. Schaefer, Benjamin Setterholm, John M. Brewer, Joe Llama, Lily L. Zhao, Andrew E. Szymkowiak, Gregory W. Henry

Abstract: To accurately characterize the planets a star may be hosting, stellar parameters must first be well-determined. $τ$ Ceti is a nearby solar analog and often a target for exoplanet searches. Uncertainties in the observed rotational velocities have made constraining $τ$ Ceti's inclination difficult. For planet candidates from radial velocity (RV) observations, this leads to substantial uncertainties… ▽ More To accurately characterize the planets a star may be hosting, stellar parameters must first be well-determined. $τ$ Ceti is a nearby solar analog and often a target for exoplanet searches. Uncertainties in the observed rotational velocities have made constraining $τ$ Ceti's inclination difficult. For planet candidates from radial velocity (RV) observations, this leads to substantial uncertainties in the planetary masses, as only the minimum mass ($m \sin i$) can be constrained with RV. In this paper, we used new long-baseline optical interferometric data from the CHARA Array with the MIRC-X beam combiner and extreme precision spectroscopic data from the Lowell Discovery Telescope with EXPRES to improve constraints on the stellar parameters of $τ$ Ceti. Additional archival data were obtained from a Tennessee State University Automatic Photometric Telescope and the Mount Wilson Observatory HK project. These new and archival data sets led to improved stellar parameter determinations, including a limb-darkened angular diameter of $2.019 \pm 0.012$ mas and rotation period of $46 \pm 4$ days. By combining parameters from our data sets, we obtained an estimate for the stellar inclination of $7\pm7^\circ$. This nearly-pole-on orientation has implications for the previously-reported exoplanets. An analysis of the system dynamics suggests that the planetary architecture described by Feng et al. (2017) may not retain long-term stability for low orbital inclinations. Additionally, the inclination of $τ$ Ceti reveals a misalignment between the inclinations of the stellar rotation axis and the previously-measured debris disk rotation axis ($i_\mathrm{disk} = 35 \pm 10^\circ$). △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 14 pages, 3 figures, 4 tables, 1 appendix, accepted for publication to AJ

EXPRES IV: Two Additional Planets Orbiting $ρ$ Coronae Borealis Reveal Uncommon System Architecture

Authors: John M. Brewer, Lily L. Zhao, Debra A. Fischer, Rachael M. Roettenbacher, Gregory W. Henry, Joe Llama, Andrew E. Szymkowiak, Samuel H. C. Cabot, Sam A. Weiss, Chris McCarthy

Abstract: Thousands of exoplanet detections have been made over the last twenty-five years using Doppler observations, transit photometry, direct imaging, and astrometry. Each of these methods is sensitive to different ranges of orbital separations and planetary radii (or masses). This makes it difficult to fully characterize exoplanet architectures and to place our solar system in context with the wealth o… ▽ More Thousands of exoplanet detections have been made over the last twenty-five years using Doppler observations, transit photometry, direct imaging, and astrometry. Each of these methods is sensitive to different ranges of orbital separations and planetary radii (or masses). This makes it difficult to fully characterize exoplanet architectures and to place our solar system in context with the wealth of discoveries that have been made. Here, we use the EXtreme PREcision Spectrograph (EXPRES) to reveal planets in previously undetectable regions of the mass-period parameter space for the star $ρ$ Coronae Borealis. We add two new planets to the previously known system with one hot Jupiter in a 39-day orbit and a warm super-Neptune in a 102-day orbit. The new detections include a temperate Neptune planet ($M{\sin{i}} \sim 20$ M$_\oplus$) in a 281.4-day orbit and a hot super-Earth ($M{\sin{i}} = 3.7$ M$_\oplus$) in a 12.95-day orbit. This result shows that details of planetary system architectures have been hiding just below our previous detection limits; this signals an exciting era for the next generation of extreme precision spectrographs. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: Accepted to AJ; 20 pages, 13 figures, 5 Tables

Detection of the Cosmological Time Dilation of High Redshift Quasars

Authors: Geraint F. Lewis, Brendon J. Brewer

Abstract: A fundamental prediction of relativistic cosmologies is that, due to the expansion of space, observations of the distant cosmos should be time dilated and appear to run slower than events in the local universe. Whilst observations of cosmological supernovae unambiguously display the expected redshift-dependent time dilation, this has not been the case for other distant sources. Here we present the… ▽ More A fundamental prediction of relativistic cosmologies is that, due to the expansion of space, observations of the distant cosmos should be time dilated and appear to run slower than events in the local universe. Whilst observations of cosmological supernovae unambiguously display the expected redshift-dependent time dilation, this has not been the case for other distant sources. Here we present the identification of cosmic time dilation in a sample of 190 quasars monitored for over two decades in multiple wavebands by assessing various hypotheses through Bayesian analysis. This detection counters previous claims that observed quasar variability lacked the expected redshift-dependent time dilation. Hence, as well as demonstrating the claim that the lack of the redshift dependence of quasar variability represents a significant challenge to the standard cosmological model, this analysis further indicates that the properties of quasars are consistent with them being truly cosmologically distant sources. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 18 pages, 3 figures. Author version of a manuscript under review at Nature Astronomy

Report number: GFL-001

What Does the Virial Coefficient of the \Hb Broad-Line Region Depend On?

Authors: Lizvette Villafaña, Peter R. Williams, Tommaso Treu, Brendon J. Brewer, Aaron J. Barth, Vivian U, Vardha N. Bennert, Hengxiao Guo, Misty C. Bentz, Gabriela Canalizo, Alexei V. Filippenko, Elinor Gates, Michael D. Joner, Matthew A. Malkan, Jong-Hak Woo, Bela Abolfathi, Thomas Bohn, K. Azalee Bostroem, Andrew Brandel, Thomas G. Brink, Sanyum Channa, Maren Cosens, Edward Donohue, Goni Halevi, Carol E. Hood , et al. (15 additional authors not shown)

Abstract: We combine our dynamical modeling black hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors, f, used in traditional reverberation mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have… ▽ More We combine our dynamical modeling black hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors, f, used in traditional reverberation mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have utilized our methods. Aiming to improve the precision of black hole mass estimates, as well as uncover any regularities in the behavior of the broad-line region (BLR), we search for correlations between f and other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficient log10(fmean,σ) and black hole mass, (ii) marginal evidence for a similar correlation between log10(frms,σ) and black hole mass, (iii) marginal evidence for an anti-correlation of BLR disk thickness with log10(fmean,FWHM)and log10(frms,FWHM), and (iv) marginal evidence for an anti-correlation of inclination angle with log10(fmean,FWHM), log10(frms,σ), and log10(fmean,σ). Lastly, we find marginal evidence for a correlation between line-profile shape, when using the root-meansquare spectrum, log10(FWHM/σ)rms, and the virial coefficient, log10(frms,σ), and investigate how BLR properties might be related to line-profile shape using caramel models. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: 17 pages, 10 figures

Testing the Cosmological Principle with CatWISE Quasars: A Bayesian Analysis of the Number-Count Dipole

Authors: Lawrence Dam, Geraint F. Lewis, Brendon J. Brewer

Abstract: The Cosmological Principle, that the Universe is homogeneous and isotropic on sufficiently large scales, underpins the standard model of cosmology. However, a recent analysis of 1.36 million infrared-selected quasars has identified a significant tension in the amplitude of the number-count dipole compared to that derived from the CMB, thus challenging the Cosmological Principle. Here we present a… ▽ More The Cosmological Principle, that the Universe is homogeneous and isotropic on sufficiently large scales, underpins the standard model of cosmology. However, a recent analysis of 1.36 million infrared-selected quasars has identified a significant tension in the amplitude of the number-count dipole compared to that derived from the CMB, thus challenging the Cosmological Principle. Here we present a Bayesian analysis of the same quasar sample, testing various hypotheses using the Bayesian evidence. We find unambiguous evidence for the presence of a dipole in the distribution of quasars with a direction that is consistent with the dipole identified in the CMB. However, the amplitude of the dipole is found to be 2.7 times larger than that expected from the conventional kinematic explanation of the CMB dipole, with a statistical significance of $5.7σ$. To compare these results with theoretical expectations, we sharpen the $Λ$CDM predictions for the probability distribution of the amplitude, taking into account a number of observational and theoretical systematics. In particular, we show that the presence of the Galactic plane mask causes a considerable loss of dipole signal due to a leakage of power into higher multipoles, exacerbating the discrepancy in the amplitude. By contrast, we show using probabilistic arguments that the source evolution of quasars improves the discrepancy, but only mildly so. These results support the original findings of an anomalously large quasar dipole, independent of the statistical methodology used. △ Less

Submitted 7 August, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 15 pages, 7 figures; accepted for publication in MNRAS

Measured Spin-Orbit Alignment of Ultra-Short Period Super-Earth 55 Cancri e

Authors: Lily L. Zhao, Vedad Kunovac, John M. Brewer, Joe Llama, Sarah C. Millholland, Christina Hedges, Andrew E. Szymkowiak, Rachael M. Roettenbacher, Samuel H. C. Cabot, Sam A. Weiss, Debra A. Fischer

Abstract: A planet's orbital alignment places important constraints on how a planet formed and consequently evolved. The dominant formation pathway of ultra-short period planets ($P<1$ day) is particularly mysterious as such planets most likely formed further out, and it is not well understood what drove their migration inwards to their current positions. Measuring the orbital alignment is difficult for sma… ▽ More A planet's orbital alignment places important constraints on how a planet formed and consequently evolved. The dominant formation pathway of ultra-short period planets ($P<1$ day) is particularly mysterious as such planets most likely formed further out, and it is not well understood what drove their migration inwards to their current positions. Measuring the orbital alignment is difficult for smaller super-Earth/sub-Neptune planets, which give rise to smaller amplitude signals. Here we present radial velocities across two transits of 55 Cancri e, an ultra-short period Super-Earth, observed with the Extreme Precision Spectrograph (EXPRES). Using the classical Rossiter-McLaughlin (RM) method, we measure 55 Cnc e's sky-projected stellar spin-orbit alignment (i.e., the projected angle between the planet's orbital axis and its host star's spin axis) to be $λ=10\substack{+17\\ -20}^{\circ}$ with an unprojected angle of $ψ=23\substack{+14\\ -12}^{\circ}$. The best-fit RM model to the EXPRES data has a radial velocity semi-amplitude of just $0.41\substack{+0.09\\ -0.10} m s^{-1}$. The spin-orbit alignment of 55 Cnc e favors dynamically gentle migration theories for ultra-short period planets, namely tidal dissipation through low-eccentricity planet-planet interactions and/or planetary obliquity tides. △ Less

Submitted 9 December, 2022; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 12 pages, 4 figures, published in Nature Astronomy

Chemo-dynamical substructure in the M31 inner halo globular clusters: Further evidence for a recent accretion event

Authors: Geraint F. Lewis, Brendon J. Brewer, Dougal Mackey, Annette M. N. Ferguson, Yuan, Li, Tim Adams

Abstract: Based upon a metallicity selection, we identify a significant sub-population of the inner halo globular clusters in the Andromeda Galaxy which we name the Dulais Structure. It is distinguished as a co-rotating group of 10-20 globular clusters which appear to be kinematically distinct from, and on average more metal-poor than, the majority of the inner halo population. Intriguingly, the orbital axi… ▽ More Based upon a metallicity selection, we identify a significant sub-population of the inner halo globular clusters in the Andromeda Galaxy which we name the Dulais Structure. It is distinguished as a co-rotating group of 10-20 globular clusters which appear to be kinematically distinct from, and on average more metal-poor than, the majority of the inner halo population. Intriguingly, the orbital axis of this Dulais Structure is closely aligned with that of the younger accretion event recently identified using a sub-population of globular clusters in the outer halo of Andromeda, and this is strongly suggestive of a causal relationship between the two. If this connection is confirmed, a natural explanation for the kinematics of the globular clusters in the Dulais Structure is that they trace the accretion of a substantial progenitor (~10^11 Msun) into the halo of Andromeda during the last few billion years, that may have occurred as part of a larger group infall. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 10 pages, 7 figures - accepted for publication in MNRAS

Report number: GFL-001

Planet Engulfment Detections are Rare According to Observations and Stellar Modeling

Authors: Aida Behmard, Fei Dai, John M. Brewer, Travis A. Berger, Andrew W. Howard

Abstract: Dynamical evolution within planetary systems can cause planets to be engulfed by their host stars. Following engulfment, the stellar photosphere abundance pattern will reflect accretion of rocky material from planets. Multi-star systems are excellent environments to search for such abundance trends because stellar companions form from the same natal gas cloud and are thus expected to share primord… ▽ More Dynamical evolution within planetary systems can cause planets to be engulfed by their host stars. Following engulfment, the stellar photosphere abundance pattern will reflect accretion of rocky material from planets. Multi-star systems are excellent environments to search for such abundance trends because stellar companions form from the same natal gas cloud and are thus expected to share primordial chemical compositions to within 0.03$-$0.05 dex. Abundance measurements have occasionally yielded rocky enhancements, but few observations targeted known planetary systems. To address this gap, we carried out a Keck-HIRES survey of 36 multi-star systems where at least one star is a known planet host. We found that only HAT-P-4 exhibits an abundance pattern suggestive of engulfment, but is more likely primordial based on its large projected separation (30,000 $\pm$ 140 AU) that exceeds typical turbulence scales in molecular clouds. To understand the lack of engulfment detections among our systems, we quantified the strength and duration of refractory enrichments in stellar photospheres using MESA stellar models. We found that observable signatures from 10 $M_{\oplus}$ engulfment events last for $\sim$90 Myr in 1 $M_{\odot}$ stars. Signatures are largest and longest lived for 1.1$-$1.2 $M_{\odot}$ stars, but are no longer observable $\sim$2 Gyr post-engulfment. This indicates that engulfment will rarely be detected in systems that are several Gyr old. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 15 pages, 12 figures; submitted to MNRAS

The Lick AGN Monitoring Project 2016: Dynamical Modeling of Velocity-Resolved H\b{eta} Lags in Luminous Seyfert Galaxies

Authors: Lizvette Villafaña, Peter R. Williams, Tommaso Treu, Brendon J. Brewer, Aaron J. Barth, Vivian U, Vardha N. Bennert, H. Alexander Vogler, Hengxiao Guo, Misty C. Bentz, Gabriela Canalizo, Alexei V. Filippenko, Elinor Gates, Frederick Hamann, Michael D. Joner, Matthew A. Malkan, Jong-Hak Woo, Bela Abolfathi, L. E. Abramson, Stephen F. Armen, Hyun-Jin Bae, Thomas Bohn, Benjamin D. Boizelle, Azalee Bostroem, Andrew Brandel , et al. (40 additional authors not shown)

Abstract: We have modeled the velocity-resolved reverberation response of the H\b{eta} broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitioring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the H\b{eta} BLR is generally… ▽ More We have modeled the velocity-resolved reverberation response of the H\b{eta} broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitioring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the H\b{eta} BLR is generally a thick disk viewed at low to moderate inclination angles. We combine our sample with prior studies and investigate line-profile shape dependence, such as log10(FWHM/σ), on BLR structure and kinematics and search for any BLR luminosity-dependent trends. We find marginal evidence for an anticorrelation between the profile shape of the broad H\b{eta} emission line and the Eddington ratio, when using the root-mean-square spectrum. However, we do not find any luminosity-dependent trends, and conclude that AGNs have diverse BLR structure and kinematics, consistent with the hypothesis of transient AGN/BLR conditions rather than systematic trends. △ Less

Submitted 28 March, 2022; originally announced March 2022.

The EXPRES Stellar Signals Project II. State of the Field in Disentangling Photospheric Velocities

Authors: Lily L. Zhao, Debra A. Fischer, Eric B. Ford, Alex Wise, Michaël Cretignier, Suzanne Aigrain, Oscar Barragan, Megan Bedell, Lars A. Buchhave, João D. Camacho, Heather M. Cegla, Jessi Cisewski-Kehe, Andrew Collier Cameron, Zoe L. de Beurs, Sally Dodson-Robinson, Xavier Dumusque, João P. Faria, Christian Gilbertson, Charlotte Haley, Justin Harrell, David W. Hogg, Parker Holzer, Ancy Anna John, Baptiste Klein, Marina Lafarga , et al. (18 additional authors not shown)

Abstract: Measured spectral shifts due to intrinsic stellar variability (e.g., pulsations, granulation) and activity (e.g., spots, plages) are the largest source of error for extreme precision radial velocity (EPRV) exoplanet detection. Several methods are designed to disentangle stellar signals from true center-of-mass shifts due to planets. The EXPRES Stellar Signals Project (ESSP) presents a self-consist… ▽ More Measured spectral shifts due to intrinsic stellar variability (e.g., pulsations, granulation) and activity (e.g., spots, plages) are the largest source of error for extreme precision radial velocity (EPRV) exoplanet detection. Several methods are designed to disentangle stellar signals from true center-of-mass shifts due to planets. The EXPRES Stellar Signals Project (ESSP) presents a self-consistent comparison of 22 different methods tested on the same extreme-precision spectroscopic data from EXPRES. Methods derived new activity indicators, constructed models for mapping an indicator to the needed RV correction, or separated out shape- and shift-driven RV components. Since no ground truth is known when using real data, relative method performance is assessed using the total and nightly scatter of returned RVs and agreement between the results of different methods. Nearly all submitted methods return a lower RV RMS than classic linear decorrelation, but no method is yet consistently reducing the RV RMS to sub-meter-per-second levels. There is a concerning lack of agreement between the RVs returned by different methods. These results suggest that continued progress in this field necessitates increased interpretability of methods, high-cadence data to capture stellar signals at all timescales, and continued tests like the ESSP using consistent data sets with more advanced metrics for method performance. Future comparisons should make use of various well-characterized data sets -- such as solar data or data with known injected planetary and/or stellar signals -- to better understand method performance and whether planetary signals are preserved. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 33 pages (+12 pages of Appendix), 10 figures, 8 tables, accepted for publication in AJ

EXPRES. III. Revealing the Stellar Activity Radial Velocity Signature of $ε$ Eridani with Photometry and Interferometry

Authors: Rachael M. Roettenbacher, Samuel H. C. Cabot, Debra A. Fischer, John D. Monnier, Gregory W. Henry, Robert O. Harmon, Heidi Korhonen, John M. Brewer, Joe Llama, Ryan R. Petersburg, Lily Zhao, Stefan Kraus, Jean-Baptiste Le Bouquin, Narsireddy Anugu, Claire L. Davies, Tyler Gardner, Cyprien Lanthermann, Gail Schaefer, Benjamin Setterholm, Catherine A. Clark, Svetlana G. Jorstad, Kyler Kuehn, Stephen Levine

Abstract: The distortions of absorption line profiles caused by photospheric brightness variations on the surfaces of cool, main-sequence stars can mimic or overwhelm radial velocity (RV) shifts due to the presence of exoplanets. The latest generation of precision RV spectrographs aims to detect velocity amplitudes $\lesssim 10$ cm s$^{-1}$, but requires mitigation of stellar signals. Statistical techniques… ▽ More The distortions of absorption line profiles caused by photospheric brightness variations on the surfaces of cool, main-sequence stars can mimic or overwhelm radial velocity (RV) shifts due to the presence of exoplanets. The latest generation of precision RV spectrographs aims to detect velocity amplitudes $\lesssim 10$ cm s$^{-1}$, but requires mitigation of stellar signals. Statistical techniques are being developed to differentiate between Keplerian and activity-related velocity perturbations. Two important challenges, however, are the interpretability of the stellar activity component as RV models become more sophisticated, and ensuring the lowest-amplitude Keplerian signatures are not inadvertently accounted for in flexible models of stellar activity. For the K2V exoplanet host $ε$ Eridani, we separately use ground-based photometry to constrain Gaussian processes for modeling RVs and TESS photometry with a light-curve inversion algorithm to reconstruct the stellar surface. From the reconstructions of TESS photometry, we produce an activity model, which reduces the rms scatter in RVs obtained with EXPRES from 4.72 m s$^{-1}$ to 1.98 m s$^{-1}$. We present a pilot study using the CHARA Array and MIRC-X beam combiner to directly image the starspots seen in the TESS photometry. With the limited phase coverage, our spot detections are marginal with current data but a future dedicated observing campaign should allow for imaging, as well as the stellar inclination and orientation with respect to its debris disk to be definitely determined. This work shows that stellar surface maps obtained with high cadence, time-series photometric and interferometric data can provide the constraints needed to accurately reduce RV scatter. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 24 pages, 12 figures, 5 tables, accepted for publication in AJ

Multi-scale photonic emissivity engineering for relativistic lightsail thermal regulation

Authors: John Brewer, Matthew F. Campbell, Pawan Kumar, Sachin Kulkarni, Deep Jariwala, Igor Bargatin, Aaswath P. Raman

Abstract: The Breakthrough Starshot Initiative aims to send a gram-scale probe to Proxima Centuri B using a laser-accelerated lightsail traveling at relativistic speeds. Thermal management is a key lightsail design objective because of the intense laser powers required but has generally been considered secondary to accelerative performance. Here, we demonstrate nanophotonic photonic crystal slab reflectors… ▽ More The Breakthrough Starshot Initiative aims to send a gram-scale probe to Proxima Centuri B using a laser-accelerated lightsail traveling at relativistic speeds. Thermal management is a key lightsail design objective because of the intense laser powers required but has generally been considered secondary to accelerative performance. Here, we demonstrate nanophotonic photonic crystal slab reflectors composed of 2H-phase molybdenum disulfide and crystalline silicon nitride, highlight the inverse relationship between the thermal band extinction coefficient and the lightsail's maximum temperature, and examine the trade-off between the acceleration distance and setting realistic sail thermal limits, ultimately realizing a thermally endurable acceleration minimum distance of 16.3~Gm. We additionally demonstrate multi-scale photonic structures featuring thermal-wavelength-scale Mie resonant geometries, and characterize their broadband Mie resonance-driven emissivity enhancement and acceleration distance reduction. Our results highlight new possibilities in simultaneously controlling optical and thermal response over broad wavelength ranges in ultralight nanophotonic structures. △ Less

Submitted 14 September, 2021; v1 submitted 3 June, 2021; originally announced June 2021.

Comments: 15 pages, 4 figures

Relativistic light sails need to billow

Authors: Matthew F. Campbell, John Brewer, Deep Jariwala, Aaswath Raman, Igor Bargatin

Abstract: We argue that light sails that are rapidly accelerated to relativistic velocities by lasers must be significantly curved in order to reduce their mechanical stresses and avoid tears. Using an integrated opto-thermo-mechanical model, we show that the diameter and radius of curvature of a circular light sail should be comparable in magnitude, both on the order of a few meters in optimal designs for… ▽ More We argue that light sails that are rapidly accelerated to relativistic velocities by lasers must be significantly curved in order to reduce their mechanical stresses and avoid tears. Using an integrated opto-thermo-mechanical model, we show that the diameter and radius of curvature of a circular light sail should be comparable in magnitude, both on the order of a few meters in optimal designs for gram-scale payloads. Moreover, when sufficient laser power is available, a sail's acceleration length decreases and its chip payload capacity increases as its curvature increases. Our findings provide guidance for emerging light sail design programs, which herald a new era of interstellar space exploration. △ Less

Submitted 22 May, 2021; originally announced May 2021.

Comments: 8 pages, 4 figures

Dynamical Modeling of the CIV Broad Line Region of the $z=2.805$ Multiply Imaged Quasar SDSS J2222+2745

Authors: Peter R. Williams, Tommaso Treu, Håkon Dahle, Stefano Valenti, Louis Abramson, Aaron J. Barth, Brendon J. Brewer, Karianne Dyrland, Michael Gladders, Keith Horne, Keren Sharon

Abstract: We present CIV BLR modeling results for the multiply imaged $z=2.805$ quasar SDSS J2222+2745. Using data covering a 5.3 year baseline after accounting for gravitational time delays, we find models that can reproduce the observed emission-line spectra and integrated CIV fluctuations. The models suggest a thick disk BLR that is inclined by $\sim$40 degrees to the observer's line of sight and with a… ▽ More We present CIV BLR modeling results for the multiply imaged $z=2.805$ quasar SDSS J2222+2745. Using data covering a 5.3 year baseline after accounting for gravitational time delays, we find models that can reproduce the observed emission-line spectra and integrated CIV fluctuations. The models suggest a thick disk BLR that is inclined by $\sim$40 degrees to the observer's line of sight and with a emissivity weighted median radius of $r_{\rm median} = 33.0^{+2.4}_{-2.1}$ light days. The kinematics are dominated by near-circular Keplerian motion with the remainder inflowing. The rest-frame lag one would measure from the models is $τ_{\rm median} = 36.4^{+1.8}_{-1.8}$ days, which is consistent with measurements based on cross-correlation. We show a possible geometry and transfer function based on the model fits and find that the model-produced velocity-resolved lags are consistent with those from cross-correlation. We measure a black hole mass of $\log_{10}(M_{\rm BH}/M_\odot) = 8.31^{+0.07}_{-0.06}$, which requires a scale factor of $\log_{10}(f_{{\rm mean},σ}) = 0.20^{+0.09}_{-0.07}$. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: 8 pages, 5 figures, submitted to ApJL

Opportunities of OO and $p$O collisions at the LHC

Authors: Jasmine Brewer, Aleksas Mazeliauskas, Wilke van der Schee

Abstract: This is the summary document of the virtual workshop "Opportunities of OO and $p$O collisions at the LHC", which took place Feb 4-10, 2021. We briefly review the presented perspectives on the physics opportunities of oxygen--oxygen and proton--oxygen collisions. The full workshop program and the recordings are available at cern.ch/OppOatLHC. This is the summary document of the virtual workshop "Opportunities of OO and $p$O collisions at the LHC", which took place Feb 4-10, 2021. We briefly review the presented perspectives on the physics opportunities of oxygen--oxygen and proton--oxygen collisions. The full workshop program and the recordings are available at cern.ch/OppOatLHC. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: Workshop summary (see http://cern.ch/OppOatLHC ), 6 pages, 9 figures

Report number: CERN-TH-2021-028

Transmission spectroscopy for the warm sub-Neptune HD3167c: evidence for molecular absorption and a possible high metallicity atmosphere

Authors: Thomas Mikal-Evans, Ian J. M. Crossfield, Bjorn Benneke, Laura Kreidberg, Julie Moses, Caroline V. Morley, Daniel Thorngren, Paul Molliere, Kevin K. Hardegree-Ullman, John Brewer, Jessie L. Christiansen, David R. Ciardi, Diana Dragomir, Courtney Dressing, Jonathan J. Fortney, Varoujan Gorjian, Thomas P. Greene, Lea A. Hirsch, Andrew W. Howard, Steve B. Howell, Howard Isaacson, Molly R. Kosiarek, Jessica Krick, John H. Livingston, Joshua D. Lothringer , et al. (4 additional authors not shown)

Abstract: We present a transmission spectrum for the warm (500-600K) sub-Neptune HD3167c obtained using the Hubble Space Telescope Wide Field Camera 3 infrared spectrograph. We combine these data, which span the 1.125-1.643 micron wavelength range, with broadband transit measurements made using Kepler/K2 (0.6-0.9 micron) and Spitzer/IRAC (4-5 micron). We find evidence for absorption by at least one of H2O,… ▽ More We present a transmission spectrum for the warm (500-600K) sub-Neptune HD3167c obtained using the Hubble Space Telescope Wide Field Camera 3 infrared spectrograph. We combine these data, which span the 1.125-1.643 micron wavelength range, with broadband transit measurements made using Kepler/K2 (0.6-0.9 micron) and Spitzer/IRAC (4-5 micron). We find evidence for absorption by at least one of H2O, HCN, CO2, and CH4 (Bayes factor 7.4; 2.5-sigma significance), although the data precision does not allow us to unambiguously discriminate between these molecules. The transmission spectrum rules out cloud-free hydrogen-dominated atmospheres with metallicities <100x solar at >5.8-sigma confidence. In contrast, good agreement with the data is obtained for cloud-free models assuming metallicities >700x solar. However, for retrieval analyses that include the effect of clouds, a much broader range of metallicities (including subsolar) is consistent with the data, due to the degeneracy with cloud-top pressure. Self-consistent chemistry models that account for photochemistry and vertical mixing are presented for the atmosphere of HD3167c. The predictions of these models are broadly consistent with our abundance constraints, although this is primarily due to the large uncertainties on the latter. Interior structure models suggest the core mass fraction is >40%, independent of a rock or water core composition, and independent of atmospheric envelope metallicity up to 1000x solar. We also report abundance measurements for fifteen elements in the host star, showing that it has a very nearly solar composition. △ Less

Submitted 6 November, 2020; originally announced November 2020.

Comments: Accepted for publication in the Astronomical Journal

EXPRES. II. Searching for Planets Around Active Stars: A Case Study of HD 101501

Authors: Samuel H. C. Cabot, Rachael M. Roettenbacher, Gregory W. Henry, Lily Zhao, Robert O. Harmon, Debra A. Fischer, John M. Brewer, Joe Llama, Ryan R. Petersburg, Andrew E. Szymkowiak

Abstract: By controlling instrumental errors to below 10 cm/s, the EXtreme PREcision Spectrograph (EXPRES) allows for a more insightful study of photospheric velocities that can mask weak Keplerian signals. Gaussian Processes (GP) have become a standard tool for modeling correlated noise in radial velocity datasets. While GPs are constrained and motivated by physical properties of the star, in some cases th… ▽ More By controlling instrumental errors to below 10 cm/s, the EXtreme PREcision Spectrograph (EXPRES) allows for a more insightful study of photospheric velocities that can mask weak Keplerian signals. Gaussian Processes (GP) have become a standard tool for modeling correlated noise in radial velocity datasets. While GPs are constrained and motivated by physical properties of the star, in some cases they are still flexible enough to absorb unresolved Keplerian signals. We apply GP regression to EXPRES radial velocity measurements of the 3.5 Gyr old chromospherically active Sun-like star, HD 101501. We obtain tight constraints on the stellar rotation period and the evolution of spot distributions using 28 seasons of ground-based photometry, as well as recent $TESS$ data. Light curve inversion was carried out on both photometry datasets to reveal the spot distribution and spot evolution timescales on the star. We find that the $> 5$ m/s rms radial velocity variations in HD 101501 are well-modeled with a GP stellar activity model without planets, yielding a residual rms scatter of 45 cm/s. We carry out simulations, injecting and recovering signals with the GP framework, to demonstrate that high-cadence observations are required to use GPs most efficiently to detect low-mass planets around active stars like HD 101501. Sparse sampling prevents GPs from learning the correlated noise structure and can allow it to absorb prospective Keplerian signals. We quantify the moderate to high-cadence monitoring that provides the necessary information to disentangle photospheric features using GPs and to detect planets around active stars. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: 25 pages, 16 figures, accepted to AJ

Space Telescope and Optical Reverberation Mapping Project. XII. Broad-Line Region Modeling of NGC 5548

Authors: P. R. Williams, A. Pancoast, T. Treu, B. J. Brewer, B. M. Peterson, A. J. Barth, M. A. Malkan, G. De Rosa, Keith Horne, G. A. Kriss, N. Arav, M. C. Bentz, E. M. Cackett, E. Dalla Bontà, M. Dehghanian, C. Done, G. J. Ferland, C. J. Grier, J. Kaastra, E. Kara, C. S. Kochanek, S. Mathur, M. Mehdipour, R. W. Pogge, D. Proga , et al. (133 additional authors not shown)

Abstract: We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$β$, C IV, and Ly$α$ broad emission lines. We find an extended disk-like H$β$ BLR with a mixture of near-circular and outflowing gas traje… ▽ More We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$β$, C IV, and Ly$α$ broad emission lines. We find an extended disk-like H$β$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$α$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$α$ emission arising at smaller radii than the H$β$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$β$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

The TESS-Keck Survey II: An Ultra-Short Period Rocky Planet and its Siblings Transiting the Galactic Thick-Disk Star TOI-561

Authors: Lauren M. Weiss, Fei Dai, Daniel Huber, John M. Brewer, Karen A. Collins, David R. Ciardi, Elisabeth C. Matthews, Carl Ziegler, Steve B. Howell, Natalie M. Batalha, Ian J. M> Crossfield, Courtney Dressing, Benjamin Fulton, Andrew W. Howard, Howard Isaacson, Stephen R. Kane, Erik A. Petigura, Paul Robertson, Arpita Roy, Ryan A. Rubenzahl, Joseph D. Twicken, Zachary R. Claytor, Keivan G. Stassun, Mason G. MacDougall, Ashley Chontos , et al. (39 additional authors not shown)

Abstract: We report the discovery of TOI-561, a multi-planet system in the galactic thick disk that contains a rocky, ultra-short period planet (USP). This bright ($V=10.2$) star hosts three small transiting planets identified in photometry from the NASA TESS mission: TOI-561 b (TOI-561.02, P=0.44 days, $R_b = 1.45\pm0.11\,R_\oplus$), c (TOI-561.01, P=10.8 days, $R_c=2.90\pm0.13\,R_\oplus$), and d (TOI-561.… ▽ More We report the discovery of TOI-561, a multi-planet system in the galactic thick disk that contains a rocky, ultra-short period planet (USP). This bright ($V=10.2$) star hosts three small transiting planets identified in photometry from the NASA TESS mission: TOI-561 b (TOI-561.02, P=0.44 days, $R_b = 1.45\pm0.11\,R_\oplus$), c (TOI-561.01, P=10.8 days, $R_c=2.90\pm0.13\,R_\oplus$), and d (TOI-561.03, P=16.3 days, $R_d=2.32\pm0.16\,R_\oplus$). The star is chemically ([Fe/H]$=-0.41\pm0.05$, [$α$/H]$=+0.23\pm0.05$) and kinematically consistent with the galactic thick disk population, making TOI-561 one of the oldest ($10\pm3\,$Gyr) and most metal-poor planetary systems discovered yet. We dynamically confirm planets b and c with radial velocities from the W. M. Keck Observatory High Resolution Echelle Spectrometer. Planet b has a mass and density of $3.2\pm0.8\,M_\oplus$ and $5.5^{+2.0}_{-1.6}\,$g$\,$cm$^{-3}$, consistent with a rocky composition. Its lower-than-average density is consistent with an iron-poor composition, although an Earth-like iron-to-silicates ratio is not ruled out. Planet c is $7.0\pm2.3\,M_\oplus$ and $1.6\pm0.6\,$g$\,$cm$^{-3}$, consistent with an interior rocky core overlaid with a low-mass volatile envelope. Several attributes of the photometry for planet d (which we did not detect dynamically) complicate the analysis, but we vet the planet with high-contrast imaging, ground-based photometric follow-up and radial velocities. TOI-561 b is the first rocky world around a galactic thick-disk star confirmed with radial velocities and one of the best rocky planets for thermal emission studies. △ Less

Submitted 14 December, 2020; v1 submitted 7 September, 2020; originally announced September 2020.

Comments: Accepted at The Astronomical Journal; 25 pages, 10 figures

Stellar Characterization of Keck HIRES Spectra with The Cannon

Authors: Malena Rice, John Brewer

Abstract: To accurately interpret the observed properties of exoplanets, it is necessary to first obtain a detailed understanding of host star properties. However, physical models that analyze stellar properties on a per-star basis can become computationally intractable for sufficiently large samples. Furthermore, these models are limited by the wavelength coverage of available spectra. We combine previousl… ▽ More To accurately interpret the observed properties of exoplanets, it is necessary to first obtain a detailed understanding of host star properties. However, physical models that analyze stellar properties on a per-star basis can become computationally intractable for sufficiently large samples. Furthermore, these models are limited by the wavelength coverage of available spectra. We combine previously derived spectral properties from the Spectroscopic Properties of Cool Stars (SPOCS) catalog (Brewer et al. 2016) with generative modeling using The Cannon to produce a model capable of deriving stellar parameters ($\log g$, $T_{\mathrm{eff}}$, and $v\sin i$) and 15 elemental abundances (C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, and Y) for stellar spectra observed with Keck Observatory's High Resolution Echelle Spectrometer (HIRES). We demonstrate the high accuracy and precision of our model, which takes just $\sim$3 seconds to classify each star, through cross-validation with pre-labeled spectra from the SPOCS sample. Our trained model, which takes continuum-normalized template spectra as its inputs, is publicly available at https://github.com/malenarice/keckspec. Finally, we interpolate our spectra and employ the same modeling scheme to recover labels for 477 stars using archival stellar spectra obtained prior to Keck's 2004 detector upgrade, demonstrating that our interpolated model can successfully predict stellar labels for different spectrographs that have (1) sufficiently similar systematics and (2) a wavelength range that substantially overlaps with that of the post-2004 HIRES spectra. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: Accepted for publication in ApJ. 27 pages, 16 figures

EXPRES I. HD~3651 an Ideal RV Benchmark

Authors: John M. Brewer, Debra A. Fischer, Ryan T. Blackman, Samuel H. C. Cabot, Allen B. Davis, Gregory Laughlin, Christopher Leet, J. M. Joel Ong, Ryan R. Petersburg, Andrew E. Szymkowiak, Lily L. Zhao, Gregory W. Henry, Joe Llama

Abstract: The next generation of exoplanet-hunting spectrographs should deliver up to an order of magnitude improvement in radial velocity precision over the standard 1 m/s state of the art. This advance is critical for enabling the detection of Earth-mass planets around Sun-like stars. New calibration techniques such as laser frequency combs and stabilized etalons ensure that the instrumental stability is… ▽ More The next generation of exoplanet-hunting spectrographs should deliver up to an order of magnitude improvement in radial velocity precision over the standard 1 m/s state of the art. This advance is critical for enabling the detection of Earth-mass planets around Sun-like stars. New calibration techniques such as laser frequency combs and stabilized etalons ensure that the instrumental stability is well characterized. However, additional sources of error include stellar noise, undetected short-period planets, and telluric contamination. To understand and ultimately mitigate error sources, the contributing terms in the error budget must be isolated to the greatest extent possible. Here, we introduce a new high cadence radial velocity program, the EXPRES 100 Earths program, which aims to identify rocky planets around bright, nearby G and K dwarfs. We also present a benchmark case: the 62-d orbit of a Saturn-mass planet orbiting the chromospherically quiet star, HD 3651. The combination of high eccentricity (0.6) and a moderately long orbital period, ensures significant dynamical clearing of any inner planets. Our Keplerian model for this planetary orbit has a residual RMS of 58 cm/s over a $\sim 6$ month time baseline. By eliminating significant contributors to the radial velocity error budget, HD 3651 serves as a standard for evaluating the long term precision of extreme precision radial velocity (EPRV) programs. △ Less

Submitted 3 June, 2020; originally announced June 2020.

Comments: 11 pages, 6 figures, accepted for publication in Astronomical Journal

Performance Verification of the EXtreme PREcision Spectrograph

Authors: Ryan T. Blackman, Debra A. Fischer, Colby A. Jurgenson, David Sawyer, Tyler M. McCracken, Andrew E. Szymkowiak, Ryan R. Petersburg, J. M. Joel Ong, John M. Brewer, Lily L. Zhao, Christopher Leet, Lars A. Buchhave, René Tronsgaard, Joe Llama, Travis Sawyer, Allen B. Davis, Samuel H. C. Cabot, Michael Shao, Russell Trahan, Bijan Nemati, Matteo Genoni, Giorgio Pariani, Marco Riva, Rafael A. Probst, Ronald Holzwarth , et al. (3 additional authors not shown)

Abstract: The EXtreme PREcision Spectrograph (EXPRES) is a new Doppler spectrograph designed to reach a radial velocity measurement precision sufficient to detect Earth-like exoplanets orbiting nearby, bright stars. We report on extensive laboratory testing and on-sky observations to quantitatively assess the instrumental radial velocity measurement precision of EXPRES, with a focused discussion of individu… ▽ More The EXtreme PREcision Spectrograph (EXPRES) is a new Doppler spectrograph designed to reach a radial velocity measurement precision sufficient to detect Earth-like exoplanets orbiting nearby, bright stars. We report on extensive laboratory testing and on-sky observations to quantitatively assess the instrumental radial velocity measurement precision of EXPRES, with a focused discussion of individual terms in the instrument error budget. We find that EXPRES can reach a single-measurement instrument calibration precision better than 10 cm/s, not including photon noise from stellar observations. We also report on the performance of the various environmental, mechanical, and optical subsystems of EXPRES, assessing any contributions to radial velocity error. For atmospheric and telescope related effects, this includes the fast tip-tilt guiding system, atmospheric dispersion compensation, and the chromatic exposure meter. For instrument calibration, this includes the laser frequency comb (LFC), flat-field light source, CCD detector, and effects in the optical fibers. Modal noise is mitigated to a negligible level via a chaotic fiber agitator, which is especially important for wavelength calibration with the LFC. Regarding detector effects, we empirically assess the impact on radial velocity precision due to pixel-position non-uniformities (PPNU) and charge transfer inefficiency (CTI). EXPRES has begun its science survey to discover exoplanets orbiting G-dwarf and K-dwarf stars, in addition to transit spectroscopy and measurements of the Rossiter-McLaughlin effect. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: 39 pages, 30 figures, accepted to AJ

An Extreme Precision Radial Velocity Pipeline: First Radial Velocities from EXPRES

Authors: Ryan R. Petersburg, J. M. Joel Ong, Lily L. Zhao, Ryan T. Blackman, John M. Brewer, Lars A. Buchhave, Samuel H. C. Cabot, Allen B. Davis, Colby A. Jurgenson, Christopher Leet, Tyler M. McCracken, David Sawyer, Mikhail Sharov, René Tronsgaard, Andrew E. Szymkowiak, Debra A. Fischer

Abstract: The EXtreme PREcision Spectrograph (EXPRES) is an environmentally stabilized, fiber-fed, $R=137,500$, optical spectrograph. It was recently commissioned at the 4.3-m Lowell Discovery Telescope (LDT) near Flagstaff, Arizona. The spectrograph was designed with a target radial-velocity (RV) precision of 30$\mathrm{~cm~s^{-1}}$. In addition to instrumental innovations, the EXPRES pipeline, presented h… ▽ More The EXtreme PREcision Spectrograph (EXPRES) is an environmentally stabilized, fiber-fed, $R=137,500$, optical spectrograph. It was recently commissioned at the 4.3-m Lowell Discovery Telescope (LDT) near Flagstaff, Arizona. The spectrograph was designed with a target radial-velocity (RV) precision of 30$\mathrm{~cm~s^{-1}}$. In addition to instrumental innovations, the EXPRES pipeline, presented here, is the first for an on-sky, optical, fiber-fed spectrograph to employ many novel techniques---including an "extended flat" fiber used for wavelength-dependent quantum efficiency characterization of the CCD, a flat-relative optimal extraction algorithm, chromatic barycentric corrections, chromatic calibration offsets, and an ultra-precise laser frequency comb for wavelength calibration. We describe the reduction, calibration, and radial-velocity analysis pipeline used for EXPRES and present an example of our current sub-meter-per-second RV measurement precision, which reaches a formal, single-measurement error of 0.3$\mathrm{~m~s^{-1}}$ for an observation with a per-pixel signal-to-noise ratio of 250. These velocities yield an orbital solution on the known exoplanet host 51 Peg that matches literature values with a residual RMS of 0.895$\mathrm{~m~s^{-1}}$. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Do Metal-Rich Stars Make Metal-Rich Planets? New Insights on Giant Planet Formation from Host Star Abundances

Authors: Johanna K. Teske, Daniel Thorngren, Jonathan J. Fortney, Natalie Hinkel, John M. Brewer

Abstract: The relationship between the compositions of giant planets and their host stars is of fundamental interest in understanding planet formation. The solar system giant planets are enhanced above solar composition in metals, both in their visible atmospheres and bulk compositions. A key question is whether the metal enrichment of giant exoplanets is correlated with that of their host stars. Thorngren… ▽ More The relationship between the compositions of giant planets and their host stars is of fundamental interest in understanding planet formation. The solar system giant planets are enhanced above solar composition in metals, both in their visible atmospheres and bulk compositions. A key question is whether the metal enrichment of giant exoplanets is correlated with that of their host stars. Thorngren et al. (2016) showed that in cool (Teq < 1000 K) giant exoplanets, the total heavy-element mass increases with total Mp and the heavy element enrichment relative to the parent star decreases with total Mp. In their work, the host star metallicity was derived from literature [Fe/H] measurements. Here we conduct a more detailed and uniform study to determine whether different host star metals (C, O, Mg, Si, Fe, and Ni) correlate with the bulk metallicity of their planets, using correlation tests and Bayesian linear fits. We present new host star abundances of 19 cool giant planet systems, and combine these with existing host star data for a total of 22 cool giant planet systems (24 planets). Surprisingly, we find no clear correlation between stellar metallicity and planetary residual metallicity (the relative amount of metal versus that expected from the planet mass alone), which is in conflict with common predictions from formation models. We also find a potential correlation between residual planet metals and stellar volatile-to-refractory element ratios. These results provide intriguing new relationships between giant planet and host star compositions for future modeling studies of planet formation. △ Less

Submitted 30 November, 2019; originally announced December 2019.

Comments: Main results in Figure 3-5. Published in AJ on 2019-11-20. 24 pages (single column style)

Journal ref: The Astronomical Journal, Volume 158, Issue 6, article id. 239, 18 pp. (2019) Link: https://iopscience.iop.org/article/10.3847/1538-3881/ab4f79

The Globular Cluster Population of NGC 1052-DF2: Evidence for Rotation

Authors: Geraint F. Lewis, Brendon J. Brewer, Zhen Wan

Abstract: Based upon the kinematics of ten globular clusters, it has recently been claimed that the ultra-diffuse galaxy, NCD 1052-DF2, lacks a significant quantity of dark matter. Dynamical analyses have generally assumed that this galaxy is pressure supported, with the relatively small velocity dispersion of the globular cluster population indicating the deficit of dark matter. However, the presence of a… ▽ More Based upon the kinematics of ten globular clusters, it has recently been claimed that the ultra-diffuse galaxy, NCD 1052-DF2, lacks a significant quantity of dark matter. Dynamical analyses have generally assumed that this galaxy is pressure supported, with the relatively small velocity dispersion of the globular cluster population indicating the deficit of dark matter. However, the presence of a significant rotation of the globular cluster population could substantially modify this conclusion. Here we present the discovery of such a signature of rotation in the kinematics of NGC 1052-DF2's globular clusters, with a velocity amplitude of $\sim12.44^{+4.40}_{-5.16}$ km/s, which, through Bayesian model comparison, represents a marginally better fit to the available kinematic data; note that this rotation is distinct from, and approximately perpendicular to, the recently identified rotation of the stellar component of NGC 1052-DF2. Assuming this truly represents an underlying rotation, it is shown that the determined mass depends upon the inclination of the rotational component and, with a moderate inclination, the resultant mass to light ratio can exceed $M/L\sim10$. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS Letters

Report number: GFL-001

Two major accretion epochs in M31 from two distinct populations of globular clusters

Authors: Dougal Mackey, Geraint F. Lewis, Brendon J. Brewer, Annette M. N. Ferguson, Jovan Veljanoski, Avon P. Huxor, Michelle L. M. Collins, Patrick Côté, Rodrigo A. Ibata, Mike J. Irwin, Nicolas Martin, Alan W. McConnachie, Jorge Peñarrubia, Nial Tanvir, Zhen Wan

Abstract: Large galaxies grow through the accumulation of dwarf galaxies. In principle it is possible to trace this growth history using the properties of a galaxy's stellar halo. Previous investigations of the galaxy M31 (Andromeda) have shown that outside a radius of 25 kpc the population of halo globular clusters is rotating in alignment with the stellar disk, as are more centrally located clusters. The… ▽ More Large galaxies grow through the accumulation of dwarf galaxies. In principle it is possible to trace this growth history using the properties of a galaxy's stellar halo. Previous investigations of the galaxy M31 (Andromeda) have shown that outside a radius of 25 kpc the population of halo globular clusters is rotating in alignment with the stellar disk, as are more centrally located clusters. The M31 halo also contains coherent stellar substructures, along with a smoothly distributed stellar component. Many of the globular clusters outside 25 kpc are associated with the most prominent substructures, while others are part of the smooth halo. Here we report a new analysis of the kinematics of these globular clusters. We find that the two distinct populations are rotating with perpendicular orientations. The rotation axis for the population associated with the smooth halo is aligned with the rotation axis for the plane of dwarf galaxies that encircles M31. We interpret these separate cluster populations as arising from two major accretion epochs, likely separated by billions of years. Stellar substructures from the first epoch are gone, but those from the more recent second epoch still remain. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: Authors' version of a Letter published in Nature on October 3rd, 2019

Towards precise stellar ages: combining isochrone fitting with empirical gyrochronology

Authors: Ruth Angus, Timothy D. Morton, Daniel Foreman-Mackey, Jennifer van Saders, Jason Curtis, Stephen R. Kane, Megan Bedell, Rocio Kiman, David W. Hogg, John Brewer

Abstract: We present a new age-dating technique that combines gyrochronology with isochrone fitting to infer ages for FGKM main-sequence and subgiant field stars. Gyrochronology and isochrone fitting are each capable of providing relatively precise ages for field stars in certain areas of the Hertzsprung-Russell diagram: gyrochronology works optimally for cool main-sequence stars, and isochrone fitting can… ▽ More We present a new age-dating technique that combines gyrochronology with isochrone fitting to infer ages for FGKM main-sequence and subgiant field stars. Gyrochronology and isochrone fitting are each capable of providing relatively precise ages for field stars in certain areas of the Hertzsprung-Russell diagram: gyrochronology works optimally for cool main-sequence stars, and isochrone fitting can provide precise ages for stars near the main-sequence turnoff. Combined, these two age-dating techniques can provide precise and accurate ages for a broader range of stellar masses and evolutionary stages than either method used in isolation. We demonstrate that the position of a star on the Hertzsprung- Russell or color-magnitude diagram can be combined with its rotation period to infer a precise age via both isochrone fitting and gyrochronology simultaneously. We show that incorporating rotation periods with 5% uncertainties into stellar evolution models improves age precision for FGK stars on the main sequence, and can, on average, provide age estimates up to three times more precise than isochrone fitting alone. In addition, we provide a new gyrochronology relation, calibrated to the Praesepe cluster and the Sun, that includes a variance model to capture the rotational behavior of stars whose rotation periods do not lengthen with the square-root of time, and parts of the Hertzsprung-Russell diagram where gyrochronology has not been calibrated. This publication is accompanied by an open source Python package, stardate, for inferring the ages of main-sequence and subgiant FGKM stars from rotation periods, spectroscopic parameters and/or apparent magnitudes and parallaxes. △ Less

Submitted 20 August, 2019; originally announced August 2019.

Comments: Accepted for publication in AJ

K2 rotation periods for low-mass Hyads and a quantitative comparison of the distribution of slow rotators in the Hyades and Praesepe

Authors: S. T. Douglas, J. L. Curtis, M. A. Agüeros, P. A. Cargile, J. M. Brewer, S. Meibom, T. Jansen

Abstract: We analyze K2 light curves for 132 low-mass ($1\ \gtrsim\ M_*\ \gtrsim\ 0.1$~${M_{\odot}}$) members of the 600--800~Myr-old Hyades cluster and measure rotation periods ($P_{rot}$) for 116 of these stars. These include 93 stars with no prior $P_{rot}$ measurement; the total number of Hyads with known $P_{rot}$ is now 232. We then combine literature binary data with Gaia DR2 photometry and astrometr… ▽ More We analyze K2 light curves for 132 low-mass ($1\ \gtrsim\ M_*\ \gtrsim\ 0.1$~${M_{\odot}}$) members of the 600--800~Myr-old Hyades cluster and measure rotation periods ($P_{rot}$) for 116 of these stars. These include 93 stars with no prior $P_{rot}$ measurement; the total number of Hyads with known $P_{rot}$ is now 232. We then combine literature binary data with Gaia DR2 photometry and astrometry to select single star sequences in the Hyades and its roughly coeval Praesepe open cluster, and derive a new reddening value of $A_V = 0.035$$\pm$$0.011$ for Praesepe. Comparing the effective temperature--$P_{rot}$ distributions for the Hyades and Praesepe, we find that solar-type Hyads rotate, on average, 0.4~d slower than their Praesepe counterparts. This $P_{rot}$ difference indicates that the Hyades is slightly older than Praesepe: we apply a new gyrochronology model tuned with Praesepe and the Sun, and find an age difference between the two clusters of 57~Myr. However, this $P_{rot}$ difference decreases and eventually disappears for lower-mass stars. This provides further evidence for stalling in the rotational evolution of these stars, and highlights the need for more detailed analysis of angular-momentum evolution for stars of different masses and ages. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: Accepted for publication in ApJ. 20 pages, 10 figures, 5 tables. CSV versions of tables 2, 3, and 4 available by request

The Mass of the White Dwarf Companion in the Self-Lensing Binary KOI-3278: Einstein vs. Newton

Authors: Daniel A. Yahalomi, Yossi Shvartzvald, Eric Agol, Avi Shporer, David W. Latham, Ethan Kruse, John M. Brewer, Lars A. Buchhave, Benjamin J. Fulton, Andrew W. Howard, Howard Isaacson, Erik A. Petigura, Samuel N. Quinn

Abstract: KOI-3278 is a self-lensing stellar binary consisting of a white-dwarf secondary orbiting a Sun-like primary star. Kruse and Agol (2014) noticed small periodic brightenings every 88.18 days in the Kepler photometry and interpreted these as the result of microlensing by a white dwarf with about 63$\%$ of the mass of the Sun. We obtained two sets of spectra for the primary that allowed us to derive t… ▽ More KOI-3278 is a self-lensing stellar binary consisting of a white-dwarf secondary orbiting a Sun-like primary star. Kruse and Agol (2014) noticed small periodic brightenings every 88.18 days in the Kepler photometry and interpreted these as the result of microlensing by a white dwarf with about 63$\%$ of the mass of the Sun. We obtained two sets of spectra for the primary that allowed us to derive three sets of spectroscopic estimates for its effective temperature, surface gravity, and metallicity for the first time. We used these values to update the Kruse and Agol (2014) Einsteinian microlensing model, resulting in a revised mass for the white dwarf of $0.539^{+0.022}_{-0.020} \, M_{\odot}$. The spectra also allowed us to determine radial velocities and derive orbital solutions, with good agreement between the two independent data sets. An independent Newtonian dynamical MCMC model of the combined velocities yielded a mass for the white dwarf of $0.5122^{+0.0057}_{-0.0058} \, M_{\odot}$. The nominal uncertainty for the Newtonian mass is about four times better than for the Einsteinian, $\pm 1.1\%$ vs. $\pm 4.1\%$ and the difference between the two mass determinations is $5.2 \%$. We then present a joint Einsteinian microlensing and Newtonian radial velocity model for KOI-3278, which yielded a mass for the white dwarf of $0.5250^{+0.0082}_{-0.0089} \, M_{\odot}$. This joint model does not rely on any white dwarf evolutionary models or assumptions on the white dwarf mass-radius relation. We discuss the benefits of a joint model of self-lensing binaries, and how future studies of these systems can provide insight into the mass-radius relation of white dwarfs. △ Less

Submitted 4 June, 2019; v1 submitted 24 April, 2019; originally announced April 2019.

Comments: ApJ Accepted; 22 Pages, 8 Figures, 6 Tables and 4 Supplementary Tables

Modeling the Echelle Spectra Continuum with Alpha Shapes and Local Regression Fitting

Authors: Xin Xu, Jessi Cisewski-Kehe, Allen B. Davis, Debra A. Fischer, John M. Brewer

Abstract: Continuum normalization of echelle spectra is an important data analysis step that is difficult to automate. Polynomial fitting requires a reasonably high order model to follow the steep slope of the blaze function. However, in the presence of deep spectral lines, a high order polynomial fit can result in ripples in the normalized continuum that increase errors in spectral analysis. Here, we prese… ▽ More Continuum normalization of echelle spectra is an important data analysis step that is difficult to automate. Polynomial fitting requires a reasonably high order model to follow the steep slope of the blaze function. However, in the presence of deep spectral lines, a high order polynomial fit can result in ripples in the normalized continuum that increase errors in spectral analysis. Here, we present two algorithms for flattening the spectrum continuum. The Alpha-shape Fitting to Spectrum algorithm (AFS) is completely data-driven, using an alpha shape to obtain an initial estimate of the blaze function. The Alpha-shape and Lab Source Fitting to Spectrum algorithm (ALSFS) incorporates a continuum constraint from a lab source reference spectrum for the blaze function estimation. These algorithms are tested on a simulated spectrum, where we demonstrate improved normalization compared to polynomial regression for continuum fitting. We show an additional application, using the algorithms for mitigation of spatially correlated quantum efficiency variations and fringing in the CCD detector of the EXtreme PREcision Spectrometer (EXPRES). △ Less

Submitted 22 April, 2019; originally announced April 2019.

The SAMI Galaxy Survey: Bayesian Inference for Gas Disk Kinematics using a Hierarchical Gaussian Mixture Model

Authors: Mathew R. Varidel, Scott M. Croom, Geraint F. Lewis, Brendon J. Brewer, Enrico M. Di Teodoro, Joss Bland-Hawthorn, Julia J. Bryant, Christoph Federrath, Caroline Foster, Karl Glazebrook, Michael Goodwin, Brent Groves, Andrew M. Hopkins, Jon S. Lawrence, Ángel R. López-Sánchez, Anne M. Medling, Matt S. Owers, Samuel N. Richards, Richard Scalzo, Nicholas Scott, Sarah M. Sweet, Dan S. Taranu, Jesse van de Sande

Abstract: We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To acco… ▽ More We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To account for beam smearing effects, we construct a modelled cube that is then convolved per wavelength slice by the seeing, before calculating the likelihood function. We show that our method can model complex gas substructure including clumps and spiral arms. We also show that kinematic asymmetries can be observed after beam smearing for regularly rotating galaxies with asymmetries only introduced in the spatial distribution of the gas. We present findings for our method applied to a sample of 20 star-forming galaxies from the SAMI Galaxy Survey. We estimate the global H$α$ gas velocity dispersion for our sample to be in the range $\barσ_v \sim $[7, 30] km s$^{-1}$. The relative difference between our approach and estimates using the single Gaussian component fits per spaxel is $Δ\barσ_v / \barσ_v = - 0.29 \pm 0.18$ for the H$α$ flux-weighted mean velocity dispersion. △ Less

Submitted 10 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 23 pages, 12 figures, accepted for MNRAS

Benchmarking Substellar Evolutionary Models Using New Age Estimates for HD 4747 B and HD 19467 B

Authors: Charlotte M. Wood, Tabetha Boyajian, Kaspar von Braun, John M. Brewer, Justin R. Crepp, Gail Schaefer, Arthur Adams, Timothy R. White

Abstract: Constraining substellar evolutionary models (SSEMs) is particularly difficult due to a degeneracy between the mass, age, and luminosity of a brown dwarf. In cases where a brown dwarf is found as a directly imaged companion to a star, as in HD 4747 and HD 19467, the mass, age, and luminosity of the brown dwarf are determined independently, making them ideal objects to use to benchmark SSEMs. Using… ▽ More Constraining substellar evolutionary models (SSEMs) is particularly difficult due to a degeneracy between the mass, age, and luminosity of a brown dwarf. In cases where a brown dwarf is found as a directly imaged companion to a star, as in HD 4747 and HD 19467, the mass, age, and luminosity of the brown dwarf are determined independently, making them ideal objects to use to benchmark SSEMs. Using the Center for High Angular Resolution Astronomy Array, we measured the angular diameters and calculated the radii of the host stars HD 4747 A and HD 19467 A. After fitting their parameters to the Dartmouth Stellar Evolution Database, MESA Isochrones and Stellar Tracks, and Yonsei-Yale isochronal models, we adopt age estimates of $10.74^{+6.75}_{-6.87}$ Gyr for HD 4747 A and $10.06^{+1.16}_{-0.82}$ Gyr for HD 19467 A. Assuming the brown dwarf companions HD 4747 B and HD 19467 B have the same ages as their host stars, we show that many of the SSEMs under-predict bolometric luminosities by $\sim$ 0.75 dex for HD 4747 B and $\sim 0.5$ dex for HD 19467 B. The discrepancies in luminosity correspond to over-predictions of the masses by $\sim$ 12\% for HD 4747 B and $\sim$ 30\% for HD 19467 B. We also show that SSEMs that take into account the effect of clouds reduce the under-prediction of luminosity to $\sim 0.6$ dex and the over-prediction of mass to $\sim 8\%$ for HD 4747 B, an L/T transition object that is cool enough to begin forming clouds. One possible explanation for the remaining discrepancies is missing physics in the models, such as the inclusion of metallicity effects. △ Less

Submitted 11 January, 2019; originally announced January 2019.

Comments: 12 pages, 6 figures, 5 tables, accepted to ApJ

Compact multi-planet systems are more common around metal poor hosts

Authors: John M. Brewer, Songhu Wang, Debra A. Fischer, Daniel Foreman-Mackey

Abstract: In systems with detected planets, hot-Jupiters and compact systems of multiple planets are nearly mutually exclusive. We compare the relative occurrence of these two architectures as a fraction of detected planetary systems to determine the role that metallicity plays in planet formation. We show that compact multi-planet systems occur more frequently around stars of increasingly lower metalliciti… ▽ More In systems with detected planets, hot-Jupiters and compact systems of multiple planets are nearly mutually exclusive. We compare the relative occurrence of these two architectures as a fraction of detected planetary systems to determine the role that metallicity plays in planet formation. We show that compact multi-planet systems occur more frequently around stars of increasingly lower metallicities using spectroscopically derived abundances for more than 700 planet hosts. At higher metallicities, compact multi-planet systems comprise a nearly constant fraction of the planet hosts despite the steep rise in the fraction of hosts containing hot and cool-Jupiters. Since metal poor stars have been underrepresented in planet searches, this implies that the occurrence rate of compact multis is higher than previously reported. Due to observational limits, radial velocity planet searches have focused mainly on high-metallicity stars where they have a higher chance of finding giant planets. New extreme-precision radial velocity instruments coming online that can detect these compact multi-planet systems can target lower metallicity stars to find them. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 5 pages, 2 figures, accepted to ApJ Letters

HD 202772A B: A Transiting Hot Jupiter Around A Bright, Mildly Evolved Star In A Visual Binary Discovered By Tess

Authors: Songhu Wang, Matias Jones, Avi Shporer, Benjamin J. Fulton, Leonardo A. Paredes, Trifon Trifonov, Diana Kossakowski, Jason Eastman, Maximilian N. Gunther, Chelsea X. Huang, Sarah Millholland, Darryl Seligman, Debra Fischer, Rafael Brahm, Xian-Yu Wang, Bryndis Cruz, Hodari-Sadiki James, Brett Addison, Todd Henry, En-Si Liang, Allen B. Davis, Rene Tronsgaard, Keduse Worku, John Brewer, Martin Kurster , et al. (30 additional authors not shown)

Abstract: We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial-velocity measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the s… ▽ More We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial-velocity measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the star is among the brightest known to host a hot Jupiter. Based on the 27days of TESS photometry, and radial velocity data from the CHIRON and HARPS spectrographs, the planet has a mass of 1.008+/-0.074 M_J and radius of 1.562+/-0.053 R_J , making it an inflated gas giant. HD 202772A b is a rare example of a transiting hot Jupiter around a quickly evolving star. It is also one of the most strongly irradiated hot Jupiters currently known. △ Less

Submitted 4 October, 2018; originally announced October 2018.

Comments: Submitted to AAS Journal

The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad-Line Region

Authors: Peter R. Williams, Anna Pancoast, Tommaso Treu, Brendon J. Brewer, Aaron J. Barth, Vardha N. Bennert, Tabitha Buehler, Gabriela Canalizo, S. Bradley Cenko, Kelsey I. Clubb, Michael C. Cooper, Alexei V. Filippenko, Elinor Gates, Sebastian F. Hoenig, Michael D. Joner, Michael T. Kandrashoff, Clifton David Laney, Mariana S. Lazarova, Weidong Li, Matthew A. Malkan, Jacob Rex, Jeffrey M. Silverman, Erik Tollerud, Jonelle L. Walsh, Jong-Hak Woo

Abstract: We present models of the H$β$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ioni… ▽ More We present models of the H$β$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ionizing source. The dynamics in our sample range from near-circular elliptical orbits to inflowing or outflowing trajectories. We measure black hole masses of $\log_{10}(M_{\rm BH}/M_\odot) = 6.48^{+0.21}_{-0.18}$ for PG 1310$-$108, $7.50^{+0.25}_{-0.18}$ for Mrk 50, $7.46^{+0.15}_{-0.21}$ for Mrk 141, $7.58^{+0.08}_{-0.08}$ for Mrk 279, $7.11^{+0.20}_{-0.17}$ for Mrk 1511, $6.65^{+0.27}_{-0.15}$ for NGC 4593, and $6.94^{+0.14}_{-0.14}$ for Zw 229$-$015. We use these black hole mass measurements along with cross-correlation time lags and line widths to recover the scale factor $f$ used in traditional reverberation mapping measurements. Combining our results with other studies that use this modeling technique, bringing our sample size to 16, we calculate a scale factor that can be used for measuring black hole masses in other reverberation mapping campaigns. When using the root-mean-square (rms) spectrum and using the line dispersion to measure the line width, we find $\log_{10}(f_{{\rm rms},σ})_{\rm pred} = 0.57 \pm 0.19$. Finally, we search for correlations between $f$ and other AGN and BLR parameters and find marginal evidence that $f$ is correlated with $M_{\rm BH}$ and the BLR inclination angle, but no significant evidence of a correlation with the AGN luminosity or Eddington ratio. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Comments: 26 pages, 14 figures. Accepted for publication in ApJ

kima: Exoplanet detection in radial velocities

Authors: J. P. Faria, N. C. Santos, P. Figueira, B. J. Brewer

Abstract: The radial-velocity (RV) method is one of the most successful in the detection of exoplanets, but is hindered by the intrinsic RV variations of the star, which can easily mimic or hide true planetary signals. kima is a package for the detection and characterization of exoplanets using RV data. It fits a sum of Keplerian curves to a timeseries of RV measurements and calculates the evidence for mode… ▽ More The radial-velocity (RV) method is one of the most successful in the detection of exoplanets, but is hindered by the intrinsic RV variations of the star, which can easily mimic or hide true planetary signals. kima is a package for the detection and characterization of exoplanets using RV data. It fits a sum of Keplerian curves to a timeseries of RV measurements and calculates the evidence for models with a fixed number Np of Keplerian signals, or after marginalising over Np. Moreover, kima can use a GP with a quasi-periodic kernel as a noise model, to deal with activity-induced signals. The hyperparameters of the GP are inferred together with the orbital parameters. The code is written in C++, but includes a helper Python package, pykima, which facilitates the analysis of the results. △ Less

Submitted 21 June, 2018; originally announced June 2018.

Comments: 3 pages; published in the Journal of Open Source Software; all code available at https://github.com/j-faria/kima

Journal ref: Journal of Open Source Software, 3(26), 487, 2018

Spectral Properties of Cool Stars: Extended Abundance Analysis of Kepler Objects of Interest

Authors: John M. Brewer, Debra A. Fischer

Abstract: Accurate stellar parameters and precise elemental abundances are vital pieces to correctly characterize discovered planetary systems, better understand planet formation, and trace galactic chemical evolution. We have performed a uniform spectroscopic analysis for 1127 stars, yielding accurate gravity, temperature, and projected rotational velocity in addition to precise abundances for 15 elements… ▽ More Accurate stellar parameters and precise elemental abundances are vital pieces to correctly characterize discovered planetary systems, better understand planet formation, and trace galactic chemical evolution. We have performed a uniform spectroscopic analysis for 1127 stars, yielding accurate gravity, temperature, and projected rotational velocity in addition to precise abundances for 15 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, and Y). Most of the stars in this sample are Kepler Objects of Interest, observed by the California-Kepler Survey (CKS) and include 1,003 stars hosting 1,562 confirmed planets. This catalog extends the uniform analysis of our previous catalog, bringing the total of homogeneously analyzed stars to almost 2,700 F, G, and K dwarfs. To ensure consistency between the catalogs, we performed an analysis of our ability to recover parameters as a function of S/N ratio and present individual uncertainties as well as functions to calculate uncertainties for parameters derived from lower S/N ratio spectra. With the updated parameters, we used isochrone fitting to derived new radii, masses and ages for the stars. Finally, we look at the Mg/Si ratios of super-Earth and sub-Neptune hosts to test whether differences in initial composition might lead to differences in planet radius. We find no differences in the Mg/Si distribution as a function of planet radius. △ Less

Submitted 30 July, 2018; v1 submitted 2 April, 2018; originally announced April 2018.

Comments: 19 pages, 12 figures, 6 tables, accepted to ApJ Supplement. Full data tables available in machine readable format in source file

Stability of the Broad Line Region Geometry and Dynamics in Arp 151 Over Seven Years

Authors: A. Pancoast, A. J. Barth, K. Horne, T. Treu, B. J. Brewer, V. N. Bennert, G. Canalizo, E. L. Gates, W. Li, M. A. Malkan, D. Sand, T. Schmidt, S. Valenti, J. H. Woo, K. I. Clubb, M. C. Cooper, S. M. Crawford, S. F. Honig, M. D. Joner, M. T. Kandrashoff, M. Lazarova, A. M. Nierenberg, E. Romero-Colmenero, D. Son, E. Tollerud , et al. (2 additional authors not shown)

Abstract: The Seyfert 1 galaxy Arp 151 was monitored as part of three reverberation mapping campaigns spanning $2008-2015$. We present modeling of these velocity-resolved reverberation mapping datasets using a geometric and dynamical model for the broad line region (BLR). By modeling each of the three datasets independently, we infer the evolution of the BLR structure in Arp 151 over a total of seven years… ▽ More The Seyfert 1 galaxy Arp 151 was monitored as part of three reverberation mapping campaigns spanning $2008-2015$. We present modeling of these velocity-resolved reverberation mapping datasets using a geometric and dynamical model for the broad line region (BLR). By modeling each of the three datasets independently, we infer the evolution of the BLR structure in Arp 151 over a total of seven years and constrain the systematic uncertainties in non-varying parameters such as the black hole mass. We find that the BLR geometry of a thick disk viewed close to face-on is stable over this time, although the size of the BLR grows by a factor of $\sim 2$. The dynamics of the BLR are dominated by inflow and the inferred black hole mass is consistent for the three datasets, despite the increase in BLR size. Combining the inference for the three datasets yields a black hole mass and statistical uncertainty of $\log_{10}($M$_{\rm BH}/\rm{M}_{\odot})=6.82^{+0.09}_{-0.09}$ with a standard deviation in individual measurements of 0.13 dex. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 21 pages, 8 figures, 2 tables. Accepted for publication in ApJ

Stellar Spin-Orbit Alignment for Kepler-9, a Multi-transiting Planetary system with Two Outer Planets Near 2:1 Resonance

Authors: Songhu Wang, Brett Addison, Debra A. Fischer, John M. Brewer, Howard Isaacson, Andrew W. Howard, Gregory Laughlin

Abstract: We present spectroscopic measurements of the Rossiter-McLaughlin effect for the planet b of Kepler-9 multi-transiting planet system. The resulting sky-projected spin-orbit angle is $λ=-13^{\circ} \pm 16^{\circ}$, which favors an aligned system and strongly disfavors highly misaligned, polar, and retrograde orbits. Including Kepler-9, there are now a total of 4 Rossiter-McLaughlin effect measuremen… ▽ More We present spectroscopic measurements of the Rossiter-McLaughlin effect for the planet b of Kepler-9 multi-transiting planet system. The resulting sky-projected spin-orbit angle is $λ=-13^{\circ} \pm 16^{\circ}$, which favors an aligned system and strongly disfavors highly misaligned, polar, and retrograde orbits. Including Kepler-9, there are now a total of 4 Rossiter-McLaughlin effect measurements for multiplanet systems, all of which are consistent with spin-orbit alignment. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: Accepted for publication in the Astronomical Journal

Planet Detectability in the Alpha Centauri System

Authors: Lily L. Zhao, Debra A. Fischer, John M. Brewer, Matt Giguere, Bárbara Rojas-Ayala

Abstract: We use more than a decade of radial velocity measurements for $α$ Cen A, B, and Proxima Centauri from HARPS, CHIRON, and UVES to identify the $M \sin i$ and orbital periods of planets that could have been detected if they existed. At each point in a mass-period grid, we sample a simulated, Keplerian signal with the precision and cadence of existing data and assess the probability that the signal c… ▽ More We use more than a decade of radial velocity measurements for $α$ Cen A, B, and Proxima Centauri from HARPS, CHIRON, and UVES to identify the $M \sin i$ and orbital periods of planets that could have been detected if they existed. At each point in a mass-period grid, we sample a simulated, Keplerian signal with the precision and cadence of existing data and assess the probability that the signal could have been produced by noise alone. Existing data places detection thresholds in the classically defined habitable zones at about $M \sin i$ of 53 M$_{\oplus}$ for $α$ Cen A, 8.4 M$_{\oplus}$ for $α$ Cen B, and 0.47 M$_{\oplus}$ for Proxima Centauri. Additionally, we examine the impact of systematic errors, or "red noise" in the data. A comparison of white- and red-noise simulations highlights quasi-periodic variability in the radial velocities that may be caused by systematic errors, photospheric velocity signals, or planetary signals. For example, the red-noise simulations show a peak above white-noise simulations at the period of Proxima Centauri b. We also carry out a spectroscopic analysis of the chemical composition of the $α$ Centauri stars. The stars have super-solar metallicity with ratios of C/O and Mg/Si that are similar to the Sun, suggesting that any small planets in the $α$ Cen system may be compositionally similar to our terrestrial planets. Although the small projected separation of $α$ Cen A and B currently hampers extreme-precision radial velocity measurements, the angular separation is now increasing. By 2019, $α$ Cen A and B will be ideal targets for renewed Doppler planet surveys. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 16 pages, 7 figures, data provided in appendix. Accepted for publication in ApJ

A physically motivated and empirically calibrated method to measure effective temperature, metallicity, and Ti abundance of M dwarfs

Authors: Mark J. Veyette, Philip S. Muirhead, Andrew W. Mann, John M. Brewer, France Allard, Derek Homeier

Abstract: The ability to perform detailed chemical analysis of Sun-like F-, G-, and K-type stars is a powerful tool with many applications including studying the chemical evolution of the Galaxy and constraining planet formation theories. Unfortunately, complications in modeling cooler stellar atmospheres hinders similar analysis of M-dwarf stars. Empirically-calibrated methods to measure M dwarf metallicit… ▽ More The ability to perform detailed chemical analysis of Sun-like F-, G-, and K-type stars is a powerful tool with many applications including studying the chemical evolution of the Galaxy and constraining planet formation theories. Unfortunately, complications in modeling cooler stellar atmospheres hinders similar analysis of M-dwarf stars. Empirically-calibrated methods to measure M dwarf metallicity from moderate-resolution spectra are currently limited to measuring overall metallicity and rely on astrophysical abundance correlations in stellar populations. We present a new, empirical calibration of synthetic M dwarf spectra that can be used to infer effective temperature, Fe abundance, and Ti abundance. We obtained high-resolution (R~25,000), Y-band (~1 micron) spectra of 29 M dwarfs with NIRSPEC on Keck II. Using the PHOENIX stellar atmosphere modeling code (version 15.5), we generated a grid of synthetic spectra covering a range of temperatures, metallicities, and alpha-enhancements. From our observed and synthetic spectra, we measured the equivalent widths of multiple Fe I and Ti I lines and a temperature-sensitive index based on the FeH bandhead. We used abundances measured from widely-separated solar-type companions to empirically calibrate transformations to the observed indices and equivalent widths that force agreement with the models. Our calibration achieves precisions in Teff, [Fe/H], and [Ti/Fe] of 60 K, 0.1 dex, and 0.05 dex, respectively and is calibrated for 3200 K < Teff < 4100 K, -0.7 < [Fe/H] < +0.3, and -0.05 < [Ti/Fe] < +0.3. This work is a step toward detailed chemical analysis of M dwarfs at a similar precision achieved for FGK stars. △ Less

Submitted 27 October, 2017; originally announced October 2017.

Comments: accepted for publication in ApJ, all synthetic spectra available at http://people.bu.edu/mveyette/phoenix/

Kronos & Krios: Evidence for accretion of a massive, rocky planetary system in a comoving pair of solar-type stars

Authors: Semyeong Oh, Adrian M. Price-Whelan, John M. Brewer, David W. Hogg, David N. Spergel, Justin Myles

Abstract: We report and discuss the discovery of a comoving pair of bright solar-type stars, HD 240430 and HD 240429, with a significant difference in their chemical abundances. The two stars have an estimated 3D separation of $\approx 0.6$ pc ($\approx 0.01$ pc projected) at a distance of $r\approx 100$ pc with nearly identical three-dimensional velocities, as inferred from Gaia TGAS parallaxes and proper… ▽ More We report and discuss the discovery of a comoving pair of bright solar-type stars, HD 240430 and HD 240429, with a significant difference in their chemical abundances. The two stars have an estimated 3D separation of $\approx 0.6$ pc ($\approx 0.01$ pc projected) at a distance of $r\approx 100$ pc with nearly identical three-dimensional velocities, as inferred from Gaia TGAS parallaxes and proper motions, and high-precision radial velocity measurements. Stellar parameters determined from high-resolution Keck HIRES spectra indicate that both stars are $\sim 4$ Gyr old. The more metal-rich of the two, HD 240430, shows an enhancement of refractory ($T_C>1200$ K) elements by $\approx 0.2$ dex and a marginal enhancement of (moderately) volatile elements ($T_C<1200$ K, C, N, O, Na, and Mn). This is the largest metallicity difference found in a wide binary pair yet. Additionally, HD 240430 shows an anomalously high surface lithium abundance ($A(\mathrm{Li})=2.75$), higher than its companion by $0.5$ dex. The proximity in phase-space and ages between the two stars suggests that they formed together with the same composition, at odds with the observed differences in metallicity and abundance patterns. We therefore suggest that the star HD~240430, "Kronos", accreted 15 $M_\oplus$ of rocky material after birth, selectively enhancing the refractory elements as well as lithium in its surface and convective envelope. △ Less

Submitted 15 September, 2017; originally announced September 2017.

Comments: Figure 9 is the main figure, submitted

Computing Entropies With Nested Sampling

Authors: Brendon J. Brewer

Abstract: The Shannon entropy, and related quantities such as mutual information, can be used to quantify uncertainty and relevance. However, in practice, it can be difficult to compute these quantities for arbitrary probability distributions, particularly if the probability mass functions or densities cannot be evaluated. This paper introduces a computational approach, based on Nested Sampling, to evaluate… ▽ More The Shannon entropy, and related quantities such as mutual information, can be used to quantify uncertainty and relevance. However, in practice, it can be difficult to compute these quantities for arbitrary probability distributions, particularly if the probability mass functions or densities cannot be evaluated. This paper introduces a computational approach, based on Nested Sampling, to evaluate entropies of probability distributions that can only be sampled. I demonstrate the method on three examples: a simple gaussian example where the key quantities are available analytically; (ii) an experimental design example about scheduling observations in order to measure the period of an oscillating signal; and (iii) predicting the future from the past in a heavy-tailed scenario. △ Less

Submitted 15 August, 2017; v1 submitted 12 July, 2017; originally announced July 2017.

Comments: Accepted for publication in Entropy. 21 pages, 3 figures. Software available at https://github.com/eggplantbren/InfoNest

Three's Company: An additional non-transiting super-Earth in the bright HD 3167 system, and masses for all three planets

Authors: Jessie L. Christiansen, Andrew Vanderburg, Jennifer Burt, B. J. Fulton, Konstantin Batygin, Björn Benneke, John M. Brewer, David Charbonneau, David R. Ciardi, Andrew Collier Cameron, Jeffrey L. Coughlin, Ian J. M. Crossfield, Courtney Dressing, Thomas P. Greene, Andrew W. Howard, David W. Latham, Emilio Molinari, Annelies Mortier, Fergal Mullally, Francesco Pepe, Ken Rice, Evan Sinukoff, Alessandro Sozzetti, Susan E. Thompson, Stéphane Udry , et al. (33 additional authors not shown)

Abstract: HD 3167 is a bright (V = 8.9), nearby K0 star observed by the NASA K2 mission (EPIC 220383386), hosting two small, short-period transiting planets. Here we present the results of a multi-site, multi-instrument radial velocity campaign to characterize the HD 3167 system. The masses of the transiting planets are 5.02+/-0.38 MEarth for HD 3167 b, a hot super-Earth with a likely rocky composition (rho… ▽ More HD 3167 is a bright (V = 8.9), nearby K0 star observed by the NASA K2 mission (EPIC 220383386), hosting two small, short-period transiting planets. Here we present the results of a multi-site, multi-instrument radial velocity campaign to characterize the HD 3167 system. The masses of the transiting planets are 5.02+/-0.38 MEarth for HD 3167 b, a hot super-Earth with a likely rocky composition (rho_b = 5.60+2.15-1.43 g/cm^3), and 9.80+1.30-1.24 MEarth for HD 3167 c, a warm sub-Neptune with a likely substantial volatile complement (rho_c = 1.97+0.94-0.59 g/cm^3). We explore the possibility of atmospheric composition analysis and determine that planet c is amenable to transmission spectroscopy measurements, and planet b is a potential thermal emission target. We detect a third, non-transiting planet, HD 3167 d, with a period of 8.509+/-0.045 d (between planets b and c) and a minimum mass of 6.90+/-0.71 MEarth. We are able to constrain the mutual inclination of planet d with planets b and c: we rule out mutual inclinations below 1.3 degrees as we do not observe transits of planet d. From 1.3-40 degrees, there are viewing geometries invoking special nodal configurations which result in planet d not transiting some fraction of the time. From 40-60 degrees, Kozai-Lidov oscillations increase the system's instability, but it can remain stable for up to 100Myr. Above 60 degrees, the system is unstable. HD 3167 promises to be a fruitful system for further study and a preview of the many exciting systems expected from the upcoming NASA TESS mission. △ Less

Submitted 6 June, 2017; originally announced June 2017.

Comments: 22 pages, 14 figures, 5 tables. Submitted to AJ March 3rd, 2017. Accepted April 28th, 2017. In press

K2-66b and K2-106b: Two extremely hot sub-Neptune-size planets with high densities

Authors: Evan Sinukoff, Andrew W. Howard, Erik A. Petigura, Benjamin J. Fulton, Ian J. M. Crossfield, Howard Isaacson, Erica Gonzales, Justin R. Crepp, John M. Brewer, Lea Hirsch, Lauren M. Weiss, David R. Ciardi, Joshua E. Schlieder, Bjoern Benneke, Jessie L. Christiansen, Courtney D. Dressing, Brad M. S. Hansen, Heather A. Knutson, Molly Kosiarek, John H. Livingston, Thomas P. Greene, Leslie A. Rogers, Sebastien Lepine

Abstract: We report precise mass and density measurements of two extremely hot sub-Neptune-size planets from the K2 mission using radial velocities, K2 photometry, and adaptive optics imaging. K2-66 harbors a close-in sub-Neptune-sized (2.49$^{+0.34}_{-0.24} R_\oplus$) planet (K2-66b) with a mass of 21.3 $\pm$ 3.6 $M_\oplus$. Because the star is evolving up the sub-giant branch, K2-66b receives a high level… ▽ More We report precise mass and density measurements of two extremely hot sub-Neptune-size planets from the K2 mission using radial velocities, K2 photometry, and adaptive optics imaging. K2-66 harbors a close-in sub-Neptune-sized (2.49$^{+0.34}_{-0.24} R_\oplus$) planet (K2-66b) with a mass of 21.3 $\pm$ 3.6 $M_\oplus$. Because the star is evolving up the sub-giant branch, K2-66b receives a high level of irradiation, roughly twice the main sequence value. K2-66b may reside within the so-called "photoevaporation desert", a domain of planet size and incident flux that is almost completely devoid of planets. Its mass and radius imply that K2-66b has, at most, a meager envelope fraction (< 5%) and perhaps no envelope at all, making it one of the largest planets without a significant envelope. K2-106 hosts an ultra-short-period planet ($P$ = 13.7 hrs) that is one of the hottest sub-Neptune-size planets discovered to date. Its radius (1.82$^{+0.20}_{-0.14} R_\oplus$) and mass (9.0 $\pm$ 1.6 $M_\oplus$) are consistent with a rocky composition, as are all other small ultra-short-period planets with well-measured masses. K2-106 also hosts a larger, longer-period planet (Rp = 2.77$^{+0.37}_{-0.23} R_\oplus$, $P$ = 13.3 days) with a mass less than 24.4 $M_\oplus$ at 99.7% confidence. K2-66b and K2-106b probe planetary physics in extreme radiation environments. Their high densities reflect the challenge of retaining a substantial gas envelope in such extreme environments. △ Less

Submitted 9 May, 2017; originally announced May 2017.

Comments: 14 pages, 8 figures, 5 tables, Accepted to AJ

The Structure of the Broad-Line Region In Active Galactic Nuclei. II. Dynamical Modeling of Data from the AGN10 Reverberation Mapping Campaign

Authors: C. J. Grier, A. Pancoast, A. J. Barth, M. M. Fausnaugh, B. J. Brewer, T. Treu, B. M. Peterson

Abstract: We present inferences on the geometry and kinematics of the broad-Hbeta line-emitting region in four active galactic nuclei monitored as a part of the fall 2010 reverberation mapping campaign at MDM Observatory led by the Ohio State University. From modeling the continuum variability and response in emission-line profile changes as a function of time, we infer the geometry of the Hbeta- emitting b… ▽ More We present inferences on the geometry and kinematics of the broad-Hbeta line-emitting region in four active galactic nuclei monitored as a part of the fall 2010 reverberation mapping campaign at MDM Observatory led by the Ohio State University. From modeling the continuum variability and response in emission-line profile changes as a function of time, we infer the geometry of the Hbeta- emitting broad line regions to be thick disks that are close to face-on to the observer with kinematics that are well-described by either elliptical orbits or inflowing gas. We measure the black hole mass to be log (MBH) = 7.25 (+/-0.10) for Mrk 335, 7.86 (+0.20, -0.17) for Mrk 1501, 7.84 (+0.14, -0.19) for 3C 120, and 6.92 (+0.24, -0.23) for PG 2130+099. These black hole mass measurements are not based on a particular assumed value of the virial scale factor f, allowing us to compute individual f factors for each target. Our results nearly double the number of targets that have been modeled in this manner, and investigate the properties of a more diverse sample by including previously modeled objects. We measure an average scale factor f in the entire sample to be log10(f) = 0.54 +/- 0.17 when the line dispersion is used to characterize the line width, which is consistent with values derived using the normalization of the MBH-sigma relation. We find that the scale factor f for individual targets is likely correlated with the black hole mass, inclination angle, and opening angle of the broad line region but we do not find any correlation with the luminosity. △ Less

Submitted 27 September, 2017; v1 submitted 5 May, 2017; originally announced May 2017.

Comments: 19 pages, 13 Figures. Accepted for publication in ApJ

The Yale-Potsdam Stellar Isochrones (YaPSI)

Authors: F. Spada, P. Demarque, Y. -C. Kim, T. S. Boyajian, J. M. Brewer

Abstract: We introduce the Yale-Potsdam Stellar Isochrones (YaPSI), a new grid of stellar evolution tracks and isochrones of solar-scaled composition. In an effort to improve the Yonsei-Yale database, special emphasis is placed on the construction of accurate low-mass models (Mstar < 0.6 Msun), and in particular of their mass-luminosity and mass-radius relations, both crucial in characterizing exoplanet-hos… ▽ More We introduce the Yale-Potsdam Stellar Isochrones (YaPSI), a new grid of stellar evolution tracks and isochrones of solar-scaled composition. In an effort to improve the Yonsei-Yale database, special emphasis is placed on the construction of accurate low-mass models (Mstar < 0.6 Msun), and in particular of their mass-luminosity and mass-radius relations, both crucial in characterizing exoplanet-host stars and, in turn, their planetary systems. The YaPSI models cover the mass range 0.15 to 5.0 Msun, densely enough to permit detailed interpolation in mass, and the metallicity and helium abundance ranges [Fe/H] = -1.5 to +0.3, and Y = 0.25 to 0.37, specified independently of each other (i.e., no fixed Delta Y/Delta Z relation is assumed). The evolutionary tracks are calculated from the pre-main sequence up to the tip of the red giant branch. The isochrones, with ages between 1 Myr and 20 Gyr, provide UBVRI colors in the Johnson-Cousins system, and JHK colors in the homogeneized Bessell & Brett system, derived from two different semi-empirical Teff-color calibrations from the literature. We also provide utility codes, such as an isochrone interpolator in age, metallicity, and helium content, and an interface of the tracks with an open-source Monte Carlo Markov-Chain tool for the analysis of individual stars. Finally, we present comparisons of the YaPSI models with the best empirical mass- luminosity and mass-radius relations available to date, as well as isochrone fitting of well-studied ste △ Less

Submitted 11 March, 2017; originally announced March 2017.

Comments: 17 pages, 14 figures; accepted for publication in the Astrophysical Journal