# Measuring Entropy Change in a Human Physiological System

Authors: Satish Boregowda; Rod Handy; Darrah Sleeth; Andrew Merryweather
Journal: Journal of Thermodynamics
 (2016)

---

## Abstract

The paper presents a novel approach involving the use of Maxwell relations to combine multiple physiological measures to provide a measure of entropy change. The physiological measures included blood pressure (BP), heart rate (HR), skin temperature (ST), electromyogram (EMG), and electrodermal response (EDR). The multiple time-series physiological data were collected from eight subjects in an experimental pilot study conducted at theHuman Engineering Laboratory of NASA Langley Research Center. The methodology included data collection during a relaxation period of eighteen minutes followed by a sixty-minute cognitive task. Two types of entropy change were computed: (a) entropy change (Δ
S
B
P) due to blood pressure, heart rate, and skin temperature and (b) entropy change (Δ
S
E
M
G) due to electromyogram, electrodermal response, and skin temperature. The results demonstrate that entropy change provides a valuable composite measure of individual physiological response to various stressors that could be valuable in the areas of medical research, diagnosis, and clinical practice.


---

## Body

## 1. Introduction

Earlier work done by several researchers has established the fact that human life processes are indeed thermodynamic in nature and hence thermodynamic laws can be used to model human physiology. It has been hypothesized by Bridgman [1] that the laws of thermodynamics are intrinsic to nature and life and are thus well positioned to characterize the physiological behavior of living systems. Aoki [2, 3] examined human thermoregulation by measuring entropy flow and production in basal and exercising conditions but did not utilize Maxwell relations to include other physiological responses. The seminal work done in the areas of biological thermodynamics by Morowitz [4, 5] does not show any application of Maxwell relations to model human physiology. Furthermore, the work done by Garby and Larsen [6] and Jou and Llebot [7] addresses mass, energy, and entropy balance of open living systems but not Maxwell relations. However, the detailed development and preliminary verification of physiological entropy change using Maxwell relations have been presented by Boregowda et al. [8, 9] and are reproduced in Appendix for the convenience of the readers. The present study utilizes time-series data collected during a cognitive task to validate and verify entropy change metrics representing two human physiological subsystems. For the convenience of the readers, the detailed development of physiological entropy change using Maxwell relations is presented in Appendix of this paper. The first one combines changes in blood pressure, heart rate, and skin temperature to provide a BP-based entropy change (Δ
S
B
P). The second one combines changes in electromyogram, electrodermal response, and skin temperature to provide an EMG-based entropy change (Δ
S
E
M
G).


## 2. Modeling and Formulation

### 2.1. Entropy Change (Δ
S
B
P: Blood Pressure, Heart Rate, and Skin Temperature)

Let us first consider the human physiological subsystem characterized by blood pressure (BP), heart rate (HR), and skin temperature (ST) as a closed thermodynamic system. For this system, the work done by the system is given by [10](1)δW=P∗dV=BP∗dHR,where P is pressure, BP is blood pressure (mmHg), V is volume, and HR is heart rate (beats/minute).The heart rate is an indirect measure of stroke volume of blood in the heart region and is, thus, used in the place of stroke volume. Please note that heart rate (in beats per minute) is equal to cardiac output (mL/min) divided by stroke volume (mL/beat). The ratio of partial changes in two physiological variables is accompanied by constancy in the third variable, as this is the basis of Maxwell relations. For example, when there is a partial change in heart rate with respect to partial change in entropy, the mean arterial pressure remains constant. As shown in Appendix, the thermodynamic potential is derived as shown starting with Helmholtz’s potential (F
=
U
-
T
S) and arriving with the following property relation: (2)dF=-BP∗dHR-S∗dST.The following Maxwell relation is obtained by invoking the exactness condition to the above property relation (2):(3)∂BP∂STV=∂S∂HRT.The partial derivative is approximated to form a modified Maxwell relation that is used in the present experimental study to compute the physiological entropy change resulting from changes in blood pressure (BP), heart rate (HR), and skin temperature (ST):(4)ΔBPSTHR=ΔSBPΔHRST,ΔSBP=ΔBP×ΔHRΔST.The skin temperature (ST) in the denominator acts as a temperature at the boundary of the closed thermodynamic system and provides a physiological reflection of emotional response. If one were to consider either internal energy or the Gibbs function to define human physiological function, then one would get a negative sign preceding the entropy change (Δ
S
B
P). The deviation from the relaxed state may be positive or negative depending on the imposed external stressor and the internal physiological condition.


### 2.2. Entropy Change (Δ
S
E
M
G: Electromyogram, Electrodermal Response, and Skin Temperature)

