---
title: 'TRICERATOPS: Bayesian Exoplanet Vetting'
url: https://www.emergentmind.com/topics/triceratops
type: topic
---

# TRICERATOPS: Bayesian Exoplanet Vetting

TRICERATOPS is a Bayesian vetting and statistical validation framework for transiting exoplanet candidates, developed for TESS Objects of Interest and designed to estimate whether a transit-like signal is more likely to be produced by a bona fide planet or by an astrophysical false positive [2002.00691]. Its defining feature is explicit multi-star modeling: it evaluates scenarios on the target star, on unresolved bound or background/foreground companions, and on resolved nearby stars identified in Gaia and the TIC, and it reports both a conventional false positive probability (FPP) and a nearby false positive probability (NFPP) [2002.00691]. This architecture is motivated by the TESS observing regime, in which pixels are \(21\) arcsec and contamination by nearby sources is common [2002.00691].

## 1. Origin, scope, and design context

TRICERATOPS was introduced as a fast Bayesian tool tailored to TESS, whose large pixels and crowded apertures make explicit contamination modeling a central part of transit vetting [2002.00691]. In its original presentation, it was described as the only TESS vetting and validation tool that models transits from nearby contaminant stars in addition to the target star, rather than assuming that nearby-star pathways have already been excluded by separate pre-vetting [2002.00691]. The framework works directly with TIC and Gaia stellar properties, aperture definitions, physical transit models, and optional follow-up constraints such as high-resolution imaging contrast curves [2002.00691].

The original benchmark applied the method to 68 previously classified TOIs and then to 384 unclassified TOIs. In that application, TRICERATOPS statistically validated 12 objects, classified 125 as likely planets, and classified 52 as likely nearby false positives [2002.00691]. Later studies extended its use from bulk triage to system-specific statistical validation, including surveys of hot Neptunes, K-star hosts, and super-Earth candidates, as well as detailed analyses of individual TOIs with extensive ground-based follow-up [2211.08490] [2301.09865] [2311.00688].

This broad usage reflects a specific niche. TRICERATOPS is not principally an instrumental false-alarm detector; rather, it is a probabilistic astrophysical scenario comparator. As a result, it is routinely embedded in larger vetting pipelines that also inspect centroid motion, odd-even depth differences, secondary eclipses, archival imaging, and follow-up photometry [2002.00691] [2211.08490].

## 2. Bayesian formulation and hypothesis space

The framework computes posterior probabilities for competing scenarios \(S_j\) using Bayes’ theorem. In the original formulation, the posterior is written as \(p(S_j\mid D)\propto p(S_j)\,p(D\mid S_j)\), with normalized scenario probabilities
\[
\mathcal{P}_j=\frac{p(S_j\mid D)}{\sum_j p(S_j\mid D)}.
\]
The marginal likelihood for each scenario is estimated by Monte Carlo sampling, typically with \(N=10^6\), and the transit-data likelihood is evaluated from modeled fluxes relative to the observed phase-folded light curve [2002.00691]. Later papers retain the same Bayesian structure, often expressing it as
\[
P(H_i\mid D)=\frac{P(D\mid H_i)\,P(H_i)}{\sum_j P(D\mid H_j)\,P(H_j)}.
\]
[2602.14840] [2508.02782]

In the original paper, TRICERATOPS enumerates 18 astrophysical scenarios [2002.00691]. These include target-star cases without unresolved companions, cases with unresolved bound companions, cases with unresolved background/foreground stars, and scenarios on resolved nearby stars. The scenario vocabulary used throughout the later literature centers on a stable subset of labels: TP for a transiting planet on the target, PTP for a planet on the primary when an unresolved bound companion is present, DTP for a diluted transiting planet on the target in the presence of a blended background/foreground star, NTP for a transiting planet on a nearby resolved star, and EB-, HEB-, and BEB-like families for eclipsing-binary false positives [2002.00691] [2209.13651] [2508.02782].

