Cepheid Research Articles

The Population of Pulsating Variable Stars in the Sextans Dwarf Spheroidal Galaxy

Abstract A large extension of the Sextans dwarf spheroidal galaxy, 7 deg2, has been surveyed for variable stars using the Dark Energy Camera at the Blanco Telescope at Cerro Tololo Inter-American Observatory, Chile. We report seven anomalous Cepheids, 199 RR Lyrae stars, and 16 dwarf Cepheids in the field. This is only the fifth extragalactic system in which dwarf Cepheids have been systematically searched. Henceforth, the new stars increase the census of stars coming from different environments that can be used to asses the advantages and limitations of using dwarf Cepheids as standard candles in populations for which the metallicity is not necessarily known. The dwarf Cepheids found in Sextans have a mean period of 0.066 day and a mean g amplitude of 0.87 mag. They are located below the horizontal branch, spanning a range of 0.8 mag: 21.9 < g < 22.7. The number of dwarf Cepheids in Sextans is low compared with other galaxies such as Carina, which has a strong intermediate-age population. On the other hand, the number and ratio of RR Lyrae stars to dwarf Cepheids are quite similar to those of Sculptor, a galaxy which, as Sextans, is dominated by an old stellar population. The dwarf Cepheid stars found in Sextans follow a well-constrained period–luminosity relationship with an rms = 0.05 mag in the g band, which was set up by anchoring to the distance modulus given by the RR Lyrae stars. Although the majority of the variable stars in Sextans are located toward the center of the galaxy, we have found two RR Lyrae stars and one anomalous Cepheid in the outskirts of the galaxy that may be extratidal stars and suggest that this galaxy may be undergoing tidal destruction. These possible extratidal variable stars share the same proper motions as Sextans, as seen by recent Gaia measurements. Two additional stars that we initially classified as foreground RR Lyrae stars may actually be other examples of Sextans extratidal anomalous Cepheids, although radial velocities are needed to prove that scenario.

Light Curve Parameters of Cepheid and RR Lyrae Variables at Multiple Wavelengths − Models vs. Observations

We present results from a comparative study of light curves of Cepheid and RR Lyrae stars in the Galaxy and the Magellanic Clouds with their theoretical models generated from the stellar pulsation codes. Fourier decomposition method is used to analyse the theoretical and the observed light curves at multiple wavelengths. In case of RR Lyrae stars, the amplitude and Fourier parameters from the models are consistent with observations in most period bins except for low metal-abundances (Z < 0.004). In case of Cepheid variables, we observe a greater offset between models and observations for both the amplitude and Fourier parameters. The theoretical amplitude parameters are typically larger than those from observations, except close to the period of 10 days. We find that these discrepancies between models and observations can be reduced if a higher convective efficiency is adopted in the pulsation codes. Our results suggest that a quantitative comparison of light curve structure is very useful to provide constraints for the input physics to the stellar pulsation models.

The Carnegie Supernova Project: Absolute Calibration and the Hubble Constant

Abstract We present an analysis of the final data release of the Carnegie Supernova Project I, focusing on the absolute calibration of the luminosity–decline rate relation for Type Ia supernovae (SNe Ia) using new intrinsic color relations with respect to the color-stretch parameter, s BV , enabling improved dust extinction corrections. We investigate to what degree the so-called fast-declining SNe Ia can be used to determine accurate extragalactic distances. We estimate the intrinsic scatter in the luminosity–decline rate relation and find it ranges from ±0.13 mag to ±0.18 mag with no obvious dependence on wavelength. Using the Cepheid variable star data from the SH0ES project, the SN Ia distance scale is calibrated and the Hubble constant is estimated using our optical and near-infrared sample, and these results are compared to those determined exclusively from a near-infrared subsample. The systematic effect of the supernova’s host galaxy mass is investigated as a function of wavelength and is found to decrease toward redder wavelengths, suggesting this effect may be due to dust properties of the host. Using estimates of the dust extinction derived from optical and near-infrared wavelengths and applying these to the H band, we derive a Hubble constant , whereas using a simple B − V color correction applied to the B band yields . Photometry of two calibrating SNe Ia from the CSP-II sample, SN 2012ht and SN 2015F, is presented and used to improve the calibration of the SN Ia distance ladder.

