Hall-Like Transversal Stress and Sandpile Criticality on Real Production Networks

Abstract: This paper develops a Hall-Sandpile model of economic instability that combines a Hall-like transversal stress mechanism with sandpile threshold dynamics on a real production-network substrate. In analogy with the physical Hall effect, where exposed flows under an external field generate stress in a transversal direction, we model economic shocks as fields that act on flow-intensive, low-redundancy, low-capacity nodes and produce systemic stress through a multiplicative conversion function. The accumulated stress drives a discrete toppling rule and an avalanche dynamics whose effective activation threshold declines with transversal exposure. The model is calibrated on annual World Input--Output Database (WIOD) production networks for 2000--2014 and simulated on the 2014 substrate (2{,}283 country--sector nodes) under three alternative propagation normalisations to avoid mechanical near-criticality from row-stochastic operators. Controlled Monte Carlo experiments over external field intensity and redundancy stress generate four ordered regimes: stable absorption, latent fragility, critical transition, and avalanche regime. Mean avalanche size and the probabilities of finite-size systemic events $\Pr(S!\geq!5)$, $\Pr(S!\geq!10)$ and $\Pr(S!\geq!20)$ rise jointly with field intensity and redundancy stress. Tail diagnostics show regime-dependent thickening of the avalanche distribution, but the estimated tail indices remain too high to interpret as evidence of universal power-law criticality. The contribution is therefore a finite-size, real-network description of how transversal stress activates structural fragility, not a claim of self-organised criticality in the global economy.

• The paper develops a Hall-Sandpile model that combines multiplicative stress conversion with dissipative network propagation, showing how field intensity, redundancy, and capacity alter effective toppling thresholds.
• Simulations on the 2,283-node WIOD production network produce four stable regimes, with mean avalanche size rising from 0.084 to 5.811 as field intensity and redundancy stress jointly increase.
• The paper finds regime-dependent tail thickening with exponents near 6, but rejects a self-organised criticality claim because finite network size, dissipation, and simultaneous driving and relaxation constrain the dynamics.

Motivation and positioning

This paper develops a "Hall-Sandpile" model of economic instability, embedding a Hall-like transversal stress mechanism into dissipative sandpile dynamics on the observed World Input-Output Database (WIOD) production network. The central conceptual move is that economic systems do not merely transmit shocks along network edges but convert shocks across dimensions: an external field $B_t$ acting on flow-intensive, low-redundancy, low-capacity nodes generates stress of a different nature than the original perturbation. This contrasts with the predominantly longitudinal contagion literature in production networks and financial cascades, where perturbations enter a node and propagate along weighted directed edges. The Hall analogy is explicitly structural rather than physical: $H_{i,t}$ is not a voltage, and the author is careful to disclaim Lorentz-force kinematics.

The paper is complementary to the author's own curvature-based "Sandpile Economics" programme, which locates structural fragility in negative-curvature concentrations; here the emphasis shifts from the shape of the substrate to the loading on it. The substrate itself is not synthetic: the simulation runs on the 2014 WIOD network of 2,283 country-sector nodes, with the 2000-2014 panel used for indicator construction.

Model structure and theoretical results

The model has four blocks. First, three alternative propagation operators are constructed from the WIOD transaction matrix: a row-share matrix (row-stochastic, spectral radius near one by construction), a leakage-adjusted dissipative matrix used as the main operator, and a max-row-normalised matrix preserving row-total heterogeneity. Second, the Hall-like stress-conversion function

$$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$

multiplicatively combines field intensity, relative flow intensity $I_{i,t}$, redundancy $D_{i,t}$ and absorptive capacity $C_{i,t}$. Third, stress accumulates with dissipation $\delta$, idiosyncratic shocks $\alpha x_{i,t}$, network propagation $\beta$, and Hall loading $\gamma$, with toppling at node-specific thresholds. Fourth, avalanche size $H_{i,t}$0 counts toppling nodes per period.

Three propositions organise the comparative statics. Proposition 1 derives the Hall-adjusted threshold $H_{i,t}$1, so transversal stress lowers the effective activation threshold without altering topology. Proposition 2 establishes field-resistance complementarity: the cross-partial $H_{i,t}$2. Proposition 3 states that avalanche regimes require joint increases in field intensity and redundancy stress, producing a positively sloped phase frontier. A lemma guarantees boundedness: under the contraction condition $H_{i,t}$3 (satisfied with $H_{i,t}$4 versus a bound of $H_{i,t}$5 given $H_{i,t}$6), the stress process has a unique stationary distribution and avalanches are finite-size bounded by $H_{i,t}$7. This corollary matters for the interpretation of all subsequent tail results: heavy tails, if any, are truncated at the network size.

Spectral diagnostics and exposure concentration

A methodological contribution is the demonstration that near-critical avalanche behaviour cannot be a mechanical artefact of normalisation. The row-share spectral radius lies in $H_{i,t}$8 across 2000-2014, but this is a property of row-stochastic construction, not of the economy. The leakage-adjusted radius is sharply dissipative at $H_{i,t}$9-$$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$0, and the max-row-normalised radius rises from 0.219 (2000) to 0.317 (2014). All avalanche results are reported to arise from dissipative propagation plus Hall loading, not a unit-radius operator.