The summary statistics are derived from these scenario posteriors. In the original paper,
\[
\mathrm{FPP}=1-(\mathcal{P}_{\mathrm{TP}}+\mathcal{P}_{\mathrm{PTP}}+\mathcal{P}_{\mathrm{DTP}})
\]
and
\[
\mathrm{NFPP}=\sum(\mathcal{P}_{\mathrm{NTP}}+\mathcal{P}_{\mathrm{NEB}}+\mathcal{P}_{\mathrm{NEBx2P}}).
\]
[2002.00691] Later applications often use the equivalent compact interpretation that FPP is the posterior mass assigned to all false-positive scenarios, while NFPP isolates the nearby-source subset [2209.13651] [2602.14840].

A recurring feature of published analyses is that low TP probability does not imply high FPP if posterior mass is concentrated in other planet-containing scenarios. This is stated explicitly for TOI-7701.01, where TP accounts for \(78.1\%\) of the posterior, but total planet-containing configurations account for \(99.7\%\), leaving \(\mathrm{FPP}=0.00191\) and \(\mathrm{NFPP}<10^{-6}\) [2606.15430].

## 3. Inputs, priors, and forward modeling

TRICERATOPS requires the transit light curve, ephemeris, an aperture definition, host-star properties, and a catalog of nearby stars. In the original implementation, users provide the TIC ID, phase-folded TESS light curve and uncertainty, the observed transit depth \(\delta_{\mathrm{obs}}\), the ephemeris, and the extraction aperture for each sector; the code then queries nearby stars within a 10-pixel radius and computes per-star flux contributions \(X_s\) using a Gaussian approximation to the PRF with \(\sigma\approx0.75\) pixels [2002.00691]. The corresponding blending relation is
\[
\delta_s=\frac{\delta_{\mathrm{obs}}}{X_s},
\]
so stars with \(\delta_s>1\) are too faint to generate the observed signal [2002.00691].

The transit and eclipse models are batman-based physical light-curve models with quadratic limb darkening, and the original implementation assumes circular orbits and fixed \(P_{\rm orb}\) for speed and tractability [2002.00691]. Likelihoods are combined with priors on orbital period, planet radius, mass ratio, background stellar population, and multiplicity. For unresolved background/foreground blends, simulated stars are drawn from TRILEGAL in a \(0.1\,\mathrm{deg}^2\) cone centered on the target [2002.00691].

High-resolution imaging enters through a weighting term on scenario likelihoods, suppressing models that require undetected companions inconsistent with the supplied contrast curve [2002.00691]. This mechanism became standard in later applications. Studies such as the validation of TOI-1221 b explicitly fed SOAR/HRCam and Gemini South/Zorro contrast curves into TRICERATOPS, while also using TIC and Gaia information, empirical background-star rates, and multiplicity priors [2209.13651]. Other works supplied Gaia DR3 source environments, propagated source positions, PRF/aperture maps, centroid constraints, and follow-up multi-band photometry to refine dilution and off-target probabilities [2606.15430] [2604.15816].

A central practical distinction in the literature is between photometry prepared for geometric parameter inference and photometry prepared for validation. In the validation of TOI-7701.01, PDCSAP was used for geometric radius inference, while re-normalized un-detrended SAP was deliberately supplied to TRICERATOPS so that field dilution and crowding would remain encoded in the data stream used for false-positive modeling [2606.15430]. This suggests that TRICERATOPS is often most informative when the photon environment is preserved rather than aggressively flattened away.

## 4. Decision metrics, thresholds, and workflow placement

The original FPP–NFPP decision plane defines three standard classes for 2-minute cadence TESS data: validated planets with \(\mathrm{FPP}<0.015\) and \(\mathrm{NFPP}<10^{-3}\), likely planets with \(\mathrm{FPP}<0.5\) and \(\mathrm{NFPP}<10^{-3}\), and likely nearby false positives with \(\mathrm{NFPP}>10^{-1}\) [2002.00691]. The same thresholds are widely reused in later target-by-target validation papers [2209.13651] [2602.14840] [2604.15816]. The original study also states that performance is best for \(\mathrm{SNR}>15\), and gives a typical runtime of about 5 minutes per target on a standard 2-core laptop [2002.00691].

