DOCSLIB.ORG

The Generalised Lomb-Scargle Periodogram a New Formalism for the floating-Mean and Keplerian Periodograms M

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
The Generalised Lomb-Scargle Periodogram a New Formalism for the floating-Mean and Keplerian Periodograms M A&A 496, 577–584 (2009) Astronomy DOI: 10.1051/0004-6361:200811296 & c ESO 2009 Astrophysics The generalised Lomb-Scargle periodogram A new formalism for the floating-mean and Keplerian periodograms M. Zechmeister and M. Kürster Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany e-mail: [email protected] Received 5 November 2008 / Accepted 19 December 2008 ABSTRACT The Lomb-Scargle periodogram is a common tool in the frequency analysis of unequally spaced data equivalent to least-squares fitting of sine waves. We give an analytic solution for the generalisation to a full sine wave fit, including an offset and weights (χ2 fitting). Compared to the Lomb-Scargle periodogram, the generalisation is superior as it provides more accurate frequencies, is less susceptible to aliasing, and gives a much better determination of the spectral intensity. Only a few modifications are required for the computation and the computational effort is similar. Our approach brings together several related methods that can be found in the literature, viz. the date-compensated discrete Fourier transform, the floating-mean periodogram, and the “spectral significance” estimator used in the SigSpec program, for which we point out some equivalences. Furthermore, we present an algorithm that implements this generalisation for the evaluation of the Keplerian periodogram that searches for the period of the best-fitting Keplerian orbit to radial velocity data. The systematic and non-random algorithm is capable of detecting eccentric orbits, which is demonstrated by two examples and can be a useful tool in searches for the orbital periods of exoplanets. Key words. methods: data analysis – methods: analytical – methods: statistical – techniques: radial velocities 1. Introduction the mean of the data was subtracted, which assumes that the mean of the data and the mean of the fitted sine function are the The Lomb-Scargle periodogram (Scargle 1982) is a widely used same. One can overcome this assumption with the introduction tool in period searches and frequency analysis of time series. It is ff y = ω + ω of an o set c, resulting in a further generalisation of this peri- equivalent to fitting sine waves of the form a cos t b sin t. odogram to the equivalent of weighted full sine wave fitting; i.e., While standard fitting procedures require the solution of a set of y = a cos ωt + b sin ωt + c. Cumming et al. (1999), who called linear equations for each sampled frequency, the Lomb-Scargle this generalisation “floating-mean periodogram”, argue that this method provides an analytic solution and is therefore both con- ffi approach is superior: “... the Lomb-Scargle periodogram fails to venient to use and e cient. The equation for the periodogram account for statistical fluctuations in the mean of a sampled si- was given by Barning (1963), and also Lomb (1976)andScargle nusoid, making it non-robust when the number of observations (1982), who furthermore investigated its statistical behaviour, is small, the sampling is uneven, or for periods comparable to especially the statistical significance of the detection of a signal. y y = or greater than the duration of the observations”. These authors For a time series (ti, i) with zero mean ( 0), the Lomb- provided a formal definition and also a sophisticated statistical Scargle periodogram is defined as (normalisation from Lomb treatment, but do not use an analytical solution for the computa- 1976): tion of this periodogram. ⎡ ⎤ ⎢ 2 2 ⎥ Basically, analytical formulae for a full sine, least-squares 1 ⎢ YCˆ τ YSˆ τ ⎥ pˆ(ω) = ⎣⎢ ˆ + ˆ ⎦⎥ (1) spectrum have already been given by Ferraz-Mello (1981), call- YYˆ CCˆ τˆ SSˆ τˆ ing this date-compensated discrete Fourier transform (DCDFT). ⎧ ⎫ ⎪ 2 2 ⎪ We prefer to adopt a notation closely related to the Lomb- 1 ⎨ i yi cos ω(ti − τˆ) i yi sin ω(ti − τˆ) ⎬ = ⎪ + ⎪(2) Scargle periodogram calling it the generalised Lomb-Scargle pe- y2 ⎩ 2 ω − τ 2 ω − τ ⎭ i i i cos (ti ˆ) i sin (ti ˆ) riodogram (GLS). Shrager (2001) tries for such an approach but did not generalise the parameterτ ˆ in Eq. (3). Moreover, our gen- where the hats are used in this paper to symbolise the classical eralised equations, which are derived in the following (Sect. 2), expressions. The parameterτ ˆ is calculated via have a comparable symmetry to the classical ones and also allow us to point out equivalences to the “spectral significance” estima- sin 2ωt tan 2ωτˆ = i i · (3) tor used in the SigSpec program by Reegen (2007) (Sect. 4). i cos 2ωti However, there are two shortcomings. First, the Lomb- 2. The generalised Lomb-Scargle periodogram Scargle periodogram does not take the measurement errors into (GLS) account. This was solved by introducing weighted sums by Gilliland & Baliunas (1987)andIrwin et al. (1989) (equiva- The analytic solution for the generalised Lomb-Scargle peri- lent to the generalisation to a χ2 fit). Second, for the analysis odogram can be obtained in a straightforward manner in the Article published by EDP Sciences 578 M. Zechmeister and M. Kürster: The generalised Lomb-Scargle periodogram same way as outlined in Lomb (1976). Let yi be the N mea- So the sums YC and YS use the weighted mean subtracted data surements of a time series at time ti and with errors σi. Fitting a and are calculated in the same way as for the Lomb-Scargle pe- full sine function (i.e. including an offset c): riodogram (but with weights). ω y(t) = a cos ωt + b sin ωt + c The generalised Lomb-Scargle periodogram p( )inEq.(4) is normalised to unity and therefore in the range of 0 ≤ p ≤ 1, 2π at given frequency ω (or period P = ω ) means to minimise with p = 0 indicating no improvement of the fit and p = 1a 2 2 the squared difference between the data yi and the model func- “perfect” fit (100% reduction of χ or χ = 0). tion y(t): As the full sine fit is time-translation invariant, there is also the possibility to introduce an arbitrary time reference N y − y 2 2 [ i (ti)] 2 τ t → t − τ CC = w 2 ω t − τ − χ = = W wi[yi − y(ti)] point ( i i ;now,e.g. i cos ( i ) σ2 w ω − τ 2 ff χ2 i=1 i ( i cos (ti )) ), which will not a ect the of the fit. If this parameter τ is chosen as where 1 1 w = = 1 w = ωτ = 2CS i W σ2 i 1 tan 2 σ2 i − W i CC SS w sin 2ωt −2 w cos ωt w sin ωt are the normalised weights1. Minimisation leads to a system of = i i i i i i (19) 2 2 (three) linear equations whose solution is derived in detail in wi cos 2ωti − ( wi cos ωti) −( wi sin ωti) Appendix A.1. Furthermore, it is shown in Appendix A.1 that 2 = w ω − the relative χ -reduction p(ω) as a function of frequency ω and the interaction term in Eq. (5) disappears, CS τ i cos (ti χ2 χ2 τ)sinω(ti − τ) − wi cos ω(ti − τ) wi sin ω(ti − τ) = 0 (proof normalised to unity by 0 (the for the weighted mean) can be written as: in Appendix A.2) and in this case we append the index τ to the τ ω χ2 − χ2 ω time dependent sums. The parameter ( ) is determined by the 0 ( ) σ ω p(ω) = (4) times ti and the measurement errors i for each frequency .So χ2 τ 0 when using as defined in Eq. (19) the periodogram in Eq. (5) 1 becomes p(ω) = SS · YC2 + CC · YS 2 − 2CS · YC · YS (5) YY · D YC2 YS 2 ω = 1 τ + τ · with: p( ) (20) YY CCτ SSτ D(ω) = CC · SS − CS 2 (6) Note that Eq. (20) has the same form as the Lomb-Scargle pe- and the following abbreviations for the sums: ff riodogram in Eq. (1) with the di erence that the errors can be weighted (weights wi in all sums) and that there is an additional Y = wiyi (7) second term in CCτ, SSτ, CS τ and tan 2ωτ (Eqs. (13)−(15) C = wi cos ωti (8) and (19), respectively) which accounts for the floating mean. The computational effort is similar as for the Lomb-Scargle = w ω S i sin ti (9) periodogram. The incorporation of the offset c requires only two ω = w ω additional sums for each frequency (namely S i sin ti = w ω ff YY = YYˆ − Y · Y YYˆ = w y2 (10) and C i cos ti or S τ and Cτ respectively). The e ort is i i even weaker when using Eq. (5) with keeping CS instead of us- τ YC(ω) = YCˆ − Y · C YCˆ = wiyi cos ωti (11) ing Eq. (20) with the parameter introduced via Eq. (19)which needs an extra preceding loop in the algorithm. If the errors are YS(ω) = YSˆ − Y · S YSˆ = wiyi sin ωti (12) taken into account as weights, also the multiplication with w i must be done. CC(ω) = CCˆ − C · C CCˆ = w cos2 ωt (13) i i For fast computation of the trigonometric sums the algorithm SS(ω) = SSˆ − S · S SSˆ = w sin2 ωt (14) of Press & Rybicki (1989) can be applied, which has advan- i i tages in the case of large data sets and/or many frequency steps. 2 CS(ω) = CSˆ − C · S CSˆ = wi cos ωti sin ωti. (15) Another possibility are trigonometric recurrences as described in Press et al. (1992). Note also that the first sum in SS can be Note that sums with hats correspond to the classical sums.
Load more

