Realistic Backtesting: Order Book, VWAP & Almgren-Chriss

We design and develop full-cycle blockchain solutions: from smart contract architecture to launching DeFi protocols, NFT marketplaces and crypto exchanges. Security audits, tokenomics, integration with existing infrastructure.
Showing 1 of 1All 1305 services
Realistic Backtesting: Order Book, VWAP & Almgren-Chriss
Complex
~3-5 days
Frequently Asked Questions

Blockchain Development Services

Blockchain Development Stages

Latest works

  • image_website-b2b-advance_0.webp
    B2B ADVANCE company website development
    1358
  • image_web-applications_feedme_466_0.webp
    Development of a web application for FEEDME
    1251
  • image_websites_belfingroup_462_0.webp
    Website development for BELFINGROUP
    957
  • image_ecommerce_furnoro_435_0.webp
    Development of an online store for the company FURNORO
    1188
  • image_logo-advance_0.webp
    B2B Advance company logo design
    646
  • image_crm_enviok_479_0.webp
    Development of a web application for Enviok
    929

A strategy often performs wonderfully on a $10K demo account, but when moving to a real $1M, results degrade sharply. For example, a $500K trade on Uniswap might slip by 0.5%, while a $5M trade slips by 2%. The cause? Ignoring liquidity. Realistic backtesting that incorporates liquidity is the most complex yet critical aspect. Without it, you risk unrealistic results. Our team, with 10+ years of experience, builds backtesting engines that account for market liquidity and order impact on price. We've helped dozens of projects avoid losses from unaccounted slippage. Accurate modeling saves up to $15,000 per month in slippage for an average trader.

What is Realistic Backtesting with Order Book, VWAP & Almgren-Chriss?

Realistic backtesting that models liquidity is a method of evaluating trading strategies accounting for the impact of large orders on market price. Without such modeling, results can be overstated by 2-3 times. Our approach uses order book snapshots, VWAP models, and the Almgren–Chriss framework for accurate prediction of slippage and market impact. Every trade carries trading costs that are underestimated without liquidity-aware backtesting.

Problems Solved by Liquidity Modeling

A typical mistake: a backtest shows great results with small capital, but when scaling, the strategy stops working. For instance, on $10K, a trade is 0.01% of daily volume; on $1M, it's 1%, which noticeably affects price. The rule: an order exceeding 0.5–1% of daily volume starts moving the market. Without modeling, this leads to overestimated returns.

Another problem is misjudging slippage. Without a liquidity model, slippage is assumed zero or constant, which is unrealistic. We use three approaches:

  • Order Book Simulation — accurate reconstruction of the order book
  • VWAP Model — approximation via candle volume
  • Almgren–Chriss Framework — optimal execution of large orders

Why Realistic Backtesting with Liquidity Is Critical?

Without accurate modeling, you cannot assess the real capacity of a strategy. Capacity is the maximum capital under which the strategy remains effective. Exceeding the threshold, market impact cancels out advantages. In one of our projects: a strategy showed a Sharpe of 2.5 on $100K but dropped to 1.1 on $5M — precisely due to unaccounted slippage. Proper modeling allows early detection of this threshold and strategy adjustment. Our clients save up to $30,000 per month in trading costs after implementing correct modeling.

When to Use Almgren–Chriss Framework?

Almgren–Chriss suits large positions where time risk is significant. If you trade 1% or more of daily volume, this model gives an optimal execution schedule. For small trades, the VWAP model is sufficient. Read more about the model in the Almgren–Chriss framework.

How We Do It: Implementation Details

Order Book Simulation

For precise modeling, we need order book snapshots. We reconstruct how much volume is available at each price level:

import numpy as np
from dataclasses import dataclass

@dataclass
class OrderBookLevel:
    price: float
    quantity: float

@dataclass
class SimulatedOrderBook:
    symbol: str
    timestamp: int
    bids: list[OrderBookLevel]
    asks: list[OrderBookLevel]

    def simulate_market_buy(self, quantity: float) -> tuple[float, float]:
        remaining = quantity
        total_cost = 0.0
        filled = 0.0

        for level in self.asks:
            if remaining <= 0:
                break

            fill_qty = min(remaining, level.quantity)
            total_cost += fill_qty * level.price
            filled += fill_qty
            remaining -= fill_qty

        if filled == 0:
            return 0.0, 0.0

        return total_cost / filled, filled

    def simulate_market_sell(self, quantity: float) -> tuple[float, float]:
        remaining = quantity
        total_proceeds = 0.0
        filled = 0.0

        for level in self.bids:
            if remaining <= 0:
                break

            fill_qty = min(remaining, level.quantity)
            total_proceeds += fill_qty * level.price
            filled += fill_qty
            remaining -= fill_qty

        if filled == 0:
            return 0.0, 0.0

        return total_proceeds / filled, filled