In practice, TRICERATOPS is rarely used in isolation. A common workflow is front-end vetting followed by probabilistic validation. One systematic hot-Neptune study first used DAVE for visual and centroid-level vetting, then ran TRICERATOPS on the surviving candidates, and later re-ran it with follow-up imaging constraints; ARIADNE provided refined stellar parameters and juliet provided physically consistent transit fits [2211.08490]. VaTEST II combined LATTE, TESS-Plot, high-resolution imaging, and TRICERATOPS to validate 11 K-star planets, while VaTEST III pre-cleared nearby stars with ground-based photometry and then set \(\mathrm{NFPP}=0\) before running TRICERATOPS on eight super-Earth candidates [2301.09865] [2311.00688].

Several studies report repeated runs to mitigate stochasticity. The hot-Neptune survey used 10 runs in its initial pass and 30 runs in its final validation stage [2211.08490]. VaTEST II used 15 iterations per target per contrast-curve file [2301.09865]. The upgraded multicolor framework reports mean values and 68% confidence intervals over 10 runs, noting the higher computational load of multiband likelihood evaluation [2508.02782]. Target-specific papers also used 20 independent runs to test robustness against MCMC stochasticity [2602.14840] [2512.10007].

The threshold values are not immutable. In the hot-Neptune desert, one study adopted a stricter validation criterion of \(\mathrm{FPP}<3\times10^{-3}\) and \(\mathrm{NFPP}<10^{-3}\), arguing that the default TRICERATOPS planet occurrence priors treat period and radius as independent and can therefore underestimate FPP in a sparse-occurrence region [2211.08490]. VaTEST III, by contrast, validated planets at \(\mathrm{FPP}\le1\%\) with \(\mathrm{NFPP}<0.1\%\) after pre-clearing nearby stars [2311.00688].

## 5. TRICERATOPS+ and the multicolor extension

TRICERATOPS+ is an upgraded version that incorporates ground-based light curves in separate bandpasses together with TESS photometry and high-resolution imaging [2508.02782]. Its core modification is the total likelihood,
\[
\ln \mathcal{L}_{\mathrm{total}}=\ln \mathcal{L}_{\mathrm{TESS}}+\sum_{\lambda}\ln \mathcal{L}_{\lambda}+\ln \mathrm{prior}_{\mathrm{companion}},
\]
which allows the framework to use chromatic information directly in model comparison [2508.02782]. The approach exploits the approximate achromaticity of a planetary transit, aside from weak limb-darkening dependence, against the intrinsically chromatic depths expected from many eclipsing-binary blends [2508.02782].

The bandpass-specific dilution relation is written as
\[
\delta_{\mathrm{obs},\lambda}=\delta_{\mathrm{true}}(\lambda)\,\frac{F_{t,\lambda}}{F_{t,\lambda}+F_{c,\lambda}},
\]
and TRICERATOPS+ supplements the original TESS-only flux-ratio machinery with bandpass-specific stellar limb darkening from ExoTiC-LD and flux-ratio splines in SDSS \(g,r,i,z\) and 2MASS \(J,H,K\), constructed using a TRILEGAL synthetic population at 100 pc [2508.02782]. In this upgraded framework, users can also supply the fractional flux of resolved contaminants in ground-based apertures [2508.02782].

The 2025 multicolor validation study applied TRICERATOPS+ to 14 TESS candidates and statistically validated six new planets in five systems: TOI-1346 b, TOI-1346 c, TOI-2719 b, TOI-4155 b, TOI-6000 b, and TOI-6324 b [2508.02782]. The same paper showed that multicolor photometry can either lower FPP and drive \(\mathrm{NFPP}\rightarrow0\), as for TOI-6000 b, or reveal chromaticity indicative of false positives, as for TOI-1254.01, TOI-1616.01, TOI-4051.01, and TOI-5706.01 [2508.02782].

