Source author record

Mikko Tuomi

Mikko Tuomi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2025arXiv

Homogeneous search for spot transits in Kepler and TESS photometry of K $-$ M-type main-sequence stars

Late-type stars are known to host numerous exoplanets, and their photometric variability, primarily caused by rotational modulation, provides a unique opportunity to study starspots. As exoplanets transit in front of their host stars, they may occult darker, spotted regions on the stellar surfaces. The monitoring of starspots from planetary transits, known as transit mapping, offers a possibility to detect small dark regions on magnetically active, late-type stars. These spots may be so small that they would be undetectable to other methods used to reconstruct stellar magnetic activity. We describe a Bayesian analysis framework on the transit light curves of planets orbiting K- and M-type main-sequence stars in search for spot occultation event candidates. We present a systematic analysis of high-precision, high-cadence light curves from Kepler and TESS to detect and characterise starspots during exoplanetary transits. According to our tests, the set of criteria applied in the analysis is robust and not prone to false positives. Our sample comprises K- and M-dwarfs hosting transiting exoplanets observed by the Kepler or TESS space telescopes at a high cadence, totalling 99 planets meeting our selection criteria. After analysing 3273 transit light curves from 99 planets, we find 105 candidates for starspot occultation events by six planets. We report new spot occultation candidates for the K-dwarfs HD 189733 and TOI-1268. The identified dark regions have a lower limit for radii between 1.6 degrees and 29.5 degrees and contrasts up to 0.69. We estimate a spot detection frequency of 3.7% and 4.2% for K- and M-dwarfs by TESS, and 37.5% for K-dwarfs by Kepler.

preprint2020arXiv

The Magellan/PFS Exoplanet Search: A 55-day period dense Neptune transiting the bright ($V=8.6$) star HD 95338

We report the detection of a transiting, dense Neptune planet candidate orbiting the bright ($V=8.6$) K0.5V star HD 95338. Detection of the 55-day periodic signal comes from the analysis of precision radial velocities from the Planet Finder Spectrograph on the Magellan II Telescope. Follow-up observations with HARPS also confirm the presence of the periodic signal in the combined data. HD 95338 was also observed by the Transiting Exoplanet Survey Satellite ({\it TESS}) where we identify a clear single transit in the photometry. A Markov Chain Monte Carlo period search on the velocities allows strong constraints on the expected transit time, matching well the epoch calculated from \tess{} data, confirming both signals describe the same companion. A joint fit model yields an absolute mass of 42.44$^{+2.22}_{-2.08} M_{\oplus}$ and a radius of 3.89$^{+0.19}_{-0.20}$ $R_{\oplus}$ which translates to a density of 3.98$^{+0.62}_{-0.64}$ \gcm\, for the planet. Given the planet mass and radius, structure models suggest it is composed of a mixture of ammonia, water, and methane. HD 95338\,b is one of the most dense Neptune planets yet detected, indicating a heavy element enrichment of $\sim$90\% ($\sim38\, M_{\oplus}$). This system presents a unique opportunity for future follow-up observations that can further constrain structure models of cool gas giant planets.

preprint2016arXiv