Recommended publications
  • Robper: an R Package to Calculate Periodograms for Light Curves Based on Robust Regression
    RobPer: An R Package to Calculate Periodograms for Light Curves Based on Robust Regression Anita M. Thieler Roland Fried Jonathan Rathjens TU Dortmund University TU Dortmund University TU Dortmund University Abstract An important task in astroparticle physics is the detection of periodicities in irregularly sampled time series, called light curves. The classic Fourier periodogram cannot deal with irregular sampling and with the measurement accuracies that are typically given for each observation of a light curve. Hence, methods to fit periodic functions using weighted regression were developed in the past to calculate periodograms. We present the R package RobPer which allows to combine different periodic functions and regression techniques to calculate periodograms. Possible regression techniques are least squares, least absolute deviations, least trimmed squares, M-, S- and τ-regression. Measurement accuracies can be taken into account including weights. Our periodogram function covers most of the approaches that have been tried earlier and provides new model-regression-combinations that have not been used before. To detect valid periods, RobPer applies an outlier search on the periodogram instead of using fixed critical values that are theoretically only justified in case of least squares regression, independent periodogram bars and a null hypothesis allowing only normal white noise. Finally, the package also includes a generator to generate artificial light curves. Keywords: periodogram, light curves, period detection, irregular sampling, robust regression, outlier detection, Cramér-von-Mises distance minimization, time series analysis, beta distri- bution, measurement accuracies, astroparticle physics, weighted regression, regression model. 1. Introduction We introduce the R (R Core Team 2015) package RobPer (Thieler, Rathjens, and Fried 2015), which can be used to calculate periodograms and detect periodicities in irregularly sampled time series.
    [Show full text]
  • Sigspec User's Manual
    Comm. in Asteroseismology - Complementary Topics Volume 1, 2010 c Austrian Academy of Sciences Preface Jeffrey D. Scargle1 1 Space Science and Astrobiology Division, NASA Ames Research Center SigSpec is a method for detecting and characterizing periodic signals in noisy data. This is an extremely common problem, not only in astronomy but in almost every branch of science and engineering. This work will be of great interest to anyone carrying out harmonic analysis employing Fourier techniques. The method is based on the definition of a quantity called spectral signif- icance – a function of Fourier phase and amplitude. Most data analysts are used to exploring only the Fourier amplitude, through the power spectrum, ignoring phase information. The Fourier phase spectrum can be estimated from data, but its interpretation is usually problematic. The spectral sig- nificance quantity conveys more information than does the conventional amplitude spectrum alone, and appears to simplify statistical issues as well as the interpretation of phase information. arXiv:1006.5081v1 [astro-ph.IM] 25 Jun 2010 Comm. in Asteroseismology - Complementary Topics Volume 1, 2010 c Austrian Academy of Sciences SigSpec User’s Manual P. Reegen1 1 Institut f¨ur Astronomie, T¨urkenschanzstraße 17, 1180 Vienna, Austria [email protected] Abstract SigSpec computes the spectral significance levels for the DFT amplitude spec- trum of a time series at arbitrarily given sampling. It is based on the analyti- cal solution for the Probability Density Function (PDF) of an amplitude level, including dependencies on frequency and phase and referring to white noise. Using a time series dataset as input, an iterative procedure including step- by-step prewhitening of the most significant signal components and MultiSine least-squares fitting is provided to determine a whole set of signal components, which makes the program a powerful tool for multi-frequency analysis.
    [Show full text]
  • A Heuristic Derivation of the Uncertainty of the Frequency Determination in Time Series Data
    Astronomy & Astrophysics manuscript no. Kallinger˙ferror c ESO 2018 November 1, 2018 A heuristic derivation of the uncertainty of the frequency determination in time series data Thomas Kallinger, Piet Reegen and Werner W. Weiss Institute for Astronomy (IfA), University of Vienna, T¨urkenschanzstrasse 17, A-1180 Vienna Accepted 22 December 2007 ABSTRACT Context. Several approaches to estimate frequency, phase and amplitude errors in time series analyses were reported in the literature, but they are either time consuming to compute, grossly overestimating the error, or are based on empirically determined criteria. Aims. A simple, but realistic estimate of the frequency uncertainty in time series analyses. Methods. Synthetic data sets with mono- and multi–periodic harmonic signals and with randomly distributed amplitude, frequency and phase were generated and white noise added. We tried to recover the input parameters with classical Fourier techniques and investigated the error as a function of the relative level of noise, signal and frequency difference. Results. We present simple formulas for the upper limit of the amplitude, frequency and phase uncertainties in time–series analyses. We also demonstrate the possibility to detect frequencies which are separated by less than the classical frequency resolution and that the realistic frequency error is at least 4 times smaller than the classical frequency resolution. Key words. methods: data analysis – methods: statistical 1. Motivation based on an analytical solution for the one–sigma error of a least–squares sinusoidal fit with a rms of σ(m). The total num- In the frequency analysis of time series, a realistic estimateof the ber of data points, the total time base of the observations, the amplitude, phase and frequency uncertainties can be of special signal amplitude, phase and frequency are N, T, a, φ and f , re- interest.
    [Show full text]
  • Sigspec User's Manual
    Comm. in Asteroseismology Volume 163, 2011 c Austrian Academy of Sciences SigSpec User’s Manual P. Reegen† Institut f¨ur Astronomie, T¨urkenschanzstraße 17, 1180 Vienna, Austria Abstract SigSpec computes the spectral significance levels for the DFT∗ amplitude spectrum of a time series at arbitrarily given sampling. It is based on the ana- lytical solution for the Probability Density Function (PDF) of an amplitude level, including dependencies on frequency and phase and referring to white noise. Us- ing a time series dataset as input, an iterative procedure including step-by-step prewhitening of the most significant signal components and MultiSine least- squares fitting is provided to determine a whole set of signal components, which makes the program a powerful tool for multi-frequency analysis. Instead of the step-by-step prewhitening of the most significant peaks, the program is also able to take into account several steps of the prewhitening sequence simultane- ously and check for the combination associated to a minimum residual scatter. This option is designed to overcome the aliasing problem caused by periodic time gaps in the dataset. SigSpec can detect non-sinusoidal periodicities in a dataset by simultaneously taking into account a fundamental frequency plus a set of harmonics. Time-resolved spectral significance analysis using a set of intervals of the time series is supported to investigate the development of eigen- frequencies over the observation time. Furthermore, an extension is available to perform the SigSpec analysis for multiple time series input files at once. In this MultiFile mode, time series may be tagged as target and comparison data.
    [Show full text]
  • Arxiv:Physics/0703160V1 [Physics.Data-An] 15 Mar 2007 Ok Nomto,Iitoueanwadubae Reliabili Unbiased and New a ( Introduce Criterion I Information, Book) W 4)
    Astronomy & Astrophysics manuscript no. AA˙2006˙6597 c ESO 2018 October 15, 2018 SigSpec I. Frequency- and Phase-Resolved Significance in Fourier Space P. Reegen1 Institut f¨ur Astronomie, Universit¨at Wien, T¨urkenschanzstraße 17, 1180 Vienna, Austria, e-mail: [email protected] Received October 19, 2006; accepted March 6, 2007 ABSTRACT Context. Identifying frequencies with low signal-to-noise ratios in time series of stellar photometry and spectroscopy, and measuring their amplitude ratios and peak widths accurately, are critical goals for asteroseismology. These are also challenges for time series with gaps or whose data are not sampled at a constant rate, even with modern Discrete Fourier Transform (DFT) software. Also the False-Alarm Probability introduced by Lomb and Scargle is an approximation which becomes less reliable in time series with longer data gaps. Aims. A rigorous statistical treatment of how to determine the significance of a peak in a DFT, called SigSpec, is presented here. SigSpec is based on an analytical solution of the probability that a DFT peak of a given amplitude does not arise from white noise in a non-equally spaced data set. Methods. The underlying Probability Density Function (PDF) of the amplitude spectrum generated by white noise can be derived explicitly if both frequency and phase are incorporated into the solution. In this paper, I define and evaluate an unbiased statistical estimator, the “spectral significance”, which depends on frequency, amplitude, and phase in the DFT, and which takes into account the time-domain sampling. Results. I also compare this estimator to results from other well established techniques and assess the advantages of SigSpec, through comparison of its analytical solutions to the results of extensive numerical calculations.
    [Show full text]
  • Spectral Analysis of Nonuniformly Sampled Data
    IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 57, NO. 3, MARCH 2009 843 Spectral Analysis of Nonuniformly Sampled Data: A New Approach Versus the Periodogram Petre Stoica, Fellow, IEEE, Jian Li, Fellow, IEEE, and Hao He, Student Member, IEEE Abstract—We begin by revisiting the periodogram to explain B. The Fourier Periodogram why arguably the plain least-squares periodogram (LSP) is prefer- able to the “classical” Fourier periodogram, from a data-fitting The “classical” Fourier transform-based periodogram (FP) viewpoint, as well as to the frequently-used form of LSP due to associated with is given by Lomb and Scargle, from a computational standpoint. Then we go on to introduce a new enhanced method for spectral analysis of nonuniformly sampled data sequences. The new method can (1) be interpreted as an iteratively weighted LSP that makes use of a data-dependent weighting matrix built from the most recent spectral estimate. Because this method is derived for the case where is the frequency variable and where, depending on the of real-valued data (which is typically more complicated to deal with in spectral analysis than the complex-valued data case), it application, the normalization factor might be different from is iterative and it makes use of an adaptive (i.e., data-dependent) (such as , see, e.g., [1] and [2]). It can be readily weighting, we refer to it as the real-valued iterative adaptive verified that can be obtained from the solution to the fol- approach (RIAA). LSP and RIAA are nonparametric methods lowing least-squares (LS) data fitting problem: that can be used for the spectral analysis of general data sequences with both continuous and discrete spectra.
    [Show full text]
  • An R Package to Calculate Periodograms for Light Curves Based on Robust Regression
    JSS Journal of Statistical Software March 2016, Volume 69, Issue 9. doi: 10.18637/jss.v069.i09 RobPer: An R Package to Calculate Periodograms for Light Curves Based on Robust Regression Anita M. Thieler Roland Fried Jonathan Rathjens TU Dortmund University TU Dortmund University TU Dortmund University Abstract An important task in astroparticle physics is the detection of periodicities in irregularly sampled time series, called light curves. The classic Fourier periodogram cannot deal with irregular sampling and with the measurement accuracies that are typically given for each observation of a light curve. Hence, methods to fit periodic functions using weighted regression were developed in the past to calculate periodograms. We present the R package RobPer which allows to combine different periodic functions and regression techniques to calculate periodograms. Possible regression techniques are least squares, least absolute deviations, least trimmed squares, M-, S- and τ-regression. Measurement accuracies can be taken into account including weights. Our periodogram function covers most of the approaches that have been tried earlier and provides new model-regression-combinations that have not been used before. To detect valid periods, RobPer applies an outlier search on the periodogram instead of using fixed critical values that are theoretically only justified in case of least squares regression, independent periodogram bars and a null hypothesis allowing only normal white noise. Finally, the package also includes a generator to generate artificial light curves. Keywords: periodogram, light curves, period detection, irregular sampling, robust regression, outlier detection, Cramér-von-Mises distance minimization, time series analysis, beta distri- bution, measurement accuracies, astroparticle physics, weighted regression, regression model.
    [Show full text]
  • Abstracts of Talks (A) Sun and Solar System
    29th ASI Meeting ASI Conference Series, 2011, Vol. 3, pp 99–189 Edited by Pushpa Khare & C. H. Ishwara-Chandra Abstracts of talks (A) Sun and solar system Modeling solar cycles using a variable meridional circulation in a flux transport dynamo model Bidya Binay Karak and Arnab Rai Choudhuri Department of Physics, Indian Institute of Science, Bangalore E-mail:bidya [email protected] Abstract. The sunspot number − a proxy of solar activity − varies roughly periodically with time. However the individual cycle durations and amplitudes are found to vary in an irregular way. An important feature of the solar cycle is the Maunder minimum during 1645-1715 when there were very few sunspots. We explore whether this irregular solar cycle can be modeled with the help of a flux transport dynamo model of the solar cycle. We model the periods of the last 23 sunspot cycles by varying the meridional circulation speed. We find that most of the cycle amplitudes also get modeled up to some extent when we model the periods. Moreover, under certain situations we are able to re- produce Maunder like grand minimum. However, we fail to reproduce these results if the value of turbulent diffusivity is reasonably low. Study of polarimetric properties of comet Levy 1990XX by mixture of compact and aggregate particles A. Suklabaidya, D. Paul , H. S. Das and A. K. Sen Department of Physics, Assam University, Silchar, Assam E-mail: abinash 80@rediffmail.com Abstract. In the present work, the observed linear polarization data of comet Levy 1990XX are studied at wavelengths 0.485 µm, 0.540 µm 100 Abstracts of talks and 0.670 µm.
    [Show full text]
  • Period Analysis of Stellar Brightness Variations Using Multichannel
    PERIOD ANALYSIS OF STELLAR BRIGHTNESS VARIATIONS USING MULTICHANNEL PHOTOMETRIC OBSERVATIONS JIANPENG YOU Division of Astronomy Department of Physical Sciences University of Oulu Finland Academic dissertation to be presented, with the permission of the Faculty of Science of the University of Oulu, for public discussion in the Auditorium Anttilansali, Linnanmaa, on 14 December, 2007, at 12 o’clock noon. REPORT SERIES IN PHYSICAL SCIENCES Report No. 45 OULU 2007 UNIVERSITY OF OULU Reviewers Prof. Werner W. Weiss, University of Vienna, Austria Prof. Nikolai Piskunov, University of Uppsala, Sweden Opponent Doc. Harry J. Lehto, University of Turku, Finland Custos Prof. Juri Poutanen, University of Oulu, Finland ISBN 978-951-42-8642-1 ISBN 978-951-42-8643-8 (PDF) ISSN 1239-4327 Oulu University Press Oulu 2007 Abstract This dissertation presents time series analysis (TSA) on the simultaneous multichannel photometric observations to study stellar differential rotation with a novel three stage, multichannel period analysis (MPA). The traditional TSA methods are typically applied only to one channel photometric observations. This PhD work devises a new MPA method, based on the traditional TSA technique, to analyze simultaneously all channel data together. The applications of MPA focus on analyzing the photometric obser- vations of brightness variations of stars, studying stellar surface spottedness and differential rotation. Similarly to the solar differential rotation, which can be studied using sunspots, stellar differential rotation can be investigated from photometric variations caused by surface spots, which can be observed instrumentally and analyzed by different TSA techniques. The PhD work shows that the MPA method has advantages both in recovering physical periodicities and in achieving higher accuracies of the final results over the traditional TSA applied to the single channel data.
    [Show full text]
  • Signature Redacted
    Asteroseismic Study of the Subgiant HD 82074 A__ _M! MASSACHUSETTS INSTfl1JE- by OF TECHNOLOGY Victoria Ashley Villar [j 5204 Submitted to the Department of Physics LIBRARIES in partial fulfillment of the requirements for the degree of Bachelor of Science in Physics at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY June 2014 @ Massachusetts Institute of Technology 2014. All rights reserved. Signature redacted Author. Department of Physics Signature redacted May 9, 2014 Certified by.. .. .. .... .. .. .... ... .. .. .. ... John A. Johnson Harvard Department of Astronomy redacted Thesis Supervisor Certified by.... Signature Joshua N. Winn Department of Physics Thesis Supervisor Signature redacted Accepted by...... ................... Nergis Mavalvala Senior Thesis Coordinator, Department of Physics 2 Asteroseismic and Interferometric Study of the Subgiant HD 82074 by Victoria Ashley Villar Submitted to the Department of Physics on May 9, 2014, in partial fulfillment of the requirements for the degree of Bachelor of Science in Physics Abstract This thesis analyzes HD 82074, a solar-mass and low-metallicity subgiant star, using ground-based asteroseismology with the CHIRON spectrometer at the Cerro Tololo Inter-American Observatory (CTIO) and interferometry with the CHARA array on Mt. Wilson. The physical parameters of subgiant stars are of particular interest in the exoplanetary field due to their importance in understanding the relationships between planet occurrence rate and stellar properties such as age, metallicity and mass. Potential systematic uncertainties in the canonical stellar models make it especially important to independently determine the masses and radii of these stars. We determine HD 82074's radius using interferometry from CHARA, and we combine this result with measurements of the spacing and frequencies of the asteroseismic oscillations of HD 82074 to determine a stellar mass.
    [Show full text]
  • Communications in Asteroseismology
    Communications in Asteroseismology Volume 163 2011 Communications in Asteroseismology Editor-in-Chief: Michel Breger, [email protected] Editorial Assistant: Isolde M¨uller, [email protected] Layout & Production Manager: Isolde M¨uller, [email protected] CoAst Editorial and Production Office T¨urkenschanzstraße 17, A - 1180 Wien, Austria http://www.oeaw.ac.at/CoAst/ [email protected] Editorial Board: Conny Aerts, Gerald Handler, Don Kurtz, Jaymie Matthews, Ennio Poretti Cover Illustration Sample Phase Distribution Diagram (For more information see page 36.) British Library Cataloguing in Publication data. A Catalogue record for this book is available from the British Library. All rights reserved ISBN 978-3-7001-7194-2 ISSN 1021-2043 Copyright c 2011 by Austrian Academy of Sciences Vienna Austrian Academy of Sciences Press A-1011 Wien, Postfach 471, Postgasse 7/4 Tel. +43-1-515 81/DW 3402-3406, +43-1-512 9050 Fax +43-1-515 81/DW 3400 http://verlag.oeaw.ac.at, e-mail: [email protected] Preface by J. D. Scargle 1 SigSpec User’s Manual by P. Reegen 3 1. What is SigSpec? ........................ 4 2. How to Run SigSpec ....................... 7 2.1. Projects . 7 2.2. Quietmode ........................ 12 3. Input................................ 12 3.1. Thetimeseriesinputfile . 12 3.2. The .ini file ....................... 12 3.3. Time series columns representing time and observable . 13 3.4. Time series columns containing statistical weights . 13 3.5. Time series columns containing subset identifiers . 15 3.6. Lowerfrequencylimit . 18 3.7. Upper frequency limit and Nyquist Coefficient . 19 3.8. Frequencyspacing and oversamplingratio .
    [Show full text]
  • THE NATURE of P-MODES and GRANULATION in PROCYON: NEW MOST1 PHOTOMETRY and NEW YALE CONVECTION MODELS D
    The Astrophysical Journal, 687:1448Y1459, 2008 November 10 # 2008. The American Astronomical Society. All rights reserved. Printed in U.S.A. THE NATURE OF p-MODES AND GRANULATION IN PROCYON: NEW MOST1 PHOTOMETRY AND NEW YALE CONVECTION MODELS D. B. Guenther Institute for Computational Astrophysics, Department of Astronomy and Physics, Saint Mary’s University, Halifax, NS B3H 3C3, Canada T. Kallinger, M. Gruberbauer, D. Huber, W. W. Weiss, R. Kuschnig Institut fu¨r Astronomie, Universita¨t Wien Tu¨rkenschanzstrasse 17, A-1180 Wien, Austria P. Demarque, F. Robinson Department of Astronomy, Yale University, New Haven CT 06511 J. M. Matthews Department of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T 1Z1, Canada A. F. J. Moffat Observatoire Astronomque du Mont Me´gantic, De´partment de Physique, Universite´ de Montre´al C. P. 6128, Succursale: Centre-Ville, Montre´al, QC H3C 3J7, Canada S. M. Rucinski Department of Astronomy and Astrophysics, University of Toronto, ON M5S 3H4, Canada D. Sasselov Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 102138 and G. A. H. Walker Department of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T 1Z1, Canada Received 2008 February 12; accepted 2008 June 30 ABSTRACT We present new photometry of Procyon, obtained by MOST during a 38 day run in 2007, and frequency analyses of those data. The long time coverage and low point-to-point scatter of the light curve yield an average noise amplitude of about 1.5Y2.0 ppm in the frequency range 500Y1500 Hz. This is half the noise level obtained from each of the previous two Procyon campaigns by MOST in 2004 and 2005.
    [Show full text]