VWAP Market Impact Model

If full order book data is unavailable, we approximate via candle volume:

class VWAPLiquidityModel:
    def __init__(
        self,
        participation_rate: float = 0.05,
        market_impact_coefficient: float = 0.1,
    ):
        self.participation_rate = participation_rate
        self.impact_coeff = market_impact_coefficient

    def estimate_execution_price(
        self,
        side: str,
        order_size_usd: float,
        candle_volume_usd: float,
        candle_close: float,
        daily_volume_usd: float,
    ) -> dict:
        max_fillable_usd = candle_volume_usd * self.participation_rate

        if order_size_usd > max_fillable_usd:
            candles_needed = int(np.ceil(order_size_usd / max_fillable_usd))
            effective_order = max_fillable_usd
            partial_fill = True
        else:
            candles_needed = 1
            effective_order = order_size_usd
            partial_fill = False

        order_fraction = effective_order / daily_volume_usd
        impact_pct = self.impact_coeff * np.sqrt(order_fraction)

        if side == 'BUY':
            execution_price = candle_close * (1 + impact_pct)
        else:
            execution_price = candle_close * (1 - impact_pct)

        return {
            'execution_price': execution_price,
            'filled_usd': effective_order,
            'partial_fill': partial_fill,
            'candles_to_complete': candles_needed,
            'slippage_pct': impact_pct * 100,
        }
VWAP Model Calibration Details

The default market_impact_coefficient is 0.1, but for a specific market, it needs calibration based on historical trades. For highly liquid pairs (BTC/USD), the coefficient may be 0.05; for low liquidity, 0.2. Accurate calibration improves slippage prediction accuracy by 15%.

Almgren–Chriss Framework

A more complex but also more accurate model for evaluating trading costs of large positions:

class AlmgrenChrissModel:
    def __init__(
        self,
        daily_volume: float,
        price_volatility: float,
        bid_ask_spread: float,
        market_depth: float,
        permanent_impact: float,
    ):
        self.V = daily_volume
        self.sigma = price_volatility
        self.epsilon = bid_ask_spread / 2
        self.eta = market_depth
        self.gamma = permanent_impact

    def optimal_schedule(
        self,
        total_size: float,
        time_horizon: int,
        risk_aversion: float = 1e-6,
    ) -> list[float]:
        T = time_horizon
        N = T

        kappa_sq = (risk_aversion * self.sigma**2) / (self.eta / self.V)
        kappa = np.sqrt(max(kappa_sq, 0))

        schedule = []
        for j in range(N):
            t = j / N
            x_j = total_size * np.sinh(kappa * (1 - t)) / np.sinh(kappa)
            if j > 0:
                trade_j = schedule[-1] - x_j if j > 0 else total_size - x_j
                schedule.append(trade_j)

        return schedule

Model Comparison

Parameter VWAP Model Almgren–Chriss
Accuracy Medium High
Data Candle volume only Full order book, volatility, spread
Speed High Medium
Time risk considered No Yes
Application Quick estimates, small orders Large positions, optimal execution

VWAP model runs 5x faster than Almgren–Chriss, but slippage accuracy is 10% lower. This trade-off should be considered when choosing.

How to Scale a Strategy Without Losing Profitability?

Important analytics for investment decisions:

def analyze_capacity(
    strategy_backtest: BacktestResult,
    volume_data: pd.DataFrame,
    participation_rate: float = 0.05,
) -> pd.DataFrame:
    capital_levels = [10_000, 50_000, 100_000, 500_000, 1_000_000, 5_000_000]
    results = []

    for capital in capital_levels:
        scale_factor = capital / strategy_backtest.initial_capital

        adjusted_returns = []
        for trade in strategy_backtest.trades:
            order_size = trade['size_usd'] * scale_factor
            avg_daily_vol = volume_data.loc[trade['date'], 'volume_usd']

            participation = order_size / avg_daily_vol
            extra_slippage = 0.1 * np.sqrt(participation)

            adjusted_pnl = trade['pnl'] * scale_factor - order_size * extra_slippage
            adjusted_returns.append(adjusted_pnl / capital)

        adjusted_sharpe = np.mean(adjusted_returns) / np.std(adjusted_returns) * np.sqrt(252)
        results.append({
            'capital': capital,
            'sharpe': adjusted_sharpe,
            'annual_return_pct': np.mean(adjusted_returns) * 252 * 100,
        })

    return pd.DataFrame(results)

The result is a capital vs return curve. The strategy remains effective up to a certain threshold, after which market impact cancels out the advantage. This threshold is the maximum strategy capacity. Our clients save up to $15,000 per month in trading costs thanks to proper modeling.