A terrestrial planet candidate in a temperate orbit around Proxima Centauri

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($α$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.

preprint2015arXiv

Benchmarking the power of amateur observatories for TTV exoplanets detection

We perform an analysis of ~80000 photometric measurements for the following 10 stars hosting transiting planets: WASP-2, -4, -5, -52, Kelt-1, CoRoT-2, XO-2, TrES-1, HD 189733, GJ 436. Our analysis includes mainly transit lightcurves from the Exoplanet Transit Database, public photometry from the literature, and some proprietary photometry privately supplied by other authors. Half of these lightcurves were obtained by amateurs. From this photometry we derive 306 transit timing measurements, as well as improved planetary transit parameters. Additionally, for 6 of these 10 stars we present a set of radial velocity measurements obtained from the spectra stored in the HARPS, HARPS-N, and SOPHIE archives using the HARPS-TERRA pipeline. Our analysis of these TTV and RV data did not reveal significant hints of additional orbiting bodies in almost all of the cases. In the WASP-4 case, we found hints of marginally significant TTV signals having amplitude 10-20 sec, although their parameters are model-dependent and uncertain, while radial velocities did not reveal statistically significant Doppler signals.

preprint2015arXiv

Comment on "Stellar activity masquerading as planets in the habitable zone of the M dwarf Gliese 581"

Robertson et al.(Reports, July 25 2014, p440-444)(1) claimed that activity-induced variability is responsible for the Doppler signal of the proposed planet candidate GJ 581d. We point out that their analysis using periodograms of residual data is incorrect, further promoting inadequate tools. Since the claim challenges the viability of the method to detect exo-Earths, we urge for more appropriate analyses (see appendix).

preprint2015arXiv

No evidence for activity correlations in the radial velocities of Kapteyn's star

Stellar activity may induce Doppler variability at the level of a few m/s which can then be confused by the Doppler signal of an exoplanet orbiting the star. To first order, linear correlations between radial velocity measurements and activity indices have been proposed to account for any such correlation. The likely presence of two super-Earths orbiting Kapteyn's star was reported in Anglada et al. (2014, MNRAS 443L, 89A), but this claim was recently challenged by Robertson et al. (2015, ApJ 805L, 22R) arguing evidence of a rotation period (143 days) at three times the orbital period of one of the proposed planets (Kapteyn's b, P=48.6 days), and the existence of strong linear correlations between its Doppler signal and activity data. By re-analyzing the data using global optimization methods and model comparison, we show that such claim is incorrect given that; 1) the choice of a rotation period at 143 days is unjustified, and 2) the presence of linear correlations is not supported by the data. We conclude that the radial velocity signals of Kapteyn's star remain more simply explained by the presence of two super-Earth candidates orbiting it. We also advocate for the use of global optimization procedures and objective arguments, instead of claims lacking of a minimal statistical support.

preprint2014arXiv

A new cold sub-Saturnian candidate planet orbiting GJ 221

We re-analyse the recently published HARPS and PFS velocities of the nearby K dwarf GJ 221 that have been reported to contain the signatures of two planets orbiting the star. Our goal is to see whether the earlier studies discussing the system fell victims of false negative detections. We perform the analyses by using an independent statistical method based on posterior samplings and model comparisons in the Bayesian framework that is known to be more sensitive to weak signals of low-mass planets. According to our analyses, we find strong evidence in favour of a third candidate planet in the system corresponding to a cold sub-Saturnian planet with an orbital period of 500 days and a minimum mass of 29 $M_{\oplus}$. Application of sub-optimal signal detection methods can leave low-amplitude signals undetected in radial velocity time-series. Our results suggest that the estimated statistical properties of low-mass planets can thus be biased because several signals corresponding to low-mass candidate planets may have gone unnoticed. This also suggests that the occurrence rates of such planets based on radial velocity surveys might be underestimated.

preprint2014arXiv

Bayesian search for low-mass planets around nearby M dwarfs. Estimates for occurrence rate based on global detectability statistics

Due to their higher planet-star mass-ratios, M dwarfs are the easiest targets for detection of low-mass planets orbiting nearby stars using Doppler spectroscopy. Furthermore, because of their low masses and luminosities, Doppler measurements enable the detection of low-mass planets in their habitable zones that correspond to closer orbits than for Solar-type stars. We re-analyse literature UVES radial velocities of 41 nearby M dwarfs in a combination with new velocities obtained from publicly available spectra from the HARPS-ESO spectrograph of these stars in an attempt to constrain any low-amplitude Keplerian signals. We apply Bayesian signal detection criteria, together with posterior sampling techniques, in combination with noise models that take into account correlations in the data and obtain estimates for the number of planet candidates in the sample. More generally, we use the estimated detection probability function to calculate the occurrence rate of low-mass planets around nearby M dwarfs. We report eight new planet candidates in the sample (orbiting GJ 27.1, GJ 160.2, GJ 180, GJ 229, GJ 422, and GJ 682), including two new multiplanet systems, and confirm two previously known candidates in the GJ 433 system based on detections of Keplerian signals in the combined UVES and HARPS radial velocity data that cannot be explained by periodic and/or quasiperiodic phenomena related to stellar activities. Finally, we use the estimated detection probability function to calculate the occurrence rate of low-mass planets around nearby M dwarfs. According to our results, M dwarfs are hosts to an abundance of low-mass planets and the occurrence rate of planets less massive than 10 M$_{\oplus}$ is of the order of one planet per star, possibly even greater. ...

