Starspot mapping with adaptive parallel tempering I: Implementation of computational code

Abstract

Starspots are thought to be regions of locally strong magnetic fields, similar to sunspots, and they can generate photometric brightness modulations. To deduce stellar and spot properties, such as spot emergence and decay rates, we implement computational code for starspot modeling. It is implemented with an adaptive parallel tempering algorithm and an importance sampling algorithm for parameter estimation and model selection in the Bayesian framework. For evaluating the performance of the code, we apply it to synthetic light curves produced with 3 spots. The light curves are specified in the spot parameters, such as the radii, intensities, latitudes, longitudes, and emergence/decay durations. The spots are circular with specified radii and intensities relative to the photosphere, and the stellar differential rotation coefficient is also included in the light curves. As a result, stellar and spot parameters are uniquely deduced. The number of spots is correctly determined: the 3-spot model is preferable because the model evidence is much greater than that of 2-spot models by orders of magnitude and more than that of 4-spot model by a more modest factor, whereas the light curves are produced to have 2 or 1 local minimum during one equatorial rotation period by adjusting the values of longitude. The spot emergence and decay rates can be estimated with error less than an order of magnitude, considering the difference of the number of spots.

arXiv comment27 pages, 14 figures, 2 tables, accepted for publication in ApJ

References in corpus (41)

• emcee: The MCMC Hammer
• Exchange Monte Carlo Method and Application to Spin Glass Simulations
• Kepler Mission Design, Realized Photometric Performance, and Early Science
• Improving the Efficiency of Markov Chain Monte Carlo for Analyzing the Orbits of Extrasolar Planets
• Superflares on Solar Type Stars Observed with Kepler I. Statistical Properties of Superflares
• Quantifying the Uncertainty in the Orbits of Extrasolar Planets
• Rotation and differential rotation of active Kepler stars
• Data analysis recipes: Using Markov Chain Monte Carlo
• Dynamic temperature selection for parallel-tempering in Markov chain Monte Carlo simulations
• Flare-productive active regions
• Stellar Limb-Darkening Coefficients for CoRot and Kepler
• Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy
• Can Superflares Occur on Our Sun?
• PASTIS: Bayesian extrasolar planet validation. I. General framework, models, and performance
• Superflares on Solar-Type Stars Observed with Kepler II. Photometric Variability of Superflare-Generating Stars : A Signature of Stellar Rotation and Starspots
• A Kepler Study of Starspot Lifetimes with Respect to Light Curve Amplitude and Spectral Type
• Do Kepler superflare stars really include slowly-rotating Sun-like stars ? - Results using APO 3.5m telescope spectroscopic observations and Gaia-DR2 data -
• Photometric Variability in Kepler Target Stars: The Sun Among Stars -- A First Look
• The Starspots of HAT-P-11: Evidence for a Solar-like Dynamo
• RUN DMC: An efficient, parallel code for analyzing Radial Velocity Observations using N-body Integrations and Differential Evolution Markov chain Monte Carlo
• Short-lived spots in solar-like stars as observed by CoRoT
• Magnetic activity and differential rotation in the very young star KIC 8429280
• Starspot activity and superflares on solar-type stars
• Superflares on solar-type stars from the first year observation of TESS
• Long Term Photometric Variability in Kepler Full Frame Images: Magnetic Cycles of Sun-Like Stars
• The Information Content in Analytic Spot Models of Broadband Precision Lightcurves
• The Differential Rotation of Kappa1 Ceti as Observed by MOST
• Lifetimes and Emergence/Decay Rates of Star Spots on Solar-type Stars Estimated by Kepler Data in Comparison with Those of Sunspots
• High Dispersion Spectroscopy of Solar-type Superflare Stars. II. Stellar Rotation, Starspots, and Chromospheric Activities
• Magnetic activity and differential rotation in the young Sun-like stars KIC 7985370 and KIC 7765135
• An Analytic Model for Rotational Modulations in the Photometry of Spotted Stars
• Measuring stellar differential rotation with high-precision space-borne photometry
• On Sunspot and Starspot Lifetimes
• Double-Dipping: A New Relation between Stellar Rotation and Starspot Activity
• Two Sun-like Superflare Stars Rotating as Slow as the Sun
• High Dispersion Spectroscopy of Solar-type Superflare Stars. I. Temperature, Surface Gravity, Metallicity, and $v\sin i$
• Absolute densities in exoplanetary systems: photodynamical modelling of Kepler-138
• Temporal Evolution of Spatially-Resolved Individual Star Spots on a Planet-Hosting Solar-type Star: Kepler 17
• Calibration of Differential Light Curves for Physical Analysis of Starspots
• High Dispersion Spectroscopy of the Superflare Star KIC6934317
• Variations of the physical parameters of the blazar Mrk 421 based on the analysis of the spectral energy distributions

Cited by in corpus (7)

• Statistical Properties of Superflares on Solar-type Stars: Results Using All of the Kepler Primary Mission Data
• Superflares, chromospheric activities and photometric variabilities of solar-type stars from the second-year observation of TESS and spectra of LAMOST
• Starspot mapping with adaptive parallel tempering. II. Application to TESS data for M-dwarf flare stars, AU Microscopii, YZ Canis Minoris, and EV Lacertae
• Characterization of starspots on a young M-dwarf K2-25: multi-band observations of stellar photometric variability and planetary transits
• Doppler imaging combined with high-cadence photometry. I. Revisiting the surface of a pre-main-sequence flare star
• Temporal Variation of Flare Occurrence Rates via the Spot Evolution on the Sun and Solar-type Stars
• PAStar: a model for stellar surface from the Sun to active stars