Imagine you are trying to predict the weather six months from now. You could pick a single forecast — "it will probably be mild, around 14°C" — and present that as your best estimate. Or you could run thousands of simulations using different starting conditions and atmospheric variables, and present a probability distribution: there is a 60% chance of mild temperatures, a 25% chance of cold, and a 15% chance of an unusually warm autumn.

The second approach is more honest, more useful, and considerably more complex. It is also, in essence, what Monte Carlo simulation does for retirement planning.

The name comes from the famous casino district in Monaco — a reference to the randomness and probability at the heart of the method. But the technique itself was born not in a casino but in a weapons laboratory in New Mexico, developed by some of the most brilliant mathematicians of the 20th century.

See Monte Carlo simulation applied to your own retirement — 5,000 simulated futures, free.

Create a free account and try odoPT

A brief history — from Los Alamos to your pension

1946 — Los Alamos, New Mexico
Stanisław Ulam and the solitaire problem
Mathematician Stanisław Ulam, recovering from illness, spent time playing solitaire.[1] He found himself wondering: what is the probability of winning a game of Canfield solitaire? Rather than calculating it analytically — which would have required an enormous combinatorial analysis — he realised it would be easier to simply play the game many times and observe how often he won. The insight was profound: for complex problems where analytical solutions are intractable, repeated random sampling can produce accurate probability estimates.
1947 — The Manhattan Project connection
John von Neumann formalises the method
Ulam shared his idea with John von Neumann, the Hungarian-American mathematician who was also working at Los Alamos. Von Neumann immediately saw the application to neutron diffusion problems in nuclear weapons design[2] — problems that were analytically impossible but tractable through repeated random simulation. He helped formalise the method and developed algorithms to generate the random numbers the simulations required. The name "Monte Carlo" was coined by physicist Nicholas Metropolis[3], a reference to Ulam's uncle who frequented the Monaco casinos.
1950s–1970s
Spread to physics, engineering, and finance
As computing power grew, Monte Carlo methods spread from nuclear physics into other domains where complex systems produced unpredictable outcomes: particle physics, fluid dynamics, aerospace engineering, and eventually financial modelling. The method was particularly useful wherever the underlying mathematics involved too many interacting variables for closed-form solutions.
1973
Black-Scholes and financial applications
The Black-Scholes options pricing model demonstrated that probabilistic mathematics could be applied to financial instruments.[4] Monte Carlo methods became a standard tool in quantitative finance for pricing complex derivatives and modelling portfolio risk — particularly for path-dependent instruments where the sequence of prices mattered, not just the final value.
1990s–2000s
Retirement planning applications
As computing became affordable and financial planning software matured, Monte Carlo simulation began appearing in retirement planning tools. The application was natural: retirement outcomes depend on a sequence of market returns over decades, making them inherently path-dependent and unsuitable for simple linear projection. Academic researchers including William Pfau began applying Monte Carlo methods to withdrawal rate research, challenging the deterministic 4% rule with probabilistic alternatives.[5][6]
Today
The gold standard for retirement modelling
Monte Carlo simulation is now considered the most rigorous method available for retirement planning. Tools like odoPT run thousands of simulations using realistic return distributions, capturing the range of outcomes that a retiree might actually experience — rather than the single projected outcome that linear calculators produce.

The mathematics — what is actually happening

At its core, Monte Carlo simulation for retirement planning involves three components: a model of how investment returns are distributed, a mechanism for generating random sequences of returns from that distribution, and a retirement scenario to run those sequences through.

Step 1: Define the return distribution

Historical data tells us that annual investment returns are not fixed — they vary, sometimes dramatically. Equity markets have historically produced a long-run average of roughly 7–8% per year in real terms[7], but individual years range from gains of 40% or more to losses of 50% or more. This variability can be characterised statistically — typically using a normal or log-normal distribution defined by its mean (the average return) and its standard deviation (how spread out the returns are).

Step 2: Generate random return sequences

Using that distribution, the simulation generates a random sequence of annual returns — perhaps 30 or 40 years long, representing a full retirement horizon. This is done using a random number generator seeded with different values for each simulation run, producing a unique sequence each time. Crucially, some sequences will start with a run of good years; others will start with a crash. The distribution of starting conditions reflects the statistical properties of the historical data.

Step 3: Run the retirement scenario

Each random return sequence is applied to the retirement scenario: starting portfolio value, annual withdrawal, asset allocation, tax treatment, State Pension income, and so on. The simulation tracks the portfolio value year by year and determines whether it survived the full planning horizon — or ran out of money, and if so, when.

Step 4: Aggregate the results

After running thousands of simulations — odoPT runs 5,000 — the results are aggregated. The percentage of simulations in which the portfolio survived is the probability of success. The distribution of portfolio values at each point in time produces the fan chart — the spread of possible futures, from pessimistic to optimistic.