preprint2014arXiv

Filtering out activity-related variations from radial velocities in a search for low-mass planets

Variations related to stellar activity and correlated noise can prevent the detections of low-amplitude signals in radial velocity data if not accounted for. This can be seen as the greatest obstacle in detecting Earth-like planets orbiting nearby stars with Doppler spectroscopy regardless of developments in instrumentation and rapidly accumulating amounts of data. We use a statistical model that is not sensitive to aperiodic and/or quasiperiodic variability of stellar origin. We demonstrate the performance of our model by re-analysing the radial velocities of the moderately active star CoRoT-7 ($\log R_{\rm HK} = -4.61$) with a transiting planet whose Doppler signal has proven rather difficult to detect. We find that the signal of the transiting planet can be robustly detected together with signals of two other planet candidates. Our results suggest that rotation periods of moderately active stars can be filtered out of the radial velocity noise, which enables the detections of low-mass planets orbiting such stars.

preprint2014arXiv

GJ 832c: A super-earth in the habitable zone

We report the detection of GJ 832c, a super-Earth orbiting near the inner edge of the habitable zone of GJ 832, an M dwarf previously known to host a Jupiter analog in a nearly-circular 9.4-year orbit. The combination of precise radial-velocity measurements from three telescopes reveals the presence of a planet with a period of 35.68+/-0.03 days and minimum mass (m sin i) of 5.4+/-1.0 Earth masses. GJ 832c moves on a low-eccentricity orbit (e=0.18+/-0.13) towards the inner edge of the habitable zone. However, given the large mass of the planet, it seems likely that it would possess a massive atmosphere, which may well render the planet inhospitable. Indeed, it is perhaps more likely that GJ 832c is a "super-Venus," featuring significant greenhouse forcing. With an outer giant planet and an interior, potentially rocky planet, the GJ 832 planetary system can be thought of as a miniature version of our own Solar system.

preprint2014arXiv

The Anglo-Australian Planet Search. XXIII. Two New Jupiter Analogs

We report the discovery of two long-period giant planets from the Anglo-Australian Planet Search. HD 154857c is in a multiple-planet system, while HD 114613b appears to be solitary. HD 114613b has an orbital period P=10.5 years, and a minimum mass m sin i of 0.48 Jupiter masses; HD 154857c has P=9.5 years and m sin i=2.6 Jupiter masses. These new data confirm the planetary nature of the previously unconstrained long-period object in the HD 154857 system. We have performed detailed dynamical stability simulations which show that the HD 154857 two-planet system is stable on timescales of at least 100 million years. These results highlight the continued importance of "legacy" surveys with long observational baselines; these ongoing campaigns are critical for determining the population of Jupiter analogs, and hence of those planetary systems with architectures most like our own Solar system.

preprint2014arXiv

The curious case of HD41248. A pair of static signals buried behind red-noise

Gaining a better understanding of the effects of stellar induced radial velocity noise is critical for the future of exoplanet studies, since the discovery of the lowest-mass planets using this method will require us to go below the intrinsic stellar noise limit. An interesting test case in this respect is that of the southern solar analogue HD41248. The radial velocity time series of this star has been proposed to contain either a pair of signals with periods of around 18 and 25 days, that could be due to a pair of resonant super-Earths, or a single and varying 25 day signal that could arise due to a complex interplay between differential rotation and modulated activity. In this letter we build-up more evidence for the former scenario, showing that the signals are still clearly significant even after more than 10 years of observations and they likely do not change in period, amplitude, or phase as a function of time, the hallmarks of static Doppler signals. We show that over the last two observing seasons this star was more intrinsically active and the noise reddened, highlighting why better noise models are needed to find the lowest amplitude signals, in particular models that consider noise correlations. This analysis shows that there is still sufficient evidence for the existence of two super-Earths on the edge of, or locked into, a 7:5 mean motion resonance orbiting HD41248.

