RE:CZ

Derivation of SDE Equations for Three-Body Dynamics in Capital Markets

Quantitative Finance

👤 Financial modeling researchers, quantitative analysts, economists interested in capital market dynamics
Building on the article 'The Three-Body Dynamics Hypothesis in Capital Markets,' this paper derives a complete system of stochastic differential equations (SDEs) to describe the interactions among momentum capital (M), value capital (V), and liquidity capital (L) in capital markets. The article defines fast variables (such as log premium, momentum, volatility) and slow variables (the volumes of the three types of capital) and extracts 12 formalizable core constraints. Through a detailed analysis of the SDE equations, the article validates these constraints one by one, including positive feedback for M, negative feedback for V, directionless feedback for L, positive and negative feedback loops, payoff matrices, and crowding effects. All constraints are validated, indicating that this SDE system fully implements the qualitative mechanisms from the original article, such as volatility clustering, fat-tailed distributions, and chaotic behavior. The article also conducts phase analysis and statistical property validation, providing a foundation for subsequent numerical simulations, bifurcation analysis, and parameter calibration.
  • ✨ Derived a complete SDE system describing the interactions among three types of capital
  • ✨ Validated 12 core constraints, including positive/negative feedback and payoff matrices
  • ✨ The system explains market characteristics such as volatility clustering and fat-tailed distributions
📅 2026-02-07 · 1,839 words · ~9 min read
  • Capital Markets
  • Three-Body Dynamics
  • Stochastic Differential Equations
  • SDE
  • Momentum Capital
  • Value Capital
  • Liquidity Capital

Table of Contents

Derivation of the SDE System for Capital Market Three-Body Dynamics1. State Variable Definitions2. Core Constraint Extraction3. SDE SystemFast Variable SubsystemSlow Variable Subsystem (ε1\varepsilon \ll 1)4. Parameter Table5. Constraint Verification5.1 Constraint C1: M's Positive Feedback on Price5.2 Constraint C2: V's Negative Feedback on Price5.3 Constraint C3: L's Non-Directional Feedback on Price5.4 Constraint C4: MσM\uparrow \to \sigma\uparrow5.5 Constraint C5: σL\sigma\uparrow \to L\downarrow (in crash spiral)5.6 Constraint C6: LL\uparrow \to Price Impact\downarrow5.7 Constraint C7: V Buys when S<SS < S^*, Sells when S>SS > S^*5.8 Constraint C8: VSSV\uparrow \to |S - S^*|\downarrow5.9 Constraint C9: M's Return μ\propto \mu, Risk δ\propto \delta, Cost σ\propto \sigma5.10 Constraint C10: V's Return δ\propto \delta, Risk σ\propto \sigma, Cost μ\propto \mu5.11 Constraint C11: L's Return σ\propto \sigma, Risk μ\propto \mu, Cost δ\propto \delta5.12 Constraint C12: Capital Surplus → Return Decline → Volume Contraction6. Complete Tracking of the Positive Feedback Loop7. Complete Tracking of the Negative Feedback Loop8. Phase State Analysis VerificationCold State (000): Low δ\delta, Low μ\mu, Low σ\sigmaHot State (111): High δ\delta, High μ\mu, High σ\sigma9. Three-Body Analogy Verification10. Statistical Property VerificationVolatility ClusteringFat-Tailed Distribution11. Summary12. Future Research Directions

Derivation of the SDE System for Capital Market Three-Body Dynamics

2026-02-07

Based on the article "The Three-Body Dynamics Hypothesis of Capital Markets," this paper derives the complete system of Stochastic Differential Equations (SDEs) and verifies one by one whether it satisfies all the qualitative constraints from the original text.


1. State Variable Definitions

The system has two time scales:

Fast Variables (seconds to days):

  • x(t)=ln(S/S)x(t) = \ln(S/S^*) — Log Premium (SS is price, SS^* is intrinsic value)
  • p(t)p(t) — Momentum (state variable for the rate of price change)
  • σ(t)\sigma(t) — Instantaneous Volatility