Let us consider the human physiological subsystem characterized by electromyogram (EMG), electrodermal response (EDR), and skin temperature (ST) as a closed thermodynamic system. For this system, the work done by the system is given by(5)δW=ℾ∗dL=EMG∗dEDR,where ℾ is tension, EMG is facial electromyogram (mV), L is length, and EDR is electrodermal response (mV).In this regard, the elements of a human physiological subsystem system are considered as equivalent to a conducting wire withℾ as tension and L as the length. The facial electromyogram (EMG) has been used in medical research and clinical practice as a measure of tension in facial muscles [11]. As a result of tension in muscles or nerves, the outer surface of the skin responds by expanding or contracting in terms of changes called electrodermal response (EDR). The skin temperature (ST) represents the thermal response to tension. The thermodynamic potential is again derived, starting with Helmholtz’s potential (F
=
U
-
T
S) arriving at the following property relation: (6)dF=-EMG∗dEDR-S∗dST.The following Maxwell relation is obtained by invoking the exactness condition to the above property relation (2):(7)∂EMG∂STEDR=∂S∂EDRST.The partial derivative is approximated to form a modified Maxwell relation that is used in the present experimental study to compute the physiological entropy change resulting from changes in blood pressure (BP), heart rate (HR), and skin temperature (ST):(8)ΔEMGSTEDR=ΔSEMGΔEDRST,ΔSEMG=ΔEMG×ΔEDRΔST.As mentioned earlier, the use of either internal energy or the Gibbs function would have resulted in a negative sign preceding the entropy change (Δ
S
P). The deviation from the relaxed state may be positive or negative, depending on the imposed external stressor and the internal physiological condition. Specifically, physiological entropy is qualitatively defined as a measure of disorder [12]. It is a kind of global measure that specifies how violent motions and reactions are occurring in nature. Hence, the entropy change in the human physiological system shows the extent of activity within the body as a whole; thus, the entropy change is a significant quantity that characterizes the human body from thermodynamic and holistic (i.e., considering a human body as a whole) viewpoints.The human physiological entropy change could be considered as a composite measure of change in the whole physiological state in response to any external stimuli or stressor. However, if a single physiological indicator alone, such as the mean arterial pressure, can provide that information, then why do we need this entropy change as a composite measure of physiological response? The answer is as follows: the physiological concepts such as stimulus response (SR) specificity, organ response (OR) specificity, individual response specificity, and autonomic balance make the human physiological response a complex phenomenon [11]. The single physiological indicators, in this regard, provide a very narrow representation of the human physiological stress response system. It is only by recognizing the interaction among human subsystems in their response to any stressor stimuli that one could build a better model of human physiology. This study makes an effort to reduce the physiological complexity in terms of a composite entropy change.


## 2.1. Entropy Change (Δ
S
B
P: Blood Pressure, Heart Rate, and Skin Temperature)

Let us first consider the human physiological subsystem characterized by blood pressure (BP), heart rate (HR), and skin temperature (ST) as a closed thermodynamic system. For this system, the work done by the system is given by [10](1)δW=P∗dV=BP∗dHR,where P is pressure, BP is blood pressure (mmHg), V is volume, and HR is heart rate (beats/minute).The heart rate is an indirect measure of stroke volume of blood in the heart region and is, thus, used in the place of stroke volume. Please note that heart rate (in beats per minute) is equal to cardiac output (mL/min) divided by stroke volume (mL/beat). The ratio of partial changes in two physiological variables is accompanied by constancy in the third variable, as this is the basis of Maxwell relations. For example, when there is a partial change in heart rate with respect to partial change in entropy, the mean arterial pressure remains constant. As shown in Appendix, the thermodynamic potential is derived as shown starting with Helmholtz’s potential (F
=
U
-
T
S) and arriving with the following property relation: (2)dF=-BP∗dHR-S∗dST.The following Maxwell relation is obtained by invoking the exactness condition to the above property relation (2):(3)∂BP∂STV=∂S∂HRT.The partial derivative is approximated to form a modified Maxwell relation that is used in the present experimental study to compute the physiological entropy change resulting from changes in blood pressure (BP), heart rate (HR), and skin temperature (ST):(4)ΔBPSTHR=ΔSBPΔHRST,ΔSBP=ΔBP×ΔHRΔST.The skin temperature (ST) in the denominator acts as a temperature at the boundary of the closed thermodynamic system and provides a physiological reflection of emotional response. If one were to consider either internal energy or the Gibbs function to define human physiological function, then one would get a negative sign preceding the entropy change (Δ
S
B
P). The deviation from the relaxed state may be positive or negative depending on the imposed external stressor and the internal physiological condition.


