Preprint
Concept Paper

This version is not peer-reviewed.

PainVector: A Multidimensional Computational Framework for the Digital Representation of Pain

Submitted:

27 September 2026

Posted:

29 September 2026

You are already at the latest version

Abstract
Pain is a multidimensional clinical phenomenon influenced by biological, psychological, behavioral, and contextual factors. Its assessment remains challenging because pain is a personal experience that cannot be directly observed and because commonly used clinical instruments, such as the Numeric Rating Scale (NRS) and Visual Analog Scale (VAS), primarily reduce pain to a unidimensional scalar value. Physiological signals such as heart rate variability (HRV), electrodermal activity (EDA), and peripheral temperature may provide complementary information, but these signals are themselves nonspecific and require contextual interpretation. We propose PainVector, a computational framework for representing pain as a structured multidimensional state rather than as a single scalar intensity value. The framework integrates self-reported pain, normalized physiological signals, and contextual variables including age, cultural context, and individual pain-expression characteristics. PainVector separates perceived pain from physiologically inferred pain and subsequently combines these components through an explicit fusion mechanism. The framework additionally incorporates temporal and system-level dimensions, including pain duration, predictability, controllability, localization, semantic interpretation, and confidence. A mathematical formulation is provided for signal normalization, contextual modulation, functional weighting, multimodal fusion, temporal prediction, and controllability analysis. A synthetic dataset is used as a proof-of-concept demonstration of the computational structure, while a two-person migraine scenario illustrates how individuals reporting identical NRS values may produce different multidimensional representations when contextual and physiological information are incorporated. PainVector is not proposed as a direct measurement of subjective consciousness, nor as a replacement for clinical pain assessment. Rather, it is designed as a computational representation layer that may support multimodal analysis, longitudinal monitoring, prediction, intervention-response evaluation, and future human–machine pain modeling.
Keywords: 
;  ;  ;  ;  ;  ;  ;  ;  ;  

1. Introduction

Pain is one of the most complex phenomena encountered in clinical medicine and neuroscience. Contemporary definitions emphasize that pain is a personal experience influenced by biological, psychological, and social factors, and that pain should not be equated with nociceptive activity alone [1].
This complexity creates a fundamental measurement problem. Traditional clinical instruments such as the NRS and VAS provide useful and clinically established estimates of perceived pain intensity. However, they compress a multidimensional phenomenon into a single scalar value. For example, two individuals may both report NRS = 7 while differing substantially in age, cultural background, behavioral expression, physiological response, pain duration, or response to intervention. Therefore, although a scalar intensity score is clinically useful, it may not fully describe the structure or dynamics of an individual’s pain state.
Physiological signals provide a potentially complementary source of information. HRV, EDA, and peripheral temperature can capture aspects of autonomic and physiological regulation [2,3]. However, these signals are not pain-specific and may also respond to stress, emotion, physical activity, environmental conditions, medication, and other physiological processes [4,5]. Consequently, the problem is not simply to replace subjective pain assessment with physiological measurement; instead, a computational framework is needed that can represent multiple heterogeneous sources of pain-related information within a common mathematical structure.
This paper introduces PainVector, a multidimensional computational representation framework designed for this purpose. The central proposition is: pain can be represented computationally as a structured multidimensional state without claiming that the representation directly captures subjective consciousness or the complete biological mechanism of pain. PainVector therefore focuses on representation, rather than claiming direct measurement of subjective experience.

3. Problem Definition

The principal challenges addressed by PainVector are:
1.
Pain is not directly observable.
2.
Pain reports are individual-dependent and context-dependent.
3.
A single scalar such as NRS cannot represent all relevant dimensions of a pain state.
4.
Physiological signals are heterogeneous and noisy.
5.
Physiological signals are not uniquely specific to pain.
6.
Individuals may differ in their expression and reporting of pain.
7.
Pain is dynamic rather than static.
8.
Prediction and response to intervention require temporal information.
9.
Multimodal pain data require a common computational representation.
The objective of PainVector is therefore: to construct a structured, machine-readable, multidimensional representation of pain-related state information that can be analyzed, compared, tracked over time, and potentially used for prediction and intervention evaluation.
The framework is not intended to: directly measure subjective consciousness; establish the presence of pain solely from physiological signals; replace clinical assessment; establish a universal biological equation for pain; or claim that cultural or behavioral coefficients can be assumed without empirical calibration. These distinctions are essential because pain and nociception are not equivalent phenomena [1].

