---
title: True Self-Avoiding Walk (TSAW) Overview
url: https://www.emergentmind.com/topics/true-self-avoiding-walk-tsaw
type: topic
---

# True Self-Avoiding Walk (TSAW) Overview

True self-avoiding walk (TSAW), also called the myopic or “true” self-avoiding random walk, is a non-Markovian nearest-neighbor walk whose transition mechanism depends on the accumulated occupation profile of the path, so that motion is biased toward less visited regions. In lattice formulations this bias is expressed through site local times or edge traversal counts; in continuum limits it is related to true self-repelling motion; and in algorithmic settings it appears as an irreversible, self-repelling sampling dynamics. The rigorous picture is strongly dimension-dependent: in \(d\ge 3\) one has diffusion for broad classes of self-interaction functions and a Brownian finite-dimensional scaling limit in a restricted quartic case, whereas in \(d=1\) the process is super-diffusive with \(t^{2/3}\) scaling and converges, after rescaling, to the Tóth–Werner true self-repelling motion [1004.4262].

## 1. Definitions and principal formulations

In the lattice model studied by Horváth, Tóth, and Vető, \(X(t)\in \mathbb Z^d\) is a continuous-time nearest-neighbor jump process with occupation-time field
\[
L(t,x)=L(0,x)+\int_0^t \mathbf 1_{\{X(s)=x\}}\,ds.
\]
A smooth rate function \(w:\mathbb R\to(0,\infty)\) is fixed, with ellipticity \(\inf_{u\in\mathbb R} w(u)>0\), and decomposed into even and odd parts
\[
s(u)=\tfrac12[w(u)+w(-u)],\qquad r(u)=\tfrac12[w(u)-w(-u)].
\]
Conditionally on \(X(t)=x\), the walker jumps to a nearest neighbor \(y\) at rate \(w(L(t,x)-L(t,y))\). In the continuum heuristic, the drift is proportional to the negative gradient of the occupation measure, so the walk is repelled from frequently visited sites [1004.4262].

The one-dimensional literature also uses a site-local-time formulation in which the jump rates out of \(x\) are
\[
\text{rate}\{x\to x\pm1\}=\exp\bigl(-\beta\,L_t(x\pm1)\bigr),
\]
or, in discrete time,
\[
P\bigl(X_{n+1}=x\pm1\mid\mathcal F_n\bigr)
=\frac{\exp\bigl(-\beta\,N_n(x\pm1)\bigr)}
{\exp\bigl(-\beta\,N_n(x-1)\bigr)+\exp\bigl(-\beta\,N_n(x+1)\bigr)},
\]
where \(N_n(y)\) is the visit count at \(y\). In the limit \(\beta\to\infty\), this becomes a strictly self-avoiding walk in which only new sites are entered when possible [2410.08694].

A distinct one-dimensional formulation uses edge local times. Writing \(\mathsf E_n^\pm(x)\) for the number of right- and left-directed traversals of adjacent edges, the transition probabilities are weighted by past edge usage, so that heavily used edges become less attractive. A related edge-repulsion model on infinite locally finite trees assigns weight \(w(n)=e^{-\beta n}\) to an edge after \(n\) traversals, and from a current vertex \(x\) the next step to a neighbor \(y\) is chosen with probability proportional to \(e^{-\beta \mathsf L_n(x,y)}\) [2502.10960].

## 2. Environment process and stationary Gibbs structure

A central device in the rigorous analysis is the environment as seen from the walker. Define
\[
\eta(t,z)=L\bigl(t,X(t)+z\bigr),\qquad z\in\mathbb Z^d.
\]
Then \(\eta(t)\) is itself a Markov process on a configuration space \(\Omega\), and its generator acts on smooth cylinder functions \(f\) by
\[
Gf(\omega)=\partial_0 f(\omega)+\sum_{|e|=1} w\bigl(\omega(0)-\omega(e)\bigr)\,[f(\tau_e\omega)-f(\omega)],
\]
where \(\tau_e\omega(x)=\omega(x+e)\). This converts the original non-Markovian walk into a Markovian environment process, at the price of enlarging the state space [1004.4262].