Process

Stage Actions
Analysis Examine your strategy, trading data, and goals
Design Select liquidity models, design architecture
Implementation Write code in Python with Foundry/Hardhat for smart contract simulation
Testing Validate on historical data, compare with real trades
Deployment Deliver code, documentation, and train your team

How to Set Up Realistic Backtesting with Liquidity

  1. Collect historical order book data (exchange snapshots).
  2. Choose a liquidity model based on data volume and accuracy.
  3. Calibrate model parameters (market impact coefficient, participation rate).
  4. Integrate the model into the backtesting engine.
  5. Run backtests with order execution simulation.
  6. Analyze capacity and scalability.

What's Included

  • Modular code supporting order book simulation, VWAP, and Almgren–Chriss
  • Integration with your existing backtesting system
  • Documentation and usage examples
  • Team training (2 hours online)
  • 30-day post-delivery support

Timeline and Cost

Development timeline: from 3 to 6 weeks depending on complexity. Cost is estimated individually after analyzing your project. Get a consultation: we will analyze your strategy and suggest the optimal liquidity model. Experience: over 50 successful projects in crypto trading and DeFi. We guarantee code quality and deadline adherence. Request a consultation to learn how we can help your strategy.

Why exchange development requires deep domain expertise

We develop exchanges — not 'chart sites,' but matching engines that process thousands of orders per second without delay, route liquidity between pools, and guarantee that no user gains access to others' funds. Teams that start with the UI and postpone the engine 'for later' end up rewriting everything in six months in 90% of cases.

Order Book vs AMM: where most projects break

Centralized exchanges (CEX) are built around an order book + matching engine. Decentralized exchanges (DEX) either also use an order book (dYdX on StarkEx, Serum/OpenBook on Solana) or an AMM with concentrated liquidity (Uniswap v3/v4, Curve, Balancer). A classic mistake when developing a CEX is implementing the matching engine on top of a relational database with transactions for each match. PostgreSQL handles ~500 RPS without special effort, but at peak loads of 5,000–10,000 orders per second, it turns into a deadlock nightmare. The correct architecture: in-memory order book (Redis Sorted Sets or custom C++/Rust structure), asynchronous writing of matches to PostgreSQL via a queue (Kafka/RabbitMQ), and a separate settlement service that finally updates balances.

For DEX, the most painful problem is sandwich attacks and MEV. A pool with a plain xy=k AMM without slippage protection becomes a target for MEV bots within hours of launch. Uniswap v2 lost hundreds of millions of dollars in user liquidity. Solutions: integration with Flashbots Protect, a commit-reveal scheme for orders, or switching to TWAMM (Time-Weighted AMM) for large trades.

Concentrated liquidity and impermanent loss

Uniswap v3 introduced concentrated liquidity – LPs choose a price range in which to provide liquidity. Capital efficiency increased 4,000x compared to v2 for stable pairs. But implementing this mechanism correctly is non-trivial. The Uniswap v3 liquidity contract uses tick-based accounting: the price space is divided into discrete ticks (tick = log₁.0001(price)), each tick stores accumulated fee growth and liquidity delta. When creating a position, the lower and upper ticks are computed, and the contract recalculates all active positions at each swap. Storage layout is critical here – incorrect variable packing in slots easily adds 40–60% to swap gas cost.

We implemented a Uniswap v3 fork for a client on Polygon with a custom fee tier system. The initial version consumed 180k gas for a swap across 2 ticks. After slot packing of variables in Tick.Info and inlining several internal calls, it dropped to 112k gas. This reduced gas costs by 38% and saved the client substantial costs on fees monthly. The techniques applied are described in the Uniswap v3 Whitepaper and confirmed by our audit experience.

How a matching engine delivers performance

A production-ready matching engine is built according to the following scheme:

  • Order ingestion layer – WebSocket gateway (Go or Rust), accepts orders, validates signature, checks balance via Redis, queues them. Latency at this level must be <1ms.
  • Matching core – single-threaded event loop (eliminates race conditions without mutexes). In memory, we hold two Sorted Sets for each trading instrument: bids and asks. FIFO matching for limit orders, immediate-or-cancel for market orders. Throughput with a proper Rust implementation – 500k–1M matches per second on a single core.
  • Settlement service – reads matches from Kafka, atomically updates balances in PostgreSQL (UPDATE accounts SET balance = balance - $1 WHERE id = $2 AND balance >= $1). Optimistic locking via row versioning.
  • Withdrawal pipeline – separate service with cold/hot wallet architecture. The hot wallet holds 5–10% of total deposits, the rest is cold storage with multi-sig (Gnosis Safe or custom HSM). Automatic withdrawals only from hot wallet, large amounts require manual authorization.