4. Scope and Conceptual Architecture

PainVector is designed as a representation layer between heterogeneous pain-related observations and downstream computational applications. The architecture consists of five conceptual layers.
Layer 1 — Input Layer. The framework receives self-reported pain, physiological signals, contextual variables, temporal information, and optional clinical information. The initial physiological implementation includes HRV, EDA, and peripheral temperature. Future extensions may include EEG, EMG, facial expression, movement, respiratory variables, neuroimaging, and wearable-derived measurements.
Layer 2 — Normalization Layer. Raw physiological measurements are transformed relative to an individual’s baseline.
Layer 3 — Contextual and Functional Modulation Layer. Individual and contextual variables modify the interpretation and weighting of the available measurements.
Layer 4 — Multimodal Fusion Layer. Perceived and physiologically inferred components are combined into a unified pain-state representation.
Layer 5 — State Representation and Analysis Layer. The resulting PainVector can be used for longitudinal monitoring, comparison, prediction, intervention-response evaluation, state-transition analysis, localization, semantic interpretation, confidence estimation, and future human–machine interaction.

5. Mathematical Representation of PainVector

PainVector is defined as
PV t = I perceived , I inferred , I final , D , P , C t l , L , M , C o n f t
where I perceived is the context-adjusted perceived pain intensity, I inferred is the physiologically inferred pain-related intensity, I final is the fused pain representation, D is duration, P is predictability, C t l is controllability, L is localization, M is semantic meaning/interpretation, and C o n f is the confidence or uncertainty associated with the representation. The vector is therefore not intended to represent pain as a single number; instead, it represents pain as a multidimensional state.

6. Methods

6.1. Signal Normalization

For each physiological feature x i , an individual baseline is estimated:
μ i = E [ x i ∣ baseline ] , σ i = SD [ x i ∣ baseline ]
The standardized deviation is
z i = x i − μ i λ i σ i + ϵ
where x i is the observed feature, μ i is the individual’s baseline mean, σ i is the baseline standard deviation, λ i is a feature-specific scaling parameter, and ϵ prevents numerical instability.
The standardized value is transformed into a bounded representation:
s i = clip f i ( z i ) , 0 , 1
where f i is a feature-specific transformation. This distinction is important because the direction of physiological change is not necessarily identical across signals: for example, a decrease in a particular HRV measure may be interpreted differently from an increase in EDA. The framework therefore does not assume a common monotonic relationship across raw physiological variables; instead s i = f i ( x i , μ i , σ i ) is calibrated empirically for each feature. The normalized physiological representation becomes
s = [ s H R V , s E D A , s T e m p , … , s n ]

6.2. Perceived Pain Representation

The basic self-reported pain value is normalized as
s self = NRS 10 , 0 ≤ s self ≤ 1
Because self-report may be influenced by contextual and behavioral factors, PainVector permits contextual modulation. Let C denote a context/cultural calibration factor and E an individual pain-expression calibration factor. Then
I perceived = clip s self C · E , 0 , 1
The parameters C and E are not assumed to be universal constants; they represent calibration parameters that would require empirical estimation from appropriate populations. In a clinically validated implementation, these parameters should be estimated from data rather than manually assigned.

6.3. Physiological Inference

The physiologically inferred component is defined as
I inferred = ∑ i = 1 n w i s i , w i ≥ 0 , ∑ i = 1 n w i = 1
For the initial implementation,
I inferred = w H R V s H R V + w E D A s E D A + w T e m p s T e m p
An illustrative configuration may use w H R V = 0.30 , w E D A = 0.40 , w T e m p = 0.30 . These values are illustrative model parameters rather than clinically validated weights; they should ultimately be estimated or optimized using empirical data.

6.4. Contextual and Functional Weighting