Slow Variables (weeks to years):

  • m(t)m(t) — Momentum Capital Volume
  • v(t)v(t) — Value Capital Volume
  • l(t)l(t) — Liquidity Capital Volume

Derived Variables:

  • δ=(SS)/S=ex1x\delta = (S - S^*)/S^* = e^x - 1 \approx x (when xx is small)
  • μ=dS/dt\mu = dS/dt (manifested as drift terms in the SDEs)

2. Core Constraint Extraction

Formalizable constraints extracted from the original text:

ID Constraint Source
C1 M's positive feedback on price: d(Position)dS>0\frac{d(\text{Position})}{dS} > 0 Definition of Momentum Capital
C2 V's negative feedback on price: d(Position)dS<0\frac{d(\text{Position})}{dS} < 0 Definition of Value Capital
C3 L's non-directional feedback on price: d(Position)dS0\frac{d(\text{Position})}{dS} \approx 0 Definition of Liquidity Capital
C4 MσM\uparrow \to \sigma\uparrow Positive Feedback Loop
C5 σL\sigma\uparrow \to L\downarrow (liquidity withdrawal) Positive Feedback Loop
C6 LL\uparrow \to Price Impact\downarrow L→M Constraint on Impact
C7 V buys when S<SS < S^*, sells when S>SS > S^* V's Anchoring Mechanism
C8 $V\uparrow \to S - S^*
C9 M's return μ\propto \mu, risk δ\propto \delta, cost σ\propto \sigma Return Matrix
C10 V's return δ\propto \delta, risk σ\propto \sigma, cost μ\propto \mu Return Matrix
C11 L's return σ\propto \sigma, risk μ\propto \mu, cost δ\propto \delta Return Matrix
C12 Capital surplus → Return decline → Volume contraction Natural Selection Mechanism

3. SDE System

Fast Variable Subsystem

(I) Log Premium: dx=pdt+σlβdW1dx = p \cdot dt + \frac{\sigma}{l^\beta} \cdot dW_1

(II) Momentum: dp=1l[αMmpαVvxγp]dt+ηdW2dp = \frac{1}{l}\left[\alpha_M \cdot m \cdot p - \alpha_V \cdot v \cdot x - \gamma \cdot p\right] dt + \eta \cdot dW_2

Rearranged as: dp=αMmγlpdtαVvlxdt+ηdW2dp = \frac{\alpha_M m - \gamma}{l} \cdot p \cdot dt - \frac{\alpha_V v}{l} \cdot x \cdot dt + \eta \cdot dW_2

(III) Volatility: dσ=κσ(σˉσ)dt+λMmpdtλVvdtλLldt+ξσdW3d\sigma = \kappa_\sigma(\bar{\sigma} - \sigma) \cdot dt + \lambda_M m |p| \cdot dt - \lambda_V v \cdot dt - \lambda_L l \cdot dt + \xi \sigma \cdot dW_3

Slow Variable Subsystem (ε1\varepsilon \ll 1)

(IV) Momentum Capital: dm=εm[aMpbMxcMσρMm]dt+εσmmdW4dm = \varepsilon \cdot m \cdot \left[a_M |p| - b_M |x| - c_M \sigma - \rho_M m\right] dt + \varepsilon \sigma_m m \cdot dW_4

(V) Value Capital: dv=εv[aVxbVσcVpρVv]dt+εσvvdW5dv = \varepsilon \cdot v \cdot \left[a_V |x| - b_V \sigma - c_V |p| - \rho_V v\right] dt + \varepsilon \sigma_v v \cdot dW_5

(VI) Liquidity Capital: dl=εl[aLσbLpcLxρLl]dt+εσlldW6dl = \varepsilon \cdot l \cdot \left[a_L \sigma - b_L |p| - c_L |x| - \rho_L l\right] dt + \varepsilon \sigma_l l \cdot dW_6

