An open firm-level model for UK VAT threshold reform
Abstract.
I present an open, reproducible firm-level microsimulation for costing UK VAT registration-threshold reforms, released with its synthetic-data generator. The threshold is a notch: once turnover crosses £85,000, VAT falls due on a firm’s whole base. I build synthetic firm records calibrated to 2023–24 HMRC and OBR statistics, and cost level, shape (a taper) and rate (a reduced-rate band) reforms statically on a common £184.8bn base: raising the threshold to £100,000 costs £753m, a taper £520m, and reduced-rate bands £484m and £242m. I characterise the notch by its dominated region—the £21,250 band above the threshold in which no firm profits from locating—and how each reform changes it. Conditional on an assumed elasticity, intensive-margin responses barely move these costings. The near-threshold shape is target-inherited: a placebo regenerated without the OBR targets returns zero excess mass, so band-calibrated data support no bunching inference.
Keywords: value-added tax, microsimulation, tax notch, policy
The United Kingdom’s Value Added Tax (VAT) applies above a registration threshold that is the highest in the OECD alongside Switzerland (Seely 2024; HM Revenue and Customs 2024). The threshold is a notch rather than a kink: once a firm’s annual taxable turnover crosses it, VAT falls due on the firm’s whole base, not only on the part above the threshold, so registration produces a discrete jump in liability rather than a marginal one. The threshold falls in a dense part of the firm distribution. Figure 1(a) shows that VAT- and PAYE-registered UK businesses are concentrated at the lower end of the turnover distribution, and Figure 1(b) shows that the VAT-registered population per £1,000 of turnover is densest at and below the threshold, so the notch applies to a large number of firms.
VAT- and/or PAYE-registered businesses (ONS, 2023–24)
VAT-registered firms (HMRC, 2023–24)
Firms respond to the notch by locating below the threshold. Tabulations of HMRC turnover data show the number of firms rising towards £84,000 and then dropping sharply at the £85,000 threshold, leaving excess mass just below it (Office for Budget Responsibility 2023). Liu et al. (2021) document this bunching in UK administrative firm data over 2004–05 to 2014–15, when the threshold rose from £58,000 to £81,000, and Liu et al. (2024) find that firms slow their turnover growth as they approach the threshold. The notch therefore distorts the size distribution and the growth decisions of firms near the margin, which places the threshold under recurrent pressure for reform.
Reform proposals change the threshold in one of three ways: its level (raising or lowering it), its shape (replacing the notch with a schedule that phases the rate in over a band), or its rate (a reduced rate for firms in a band above the threshold). The three are not equivalent. Changing the level moves the notch along the turnover axis; changing the shape or the rate alters the size of the discrete jump itself. Costing the three on a common basis is the practical difficulty this paper addresses. A change in the level can be costed by mechanical accounting, but a change in the shape of the schedule cannot be evaluated with the reduced-form bunching methods used to study the threshold, because a continuous schedule leaves no single point of excess mass to exploit. The firm-level data needed to study turnover near the threshold is, in addition, not available in open form.
In this paper I build an open, reproducible firm-level microsimulation for costing UK VAT threshold reforms on a common static basis. I construct synthetic firm microdata calibrated to published official aggregates—HMRC registered-firm counts by turnover band and by sector, net VAT liability by turnover band, and ONS business counts and employment bands—for the 2023–24 tax year, when the threshold was £85,000, and separately for 2024–25, when it rose to £90,000. Each firm’s net VAT liability is the standard rate applied to its value added. The synthetic data can be released openly and regenerated under any counterfactual schedule.
On these data I do four things. First, I cost level reforms and schedule reforms—a graduated taper and a reduced-rate band—statically, by applying each counterfactual schedule mechanically to the firm distribution. Second, I characterise the notch’s distortion through its dominated region—the exact, elasticity-free range of turnover above the threshold in which no firm profits from locating—and how each reform changes it. Third, I price the level and rate reforms behaviourally: each firm re-optimises its turnover within its schedule region, conditional on an assumed turnover elasticity \(e\) that the data do not identify, reported as an \(e\)-sensitivity range that nests the static costing exactly in its \(e\to0\) limit. Fourth, I show that the synthetic data’s near-threshold structure is target-inherited and cannot support bunching inference: the £85,000 profile is calibrated to the OBR’s published £1,000-band data and the estimator reads exactly that back, a placebo without the fine targets returns zero, and the uncalibrated £90,000 vintage shows a spurious band-edge signal.
Four findings follow. First, the static pipeline prices all the reforms on a common basis, and its anchor costing tracks the Government’s published number once the deregistration convention is stated. My static costing of the April 2024 increase from £85,000 to £90,000 depends on what released firms are assumed to do: if every firm in the band deregisters, the 2025–26 cost is £323m, overshooting HMRC’s published £185m (HM Treasury 2024) in line with the official costing’s phased registration response; if the roughly 43% voluntary-registration share that Liu et al. (2021) document is retained, it is £184m, in line with the published figure—agreement finer than the model’s liability-level calibration error supports reading as precision (Section 4). Both conventions reproduce the official profile’s narrowing and its 2028–29 sign change. Forward moves from the £90,000 base range from a gain of about £1,040m at £70,000 to a cost of about £1,450m at £120,000, against a base of about £200.9bn in 2025–26.
Second, the shape and rate reforms act on the notch’s distortion in a way that a change in level does not. The notch creates a dominated region of width \(a = T^{*}\tau/(1-\tau) = \unicode{xA3} 21{,}250\) above the threshold (£85,000–£106,250), a range of turnover in which no firm can locate without being strictly worse off than at the threshold. This region is an exact property of the schedule: under the value-added formulation of Section 7 its width is the same for every firm with positive value added, whatever its input share, and it requires no behavioural calibration. It also spans a populous part of the firm distribution—about 148,000 weighted firms in the calibrated data—so its width is an economically meaningful distortion index. Raising the threshold relocates this region without removing it. A reduced rate does not remove it: the jump at the threshold falls (from £17,000 to £12,750 at a 15% band), but the band reverts to the standard rate at its £105,000 top, introducing a second notch and a second dominated region, so the total dominated width is essentially unchanged—£21,562 at 15% and £22,569 at 10%, both marginally above the £21,250 baseline. Only the graduated taper, being continuous, removes the dominated region. On the common £85,000 data-year base, raising the threshold to £100,000 costs £753m and removes none of the dominated region; the taper removes all of it for £520m; and the 10% and 15% reduced rates cost £484m and £242m while splitting the dominated band in two. On the dominated-region metric the level move is the most expensive option per unit of distortion removed—a comparison that indexes the misallocation distortion only and omits the compliance-cost saving a higher threshold delivers (Section 8).
Third, conditional on the assumed elasticity, intensive-margin responses barely move these costings. A pure level rise has exactly zero intensive-margin revenue offset at every \(e\): released firms leave the VAT base, so their expansion is untaxed, and no other firm’s effective rate changes—a result the simulator verifies to £0.1m. For the reduced-rate bands the offset is second order: at \(e=0.17\) the 10% band’s cost falls from £484m to £475m and the 15% band’s from £242m to £235m, under 6% even at \(e=0.32\). I do not report a behavioural cost for the graduated taper: its rate varies continuously with turnover, so its marginal-rate channel is outside the flat-rate machinery, and only its static cost and dominated-region removal are reported. The sweep values \(e \in \{0.05, 0.17, 0.32\}\) are assumptions: the lower end has an external anchor (Kleven and Waseem 2013), and the published UK VAT-base estimates of Liu et al. (2021) fall between it and the midpoint. The consequence of these results runs opposite to the usual worry: the static costings are robust to the intensive margin, and the behavioural uncertainty that matters—registration and location choice at the notch—is exactly the margin the synthetic data cannot identify.
Fourth, the open data carry a transparency lesson about bunching: every sub-band feature of band-calibrated synthetic data is exactly what its targets put there, and this paper demonstrates the point in three directions. The 2023–24 population’s near-threshold density is calibrated to the OBR’s published £1,000-band profile (Office for Budget Responsibility 2023), so it reproduces the administratively observed bunching shape—and the estimator reads that shape back as excess mass (\(E=7{,}933\)), which is target inheritance, not behaviour. A placebo that regenerates the population without the fine targets returns \(E=0\): remove the targets and the “bunching” vanishes. And the 2024–25 vintage, for which no fine bands are published, shows what the trap looks like from the other side: the coarse £90,000 band edge alone produces an excess-bunching statistic within 10% of the administrative estimate of Liu et al. (2021)—a coincidence of band geometry, four-fold unstable across estimator specifications, with no behaviour behind it. An earlier draft of this paper reported “reproducing” the administrative step when a liability mis-scaling in the generator manufactured it; the correction removed it, and the OBR targets now supply the shape from published data instead. The lesson is symmetric and general: synthetic bunching is not evidence of behaviour, its absence is not evidence of no behaviour, and any sub-band feature of calibrated data should be traced to its targets before it is read as economics. The behavioural evidence on UK VAT bunching remains Liu et al. (2021)’s administrative estimate.
The paper draws on four strands of work. The first is the methodology of bunching and notches (Chetty et al. 2011; Saez 2010; Kleven and Waseem 2013; Kleven 2016; Best et al. 2015), which provides the analytic notion of a dominated region; I use the dominated region to characterise a distortion that leaves no single point of excess mass to exploit.
The second is the empirical literature on VAT registration thresholds (Liu et al. 2021, 2024; Onji 2009; Harju et al. 2019; Asatryan and Peichl 2017; Nandi and Warwick 2020; Bashir et al. 2024). Liu et al. (2021) is closest: it estimates bunching at the UK notch on administrative data over 2004–2014 and documents voluntary registration; I take its estimates as the established behavioural facts and cost shape and rate reforms statically on a common basis, which that literature does not do. The international evidence (Harju et al. 2019; Onji 2009; Asatryan and Peichl 2017; Nandi and Warwick 2020) notes that bunching at a VAT threshold can reflect compliance costs, firm-splitting, or reporting rather than real turnover responses.
The third strand questions whether bunching identifies a structural elasticity at all: Blomquist et al. (2021) show that the elasticity is not non-parametrically identified from bunching at a kink or a notch without strong functional-form assumptions, and Bertanha et al. (2023) develop the notch case. This literature motivates my refusal to read an elasticity off the synthetic data: even setting aside the placebo evidence that any synthetic step is mechanical, the mapping from excess mass to an elasticity requires assumptions I am unwilling to impose. The fourth strand is the structural literature on size-based thresholds (Garicano et al. 2016; Gourio and Roys 2014), which identifies policy-invariant primitives at the French 50-employee threshold and computes welfare from them; the present paper makes no such structural claim.
This paper makes three contributions. First, it provides an open, reproducible firm-level microsimulation for costing UK VAT threshold reforms, where the existing firm-level evidence rests on restricted administrative data; the released pipeline regenerates every number in the paper from public inputs. Second, it characterises the notch’s distortion through its exact dominated region and costs two schedule reforms—a graduated taper and a reduced-rate band—on the same static basis as level reforms. Third, it demonstrates, by construction and by correction, that band-calibrated synthetic data carry no sub-band behavioural information beyond their targets: the paper documents how its own earlier build manufactured a plausible-looking density step, how the correction removed it, how calibrating to the OBR’s published fine bands restores the administrative shape as an explicit target, and how the same band-edge mechanism manufactures a spurious signal where no fine targets exist. Identifying the turnover elasticity from administrative micro-data, a welfare accounting, revenue-neutral designs, and an explicit model of voluntary registration are left for future work.
The remainder of the paper proceeds as follows. Section 2 describes the institutional background of the UK VAT threshold. Section 3 sets out the synthetic firm microdata and its calibration to official aggregates. Section 4 sets out the static costing of threshold changes and the schedule reforms. Section 5 shows that the synthetic data contain no usable bunching signal and validates the estimator’s power on injected signals. Section 6 develops the notch and its exact dominated region, and shows how each reform changes it. Section 7 prices the reforms behaviourally, conditional on an assumed elasticity, and derives the level-rise invariance. Section 8 concludes.
The United Kingdom levies a value-added tax at a standard rate of 20% on most goods and services. A business must register once its taxable turnover over any rolling twelve-month period exceeds the registration threshold, after which it charges VAT on its taxable sales and reclaims VAT paid on inputs. Because registration makes VAT due on the firm’s entire taxable turnover rather than only the amount above the cutoff, the threshold is a notch—a discrete jump in liability—and not a kink. The cost of crossing is largest for consumer-facing firms with limited pass-through and small reclaimable input shares, for whom the incentive is to restrain turnover and locate just below the threshold (Liu et al. 2021; Kleven and Waseem 2013).
Real VAT is a tax on value added, not on turnover: a registered firm charges output VAT on its sales but reclaims input VAT on its purchases. The model reflects this: each synthetic firm’s net liability is the standard rate applied to its value added (Section 3), and under the value-added formulation—in which unregistered firms bear no unreclaimed input VAT—the notch’s dominated region is invariant to the firm’s input share (Section 6). Two features of the real system remain outside the model’s scope. A large share of below-threshold firms register voluntarily—roughly \(43\%\) of eligible below-threshold firms (Liu et al. 2021), driven by input-VAT reclaim and business-to-business sales—and for them the threshold carries no bunching incentive; voluntary registration enters the model only as a registration-rate parameter and a retention sensitivity, not as a choice. And net-repayment traders (zero-rated exporters and similar, HMRC’s negative-liability column) are not represented, because the public aggregates do not resolve rated structure by turnover band.
The registration threshold was held at £85,000 from 1 April 2017 to 31 March 2024—a seven-year nominal freeze, the longest in the history of the tax—before being raised to £90,000 on 1 April 2024, with the deregistration threshold rising from £83,000 to £88,000 (HM Revenue and Customs 2024; Seely 2024). Because nominal turnover grows with inflation and real activity, freezing the threshold lowers it in real terms and draws an increasing number of firms towards the registration margin; this fiscal-drag mechanism was a central motivation for the April 2024 increase. My data year is 2023–24, so the prevailing statutory threshold throughout my analysis is £85,000. Figure 2 shows the pattern in the turnover distribution: the number of firms rises towards £84,000 and falls sharply above £85,000, and the projected 2025–26 distribution is more sharply bunched than the 2019–20 outturn as fiscal drag draws more firms towards the frozen threshold.
The UK threshold is the highest in the OECD alongside Switzerland (HM Revenue and Customs 2024), and there is debate over whether to raise it—to reduce the compliance burden on small firms—or lower it, to broaden the base and reduce the growth distortion at the margin. Three families of reform feature, and I evaluate all three: a change in the threshold’s level (the £85,000 to £90,000 increase of April 2024 and further moves around it); a change in its shape (a graduated taper that phases the effective rate in over a band, replacing the notch with a continuous profile); and a change in its rate (a reduced rate for small firms in a band above the threshold, with the standard 20% rate retained above). The taper and reduced-rate designs follow proposals in the UK policy debate: Office of Tax Simplification (2017) canvassed a financial incentive or reduced effective rate to smooth the cliff-edge at the threshold, and HM Treasury’s call for evidence on the registration threshold examined smoothing mechanisms in the same family (HM Treasury 2018). The European Union ran a closely related “graduated relief” scheme for small enterprises and abolished it from January 2025 on complexity grounds, in favour of a simple exemption with a transition tolerance (Council of the European Union 2020)—relevant precedent for the administrability of taper-style designs. A banded reduced rate introduces a smaller secondary notch at the band top, which the continuous taper avoids. I cost all three by the same machinery in Section 4.
Studying firm responses at the threshold needs firm-level turnover near the cutoff, which the official statistics publish only in coarse turnover bands. I therefore build a synthetic firm population calibrated to those aggregates: it reproduces the official totals while resolving the individual firm that the authorities do not release, opening the analysis to anyone. The same calibrate-to-published-aggregates approach underlies open household-level tax-benefit microsimulation (PolicyEngine 2025); here I apply it at the firm level. The generator takes the registration threshold as a parameter and is documented in Appendix 9.1; I build one population for 2023–24 (threshold £85,000), the paper’s primary data year, and one for 2024–25 (threshold £90,000), used for the forward sweep of Section 4.
The archived CSV inputs in this repository remain the reproduction source for the paper’s reported numbers; a pinned migration snapshot also represents the same numeric target surface through PolicyEngine Ledger and the experimental Populace firm generator, and Appendix 9.1 reports the parity check between that shared source-of-truth path and the paper inputs.
Index firms by \(i\). For each sector \(s\) and ONS turnover band \(b=[\underline{b},\overline{b}]\), the UK Business: Activity, Size and Location table gives a business count \(N_{s,b}\), and I draw \(N_{s,b}\) firms with within-band turnover \[y_i \;=\; \underline{b} + (\overline{b}-\underline{b})\,u_i + \varepsilon_i, \qquad u_i \sim \mathcal{U}(0,1), \;\; \varepsilon_i \sim \mathcal{N}(0,\sigma_b^2),\] truncated below at zero, where \(\varepsilon_i\) smooths the within-band distribution. The ONS register covers businesses registered for VAT and/or PAYE—about 2.7 million—so the synthetic population represents that registered-business frame, not the full business population of roughly 5.6 million in the DBT business population estimates; the difference is concentrated below the threshold and matters for threshold cuts, which I flag where relevant. Each firm is given input expenditure \(x_i = \rho_i\,y_i\), with the input–output ratio \(\rho_i\) drawn from a rescaled Beta distribution with sector-specific shifts, clamped to \(\rho_i \in [0.1, 0.95]\) (mean about 0.6, so value added averages about 40% of turnover, in line with the UK non-financial business economy). Its net VAT liability is the standard rate applied to its value added, \[v_i \;=\; \tau\,(y_i - x_i), \qquad \tau = 0.20,\] so every firm’s net rate \(v_i/y_i = \tau(1-\rho_i)\) lies in \([1\%, 18\%]\), strictly positive and below the statutory cap of 20%. The model therefore abstracts from net-repayment traders (firms with zero-rated outputs whose input reclaim exceeds output VAT—HMRC’s negative-liability column); representing them requires rated-structure detail (zero, reduced, and exempt output shares) that the public aggregates do not resolve by band. Employment is assigned from the ONS employment-band shares.1
The base population reproduces the ONS structure but not the HMRC VAT-registered totals, so I re-weight it. Each firm receives a strictly positive weight \(w_i = e^{\theta_i}\), and for a vector of official targets \(T=(T_k)_{k=1}^{K}\) the weighted synthetic counterpart of target \(k\) is \[\widehat{T}_k(\theta) \;=\; \sum_i a_{ki}\, w_i ,\] where \(a_{ki}\) is firm \(i\)’s contribution to target \(k\): an indicator of band or sector membership for a count target, and the firm’s net liability \(v_i\) for a VAT-liability target. I choose \(\theta\) to minimise a symmetric relative-error loss with per-target importance weights \(\lambda_k\), \[L(\theta) \;=\; \frac{1}{K}\sum_{k} \lambda_k \min\!\Big\{ \big(\widehat{T}_k/T_k - 1\big)^2,\; \big(T_k/\widehat{T}_k - 1\big)^2 \Big\} \;+\; \frac{\alpha}{N}\sum_i \lvert \theta_i \rvert ,\] the final term a mean absolute penalty on the log-weights. The targets are the HMRC VAT-registered counts by turnover band and by sector, the HMRC net VAT liability by turnover band, the ONS employment-band totals, and—for the 2023–24 vintage—the OBR’s published near-threshold profile: the March 2023 Economic and Fiscal Outlook Chart C data give HMRC counts of businesses by £1,000 turnover band over £65,000–£90,000 (Office for Budget Responsibility 2023), which I interpolate to the 2023–24 data year between the 2019–20 outturn and the OBR’s 2025–26 frozen-threshold projection (the chart’s final row, the £90,000 bin, lies beyond the fine window and is not used). All twenty-five bins enter as shape targets, applied side-consistently: each side of the threshold takes the OBR profile’s within-side shape, normalised over its own window (£65,000–£85,000 below; £85,000–£90,000 above) and scaled to the synthetic population’s own mass there. The chart counts HMRC traders—a different unit and, below the threshold, a substantially narrower universe than the ONS business frame—so its levels are not imported on either side. An earlier build that mixed conventions (direct counts above the threshold, frame-scaled shape below) placed more mass just above the threshold than just below it, inverting the administrative profile’s cross-threshold ordering; under side-consistent scaling the cross-threshold step comes from the frame’s own band structure and the OBR data supply only the within-side geometry—the rise into the threshold and the decline beyond it. These fine targets pin the near-threshold density to the OBR profile’s within-side geometry at frame-implied levels—the region that coarse bands alone leave to the weight optimiser (Section 5); the resulting cross-threshold step (\(-25\%\)) is about half the OBR chart’s own (\(-60\%\)), because the chart’s below-threshold universe is narrower than the frame. No fine bands are published for the £90,000 threshold era, so the 2024–25 vintage keeps coarse bands only. Turnover-band and near-threshold targets carry the largest \(\lambda_k\), followed by VAT liability by band. Two liability dimensions are reported as informational diagnostics but excluded from the optimizer. VAT liability by sector is excluded because the generator draws a common input-share distribution with sector shifts rather than calibrating the input/output tax structure by sector. The below-threshold liability band (£1–threshold, £1.46bn in 2023–24) is excluded because that total is remitted by voluntary registrants—input-reclaim traders averaging about £2,150 of net remittance—whom the model’s standard-rate-on-value-added liability does not represent (Section 2); imposing it on the whole below-threshold population would force the optimiser to suppress weights in the £50,000–£85,000 region, distorting the near-threshold shape the fine targets exist to pin. The problem is solved by gradient descent, and the symmetric loss keeps targets of very different magnitudes on a common scale.
A firm is registered with certainty when its data-year turnover exceeds the threshold \(T^{*}\). This is an approximation to the statutory test, which is a rolling 12-month test on taxable turnover with a roughly two-month registration lag, a forward-looking 30-day test, and an exception for temporary breaches (Section 2); the model also registers firms whose turnover is largely exempt, whom the law would not require to register. For a single-year cross-section calibrated to registered-population totals these timing and composition wedges are second order for the aggregates, but they are worth naming: the HMRC band tables classify traders by total declared outputs (VAT-return Box 6), and the ONS turnover measure is likewise total turnover, so the running variable throughout is total rather than taxable turnover. Below the threshold a firm registers voluntarily with probability \(r = M / n_{<}\), where \(M\) is the HMRC count of registered firms in the £1–threshold band and \(n_{<}\) is the unweighted number of synthetic below-threshold firm rows. This sets the expected row-level voluntary-registration rate; I do not treat below-threshold voluntary registration as a separate weighted calibration target. The resulting dataset comprises approximately \(2.94\) million firm records per vintage, weighted to represent roughly \(2.5\) million registered businesses. Checked against the official sources along every calibrated dimension, the vintages attain overall calibration accuracies of 89.2% (2023–24) and 92.8% (2024–25), with dimension scores ranging from 77.9% (2023–24 employment bands) to 94.2%; Appendix 9.1 reports the full table. The headline is the mean of five calibrated-dimension scores, each clipped at zero as \(\max\{0,1-\lvert \widehat{T}-T\rvert/\lvert T\rvert\}\); it excludes the VAT-liability-by-sector diagnostic described above.
My primary data year is 2023–24, when the statutory threshold was £85,000. The HMRC band targets step down across the threshold—678,350 firms in the £1–£85,000 band against 305,320 in £85,000–£150,000, or roughly \(8{,}000\) against \(4{,}700\) firms per £1,000 of band width—and the calibrated population necessarily inherits that band-average step. The local shape within a few thousand pounds of the threshold is pinned by the OBR near-threshold targets described above: the calibrated density rises into the threshold (about \(13{,}800\) firms per £1,000 in the £84,000–£85,000 bin), steps down across it (about \(10{,}300\) just above), and declines through the band—the administratively observed bunching geometry (Office for Budget Responsibility 2023; Liu et al. 2021), present in the file because it is targeted, not because any synthetic firm chooses it. In an earlier build without these targets the local shape was an unpinned optimiser equilibrium (and, under the pre-correction liability scale, a spurious bunching look-alike); the central data caveat of the paper stands in all three configurations: no sub-band feature of band-calibrated synthetic data carries behavioural information beyond what its targets put there (Section 5 quantifies this with placebo and recovery tests). Counterfactual policies—including the £85,000 to £90,000 reform and the forward threshold sweep of Section 4—are evaluated on these firms by re-applying the registration rule to the microdata for the static costings, holding the empirical turnover distribution fixed. A behavioural layer that re-prices firms under counterfactual schedules, conditional on an assumed turnover elasticity, is developed in Section 7.
A static costing measures the first-order revenue effect of a reform by re-applying the tax rule to a fixed firm population, with no behavioural response. Let firm \(i\) have turnover \(y_i\), calibrated net VAT liability \(v_i = \tau(y_i - x_i)\) (Section 3), and weight \(w_i\). A firm remits VAT only once registered, which the model equates with turnover above the threshold \(T\), so total VAT revenue at threshold \(T\) is the weighted sum of liabilities over registered firms, \[R(T) \;=\; \sum_i w_i\, v_i \,\mathbf{1}\{y_i \ge T\}.\] Two scope conventions follow from this formula and are worth stating. First, below-threshold voluntary registrants’ remittances (a £1.46bn HMRC target in 2023–24) are outside \(R(T)\), so a threshold cut credits the full liability of newly covered firms even though roughly 43% of below-threshold firms would already be registered voluntarily (Liu et al. 2021). Second, and symmetrically, the headline treats every firm released by a threshold rise as deregistering; the anchor comparison below reports the voluntary-retention sensitivity. Because turnover is held fixed, moving the threshold from \(T_0\) to \(T\) changes revenue only through the firms whose registration status flips: \[\Delta R \;=\; R(T) - R(T_0) \;=\; -\,\operatorname{sign}(T - T_0) \sum_i w_i\, v_i \,\mathbf{1}\{\,T_{\min} \le y_i < T_{\max}\,\},\] where \(T_{\min}=\min(T,T_0)\) and \(T_{\max}=\max(T,T_0)\): lowering the threshold draws firms into the net and raises revenue, while raising it releases firms and loses revenue. The estimate is the mechanical reclassification effect and excludes the behavioural response of Section 7. I age the microdata to a given fiscal year by a cumulative nominal-growth factor; band membership is evaluated on data-year turnover against each year’s nominal thresholds, and only liabilities are aged—an approximation stated because a fully aged band would shift membership toward the pinned near-threshold density. At the £90,000 threshold this puts total VAT revenue at £200.9bn in 2025–26 and £206.1bn in 2026–27.2
I compare the static pipeline against HMRC’s costing of the April 2024 increase from £85,000 to £90,000, announced and costed at Spring Budget 2024 (HM Treasury 2024). The comparison requires the costing baseline convention: the official pre-measure baseline is not a frozen £85,000 threshold but a counterfactual path in which the freeze legislated to March 2026 ends and the threshold resumes uprating (85, 85, 87, 89, 92 thousand pounds over 2024–25 to 2028–29 in the model’s implementation), so the measure’s impact shrinks over the forecast and changes sign in 2028–29, when the counterfactual threshold overtakes £90,000. The official profile is \(-\)150, \(-\)185, \(-\)125, \(-\)50, \(+\)65 £m. My mechanical costing depends on the deregistration convention. If every released firm deregisters, the model gives \(-\)317, \(-\)323, \(-\)195, \(-\)66, \(+\)103: it reproduces the sign, the narrowing profile, and the 2028–29 turning point, but overshoots the official magnitudes throughout (\(-\)317 against \(-\)150 in 2024–25; \(-\)323 against \(-\)185 in 2025–26). If instead the Liu et al. (2021) voluntary-registration share (43%) of released-firm liability is retained, the model gives \(-\)180, \(-\)184, \(-\)111, \(-\)38, \(+\)59: above the official number in the first year, close to it in 2025–26 (\(-\)184.1 against \(-\)185), and modestly below it thereafter (Figure 4). That is the honest reading of the anchor: the official costing embeds a phased registration response (HMRC expected roughly 28,000 fewer registrants in the first year, against about 48,000 firms in the released band here, now pinned by the OBR fine-band targets), so the full-deregistration convention overshoots throughout, and the retention convention crosses the official path—above it in the first year, below it from 2026–27. Because my microdata are calibrated to HMRC VAT aggregates, even a close match would confirm internal consistency rather than out-of-sample accuracy. The sub-£1m agreement in 2025–26 is, in particular, finer than the model’s own liability level supports: the £85,000–£150,000 band’s calibrated liability sits \(9.8\%\) above its HMRC target total (footnote above), and rescaling the released band to target puts the retention figure near \(-\unicode{xA3}168\)m—the year-2 agreement is a coincidence within that error band, not model precision. The retention share is itself a convention: LLAT’s roughly \(43\%\) is a firm share among all below-threshold eligibles, applied here as a liability share among just-released firms.
I next cost a menu of further threshold locations, taking the post-reform £90,000 threshold as the baseline. This sweep runs on the 2024–25 vintage population (Section 3), whose liability-by-band calibration is the marginally stronger of the two (92.0% against 91.3%; both vintages score about 90% overall). Each row of Table 1 costs a move of the threshold to a new location for 2025–26, holding the turnover distribution fixed, so the revenue change reflects only the reclassification of firms into or out of the VAT net. Because the optimiser determines the within-band shape near the threshold (Section 3), I evaluate the reclassified band from a smooth extrapolation of the above-threshold profile rather than the raw local mass. Lowering the threshold draws firms in and raises revenue; raising it loses revenue; the threshold-cut rows inherit the registered-business frame, so they omit currently unregistered businesses that a lower threshold would newly cover. The schedule is approximately linear with a mild concavity: each £5,000 step changes revenue by roughly £230m–£265m, slightly less per step at higher thresholds because the number of firms reclassified per step falls (from about \(38{,}000\) to about \(26{,}000\)) faster than per-firm liability rises. A £30,000 rise to £120,000 costs about £1.4bn, 0.7% of the £200.9bn base, because the threshold governs many small firms whose individual liabilities are small. Figure 5 shows the sweep graphically. These are static results: they bound the first-order fiscal effect of a threshold change but abstract from the behavioural response of Section 7.
| Threshold | Revenue change (£m) | Change (%) | VAT-paying firms (000s) |
|---|---|---|---|
| £70,000 | \(+1{,}042.8\) | \(+0.52\) | \(+145.7\) |
| £75,000 | \(+776.2\) | \(+0.39\) | \(+107.2\) |
| £80,000 | \(+513.6\) | \(+0.26\) | \(+70.1\) |
| £85,000 | \(+254.8\) | \(+0.13\) | \(+34.4\) |
| £90,000 (baseline) | \(0\) | \(0\) | \(0\) |
| £95,000 | \(-250.9\) | \(-0.12\) | \(-33.0\) |
| £100,000 | \(-498.0\) | \(-0.25\) | \(-64.6\) |
| £105,000 | \(-741.1\) | \(-0.37\) | \(-94.8\) |
| £110,000 | \(-980.4\) | \(-0.49\) | \(-123.7\) |
| £115,000 | \(-1{,}215.7\) | \(-0.61\) | \(-151.2\) |
| £120,000 | \(-1{,}447.1\) | \(-0.72\) | \(-177.3\) |
Revenue impact (£m)
Change in VAT-paying firms (000s)
The threshold sweep above moves the location of a single cutoff. The same static principle—re-applying a tax rule to the fixed turnover distribution, with no behavioural response—also costs reforms that change the shape of the schedule rather than merely its position. I consider four reforms and report their first-order revenue effects against a common baseline: the £85,000-notch / £184.8bn 2023–24 base. (This 2023–24 data-year base differs from the £200.9bn figure in the sweep above, which runs on the 2024–25 vintage aged to 2025–26 at the post-reform £90,000 threshold; the schedule reforms are costed in the data year to match the bunching and dominated-region analysis. Both bases are theoretical-liability (VTTL-style) aggregates in the sense of the footnote above, not cash-receipts lines.) For each reform I recompute the weighted sum of liabilities \(R = \sum_i w_i\, v_i\) under the new schedule, holding each firm’s turnover fixed, and difference it against the baseline.
Raising the threshold from £85,000 to £100,000 releases the firms with turnover in \([\unicode{xA3}85{,}000, \unicode{xA3}100{,}000)\)—about \(111{,}900\) firms—from the VAT net, for a static revenue change of \(-\)£753m on the common base. Because the released band lies above the data-year threshold, its liabilities are directly calibrated (and its lower reach pinned by the OBR fine-band targets), so the band is summed directly; the sweep’s smooth-counterfactual method, which irons out the near-threshold concentration by construction, gives \(-\)£697m and \(100{,}700\) firms, 7% and 10% below the direct sum—a gap that now measures the real, target-pinned concentration of firms just above the threshold rather than an optimiser artifact. On this common base the level move carries the largest static cost of the menu and leaves the dominated region in place, relocating it to the new cutoff (Section 6); the taper removes the dominated region at a static cost of \(-\)£520m. Which matters more depends on what the reader wants from the reform—the dominated-region metric indexes misallocation only, and Section 8 discusses what it omits.
A graduated taper replaces the discrete notch with an effective rate that phases linearly from \(0\%\) at £85,000 to the full \(20\%\) at £105,000. Firms in the band remit a tapered rather than a full rate on their (fixed) turnover; mechanically this reduces revenue by \(-\)£520m, with about \(140{,}800\) firms in the £85,000–£105,000 band facing a lower effective rate.
A banded reduced rate charges firms in the £85,000–£105,000 band a reduced rate \(\tau_{\text{low}}\in\{10\%,15\%\}\) in place of the standard \(20\%\), with \(20\%\) retained above £105,000. Statically this reduces revenue by \(-\)£484m at \(\tau_{\text{low}}=10\%\) and by \(-\)£242m at \(\tau_{\text{low}}=15\%\), again with about \(140{,}800\) firms in the band facing the lower rate.
| Reform | Lever | Static revenue (£m) |
|---|---|---|
| Raise threshold to £100,000 | Location | \(-753\) |
| Graduated taper $[85,105 | {k}]$ Shape ( | phase-in) \(-520\) |
| Reduced rate \(10\%\) $[85,105\ | text{k}]$ Rate (s | tep) \(-484\) |
| Reduced rate \(15\%\) $[85,105\ | text{k}]$ Rate (s | tep) \(-242\) |
These static figures bound the first-order fiscal effect of each reform but abstract from the behavioural response of Section 7. One feature of the shape reforms is not visible in a static costing: unlike a threshold move, which transports the notch to a new cutoff, the taper and the reduced rate change the size of the notch itself—the discrete drop in net revenue at the threshold and the empty band of turnover it induces. Quantifying that change requires the structural object I develop in Section 6, where I derive the notch’s dominated region and show analytically how each of these reforms alters it.
This section applies the polynomial counterfactual-density estimator of Chetty et al. (2011) and Kleven and Waseem (2013)—as applied to UK VAT data by Liu et al. (2021)—to the synthetic firm population. Its finding is about the data, not about firm behaviour: every sub-band feature of a band-calibrated synthetic population is exactly what its targets put there, and this section demonstrates the point in three directions. The 2023–24 population’s near-threshold density is calibrated to the OBR’s published £1,000-band profile (Section 3), and the estimator reads that profile back as excess mass—target inheritance, not behaviour. A placebo that regenerates the population without those targets returns zero. And the 2024–25 vintage, for which no fine bands are published, produces a very large spurious excess-mass signal sitting exactly on an HMRC calibration-band edge. Together with a power (recovery) test, these results establish the section’s transparency claim: band-calibrated synthetic data support no bunching inference in either direction, and any sub-band feature of such data should be traced to its calibration targets before it is read as economics. The behavioural evidence on UK VAT bunching remains the administrative estimate of Liu et al. (2021).
The section also documents a correction. An earlier draft of this paper reported an excess mass of \(8{,}712\) firms at £85,000 that appeared to “reproduce” the density step documented in administrative data. That excess mass was an artifact of a liability mis-scaling in the generator (net liability set to value added itself rather than the standard rate applied to value added, inflating per-firm liabilities roughly fivefold; Section 3, footnote), which distorted how the weight optimiser distributed mass near the threshold. Correcting the scaling removed the step entirely; the near-threshold profile is now supplied by the OBR targets instead, as published data rather than an optimiser side-effect.3 A plausible-looking, literature-adjacent bunching statistic was manufactured by an accounting bug in the data layer and survived until the liability accounting was audited. I report this openly because it is the sharpest available demonstration of the section’s thesis.
Following Liu et al. (2021), I group firms into fine turnover bins of £1,000 and estimate a counterfactual density—the distribution that would obtain absent the registration notch—by polynomial regression on bin counts, excluding a manipulation window around the threshold \(y^{*}\): \[c_j \;=\; \sum_{l=0}^{q} \beta_l \, (y_j)^{l} \;+\; \sum_{i \in [y^{*}_{-},\, y^{*}_{+}]} \gamma_i \, \mathbb{1}\{j = i\} \;+\; \varepsilon_j ,\] where \(c_j\) is the number of firms in bin \(j\), \(y_j\) is the distance of bin \(j\) from the threshold, \(q\) is the polynomial order, and \([y^{*}_{-},\,y^{*}_{+}]\) is the excluded window. The fitted polynomial is then projected into the excluded region to recover the counterfactual density \(g(y)\), with a mass-conservation restriction that locates the upper edge of the manipulation window endogenously (excess mass below the threshold must equal missing mass above it; Appendix 9.2). Excess bunching is the integrated difference between the observed and counterfactual densities just below the threshold, \[B(y^{*}) \;=\; \int_{y^{*}-\Delta y^{*}}^{y^{*}} g(y)\,\mathrm{d}y \;\approx\; g(y^{*})\,\Delta y^{*},\] the approximation holding when \(g\) is locally smooth. The bunching ratio \(b(y^{*}) = B(y^{*})/g(y^{*}) \approx \Delta y^{*}/y^{*}\) is the width of the bunching region expressed as a fraction of the threshold—a normalised ratio, not an elasticity. The estimator delivers the raw geometry of any density step at the threshold; whether that step reflects firm behaviour is a separate question. The estimator reports geometry only (\(b\), excess mass \(E\), and the mass-conservation objects): no elasticity is computed from it, both because the identification critique of bunching designs applies (Blomquist et al. 2021; Bertanha et al. 2023) and because on these data there is no behavioural mass to convert.
Applied to the 2023–24 population at the £85,000 threshold, the estimator returns an excess mass of \(E = 7{,}933\) weighted firms: it reads back the bunching profile that the OBR £1,000-band targets place in the file (Section 3). This is target inheritance, not behaviour—the generator contains no firm location-choice—and the placebo below proves it by deletion. On the Liu et al. (2021) normalisation the target-inherited excess is \(b_{\mathrm{LLAT}} = 1.91\), about 1.4 times their administrative \(1.361\); the two are not comparable estimates of the same object—LLAT’s sample ends in 2014–15, the profile here is interpolated OBR data whose cross-threshold level is frame-implied (Section 3), and the readback varies severalfold across estimator specifications (Appendix 9.2)—so the distance from \(1.361\) carries no information. Mass conservation fails to bind at £85,000: the marginal-buncher walk is censored at the search bound (\(y_R = \unicode{xA3}100\)k) with missing mass \(\Delta_R = 5{,}761 < E\), all of it accruing beyond the £90,000 fine-target boundary—a calibration-band edge inside the estimator’s own window, not administrative-style missing mass (Appendix 9.2).
Applied to the 2024–25 population at the £90,000 threshold—which has no fine targets—the same estimator returns \(E = 11{,}072\) and \(b_{\mathrm{LLAT}} = 1.23\): within ten per cent of Liu et al. (2021)’s administrative \(1.361\). Nothing behavioural produced it. It is the coarse £90,000 band edge, where the £1–£90,000 and £90,000–£150,000 targets meet, read through the estimator; its magnitude swings four-fold across estimator specifications (\(E\) between roughly \(8{,}600\) and \(36{,}700\) over the degree–window grid, Appendix 9.2) and the manipulation-window edge sits at the search bound. The pair of results is the section’s central exhibit (Figure 6): where fine targets exist, the estimator recovers exactly the shape they encode; where they do not, band geometry alone can manufacture a bunching statistic in the publishable range of the literature—a researcher looking to “validate” synthetic data against Liu et al. (2021) would have found agreement here, and it would have meant nothing.
2023–24 population, £85,000: target-inherited profile
2024–25 population, £90,000: band-edge artifact
The placebo diagnoses where any threshold-adjacent structure comes from by manipulating the calibration targets directly. HMRC publishes turnover counts only in coarse bands across the threshold: in 2023–24, \(678{,}350\) firms at or below £85,000 (about \(7{,}981\) per £1,000) against \(305{,}320\) between £85,000 and £150,000 (about \(4{,}697\) per £1,000), so a band-average step is present in the target before any firm is generated. In Placebo A I reweight the synthetic firms so the per-£1,000 density follows a single smooth trend across the threshold, holding total mass fixed. In Placebo B I regenerate the entire population and re-calibrate from scratch with no fine structure in the targets at all: the OBR near-threshold bins are removed and the two HMRC band counts straddling the threshold are replaced by a split implied by one common density. The results track the targets exactly. The control—calibrated to the OBR profile—returns \(E = 7{,}933\); Placebo A returns \(E = 0\) (\(b = -0.059\)); Placebo B returns \(E = 0\) (\(b = -0.073\)). Deleting the fine targets deletes the “bunching”; and where a coarse step is present in the targets—as at the 2024–25 vintage’s £90,000 band edge—the estimator reads it as (artifactual) mass. Inserting or removing structure in the targets inserts or removes “bunching” in the estimates.
A null result is only informative if the estimator could have detected a signal. I therefore inject a Kleven–Waseem-consistent behavioural response of known magnitude into the step-free (Placebo A) population: firms relocate to just below the threshold from a donor region extending to a true marginal buncher at £112,000—beyond the estimator’s exclusion window, as with genuine bunching—with relocation probability declining linearly in distance from the threshold. Scored by the headline estimator with no reference to the known null, injected excess masses of \(5{,}000\) and \(8{,}000\) weighted firms are recovered at \(49\%\) and \(62\%\) respectively, while a small \(2{,}000\)-firm injection sits near the detection floor (\(15\%\)); the estimated marginal-buncher location is attenuated toward the threshold throughout. The attenuation has a measured mechanical source: part of the missing-mass region lies outside the exclusion window, where it pulls the fitted counterfactual down. The estimator therefore under-states a true response, more severely for small ones, and does not over-state one—so the target-inherited \(E=7{,}933\) at £85,000 is, if anything, an attenuated reading of the targeted profile.4
The behavioural fact in the literature is administrative. Liu et al. (2021) (hereafter LLAT) document genuine bunching at the UK VAT registration notch in HMRC micro-data over 2004–05 to 2014–15, during which the threshold rose from £58,000 to £81,000, with an excess-bunching ratio of \(b = 1.361\) (SE \(0.202\)) around the year-specific threshold and a manipulation window of \(-\unicode{xA3}14{,}000\) to \(+\unicode{xA3}24{,}000\). My data’s \(b_{\mathrm{LLAT}} = 1.91\) is not a second estimate of their parameter: it is the OBR target profile read back through their normalisation; and the uncalibrated 2024–25 vintage’s \(b_{\mathrm{LLAT}} = 1.23\), despite sitting within 10% of their estimate, is not one either—it is band geometry.5 What does carry across is institutional context: LLAT report that almost half (roughly \(43\%\)) of eligible below-threshold firms register voluntarily, which dampens observed bunching relative to a pure notch and which motivates the voluntary-retention sensitivity of Section 4. For wider order-of-magnitude context, Office for Budget Responsibility (2023) count firms holding turnover below the frozen threshold over a wide band (rising from \(23{,}000\) in 2017–18 to about \(44{,}000\) in 2025–26, with foregone turnover of £110m rising to £350m); this is a different object from a near-threshold spike and is cited only as context.
These results sharpen, rather than weaken, the case for grounding the policy analysis in structure. The reduced-form density around the threshold cannot price reforms—it has no single point of excess mass to forward-map onto a change in the shape of the schedule—and on band-calibrated synthetic data it cannot be read behaviourally at all, in either direction. What the structure of Section 6 does deliver exactly, for every reform considered, is the notch’s dominated region and how each schedule change alters it. It is precisely because the reduced-form bunching is non-identifying here that the exact dominated region carries the distortion characterisation. The behavioural layer that re-prices each firm’s turnover under counterfactual schedules—conditional on an assumed turnover elasticity—is developed in Section 7.
The behavioural response to the registration notch is documented in administrative data by Liu et al. (2021); the synthetic data contain no bunching signal to exploit (Section 5). To characterise the notch precisely I therefore turn to its analytic structure, which yields an exact distortion measure that needs no behavioural calibration. I set out the firm-level structure following Kleven and Waseem (2013), Saez (2010), and Kleven (2016), whose bunching framework Liu et al. (2021) adapt to UK VAT. The UK registration threshold is a notch, not a kink: once turnover crosses the threshold the firm becomes liable for VAT on its entire base, not merely on the increment above it. The section’s central deliverable is the notch’s dominated region—an exact, elasticity-free property of the schedule that indexes the distortion the threshold imposes (Section 6.2)—and its analytic response to each schedule reform, complementing the static reform costs of Section 4.1. Section 7 prices the intensive-margin response to each reform conditional on an assumed turnover elasticity; nothing in this section depends on that assumption.
The model prices the registration margin for firms that remit positive net VAT. Two populations are outside its scope and are treated as such throughout: below-threshold voluntary registrants (roughly 43% of eligible firms in (Liu et al. 2021)), who have chosen registration and so face no notch at the threshold, and net-repayment traders (zero-rated exporters and similar), whose remittance is negative for rated-structure reasons the generator does not represent (Section 3).
Each firm \(i\) chooses turnover \(y_i\) to maximise after-tax profit. Let \(c_i(y_i)\) denote the cost of producing turnover \(y_i\), \(T^{*}\) the registration threshold, and \(\tau\) the standard VAT rate. In the simplest statement VAT applies to the firm’s entire turnover once it registers, so after-tax revenue is discontinuous and profit is
\[ \pi_i(y_i) \;=\; \begin{cases} \;y_i - c_i(y_i), & y_i < T^{*} \quad \text{(unregistered)},\\[4pt] \;(1-\tau)\,y_i - c_i(y_i), & y_i \ge T^{*} \quad \text{(registered)}. \end{cases} \tag{1}\]
At \(T^{*}\) after-tax revenue falls discretely by \(\tau\,T^{*}\): a firm that nudges its turnover from just below to just above the threshold loses \(\tau\,T^{*}\) of net revenue on turnover it was already earning. With the standard rate \(\tau=0.20\) and threshold \(T^{*}=\unicode{xA3}85{,}000\) this discrete drop is \(\tau\,T^{*}=\unicode{xA3}17{,}000\) on whole turnover. Real VAT taxes value added, and the corrected data give each firm net liability \(v_i=\tau(y_i-x_i)\) with deductible inputs \(x_i=\delta_i y_i\) (Section 3); profit is then \((1-\delta_i)y_i - \tilde c_i(y_i)\) unregistered and \((1-\delta_i)(1-\tau)y_i - \tilde c_i(y_i)\) registered—the deductible share multiplies both branches, so the notch geometry of Equation 1 carries over firm by firm, with the cash size of the jump scaled to \(\tau(1-\delta_i)T^{*}\) (about £6,800 at the mean value-added share of 40%). This drop—absent from a kink, where only the marginal rate changes at the threshold while average liability moves continuously—is the defining feature of a notch: it makes locating just below the cutoff strictly attractive and creates the empty band derived next.
The notch creates a range of turnover just above \(T^{*}\) that no firm chooses, because every point in it yields strictly lower profit than locating at \(T^{*}\). A firm at the threshold earns net revenue \(T^{*}\); a registered firm at \(T^{*}+a\) earns \((1-\tau)(T^{*}+a)\). Equating the two—the cost difference across this narrow band is negligible—the upper edge of the dominated region is
\[ (1-\tau)\,(T^{*}+a) \;=\; T^{*} \quad\Longrightarrow\quad a \;=\; T^{*}\,\frac{\tau}{1-\tau} . \tag{2}\]
This is the Kleven–Waseem dominated-region width (Kleven and Waseem 2013). At my parameters it gives \(a = \unicode{xA3}85{,}000\times 0.20/0.80 = \unicode{xA3}21{,}250\), so the dominated region is \((T^{*},\,T^{*}+a)=(\unicode{xA3}85{,}000,\ \unicode{xA3}106{,}250)\), shown in Figure 7.
shape. It is the width \(a\), not the count, that is exact.
Because the width \(a=T^{*}\tau/(1-\tau)\) depends on the threshold and the rate alone, the effect of each schedule reform of Section 4.1 on the dominated region is exact and follows directly from the notch arithmetic of Equation 2 and Equation 3, with no dependence on the cost specification, any behavioural elasticity, or any simulated re-optimisation. Raising the threshold to £100,000 leaves the rate unchanged and therefore relocates the region without altering its width: it is the same £21,250 band, transported to the new cutoff. A banded reduced rate does not remove the distortion. It reduces the largest single notch—the jump at \(T^{*}\) falls from \(\tau T^{*}=\unicode{xA3}17{,}000\) to \(r\,T^{*}=\unicode{xA3}12{,}750\) at a \(15\%\) band rate—but, because the band reverts to the standard rate \(\tau\) at its upper edge \(T_{1}=\unicode{xA3}105{,}000\), it introduces a secondary notch there, and hence a second dominated region the simple “\(a\) shrinks with \(\tau\)” reasoning omits. The primary region at \(T^{*}\) (the band rate \(r\) on whole turnover) has the usual width \(a=T^{*}r/(1-r)\): \(\unicode{xA3}85{,}000\times0.15/0.85=\unicode{xA3}15{,}000\) at \(15\%\) and \(\unicode{xA3}85{,}000\times0.10/0.90=\unicode{xA3}9{,}444\) at \(10\%\). The secondary region at the band top—where the schedule steps from \(r\) back to \(\tau\) on whole turnover—has width
\[ a' \;=\; T_{1}\,\frac{\tau-r}{1-\tau}, \tag{3}\]
giving \(\unicode{xA3}105{,}000\times0.05/0.80=\unicode{xA3}6{,}562\) at \(15\%\) and \(\unicode{xA3}105{,}000\times0.10/0.80=\unicode{xA3}13{,}125\) at \(10\%\). The total dominated-turnover width is therefore \(a+a'=\unicode{xA3}21{,}562\) at \(15\%\) and \(\unicode{xA3}22{,}569\) at \(10\%\)—both marginally above the \(\unicode{xA3}21{,}250\) baseline, not below it. The reduced rate splits one wide empty band into two narrower ones rather than shrinking the total: it lowers the largest single notch (to \(\unicode{xA3}12{,}750\) at \(15\%\)) at the price of adding a secondary notch of \((\tau-r)\,T_{1}=\unicode{xA3}5{,}250\). The firm counts confirm this (results/dominated_region_mass.txt): the primary regions hold about \(111{,}900\) (\(15\%\)) and \(77{,}300\) (\(10\%\)) weighted firms, the secondary regions add about \(34{,}200\) and \(65{,}500\), and the totals—\(146{,}100\) and \(142{,}900\)—sit within \(3.3\%\) of the \(147{,}700\)-firm baseline band. A graduated taper, by contrast, has no notch and hence no dominated region: marginal net revenue declines smoothly through the band rather than off a cliff, so the band of dominated turnover—and with it the discrete incentive to suppress turnover—is removed entirely. On the dominated-region metric, then, only the taper removes the distortion; the reduced rate splits it and the level move relocates it. The metric indexes the misallocation distortion only—what it omits, and how the reforms compare on cost, is taken up in Section 8. The dominated region is the paper’s exact, elasticity-free measure of the distortion the notch imposes, and the analytic statement above is exactly how much of that distortion each reform of Section 4.1 removes.
The reforms of Section 4 are costed statically: each schedule is re-applied to a fixed turnover distribution, and the revenue change is the mechanical reclassification of firms with no response in turnover. Section 5 shows that the synthetic data contain no behavioural signal from which a response could be estimated, and Section 6 supplies the one object that needs no behavioural calibration at all: the exact dominated region. What is missing from this sequence is the behavioural layer itself—a costing that lets firms re-optimise turnover when a reform changes the effective rate they face. This section supplies that layer for an assumed turnover elasticity \(e\), reported as an \(e\)-sensitivity range rather than a point estimate, and nesting the static costing exactly in its \(e\to0\) limit. Its headline is deflationary, and usefully so: on the intensive margin the static costings are robust. A pure threshold rise has exactly zero intensive-margin revenue offset at every \(e\), and the reduced-rate bands’ offsets stay under \(6\%\) of the static cost across the sweep.
I adapt the iso-elastic model of Kleven and Waseem (2013), in the sufficient-statistics tradition of Saez (2010), writing the tax on value added rather than on turnover. (Liu et al. (2021) caution that turnover-tax bunching formulae do not carry to VAT unmodified, because input costs intervene; the value-added statement below is the adaptation this requires.) In this statement—formulation A—a firm of ability \(n_i\) chooses turnover \(y\) to maximise
\[ \pi_i(y) \;=\; (1-\delta_i)\,\bigl(1-\tau\,f(y)\bigr)\,y \;-\; C(y;n_i,e), \qquad C(y;n,e)=\frac{n}{1+1/e}\left(\frac{y}{n}\right)^{1+1/e}, \tag{4}\]
where \(\delta_i\in[0,1)\) is the firm’s deductible-input share of turnover, so that value added is \((1-\delta_i)\,y\); \(\tau\) is the statutory VAT rate; and \(f(y)\in[0,1]\) is the fraction of \(\tau\) that the schedule levies on the firm once registered (\(f=0\) below the threshold, \(f=1\) for a fully registered firm under the standard rate). Here \(e>0\) is the elasticity of turnover with respect to the net-of-tax rate.
It is essential to read \(C\) correctly. Because the deductible inputs \(\delta_i\,y\) are already netted out of revenue through the \((1-\delta_i)\) factor, \(C(y;n_i,e)\) must be understood as the firm’s own-factor production cost—the labour and capital effort the firm itself supplies—and not its total cost of production. Reading \(C\) as total cost would double-count the bought-in inputs, which the value-added factor has already removed. Marginal own-factor cost is \((y/n)^{1/e}\), so the undistorted optimum is \(y=n\): ability \(n\) is the firm’s frictionless turnover. A firm registered under a flat rate (\(f\equiv1\)) chooses
\[ y^{*} \;=\; n\,\bigl[(1-\delta)(1-\tau)\bigr]^{e}, \tag{5}\]
so the retained wedge facing the firm is the product \((1-\delta)(1-\tau)\) of the value-added share and the statutory net-of-tax rate, and \(\mathrm{d}\ln y^{*}/\mathrm{d}\ln\!\bigl[(1-\delta)(1-\tau)\bigr]=e\): \(e\) is exactly the knob that governs the intensive response. Figure 8 traces the formulation-A profit and its optima across the notch and taper schedules for three deductible shares.
The reform costing needs each firm’s optimal turnover under a counterfactual schedule, given that its observed turnover is optimal under the baseline notch. Formulation A makes this a one-line object. Inverting Equation 5 at the baseline (fraction \(f_0\)) and re-solving at the reform (fraction \(f_1\) on the firm’s region of the schedule), the ability and the deductible share cancel in the ratio:
\[ y^{*}_i \;=\; y_{\mathrm{obs},i}\, \left(\frac{1-\tau f_1}{1-\tau f_0}\right)^{e}. \tag{6}\]
The intensive response depends only on \(e\), the statutory rate, and the schedule fractions—not on \(\delta_i\), not on the firm’s net remittance rate, and not on the recovered ability level. Ability recovery (\(n_i=y_{\mathrm{obs},i}\) below the threshold, \(n_i=y_{\mathrm{obs},i}/[(1-\delta_i)(1-\tau)]^{e}\) above) is an accounting anchor—it rationalises the observed allocation given \(e\) and reproduces the baseline by construction—but no reform number below depends on it.
The margin priced here is purely intensive: each firm re-optimises within the region of the reform schedule that contains its observed turnover, per Equation 6, clipped at its region’s edges. Relocation across a notch—bunching below a threshold from above, or de-registering past a band edge—is the extensive margin, which is the separate analytic object of Section 6 (the dominated region and the marginal buncher) and is deliberately out of scope: the synthetic data contain no information about it (Section 5), and pricing it would require assumptions about re-bunching dynamics the paper has declined to make. Confining the response to the firm’s own region also makes the \(e\to0\) limit collapse onto the static costing exactly, with no solver tolerance.6
The elasticity \(e\) is not identified from these data, so every behavioural magnitude below is conditional on the assumed \(e\), and the sweep is a sensitivity analysis rather than a set of forecasts. The sweep values \(e\in\{0.05,\,0.17,\,0.32\}\) are assumptions: \(0.05\) is the external anchor—Kleven and Waseem (2013) report lower-bound structural elasticities mostly in \([0.05, 0.15]\)—and \(0.17\) and \(0.32\) are chosen to span the range used in the threshold literature; the published UK VAT-base elasticities of Liu et al. (2021) (roughly \(0.09\)–\(0.14\)) fall between the anchor and the midpoint. A reader may substitute any preferred value directly. The robust, \(e\)-free headline objects remain the static costs of Section 4 and the dominated region of Section 6.
Table 3 reports the behavioural cost of each of the three flat-rate reforms—the raised threshold and the two reduced-rate bands—measured as the change in revenue relative to the £85,000-notch baseline, alongside the static figure of Table 2, for each \(e\). For context, the marginal buncher \(n_H(e)\)—the highest ability that still finds bunching at the threshold optimal, the extensive-margin object held fixed here—sits at £112,795, £127,382 and £143,527 for \(e=0.05, 0.17, 0.32\).
| Reform | Static | \(e=0.05\) | \(e=0.17\) | \(e=0.32\) | Firms (\(e{=}0.17\)) |
|---|---|---|---|---|---|
| Raise threshold to £100,000 | $- | 753$m $- | e{xA3}753$m $-\ | unicode{xA3}753$m | \(-\unicode{xA3}753\)m \(108{,}973\) |
| Reduced rate \(10\%\) $[85,105\ | text{k}]$ $-\ | unicode{xA3}484$m | $-\unicode{xA3 | }481$m $- | de{xA3}475$m \(-\unicode{xA3}467\)m \(137{,}944\) |
| Reduced rate \(15\%\) $[85,105\ | text{k}]$ $-\ | unicode{xA3}242$m | $-\unicode{xA3 | }240$m $- | de{xA3}235$m \(-\unicode{xA3}228\)m \(137{,}944\) |
Note. The graduated taper is omitted: its rate varies continuously with turnover, so its intensive response is outside the region-confined solve, and its marginal-rate channel is not represented by the flat-rate first-order condition. Its \(e\)-free figures—the static cost of \(-\)£520m and the exact removal of the dominated region—are in Sections 4 and 6.
The raise-to- £100,000 row is identical to its static cost at every \(e\), and this is a result, not a failure to respond. Partition the firms. Those with observed turnover below £85,000 face \(f=0\) before and after; nothing changes. Those released by the raise (\(y_{\mathrm{obs}}\in[\unicode{xA3}85\text{k}, \unicode{xA3}100\text{k})\)) see their rate fall to zero and expand per Equation 6 by the factor \((1-\tau)^{-e}\) (about \(+3.9\%\) at \(e=0.17\))—but their expansion is untaxed: they have left the VAT base, so the revenue loss is their full baseline remittance, exactly the static figure. Those above £100,000 face \(f=1\) before and after; nothing changes. The only candidates for a taxed response are released firms whose expansion would carry them across the new threshold, and for every such firm the new notch’s bunching condition binds (its ability lies strictly below the new marginal buncher), so it optimally stops at the threshold and remits nothing. Base broadening from a threshold rise is therefore real economic activity but generates no VAT: conditional on this model, the static costing of a level move is not an upper bound that behaviour erodes—it is the behavioural costing on the intensive margin. The simulator asserts this invariance programmatically (to £0.1m at each swept \(e\)).
A reduced-rate band lowers the rate on band firms from \(\tau\) to \(r\), so each scales up by Equation 6 with \(f_1 = r/\tau\): about \(+2.0\%\) at the \(10\%\) rate and \(e=0.17\), and about \(+1.0\%\) at \(15\%\). Remittance scales with turnover, so the offset is roughly the scale gain times the band’s (rate-reduced) remittance—£9m against a static £484m at \(e=0.17\) for the \(10\%\) band, reaching £17m (\(3.5\%\)) only at \(e=0.32\). The direction is uniform—a larger \(e\) makes the bands slightly cheaper—but the magnitude never approaches the static term. Figure 9 traces the costs across the full range of \(e\): the raise is exactly flat, the bands nearly so.
Three properties are asserted programmatically in the released code on every run: the baseline solve reproduces each firm’s observed turnover exactly (the accounting anchor is internally consistent); the \(e\to0\) behavioural cost equals the static cost to £0.1m for every priced reform; and the solved response satisfies \(\mathrm{d}\ln y^{*}/\mathrm{d}\ln(1-\tau f)=e\) to machine precision. The static costing is thus the \(e=0\) corner of the same machine, not a separate calculation.
This section completes the costing framework, and its contribution is a robustness statement. Where Section 4 prices each reform mechanically and Section 6 measures the distortion each removes, the intensive-margin layer shows that naming an elasticity barely moves the static numbers: not at all for a level move, and by under \(6\%\) for the reduced-rate bands. The behavioural uncertainty that could move these costings materially sits on the extensive margin—registration, de-registration, and location choice at the notch—which the synthetic data cannot identify and which this paper prices only through its exact analytic geometry.
This paper contributes an open, reproducible firm-level microsimulation for analysing the UK VAT registration threshold, released in full—including its synthetic-data generator—for replication. Its value is as usable policy infrastructure rather than a structural identification result: I make no claim to recover policy-invariant firm primitives, unlike Garicano et al. (2016) and Gourio and Roys (2014). The analysis is anchored to the 2023–24 data year, in which the statutory threshold was £85,000, with a second 2024–25 population (threshold £90,000) for the forward sweep. Table 4 collects the headline objects.
| Object | Symbol / definition | Value | Source |
|---|---|---|---|
| Static cost, £85k to £90k anchor | dereg. / retention conventions | \(-\unicode{xA3}323\)m / $-\unicode{ | xA3}184$m Section 4 |
| Static sweep, £90k baseline | per £5k step | £230m–£265m | Table 1 |
| Dominated region | \(a=T^{*}\tau/(1-\tau)\) | £21,250 | Section 6.2 |
| Firms in the dominated band (weighted) | — | \(\approx148{,}000\) | Section 6.2 |
| Raise to £100k, static and behavioural | invariant in \(e\) | \(-\unicode{xA3}753\)m | Sections 4, 7 |
| Graduated taper, static | — | \(-\unicode{xA3}520\)m | Section 4.1 |
Four results stand out. First, the static pipeline prices level and schedule reforms on a common £85,000 / £184.8bn base (\(-\)£753m to raise the threshold to £100,000, \(-\)£520m for a graduated taper, \(-\)£484m and \(-\)£242m for \(10\%\) and \(15\%\) banded reduced rates), and its anchor costing of the £85,000 to £90,000 move tracks HMRC’s published figure once the deregistration convention is stated: the full-deregistration convention overshoots (\(-\)£323m against the published \(-\)£185m in 2025–26, consistent with the official costing’s phased registration response), the retention convention runs close to the published 2025–26 figure (\(-\)£184.1m against \(-\)£185m; agreement finer than the calibrated liability level supports, Section 4), and both conventions reproduce the official profile’s narrowing and its 2028–29 sign change under the uprated counterfactual threshold path.
Second, the notch’s distortion is captured exactly by its dominated region, the Kleven and Waseem (2013) width \(a=T^{*}\tau/(1-\tau)=\unicode{xA3}21{,}250\)—an accounting identity rather than a new structural result, but one that is invariant to the firm’s input share under the value-added formulation and spans a populous range of roughly \(148{,}000\) weighted firms. Setting cost against distortion on this metric: raising the threshold carries the largest static cost and relocates the region without changing its width; a banded reduced rate lowers the largest single notch (from £17,000 to £12,750 at \(15\%\)) but adds a secondary notch at the band top, leaving the total dominated width essentially unchanged (£21,562 at \(15\%\), £22,569 at \(10\%\)); only the graduated taper removes the dominated region, at a static cost of \(-\)£520m against the \(10\%\) band’s \(-\)£484m, which removes none of it. The metric indexes misallocation only; the limitations paragraph below states what it omits.
Third, the synthetic data carry a transparency lesson: sub-band structure in band-calibrated microdata is exactly what its targets put there, so it supports no bunching inference in either direction. The 2023–24 population’s near-threshold profile is calibrated to the OBR’s published £1,000-band data, and the estimator reads it back as target-inherited excess mass (\(E=7{,}933\)); a placebo that regenerates without the fine targets returns zero; and the 2024–25 vintage, which has no fine targets, shows the coarse £90,000 band edge alone yielding an excess-bunching statistic within 10% of the administrative estimate—coincidental, specification-unstable, and behaviour-free. A redesigned power test confirms the estimator detects moderate genuine signals (recovering half to two-thirds of 5,000–8,000-firm injections, attenuated, never inflated). An earlier draft of this paper reported “reproducing” the administrative density step at £85,000 when a liability mis-scaling manufactured it; the paper now documents the artifact, its correction, and the replacement of optimiser side-effects with published fine-band targets as the sharpest demonstration of the lesson. The behavioural evidence on threshold bunching remains the administrative estimate of Liu et al. (2021).
Fourth, the behavioural layer’s contribution is a robustness statement: on the intensive margin, naming an elasticity barely moves the static costings. A pure level rise has exactly zero intensive-margin offset at every \(e\)—released firms leave the VAT base, so their expansion is untaxed, and boundary-crossers optimally bunch below the new threshold—and the reduced-rate bands’ offsets stay under \(6\%\) of the static cost across \(e\in\{0.05,0.17,0.32\}\) (the \(10\%\) band falls from £484m to £475m at \(e=0.17\)). The graduated taper is priced only statically: its continuously varying rate is outside the flat-rate machinery. The behavioural uncertainty that could move these costings materially sits on the extensive margin—registration and location choice at the notch—which these data cannot identify and the paper prices only through its exact analytic geometry.
The model abstracts from several margins, each stated where it binds. It has no firm location choice and no extensive-margin response; voluntary registration enters only as a registration-rate parameter and a retention sensitivity, not a choice model; net-repayment traders are outside the liability model (Section 3); the statutory rolling registration test is approximated by a data-year turnover test on total rather than taxable turnover; the near-threshold profile rests on two stated assumptions (the 2023–24 interpolation between the OBR’s 2019–20 outturn and 2025–26 frozen-threshold projection, and the shape-only, side-consistent treatment of the OBR bins, whose chart universe differs from the ONS frame); and the calibrated revenue bases sit about \(5.8\%\) above the HMRC liability-by-band target sums, with liability-by-band calibration scores of \(91.3\%\) and \(92.0\%\) (Section 4, footnote). The dominated-region width indexes the misallocation distortion only and omits the compliance and administrative-cost saving that a higher threshold delivers by removing small firms from VAT—the classic rationale for a high registration threshold (Keen and Mintz 2004)—so the efficiency comparison is partial; it likewise books no compliance saving to the taper or reduced rate, which keep every band firm inside the VAT system, and a taper’s administrability (an effective rate that depends on annual turnover cannot be charged on invoices, so a real design would operate as remittance-side relief) is a design question the costing does not address. The menu reforms are costed only in aggregate: because VAT liability by sector is excluded from the calibration optimiser (Section 3), the framework does not yet report the sectoral incidence of a taper or reduced rate, a cut of obvious policy interest, and the split of band-total liability across sub-band turnover—which the anchor and relocation costings use—rests on the generator’s common input-share distribution rather than on sector-composition data. The analysis is partial-equilibrium, abstracting from pass-through into prices (Benedek et al. 2015; Copestake et al. 2026). The chief outstanding empirical task is to identify the turnover elasticity from administrative micro-data, which would turn the \(e\)-sensitivity range into a point estimate; a welfare/MVPF accounting would net compliance costs against the misallocation distortion and support revenue-neutral schedule designs. Above all, validation against administrative firm-level data with genuine location choice would replace the present internal-consistency checks with out-of-sample tests.
This appendix collects supporting material in two parts: data construction and calibration robustness (Appendix 9.1); and statistical inference for the reduced-form bunching results (Appendix 9.2), including the placebo and power tests establishing that band-calibrated synthetic data cannot support bunching inference in either direction.
Section 3 describes the two-stage synthetic-population pipeline, parameterised by the registration threshold, and its calibration to HMRC and ONS aggregates. This appendix records the implementation detail deferred from that account.
Each of the approximately \(2.94\) million firm rows carries a turnover, an input expenditure, a net VAT liability (\(v_i = 0.20\,(y_i - x_i)\)), a calibration weight, and a registration status; the weights are chosen so the weighted totals reproduce the official targets while the population sums to about \(2.5\) million firms. Because the population is calibrated to the HMRC aggregates, agreement between my weighted outputs and those aggregates is a check of internal consistency, not external validation: any target used in calibration is reproduced by construction. The headline calibration accuracies (89.2% and 92.8% for the two vintages) and the anchor comparison of Section 4 are therefore consistency checks on the pipeline, not out-of-sample tests of firm behaviour.
The reported calibration score is a bounded summary statistic. For each dimension I compute \(\max\{0,1-\lvert \widehat{T}-T\rvert/\lvert T\rvert\}\), so a dimension missed by more than 100 percent receives zero rather than a negative score. The headline “overall” number is the unweighted mean of the five dimensions used as calibration targets—HMRC turnover bands, ONS population, employment bands, HMRC sector counts, and VAT liability by turnover band. VAT liability by sector is excluded from that mean and reported separately as an informational diagnostic, because the current generator does not calibrate sector-specific input/output VAT structure.
The official target tables are being moved into PolicyEngine Ledger, with Populace providing the shared synthetic-population generator. I therefore keep the paper’s archived processed CSVs as the reproduction source for the reported results, and treat the pinned Populace/Ledger path as an auditable migration check rather than a silent replacement. The snapshot used here is PolicyEngine/ledger pull request 67 at commit cd98b5c and PolicyEngine/populace pull request 223 at commit fa20daf; both upstream pull requests merged on June 30, 2026. For the 2024–25 vintage, the Ledger-backed targets match the paper’s processed numeric inputs exactly after dropping presentation-only labels, totals, and the HMRC “Unknown” column that the generator does not calibrate: six normalized source tables checked, zero mismatches, and maximum numeric difference zero. (The “2024–25 vintage” here denotes the HMRC/ONS statistical release labelled 2024–25, the same target tables the paper’s 2024–25 population calibrates to.) A 1,000-iteration Populace run from those targets generated 2,946,015 firm rows; its validator reported a 93.8 percent overall score against the paper validator’s then-90.5 percent, a comparison that is not like-for-like (different band sets and sector-target definitions). Two caveats now attach to that pinned snapshot. First, it predates this paper’s net-liability correction (Section 3, footnote): both generators at the pin date shared the mis-scaled \(v_i = y_i - x_i\) convention, so the experimental Populace firm generator inherits the same defect, which I have reported upstream; the target-parity result (zero mismatches) is unaffected, since it concerns the input tables, not the generated population. Second, the paper validator’s scores on the corrected build are 89.2 percent (2023–24) and 92.8 percent (2024–25). The comparison table and structured provenance are reproduced by firm-microsim-populace-ledger and checked into results/populace_ledger_comparison.txt and results/populace_ledger_provenance.json.
The no-VAT counterfactual density of Section 5 is fitted by polynomial regression on the observed density outside a manipulation window around the threshold. The baseline excludes \([T^{*}-15\text{k}, T^{*}+15\text{k}]\) and fits a degree-7 polynomial. Widening the window guards against contamination of the counterfactual by the behavioural response but reduces the data available for the fit; raising the degree improves in-sample fit but risks over-fitting the tails. The sensitivity of the excess-mass and bunching-ratio estimates to window width and degree, with associated uncertainty, is reported in Appendix 9.2.
This appendix documents the inference procedure for the reduced-form bunching results of Section 5: the mass-conservation constraint that disciplines the counterfactual density, why no standard errors are reported, sensitivity to the polynomial degree, the exclusion window, and the generator seed, the placebo tests, and the power (recovery) test. All numbers are produced by the released estimator run on the regenerated weighted synthetic populations for both vintages; the machine-readable artifacts are results/bunching_inference.txt, results/placebo_bunching.txt, results/recovery_bunching.txt and results/seed_sensitivity.txt. The estimator reports the geometry only—the bunching ratio \(b\), its Liu et al. (2021)-normalised counterpart \(b_{\mathrm{LLAT}}\), the excess mass \(E\), and the mass-conservation objects (\(\Delta_R\), \(y_R\))—and computes no elasticities: the mapping from excess mass to an elasticity is not identified without functional-form assumptions (Blomquist et al. 2021; Bertanha et al. 2023), and on these data there is no behavioural mass to convert in the first place.7
The canonical bunching estimator of Chetty et al. (2011) and Kleven and Waseem (2013) imposes that the counterfactual density conserves mass: every firm removed from above the threshold reappears as excess mass below it. An unconstrained polynomial fit on bins outside a hand-set window computes the excess mass \(E\) (below \(T^{*}\)) and the missing mass \(\Delta_{R}\) (above \(T^{*}\)) independently, with no guarantee that \(E=\Delta_{R}\) and no endogenous location of the marginal buncher. I impose the constraint in two ways. First, the counterfactual polynomial in \((y-T^{*})\) is fit by an iterated procedure in which the density of the bins above the excluded region is rescaled until the fitted counterfactual integrates to the observed total mass over the estimation range \([\unicode{xA3}20\text{k},\unicode{xA3}140\text{k}]\) (the regressor is normalised to \([-1,1]\) for numerical stability of the high-degree basis). Second, the upper edge of the manipulation region—the marginal buncher \(y_{R}\)—is located endogenously by integrating missing mass upward from \(T^{*}\) until it equals the excess mass below, \[\int_{T^{*}}^{y_{R}} \bigl[f^{\mathrm{cf}}(y)-f^{\mathrm{obs}}(y)\bigr]_{+}\,dy \;=\; \int_{T^{*}-W}^{T^{*}} \bigl[f^{\mathrm{obs}}(y)-f^{\mathrm{cf}}(y)\bigr]_{+}\,dy \;=\; E .\] The constraint’s behaviour differs by vintage, diagnostically so. On the 2024–25 population it is satisfied in the interior (\(y_R = \unicode{xA3}101\)k, \(\Delta_R = 11{,}452 \ge E\)): the band-edge artifact is not a local spike. On the 2023–24 population it fails to bind: the walk is censored at the search cap (\(y_R = \unicode{xA3}100\text{k} = T^{*}+W\)) with \(\Delta_R = 5{,}761 < E = 7{,}933\), and decomposing the missing mass locates all of it in \((\unicode{xA3}90\text{k},\unicode{xA3}100\text{k}]\) and none in the OBR-pinned \((\unicode{xA3}85\text{k},\unicode{xA3}90\text{k}]\) segment. The counterfactual exceeds the observed density only beyond the £90,000 fine-target boundary—the same calibration-band-edge mechanism as the 2024–25 exhibit, operating inside the headline estimator’s own window—so the reported \(y_R\) is a property of the target boundary, not evidence of administrative-style missing mass.
Earlier drafts attached bootstrap standard errors to the bunching statistics (resampling firm rows with replacement and re-running the estimator). Those numbers are not reported anywhere in this paper, because the exercise has no estimand: the synthetic file is a deterministic construction from published targets, not a sample from a population, and row-resampling dispersion scales as \(1/\sqrt{n}\) in a row count the analyst chooses. A standard error that can be made arbitrarily small by generating more synthetic rows conveys precision that does not exist. Two reportable quantities replace them. First, the specification grids below, which show how the point estimates move with the estimator’s own knobs. Second, generator-seed sensitivity (results/seed_sensitivity.txt): regenerating the full 2023–24 population under three seeds moves \(E\) by \(\pm161\) firms (\(\pm2\%\)), \(b_{\mathrm{LLAT}}\) by \(\pm0.02\), the raise-to-£100k cost by \(\pm\unicode{xA3}2.1\)m, and the registered base by \(\pm\unicode{xA3}0.08\)bn—the pipeline’s actual Monte Carlo error. Table 5 reports the point estimates. Two conventions matter for reading them. Bins are centred on integer £k values (edges at half-integers), so the bin containing the threshold spans £84.5k–£85.5k and its mass is attributed to the above side, attenuating \(E\) and \(\Delta_R\) relative to edge-aligned binning. And \(E\) sums the positive part of the observed-minus-counterfactual gap in the excluded window below the threshold; the net gap there is \(1{,}570\), so the statistic reads gross rather than net excess.
| Vintage | \(E\) | \(\Delta_R\) | \(b_{\mathrm{LLAT}}\) | \(y_R\) (£k) |
|---|---|---|---|---|
| 2023–24 (£85k, OBR fine targets) | \(7{,}933\) | \(5{,}761\) | \(1.91\) | \(100.0\) |
| 2024–25 (£90k, coarse bands only) | \(11{,}072\) | \(11{,}452\) | \(1.23\) | \(101.0\) |
Table entries for a grid of polynomial degree \(\in\{5,6,7,8\}\) and symmetric exclusion window \(\in\{\unicode{xA3}10\text{k},\ldots,\unicode{xA3}25\text{k}\}\) are in results/bunching_inference.txt; the pattern matters more than the cells. On the 2023–24 population the target-inherited excess ranges over roughly \(6{,}700\)–\(9{,}800\) in the interior cells (degrees 6–7, windows £10–15k; 4 of the 16 grid cells), inflating in the corner cells where the polynomial confuses the profile’s shoulders with its counterfactual. On the 2024–25 population the band-edge artifact swings four-fold across the grid (\(E\) between roughly \(8{,}600\) and \(36{,}700\); \(b_{\mathrm{LLAT}}\) between \(0.94\) and \(4.46\)) while its headline cell lands within 10% of the administrative estimate: neither specification stability nor literature-proximity certifies a synthetic bunching statistic—a caution for any bunching exercise on calibrated synthetic data.
The placebo makes the target-inheritance mechanism precise by manipulating the calibration targets and re-running the estimator (Table 6; Section 5.3 describes the constructions). The estimates track the targets exactly: the control, calibrated to the OBR fine-band profile, returns \(E=7{,}933\); Placebo A (smooth reweighting) returns \(E=0\); and Placebo B (no fine targets, un-stepped coarse bands, full regeneration and recalibration) returns \(E=0\). Removing structure from the targets removes structure from the estimates because the targets are the only place sub-band structure can come from: the generator has no firm location-choice.
| Population | \(b_{\mathrm{LLAT}}\) | \(E\) |
|---|---|---|
| Actual (control, OBR fine-band targets) | \(1.911\) | \(7{,}933\) |
| Placebo A (reweight density smooth) | \(0.000\) | \(0\) |
| Placebo B (no fine targets, un-stepped) | \(0.395\) | \(0\) |
The placebo establishes specificity: no false positive on step-free data. The complementary property is power. Starting from the step-free Placebo A population, I inject a Kleven–Waseem-consistent behavioural signal of known magnitude: firms in the donor region \((T^{*}, \unicode{xA3}112\text{k}]\) relocate to \([\unicode{xA3}80\text{k}, T^{*})\) with probability declining linearly from the threshold to zero at the true marginal buncher \(y_R^{\mathrm{true}}=\unicode{xA3}112\)k, so the missing-mass region extends £12k beyond the exclusion window—as with genuine bunching—and the counterfactual fit confronts genuinely depressed above-window bins. Recovery is scored with the headline estimator alone (Table 7).8
| \(E_{\mathrm{true}}\) (injected) | \(E_{\mathrm{hat}}\) | Recovery | \(\hat y_R\) (true 112) |
|---|---|---|---|
| \(2{,}000\) | \(298\) | \(14.9\%\) | \(85.8\) |
| \(5{,}000\) | \(2{,}445\) | \(48.9\%\) | \(91.9\) |
| \(8{,}000\) | \(4{,}918\) | \(61.5\%\) | \(100.0\) |
Recovery rises with signal size—half to two-thirds of moderate (5,000–8,000-firm) injections, while a small 2,000-firm injection sits near the detection floor—and the estimator under-states the truth rather than inflating it. The attenuation has a measured mechanical source: the depressed bins between the window top (£100k) and the true marginal buncher (£112k) enter the counterfactual fit and pull it down, truncating measured missing mass and drawing \(\hat y_R\) toward the threshold. The implication is the one Section 5 turns on: the target-inherited excess on the 2023–24 data is read with the same attenuation a genuine signal would face, so the estimator, if anything, under-states the targeted profile; and no attenuation argument can turn the 2024–25 band-edge artifact into behaviour. Both tests are part of the released code and are fully reproducible.
An earlier draft of this generator set \(v_i = y_i - x_i\) with \(\rho_i \in [0.6, 1.5]\)—omitting the factor \(\tau\), so per-firm liabilities were about five times too large and a large minority of firms carried negative value added. The weight optimiser compensated on the band totals, leaving aggregate calibration scores similar but distorting the near-threshold density; Section 5 documents the consequence and the correction.↩︎
These aggregates are weighted sums of the synthetic registered population’s net VAT liabilities aged forward—a theoretical-liability (VTTL-style) base, not the published cash-receipts line (around £169bn in 2023–24; OBR forecast near £180bn for 2025–26), which is reduced by refunds, the VAT gap, and collection timing. The unaged bases also sit about 5.8% above the HMRC liability-by-band target sums they are calibrated toward (£184.8bn against £174.6bn above-threshold in 2023–24; £190.7bn against £180.4bn in 2024–25, and by \(9.8\%\) in the £85,000–£150,000 band specifically): the multi-target optimiser trades the liability totals off against the count targets, and the liability-by-band calibration scores (91.3% and 92.0%) bound the shortfall. Band-to-band differences, which drive every costing in the paper, are less exposed than the level. The growth factors are cumulative from 2023–24 and are applied to both vintages; to the extent the 2024–25 vintage already embodies 2024–25 nominal levels, the aged 2024–25 aggregates overstate nominal growth by up to one year’s factor (\(\approx\)3%), within the calibration error just described.↩︎
The defect and its diagnosis are recorded in the replication repository (issue #15); the near-threshold targets are issue #23. A second correction followed: the below-threshold HMRC liability total, which belongs to voluntary registrants the liability model does not represent, was initially calibrated against the whole below-threshold population, draining weights around the fine-target window; on an intermediate build it also manufactured a spectacular (\(E\approx196{,}000\)), specification-robust artifact at the £90,000 band edge. It is now an informational diagnostic (Section 3). Aggregate calibration scores barely moved under any of these changes—the optimiser hit the band totals every time—which is itself part of the lesson: aggregate fit does not discipline sub-band shape.↩︎
An earlier version of this exercise relocated mass entirely within the exclusion window, which the counterfactual fit cannot see by construction; its reported “\(\sim\)90% recovery” measured deposit bookkeeping, not estimator power. The redesigned test is in analysis/recovery_bunching.py.↩︎
The earlier draft’s \(b_{\mathrm{LLAT}} = 0.94\), read against LLAT’s \(1.361\), illustrated how a band-inherited artifact can land in the publishable range of the literature; that artifact’s disappearance under the liability correction, and its replacement by an explicit published-data target, are reported above.↩︎
An earlier version of this simulator computed the registered optimum by damped fixed-point iteration on the discontinuous schedule. For any firm whose ability straddles a reform notch that iteration has no fixed point—the iterate oscillates across the notch and the reported turnover is an artifact of the iteration count. The reported behavioural offsets in the earlier draft (for example \(-\unicode{xA3}292\)m against a static \(-\unicode{xA3}508\)m for the raise at \(e=0.17\)) were dominated by this artifact. The corrected, region-confined solve is closed form; the replication repository documents the change.↩︎
Earlier versions of the estimator also printed a CES substitution elasticity and two wedge-normalised turnover elasticities, the latter inheriting an ad hoc \(\tau/2\) normalisation from a deleted smoothed-schedule model. These outputs were labelled non-identified but were still reported; they have been removed from the code and the paper.↩︎
An earlier version of this exercise relocated mass entirely within the exclusion window. The fitted counterfactual is identical across such injections by construction (it is fit outside the window and rescaled by total mass, which relocation preserves), so that version measured deposit bookkeeping, not power, and its stated attenuation mechanism was wrong. The redesigned test supersedes it.↩︎