Gaia Data Release 2

Context. The second data release of the Gaia mission (DR2) includes an advance catalogue of variable stars. The classifications of these stars are based on sparse photometry from the first 22 months of the mission. Aims. We set out to investigate the purity and completeness of the all-sky Gaia classification results with the help of the continuous light curves of the observed targets from the Kepler and K2 missions, focusing specifically on RR Lyrae and Cepheid pulsators, outside the Galactic bulge region. Methods. We cross-matched the Gaia identifications with the observations collected by the Kepler space telescope. We inspected the light curves visually, then calculated the relative Fourier coefficients and period ratios for the single- and double-mode K2 RR Lyrae stars to further classify them. Results. We identified 1443 and 41 stars classified as RR Lyrae or Cepheid variables in Gaia DR2 in the targeted observations of the two missions and 263 more RR Lyre targets in the full-frame images (FFI) of the original mission. We provide the cross-match of these sources. We conclude that the RR Lyrae catalogue has a completeness between 70–78%, and provide a purity estimate of between 92 and 98% (targeted observations) with lower limits of 75% (FFI stars) and 51% (K2 worst-case scenario). The low number of Cepheids prevents us from drawing detailed conclusions, but the purity of the DR2 sample is estimated to be about 66%.

Variability in Proto-planetary Nebulae. V. Velocity and Light Curve Analysis of IRAS 17436+5003, 18095+2704, and 19475+3119

Abstract We have obtained contemporaneous light, color, and radial velocity data for three proto-planetary nebulae (PPNe) over the years 2007 to 2015. The light and velocity curves of each show similar periods of pulsation, with photometric periods of 42 and 50 days for IRAS 17436+5003, 102 days for IRAS 18095+2704, and 35 days for IRAS 19475+3119. The light and velocity curves are complex with multiple periods and small, variable amplitudes. Nevertheless, at least over limited time intervals, we were able to identify dominant periods in the light, color, and velocity curves and compare the phasing of each. The color curves appear to peak with or slightly after the light curves while the radial velocity curves peak about a quarter of a cycle before the light curves. Similar results were found previously for two other PPNe, although for them the light and color appeared to be in phase. Thus, it appears that PPNe are brightest when smallest and hottest. These phase results differ from those found for classical Cepheid variables, where the light and velocity differ by half a cycle, and are hottest at about average size and expanding. However, they do appear to have similar phasing to the larger-amplitude pulsations seen in RV Tauri variables. Presently, few pulsation models exist for PPNe, and these do not fit the observations well, especially the longer periods observed. Model fits to these new light and velocity curves would allow masses to be determined for these post-AGB objects, and thereby provide important constraints to post-AGB stellar evolution models of low- and intermediate-mass stars.

Old-Aged Primary Distance Indicators

Old-aged stellar distance indicators are present in all Galactic structures (halo, bulge, disk) and in galaxies of all Hubble types and, thus, are immensely powerful tools for understanding our Universe. Here we present a comprehensive review for three primary standard candles from Population II: (i) RR Lyrae type variables (RRL), (ii) type II Cepheid variables (T2C), and (iii) the tip of the red giant branch (TRGB). The discovery and use of these distance indicators is placed in historical context before describing their theoretical foundations and demonstrating their observational applications across multiple wavelengths. The methods used to establish the absolute scale for each standard candle is described with a discussion of the observational systematics. We conclude by looking forward to the suite of new observational facilities anticipated over the next decade; these have both a broader wavelength coverage and larger apertures than current facilities. We anticipate future advancements in our theoretical understanding and observational application of these stellar populations as they apply to the Galactic and extragalactic distance scale.

Trapping of plasma enabled by pycnoclinic acoustic force