preprint2014arXiv

Two planets around Kapteyn's star : a cold and a temperate super-Earth orbiting the nearest halo red-dwarf

Exoplanets of a few Earth masses can be now detected around nearby low-mass stars using Doppler spectroscopy. In this paper, we investigate the radial velocity variations of Kapteyn's star, which is both a sub-dwarf M-star and the nearest halo object to the Sun. The observations comprise archival and new HARPS, HIRES and PFS Doppler measurements. Two Doppler signals are detected at periods of 48 and 120 days using likelihood periodograms and a Bayesian analysis of the data. Using the same techniques, the activity indicies and archival ASAS-3 photometry show evidence for low-level activity periodicities of the order of several hundred days. However, there are no significant correlations with the radial velocity variations on the same time-scales. The inclusion of planetary Keplerian signals in the model results in levels of correlated and excess white noise that are remarkably low compared to younger G, K and M dwarfs. We conclude that Kapteyn's star is most probably orbited by two super-Earth mass planets, one of which is orbiting in its circumstellar habitable zone, becoming the oldest potentially habitable planet known to date. The presence and long-term survival of a planetary system seems a remarkable feat given the peculiar origin and kinematic history of Kapteyn's star. The detection of super-Earth mass planets around halo stars provides important insights into planet-formation processes in the early days of the Milky Way.

preprint2013arXiv

A dynamically-packed planetary system around GJ667C with three super-Earths in its habitable zone

Since low-mass stars have low luminosities, orbits at which liquid water can exist on Earth-sized planets are relatively close-in, which produces Doppler signals that are detectable using state-of-the-art Doppler spectroscopy. GJ 667C is already known to be orbited by two super-Earth candidates. We investigate whether the data supports the presence of additional companions. We obtain new Doppler measurements from HARPS extracted spectra and combined them with those obtained from the PFS and HIRES spectrographs. We used Bayesian and periodogram-based methods to re-assess the number of candidates and evaluated the confidence of each detection. Among other tests, we validated the planet candidates by analyzing correlations of each Doppler signal activity indices and investigate quasi-periodicity. Doppler measurements of GJ 667C are described better by six Keplerian-like signals: the two known candidates (b and c); three additional few-Earth mass candidates with periods of 92, 62 and 39 days (d, e and f); a cold super-Earth in a 260-day orbit (g) and tantalizing evidence of a $\sim$ 1 M$_\oplus$ object in a close-in orbit of 17 days (h). We explore whether long-term stable orbits are compatible with the data by integrating 8$\times 10^4$ solutions derived from the Bayesian samplings. The system consisting of six planets is compatible with dynamically stable configurations. As for the solar system, the most stable solutions do not contain mean-motion resonances and are described well by analytic Laplace-Lagrange solutions. The presence of a seventh planet (h) is supported by the fact that it appears squarely centered on the only island of stability left in the six-planet solution. Habitability assessments accounting for the stellar flux, as well as tidal dissipation effects, indicate that three (maybe four) planets are potentially habitable...

preprint2013arXiv

Signals embedded in the radial velocity noise

Bayesian data analysis techniques, together with suitable statistical models, can be used to obtain much more information from noisy data than the traditional frequentist methods. For instance, when searching for periodic signals in noisy data, the Bayesian techniques can be used to define exact detection criteria for low-amplitude signals - the most interesting signals that might correspond to habitable planets. We present an overview of Bayesian techniques and present detailed analyses of the HARPS-TERRA velocities of HD 40307, a nearby star observed to host a candidate habitable planet, to demonstrate in practice the applicability of Bayes' rule to astronomical data.

preprint2013arXiv

