Engineering Papers⌕ Search

Engineering topics

Mighell, Kenneth J.

Publications and source records attributed to Mighell, Kenneth J..

Kepler and the long-period variables

High-precision Kepler photometry is used to explore the details of asymptotic giant branch (AGB) light curves. Since AGB variability has a typical timescale on the order of a year, we discuss at length the removal of long-term trends and quarterly changes in Kepler data. Photometry for a small sample of nine semi-regular (SR) AGB stars is examined using a 30 minute cadence over a period of 45 months. While undergoing long-period variations of many magnitudes, the light curves are shown to be smooth at the millimagnitude level over much shorter time intervals. No flares or other rapid events were detected on a sub-day timescale. The shortest AGB period detected is on the order of 100 days. All the SR variables in our sample are shown to have multiple modes. This is always the first overtone, typically combined with the fundamental. A second common characteristic of SR variables is shown to be the simultaneous excitation of multiple closely separated periods for the same overtone mode. Approximately half the sample had a much longer variation in the light curve, likely a long secondary period (LSP). The light curves were all well represented by a combination of sinusoids. However, the properties of the sinusoids are time variable, with irregular variations present at low levels. No non-radial pulsations were detected. It is argued that the LSP variation seen in many SR variables is intrinsic to the star and linked to multiple mode pulsation.

79 ASTRONOMY AND ASTROPHYSICS↗

The variable sky of deep synoptic surveys

The discovery of variable and transient sources is an essential product of synoptic surveys. The alert stream will require filtering for personalized criteria—a process managed by a functionality commonly described as a Broker. In order to understand quantitatively the magnitude of the alert generation and Broker tasks, we have undertaken an analysis of the most numerous types of variable targets in the sky—Galactic stars, quasi-stellar objects (QSOs), active galactic nuclei (AGNs), and asteroids. It is found that the Large Synoptic Survey Telescope (LSST) will be capable of discovering ∼10{sup 5} high latitude (|b| > 20°) variable stars per night at the beginning of the survey. (The corresponding number for |b| < 20° is orders of magnitude larger, but subject to caveats concerning extinction and crowding.) However, the number of new discoveries may well drop below 100 per night within less than one year. The same analysis applied to GAIA clarifies the complementarity of the GAIA and LSST surveys. Discovery of AGNs and QSOs are each predicted to begin at ∼3000 per night and decrease by 50 times over four years. Supernovae are expected at ∼1100 per night, and after several survey years will dominate the new variable discovery rate. LSST asteroid discoveries will start at >10{sup 5} per night, and if orbital determination has a 50% success rate per epoch, they will drop below 1000 per night within two years.

79 ASTRONOMY AND ASTROPHYSICS↗

PERIOD ERROR ESTIMATION FOR THE KEPLER ECLIPSING BINARY CATALOG

The Kepler Eclipsing Binary Catalog (KEBC) describes 2165 eclipsing binaries identified in the 115 deg{sup 2} Kepler Field based on observations from Kepler quarters Q0, Q1, and Q2. The periods in the KEBC are given in units of days out to six decimal places but no period errors are provided. We present the PEC (Period Error Calculator) algorithm, which can be used to estimate the period errors of strictly periodic variables observed by the Kepler Mission. The PEC algorithm is based on propagation of error theory and assumes that observation of every light curve peak/minimum in a long time-series observation can be unambiguously identified. The PEC algorithm can be efficiently programmed using just a few lines of C computer language code. The PEC algorithm was used to develop a simple model that provides period error estimates for eclipsing binaries in the KEBC with periods less than 62.5 days: log {sigma}{sub P} Almost-Equal-To - 5.8908 + 1.4425(1 + log P), where P is the period of an eclipsing binary in the KEBC in units of days. KEBC systems with periods {>=}62.5 days have KEBC period errors of {approx}0.0144 days. Periods and period errors of seven eclipsing binary systems in the KEBC were measured using the NASA Exoplanet Archive Periodogram Service and compared to period errors estimated using the PEC algorithm.

79 ASTRONOMY AND ASTROPHYSICS↗