Sound can hold partially ionized sulfur at the center of a spherical bulb. We use the sulfur plasma itself to drive a 180 dB re 20 $\ensuremath{\mu}\mathrm{Pa}$ sound wave by periodically heating it with microwave pulses at a frequency that matches the lowest order, spherically symmetric, acoustic resonance of the bulb. To clarify the trapping mechanism, we generalize acoustic radiation pressure theory to include gaseous inhomogeneities and find an interaction of high-amplitude sound with density gradients in the gas through which it propagates. This is the pycnoclinic acoustic force (PAF). Though generated by rapidly oscillating sound waves, it has a finite time average and manipulates the plasma through density gradients at its boundary. The PAF is essential for the description of the trap holding a plasma against its own buoyancy as well as understanding convection in the region outside the plasma. It has implications for pulse tubes, thermoacoustic engines, thermal vibrational convection in microgravity, combustion in the presence of sound, and the modeling of Cepheid variable stars.

Photometric study of the SMCNOD using variable stars from the OGLE-IV survey

We present a study of a recently discovered stellar overdensity near the northern edge of the Small Magellanic Cloud (SMCNOD). We exploited variable stars from the fourth release of the Optical Gravitational Lensing Experiment (OGLE). We used mainly pulsating variable stars and investigated their potential association with the SMCNOD using their spatial distribution and distances. We found four rather spatially concentrated anomalous Cepheids and eight evenly dispersed RR~Lyrae stars to be most likely members of this overdensity. The anomalous Cepheids inside the SMCNOD trace possible intermediate-age population with ages ranging between $2-4.5$ Gyr. The age distribution of anomalous Cepheids seems to be in a good agreement with the age distribution of anomalous Cepheids in the SMC. Using empirical relations for RR~Lyrae stars we determined the median metallicity for a possible old population in the SMCNOD to be $\rm [Fe/H]_{\rm SMCNOD}=-1.71\pm0.21$ dex, which is in agreement with median metallicity of the old SMC population. The density profile for anomalous Cepheids shows a small anomaly at the position of the SMCNOD, on the other hand, RR~Lyrae variables show no such deviation. The probability of finding the observed number of variable stars at the location of the SMCNOD by chance is very low for anomalous Cepheids (0.7 %) but high for RR~Lyrae stars (13.0 %). Based on its variable stars content, we thus confirm the presence of a modest overdensity in intermediate-age stars in the SMCNOD and conclude that it probably has its origin in the SMC rather than to be the remnant of an accreted dwarf galaxy.

Astrophysical tests of modified gravity: Stellar and gaseous rotation curves in dwarf galaxies

Chameleon theories of gravity predict that the gaseous component of isolated dwarf galaxies rotates with a faster velocity than the stellar component. In this paper, we exploit this effect to obtain new constraints on the model parameters using the measured rotation curves of six low surface brightness galaxies. For $f(R)$ theories, we rule out values of $f_{R0}>10^{-6}$. For more general theories, we find that the constraints from Cepheid variable stars are currently more competitive than the bounds we obtain here but we are able to rule out self-screening parameters $\chi_c>10^{-6}$ for fifth-force strengths (coupling of the scalar to matter) as low as $0.05$ the Newtonian force. This region of parameter space has hitherto been inaccessible to astrophysical probes. We discuss the future prospects for improving these bounds.

New Cepheid variables in the young open clusters Berkeley 51 and Berkeley 55

As part of a wider investigation of evolved massive stars in Galactic open clusters, we have spectroscopically identified three candidate classical Cepheids in the little-studied clusters Berkeley 51, Berkeley 55 and NGC 6603. Using new multi-epoch photometry, we confirm that Be 51 #162 and Be 55 #107 are bona fide Cepheids, with pulsation periods of 9.83+/-0.01 d and 5.850+/-0.005 d respectively, while NGC 6603 star W2249 does not show significant photometric variability. Using the period-luminosity relationship for Cepheid variables, we determine a distance to Be 51 of 5.3(+1.0,-0.8) kpc and an age of 44(+9,-8) Myr, placing it in a sparsely-attested region of the Perseus arm. For Be 55, we find a distance of 2.2+/-0.3 kpc and age of 63(+12,-11) Myr, locating the cluster in the Local arm. Taken together with our recent discovery of a long-period Cepheid in the starburst cluster VdBH222, these represent an important increase in the number of young, massive Cepheids known in Galactic open clusters. We also consider new Gaia (data release 2) parallaxes and proper motions for members of Be 51 and Be 55; the uncertainties on the parallaxes do not allow us to refine our distance estimates to these clusters, but the well-constrained proper motion measurements furnish further confirmation of cluster membership. However, future final Gaia parallaxes for such objects should provide valuable independent distance measurements, improving the calibration of the period-luminosity relationship, with implications for the distance ladder out to cosmological scales.