Status of the Calan-Hertfordshire Extrasolar Planet Search

In these proceedings we give a status update of the Calan-Hertfordshire Extrasolar Planet Search, an international collaboration led from Chile that aims to discover more planets around super metal-rich and Sun-like stars, and then follow these up with precision photometry to hunt for new bright transit planets. We highlight some results from this program, including exoplanet and brown dwarf discoveries, and a possible correlation between metallicity and planetary minimum mass at the lowest planetary masses detectable. Finally we discuss the short-term and long-term future pathways this program can take.

preprint2013arXiv

Two Super-Earths Orbiting the Solar Analogue HD41248 on the edge of a 7:5 Mean Motion Resonance

The number of multi-planet systems known to be orbiting their host stars with orbital periods that place them in mean motion resonances is growing. For the most part, these systems are in first-order resonances and dynamical studies have focused their efforts towards understanding the origin and evolution of such dynamically resonant commensurabilities. We report here the discovery of two super-Earths that are close to a second-order dynamical resonance, orbiting the metal-poor ([Fe/H]=-0.43 dex) and inactive G2V star HD41248. We analysed 62 HARPS archival radial velocities for this star, that until now, had exhibited no evidence for planetary companions. Using our new Bayesian Doppler signal detection algorithm, we find two significant signals in the data, with periods of 18.357 days and 25.648 days, indicating they could be part of a 7:5 second-order mean motion resonance. Both semi-amplitudes are below 3m/s and the minimum masses of the pair are 12.3 and 8.6Mearth, respectively. Our simulations found that apsidal alignment stabilizes the system, and even though libration of the resonant angles was not seen, the system is affected by the presence of the resonance and could yet occupy the 7:5 commensurability, which would be the first planetary configuration in such a dynamical resonance. Given the multitude of low-mass multiplanet systems that will be discovered in the coming years, we expect more of these second-order resonant configurations will emerge from the data, highlighting the need for a better understanding of the dynamical interactions between forming planetesimals.

preprint2013arXiv

Up to four planets around the M dwarf GJ 163. Sensitivity of Bayesian planet detection criteria to prior choice

Exoplanet Doppler surveys are currently the most efficient means to detect low-mass companions to nearby stars. Among these stars, the light M dwarfs provide the highest sensitivity to detect low-mass exoplanet candidates. Evidence is accumulating that a substantial fraction of these low-mass planets are found in high-multiplicity planetary systems. GJ 163 is a nearby inactive M dwarf with abundant public observations obtained using the HARPS spectrograph. We obtain and analyse radial velocities from the HARPS public spectra of GJ 163 and investigate the presence of a planetary companions orbiting it. The number of planet candidates detected might depend on some prior assumptions. Since the impact of prior choice has not been investigated throughly previously, we study the effects of different prior densities on the detectability of planet candidates around GJ 163. We use Bayesian tools, i.e. posterior samplings and model comparisons, when analysing the GJ 163 velocities. We consider models accounting for the possible correlations of subsequent measurements. We also search for activity-related counterparts of the signals we observe and test the dynamical stability of the planetary systems corresponding to our solutions using direct numerical integrations of the orbits. We find that there are at least three planet candidates orbiting GJ 163. The existence of a fourth planet is supported by the data but the evidence in favor of the corresponding model is not yet conclusive. The second innermost planet candidate in the system with an orbital period of 25.6 days and a minimum mass of 8.7 Me is inside the liquid-water habitable zone of the star.

preprint2012arXiv

A Hot Uranus Orbiting the Super Metal-rich Star HD77338 and the Metallicity - Mass Connection