## 2.2. Entropy Change (Δ
S
E
M
G: Electromyogram, Electrodermal Response, and Skin Temperature)

Let us consider the human physiological subsystem characterized by electromyogram (EMG), electrodermal response (EDR), and skin temperature (ST) as a closed thermodynamic system. For this system, the work done by the system is given by(5)δW=ℾ∗dL=EMG∗dEDR,where ℾ is tension, EMG is facial electromyogram (mV), L is length, and EDR is electrodermal response (mV).In this regard, the elements of a human physiological subsystem system are considered as equivalent to a conducting wire withℾ as tension and L as the length. The facial electromyogram (EMG) has been used in medical research and clinical practice as a measure of tension in facial muscles [11]. As a result of tension in muscles or nerves, the outer surface of the skin responds by expanding or contracting in terms of changes called electrodermal response (EDR). The skin temperature (ST) represents the thermal response to tension. The thermodynamic potential is again derived, starting with Helmholtz’s potential (F
=
U
-
T
S) arriving at the following property relation: (6)dF=-EMG∗dEDR-S∗dST.The following Maxwell relation is obtained by invoking the exactness condition to the above property relation (2):(7)∂EMG∂STEDR=∂S∂EDRST.The partial derivative is approximated to form a modified Maxwell relation that is used in the present experimental study to compute the physiological entropy change resulting from changes in blood pressure (BP), heart rate (HR), and skin temperature (ST):(8)ΔEMGSTEDR=ΔSEMGΔEDRST,ΔSEMG=ΔEMG×ΔEDRΔST.As mentioned earlier, the use of either internal energy or the Gibbs function would have resulted in a negative sign preceding the entropy change (Δ
S
P). The deviation from the relaxed state may be positive or negative, depending on the imposed external stressor and the internal physiological condition. Specifically, physiological entropy is qualitatively defined as a measure of disorder [12]. It is a kind of global measure that specifies how violent motions and reactions are occurring in nature. Hence, the entropy change in the human physiological system shows the extent of activity within the body as a whole; thus, the entropy change is a significant quantity that characterizes the human body from thermodynamic and holistic (i.e., considering a human body as a whole) viewpoints.The human physiological entropy change could be considered as a composite measure of change in the whole physiological state in response to any external stimuli or stressor. However, if a single physiological indicator alone, such as the mean arterial pressure, can provide that information, then why do we need this entropy change as a composite measure of physiological response? The answer is as follows: the physiological concepts such as stimulus response (SR) specificity, organ response (OR) specificity, individual response specificity, and autonomic balance make the human physiological response a complex phenomenon [11]. The single physiological indicators, in this regard, provide a very narrow representation of the human physiological stress response system. It is only by recognizing the interaction among human subsystems in their response to any stressor stimuli that one could build a better model of human physiology. This study makes an effort to reduce the physiological complexity in terms of a composite entropy change.


## 3. Methodology

### 3.1. Subjects

The data in the study were collected on eight subjects (5 male, 3 female). These subjects completed a standard psychophysiological stress profile procedure routinely used for assessment in the Human Engineering Laboratory at NASA Langley Research Center. The participants were all healthy (i.e., without any major health problems).


### 3.2. Data Collection

The physiological data were collected by a BioPac system via National Instrument LabView. The multiple physiological responses were collected from eight subjects who completed the 78-minutePhysiological Stress Profile and included the two following conditions.Condition 1 (relaxation period).
Subjects relaxed in a semireclining position with eyes open for eighteen minutes listening to a guided relaxation tape. The physiological stress response data were collected. These eighteen-minute time-series data,B
P
O, H
R
O, S
T
O, E
M
G
O, and E
D
R
O, were averaged and used in the physiological entropy change calculations as baseline data points.Condition 2 (task period).
Subjects completed a series of cognitive tasks presented by an oil refinery simulation program calledDexter. The primary goal of this simulation program was to examine the human physiological response to activities that induce boredom, drift in attention, and other negative effects on performance. The task lasted for sixty minutes and two data samples per minute were collected. Each one of these physiological measures—B
P
T
a
s
k, H
R
T
a
s
k, S
T
T
a
s
k, E
M
G
T
a
s
k, and E
D
R
T
a
s
k—was used in the calculations as shown in the next subsection.


### 3.3. Illustrative Example