PainVector allows physiological components to receive different functional weights according to individual and contextual conditions:
w ˜ i = w i 0 · A i · g i ( C , E ) · h i ( B i ) ∑ j = 1 n w j 0 · A j · g j ( C , E ) · h j ( B j )
where w i 0 is the initial weight, A i represents age-dependent modulation, g i ( C , E ) represents contextual modulation, h i ( B i ) represents baseline-deviation modulation, and B i represents the magnitude of deviation from the expected baseline state.
Age-dependent modulation is represented generally as A i = A i ( age ) rather than assuming a universal linear relationship. For a simple illustrative model,
A i ( age ) = 1 + α i age − age 0 100
where α i controls the magnitude of age-dependent modulation and age 0 is a reference age. This formulation preserves the age-dependent weighting concept while avoiding a common multiplier that would cancel during normalization.

6.5. Baseline Deviation

The deviation of each component from its expected state is
B i = | s i − μ i * |
where μ i * represents the expected normalized baseline value. The functional modulation is then
h i ( B i ) = 1 + κ i B i , κ i ≥ 0
where κ i controls the sensitivity of the model to deviations from baseline. This allows the framework to distinguish an individual’s absolute physiological value from the individual’s deviation from their own expected state, which is central to the personalized nature of PainVector.

6.6. Final Pain Representation

The perceived and inferred components are combined through a fusion coefficient:
I final = η I perceived + ( 1 − η ) I inferred , 0 ≤ η ≤ 1
The parameter η controls the relative contribution of self-report and physiological information; η = 0.5 produces equal weighting. η should ultimately be calibrated empirically rather than assumed to be universally optimal. I final is a computationally fused representation and should not be interpreted as a direct measurement of subjective pain.

6.7. Temporal Dimension: Pain Duration

If M pain episodes are observed and the duration of episode m is τ m , the mean episode duration is
D = 1 M ∑ m = 1 M τ m
Alternatively, the framework can preserve the complete temporal sequence, D t = τ t , allowing PainVector to represent dynamic pain trajectories rather than only an average duration.

6.8. Predictability

Given an observed future pain state I t + 1 and a prediction I ^ t + 1 | t , the prediction error is
MSE = 1 T ∑ t = 1 T I t + 1 − I ^ t + 1 | t 2
A normalized predictability score is defined as
P = 1 − MSE Var ( I t + 1 ) + ϵ
where ϵ prevents division by zero. This formulation is analogous to a variance-normalized prediction score. Because the resulting value can become negative when predictions perform worse than a simple variance-based reference, the model should not automatically assume 0 ≤ P ≤ 1 unless an additional bounding transformation is applied; if a bounded representation is required, P bounded = clip ( P , 0 , 1 ) can be used.

6.9. Controllability

Pain can also be represented as a dynamic system influenced by interventions. A simplified state-transition model is
I t = β 0 + β 1 I t − 1 + β 2 a t − 1 + ϵ t
where I t − 1 is the previous state, a t − 1 represents an intervention or action, β 0 , β 1 , β 2 are model parameters, and ϵ t represents unexplained variation. A conceptual controllability index is
C t l = Var ( I t ∣ observed ) − Var ( I t ∣ optimal action ) Var ( I t ∣ observed ) + ϵ
Controllability reflects the potential reduction in state variability achievable through an intervention policy. This component should be interpreted as a model-based property, not as a direct physiological measurement; its clinical validity requires longitudinal intervention data.

6.10. Localization

Pain location can be represented probabilistically:
L = arg max l P ( location = l ∣ F )
where l represents a candidate anatomical region and F represents available features. A future implementation may represent localization as a full distribution L = [ P ( l 1 ) , P ( l 2 ) , … , P ( l k ) ] rather than retaining only the maximum-probability location.

6.11. Semantic Meaning

Pain can also be associated with contextual or semantic interpretations:
M = arg max m P ( meaning = m ∣ F )
where m represents a semantic category (e.g., acute, chronic, nociceptive, neuropathic, musculoskeletal, procedure-related, stress-associated). This component should not be interpreted as a definitive clinical diagnosis; it represents a future computational semantic layer that requires supervised clinical validation.

6.12. Confidence and Uncertainty

PainVector includes an explicit confidence dimension:
C o n f = λ C o n f var + ( 1 − λ ) C o n f entropy
where C o n f var represents confidence derived from model variance, C o n f entropy represents confidence derived from predictive uncertainty or entropy, and λ controls their relative contribution. Alternatively, uncertainty may be represented directly as U = 1 − C o n f . This allows the system to distinguish a high pain estimate with high confidence from a high pain estimate with substantial uncertainty, a distinction that is important for clinical decision-support systems.

