Scaling Therapy Supervision: Assessing Psychotherapy Fidelity Using Large Language Models
Authors: Arvind Pillai, Alexis Roach, Jason Holden, Varsha D Badal, Subigya Nepal, Yu Yvonne Wu, Yuliang Chen, Sudarshan Regmi, Sarah Masud Preum, Nicholas C Jacobson, Nikhil Singh, Eric Granholm, Colin Depp, Andrew Campbell
Organizations: Dartmouth College · University of California, San Diego · University of Virginia
Psychotherapy fidelity assessment measures how faithfully a therapist adheres to a treatment protocol, a critical determinant of positive patient outcomes. However, fidelity assessment is bottlenecked by manual expert review, rendering the process labor-intensive and unscalable. To scale therapy supervision and support clinician self-learning, we present an integrated application of large language models (LLMs) that automates psychotherapy fidelity workflows and generates actionable supervisory feedback from authentic, long-form clinical transcripts. We evaluate three open-weight models across multiple task decomposition configurations and reasoning strategies, demonstrating that evidence-first reasoning improves fidelity estimation. Furthermore, we establish the framework's external construct validity on an independent randomized controlled trial (RCT), where it reliably distinguishes Cognitive Behavioral Therapy (CBT) from non-CBT sessions across core therapeutic skills. Through a user study that evaluates the generated feedback, mental health experts indicate that the pipeline synthesizes actionable clinical feedback. Ultimately, our work offers a practical, scalable solution toward bridging the community mental health supervision gap.
Figures & tables
Figure 1: Scaling therapy supervision. Left: Manual review of a 60-minute session is highly time-intensive, limiting scalability. Right: The proposed automated pipeline leverages LLM reasoning over de-identified transcripts to automatically generate fidelity ratings and actionable feedback.
Decomposition
Metric
Qwen3.5-9B
Gemma4-E4B-it
Gemma4-26B-A4B-it
RZ
EV
RF
RZ
EV
RF
RZ
EV
RF
Category
BA
0.67 [0.64–0.69]
0.72 [0.70–0.74]
0.65 [0.62–0.67]
0.62 [0.60–0.64]
0.61 [0.59–0.63]
0.57 [0.56–0.59]
0.58 [0.56–0.60]
0.67 [0.65–0.69]
0.67 [0.64–0.69]
MCC
0.38 [0.33–0.42]
0.45 [0.41–0.49]
0.29 [0.24–0.33]
0.34 [0.29–0.39]
0.35 [0.30–0.39]
0.22 [0.17–0.27]
0.32 [0.29–0.36]
0.43 [0.39–0.48]
0.39 [0.34–0.43]
Clinical
BA
0.64 [0.62–0.66]
0.64 [0.61–0.67]
0.61 [0.59–0.65]
0.64 [0.63–0.67]
0.66 [0.64–0.68]
0.65 [0.63–0.67]
0.56 [0.55–0.57]
0.69 [0.67–0.71]
0.70 [0.68–0.72]
MCC
0.35 [0.30–0.39]
0.31 [0.24–0.37]
0.24 [0.20–0.33]
0.38 [0.33–0.46]
0.39 [0.34–0.44]
0.38 [0.32–0.42]
0.28 [0.23–0.31]
0.46 [0.42–0.50]
0.45 [0.40–0.50]
Section
BA
0.65 [0.63–0.67]
0.65 [0.63–0.67]
0.65 [0.62–0.67]
0.67 [0.65–0.69]
0.68 [0.66–0.70]
0.67 [0.65–0.69]
0.53 [0.52–0.54]
0.63 [0.61–0.65]
0.68 [0.65–0.70]
Table 1: Evaluation of reasoning and decomposition strategies for fidelity estimation (T1). Performance metrics report Balanced Accuracy (BA) and Matthews Correlation Coefficient (MCC), with 95% confidence intervals provided in brackets. Reasoning strategies: Rationalization (RZ), Evidential Estimation (EV), and Reflection (RF).
Figure 2: Comparing LLM fidelity estimation performance vs. human performance (T1). Performance reports Cohen’s κ .
Figure 3: Prediction Powered Inference of overall session fidelity prevalence using Qwen3.5-9B (T2). Red crosses denote the biased, uncorrected raw population averages ( μunlabeled ) produced by each model configuration. The colored dots and corresponding horizontal brackets represent the unbiased PPI rectified means ( μ^PPI ) and their 95% confidence intervals. Grey dashed lines illustrate the direction and magnitude of the statistical rectifier, defined as Δ=μ^PPI−μunlabeled .
Figure 4: Comparison of specific fidelity skill scores between control and intervention groups (T3). The framework’s ability to distinguish between the control group (Non-CBT, μC ) and the intervention group (CBT, μI ) is evaluated across five specific section skills. Statistical significant indicated as ∗∗:p<0.01 , ∗:p<0.05 .
Figure 5: Absolute quality evaluation of generated feedback (T6 ). Mental health experts rated feedback on a 1–5 Likert scale across actionability, clarity, and usefulness.
Appendix figures & tables14 assets
Supplementary material from the paper’s appendix.
Appendix
Category
Item
Description
a) Agenda
1
Therapist noted client’s current emotional state regarding agenda setting.
2
Therapist and client established agenda for session.
3
Priorities for agenda items were established.
4
Agenda was appropriate for time allotment (neither too ambitious nor too limited).
5
Client discussed salient events or problems occurring since last session.
6
The agenda was adhered to during the session.
Appendix
Table 2: CTS-Psy Section I: General Items
Category
Item
Description
f) Guided
1
Used questions to determine meaning client attached to an event.
Discovery
2
Used questions to show incongruities without demeaning the person.
3
Used questions to help client explore various facets of a problem.
4
Used questions to examine client’s arbitrary conclusions or assumptions.
5
Used questions to elicit alternative ways of solving a problem.
6
Used questions to consider alternative explanations.
Appendix
Table 3: CTS-Psy Section II: Specific Items. Note that we do not use the overall score category j).
Study
#Dialogues
#Turns (Form)
Type
#Ratings
Cactus ( Lee et al., 2024 )
31577
16.6 (Short)
Synthetic
0 (unrated)
TherapyGym ( Huang et al., 2026 )
116
10 (Short)
Synthetic
1276
Ours
30 (420 full)
∼ 180 ± 81 (Long)
Real
1800
Appendix
Table 4: Comparison of prior evaluations and data for dialogue-based therapy.
Decomposition
Qwen3.5-9B
Gemma4-E4B-it
Gemma4-26B-A4B-it
RZ
EV
RF
RZ
EV
RF
RZ
EV
RF
Category
0.31
0.46
0.26
0.26
0.29
0.15
0.17
0.44
0.35
Clinical
0.32
0.33
0.32
0.32
0.37
0.36
0.20
0.47
0.42
Section
0.32
0.28
0.36
0.36
0.36
0.36
0.14
0.33
0.43
Obs
0.36
0.31
0.30
0.30
0.38
0.28
0.09
0.45
0.47
Full
0.26
0.26
0.36
0.36
0.33
0.33
0.16
0.44
0.50
Appendix
Table 5: Comparing LLM fidelity estimation performance vs. human performance (T1). Performance reports Cohen’s κ . Human performance upper bound is κ=0.71 . Reasoning strategies evaluated include Rationalization (RZ), Evidential Estimation (EV), and Reflection (RF).
Category
BA
MCC
Agenda
0.68 [0.61–0.75]
0.39 [0.24–0.54]
Feedback
0.54 [0.47–0.62]
0.08 [-0.06–0.23]
Understanding
0.63 [0.48–0.78]
0.26 [-0.04–0.54]
Interpersonal
0.78 [0.44–0.96]
0.22 [-0.04–0.44]
Collaboration
0.53 [0.48–0.59]
0.11 [-0.06–0.28]
Discovery
0.55 [0.46–0.63]
0.08 [-0.06–0.21]
Appendix
Table 6: Qwen3.5-9B with evidential reasoning and category level decomposition performance across categories.
Figure 6: Prediction Powered Inference of overall session fidelity prevalence using Gemma4-E4B-it (T2). Red crosses denote the biased, uncorrected raw population averages ( μunlabeled ) produced by each model configuration. The colored dots and corresponding horizontal brackets represent the unbiased PPI rectified means ( μ^PPI ) and their 95% confidence intervals. Grey dashed lines illustrate the direction and magnitude of the statistical rectifier, defined as Δ=μ^PPI−μunlabeled .
Figure 7: Prediction Powered Inference of overall session fidelity prevalence using Gemma4-26B-A4B-it (T2). Red crosses denote the biased, uncorrected raw population averages ( μunlabeled ) produced by each model configuration. The colored dots and corresponding horizontal brackets represent the unbiased PPI rectified means ( μ^PPI ) and their 95% confidence intervals. Grey dashed lines illustrate the direction and magnitude of the statistical rectifier, defined as Δ=μ^PPI−μunlabeled .
#
Role
Degree
SMI
CBT
Cog.
1
Therapist
Master’s
Extremely Familiar
Extremely Familiar
Extremely Familiar
2
Clinical Psychologist
Doctorate
Extremely Familiar
Extremely Familiar
Extremely Familiar
3
Therapist
Doctorate
Very Familiar
Extremely Familiar
Extremely Familiar
Appendix
Table 7: Expertise information of mental health experts.
Figure 8: Post-Hoc Estimation (Direct-then-Rationalize). The LLM immediately assigns a binary score yi per item, then provides a brief rationale as post-hoc evidence. The user turn injects the transcript ⟨ T ⟩ , the decomposition-specific instruction (Table 8 ), and the CTS-Psy item set Q .
Figure 9: Evidential Estimation (Extract-Assess-Predict). A grounded three-step chain: the model first extracts verbatim relevant_spans from the transcript, then assesses them against the criterion’s operational definition, and only then assigns the binary score yi . This enforces an evidence-first reasoning order.
Figure 10: Self-Reflection (Critic-in-the-Loop). After an initial Evidential Estimation pass (Step 1), a separate reflection LLM instantiated as an expert clinical supervisor critiques each rating across rubric adherence, evidence coverage, and confidence calibration (Step 2). The consolidated critique is fed back to the primary generator, which then revises its binary decisions and supporting evidence (Step 3). This loop may be applied for k≥1 iterations.
Decomposition
Partition
∣Q∣
Task Instruction
Full
—
60
Rate the transcript based on all 60 CTS-Psy items below, organized by section:
Observability
Explicit
24
Rate the explicitly observable items below (Agenda, Feedback, Homework, Social):
Rate the 30 General CTS-Psy items below (Agenda, Feedback, Understanding, Interpersonal Effectiveness, Collaboration):
Specific
30
Rate the 30 Specific CTS-Psy items below (Discovery, Cognition, Intervention, Homework, Social Skills Training):
Clinical
Relationship
18
Rate the Therapeutic Relationship items below (Understanding, Interpersonal Effectiveness, Collaboration):
Appendix
Table 8: Decomposition-specific task instructions injected into the user turn of each prompt (Figures 8 – 10 ). Category-level decomposition issues 10 independent inference passes, one per CTS-Psy category (Agenda, Feedback, Understanding, Interpersonal Effectiveness, Collaboration, Discovery, Cognition, Intervention, Homework, Social).
Figure 11: Feedback Generation. Given the three lowest-scoring CBT competency categories from the fidelity-estimation stage (§4.2), the LLM produces one positive observation and one actionable, session-specific suggestion per weak competency area. The user turn injects ⟨ bottom-3 fidelity evidence ⟩ : the per-item fidelity statement, score, relevant transcript spans, and criterion assessment for each of the 18 items spanning the three weakest categories. In the RAG condition, retrieved therapy-manual passages relevant to this evidence are additionally prepended to the user turn before generation.
Figure 12: LLM-as-a-Judge (Paired Comparison). For each dimension (Actionability, Clarity), a judge LLM independently compares RAG-grounded and non-RAG feedback generated for the same session (Figure 11 ), presented order-randomized as Feedback A/B, and returns a single comparative score ∈{−1,0,+1} . Each dimension is scored by a separate model call so that no single judge call sees more than one dimension at once.
Department of Computer Science, Stanford University, Stanford, CA, USA · Department of Psychiatry and Behavioral Sciences, Stanford University, Stanford, CA, USA · The University of Hong Kong, Hong Kong +2
CogNosco Lab, Department of Psychology and Cognitive Science, University of Trento, Italy. · Department of Education, Psychology and Communication Sciences, University of Bari, Italy.