The rapid growth of the large language model (LLM) ecosystem raises a critical question: are seemingly diverse models truly independent? Shared pretraining data, distillation, and alignment pipelines can induce hidden behavioral dependencies, or latent entanglement, that undermine multi-model systems such as LLM-as-a-judge pipelines and ensemble verification, which implicitly assume independent signals. In practice, this manifests as correlated reasoning patterns and synchronized failures, where apparent agreement reflects shared error modes rather than independent validation. To address this, we develop a statistical framework for auditing behavioral entanglement among black-box LLMs. Our approach introduces a multi-resolution hierarchy that characterizes the joint failure manifold through two information-theoretic metrics: (i) a Difficulty-Weighted Behavioral Entanglement Index (BEI), which amplifies synchronized failures on easy tasks, and (ii) a Cumulative Information Gain (CIG) metric, which captures directional alignment in erroneous responses. Through experiments on 18 LLMs from six model families, we identify statistically significant behavioral entanglement. Such behavioral dependence is associated with judge over-endorsement bias on a disjoint MMLU-Pro evaluation set (rho = 0.508 for BEI and rho = 0.520 for CIG; p < 0.01). The association further transfers to the MATH-500 benchmark (rho = 0.441 for BEI and rho = 0.457 for CIG; p < 0.05), providing cross-benchmark evidence that the identified dependency structure generalizes beyond the data and response format used for its estimation. Finally, we demonstrate a practical use case of entanglement through de-entangled verifier ensemble reweighting, achieving 3.5 and 2.6 percentage-point gains in accuracy and precision, respectively, over majority voting.
Figures & tables
Figure 1: A Multi-Resolution Entanglement Hierarchy for Behavioral Dependence
Figure 2: Behavioral entanglement graph of LLMs on MMLU-Pro: (a) co-failure entanglement (BEI); and (b) directional entanglement (CIG). Reported are only statistically significant BEI and CIG values and FDR corrected relationships (q-value <0.05 ).
Metric
MMLU-pro
MATH-500
ρ
p -value
ρ
p -value
Co-failure audit
Corr.
0.063
0.724
0.187
0.360
MI
0.113
0.525
0.208
0.309
Jaccard
0.428
0.012 ∗
0.275
0.175
BEI
0.508
0.002 ∗∗
0.441
0.024 ∗
Table 1: Metrics and their Spearman correlations with ΔPrec on MMLU-pro and MATH-500.
w
MMLU-pro ρ
MATH-500 ρ
0
0.289
0.032
1
0.508 ∗∗
0.441 ∗
2
0.490 ∗∗
0.275
3
0.504 ∗∗
0.231
4
0.465 ∗∗
0.236
5
0.363 ∗
0.180
Table 2: Sensitivity of BEI correlations to the difficulty-weighting exponent w .
Aggregation
Acc
F1
Precision
Δ Acc
Majority Vote
0.847
0.901
0.870
–
Accuracy-based Reweight
0.881
0.919
0.891
+0.034
Entangle-based Reweight
0.882
0.922
0.896
+0.035
Table 3: Verifier aggregation performance on MMLU-Pro.
Figure 3: Calibrated verifier ensemble weights.
Appendix figures & tables5 assets
Supplementary material from the paper’s appendix.
Appendix
Metric
Mean
Std.
Range
AUC
0.865
0.050
[0.752, 0.934]
ECE
0.041
0.012
[0.022, 0.077]
Brier
0.130
0.039
[0.065, 0.192]
Appendix
Table C.1: Calibration diagnostics for the difficulty-response models used in BEI.
Figure C.1: Reliability curve of the difficulty-response models used for BEI calibration.
Model 1
Model 2
BEI
p -value
q -value
Llama-3-70B
Llama-3_1-70B
0.0525
5.00E-05
9.50E-04
Llama-2-70b-hf
Qwen1.5-110B
0.0398
5.00E-05
9.50E-04
Llama-3-70B
Qwen1.5-110B
0.0352
5.00E-05
9.50E-04
Qwen1.5-14B-Chat
Qwen1.5-72B-Chat
0.0333
5.00E-05
9.50E-04
Llama-2-70b-hf
Llama-3-70B
0.0326
5.00E-05
9.50E-04
Llama-3_1-70B
Qwen1.5-110B
0.0325
5.00E-05
9.50E-04
Appendix
Table C.2: Top 10 pairwise BEI values with raw p -values and FDR-adjusted q -values.
Model 1
Model 2
CIG
p -value
q -value
Llama-3-70B
Llama-3_1-70B
0.1332
5.00E-05
1.09E-03
Claude-3-5-sonnet-20241022
GPT-4o-2024-08-06
0.0525
5.00E-05
1.09E-03
GPT-4o-2024-08-06
GPT-4o-mini
0.0523
5.00E-05
1.09E-03
DeepSeek-chat-v2_5
Gemini-1.5-flash-002
0.0502
5.00E-05
1.09E-03
Claude-4-6-Sonnet
Gpt-5
0.0471
5.00E-05
1.09E-03
Gemini-1.5-flash-002
Gemini-1.5-pro-002
0.0454
5.00E-05
1.09E-03
Appendix
Table C.3: Top pairwise CIG values with raw p -values and FDR-adjusted q -values.
Target
Test
Null hypothesis
p -value
Interpretation
ΔPrec
Paired t -test
E[pt(V2)−pt(V1)]=0
0.421
No mean shift
Ranking
Equivalence test
ρ≤1−ϵ,ϵ=0.05
7.6×10−5
Ranking consistent
Appendix
Table C.4: Joint-prompt versus independent-prompt evaluation ablation.