Robust Period Estimation Using Mutual Information for Multiband Light Curves in the Synoptic Survey Era

The Large Synoptic Survey Telescope (LSST) will produce an unprecedented amount of light curves using six optical bands. Robust and efficient methods that can aggregate data from multidimensional sparsely-sampled time series are needed. In this paper we present a new method for light curve period estimation based on the quadratic mutual information (QMI). The proposed method does not assume a particular model for the light curve nor its underlying probability density and it is robust to non-Gaussian noise and outliers. By combining the QMI from several bands the true period can be estimated even when no single-band QMI yields the period. Period recovery performance as a function of average magnitude and sample size is measured using 30,000 synthetic multi-band light curves of RR Lyrae and Cepheid variables generated by the LSST Operations and Catalog simulators. The results show that aggregating information from several bands is highly beneficial in LSST sparsely-sampled time series, obtaining an absolute increase in period recovery rate up to 50%. We also show that the QMI is more robust to noise and light curve length (sample size) than the multiband generalizations of the Lomb Scargle and Analysis of Variance periodograms, recovering the true period in 10-30% more cases than its competitors. A python package containing efficient Cython implementations of the QMI and other methods is provided.

Interferometry in the era of time-domain astronomy

The physical nature of time variable objects is often inferred from photometric light-curves and spectroscopic variations. Long-baseline optical interferometry has the power to resolve the spatial structure of time variable sources directly in order to measure their physical properties and test the physics of the underlying models. Recent interferometric studies of variable objects include measuring the angular expansion and spatial structure during the early stages of novae outbursts, studying the transits and tidal distortions of the components in eclipsing and interacting binaries, measuring the radial pulsations in Cepheid variables, monitoring changes in the circumstellar discs around rapidly rotating massive stars, and imaging starspots. Future applications include measuring the image size and centroid displacements in gravitational microlensing events, and imaging the transits of exoplanets. Ongoing and upcoming photometric surveys will dramatically increase the number of time-variable objects detected each year, providing many potential targets to observe interferometrically. For short-lived transient events, it is critical for interferometric arrays to have the flexibility to respond rapidly to targets of opportunity and optimize the selection of baselines and beam combiners to provide the necessary resolution and sensitivity to resolve the source as its brightness and size change. We discuss the science opportunities made possible by resolving variable sources using long baseline optical interferometry.

Is the Hubble Constant Scale-Dependent?

An exact determination of the Hubble constant remains one of key problems in cosmology for almost a century. However, its modern values derived by various methods still disagree from each other by almost 10%; the greater values being obtained by measurements at the relatively small distances (e.g., by Cepheid stars as the standard candles), while the smaller values being characteristic of the methods associated with huge spatial scales (e.g., the analysis of the cosmic microwave background fluctuations). A reasonable way to resolve this puzzle is to assume that the Hubble constant is inherently scale-dependent. This idea seems to be particularly attractive in light of the latest observational results on the early-type galaxies, where the dark-matter haloes are almost absent. Therefore, an average contribution of the irregularly-distributed dark matter to the rate of the cosmological expansion should be substantially different at various spatial scales. As follows from the rough estimates, the corresponding variation of the Hubble constant can be 10% and even more, which well explains the spread in its values obtained by the various methods.

Magellanic Clouds Cepheids: Thorium Abundances

The analysis of the high-resolution spectra of 31 Magellanic Clouds Cepheid variables enabled the identification of thorium lines. The abundances of thorium were found with spectrum synthesis method. The calculated thorium abundances exhibit correlations with the abundances of other chemical elements and atmospheric parameters of the program stars. These correlations are similar for both Clouds. The correlations of iron abundances of thorium, europium, neodymium, and yttrium relative to the pulsational periods are different in the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), namely the correlations are negative for LMC and positive or close to zero for SMC. One of the possible explanations can be the higher activity of nucleosynthesis in SMC with respect to LMC in the recent several hundred million years.