The invariant law is a translation-invariant Gibbs measure \(\pi\). Writing
\[
R(u)=\int_0^u r(v)\,dv,
\]
the finite-volume conditional density of \(\pi\) on \(\omega_\Lambda\) given \(\omega_{\Lambda^c}\) is proportional to an exponential involving nearest-neighbor differences \(R(\omega(x)-\omega(y))\). In \(d\ge 3\), strict convexity of \(R\), expressed in the cited work by \(r'(\cdot)>c>0\), guarantees existence and uniqueness of \(\pi\), and \(\pi\) is stationary and ergodic for the environment process [1004.4262].

The special linear case \(r(u)=u\), for which \(R(u)=u^2/2\), is structurally distinguished. Then \(\pi\) becomes the massless free Gaussian field on \(\mathbb Z^d\), with covariance
\[
E_\pi[\omega(x)\omega(y)]=(-\Delta)^{-1}_{x,y}=:b(x-y).
\]
This Gaussian setting is the one in which the most complete central limit theorem is obtained [1009.0401].

## 3. Diffusion in \(d\ge 3\)

For the non-recurrent dimensions \(d\ge 3\), the stationary environment immediately yields zero asymptotic speed: for \(\pi\)-almost every initial environment, almost surely
\[
X(t)/t\to 0.
\]
The displacement nevertheless has diffusive order. For every unit vector \(e\in\mathbb R^d\), there exist constants \(0<\gamma\le \Gamma<\infty\) such that
\[
\gamma \le \liminf_{t\to\infty} t^{-1}E[(e\cdot X(t))^2]
\le \limsup_{t\to\infty} t^{-1}E[(e\cdot X(t))^2]
\le \Gamma.
\]
The proof proceeds by decomposing the coordinate process into orthogonal martingales and compensators, then controlling the compensators through variational estimates, infrared bounds, and Brascamp–Lieb inequalities [1004.4262].

A more explicit hypothesis class is given by the conditions
\[
r'(u)\ge c>0,\qquad s(u)\le C\exp\{(c-\varepsilon)u^2/2\},
\]
together with an entire-function growth condition on \(r\). Under these assumptions, the cited work proves finite, positive upper and lower diffusive bounds for every coordinate direction [1009.0401].

The full finite-dimensional central limit theorem is proved in a restricted quartic case:
\[
r(u)=u,\qquad s(u)=s_4u^4+s_2u^2+s_0,
\]
with \(|s_4|\) sufficiently small. Then the asymptotic covariance matrix
\[
\Sigma_{ij}=\lim_{t\to\infty} t^{-1}E[X_i(t)X_j(t)]
\]
exists and is non-degenerate, and the rescaled process
\[
X^N(t)=N^{-1/2}X(Nt)
\]
has finite-dimensional distributions converging to those of a \(d\)-dimensional Brownian motion with covariance \(\Sigma\) [1004.4262].

These results rigorously confirm part of the renormalization-group predictions of Amit, Parisi, and Peliti for \(d\ge 3\): standard diffusive scaling, Gaussian finite-dimensional limits in the quartic subcase, and \(t\)-linear growth of the mean-square displacement under mild assumptions on the rate function [1009.0401].

## 4. Non-reversible Kipnis–Varadhan theory and the graded-sector method

The rigorous route to the central limit theorem is non-reversible Kipnis–Varadhan theory applied to the environment process. In \(L^2(\Omega,\pi)\), the generator is split as
\[
S=-\tfrac12(G+G^*),\qquad A=\tfrac12(G-G^*),
\]
so that \(G=-S+A\), with \(S\ge 0\) self-adjoint and \(A\) skew-adjoint. The reversible part is carried by \(-S\), while the asymmetry of TSAW is encoded by \(A\) [1004.4262].

For a centered observable \(f\), let \(u_\lambda=(\lambda-G)^{-1}f\). The resolvent form of Kipnis–Varadhan theory requires
\[
\lim_{\lambda\to 0}\lambda^{1/2}\|u_\lambda\|=0,\qquad
\lim_{\lambda\to 0}\|S^{1/2}u_\lambda\|<\infty.
\]
When these conditions hold, the additive functional \(\int_0^T f(\eta(s))\,ds\) admits an efficient martingale approximation, and a central limit theorem follows [1004.4262].

The main technical difficulty is verifying the resolvent criteria in a non-reversible, graded state space. The cited analysis decomposes \(L^2(\pi)=\bigoplus_{n=0}^\infty H_n\) according to the degree of the environment observable and introduces a diagonal operator \(D\) that acts approximately as \(n\) on \(H_n\). The graded-sector condition then requires operator-norm control on grade-preserving and grade-changing components such as
\[
\|D^{-1/2}(S_{n,n}+A_{n,n})D^{-1/2}\|=O(n^\kappa),\qquad
\|D^{-1/2}A_{n,n\pm j}D^{-1/2}\|=O(n^\kappa).
\]
Horváth, Tóth, and Vető slightly weaken the original Sethuraman–Varadhan–Yau graded-sector condition, and in the Gaussian subcase explicit Fock-space operator estimates make the program verifiable [1004.4262].

This framework is significant because it isolates the essential analytic structure of TSAW in high dimension: the process is non-reversible, but it is still sufficiently controlled by the reversible part of the environment dynamics to permit a Brownian limit theorem.

## 5. One-dimensional scaling and the true self-repelling motion

In one dimension, TSAW is super-diffusive rather than diffusive. Rigorous results summarized by Maggs state that as \(t\to\infty\), the displacement \(X(t)\) and endpoint local time \(H(t)=N_t(X(t))\) satisfy
\[
\rho_1(t,x)=\mathbb P(X(t)=x)=t^{-2/3}\nu_1(u)\big|_{u=x\,t^{-2/3}},
\]
\[
\rho_2(t,h)=\mathbb P(H(t)=h)=t^{-1/3}\nu_2(v)\big|_{v=h\,t^{-1/3}},
\]
where \(\nu_1\) and \(\nu_2\) have explicit Airy-type integral representations. In particular,
\[
\langle X(t)^2\rangle^{1/2}\sim C\,t^{2/3},\qquad
\langle H(t)\rangle\sim C' t^{1/3},
\]
and
\[
X(t)/t^{2/3}\xrightarrow[t\to\infty]{\rm law} U,\qquad
H(t)/t^{1/3}\xrightarrow[t\to\infty]{\rm law} V,
\]
with densities \(\nu_1\) and \(\nu_2\). In the associated dynamic-exponent notation \(t\sim |x|^{1+z}\), one has \(z=1/2\). The cited discussion places the model in a super-diffusive universality class with several KPZ-like features, while emphasizing that its scaling functions differ from the GUE Tracy–Widom form [2410.08694].

A complete functional scaling limit is now available. Kosygina and Peterson prove that
\[
Y_t^n=n^{-2/3}X_{\lfloor nt\rfloor},\qquad
H_t^n=n^{-1/3}L\bigl(\lfloor nt\rfloor,X_{\lfloor nt\rfloor}\bigr)
\]
converge jointly in \(D([0,\infty),\mathbb R^2)\), with the Skorokhod \(J_1\)-topology, to the Tóth–Werner true self-repelling motion together with its local time at the current position. Their proof is organized around a joint generalized Ray–Knight theorem for rescaled local-time curves, absorption points, and merge points, followed by tightness and finite-dimensional convergence arguments [2502.10960].

This one-dimensional theory shows that the phrase “diffusive limit” is dimension-sensitive for TSAW. In \(d=1\), the correct scaling is \(n^{2/3}\), the limit is non-Gaussian, and the local-time field is part of the limiting object rather than a secondary observable.

## 6. Dimensional dependence, geometric phase structure, and rare-event geometry

The asymptotic regime depends sharply on the underlying geometry.

| Regime | Main behavior | Source |
|---|---|---|
| \(\mathbb Z^d\), \(d\ge 3\) | Diffusive bounds; Brownian finite-dimensional limit in a restricted quartic case | [1004.4262] |
| \(\mathbb Z^2\) | Expected \(E|X(t)|^2\sim t\log t\); \(\pi\) does not exist without an infrared cutoff; only partial non-rigorous and anisotropic results are known | [1004.4262] |
| \(\mathbb Z\) | Super-diffusive \(t^{2/3}\) scaling; convergence to true self-repelling motion after rescaling | [2502.10960] |
| Infinite locally finite trees | Recurrent if \(\mathrm{br}_r(\mathcal T)<1/2\), transient if \(\mathrm{br}_r(\mathcal T)>1/2\) | [2604.24389] |

For trees, the asymptotic criterion is formulated through the branching-ruin number
\[
\mathrm{br}_r(\mathcal T)=\sup\Bigl\{\gamma>0:\inf_{\pi\in\Pi}\sum_{e\in\pi}|e|^{-\gamma}>0\Bigr\},
\]
where \(\Pi\) is the collection of minimal cut-sets. The boundary \(\partial\mathcal T\), equipped with the metric
\[
d(\xi,\eta)=1/\bigl|\hbox{last common edge of }\xi\hbox{ and }\eta\bigr|,
\]
has Hausdorff dimension equal to \(\mathrm{br}_r(\mathcal T)\). The cited theorem shows a sharp phase transition: the TSAW on \(\mathcal T\) is recurrent when \(\mathrm{br}_r(\mathcal T)<1/2\) and transient when \(\mathrm{br}_r(\mathcal T)>1/2\). The proof uses Rubin’s construction, one-dimensional ruin probabilities with asymptotic \(r_n\sim Cn^{-1/2}\), and a quasi-independent ruin percolation argument in the sense of Lyons [2604.24389].

Two-dimensional rare-event geometry reveals an additional distinction between bulk and tail behavior. Schawe and Hartmann study the convex-hull area and perimeter of a TSAW on the square lattice with transition probabilities
\[
p_i=\frac{\exp(-\beta n_i)}{\sum_{j\in\mathcal N}\exp(-\beta n_j)}.
\]
Their simulations confirm bulk scaling with
\[
r\sim T^\nu(\ln T)^\alpha,\qquad \nu=\tfrac12,\ \alpha=\tfrac14,
\]
but the large-deviation rate function for the convex-hull area behaves as
\[
\Phi(s)\propto s^\kappa,\qquad \kappa\approx 1.17\pm 0.03,
\]
which is inconsistent with the naïve prediction \(\kappa=1\) and instead agrees with \(7/6\), the value associated in the paper with the Smart Kinetic SAW exponent \(\nu=4/7\). The qualitative explanation offered there is that in the far-right tail, traps and loops are strongly suppressed, so the walk behaves effectively like a loop-free growth process even though the bulk is governed by the TSAW exponent \(\nu=1/2\) [1808.10698].

## 7. Algorithmic realizations and sampling applications

TSAW has also become a mechanism for irreversible sampling algorithms. In the one-dimensional event-chain Monte Carlo (ECMC) setting, Maggs and collaborators identify a microscopic dynamics whose large-scale law coincides with the TSAW. The algorithm uses factorized Metropolis updates: one active particle moves continuously to the right at unit speed until it collides with a neighbor, collision times are determined by
\[
\Delta V^+(r)=-T\ln(\mathrm{Uniform}(0,1)),
\]
and activity is then passed to the neighbor. No rejections occur, detailed balance is broken, and only global balance is preserved. Adding a factor field
\[
V(r)\mapsto V(r)+p\,r
\]
cancels the mean drift and exposes the TSAW scaling. For harmonic chains at zero net tension, the resulting dynamics reproduces the \(t^{2/3}\) and \(t^{1/3}\) laws for displacement and endpoint local time, and yields an autocorrelation decay exponent \(z=1/2\), compared in the cited summary with \(z\approx 1\) for classic ECMC and \(z=2\) for reversible Monte Carlo [2410.08694].

A different algorithmic extension places TSAW on a finite state space \(V\) with an irreducible base kernel \(P\) and stationary distribution \(\pi\). If
\[
N_t(i,j)=\sum_{s=0}^{t-1}\mathbf 1\{X_s=i,\ X_{s+1}=j\},\qquad
L_t(i)=\sum_{s=0}^{t-1}\mathbf 1\{X_s=i\},
\]
and
\[
\epsilon_t(i,j)=N_t(i,j)-L_t(i)P_{ij},
\]
then the adaptive transition rule is
\[
Q_t(i,j)=\frac{P_{ij}\exp\{-\lambda\,\epsilon_t(i,j)\}}
{\sum_{k\in V}P_{ik}\exp\{-\lambda\,\epsilon_t(i,k)\}}.
\]
For every state \(i\) and every edge \((i,j)\) with \(P_{ij}>0\), the empirical counts satisfy
\[
L_t(i)-t\pi_i=O(\sqrt{\log t}),\qquad
N_t(i,j)-t\pi_iP_{ij}=O(\sqrt{\log t})
\quad\text{almost surely}.
\]
Consequently, for every bounded \(f:V\to\mathbb R\),
\[
\left|\frac1t\sum_{s=0}^{t-1}f(X_s)-\sum_{i\in V}\pi_i f(i)\right|
=O\!\left(\frac{\sqrt{\log t}}{t}\right)
\quad\text{almost surely}.
\]
The cited paper contrasts this with the usual \(O(t^{-1/2})\) error scaling for empirical averages under standard random-walk-based MCMC and attributes the improvement to the self-avoidance mechanism, which suppresses local over-sampling and enforces near-balancing of empirical flow along each row. It also notes the practical cost of storing the full matrix of transition counts and identifies extensions to continuous or very large state spaces as open [2605.30532].

Across these settings, TSAW functions both as a probabilistic model of long-memory self-repulsion and as a design principle for non-reversible algorithms. The shared structural theme is the same: empirical overuse of sites or edges feeds back into future motion, but the resulting asymptotics range from Brownian diffusion in high-dimensional lattices to non-Gaussian \(t^{2/3}\) scaling in one dimension, geometry-controlled phase transitions on trees, and sharpened empirical-integration bounds on finite state spaces.

Source: https://www.emergentmind.com/topics/true-self-avoiding-walk-tsaw