The relative Hall-like exposure $$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$1 is extremely concentrated. The mean stays near $$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$2 while the 2014 maximum reaches 0.317 — a single node (CHN_E36, Chinese water collection and supply) carries nearly one third of total relative exposure, with the top fifteen nodes accounting for roughly 90%. The upper tail is dominated by Chinese sectors and water-utility nodes whose measured zero outgoing redundancy inflates structural resistance ($$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$3). The author concedes plainly that this concentration is a calibration feature of the WIOD accounting structure — large inflows, no recorded intermediate outflows — and should not be read as a sectoral vulnerability diagnosis. This is a substantive limitation of the exposure index, not merely a caveat.

Regime structure and phase diagrams

Baseline Monte Carlo experiments (15,000 periods per scenario, 100 replications) produce four ordered stationary regimes:

Regime	$$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$4	$$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$5	Mean $$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$6	$$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$7	$$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$8	$$H_{i,t} = \frac{B_t\, I_{i,t}}{D_{i,t}\, C_{i,t} + \varepsilon}$$9	$I_{i,t}$0
Stable absorption	0.45	0.7	0.084	0.084	0.000	0.000	0.000
Latent fragility	0.70	1.4	0.489	0.465	0.000	0.000	0.000
Critical transition	1.00	1.8	2.036	0.807	0.099	0.001	0.000
Avalanche	1.35	2.3	5.811	0.953	0.585	0.170	0.001

Mean avalanche size spans nearly two orders of magnitude, and each regime converges to a distinct stationary level — the dynamics is dissipative and non-explosive. The phase grid over $I_{i,t}$1 and $I_{i,t}$2 confirms Proposition 3: the absorption-to-avalanche frontier is positively sloped, with $I_{i,t}$3 activating first, $I_{i,t}$4 requiring stronger joint loading, and $I_{i,t}$5 never exceeding roughly 7% even at the most aggressive corner. The joint nature of the transition carries a direct diagnostic implication: single-indicator early-warning systems can miss the regime transition because the orthogonal axis has not yet moved.

Tail behaviour and the explicit rejection of SOC claims

The most disciplined part of the paper is its treatment of tails. Maximum-likelihood power-law estimation (Clauset-Shalizi-Newman methodology) yields $I_{i,t}$6 in the critical-transition regime and $I_{i,t}$7 in the avalanche regime ($I_{i,t}$8), with latent fragility returning a thin tail ($I_{i,t}$9) and stable absorption uninformative. These exponents are far steeper than the canonical SOC range of 1-2 (often near $D_{i,t}$0 in Bak-Tang-Wiesenfeld). The author therefore reports regime-dependent tail thickening but explicitly declines to claim self-organised criticality, aligning with the Stumpf-Porter critique of hasty power-law labelling.

Three factors are offered to interpret the steep exponents without treating them as falsification of SOC: (i) finite size, since $D_{i,t}$1 and maximum observed avalanches ($D_{i,t}$2) are close enough to system size for finite-size corrections to steepen the apparent slope; (ii) the subcritical dissipative operator ($D_{i,t}$3), which precludes the long-range correlations a true critical point requires; and (iii) driving at the same time scale as relaxation, rather than the infinite time-scale separation of canonical SOC. The honest conclusion is that $D_{i,t}$4 is a finite-size, finite-driving-rate diagnostic. Whether the steepening is a genuine finite-size artefact is left open, as the paper cannot vary $D_{i,t}$5 within WIOD.

Robustness checks strengthen the design: the four-regime ordering survives all three propagation operators, with tail exponents in the avalanche regime confined to $D_{i,t}$6, confirming that the results are not artefacts of row-stochastic normalisation. Sensitivity to $D_{i,t}$7, burn-in length and replication count is minimal.

Policy mechanics and limitations

The framework implies that policy leverage operates on the denominator of the stress-conversion function rather than on shocks themselves. Since $D_{i,t}$8 and $D_{i,t}$9 enter multiplicatively, simultaneous modest improvements in redundancy (alternative suppliers, modular standards, inventories) and capacity (liquidity buffers, insurance) yield disproportionate reductions in transversal stress. The extreme concentration of the exposure upper tail makes targeted interventions at the top of the distribution far more effective than median-node policies — a network analogue of granular intervention. The phase-frontier geometry further implies that resilience investment is most valuable anticipatorily, before $C_{i,t}$0 rises.

The limitations are stated with unusual candour. The calibration ends at WIOD 2014 and depends on harmonisation choices. The redundancy and capacity proxies use only network-level flow information, ignoring inventories, balance-sheet liquidity and alternative-supplier data. The external fields are controlled experiments, not measured shocks, so the results say nothing about the probability that any specific field materialises. The connection to historical crises (2008-09, the 2011 Tohoku earthquake, 2020, 2022) is explicitly a re-description, not an identification. And the tail exponents, being finite-size diagnostics on a 2,283-node substrate with rare large events, are statistically noisy.

Conclusion

The paper's contribution is a finite-size, real-network description of how Hall-like transversal stress activates structural fragility, deliberately framed as not a claim of self-organised criticality in the global economy. Three ingredients — a multiplicative stress-conversion function, dissipative sandpile dynamics, and a controlled field grid — suffice on the observed WIOD topology to generate an ordered four-regime structure with a positively sloped, normalisation-robust phase frontier. The open questions are specific: how $C_{i,t}$1 scales with network size and driving-rate ratio, whether measured shock proxies can populate the field $C_{i,t}$2 to enable historical identification, and how the mechanism integrates with curvature-based fragility measures. The paper is most valuable precisely where it is most restrained: it demonstrates that loading-induced regime transitions and tail thickening are robust features of a real dissipative production network, while refusing to over-interpret them as evidence of universal criticality.