7. Complete PainVector Definition

The complete representation is therefore
PV t = I perceived , I inferred , I final , D , P , C t l , L , M , C o n f t
which can be interpreted as the pipeline
Pain → Multimodal Observations → Normalized State → Contextual Integration → Fusion → Multidimensional Representation

8. Computational Pipeline

The complete computational pipeline proceeds as follows: (1) baseline estimation ( μ i , σ i ); (2) signal normalization ( x i → s i ); (3) perceived pain representation ( NRS → I perceived ); (4) physiological inference ( s → I inferred ); (5) contextual weighting ( w i → w ˜ i ); (6) multimodal fusion ( I perceived , I inferred → I final ); (7) dynamic representation ( D , P , C t l ); (8) extended representation ( L , M , C o n f ); producing the final output PV t .

9. Synthetic Proof-of-Concept Dataset

To demonstrate the internal computational structure, a synthetic dataset of ten simulated individuals is generated using the open-source script accompanying this paper (painvector_demo.py; fixed random seed for reproducibility). The data are synthetically generated and are not clinical patient data. An example run is shown in Table 1.
These data are intended only to demonstrate how heterogeneous inputs can be transformed into a common computational representation. They do not demonstrate clinical accuracy, diagnostic performance, or physiological validity.

10. Migraine Demonstration Scenario

To illustrate the multidimensional nature of the framework, consider two hypothetical individuals with migraine, both reporting NRS = 7 , so s self = 0.7 for both. The individuals differ in contextual and physiological parameters (Table 2).
Thus, despite identical NRS values, the computational representations differ substantially ( I final , A = 0.811 vs. I final , B = 0.535 ). This example demonstrates the mechanism of the framework, not evidence that the resulting difference accurately reflects actual clinical pain; the difference arises purely from the combination of contextual calibration and physiological inputs as defined in Section 6.2, Section 6.3, and Section 6.6.

11. Validation Framework

The present work establishes a computational framework rather than a clinically validated diagnostic system. Future empirical validation should occur in multiple stages: (1) construct validation of individual dimensions (e.g., perceived pain versus established clinical scales, physiological components versus experimentally induced pain, duration versus longitudinal records, localization versus documented anatomical location); (2) convergent validation against established instruments (NRS, VAS, multidimensional pain questionnaires, clinical outcomes); (3) predictive validation of subsequent pain intensity, trajectory, treatment response, exacerbation, and recovery using prospective longitudinal datasets; (4) sensitivity analysis over η , w i , λ i , κ i , and contextual parameters; and (5) external validation across age groups, sexes, cultural populations, pain conditions, clinical environments, wearable devices, and measurement protocols.

12. Discussion

PainVector proposes a shift from scalar pain scoring toward multidimensional computational representation. Traditional instruments answer a relatively narrow question — how intense is the pain according to the patient’s report — whereas PainVector attempts to represent a broader state: the current computationally observable configuration of pain-related information for a given individual. This distinction is central to the framework.
The model does not assume that physiological signals constitute pain itself; physiological measurements are treated as complementary observations that may contribute information about the state of the individual. This distinction is particularly important because physiological responses can be influenced by many factors unrelated to pain. As reviewed in Section 2, contemporary machine-learning approaches to pain assessment continue to investigate physiological and multimodal markers while highlighting methodological heterogeneity and the need for stronger external and clinical validation [6,10]. PainVector therefore positions multimodal physiological information as part of a representation architecture rather than as an objective replacement for subjective reporting, and, unlike the direct classifiers surveyed above, is explicitly designed as an intermediate, interpretable state onto which such estimators could in principle be mapped.
Another important property of the framework is individual baseline normalization: rather than assuming that the same physiological value has the same meaning for every person, PainVector compares measurements against an individual’s own baseline, creating a pathway toward personalized computational pain representation. The framework also introduces temporal dimensions, modeling Pain ( t 1 ) , Pain ( t 2 ) , Pain ( t 3 ) , … rather than an isolated measurement, enabling trajectory → prediction → response → control analysis and thus the possibility of moving from static assessment toward dynamic pain-state modeling.