The fan chart is the visual signature of Monte Carlo simulation. Each line represents one simulated retirement — one possible future. The width of the band at any point in time represents the uncertainty in the outcome. A narrow band means most futures lead to similar results; a wide band means the outcome is highly dependent on which sequence of returns you happen to experience.

Why Monte Carlo is genuinely good for retirement planning

Strengths
  • Models the full range of possible futures, not just one
  • Captures sequence of returns risk — the order of returns matters
  • Produces a probability of success rather than a false certainty
  • Allows sensitivity analysis — change one input and see the impact
  • Models complex interactions between income sources, tax wrappers, and spending
  • Handles path-dependent strategies like cash buffers and phased withdrawals
  • Honest about uncertainty rather than hiding it
Limitations
  • Results are only as good as the input assumptions
  • Standard implementations assume returns are independently distributed — ignoring crash clustering
  • Historical return distributions may not predict future returns
  • Does not model structural changes — demographic shifts, policy changes
  • Can produce false confidence if limitations are not understood
  • A 90% success rate still means a 10% chance of failure — which is real
  • Standard models do not capture tail risks well

The sequence problem — what linear calculators miss

The most important strength of Monte Carlo simulation for retirement planning is its ability to capture sequence of returns risk. A linear calculator that assumes a fixed 6% return every year will produce the same projected outcome regardless of whether good returns come early or late in retirement. Monte Carlo does not make this assumption — it tests every possible ordering, and the results show clearly that the order matters enormously.

Two retirees with identical average returns over 30 years can have radically different outcomes depending on when the bad years arrive. Only probabilistic simulation can reveal this — and reveal how much it matters for your specific situation.

Honest uncertainty

Perhaps the most important thing Monte Carlo does is refuse to pretend that the future is knowable. A linear calculator says "your money will last until you are 87." Monte Carlo says "in 82% of simulated futures, your money lasts 30 years. In 18%, it does not." The second answer is less comfortable. It is also far more useful, because it tells you there is work to do — and gives you the tools to do it.

Where Monte Carlo can mislead you

Monte Carlo simulation is more honest than linear projection. It is not infallible. There are several significant limitations that users of any Monte Carlo tool — including odoPT — should understand.

The independence assumption

Standard Monte Carlo implementations assume that each year's return is statistically independent of the previous year's return. In reality, markets exhibit clustering[8] — bad years tend to follow bad years, and crises tend to be more severe and prolonged than random sampling from a normal distribution would suggest. A standard Monte Carlo model draws each year's return independently, which can underestimate the probability and severity of extended downturns.

This is the most technically significant limitation of standard Monte Carlo for retirement planning. A model that draws returns independently cannot fully capture a scenario like the 2008 financial crisis, where losses were severe, sustained, and correlated across asset classes. The model will generate some runs that look like 2008 — but fewer of them, and with less severity, than actually occurred.

The input assumption problem

Monte Carlo simulation is only as good as its inputs. The mean return and standard deviation used to define the return distribution are typically derived from historical data. But past returns do not guarantee future returns. A model calibrated on 20th century US equity returns — which were exceptionally strong by global historical standards — may overstate the probability of success for future retirees facing different return environments.[10]

This is particularly relevant for UK investors. UK equity market returns have historically been lower than US returns, and gilt yields have followed a different path. Using US-calibrated assumptions for a UK retirement plan will produce systematically optimistic results.

The false precision problem

Seeing "82% probability of success" feels precise. It is not. The actual probability depends entirely on the assumptions built into the model — return distribution, inflation, longevity, tax rates. Change any of those assumptions and the number changes. An 82% result from a generous set of assumptions might correspond to a 68% result under more conservative ones. The number should be treated as an indicator, not a measurement.

The 10% tail problem

A 90% success rate is often treated as "safe." But a 10% probability of failure is not negligible — it means that in one in ten simulated futures, the portfolio runs out of money. For a retiree facing a 30-year retirement, a 10% chance of financial ruin in old age is a risk that deserves serious attention, not a footnote. Good Monte Carlo tools show the downside scenarios explicitly, not just the headline success rate.

How odoPT is designed to address these limitations

odoPT is built with an awareness of Monte Carlo's limitations. Here is how the tool addresses each of the major concerns:

LimitationHow odoPT addresses it
Independence assumption underestimates crises Crisis scenario modelling. Rather than relying solely on independently-drawn random returns, odoPT allows you to specify structured crisis scenarios — based on historical events like the 2008 GFC, the dot-com crash, or the 2020 COVID shock. These scenarios inject a defined crash profile into the simulation, capturing the severity and duration that standard random sampling may underestimate. You can stack multiple crises and test combined shock scenarios.
Historical return assumptions may be optimistic UK-calibrated return assumptions. odoPT uses return distributions calibrated to UK market history across multiple asset classes, not US data. The 32 available asset classes each carry historically-grounded return and volatility assumptions specific to the UK investment environment.
Single withdrawal rate ignores flexibility Multiple withdrawal strategies. Rather than assuming a fixed withdrawal amount, odoPT supports seven distinct withdrawal strategies including cash buffer band, floor-plus-upside, and expense-matched approaches. These strategies allow the simulation to model realistic spending flexibility — drawing less in bad years, more in good ones. See the withdrawal strategies guide for a full comparison.
False precision in the success rate Fan chart and percentile display. odoPT shows the full probability envelope — P5, P10, P25, P50, P75, P90, P95 — not just the headline success rate. The fan chart makes the uncertainty visible. The P10 outcome (the 10th percentile — one of the worst simulated futures) is displayed alongside the median, ensuring the downside is never hidden.
Tail risks underrepresented 5,000 simulations + crisis stacking. Running 5,000 simulations rather than a few hundred increases the representation of tail outcomes. Combined with explicit crisis scenario modelling, this gives a more complete picture of the severe downside scenarios that standard random sampling may undercount.
Does not model tax and vehicle interactions Full tax overlay. odoPT models the tax treatment of pension drawdown, ISA withdrawals, and GIA disposals — including income tax, capital gains tax, and the interaction with the personal allowance and State Pension income. The simulation tracks after-tax income rather than gross withdrawals, giving a more accurate picture of what you actually have to spend.

The crisis sequence approach — why it matters

The single most important enhancement odoPT makes to standard Monte Carlo is the ability to model structured crisis scenarios rather than relying entirely on randomly generated return sequences.

Here is why this matters. Standard Monte Carlo generates crashes — some simulated sequences will include severe years, and the distribution will produce some very bad outcomes. But the crashes generated by random sampling tend to be short, sharp, and followed by rapid recovery. Real market crises are often different: they are prolonged, they correlate across asset classes, and they arrive with particular timing patterns that random sampling does not fully replicate.

The 2008 global financial crisis took approximately 55 months to reach its peak drawdown and roughly 50 months to recover.[9] The dot-com crash took longer. These are not the short, sharp shocks that a normal distribution generates most naturally — they are sustained, structural events with specific profiles.

By modelling the 2008 GFC as a structured scenario — with its specific run-up, peak drawdown of 57%, duration to trough, and recovery profile — odoPT allows you to ask the direct question: "If I retire and the next five years look like 2008, what happens to my plan?" That is a more useful question than "what happens in a randomly bad sequence?"

You can also stack crises — modelling a scenario where a crash is followed by a period of elevated inflation, or where two crises arrive within a few years of each other. These combined shock scenarios test the robustness of a retirement plan in a way that standard Monte Carlo cannot replicate.

Withdrawal strategies — the other half of the equation

Monte Carlo simulation models what markets do to your portfolio. Withdrawal strategies model what you do in response — and the interaction between the two determines the outcome.

A fixed withdrawal strategy — taking the same amount every year regardless of market conditions — is the simplest approach and the one most basic retirement calculators assume. It is also the most rigid, and the most vulnerable to bad sequences. If markets fall 30% and you continue withdrawing the same fixed amount, you are selling units at the worst possible time with no flexibility.

More sophisticated strategies — including the cash buffer band approach, floor-plus-upside, and expense-matched withdrawal — allow spending to respond dynamically to portfolio performance. In down years, spending falls slightly. In good years, it rises. This flexibility can dramatically improve retirement sustainability without requiring a lower average withdrawal rate.

odoPT models seven distinct withdrawal strategies and shows you how each one performs across the full range of simulated futures. The differences between strategies can be striking — and the right choice depends on your specific income needs, flexibility, and risk tolerance.

For a full comparison of the withdrawal strategies available in odoPT and guidance on which is most appropriate for different situations, see the withdrawal strategies guide.

What Monte Carlo cannot do — and what that means for you

Even the most sophisticated Monte Carlo implementation cannot predict the future. It can tell you the probability of different outcomes under a defined set of assumptions — it cannot tell you which future you will actually live in.

This is not a failure of the method. It is a fundamental property of complex systems. Markets are influenced by events that are not in the historical data — pandemics, wars, technological disruptions, policy changes — and no model can fully anticipate them.

What Monte Carlo can do is give you a structured way to think about uncertainty, to stress-test your plan against a range of plausible scenarios, and to make decisions that are robust across those scenarios rather than optimised for a single projected future.