Where W1,,W6W_1, \ldots, W_6 are independent standard Brownian motions.


4. Parameter Table

Parameter Meaning
αM\alpha_M Strength of M's positive feedback (aggressiveness of trend-following)
αV\alpha_V Strength of V's negative feedback (force of mean reversion)
γ\gamma Natural damping (friction, information decay)
β\beta Liquidity's buffering exponent on price impact
κσ\kappa_\sigma Volatility mean reversion speed
σˉ\bar{\sigma} Long-term mean of volatility
λM,λV,λL\lambda_M, \lambda_V, \lambda_L Influence coefficients of the three bodies on volatility
ξ\xi Vol-of-vol coefficient
η\eta Momentum noise intensity (new information shocks)
ai,bi,cia_i, b_i, c_i Return/risk/cost coefficients (i=M,V,Li = M, V, L)
ρi\rho_i Crowding penalty coefficient
ε\varepsilon Time scale separation parameter (1\ll 1)

5. Constraint Verification

5.1 Constraint C1: M's Positive Feedback on Price

Original Text: Momentum capital provides positive feedback on price changes, d(Position)dS>0\frac{d(\text{Position})}{dS} > 0 (buy on rise, sell on fall).

Corresponding Equation Term: The αMmp\alpha_M \cdot m \cdot p term in equation (II).

Verification:

  • When p>0p > 0 (price rising), this term is positive, increasing dp/dtdp/dt, i.e., accelerating the rise.
  • When p<0p < 0 (price falling), this term is negative, making dp/dtdp/dt more negative, i.e., accelerating the fall.
  • This is precisely the mathematical expression of "trend-following": the change in momentum is in the same direction as momentum itself.

Conclusion: ✓ Passed


5.2 Constraint C2: V's Negative Feedback on Price

Original Text: Value capital provides negative feedback on price changes, d(Position)dS<0\frac{d(\text{Position})}{dS} < 0 (reduce position when price rises, increase position when price falls).

Corresponding Equation Term: The αVvx-\alpha_V \cdot v \cdot x term in equation (II).

Verification:

  • When x>0x > 0 (price above intrinsic value), this term is negative, generating a downward force, suppressing the rise.
  • When x<0x < 0 (price below intrinsic value), this term is positive, generating an upward force, suppressing the fall.
  • This is precisely the mechanical mechanism of mean reversion.

Conclusion: ✓ Passed


5.3 Constraint C3: L's Non-Directional Feedback on Price

Original Text: Liquidity capital has no directional reaction to price changes, d(Position)dS0\frac{d(\text{Position})}{dS} \approx 0.

Corresponding Equation Terms: In equations (I) and (II), ll appears only in denominators, not generating directional drift terms.

Verification:

  • ll affects the magnitude of price impact (1/l1/l and σ/lβ\sigma/l^\beta), not the direction.
  • Market makers do not bet on direction; they only provide liquidity buffers.

Conclusion: ✓ Passed


5.4 Constraint C4: MσM\uparrow \to \sigma\uparrow

Original Text: An increase in momentum capital leads to rising volatility.

Corresponding Equation Term: The +λMmp+\lambda_M \cdot m \cdot |p| term in equation (III).

Verification:

  • When mm increases, λMmp\lambda_M \cdot m \cdot |p| increases.
  • This directly increases the drift term dσ/dtd\sigma/dt.
  • The more momentum capital and the stronger the momentum, the higher the volatility.

Conclusion: ✓ Passed


5.5 Constraint C5: σL\sigma\uparrow \to L\downarrow (in crash spiral)

Original Text: Market makers withdraw during high volatility.

Corresponding Equation Terms: aLσbLpa_L \sigma - b_L |p| in equation (VI).