13. Relationship to AI and Human–Machine Systems

PainVector can function as an intermediate representation layer for future AI systems. Rather than feeding raw heterogeneous signals directly into a black-box model, a computational system could receive PV t as a structured state representation, yielding the conceptual pipeline
Raw Data → PainVector → AI Model → Prediction / Decision
Potential applications include clinical decision support, longitudinal monitoring, personalized pain management, digital therapeutics, wearable systems, rehabilitation, and human–machine interaction. However, such future applications should be distinguished from claims about machine consciousness or artificial subjective experience: PainVector itself makes no claim that a machine possessing a PainVector necessarily experiences pain.

14. Limitations

Several limitations must be acknowledged. (1) Synthetic demonstration data: the present proof-of-concept dataset is synthetic and therefore cannot establish clinical validity. (2) Physiological non-specificity: HRV, EDA, and temperature are influenced by many physiological and psychological processes and are not specific biomarkers of pain. (3) Parameter calibration: parameters such as C , E , η , w i , λ i , κ i currently represent framework parameters and require empirical estimation. (4) Cultural modeling: cultural effects on pain expression are complex and should not be reduced to simplistic universal coefficients. (5) Individual expression: pain expression varies across individuals, making calibration essential. (6) Simplified dynamic model: the current controllability and prediction equations are intentionally simplified representations of potentially much more complex nonlinear systems. (7) Localization and meaning: these dimensions require appropriate supervised datasets and clinical annotation. (8) No claim regarding consciousness: PainVector does not model subjective consciousness, suffering, phenomenology, or artificial sentience. (9) Clinical validation: remains necessary before the framework can be considered for clinical decision-making.

15. Future Work

Future development of PainVector may include clinical data collection; longitudinal datasets; adaptive individual baseline estimation; machine-learning-based parameter estimation; EEG, EMG, and respiratory-signal integration; facial-expression analysis; wearable-device integration; personalized calibration; nonlinear dynamical modeling; reinforcement-learning-based intervention optimization; multicenter and cross-cultural validation; and integration with personalized medicine platforms. A particularly important future direction is an adaptive PainVector,
P V t + 1 = F ( P V t , X t , A t )
where P V t represents the current pain state, X t represents new multimodal observations, and A t represents an intervention or action — transforming PainVector from a static representation into a dynamic state-space model.

16. Conclusion

PainVector introduces a multidimensional computational framework for representing pain-related state information. Instead of reducing pain to a single scalar value, the framework integrates self-report, physiological signals, context, and temporal dynamics into the structured representation PV t defined above. The framework separates perceived pain from physiologically inferred information and provides an explicit multimodal fusion mechanism, further introducing temporal, predictive, controllability, localization, semantic, and uncertainty dimensions. The current work should be considered a theoretical and computational framework with synthetic proof-of-concept demonstration, rather than a clinically validated pain measurement system. Its principal contribution is therefore not the claim that pain has been solved computationally, but the proposal of a structured mathematical representation through which heterogeneous pain-related observations can be integrated, analyzed, compared, and potentially modeled dynamically. Future empirical validation will determine whether the proposed representation provides clinically meaningful advantages over existing pain assessment approaches.

Funding

This research received no external funding.

Data Availability Statement

The synthetic dataset generator and migraine demonstration script (painvector_demo.py) are provided alongside this preprint to allow full reproduction of Table 1 and Table 2.

Conflicts of Interest

The author declares no conflict of interest.

Use of Artificial Intelligence

The author used an AI assistant (Claude, Anthropic) for language translation/editing, LaTeX formatting, code implementation, and literature search assistance during manuscript preparation. The scientific framework, mathematical formulation, and conclusions are the author’s own.