The value of Monte Carlo simulation is not that it tells you what will happen. It is that it tells you how your plan performs across many possible versions of what might happen — and gives you the information you need to improve it.

A plan that survives 85% of simulated futures is more robust than one that survives 65% — even if neither can guarantee survival. The goal is not certainty, which is unavailable. The goal is a plan that is resilient enough across enough scenarios that you can retire with justified confidence rather than false certainty or unnecessary anxiety.

Monte Carlo simulation, used honestly and with an understanding of its limitations, is the best tool currently available for building that kind of plan. It is imperfect. It is still far better than the alternative.

Run Monte Carlo simulation on your own retirement — 5,000 futures, UK-calibrated, free.

Create a free account and try odoPT
Academic & Primary Sources
  1. [1] Ulam, S. M. (1991). Adventures of a Mathematician. University of California Press. Ulam's own account of the solitaire problem and the development of the Monte Carlo method at Los Alamos. Book
  2. [2] Metropolis, N., & Ulam, S. (1949). The Monte Carlo Method. Journal of the American Statistical Association, 44(247), 335–341. The original published description of the method. Academic
    doi:10.2307/2280232
  3. [3] Metropolis, N. (1987). The Beginning of the Monte Carlo Method. Los Alamos Science, Special Issue, 125–130. Metropolis's first-hand account of naming the method and its development. Academic
  4. [4] Black, F., & Scholes, M. (1973). The Pricing of Options and Corporate Liabilities. Journal of Political Economy, 81(3), 637–654. Academic
    doi:10.1086/260062
  5. [5] Bengen, W. P. (1994). Determining Withdrawal Rates Using Historical Data. Journal of Financial Planning, 7(4), 171–180. The original paper establishing the 4% rule based on US historical market data. Academic
  6. [6] Pfau, W. D. (2011). Safe Savings Rates: A New Approach to Retirement Planning over the Lifecycle. Journal of Financial Planning, 24(5), 42–50. Pfau's probabilistic challenge to the deterministic 4% rule. Academic
  7. [7] Dimson, E., Marsh, P., & Staunton, M. (2002). Triumph of the Optimists: 101 Years of Global Investment Returns. Princeton University Press. The definitive long-run study of global equity returns including UK market history. Book
  8. [8] Mandelbrot, B., & Hudson, R. L. (2004). The (Mis)Behaviour of Markets: A Fractal View of Risk, Ruin and Reward. Profile Books. Foundational work on fat tails, volatility clustering, and why standard normal distributions underestimate financial risk. Book
  9. [9] National Bureau of Economic Research. (2010). US Business Cycle Expansions and Contractions. NBER Business Cycle Dating Committee. Used for GFC timeline (peak December 2007, trough June 2009; equity recovery to prior peak approximately 2013). Data
    nber.org
  10. [10] Clare, A., Seaton, J., Smith, P. N., & Thomas, S. (2017). Sequence of Returns Risk, Pension Drawdown Risk and Sustainable Withdrawal Rates. Journal of Retirement, 5(1), 37–53. UK-specific analysis finding sustainable withdrawal rates materially lower than US-derived 4% rule. Academic
  11. [11] Guyton, J. T., & Klinger, W. J. (2006). Decision Rules and Maximum Initial Withdrawal Rates. Journal of Financial Planning, 19(3), 48–58. Research demonstrating that dynamic withdrawal rules — reducing spending in down markets — substantially improve portfolio sustainability. Academic
TermMeaning
Monte Carlo simulationA computational method that uses repeated random sampling to model the probability of different outcomes in a system with inherent uncertainty. In retirement planning, it runs thousands of simulations using different sequences of market returns to produce a probability distribution of outcomes.
Return distributionA statistical description of how investment returns are spread — characterised by a mean (average return) and standard deviation (variability). Used to generate random return sequences in Monte Carlo simulation.
Standard deviationA measure of how spread out a set of values is around the mean. A higher standard deviation means returns are more variable — higher highs and lower lows.
Fan chartA visual representation of Monte Carlo output showing the spread of simulated portfolio values over time. The shaded band shows where the bulk of outcomes cluster; the edges show the optimistic and pessimistic extremes.
Probability of successThe percentage of simulated scenarios in which the portfolio does not run out of money before the end of the planning horizon.
Path-dependentA property of systems where the sequence of events matters, not just the final state. Retirement planning is path-dependent because the order of investment returns — not just their average — determines whether the portfolio survives.
Crisis scenarioIn odoPT, a structured shock event — based on historical market crises — applied to a simulation path. Captures the severity and duration of real-world crashes more accurately than random sampling alone.
Tail riskThe risk of extreme outcomes — events in the far ends of a probability distribution. Standard Monte Carlo implementations can underestimate tail risk because extreme events cluster in ways that random sampling does not fully replicate.