Automated essay scoring (AES) with neural network technology, such as large language models (LLMs), is of interest in language assessment for minimizing hallucinations and increasing score validity. English language testing companies either benchmark by agreement with human raters or implement LLMs in a decomposed architecture to minimize the impact of the “black box phenomenon.” However, high-stakes testing needs scoring that can be controlled and verified. We investigate whether scoring rubrics are represented inside an LLM-based scorer and whether they can be controlled. We probe and steer 25 open-weight base LLMs (0.13B to 10.31B parameters, 14 families) on human-scored BESTEP English writing responses across three task types and six scoring rubrics. We find holistic scores are generally decodable from the models’ activations. Furthermore, scoring rubrics occupy partially distinguishable directions in representation space for longer answers (Parts 2 and 3), but not the shortest (Part 1). Finally, steering along a rubric direction changes scores for five of six scoring rubrics. Among the four reliably measured rubrics, applying the matching intervention (e.g., steering toward ‘better grammar’ on a grammar-degraded essay) offered no detectable advantage over applying a mismatched one (e.g., steering toward ‘better vocabulary’) in 11 of 12 head-to-head comparisons. The models behave as if every scoring rubric shares a single control over scores. This may suggest the models collapse the rubrics internally, or they reproduce the halo effects observed in human raters. Our contribution is showing whether scoring rubrics can be controlled independently, which reveals evidence that score agreement alone cannot provide.
Keywords: Automated Essay Scoring, Large Language Models, Rubric-Based Assessment, Interpretability, Construct Validity
Automated essay scoring (AES) has recently started implementing neural network models, e.g., large language models (LLMs), in the grading process (Johnsi and Kumar, 2026). Agreement with human raters rose along the way, but the models remain opaque to inspection. Transparency and interpretability remain open priorities in AES investigations (Johnsi and Kumar, 2026). Control and verifiability are hard to explain in LLM-based AES systems, which is one of the biggest challenges when implementing LLM technologies in high-stakes language assessments. Test providers respond in different ways. Some rely on mixed models where human raters have the final score. Some develop feature-based machine learning models in which humans define what is measured. Others, such as the Duolingo English Test (DET), decompose the writing construct into subconstructs and constrain any LLM-derived feature score to a single subconstruct, so that the model cannot influence the holistic score except through that one slot. As stated in the report, scoring directly with deep learning “cedes both control and verifiability” (Madnani et al., 2026).1
Although AES designs such as the DET minimize the black box impact on the score, the core question remains unsolved: what exactly an LLM scorer does on the inside. As LLMs become more and more ubiquitous in language assessment, we need to know whether a model’s internal representations contain the scoring rubrics at all, and whether each scoring rubric is controllable and verifiable separately. Current score agreement benchmark tests cannot answer this. Prior interpretability work on AES examines what a scorer has learned or which features it weights (Fiacco et al., 2023; Wang et al., 2026). Given that, we approach the topic with two AES principles in mind, control and verifiability, by adding causal intervention on the representations themselves, mirroring how these systems operate in real-world assessments.
In this study, we audit LLM-based AES from the inside by probing and steering 25 open-weight base LLMs (0.13B to 10.31B parameters, 14 families), while they score the BESTEP (Taiwan’s Ministry-of-Education English exam) responses across three task types and six analytic scoring rubrics. We carry out the investigation guided by three questions. RQ1: Is the holistic score decodable from a model’s activations, and where does that signal live? RQ2: Do the six scoring rubrics occupy distinguishable directions in representation space? RQ3: Does steering a rubric’s direction move the score specifically for that scoring rubric, or does any rubric-adjacent push move it about equally? In this paper, we choose base models as a baseline, so that fine-tuned variants (instruction-tuned, chat) can be compared against them in future studies. Our contributions are a white-box audit protocol for LLM-based AES that separates a rubric’s decodability, separability, and causal control. Moreover, the investigation provides evidence that, in these base models, decodability does not imply controllability. The code and the database behind every number are public for verification at https://github.com/lopentu/inside-the-scorer-paclic2026.
AES and LLM-based scoring. Automated essay scoring (surveyed in Johnsi and Kumar (2026)) first moved from feature-engineered methods to BERT-based deep learning: domain-adaptive training, multi-scale essay representations, and hierarchical transfer learning each closed part of the gap to human-rater agreement (Ramesh and Sanampudi, 2022; Wang et al., 2022; Cao et al., 2020; Xue et al., 2021; Settles et al., 2020). A second wave prompts general-purpose LLMs directly as scorers or raters, with a mixed record. Some studies report LLM ratings holding up under real psychometric scrutiny, generalizability-theory and many-facet Rasch analyses across thousands of essays and multiple genres (Song et al., 2025; Wang et al., 2025a), and case studies where an LLM produces both a score and trait-level comments across several analytic scoring rubrics at once (Wang et al., 2025b; Cho et al., 2024). Others find the opposite: robustness failures under adversarial stress-testing and demographic bias in ratings (Kabra et al., 2022; Yamashita, 2025), systematic disagreement between LLM and human grading behavior (Mathew et al., 2026), and validity that shifts by language, genre, and proficiency level rather than holding uniformly (Mizumoto and Eguchi, 2023; Pack et al., 2024; Mansour et al., 2024; Gaggioli et al., 2025; Barbosa et al., 2025; Yavuz et al., 2025; Yancey et al., 2023). A related strand looks past scoring accuracy to classroom use, from ChatGPT’s promise and pitfalls as a writing tutor to its measurable effect on students’ own writing when used for feedback (Barrot, 2023; Ding et al., 2025).
Control and verifiability in operational scoring. High-stakes providers mostly manage opacity through design. The Duolingo English Test partitions features by subconstruct, so a deep-learning feature, including a fine-tuned LLM rating, reaches the score only through its own subconstruct, and the contribution of each subconstruct can be measured and explained (Madnani et al., 2026). The same report notes that explaining why an end-to-end neural scorer assigned a score remains substantially harder. Prior analyses of neural AES scorers are observational: Wang et al. (2026) analyze how LLM scorers weight textual features, and Fiacco et al. (2023) align components of a transformer scorer’s activations with human-understandable feature groups. We intervene causally on the representations.
Rubric-aware and rationale-based scoring. A separate line augments scores with rubric-level explanations or rationales (Bueno et al., 2026; Chu et al., 2025; Jiao et al., 2025; Kucia et al., 2026). These works evaluate rubric outputs; we evaluate rubric representations.
Mechanistic interpretability and faithfulness. Our probing and steering methodology builds on activation-patching practice (Zhang and Nanda, 2023), causal concept-based explanation (Wu et al., 2023b), and work questioning whether model-generated explanations are faithful to internal computation (Atanasova et al., 2023; Jacovi and Goldberg, 2021; Chuang et al., 2026; Ferrando et al., 2023; Demircan et al., 2024), a concern that applies to AES, where a plausible-sounding rationale is not evidence of rubric-specific control. The probing literature motivates our two checks. Elazar et al. (2021) show that a property’s decodability from a representation does not mean the model actually uses that property, only that it is present, and Hewitt and Liang (2019) show that a probe’s accuracy is uninterpretable without a selectivity check against a control task the probe could only solve by memorization rather than by reading the representation. RQ1’s decodability result and RQ3’s steering-selectivity test correspond to these two checks, applied to scoring rubrics rather than linguistic properties.
Rubric-adjacent NLP. For the grammar and vocabulary scoring rubrics specifically, we draw on grammatical error correction (Bryant et al., 2023; Wu et al., 2023a; Fang et al., 2023) and lexical and proficiency assessment (Aleksandrova and Pouliot, 2023; Bannò et al., 2024; Bannò et al., 2025) literatures, and for inter-rater reliability, the standard kappa statistic (Cohen, 1960).
Models. We probe and steer 25 open-weight, base (non-instruction-tuned) LLMs spanning 0.13B to 10.31B parameters and 16 to 42 layers across 14 model families (Llama, Gemma, Falcon3, SmolLM2/3, Granite, Mistral, OLMo, TAIDE, Apertus, and others; see Appendix C, Table 6 for the full catalog), selected for complete output coverage across every phase of the pipeline (probing, geometry, and steering). All 25 are base models, by design. Operational scorers are often fine-tuned (the Duolingo English Test, for instance, uses LLM features fine-tuned on human ratings (Madnani et al., 2026)), but a base model carries no instruction or format tuning that could confound the score readout, so it gives a clean baseline against which instruction-tuned, chat, and other fine-tuned variants can be compared in future work. The parameter range also reflects a deployment constraint raised in the language-testing literature: local institutions weighing AI-assisted assessment need models practical to host on-premises, keeping response data local, rather than only the largest cloud-hosted systems (Chuang and Yan, 2025).
Dataset. All models are probed and scored on the same fixed set of 75 human-verified BESTEP writing responses (the gold-75), spanning two test forms and three task parts (Table 1): Part 1 (Answering Questions, a short factual-response task), Part 2 (Expressing Opinions), and Part 3 (Integrated Writing, the longest and most demanding task). BESTEP assesses six analytic scoring rubrics, but not every scoring rubric applies to every part (Table 3): grammar and vocabulary apply throughout, spelling and punctuation is scored only in the short-answer part (Part 1), and argument relevance only in the integrated-essay part (Part 3).
| Form | Part 1 | Part 2 | Part 3 | Total |
|---|---|---|---|---|
| Form 1 | 11 | 12 | 15 | 38 |
| Form 2 | 10 | 12 | 15 | 37 |
| Total | 21 | 24 | 30 | 75 |
Activation extraction. For every model, we extract final-token residual-stream activations at every layer for each of the 75 responses; these activations are the shared input to the probing (RQ1), contrastive-direction (RQ2), and steering (RQ3) procedures described in the next three sections.
Code and Data Availability. Code, the analysis scripts, and the database slice behind every reported number are available at https://github.com/lopentu/inside-the-scorer-paclic2026. They are released for noncommercial research use under the PolyForm Noncommercial License 1.0.0. The degraded essay pairs (Appendix A) contain exam responses and are not yet public; they are available from the authors for research purposes on request.
Method. For each (model, layer, part) cell we fit a Ridge regression probe on final-token residual-stream activations to predict the response’s holistic band score, using leave-one-out cross-validation (LOOCV; Figure 1): each of the 21 (Part 1), 24 (Part 2), or 30 (Part 3) responses is held out once, the probe is trained on the rest, and every held-out prediction is pooled to compute one R2 per cell. We report the peak-R2 layer, the peak R2 itself, and the peak-to-final-layer drop, using the pooled-form scope (both test forms combined) as the headline result. To check that a peak R2 exceeds what the probing procedure produces on its own, each cell also gets a label-shuffle permutation test (1,000 shuffles of the scores). Because the peak is the best of 16 to 42 layers, we report the peak layer’s p-value both raw and Bonferroni-corrected over that model’s layers.
Results. Score information is decodable in most models, but not equally for every part, and not in a simple “easy parts early, hard parts late” order. Against the permutation null, the peak layer beats chance in all 25 models for Part 3 and in 22 of 25 for Part 2 after the layer correction, but in only 9 of 25 for Part 1 (24 of 25 uncorrected). The one consistent timing pattern: Part 1 (the shortest, short-answer task) peaks latest or second-latest among the three parts in 23 of 25 models (Figure 2). A response-length explanation for this timing does not hold up: relative peak-to-final signal drop ranges from 0.0% to 2611.3% across models, and the part with the largest drop varies by model, so timing alone cannot be read as a complexity gradient.
Peak timing and decodability strength diverge. Part 1 has the lowest peak R2 in 24 of 25 models and is never the highest; across models, mean peak R2 is 0.432 for Part 1, 0.633 for Part 2, and 0.773 for Part 3 (Table 2, Figure 3). Part 1 is also the least stable across models (peak-R2 variance 0.0338, versus 0.0095 for Part 2 and 0.0038 for Part 3). Part 1 has the fewest gold responses (21), the shortest responses (about 47 words on average, versus about 117 for Part 2 and about 214 for Part 3), and the coarsest scale, which plausibly explains why its signal is both weaker and more model-dependent: less textual evidence and fewer observations for the probe to exploit.
| Part | Mean peak R2 | Median peak R2 | Peak-R2 variance | Min | Max |
|---|---|---|---|---|---|
| Part 1 | 0.432 | 0.485 | 0.0338 | 0.007 | 0.702 |
| Part 2 | 0.633 | 0.631 | 0.0095 | 0.413 | 0.877 |
| Part 3 | 0.773 | 0.773 | 0.0038 | 0.605 | 0.858 |
For the two longer tasks, models carry information about essay quality that a label-shuffle null cannot explain. For BESTEP’s shortest task type the signal peaks later, is weaker and less stable, and is above the corrected null in only 9 of 25 models. We cannot explain this pattern mechanistically here, but RQ2 and RQ3 return to it. The Bonferroni correction is conservative because adjacent layers are correlated, and 1,000 shuffles bound how small a p-value can get (Limitations).
Method. RQ2 uses the same 75 gold-standard responses plus rubric-targeted clean-versus-degraded mutations, since holistic scores blend scoring rubrics and cannot support construct-geometry claims on their own. Mutations are generated by prompting llama3.3:70b (via Ollama) with one of 12 rubric-quoting specs, one per (part, scoring rubric) combination the rubric scores (Table 3), each spec quoting the official BESTEP band descriptors for that scoring rubric and instructing the model to degrade the essay to the low band on that dimension only, holding every other dimension, and meaning, constant (length preserved to within 10%). Two trained BESTEP graders verified each mutated essay against its intended target dimension before use (Appendix A gives the procedure and worked examples). For each (model, layer, part, scoring rubric) cell we estimate a contrastive direction as the normalized difference between clean and mutated-essay activation centroids, and its split-half reliability by recomputing the direction from random halves of the essay pairs. We compare same-rubric to different-rubric cosine similarity within the same essay part (a same-part, cross-rubric comparison avoids confounding “different construct” with “different essay”), retaining only reliability-gated rows (split-half reliability at least 0.6).
| Scoring rubric | Part 1 | Part 2 | Part 3 |
|---|---|---|---|
| Grammar | v | v | v |
| Vocabulary | v | v | v |
| Cohesion | v | v | |
| Spelling/Punct. | v | ||
| Task Relevance | v | v | |
| Argument Relevance | v |
Results.
Separation is uneven across the assessment. Averaged across 1,470 reliability-gated layer-part cells, the mean same-versus-different geometry gap is 0.334 (74.5% of cells positive), but this is not stable across task parts: Part 1’s four scoring rubrics are, on average, negative (mean gap -0.220, 6.2% of rows positive), no more separated from each other than a single scoring rubric is from itself across parts, while Parts 2 and 3 separate cleanly (89.7% and 99.2% of cells positive): reliable, rubric-specific directions are broadly available in Parts 2 and 3 but sparse in Part 1.
Pair by pair, the range is wide. Pooling all reliability-gated rows
by literal rubric pair, using disattenuated cosine similarity (a
correction for each direction’s own measurement noise, so a value of 1
means statistically indistinguishable from the same construct and 0
means unrelated), the 12 comparable pairs range from anti-correlated to
statistically identical (three of the 15 possible pairs never co-occur
in the same part and are not comparable at all; Table 4, alongside RQ3’s
steering results below). argument_relevance paired with
grammar is the most separated pair (-0.102, anti-correlated
beyond noise), and argument_relevance, a Part-3-only
scoring rubric, carries most of the separation signal in the whole
rubric. At the opposite end, grammar paired with
spelling_punctuation sits at the 1.0 ceiling: the two are
statistically indistinguishable from the same construct, and
spelling_punctuation paired with vocabulary is
nearly as entangled (0.984). Both of the most-entangled pairs are
Part-1-only comparisons, consistent with Part 1’s negative geometry gap
above.
So separation is broad in the two longer, more integrated task types, concentrated in one scoring rubric (argument relevance) beyond that, and absent in the shortest task type.
Method. RQ3 reuses RQ2’s rubric directions and tests them causally rather than only representationally. A direction is scaled to unit length and added to the residual stream at one layer, at every token position of the grading prompt, with magnitude alpha times that layer’s mean residual-stream norm (alpha = 0.05, 0.1, 0.25, 0.5), so strength is comparable across layers and models. The outcome is the change in expected band score: the probability-weighted mean over the digit tokens at the answer position of the same grading prompt (no learned readout), minus the unsteered baseline. For degraded essays with a known target scoring rubric we inject four kinds of direction: the matching rubric direction (self), a mismatched rubric direction (other), an orthogonal control (a rubric direction with the axis shared by all scoring rubrics in the part removed, which tests whether overlap between directions hides specificity), and a generic holistic direction (the direction the RQ1 probe uses to predict the band score, which tests whether any score-raising push works). Self and other effects are paired at the finest valid unit, (model, layer, part, form, degradation, alpha), giving 74,304 paired cells across all 25 models, and we test whether the matching direction beats a mismatched one by a practically meaningful margin.
Results.
Steering works for five of six scoring rubrics. Grammar, vocabulary,
cohesion, spelling/punctuation, and task relevance each show a clearly
positive self effect (0.168 to 0.359 score points). Only
argument_relevance’s own direction produces no detectable
effect: it neither raises nor lowers the score.
Steering is not selective. Pooled across all paired cells, the matching direction beats a mismatched one by only 0.0074 score points on average (95% CI [0.0036, 0.0117]; 1.74% of the absolute self effect, 95% CI [0.85%, 2.74%]), and this holds under every robustness check we ran: the gap holds at the cell level (0.0085, 95% CI [0.0051, 0.0119]), never exceeds 0.0083 under leave-one-model/family-out, and a cluster bootstrap that resamples whole model families rather than individual models, accounting for shared training lineage such as the three Llama-3.x checkpoints, gives an almost identical estimate (0.0073, 95% CI [0.0035, 0.0108]).
The pooled number combines two different patterns. Broken out pair by
pair (Figure 4), among the four scoring rubrics that both work and are
reliably measured (grammar, vocabulary, cohesion, spelling/punctuation),
the matching direction shows no detectable advantage over a mismatched
one in 11 of their 12 head-to-head comparisons; the one exception is
vocabulary losing to spelling/punctuation, not because vocabulary is
weak, but because spelling/punctuation happens to be the single
strongest general-purpose lever of the six.
argument_relevance loses every comparison it has data for,
for the different reason that its own direction never moves the score.
Re-run separately within each task part (36 within-part comparisons
among comparably-strong scoring rubrics), the same pattern holds in
every part: not one instance of a scoring rubric beating a comparably
strong rival specifically because it matched the essay’s
degradation.
| Rubric pair | n | Sim. | 95% CI |
|---|---|---|---|
| Arg. Rel. x Gram. | 564 | -0.10 | [-0.13, -0.08] |
| Spell. x Task Rel. | 6 | 0.01 | [-0.36, 0.37] |
| Arg. Rel. x Cohes. | 586 | 0.08 | [0.06, 0.11] |
| Gram. x Task Rel. | 416 | 0.28 | [0.25, 0.32] |
| Arg. Rel. x Vocab. | 419 | 0.31 | [0.29, 0.34] |
| Cohes. x Gram. | 800 | 0.39 | [0.35, 0.42] |
| Gram. x Vocab. | 823 | 0.49 | [0.46, 0.52] |
| Task Rel. x Vocab. | 336 | 0.51 | [0.48, 0.55] |
| Cohes. x Task Rel. | 334 | 0.52 | [0.48, 0.55] |
| Cohes. x Vocab. | 691 | 0.61 | [0.58, 0.65] |
| Spell. x Vocab. | 192 | 0.98 | [0.96, 1.00] |
| Gram. x Spell. | 328 | 1.00 | [1.00, 1.00] |
Two controls calibrate these effects (Table 5). The orthogonal control does not behave like a flat baseline: it lowers the score slightly and more so at larger alpha, a small movement with the opposite sign of any working self effect, so the model is not indifferent to an injected direction, only to which scoring rubric it targets. The generic holistic direction raises the score less than self or other at every alpha tested, which argues against a pure “generic tampering” account of the specificity gap above; its alpha slope falls inside the self/other range, so it behaves like a weaker member of the same family of levers. At the strongest alpha (0.5) the average self effect is 0.361 score points on the 1-6 band scale, about 6% of the band range, a moderate intervention.
| Metric | Orthogonal | Band probe |
|---|---|---|
| Delta | -0.097 | 0.130 |
| Delta 95% CI | [-0.192, -0.020] | [0.024, 0.225] |
| Slope | -0.288 | 0.550 |
| Slope 95% CI | [-0.519, -0.063] | [0.251, 0.846] |
Tying back to RQ2, a pair’s representational separation (Table 4, above) does predict how large a causal difference that pair shows under steering (r=0.717 across all comparable pairs), but it does not predict which direction wins.
RQ2’s directions are reliability-gated, but RQ3 injects the raw, non-disattenuated direction regardless of noise, and different scoring rubrics’s directions are already substantially correlated (Table 4). Noise and this overlap could be masking a larger true specificity, and three checks address that. The first restricts to progressively more reliable directions, which raises the specificity ratio from 1.8% (all cells) to 5.4% (RQ2 geometry available) to 3.1% at reliability at least 0.7 to 8.3% at the top reliability quartile (self 0.670 vs. other 0.605, gap 95% CI [0.057, 0.072]): less reliable directions do understate true specificity, but even the best-measured quarter shows self only modestly outperforming other. A second check removes the axis shared by all scoring rubrics in a part from each direction (the orthogonal-control condition above) and repeats self vs. other on these overlap-corrected directions, using data already collected. If overlap were concealing real independent control, this should reveal it, and it does not. Both self_orth and other_orth average slightly negative (-0.093, -0.100), and their gap (0.0064, 95% CI [0.0003, 0.0125], ratio 4.82%) matches the reliability-gated result: once the shared component is removed, the rubric-specific residual is not a strong positive lever for either scoring rubric, so most of the raw effect was the shared component, not overlap hiding independent control.
A third check isolates what a direction shares only with its specific rival. For each pair of scoring rubrics we strip from each direction its projection onto the other, so that the matching and mismatched directions lose exactly their mutual overlap (pairs too collinear to leave a residual are skipped), and inject the two residuals next to a norm-matched random direction as a noise floor. Across all 25 models (884 pair cells), both residuals raise the score above the noise floor, by 0.194 (matching) and 0.181 (mismatched) points, so rubric-specific residual structure exists and moves the score. The matching residual still does not win: the gap is 0.013 points (95% CI [-0.015, 0.042]), about 4.5% of the average absolute matching effect, and the interval includes zero. Pairs below and above a cosine of 0.6 give the same picture (gaps of 0.012 and 0.014). Removing the overlap pair by pair leaves the levers interchangeable.
Taken together, the three research questions give one account of what a rubric-based LLM scorer does inside. The model carries decodable information about essay quality (RQ1), strongly for the two longer tasks and weakly for the shortest. It represents at least some scoring rubrics as distinguishable directions rather than one undifferentiated quality signal (RQ2). But distinguishable representation does not translate into independently controllable behavior (RQ3). Steering a rubric direction reliably moves the score, yet pushing the matching direction is, for the scoring rubrics that are both effective and reliably measured, essentially interchangeable with pushing a different one. The model appears to exercise one shared control with several equally effective handles, not six independent, rubric-specific ones.
This matters for how neural AES is justified. Operational providers that need control and verifiability get it by construction: in a decomposed design each subconstruct has its own slot, and a neural feature can influence the score only through its slot (Madnani et al., 2026). An end-to-end LLM scorer offers no such slot, and our audit asks whether it contains something equivalent on its own. In these base models it does not: scoring rubric information is present and partly separable, but it is not exposed as separate controls. High agreement with human raters, or the mere decodability of scoring rubrics, is therefore not evidence that a model uses those scoring rubrics as separately controllable scoring levers. Controllability is a diagnostic that score agreement cannot supply.
What the diagnostic means is open, because shared control has two readings. One is that the model collapses scoring rubrics that the test keeps apart. The other is that the model faithfully reproduces the structure of human analytic ratings, in which dimensions are known to move together: halo, the carry-over of a rater’s overall impression into separate trait scores, is a documented rater effect in the scoring of written composition (Engelhard Jr, 1994; Myford and Wolfe, 2004), and its relation to rating accuracy is itself contested (Fisicaro, 1988). Our gold set carries only holistic band scores, so we cannot separate the two readings here. Either way, a provider cannot assume that rubric labels from a neural scorer are separately controllable without checking.
The evidence for separation and for nonselective steering also comes
from partly different scoring rubrics, and the two should be read
together. argument_relevance carries most of the separation
signal in RQ2 but its steering effect is undetectable, whereas the four
scoring rubrics that steer and are reliably measured have moderately to
highly similar directions (Table 4). Nonselective steering among
overlapping directions is what overlap alone would predict, which is why
the pairwise check in RQ3 matters: with each pair’s overlap removed,
both residuals still move the score and still do so about equally.
The pattern is not uniform across the assessment design. Part 1, BESTEP’s shortest task, shows the weakest and least stable score signal (RQ1) and no meaningful rubric separation (RQ2), while the two longer, more integrated parts support a stronger signal and real, if not independently controllable, rubric structure. Validity claims should be checked at that grain, not assumed to generalize across task format.
All 25 models are base (non-instruction-tuned). We chose them as a
baseline, but it means we cannot say how instruction-tuned, chat, or
task-fine-tuned variants, which resemble deployed scorers more closely,
would behave. The gold set is 75 essays from one assessment (BESTEP),
which bounds essay-level diversity. Model coverage tops out at 10.31B
parameters, so results may not hold for substantially larger models.
argument_relevance’s own steering direction never reaches
significance; whether that reflects a representational gap or a scoring
rubric that needs more than surface language features to steer is open.
More models, other writing assessments, and non-English rubric-based
assessments could test whether the findings travel.
Two limits concern the degraded essays that define every rubric direction. They were generated by one LLM, and two trained BESTEP graders then checked every one of the 104: each, separately, saw an original and its degraded version, was told the essay had been degraded on one scoring rubric, and had to name which, and both confirmed the intended scoring rubric every time (Appendix A). That shows the target scoring rubric is recognizable, not that nothing else changed. The grammar/spelling_punctuation similarity ceiling in Table 4, and the near-ceiling spelling/vocabulary pair, can arise because errors in one dimension leave traces in another even in “pure” stimuli (Appendix A gives an example), so we do not read those pairs as evidence about the model alone. Cleaner single-dimension manipulations are future work.
Because the gold set carries no rubric-level human ratings, we cannot test whether shared control reflects the models or the halo structure of human analytic ratings (Discussion); a disattenuated human inter-scoring rubric correlation matrix would let a future study separate the two.
Three statistical limits apply. First, the RQ1 correction for layer selection is Bonferroni, which is conservative because adjacent layers are correlated, and 1,000 label shuffles set a floor on each p-value; a max-over-layers permutation test would be exact. Second, the paired-cell analysis is not a simple random sample: layers, alpha values, and models sharing a training lineage are correlated. The primary inference averages within model before pooling, and the family-clustered bootstrap gives an essentially unchanged estimate, but a full hierarchical model across family, essay, scoring rubric, and layer would be a stronger test. Third, the pairwise check removes overlap between two directions at a time, and pairs too collinear to leave a residual are skipped, so the most entangled pairs are the ones it cannot speak to. How best to isolate one direction’s construct-specific component from a correlated one is itself an active methodological question (Elazar et al., 2021).
LLM-based essay scoring is attractive but opaque, and high-stakes testing asks for control and verifiability. We examined the scorer from the inside, asking whether an LLM scorer’s scoring rubrics are represented and independently controllable. Across 25 open-weight base models, holistic scores are decodable (clearly for the two longer task types, weakly for the shortest) and scoring rubrics are partly separable, but neither property implies independent causal control: steering a rubric direction moves the score as a general lever for the great majority of the scoring rubrics this design can measure reliably, and removing each pair’s overlap does not change that. This may reflect the models, the halo structure of the human ratings they imitate, or both; rubric-level human ratings would be needed to tell them apart. Independent controllability is a diagnostic that any rubric-based LLM scorer should be examined with before its rubric labels are treated as separately meaningful, and on this evidence, base LLMs do not yet offer it. Instruction-tuned and fine-tuned scorers are the natural next comparison.
RQ2’s contrastive directions are estimated from clean-versus-degraded essay pairs (§RQ2 Method). Each degradation is produced by prompting llama3.3:70b with a spec that quotes the official BESTEP low-band descriptor for one scoring rubric and instructs the model to hold every other dimension constant. The first example below degrades one source essay on two adjacent scoring rubrics, and the second shows a cohesion degradation, to make the intended separation (and its limits) concrete.
Part 1, grammar (target: verb form, agreement, word order only).
The grammar mutant changes verb forms, subject-verb agreement, and word order while leaving spelling intact; the spelling/punctuation mutant misspells content words and drops punctuation while leaving grammar intact, but it also drops sentence-initial capitalization (“the place”, “i visit”), a surface feature a grammar-sensitive signal could partly pick up on. This is a concrete illustration of the limitation discussed below: even a mutation that targets one dimension by instruction can leave a residual trace in an adjacent one.
Part 2, cohesion (target: connectives and sentence-joining only).
Here the degradation removes connectives (“First of all”, “In this way”) and splits joined sentences into flat, simple ones, with no change to individual sentences’ grammar or vocabulary, a cleaner separation than the grammar/spelling pair above, consistent with cohesion’s lower same-rubric entanglement in Table 4.
Verification. Every degraded essay was checked by two trained BESTEP essay graders before use. Working separately, each grader was shown an original essay together with its degraded version, was told that the degraded version had been downgraded on one scoring rubric, and had to identify which. All 104 degraded essays (20 for Part 1, 36 for Part 2, and 48 for Part 3, that is 5, 9, or 12 for each part and scoring rubric) were checked this way. Both graders confirmed the intended scoring rubric for every one, a consensus of 100%, and the check was run as a random double sample and as a cross-comparison between the graders. This establishes that the intended scoring rubric is recognizable in each degraded essay. It does not establish that no other dimension changed, and the grammar/spelling example above shows that adjacent dimensions can carry a trace. The degraded essay pairs contain exam responses and are not yet public; they are available from the authors for research purposes on request.
Figure 5 shows the RQ1 depth-profile curve (Figure 2’s pooled version, disaggregated) for each of the 25 models individually, for readers who want to check whether the pooled pattern in Figure 2 holds for a specific model rather than only in aggregate.
Table 6 (see page 10) lists the 25 open-weight base models probed and steered in this study, sorted by parameter count.
| Model | Org | Params (B) | Layers | Published |
|---|---|---|---|---|
| SmolLM2-135M | Hugging Face | 0.13 | 30 | 2024-10 |
| gemma-3-270m | 0.27 | 18 | 2025-08 | |
| SmolLM2-360M | Hugging Face | 0.36 | 32 | 2024-10 |
| gemma-3-1b-pt | 1.00 | 26 | 2025-03 | |
| Llama-3.2-1B | Meta | 1.24 | 16 | 2024-09 |
| OLMo-2-0425-1B | AI2 | 1.48 | 16 | 2025-04 |
| Falcon3-1B-Base | TII | 1.67 | 18 | 2024-12 |
| SmolLM2-1.7B | Hugging Face | 1.71 | 24 | 2024-10 |
| granite-3.3-2b-base | IBM | 2.53 | 40 | 2025-04 |
| gemma-2-2b | 2.61 | 26 | 2024-07 | |
| SmolLM3-3B-Base | Hugging Face | 3.08 | 36 | 2025-07 |
| Llama-3.2-3B | Meta | 3.21 | 28 | 2024-09 |
| granite-4.1-3b-base | IBM | 3.40 | 40 | 2026-04 |
| gemma-3-4b-pt | 4.30 | 34 | 2025-03 | |
| Llama-2-7b-hf | Meta | 6.74 | 32 | 2023-07 |
| TAIDE-LX-7B | TAIDE (NARLabs/NCHC) | 6.94 | 32 | 2024-05 |
| Mistral-7B-v0.3 | Mistral AI | 7.25 | 32 | 2024-05 |
| Falcon3-7B-Base | TII | 7.46 | 28 | 2024-12 |
| Apertus-8B-2509 | Swiss AI Initiative | 8.00 | 32 | 2025-09 |
| Llama-3.1-8B | Meta | 8.03 | 32 | 2024-07 |
| granite-3.3-8b-base | IBM | 8.17 | 40 | 2025-04 |
| granite-4.1-8b-base | IBM | 8.38 | 40 | 2026-04 |
| Mistral-NeMo-Minitron-8B-Base | NVIDIA | 8.41 | 40 | 2024-08 |
| gemma-2-9b | 9.24 | 42 | 2024-06 | |
| Falcon3-10B-Base | TII | 10.31 | 40 | 2024-12 |
Duolingo Research Report DRR-26-09, available at https://englishtest.duolingo.com/research.↩︎