ORBITAL PHASE VARIATIONS OF THE ECCENTRIC GIANT PLANET HAT-P-2b

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 {mu}m bands of the Spitzer Space Telescope. The 3.6 and 4.5 {mu}m data sets span an entire orbital period of HAT-P-2b (P = 5.6334729 d), making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 {mu}m that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39 {+-} 0.28, 5.84 {+-} 0.39, and 4.68 {+-} 0.37 hr after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138% {+-} 0.0089%, 0.1162% {+-} 0.0080%, and 0.1888% {+-} 0.0072% in the 3.6, 4.5, and 8.0 {mu}m bands, respectively. Our measured secondary eclipse depths of 0.0996% {+-} 0.0072%, 0.1031% {+-} 0.0061%, 0.071%{sub -0.013%}{sup +0.029,} and 0.1392% {+-} 0.0095% in the 3.6, 4.5, 5.8, and 8.0 {mu}m bands, respectively, indicate that the planet cools significantly from its peak temperature before we measure the dayside flux during secondary eclipse. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity (e = 0.50910 {+-} 0.00048) and argument of periapse ({omega} = 188. Degree-Sign 09 {+-} 0. Degree-Sign 39) of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

Quick and Dirty WFPC2 Stellar Photometry

The latest release of the IRAF (Image Reduction and Analysis Facility) MXTOOLS package includes the new tasks QDWFPC2 (Quick and Dirty Wide Field and Planetary Camera, which does quick CCD (charge-coupled device) stellar photometry on two Hubble Space Telescope (HST) WFPC2 (Wide Field and Planetary Camera) observations: WFPC2COLOR, which converts HST WFPC2 instrumental magnitudes to standard colors using the Holtzman et al. color equations, and QDCMD (Quick and Dirty Color Magnitude Diagram), which reads the output of WFPC2COLOR and displays a color-magnitude diagram on a user-chosen graphics device.

Mighell, Kenneth J.↗

Astrometry of the omega Centauri Hubble Space Telescope Calibration Field

Astrometry, on the International Celestial Reference Frame (epoch J2000.0), is presented for the Walker (1994, PASP, 106, 828) stars in the omega Centauri (=NGC 5139=C 1323-1472) Hubble Space Telescope Wide Field/Planetary Camera (WF/PC) calibration field of Harris et al. (1993, AJ, 105, 1196). Harris et al. stars were first identified on a WFPC2 observation of the omega Cen HST calibration field. Relative astrometry of the Walker stars in this field was then obtained using Walker's CCD positions and astrometry derived using the STSDAS METRIC task on the positions of the Harris et al. stars on the WFPC2 observation. Finally, the relative astrometry, which was based on the HST Guide Star Catalog, is placed on the International Celestial Reference Frame with astrometry from the USNO-A2.0 catalog. An ASCII text version of the astrometric data of the Walker stars in the omega Cen HST calibration field is available electronically in the online version of the article.

Mighell, Kenneth J.↗

WFPC2 Observations of the URSA Minor Dwarf Spheroidal Galaxy

We present our analysis of archival Hubble Space Telescope Wide Field Planetary Camera 2 (WFPC2) observations in F555W (approximately V) and F814W (approximately I) of the central region of the Ursa Minor dwarf spheroidal galaxy. The V versus V - I color-magnitude diagram features a sparsely populated blue horizontal branch, a steep thin red giant branch, and a narrow subgiant branch. The main sequence reaches approximately 2 magnitudes below the main-sequence turnoff (V(sup UMi, sub TO) approximately equals 23.27 +/- 0.11 mag) of the median stellar population. We compare the fiducial sequence of the Galactic globular cluster M92 (NGC 6341). The excellent match between Ursa Minor and M92 confirms that the median stellar population of the UMi dSph galaxy is metal poor ([Fe/H](sub UMi) approximately equals [Fe/H](sub M92) approximately equals -2.2 dex) and ancient (age(sub UMi)approximately equalsage(sub M92) approximately equals 14 Gyr). The B - V reddening and the absorption in V are estimated to be E(B - V) = 0.03 +/- 0.01 mag and A(sup UMi, sub V) = 0.09 +/- 0.03 mag. A new estimate of the distance modulus of Ursa Minor, (m - M)(sup UMi, sub 0) = 19.18 +/- 0.12 mag, has been derived based on fiducial-sequence fitting M92 [DELTA.V(sub UMi - M92) = 4.60 +/- 0.03 mag and DELTA(V - I)(sub UMi - M92) = 0.010 +/- 0.005 mag] and the adoption of the apparent V distance modulus for M92 of (m - M)(sup M92, sub V) = 14.67 +/- 0.08 mag (Pont et al. 1998, A&A, 329, 87). The Ursa Minor dwarf spheroidal galaxy is then at a distance of 69 +/- 4 kpc from the Sun. These HST observations indicate that Ursa Minor has had a very simple star formation history consisting mainly of a single major burst of star formation about 14 Gyr ago which lasted approximately < 2 Gyr. While we may have missed minor younger stellar populations due to the small field-of-view of the WFPC2 instrument, these observations clearly show that most of the stars in the central region Ursa Minor dwarf spheroidal galaxy are ancient. If the ancient Galactic globular clusters, like M92, formed concurrently with the early formation of the Milky Way galaxy itself, then the Ursa Minor dwarf spheroidal is probably as old as the Milky Way.

Mighell, Kenneth J.↗

Parameter Estimation in Astronomy with Poisson-Distributed Data: The (CHI)2(gamma) Statistic - 1

Applying the standard weighted mean formula, [Sigma (sub i)n(sub i)ssigma(sub i, sup -2)], to determine the weighted mean of data, n(sub i), drawn from a Poisson distribution, will, on average, underestimate the true mean by approx. 1 for all true mean values larger than approx.3 when the common assumption is made that the error of the i th observation is sigma(sub i) = max square root of n(sub i), 1).This small, but statistically significant offset, explains the long-known observation that chi-square minimization techniques which use the modified Neyman'chi(sub 2) statistic, chi(sup 2, sub N) equivalent Sigma(sub i)((n(sub i) - y(sub i)(exp 2)) / max(n(sub i), 1), to compare Poisson - distributed data with model values, y(sub i), will typically predict a total number of counts that underestimates the true total by about 1 count per bin. Based on my finding that weighted mean of data drawn from a Poisson distribution can be determined using the formula [Sigma(sub i)[n(sub i) + min(n(sub i), 1)](n(sub i) + 1)(exp -1)] / [Sigma(sub i)(n(sub i) + 1)(exp -1))], I propose that a new chi(sub 2) statistic, chi(sup 2, sub gamma) equivalent, should always be used to analyze Poisson- distributed data in preference to the modified Neyman's chi(exp 2) statistic. I demonstrated the power and usefulness of,chi(sub gamma, sup 2) minimization by using two statistical fitting techniques and five chi(exp 2) statistics to analyze simulated X-ray power - low 15 - channel spectra with large and small counts per bin. I show that chi(sub gamma, sup 2) minimization with the Levenberg - Marquardt or Powell's method can produce excellent results (mean slope errors approx. less than 3%) with spectra having as few as 25 total counts.

Mighell, Kenneth J.↗

WFPC2 Observations of Star Clusters in the Magellanic Clouds: The Oldest Star Clusters in the Small Magellanic Cloud - Report 2

We present our analysis of archival Hubble Space Telescope Wide Field Planetary Camera 2 (WFPC2) observations in F45OW ( approximately B) and F555W (approximately V) of the intermediate-age populous star clusters NGC 121, NGC 339, NGC 361, NGC 416, and Kron 3 in the Small Magellanic Cloud. We use published photometry of two other SMC populous star clusters, Lindsay 1 and Lindsay 113, to investigate the age sequence of these seven populous star clusters in order to improve our understanding of the formation chronology of the SMC. We analyzed the V vs B-V and M(sub V) vs (B-V)(sub 0) color-magnitude diagrams of these populous Small Magellanic Cloud star clusters using a variety of techniques and determined their ages, metallicities, and reddenings. These new data enable us to improve the age-metallicity relation of star clusters in the Small Magellanic Cloud. In particular, we find that a closed-box continuous star-formation model does not reproduce the age-metallicity relation adequately. However, a theoretical model punctuated by bursts of star formation is in better agreement with the observational data presented herein.

Mighell, Kenneth J.↗