---
title: First Attempt Time Algorithm (FATA)
url: https://www.emergentmind.com/topics/first-attempt-time-algorithm-fata
type: topic
---

# First Attempt Time Algorithm (FATA)

The First Attempt Time Algorithm (FATA) is an acceptance–rejection approach for unbiased simulation of the first-passage time (FPT) in stochastic processes, particularly for Itô diffusions and fractional Brownian motion, to constant or time-dependent boundaries. FATA eliminates the bias associated with time-discretization by simulating candidate first-passage events under a simpler reference process (standard Brownian motion) and correcting via a likelihood ratio derived from Girsanov's theorem. Extensions exist for non-Markovian systems such as fractional Brownian motion, employing adaptive bisection and conditional sampling to control error and computational cost. FATA is exact, flexible, and particularly well-suited when high numerical precision is required, as in computational neuroscience or fine-grained stochastic modeling [2412.13060, 1908.11634].

## 1. Problem Formulation and Mathematical Framework

FATA addresses the exact simulation of the first-passage time
\[
\tau = \inf \{ t \geq 0 : X_t = L(t) \}
\]
for a one-dimensional diffusion $X$ governed by the Itô SDE
\[
dX_t = \alpha(X_t)\,dt + dW_t, \quad X_0 = x_0,
\]
to a (possibly time-dependent) boundary $L : [0,\infty) \to \mathbb{R}$. For many systems, discretizing either the path $X$ or the boundary $L$ introduces bias in the estimate of $\tau$. FATA circumvents this by transforming the simulation problem into an acceptance–rejection process over candidate FPTs drawn from the law of a reference process (most often standard Brownian motion) to the target boundary [2412.13060].

The algorithm can also be adapted to fractional Brownian motion (fBm) and general Gaussian processes, where the Markov property no longer holds. In these systems, adaptive refinement of the simulation mesh is used to retain accuracy and efficiency [1908.11634].

## 2. Core Algorithmic Structure

FATA for Itô diffusions combines candidate generation and correction:

1. **Candidate generation**: Sample a candidate $\tau_B$ from the known first-passage distribution of a standard Brownian motion $W_t$ to $L(t)$. For constant and linear boundaries, closed forms exist, e.g., for $L(t)\equiv \ell$,
   \[
   f_{\tau_B}(t) = \ell (2\pi t^3)^{-1/2} \exp(-\ell^2/2t),\quad t>0.
   \]
   For general $L(t)$, iterative schemes such as the Herrmann–Tanré approach are employed [2412.13060].
   