References

  1. Raja, S.N.; Carr, D.B.; Cohen, M.; Finnerup, N.B.; Flor, H.; Gibson, S.; Keefe, F.J.; Mogil, J.S.; Ringkamp, M.; Sluka, K.A.; et al. The revised International Association for the Study of Pain definition of pain: concepts, challenges, and compromises. Pain 2020, 161, 1976–1982. [Google Scholar] [CrossRef] [PubMed]
  2. Shaffer, F.; Ginsberg, J.P. An Overview of Heart Rate Variability Metrics and Norms. Front. Public Health 2017, 5, 258. [Google Scholar] [CrossRef] [PubMed]
  3. Boucsein, W. Electrodermal Activity, 2nd ed.; Springer, 2012. [Google Scholar]
  4. Forte, G.; Troisi, G.; Pazzaglia, M.; De Pascalis, V.; Casagrande, M. Heart Rate Variability and Pain: A Systematic Review. Brain Sci. 2022, 12, 153. [Google Scholar] [CrossRef] [PubMed]
  5. Kong, Y.; Chon, K.H. Electrodermal activity in pain assessment and its clinical applications. Appl. Phys. Rev. 2024, 11, 031316. [Google Scholar] [CrossRef]
  6. Gkikas, S.; Tsiknakis, M. Automatic assessment of pain based on deep learning methods: A systematic review. Comput. Methods Programs Biomed. 2023, 231, 107365. [Google Scholar] [CrossRef] [PubMed]
  7. Pouromran, F.; Radhakrishnan, S.; Kamarthi, S. Exploration of physiological sensors, features, and machine learning models for pain intensity estimation. PLoS ONE 2021, 16, e0254108. [Google Scholar] [CrossRef] [PubMed]
  8. Fernandez-Rojas, R.; Hirachan, N.; Brown, N.; Waddington, G.; Murtagh, L.; Seymour, B.; Goecke, R. Multimodal physiological sensing for the assessment of acute pain. Front. Pain Res. 2023, 4, 1150264. [Google Scholar] [CrossRef] [PubMed]
  9. Gozzi, N.; Preatoni, G.; Ciotti, F.; Hubli, M.; Schweinhardt, P.; Curt, A.; Raspopovic, S. Unraveling the physiological and psychosocial signatures of pain by machine learning. Med 2024, 5, 1495–1509. [Google Scholar] [CrossRef] [PubMed]
  10. Cascella, M.; Schiavo, D.; Cuomo, A.; et al. Artificial Intelligence for Automatic Pain Assessment: Research Methods and Perspectives. Pain Res. Manag. 2023. [Google Scholar] [CrossRef] [PubMed]
Table 1. Synthetic proof-of-concept dataset ( n = 10 ), generated by painvector_demo.py with a fixed random seed. HRV, EDA, and Temp are pre-normalized signals in [ 0 , 1 ] ; I perceived , I inferred , and I final follow Eqs. (7), (9), and (14) respectively.
Table 1. Synthetic proof-of-concept dataset ( n = 10 ), generated by painvector_demo.py with a fixed random seed. HRV, EDA, and Temp are pre-normalized signals in [ 0 , 1 ] ; I perceived , I inferred , and I final follow Eqs. (7), (9), and (14) respectively.
Ind. HRV EDA Temp NRS C E I perceived I inferred
1 0.64 0.03 0.28 4 0.86 0.84 0.554 0.288
2 0.74 0.55 0.59 3 0.81 0.89 0.416 0.619
3 0.51 0.03 0.20 8 1.08 0.97 0.764 0.225
4 0.45 0.28 0.87 9 1.12 1.08 0.744 0.508
5 0.34 0.16 0.96 5 0.84 0.95 0.627 0.454
6 0.36 0.34 0.26 3 1.09 1.01 0.273 0.322
7 0.97 0.38 0.55 9 1.05 1.15 0.745 0.608
8 0.36 0.19 0.07 8 0.89 0.92 0.977 0.205
9 0.08 0.23 0.10 5 0.98 1.13 0.452 0.146
10 0.16 0.36 0.67 8 1.17 1.06 0.645 0.393
Table 2. Migraine demonstration: identical self-reported intensity, different contextual and physiological profiles, different final representations. Values reproduced by painvector_demo.py.
Table 2. Migraine demonstration: identical self-reported intensity, different contextual and physiological profiles, different final representations. Values reproduced by painvector_demo.py.
Subject Age C E HRV EDA Temp I perceived I final
A 25 0.8 0.9 0.7 0.8 0.4 0.972 0.811
B 60 1.2 1.1 0.5 0.6 0.5 0.530 0.535
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.