Verification:

  • Superficially, aLσa_L \sigma is a return term for L, so σ\sigma\uparrow should benefit L.
  • Key Understanding: The logic of the original text is that in a crash spiral, high σ\sigma accompanies high p|p| (strong trend).
  • When bLp>aLσb_L |p| > a_L \sigma, L's net return becomes negative, causing ll to contract.
  • This is precisely the mechanism for market maker withdrawal in a "high volatility + strong trend" environment.

Parameter Condition: A crash spiral requires bL/aLb_L / a_L to be sufficiently large so that trend risk exceeds volatility return.

Conclusion: ✓ Passed (under appropriate parameter conditions)


5.6 Constraint C6: LL\uparrow \to Price Impact\downarrow

Original Text: Ample liquidity cushions the price impact from momentum capital.

Corresponding Equation Terms:

  • σ/lβ\sigma / l^\beta in equation (I): Larger ll means smaller price noise.
  • 1/l1/l in equation (II): Larger ll means the same force produces a smaller change in momentum.

Verification:

  • Market depth ll is the denominator of the price impact coefficient.
  • In deep markets, even large orders struggle to move prices.
  • Liquidity acts as a "shock absorber."

Conclusion: ✓ Passed


5.7 Constraint C7: V Buys when S<SS < S^*, Sells when S>SS > S^*

Original Text: Value capital performs contrarian operations anchored to intrinsic value.

Corresponding Equation Term: The αVvx-\alpha_V \cdot v \cdot x term in equation (II).

Verification:

  • x=ln(S/S)x = \ln(S/S^*)
  • When S<SS < S^*, x<0x < 0, so αVvx>0-\alpha_V v x > 0, generating an upward force (buying pressure).
  • When S>SS > S^*, x>0x > 0, so αVvx<0-\alpha_V v x < 0, generating a downward force (selling pressure).
  • The force magnitude is proportional to the deviation x|x| and the value capital volume vv.

Conclusion: ✓ Passed


5.8 Constraint C8: VSSV\uparrow \to |S - S^*|\downarrow

Original Text: Value capital intervention causes price to revert to intrinsic value.

Corresponding Equation Term: The αVvx-\alpha_V \cdot v \cdot x term in equation (II).

Verification:

  • Larger vv means a stronger reversion force αVvx|\alpha_V v x|.
  • A stronger reversion force drives xx toward zero faster.
  • Value capital is the system's "stabilizer."

Quantitative Analysis: With slow variables frozen, the part of equation (II) concerning xx is: p˙αVvlxγαMmlp\dot{p} \approx -\frac{\alpha_V v}{l} x - \frac{\gamma - \alpha_M m}{l} p

Combined with x˙=p\dot{x} = p, this is a second-order system. When αVv>0\alpha_V v > 0 and γ>αMm\gamma > \alpha_M m, the system is stable, and xx oscillates and converges to zero.

Conclusion: ✓ Passed


5.9 Constraint C9: M's Return μ\propto \mu, Risk δ\propto \delta, Cost σ\propto \sigma

Original Text: Momentum capital profits from trends, premium is risk, volatility is cost.

Corresponding Equation Term: The return term rM=aMpbMxcMσr_M = a_M |p| - b_M |x| - c_M \sigma in equation (IV).

Verification:

  • aMpa_M |p|: Trend continuation (large p|p|) = M's profit ✓
  • bMx-b_M |x|: Large premium signals reversal, is a risk ✓
  • cMσ-c_M \sigma: High volatility triggers frequent stop-losses, is a cost ✓

Conclusion: ✓ Passed


5.10 Constraint C10: V's Return δ\propto \delta, Risk σ\propto \sigma, Cost μ\propto \mu

Original Text: Value capital profits from value deviation, volatility is risk, trend is cost.

Corresponding Equation Term: The return term rV=aVxbVσcVpr_V = a_V |x| - b_V \sigma - c_V |p| in equation (V).