We announce the discovery of a low-mass planet orbiting the super metal-rich K0V star HD77338 as part of our on-going Calan-Hertfordshire Extrasolar Planet Search. The best fit planet solution has an orbital period of 5.7361\pm0.0015 days and with a radial velocity semi-amplitude of only 5.96\pm1.74 m/s, we find a minimum mass of 15.9+4.7-5.3 Me. The best fit eccentricity from this solution is 0.09+0.25-0.09, and we find agreement for this data set using a Bayesian analysis and a periodogram analysis. We measure a metallicity for the star of +0.35\pm0.06 dex, whereas another recent work (Trevisan et al. 2011) finds +0.47\pm0.05 dex. Thus HD77338b is one of the most metal-rich planet host stars known and the most metal-rich star hosting a sub-Neptune mass planet. We searched for a transit signature of HD77338b but none was detected. We also highlight an emerging trend where metallicity and mass seem to correlate at very low masses, a discovery that would be in agreement with the core accretion model of planet formation. The trend appears to show that for Neptune-mass planets and below, higher masses are preferred when the host star is more metal-rich. Also a lower boundary is apparent in the super metal-rich regime where there are no very low-mass planets yet discovered in comparison to the sub-solar metallicity regime. A Monte Carlo analysis shows that this, low-mass planet desert, is statistically significant with the current sample of 36 planets at around the 4.5σ level. In addition, results from Kepler strengthen the claim for this paucity of the lowest-mass planets in super metal-rich systems. Finally, this discovery adds to the growing population of low-mass planets around low-mass and metal-rich stars and shows that very low-mass planets can now be discovered with a relatively small number of data points using stable instrumentation.

preprint2012arXiv

A planetary system with gas giants and super-Earths around the nearby M dwarf GJ 676A. Optimizing data analysis techniques for the detection of multi-planetary systems

Several M dwarfs are targets of systematical monitoring in searches for Doppler signals caused by low-mass exoplanet companions. As a result, an emerging population of high-multiplicity planetary systems around low-mass stars are being detected as well. We optimize classic data analysis methods and develop new ones to enhance sensitivity towards low-amplitude planets in high-multiplicity systems. We apply these methods to the public HARPS observations of GJ 676A, a nearby M dwarf with one reported gas giant companion. We re-derived Doppler measurements using the template matching method (HARPS-TERRA). We used refined versions of periodograms to assess the presence of additional low-mass companions. We analyse the same dataset with Bayesian tools and compared the performance of both approaches. We confirm the reported massive gas giant candidate and a long period trend, whose curvature is now well detected. We also find very secure evidence of two new candidates in close-in orbits and masses in the super-Earth mass regime. Despite the increased sensitivity of the new periodogram tools, we find that Bayesian methods are more sensitive in the early detection of candidate signals. While hardware development is important, development of data analysis techniques can help to reveal new results from existing data sets with significantly fewer resources. This new system holds the record of minimum-mass range (from Msin i = 4.5 M_Earth to 5 M_Jup) and period range (from 3.6 days to more than 10 years). Although all the planet candidates are substantially more massive, it is the first exoplanetary system with a general architecture similar to our solar system. GJ 676A can be happily added to the family of high-multiplicity planetary systems around M dwarfs.

preprint2012arXiv

Counting the number of planets around GJ 581. False positive rate of Bayesian signal detection methods

The four-planet system around GJ 581 has received attention because it has been claimed that there are possibly two additional low-mass companions as well - one of them being a planet in the middle of the stellar habitable zone. We re-analyse the available HARPS and HIRES Doppler data in an attempt to determine the false positive rate of our Bayesian data analysis techniques and to count the number of Keplerian signals in the GJ 581 data. We apply the common Lomb-Scargle periodograms and posterior sampling techniques in the Bayesian framework to estimate the number of signals in the radial velocities. We also analyse the HARPS velocities sequentially after each full observing period to compare the sensitivities and false positive rates of the two signal detection techniques. By relaxing the assumption that the radial velocity noise is white, we also demonstrate the consequences that noise correlations have on the obtained results and the significances of the signals. According to our analyses, the number of Keplerian signals favoured by the publicly available HARPS and HIRES radial velocity data of GJ 581 is four. This result relies on the sensitivity of the Bayesian statistical analysis techniques but also depends on the assumed noise model. We also show that the radial velocity noise is actually not white and that this feature has to be accounted for when analysing radial velocities in a search for low-amplitude signals corresponding to low-mass planets. ...