2. **Girsanov correction**: Accept $\tau_B$ as the FPT for the original drifted process with probability
   \[
   \alpha(\tau_B) = \eta(\tau_B) = \mathbb{E}_Q\left[\exp\left\{-\int_0^{\tau_B} [\gamma_1(s) + \gamma_2(X^Q_s)] ds\right\}\,\bigg|\, \tau_B\right],
   \]
   where $Q$ denotes the Wiener measure and
   \begin{align*}
   \gamma_1(s) &= -A'(L(s)) L'(s),\\
   \gamma_2(x) &= \frac{1}{2}[\alpha'(x)+\alpha^2(x)],
   \end{align*}
   with $A'(x) = \alpha(x)$. $X^Q$ is a Brownian bridge conditioned to start at $x_0$ and hit $L(\tau_B)$ at time $\tau_B$ [2412.13060].

For fBm and other Gaussian processes, FATA (also termed Adaptive Bisection Sectioning, "ABSec") uses a two-phase procedure: initialize on a coarse grid via the Davies–Harte algorithm, truncate beyond the first crossing, and recursively bisect "critical" intervals by sampling conditional midpoints. Refinement continues until the probability of missing a crossing on any interval is bounded by a user-specified $\epsilon'$ [1908.11634].

## 3. Extensions to Time-Dependent and Nonlinear Boundaries

FATA naturally handles time-dependent boundaries by updating the candidate generation step to sample from the FPT of a Brownian motion to a moving boundary $L(t)$. For boundaries that are not piecewise linear or constant, iterative inverse-Gaussian schemes are applied. The acceptance probability and related Girsanov correction terms are recomputed, with the only requirement being that $L(s)$ and its derivative $L'(s)$ are available for computing $\gamma_1(s)$ [2412.13060].

For fBm, the approach adapts to the loss of the Markov property by replacing candidate generation and acceptance steps with adaptive mesh refinement. Bridges that could conceal a crossing are recursively bisected, and at each bisection, conditional distributions are updated and crossing probabilities are re-evaluated. Error per bridge is at most $\epsilon'$, and overall false-miss probability $P^{\rm tot}_{\rm error}$ can be set arbitrarily small with complexity increasing only logarithmically in $1/\epsilon'$ [1908.11634].

## 4. Computational Complexity and Efficiency

For Itô diffusions, the expected number of FATA candidate draws needed for acceptance is
\[
\mathbb{E}[I] = \exp\{-A(x_0)+A(L(0))\} = 1/\mathbb{E}_Q[\eta(\tau_B)],
\]
and is exponentially sensitive to the worst-case value $\kappa = \sup_{0\leq s\leq\tau_B}[\gamma_1(s)+\gamma_2(\cdot)]$. An exponential upper bound follows:
\[
\mathbb{E}[I] \leq \mathbb{E}_Q[\exp(\kappa\tau_B)].
\]
To reduce cost when $\kappa$ is large, techniques such as shifting $\gamma_1,\gamma_2$ by their mininum or "space splitting"—dividing the boundary into $k$ sublevels—can bring costs to $O(k)$ per run [2412.13060].

For fBm, FATA attains an effective computational complexity of $O((\ln N_{\text{eff}})^3)$ for time resolution $N_{\text{eff}}$, compared to $O(N\ln N)$ for the Davies–Harte method on a uniform grid. Empirically, for $H=0.33$ and $N_{\text{eff}}=2^{32}$, FATA is $5,000\times$ faster and $10^4\times$ more memory efficient. With $H=0.25$ and $N_{\text{eff}}=2^{42}$, the speedup is $3\cdot 10^5$ and memory cost $10^6$ times smaller. The complexity increases only logarithmically with decreasing target global error $P^{\rm tot}_{\rm error}$ [1908.11634].

| Algorithm     | Time Complexity               | Memory Complexity | Bias Control         |
|---------------|------------------------------|-------------------|---------------------|
| FATA (diffusion)    | Random (potentially $\exp(\kappa)$) | Problem-dependent  | No bias (exact law) |
| FATA (fBm/ABSec)    | $O((\ln N_{\text{eff}})^3)$         | $O(\ln N_{\text{eff}})$ | Tunable, bounded error |
| Davies–Harte  | $O(N\ln N)$                  | $O(N)$            | Discretization bias |

## 5. Numerical Results and Empirical Behavior

Representative numerical experiments for FATA have been reported for diverse SDEs and boundary types:

- **Sinusoidal drift SDE with linear boundary**: FATA produces unbiased FPT samples. Euler–Maruyama and improved EM methods exhibit $O(1)$ or $O(\sqrt{\Delta})$ bias that persists unless exceedingly fine discretization is used; FATA achieves exactness regardless of $\Delta$ [2412.13060].
- **SDE with exponential boundary**: When explicit sampling is unavailable, iterative inverse-Gaussian schemes allow FATA to converge to the true density as tolerance $\epsilon\to 0$, with run-time scaling moderately as $\epsilon$ decreases [2412.13060].
- **Neuronal spiking models**: For a quadratic leaky integrate-and-fire model with adaptive threshold $L(t)$, FATA produces exact spike-time samples, including under nonlinear SDEs for which EM can fail (blow-up). Reported performance is $O(1\,\mathrm{s})$ per $10^4$ spikes, with first and second moments matching theoretical predictions [2412.13060].
- **Fractional Brownian motion**: For $H=0.33$ and effective grid size $2^{32}$, Davies–Harte is infeasible (beyond $2^{24}$), while FATA requires $2\,\mathrm{ms}$ per sample and $30\,\mathrm{MB}$ RAM. Similar drastic improvements are reported for more subdiffusive $H=0.25$ and even higher resolutions [1908.11634].

## 6. Limitations and Remedies

The principal limitation of FATA for drifted diffusions is the potentially large exponential cost with respect to the maximal Girsanov correction $\kappa$ and large candidate FPTs $\tau_B$. Remedies include recentering correction terms or partitioning the boundary range to reduce per-step rejection rates [2412.13060]. For fBm and adaptive bisection, double-precision subtraction when updating variances leads to floating-point limitations; this manifests as practical bounds $N_{\rm eff} \lesssim 2^{10.5/H}$, with extended precision arithmetic offering a solution [1908.11634].

## 7. Comparison with Alternative Approaches

Conventional methods such as Euler–Maruyama, improved EM, and Davies–Harte are fundamentally limited by time-discretization bias and computational cost scaling linearly with target resolution or step size. In contrast, FATA (for both Markovian and non-Markovian systems) achieves either unbiased or provably bounded-error FPT estimates with substantially reduced computational resource requirements at fine resolutions. For applications demanding extreme precision or unbiased estimates—such as theoretical validation, rare event simulation, or computational neuroscience—FATA provides strong advantages, provided the potential for exponential algorithmic cost is managed for the specific problem instance [2412.13060, 1908.11634].

Source: https://www.emergentmind.com/topics/first-attempt-time-algorithm-fata