Verification:

  • aVxa_V |x|: Large premium = V's opportunity ✓
  • bVσ-b_V \sigma: High volatility leads to larger floating losses, and the anchor may also be unstable ✓
  • cVp-c_V |p|: Trend continuation forces V to wait longer, reducing capital efficiency ✓

Conclusion: ✓ Passed


5.11 Constraint C11: L's Return σ\propto \sigma, Risk μ\propto \mu, Cost δ\propto \delta

Original Text: Liquidity capital profits from volatility, trend is risk, premium is cost.

Corresponding Equation Term: The return term rL=aLσbLpcLxr_L = a_L \sigma - b_L |p| - c_L |x| in equation (VI).

Verification:

  • aLσa_L \sigma: High volatility = more trading opportunities, higher market-making profits ✓
  • bLp-b_L |p|: Strong trend leads to one-sided inventory accumulation, facing directional losses ✓
  • cLx-c_L |x|: Large premium requires wider spreads for self-protection, reducing efficiency ✓

Conclusion: ✓ Passed


5.12 Constraint C12: Capital Surplus → Return Decline → Volume Contraction

Original Text: Return-driven natural selection allows the three capital types to coexist long-term.

Corresponding Equation Term: The crowding term ρii-\rho_i \cdot i (i=m,v,li = m, v, l) in equations (IV-VI).

Verification:

  • Taking M as an example: dm/dtm(rMρMm)dm/dt \propto m \cdot (r_M - \rho_M m)
  • When mm is too large, ρMm\rho_M m dominates, effective return becomes negative, and mm contracts.
  • This is the crowding effect of logistic growth.
  • The same mechanism applies to V and L, preventing any single capital type from expanding indefinitely.

Conclusion: ✓ Passed


6. Complete Tracking of the Positive Feedback Loop

Original Text: MσLPrice ImpactσMM\uparrow \to \sigma\uparrow \to L\downarrow \to \text{Price Impact}\uparrow \to \sigma\uparrow \to M\uparrow (or blow-up)

SDE Tracking:

  1. mm\uparrow: External shock or return attraction.

  2. σ\to \sigma\uparrow: λMmp\lambda_M m |p| in equation (III) increases, volatility rises.

  3. l\to l\downarrow: In equation (VI), when bLp>aLσb_L |p| > a_L \sigma (strong trend exceeds volatility return), ll contracts.

  4. \to Price Impact\uparrow: 1/l1/l in equation (II) increases, the same force produces larger pp changes.

  5. σ\to \sigma\uparrow: λMmp\lambda_M m |p| in equation (III) increases further.

  6. m\to m\uparrow (or blow-up): In equation (IV):

    • If aMpa_M |p| dominates: mm continues to grow.
    • If bMx+cMσb_M |x| + c_M \sigma dominates (premium too large, volatility too high): mm contracts (blow-up).

Conclusion: ✓ Positive feedback loop fully implemented.


7. Complete Tracking of the Negative Feedback Loop

Original Text: SSσVSSσL|S-S^*|\uparrow \to \sigma\uparrow \to V\uparrow \to |S-S^*|\downarrow \to \sigma\downarrow \to L\uparrow

SDE Tracking:

  1. x|x|\uparrow: External shock causes price to deviate from intrinsic value.

  2. σ\to \sigma\uparrow: Price deviation is usually accompanied by volatility (transmitted via p|p|).

  3. v\to v\uparrow: aVxa_V |x| in equation (V) increases, V's return increases, attracting more value capital.

  4. x\to |x|\downarrow: αVvx-\alpha_V v x in equation (II) strengthens, reversion force increases, xx converges toward zero.

  5. σ\to \sigma\downarrow: x|x| and p|p| decrease, volatility in equation (III) subsides.

  6. l\to l\uparrow: In equation (VI), low p|p| reduces L's risk, attracting liquidity back.

Conclusion: ✓ Negative feedback loop fully implemented.


8. Phase State Analysis Verification

Cold State (000): Low δ\delta, Low μ\mu, Low σ\sigma