preprint2012arXiv

Evidence for 9 planets in the HD 10180 system

We re-analyse the HARPS radial velocities of HD 10180 and calculate the probabilities of models with differing numbers of periodic signals in the data. We test the significance of the seven signals, corresponding to seven exoplanets orbiting the star, in the Bayesian framework and perform comparisons of models with up to nine periodicities. We use posterior samplings and Bayesian model probabilities in our analyses together with suitable prior probability densities and prior model probabilities to extract all the significant signals from the data and to receive reliable uncertainties for the orbital parameters of the six, possibly seven, known exoplanets in the system. According to our results, there is evidence for up to nine planets orbiting HD 10180, which would make this this star a record holder in having more planets in its orbits than there are in the Solar system. We revise the uncertainties of the previously reported six planets in the system, verify the existence of the seventh signal, and announce the detection of two additional statistically significant signals in the data. If of planetary origin, these two additional signals would correspond to planets with minimum masses of 5.1$^{+3.1}_{-3.2}$ and 1.9$^{+1.6}_{-1.8}$ M$_{\oplus}$ on orbits with 67.55$^{+0.68}_{-0.88}$ and 9.655$^{+0.022}_{-0.072}$ days periods (denoted using the 99% credibility intervals), respectively.

preprint2012arXiv

Habitable-zone super-Earth candidate in a six-planet system around the K2.5V star HD 40307

The K2.5 dwarf HD 40307 has been reported to host three super-Earths. The system lacks massive planets and is therefore a potential candidate for having additional low-mass planetary companions. We re-derive Doppler measurements from public HARPS spectra of HD 40307 to confirm the significance of the reported signals using independent data analysis methods. We also investigate these measurements for additional low-amplitude signals. We used Bayesian analysis of our radial velocities to estimate the probability densities of different model parameters. We also estimated the relative probabilities of models with differing numbers of Keplerian signals and verified their significance using periodogram analyses. We investigated the relation of the detected signals with the chromospheric emission of the star. As previously reported for other objects, we found that radial velocity signals correlated with the S-index are strongly wavelength dependent. We identify two additional clear signals with periods of 34 and 51 days, both corresponding to planet candidates with minimum masses a few times that of the Earth. An additional sixth candidate is initially found at a period of 320 days. However, this signal correlates strongly with the chromospheric emission from the star and is also strongly wavelength dependent. When analysing the red half of the spectra only, the five putative planetary signals are recovered together with a very significant periodicity at about 200 days. This signal has a similar amplitude as the other new signals reported in the current work and corresponds to a planet candidate with M sin i = 7 Me (HD 40307 g). ...

preprint2012arXiv

Probabilities of exoplanet signals from posterior samplings

Estimating the marginal likelihoods is an essential feature of model selection in the Bayesian context. It is especially crucial to have good estimates when assessing the number of planets orbiting stars when the models explain the noisy data with different numbers of Keplerian signals. We introduce a simple method for approximating the marginal likelihoods in practice when a statistically representative sample from the parameter posterior density is available. We use our truncated posterior mixture estimate to receive accurate model probabilities for models with differing number of Keplerian signals in radial velocity data. We test this estimate in simple scenarios to assess its accuracy and rate of convergence in practice when the corresponding estimates calculated using deviance information criterion can be applied to receive trustworthy results for reliable comparison. As a test case, we determine the posterior probability of a planet orbiting HD 3651 given Lick and Keck radial velocity data. The posterior mixture estimate appears to be a simple and an accurate way of calculating marginal integrals from posterior samples. We show, that it can be used to estimate the marginal integrals reliably in practice, given a suitable selection of parameter λ, that controls its accuracy and convergence rate. It is also more accurate than the one block Metropolis-Hastings estimate and can be used in any application because it is not based on assumptions on the nature of the posterior density nor the amount of data or parameters in the statistical model.

preprint2012arXiv

Signals embedded in the radial velocity noise. Periodic variations in the tau Ceti velocities