Let us consider subject #1; the eighteen-minute time-series data is averaged to get the following baseline values. Note thatB
P
P
,
O and V
H
R
,
O are averaged and used in the physiological entropy generation calculations as baseline data points:B
P
O (mmHg): 118.74.
H
R
O (bpm): 68.49.
S
T
O  (°F): 93.14.
E
M
G
O (mV): 334.17.
E
D
R
O (μmho): 557.87.Each one of the physiological measures, B
P
T
a
s
k, H
R
T
a
s
k, S
T
T
a
s
k, E
M
G
T
a
s
k, and E
D
R
T
a
s
k, from the time-series is used to find the entropy change. Let us consider the values at the 5th minute during the task. A sample calculation of physiological entropy change at the 5th minute for subject #1 is shown as follows:(9)ΔSBP=BPTask-BPO∗HRTask-HROSTTask-STO/1.8,ΔSBP=136.05-118.74∗81.97-68.4992.59-93.14/1.8-D=-178.33mmHg·bpm/K. The EMG-based entropy change is calculated as follows:(10)ΔSEMG=EMGTask-EMGO∗EDRTask-EDROSTTask-STO/1.8±D,ΔSEMG=335.25-334.17∗678.93-557.8792.59-93.14/1.8±D=-100.21mV·μmho/K.Note the following:(1)
Divide temperature differenceΔ
S
T by 1.8 to get the units in Kelvin (K).
(2)
The factorD is used to deal with zeros and small temperature difference that would cause huge entropy change numbers. In this case, D
=
-
1 since Δ
S
T is negative. If it is positive then D
=
+
1. This keeps the magnitude of the entropy change to within manageable limits while keeping the direction of temperature change congruent.The above sample calculations are done for the entire task period of sixty minutes for each subject and repeated for all the eight subjects as illustrated in the next section.


## 3.1. Subjects

The data in the study were collected on eight subjects (5 male, 3 female). These subjects completed a standard psychophysiological stress profile procedure routinely used for assessment in the Human Engineering Laboratory at NASA Langley Research Center. The participants were all healthy (i.e., without any major health problems).


## 3.2. Data Collection

The physiological data were collected by a BioPac system via National Instrument LabView. The multiple physiological responses were collected from eight subjects who completed the 78-minutePhysiological Stress Profile and included the two following conditions.Condition 1 (relaxation period).
Subjects relaxed in a semireclining position with eyes open for eighteen minutes listening to a guided relaxation tape. The physiological stress response data were collected. These eighteen-minute time-series data,B
P
O, H
R
O, S
T
O, E
M
G
O, and E
D
R
O, were averaged and used in the physiological entropy change calculations as baseline data points.Condition 2 (task period).
Subjects completed a series of cognitive tasks presented by an oil refinery simulation program calledDexter. The primary goal of this simulation program was to examine the human physiological response to activities that induce boredom, drift in attention, and other negative effects on performance. The task lasted for sixty minutes and two data samples per minute were collected. Each one of these physiological measures—B
P
T
a
s
k, H
R
T
a
s
k, S
T
T
a
s
k, E
M
G
T
a
s
k, and E
D
R
T
a
s
k—was used in the calculations as shown in the next subsection.


## 3.3. Illustrative Example

Let us consider subject #1; the eighteen-minute time-series data is averaged to get the following baseline values. Note thatB
P
P
,
O and V
H
R
,
O are averaged and used in the physiological entropy generation calculations as baseline data points:B
P
O (mmHg): 118.74.
H
R
O (bpm): 68.49.
S
T
O  (°F): 93.14.
E
M
G
O (mV): 334.17.
E
D
R
O (μmho): 557.87.Each one of the physiological measures, B
P
T
a
s
k, H
R
T
a
s
k, S
T
T
a
s
k, E
M
G
T
a
s
k, and E
D
R
T
a
s
k, from the time-series is used to find the entropy change. Let us consider the values at the 5th minute during the task. A sample calculation of physiological entropy change at the 5th minute for subject #1 is shown as follows:(9)ΔSBP=BPTask-BPO∗HRTask-HROSTTask-STO/1.8,ΔSBP=136.05-118.74∗81.97-68.4992.59-93.14/1.8-D=-178.33mmHg·bpm/K. The EMG-based entropy change is calculated as follows:(10)ΔSEMG=EMGTask-EMGO∗EDRTask-EDROSTTask-STO/1.8±D,ΔSEMG=335.25-334.17∗678.93-557.8792.59-93.14/1.8±D=-100.21mV·μmho/K.Note the following:(1)
Divide temperature differenceΔ
S
T by 1.8 to get the units in Kelvin (K).
(2)
The factorD is used to deal with zeros and small temperature difference that would cause huge entropy change numbers. In this case, D
=
-
1 since Δ
S
T is negative. If it is positive then D
=
+
1. This keeps the magnitude of the entropy change to within manageable limits while keeping the direction of temperature change congruent.The above sample calculations are done for the entire task period of sixty minutes for each subject and repeated for all the eight subjects as illustrated in the next section.