New Parallaxes of Galactic Cepheids from Spatially Scanning the Hubble Space Telescope: Implications for the Hubble Constant

Abstract We present new measurements of the parallax of seven long-period (≥10 days) Milky Way (MW) Cepheid variables (SS CMa, XY Car, VY Car, VX Per, WZ Sgr, X Pup, and S Vul) using one-dimensional astrometric measurements from spatial scanning of Wide-Field Camera 3 on the Hubble Space Telescope (HST). The observations were obtained at ∼6 month intervals over 4 years. The distances are 1.7–3.6 kpc, with a mean precision of 45 μas (signal-to-noise ratio (S/N) ≈ 10) and a best precision of 29 μas (S/N = 14). The accuracy of the parallaxes is demonstrated through independent analyses of >100 reference stars. This raises to 10 the number of long-period Cepheids with significant parallax measurements, 8 obtained from this program. We also present high-precision mean F555W, F814W, and F160W magnitudes of these Cepheids, allowing a direct, zeropoint-independent comparison to >1800 extragalactic Cepheids in the hosts of 19 SNe Ia. This sample addresses two outstanding systematic uncertainties affecting prior comparisons of MW and extragalactic Cepheids used to calibrate the Hubble constant (H 0): their dissimilarity of periods and photometric systems. Comparing the new parallaxes to their predicted values derived from reversing the distance ladder gives a ratio (or independent scale for H 0) of 1.037 ± 0.036, consistent with no change and inconsistent at the 3.5σ level with a ratio of 0.91 needed to match the value predicted by Planck cosmic microwave background data in concert with ΛCDM. Using these data instead to augment the Riess et al. measurement of H 0 improves the precision to 2.3%, yielding 73.48 ± 1.66 km s−1 Mpc−1, and the tension with Planck + ΛCDM increases to 3.7σ. The future combination of Gaia parallaxes and HST spatial scanning photometry of 50 MW Cepheids can support a <1% calibration of H 0.

Spectroscopic properties of a two-dimensional time-dependent Cepheid model

Context. Standard spectroscopic analyses of variable stars are based on hydrostatic 1D model atmospheres. This quasi-static approach has not been theoretically validated. Aim. We aim at investigating the validity of the quasi-static approximation for Cepheid variables. We focus on the spectroscopic determination of the effective temperature Teff, surface gravity log g, microturbulent velocity ξt, and a generic metal abundance log A, here taken as iron.Methods. We calculated a grid of 1D hydrostatic plane-parallel models covering the ranges in effective temperature and gravity that are encountered during the evolution of a 2D time-dependent envelope model of a Cepheid computed with the radiation-hydrodynamics code CO5BOLD. We performed 1D spectral syntheses for artificial iron lines in local thermodynamic equilibrium by varying the microturbulent velocity and abundance. We fit the resulting equivalent widths to corresponding values obtained from our dynamical model for 150 instances in time, covering six pulsational cycles. In addition, we considered 99 instances during the initial non-pulsating stage of the temporal evolution of the 2D model. In the most general case, we treated Teff, log g, ξt, and log A as free parameters, and in two more limited cases, we fixed Teff and log g by independent constraints. We argue analytically that our approach of fitting equivalent widths is closely related to current standard procedures focusing on line-by-line abundances.Results. For the four-parametric case, the stellar parameters are typically underestimated and exhibit a bias in the iron abundance of ≈−0.2 dex. To avoid biases of this type, it is favorable to restrict the spectroscopic analysis to photometric phases ϕph ≈ 0.3…0.65 using additional information to fix the effective temperature and surface gravity.Conclusions. Hydrostatic 1D model atmospheres can provide unbiased estimates of stellar parameters and abundances of Cepheid variables for particular phases of their pulsations. We identified convective inhomogeneities as the main driver behind potential biases. To obtain a complete view on the effects when determining stellar parameters with 1D models, multidimensional Cepheid atmosphere models are necessary for variables of longer period than investigated here.

Magellanic Clouds Cepheids: Thorium Abundances