The abilities of radial velocity exoplanet surveys to detect the lowest-mass extra-solar planets are currently limited by a combination of instrument precision, lack of data, and "jitter". Jitter is a general term for any unknown features in the noise, and reflects a lack of detailed knowledge of stellar physics (asteroseismology, starspots, magnetic cycles, granulation, and other stellar surface phenomena), as well as the possible underestimation of instrument noise. We study an extensive set of radial velocities for the star HD 10700 ($τ$ Ceti) to determine the properties of the jitter arising from stellar surface inhomogeneities, activity, and telescope-instrument systems, and perform a comprehensive search for planetary signals in the radial velocities. We perform Bayesian comparisons of statistical models describing the radial velocity data to quantify the number of significant signals and the magnitude and properties of the excess noise in the data. We reach our goal by adding artificial signals to the "flat" radial velocity data of HD 10700 and by seeing which one of our statistical noise models receives the greatest posterior probabilities while still being able to extract the artificial signals correctly from the data. We utilise various noise components to assess properties of the noise in the data and analyse the HARPS, AAPS, and HIRES data for HD 10700 to quantify these properties and search for previously unknown low-amplitude Keplerian signals. ...

preprint2011arXiv

Application of Bayesian model inadequacy criterion for multiple data sets to radial velocity models of exoplanet systems

We present a simple mathematical criterion for determining whether a given statistical model does not describe several independent sets of measurements, or data modes, adequately. We derive this criterion for two data sets and generalise it to several sets by using the Bayesian updating of the posterior probability density. To demonstrate the usage of the criterion, we apply it to observations of exoplanet host stars by re-analysing the radial velocities of HD 217107, Gliese 581, and \u{psion} Andromedae and show that the currently used models are not necessarily adequate in describing the properties of these measurements. We show that while the two data sets of Gliese 581 can be modelled reasonably well, the noise model of HD 217107 needs to be revised. We also reveal some biases in the radial velocities of \u{psion} Andromedae and report updated orbital parameters for the recently proposed 4-planet model. Because of the generality of our criterion, no assumptions are needed on the nature of the measurements, models, or model parameters. The method we propose can be applied to any astronomical problems, as well as outside the field of astronomy, because it is a simple consequence of the Bayes' rule of conditional probabilities.

preprint2011arXiv

Bayesian re-analysis of the radial velocities of Gliese 581. Evidence in favour of only four planetary companions

The Gliese 581 planetary system has received attention because it has been proposed to host a low-mass planet in its habitable zone. We re-analyse the radial velocity measurements reported to contain six planetary signals to see whether these conclusions remain valid when the analyses are made using Bayesian tools instead of the common periodogram analyses. We analyse the combined radial velocity data set obtained using the HARPS and HIRES spectrographs using posterior sampling techniques and computation of the posterior probabilities of models with differing numbers of Keplerian signals. We do not fix the orbital eccentricities and stellar jitter to certain values but treat these as free parameters of our statistical models. Hence, we can take the uncertainties of these parameters into account when assessing the number of planetary signals present in the data, the point estimates of all of the model parameters, and the uncertainties of these parameters. We conclude that based on the Bayesian model probabilities and the nature of the posterior densities of the different models, there is evidence in favour of four planets orbiting GJ 581. The HARPS and HIRES data do not imply the conclusion that there are two additional companions orbiting GJ 581. We also revise the orbital parameters of the four companions in the system. Especially, according to our results, the eccentricities of all the companions in the system are consistent with zero.

preprint2011arXiv

Soil carbon model Yasso07 graphical user interface

In this article, we present a graphical user interface software for the litter decomposition and soil carbon model Yasso07 and an overview of the principles and formulae it is based on. The software can be used to test the model and use it in simple applications. Yasso07 is applicable to upland soils of different ecosystems worldwide, because it has been developed using data covering the global climate conditions and representing various ecosystem types. As input information, Yasso07 requires data on litter input to soil, climate conditions, and land-use change if any. The model predictions are given as probability densities representing the uncertainties in the parameter values of the model and those in the input data - the user interface calculates these densities using a built-in Monte Carlo simulation.