## 4. Results and Analysis

The baseline data were obtained by averaging the eighteen-minute relaxation period and are provided in Table1 for the eight subjects. The Δ
S
B
P and Δ
S
E
M
G results are presented in Figures 1 and 2.Table 1
Baseline values from the relaxation period for eight subjects.
Subjects
B
P
O  (mmHg)
H
R
O (bpm)
S
T
O (°F)
EMG
O (mV)
EDR
O (μmho)
1
118.74
68.49
93.14
334.17
557.87
2
120.36
68.30
90.26
903.86
100.86
3
144.37
44.64
83.72
418.81
52.69
4
145.67
45.58
91.34
372.77
40.54
5
135.90
67.85
90.92
897.79
20.00
6
139.43
47.44
87.73
477.24
205.90
7
140.17
73.16
75.64
373.94
381.62
8
122.73
65.70
91.46
886.80
11.08Figure 1
Entropy change (Δ
S
B
P) across eight subjects: 1, 6, and 7 show significantly different patterns of variations compared to others.Figure 2
Entropy change (Δ
S
E
M
G) across eight subjects: 2, 5, and 8 show significantly different patterns of variations while 7 exhibits a unique pattern.Figures1 and 2 quantitatively demonstrate key psychophysiological important concepts such as stimulus response (SR) specificity, organ response (OR) specificity, and individual (IR) response specificity [11]. For each subject, the response to the task varies within the sixty-minute period. For instance, in Figure 1, subject 4 shows normal variability with Δ
S
B
P throughout the task. This is a demonstration of stimulus response (SR) specificity. In Figure 2, the same subject #4 shows greater variability with Δ
S
E
M
G in the beginning of the task while learning to settle down as the task progresses. This is a demonstration of organ response (OR) specificity in which the subject is exhibiting higher reactivity via EMG and EDR combinations. Also, one can observe wide variations between individual subjects in both Figures 1 and 2, which demonstrates individual response (IR) specificity. Furthermore, the average entropy change value over a 60-minute time period was obtained for all eight subjects and shown graphically in Figure 3. It is observed that there is a wide variation among the subjects in terms of average entropy change response to the task demonstrating interindividual differences. These results demonstrate the importance of a thermodynamics-based approach to conducting any medical and/or clinical research involving human subjects.Figure 3
Average entropy change response in eight subjects.It is clear from the results shown in Figure3 that subjects show greater reactivity via entropy change (Δ
S
E
M
G) induced by EMG and EDR than that by entropy change (Δ
S
B
P) induced by blood pressure and heart rate. Both have temperature change as common denominators. Because of higher reactivity, Δ
S
E
M
G involving EMG and EDR measures need to be included in any comprehensive medical stress studies.


## 5. Conclusion

The study was based on the premise that Maxwell relations could be used to compute entropy change in terms of measurable physical variables. The experimental study involved eight subjects and was conducted in the Human Engineering Laboratory at NASA Langley Research Center. The physiological variables were measured at two states, relaxed and stressor task. The average physiological measures taken during relaxation were taken as the baseline values. The physiological changes in blood pressure, heart rate, electromyogram, electrodermal response, and skin temperature from the reference (State 1) to stressor (State 2) were used in the Maxwell relations obtaining entropy changes,Δ
S
B
P, and Δ
S
E
M
G, respectively. The results validate and verify key psychophysiological concepts such as individual response (IR) specificity, organ response (OR) specificity, and stimulus response (IR) specificity.The main purpose of this study was to use a novel thermodynamics-based approach to provide quantitative framework to assess human physiological stress responses. It is clear from the result that the modified Maxwell relations indeed provide something valuable to the biomedical engineering, medical, and health community, a quantitative measure of disorder or stress in human physiology. The entropy change could be postulated as disorder (probability) and incorporated in Shannon’s entropy formulation [13] to link physical and informational aspects of human physiology. As this is beyond the scope of this paper, it could be valuable to the biomedical engineering community to utilize this concept for future applications. This fundamental concept could also lead to development of noninvasive, wearable sensors in order to gain insight into the prevention of human errors and occupational injuries and diseases. Future efforts will target field applications of this methodology. Specifically, the researchers plan to evaluate its usability and merit in occupational settings where workers are exposed to either physical (e.g., extreme temperatures, radiation), chemical (e.g., particulate matter, organics solvents), or biological stressors (e.g., mold, fungi).



---

*Source: 4932710-2016-02-22.xml*