The analysis of the high-resolution spectra of 31 Magellanic Clouds Cepheid variables enabled the identification of thorium lines. The abundances of thorium were found with spectrum synthesis method. The calculated thorium abundances exhibit correlations with the abundances of other chemical elements and atmospheric parameters of the program stars. These correlations are similar for both Clouds. The correlations of iron abundances of thorium, europium, neodymium, and yttrium relative to the pulsational periods are different in the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), namely the correlations are negative for LMC and positive or close to zero for SMC. One of the possible explanations can be the higher activity of nucleosynthesis in SMC with respect to LMC in the recent several hundred million years.

First cosmological constraints combining Planck with the recent gravitational-wave standard siren measurement of the Hubble constant

The recent observations of gravitational-wave and electromagnetic emission produced by the merger of the binary neutron-star system GW170817 have opened the possibility of using standard sirens to constrain the value of the Hubble constant. While the reported bound of $H_0=70_{-8}^{+12}$ at $68 \%$ C.L. is significantly weaker than those recently derived by observations of Cepheid variables, it does not require any form of cosmic distance ladder and can be considered as complementary and, in principle, more conservative. Here we combine, for the first time, the new measurement with the Planck Cosmic Microwave Background observations in a $12$ parameters extended $\Lambda$CDM scenario, where the Hubble constant is weakly constrained from CMB data alone and bound to a low value $H_0=55^{+7}_{-20}$ km/s/Mpc at $68 \%$ C.L. We point out that the non-Gaussian shape of the GW170817 bound makes lower values of the Hubble constant in worst agreement with observations than what expected from a Gaussian form. The inclusion of the new GW170817 Hubble constant measurement therefore significantly reduces the allowed parameter space, improving the cosmological bounds on several parameters as the neutrino mass, curvature and the dark energy equation of state.

Cracking the Conundrum of F-supergiant Coronae

Abstract Chandra X-ray and HST far-ultraviolet (FUV) observations of three early-F supergiants have shed new light on a previous puzzle involving a prominent member of the class: α Persei (HD 20902: F5 Ib). The warm supergiant is a moderately strong, hard coronal ( K) X-ray source, but has 10 times weaker “subcoronal” Si iv 1393 Å ( K) emissions than early-G supergiants of similar high-energy properties. The α Per X-ray excess was speculatively ascribed to a close-in hyperactive G-dwarf companion, which could have escaped previous notice, lost in the glare of the bright star. However, a subsequent dedicated multi-wavelength imaging campaign failed to find any evidence for a resolved secondary. The origin of the α Per high-energy dichotomy then devolved to (1) an unresolved companion or (2) intrinsic coronal behavior. Exploring the second possibility, the present program has found that early-F supergiants do appear to belong to a distinct coronal class, characterized by elevated X-ray/FUV ratios, although sharing some similarities with Cepheid variables in their transitory X-ray “high states.” Remarkably, the early-F supergiants now are seen to align with the low-activity end of the X-ray/FUV sequence defined by late-type dwarfs, suggesting that the disjoint behavior relative to the G supergiants might be attributed to thinner outer atmospheres on the F types, as in dwarfs, but in this case perhaps caused by a weakened “ionization valve” effect due to overly warm photospheres.

Gaia DR1 Evidence of Disrupting the Perseus Arm

Abstract We have discovered a clear sign of the disruption phase of the Perseus arm in the Milky Way using Cepheid variables, taking advantage of the accurately measured distances of Cepheids and the proper motions from Gaia Data Release 1. Both the Galactocentric radial and rotation velocities of 77 Cepheids within 1.5 kpc of the Perseus arm are correlated with their distances from the locus of the Perseus arm, as the trailing side is rotating faster and moving inward compared to the leading side. We also found a negative vertex deviation for the Cepheids on the trailing side, −27.°6 ± 2.°4, in contrast to the positive vertex deviation in the solar neighborhood. This is, to our knowledge, the first direct evidence that the vertex deviation around the Perseus arm is affected by the spiral arm. We compared these observational trends with our N-body/hydrodynamics simulations based on a static density-wave spiral scenario and with those based on a transient dynamic spiral scenario. Although our comparisons are limited to qualitative trends, they strongly favor the conclusion that the Perseus arm is in the disruption phase of a transient arm.