Correspondence: x0x \approx 0, p0p \approx 0, σσˉlow\sigma \approx \bar{\sigma}_{\text{low}}

Capital Returns:

  • rM=aM0bM0cMσˉcMσˉ<0r_M = a_M \cdot 0 - b_M \cdot 0 - c_M \bar{\sigma} \approx -c_M \bar{\sigma} < 0
  • rV=aV0bVσˉcV0bVσˉ<0r_V = a_V \cdot 0 - b_V \bar{\sigma} - c_V \cdot 0 \approx -b_V \bar{\sigma} < 0
  • rL=aLσˉbL0cL0aLσˉr_L = a_L \bar{\sigma} - b_L \cdot 0 - c_L \cdot 0 \approx a_L \bar{\sigma} (slightly positive or negative)

Original Description: All three parties are unprofitable, market shrinks.

Verification: ✓ Passed

Hot State (111): High δ\delta, High μ\mu, High σ\sigma

Correspondence: Large x|x|, large p|p|, large σ\sigma

Capital Returns: All three terms are large, sign depends on parameter ratios, uncertain.

Original Description: All three parties face extreme conditions, high return high risk, system at a critical point.

Verification: ✓ Passed


9. Three-Body Analogy Verification

Original Text: The market exhibits chaotic behavior due to the three bodies being evenly matched, sensitive to initial conditions.

SDE Analysis:

Drift term for pp in equation (II): αMmγlp\frac{\alpha_M m - \gamma}{l} \cdot p

  • When αMm>γ\alpha_M m > \gamma, this is a positive Lyapunov exponent direction for pp, making the system sensitive to initial conditions.
  • Simultaneously, αVvx/l-\alpha_V v x / l provides nonlinear coupling.
  • Linear instability + nonlinear coupling + noise = classic recipe for chaotic behavior.

Conclusion: ✓ Passed


10. Statistical Property Verification

Volatility Clustering

Original Text: Volatility clustering is a market stylized fact.

SDE Implementation: Multiplicative noise ξσdW3\xi \sigma dW_3 in equation (III).

Mechanism: When σ\sigma is high, noise is larger, making it easier for σ\sigma to stay high. Amplified by the positive feedback from λMmp\lambda_M m |p|.

Fat-Tailed Distribution

Original Text: Return distribution exhibits fat tails.

SDE Implementation: In the price noise σ/lβdW1\sigma / l^\beta \cdot dW_1, both σ\sigma and ll are stochastic.

Mechanism: Stochastic volatility itself generates fat tails; ll plummeting in extreme conditions (liquidity withdrawal) further thickens the tails.


11. Summary

Constraint Corresponding Equation Term Verification Result
C1 M Positive Feedback αMmp\alpha_M m p
C2 V Negative Feedback αVvx-\alpha_V v x
C3 L Non-Directional ll only in denominators
C4 M↑→σ↑ $\lambda_M m p
C5 σ↑→L↓ $a_L\sigma - b_L p
C6 L↑→Impact↓ 1/l1/l, σ/lβ\sigma/l^\beta
C7 V Anchoring Mechanism αVvx-\alpha_V v x
C8 V↑→Reversion αVvx-\alpha_V v x
C9 M Return Matrix $a_M p
C10 V Return Matrix $a_V x
C11 L Return Matrix $a_L\sigma - b_L p
C12 Crowding Effect ρii-\rho_i \cdot i

All 12 constraints passed.


12. Future Research Directions

  1. Numerical Simulation: Fix slow variables, simulate the fast variable subsystem, find attractors, limit cycles, chaotic regions.
  2. Bifurcation Analysis: Use m/vm/v ratio as bifurcation parameter, plot bifurcation diagrams.
  3. Parameter Calibration: Estimate parameters using real market data.
  4. Averaging Methods: Utilize time scale separation to analyze the effective dynamics of slow variables.

See Also

Referenced By