Component Technology Latency / Throughput
Order gateway Go + WebSocket <1ms p99
Matching engine Rust (in-memory) 500k+ orders/sec
Balance store Redis (write-through) <0.5ms
Settlement DB PostgreSQL 14+ ~50k TPS with partitioning
Event streaming Apache Kafka 1M+ events/sec
Blockchain node Geth / Solana validator depends on chain

How our exchange development process ensures reliability

Smart contracts and gas optimization

For EVM-based DEX (Ethereum, Arbitrum, Optimism, Polygon), the entire critical path lives in Solidity. Main contracts: Pool, Factory, Router, PositionManager (for v3-like), and Quoter for off-chain calculations. Typical mistakes we see in audits:

Reentrancy via callback. Uniswap v3 uses flash swap with a callback (uniswapV3SwapCallback). If your router lacks a nonReentrant guard and you don't check msg.sender == pool, the contract gets drained via a nested call. This is not hypothetical – several v3 forks lost funds this way.

Oracle manipulation in AMM. If your contract uses the spot price from the pool for collateral calculation, it is front-runnable. Correct: TWAP over 30+ minutes (Uniswap v3 OracleLib) or an external oracle (Chainlink).

Unbounded loops in liquidity range. If a swap crosses many ticks in a row (price impact 80%+), gas may exceed the block limit. Need MAX_TICKS_CROSSED with partial fill and returning the remainder.

For Solana DEX (Anchor framework, Rust), the architecture is fundamentally different: account-based model, Program Derived Addresses (PDA) instead of storage, Cross-Program Invocations instead of internal calls. Solana's throughput (~3,000–4,000 TPS vs 15–30 on Ethereum mainnet) allows building on-chain order books – exactly what Phoenix DEX does.

Liquidity bootstrapping and aggregator integration

Launching a pool is not enough – you need to ensure liquidity at launch. Practical mechanisms:

  • Liquidity Bootstrapping Pool (LBP) – initial price is high, asset weights dynamically shift, creating selling pressure and even token distribution. Implemented in Balancer v2.
  • Initial Liquidity Offering via Uniswap v3 – adding liquidity in a narrow range around the initial price, then gradually expanding as volume grows. Requires active liquidity management or integration with Arrakis/Gamma.
  • Integration with 1inch, Paraswap, Li.Fi – aggregators bring traffic but require standard compliance: the pool must have correct getAmountsOut, support ERC-20 approval/permit, and not have custom transfer hooks that break the aggregator's routing.

Development process and deliverables

Analytics and design begin with choosing the architectural model: CEX with custodial storage, non-custodial DEX, or hybrid (off-chain order book + on-chain settlement, like dYdX v3). This decision determines everything – regulatory load, tech stack, team.

Development proceeds in layers: first smart contracts with full Foundry coverage (fuzzing, invariant testing), then backend services, then integration layer, and finally frontend. Testing includes fork testing on mainnet via Foundry – we reproduce real liquidity conditions, not synthetic ones.

Audit is mandatory before mainnet deployment. For DEX contracts, minimally one firm with manual review (Trail of Bits, Spearbit, Code4rena contest). For CEX custody, audit of key storage processes. We guarantee all contracts undergo formal verification and fuzzing testing (Echidna, Foundry invariant).

Estimated timelines

Exchange type Timeframe
DEX (AMM, xy=k) 3 to 5 months
DEX with concentrated liquidity (v3-like) 6 to 10 months
CEX (matching engine + custody + trading UI) 8 to 14 months
Integration with existing protocol 4 to 8 weeks

Cost is calculated individually after a technical briefing: chain selection, throughput requirements, custodial model. Our certified engineers with 10+ years of experience will help you choose the optimal architecture and avoid common pitfalls. Contact our team for a detailed proposal.

Pitfalls to avoid at launch

  • Forgetting the price oracle in AMM. Spot price can be manipulated with a flash loan in one transaction. If your lending protocol uses the spot price from its own pool, that's a bug.
  • Hot wallet without limits. A CEX without daily limits on automatic withdrawals is an invitation for attackers. Compromising one key should lose at most 10% of total funds.
  • Absence of circuit breaker. A 40% price drop in 5 minutes should halt automatic liquidations or withdrawals until manual review. Without this, a cascading liquidation spiral destroys all TVL.
  • Incorrect decimal handling. USDC uses 6 decimals, WBTC – 8, most tokens – 18. Mixing without normalization leads to either precision loss or overflow. Solidity has no float; we work with fixed-point using FullMath (mulDiv with overflow protection).

Want to avoid these problems? Get a consultation — we will select the architecture for your project and provide exact timelines. Order exchange development with quality guarantee and ongoing support.