Distinct dynamics of conceptual and referential disruptions in human reading and large language model processing
Authors: Rui He, Nihal Altay, Wolfram Hinzen
Organizations: Grammar and Cognition Lab, Department of Translation & Language Sciences, Universitat Pompeu Fabra, Barcelona, Spain · Institut Català de Recerca i Estudis Avançats (ICREA), Barcelona, Spain
Linguistic meaning is grounded in conceptual content, from which reference to particular entities emerges as words enter discourse. To examine the processing dynamics associated with these two dimensions of meaning, we selectively disrupted conceptual or referential information in short narratives and traced the resulting effects in human self-paced reading and in the predictive and representational processing of large language models. In human reading, conceptual disruptions produced a strong but localized processing cost, emerging immediately after the distorted word, reaching an early maximum, and then declining rapidly. Referential disruptions produced weaker effects, which decreased more gradually across subsequent words, and were more strongly modulated by sentence boundaries. In the language model, both disruptions emerged immediately at the manipulated word. Contextual model surprisal showed a pattern closely paralleling human reading: conceptual disruption produced a larger, more locally concentrated effect that decayed rapidly, whereas referential disruption produced a smaller and more gradual downstream effect. Output-layer representations showed a different pattern: referential disruption produced a larger initial displacement, while both distortions were subsequently characterized by power-law decay. Together, these results provide convergent evidence for distinguishable processing dynamics of two types of meaning: conceptual information imposes a more locally concentrated integration cost, whereas referential information engages a more distributed process of maintaining discourse-level identity.
Figures & tables
Figure 1: Comparability of stimulus versions on whole-story and lexical characteristics. Distributions of (A) whole-story perplexity, (B) mean lexical frequency (Zipf scale), and (C) mean word length across Original, conceptually distorted (CD), and referentially distorted (RD) versions of the 21 stories. Points represent individual stories, with grey connecting lines linking the three matched versions of each story. Pairwise differences were evaluated using Siegel–Castellan comparisons following the corresponding Friedman test. Significance labels indicate FDR-adjusted p values: n.s., p≥0.05 ; *, p<0.05 ; **, p<0.01 ; ***, p<0.001 .
Figure 2: Local effects of conceptual and referential distortions in human reading and LLM processing. (A) Estimated effects of conceptual distortion (CD) and referential distortion (RD), relative to the matched Original condition, on residualized log reading time. (B) Changes in contextual surprisal produced by CD and RD in Qwen3-4B-Base, calculated relative to the corresponding Original words. (C) Output-layer representation distance between distorted and Original versions, quantified as cosine distance between the corresponding final-layer hidden states. K denotes the final distorted token; K-3 to K-1 and K+1 to K+3 denote positions relative to K, and K_END denotes the sentence-final word. Points show estimated effects and error bars indicate 95% confidence intervals. The horizontal dashed line indicates zero effect. Asterisks indicate effects that remained significant after FDR correction across analysis positions, applied separately for each distortion type within each outcome (FDR-adjusted p<0.05 ).
Figure 3: Downstream reading-time dynamics following conceptual and referential distortions. (A–B) Distance-specific effects of conceptual distortion (CD) and referential distortion (RD), relative to matched Original observations, on residualized log reading time from K+1 to K+10. Points show mixed-effects-model estimates with 95% confidence intervals, and solid lines show the best-fitting functional forms selected by AICc. CD was best described by a log-normal trajectory, whereas RD was best described by a linear trajectory. The trajectories differed significantly across distance, χ2(9)=66.58,p<0.001 . (C–D) AICc values for the 11 candidate functional forms fitted to the CD and RD trajectories. Lower values indicate better model support, and stars mark the best-supported model. (E) Additional distortion effect associated with reaching the sentence boundary. Boundary modulation differed significantly between CD and RD ( p=0.035 ). Displayed p values above each estimate indicate the within-distortion boundary effects. (F) Model-estimated distortion effects averaged over an early downstream window (K+1 to K+3) and a late window (K+8 to K+10). The distance-window × distortion-type interaction was significant ( p=0.042 ), indicating that the early-to-late change differed between CD and RD. Error bars indicate 95% CI. Dashed lines indicate zero effect.
Figure 4: Downstream contextual surprisal dynamics following conceptual and referential distortions. (A–B) Distance-specific effects of conceptual distortion (CD) and referential distortion (RD), relative to matched Original observations, on contextual surprisal from K to K+10. Points show clustered-OLS estimates with 95% confidence intervals, and solid lines show the best-fitting functional forms selected by AICc. CD was best described by a log-Cauchy trajectory, whereas RD was best described by an exponential trajectory. The trajectories differed significantly across distance, F(10,122)=4.73,p<0.001 . (C–D) AICc values for the 11 candidate functional forms fitted to the CD and RD trajectories. Lower values indicate better model support, and stars mark the best-supported model. (E) Additional surprisal effect associated with reaching the sentence boundary. Boundary modulation differed significantly between CD and RD ( p=0.019 ). Displayed p values above each estimate indicate the within-distortion boundary effects. (F) Model-estimated surprisal effects averaged over an early downstream window (K+1 to K+3) and a late window (K+8–K+10). The distance-window × distortion-type interaction was significant ( p=0.017 ), indicating that the early-to-late change differed between CD and RD. Error bars indicate 95% CI. Dashed lines indicate zero effect.
Figure 5: Downstream representational dynamics following conceptual and referential distortions. (A–B) Distance-specific effects of conceptual distortion (CD) and referential distortion (RD), relative to matched Original representations, on output-layer representation distance from K to K+10. Representation distance was quantified as cosine distance between the corresponding final-layer hidden states. Points show clustered-OLS estimates with 95% confidence intervals, and solid lines show the best-fitting functional forms selected by AICc. Both CD and RD were best described by power-law trajectories. The trajectories differed significantly across distance, F(10,122)=7.70,p<0.001 . (C–D) AICc values for the 11 candidate functional forms fitted to the CD and RD trajectories. Lower values indicate better model support, and stars mark the best-supported model. (E) Additional representation-distance effect associated with reaching the sentence boundary. Boundary modulation did not differ significantly between CD and RD ( p=0.142 ). Displayed p values above each estimate indicate the within-distortion boundary effects. (F) Model-estimated representation distances averaged over an early downstream window (K+1 to K+3) and a late window (K+8 to K+10). The distance-window × distortion-type interaction was not significant ( p=0.443 ), indicating no evidence that the early-to-late change differed between CD and RD. Error bars indicate 95% CI. Dashed lines indicate zero effect.
Appendix figures & tables3 assets
Supplementary material from the paper’s appendix.
Appendix
Figure S1: Local effects of surprisal and representational distance with Llama 3.2 3B.
Figure S2: Effect decay of Llama 3.2 surprisal.
Figure S3: Effect decay of Llama 3.2 representational distance.
Language models (LMs) are often used as a tool to model human language processing. Recent studies suggest that severely restricting LMs' context window improves their fit to human psycholinguistic data by simulating human working memory constraints. However, it is possible that this strict memory-decay approach overlooks humans' reliance on long-range structural representations, such as discourse structre. In this work, we systematically vary the context window size of GPT-2 across four large-scale naturalistic English reading-time datasets and observe a U-shaped relationship: Although restricted contexts (< 20 tokens) successfully capture local memory limitations, expanded contexts (500--1,000 tokens) ultimately yield the highest overall psycholinguistic fit. To investigate the mechanism driving this benefit, we conduct a counterfactual inference-time experiment that disrupts cross-sentential entity chains by pronominalizing repeated discourse entities. Obscuring these structural linkages significantly degrades the predictive power of larger context windows by 20% to 40%. Our experiments demonstrate that tracking long-range coreference relations is one important factor for the alignment between LM surprisal and human reading behavior, and approximate the extent to which human comprehenders use global discourse relations during language processing.
Kohei Kajikawa, Lin Ai, Tatsuki Kuribayashi +1
Department of Linguistics, Georgetown University, USA · Computing and Mathematical Sciences Division, MBZUAI, UAE · Center for Language AI Research, Tohoku University, Japan
A recent study (Kuribayashi et al., 2025) has shown that human sentence processing behavior, typically measured on syntactically unchallenging constructions, can be effectively modeled using surprisal from early layers of large language models (LLMs). This raises the question of whether such advantages of internal layers extend to more syntactically challenging constructions, where surprisal has been reported to underestimate human cognitive effort. In this paper, we begin by exploring internal layers that better estimate human cognitive effort observed in syntactic ambiguity processing in English. Our experiments show that, in contrast to naturalistic reading, later layers better estimate such a cognitive effort, but still underestimate the human data. This dual alignment sheds light on different modes of sentence processing in humans and LMs: naturalistic reading employs a somewhat weak prediction akin to earlier layers of LMs, while syntactically challenging processing requires more fully-contextualized representations, better modeled by later layers of LMs. Motivated by these findings, we also explore several probability-update measures using shallow and deep layers of LMs, showing a complementary advantage to single-layer's surprisal in reading time modeling.
Tatsuki Kuribayashi, Alex Warstadt, Yohei Oseki +1
Human language comprehension unfolds sequentially: each word is processed in the context of those that came before, and the interpretation builds incrementally over time. Surprisal, the negative log probability of a word given its context, has been the dominant predictor of incremental processing cost. But surprisal reduces rich sequential representations to a single scalar at each word, discarding information about the direction in which the interpretation has been evolving. Dynamical-systems approaches suggest that the trajectory of the evolving interpretive state, not just its position at each moment,should shape processing, and language itself may have local momentum, since speakers plan utterances a few words at a time. We introduce trajectory extrapolation error: at each word, we fit a linear trajectory to the preceding hidden states of a transformer language model and measure deviation from the extrapolated path. On the Natural Stories corpus, this measure is nearly orthogonal to surprisal (r = .044) and independently predicts self-paced reading times. The effect is especially pronounced in garden-path sentences, strengthens with model scale (GPT-2 Small to Large), and replicates across architectures with different positional encoding schemes (GPT-2 vs. Pythia/RoPE). A displacement control shows the effect is not reducible to representational change magnitude: displacement and extrapolation error predict in opposite directions. These findings reveal two dissociable components of processing cost: word-level prediction error (surprisal) and sensitivity to the local momentum of the unfolding interpretation (trajectory extrapolation error).
Elan Barenholtz
Machine Perception & Cognitive Robotics Laboratory Department of Psychology / Center for Complex Systems Florida Atlantic University