A system-specific demonstration appears in the validation of TOI-2267 d in a \(0.\!^{\prime\prime}384\) M-dwarf binary. There, TRICERATOPS+ combined a stacked TESS profile, two diffuser-assisted Hale/Palomar WIRC \(J\)-band transits, and speckle-imaging contrast curves. The reported result was \(\mathrm{FPP}=4.7\times10^{-6}\) with a 68% range of \((4.5\times10^{-7},5.0\times10^{-5})\) and \(\mathrm{NFPP}=0\), while the host-star assignment remained ambiguous despite transit-shape stellar-density analysis [2512.10007].

## 6. Empirical behavior, strengths, and limitations

Published applications illustrate how strongly TRICERATOPS responds to contamination information. In TOI-4616 b, an initial TESS-only run on detrended data gave \(\mathrm{FPP}=0.823\) and \(\mathrm{NFPP}=0.74\), dominated by nearby eclipsing-binary and nearby transiting-planet scenarios associated with a neighbor in the TESS aperture. Repeating the analysis with uncontaminated ground-based photometry and a NESSI contrast curve removed those pathways, yielding a final \(\mathrm{FPP}=0.0135\) and \(\mathrm{NFPP}=0\), sufficient for validation under the standard criterion [2603.10905]. In TOI-1221 b, the framework returned \(\mathrm{FPP}=0.0014\pm0.0003\) and \(\mathrm{NFPP}=0.0012\pm0.0004\); the authors emphasized that these values were conservative because TRICERATOPS did not ingest reconnaissance RV constraints or nearby-eclipsing-binary clearances from ground-based photometry [2209.13651].

Other studies show the same pattern in more favorable environments. TOI-1752 b and c were both validated with low FPP and essentially zero NFPP after combining TESS photometry, Gaia DR3 information, AO contrast curves, SPOC centroid diagnostics, and multicolor follow-up [2604.15816]. For TOI 7475.01, a workflow based on natural-flux-preserving light-curve preparation, BLS detection, Gaia DR3 contamination checks, and centroid stability produced 20 TRICERATOPS runs with mean \(\mathrm{FPP}=0.000000\pm10^{-6}\) and \(\mathrm{NFPP}=0.000000\pm0.0\); the posterior mass was concentrated entirely in TP, PTP, and DTP, giving a cumulative planetary probability in the TIC 376866659 system of approximately \(100\%\) [2602.14840]. For TOI-7701.01, ensemble TRICERATOPS analysis on re-normalized SAP yielded \(\mathrm{FPP}=0.00191\pm0.00247\) and \(\mathrm{NFPP}<10^{-6}\), while the dominant TP scenario converged to \(R_p=7.86\,R_\oplus\), close to the PDCSAP-derived geometric radius of \(8.07\,R_\oplus\) [2606.15430].

The strengths of the framework are therefore clear in the literature: explicit nearby-source modeling, direct use of Gaia/TIC environments, natural integration of contrast curves, and a quantitative NFPP that is especially relevant for TESS’s large pixels [2002.00691]. The multicolor extension adds a direct chromaticity test that is often decisive for crowded or blended systems [2508.02782].

The limitations are equally consistent across studies. TRICERATOPS focuses on astrophysical scenarios rather than instrumental/systematic false alarms [2002.00691] [2508.02782]. The original implementation assumes circular orbits and depends on accurate stellar parameters, catalog completeness, and the Gaussian PRF approximation [2002.00691]. Several papers note that results depend on priors and on the completeness of the background-star model, and that some external constraints—especially RV limits and certain neighbor clearances—may not be ingested directly, which can make the reported FPP conservative or, in occurrence-sparse regimes, optimistic unless stricter thresholds are imposed [2209.13651] [2211.08490] [2602.14840]. This suggests that TRICERATOPS is most robust when used as the probabilistic core of a broader vetting architecture that includes centroid diagnostics, high-resolution imaging, and targeted ground-based photometry.

Source: https://www.emergentmind.com/topics/triceratops