Skip to main content


New concept for quantification of similarity relates entropy and energy of objects: First and Second Law entangled, group behavior of micro black holes expected


When the free energy of similar but distinct molecule-sized objects is plotted against the temperature at which their energy and entropy contributions cancel, a highly significant linear dependence results from which the degree of similarity between the distinctly different members within the group of objects can be quantified and a relationship between energy and entropy is derived. This energy-entropy relationship entirely reflects the mathematical structure of thermodynamic equations, is in this sense fundamental and therefore does probably not dependent on material nor scale. The energy-entropy relationship is likely to be of general interest in molecular biology, population biology, synthetic biology, biophysics, chemical thermodynamics, systems chemistry and physics, most notably in particle physics and cosmology. In physics we predict a consistent and perhaps testable way of classifying micro black holes, to be generated in future Large Hadron Collider experiments, by their gravitational energy and area entropy.


The larger the physical scale is, the less frequently the term 'energy' and the more frequently the term 'entropy' is used in physics discussions. Energy, in the sense of 'bound' or 'inner' energy, is an entity that is usually measured experimentally in some more or less direct way. Entropy is an entity impossible to measure directly; it can only be determined either in conjunction with measured energy and another measured experimental parameter, free energy for instance, or it is calculated or counted using statistical mechanics or some other theory on the degeneracy of microstates. Since, owing to their distance from the observer, very large-scale physical objects are difficult to measure directly, the preferential use of entropy and the Second Law of thermodynamics is not astonishing in cosmology, neither is the preferential use of energy in quantum physics, in particular, strict energy conservation as expressed through the First Law of thermodynamics. Of course both laws apply a priori to all scales and physics, and of course the above statements are not based on statistical analyses or other objective grounds but on the subjective impression of the author to whom correspondence should be addressed.

In this article we present very briefly the results of a comprehensive analysis of published experimental thermodynamic data on the unfolding of many hundreds of proteins and nucleic acids, on molecular associations in host-guest complexes, on the stability of ab initio (quantum mechanically) calculated water clusters and the semi-empirically (force field) calculated formation thermodynamics of small organic molecules from their elements. We then mainly discuss the consequences when i) these numerical results are first grouped into families that distinguish ensembles of evidently similar objects, ii) the grouped results are correlated in a specific two-dimensional projection of a five-dimensional parameter space and, ultimately, iii) the results are detached from the molecular scale.

The discussion begins with deriving an equation that relates energy changes to entropy changes of the same objects without usage of additional empirical parameters or functions that are not explained from the fundamentals. The only new 'entity' or 'information' is the fact that the objects are grouped into families of obviously similar characteristics. Protein mutants and nucleic acid variants are macromolecules that usually differ only very little in overall shape and folding potential - only one or two in dozens or hundreds of 'chain links' are different within the same group - but may differ rather heavily in measured energy and entropy of folding. It is known since 1970 that in many very different chemical and biological systems large entropy and energy contributions compensate one another, to give small resulting free energy changes, that is, small net effects. We do not discuss this here - our studies on the compensation effect and statistical significance of the utilised linear regressions are described in full detail to be published elsewhere - but rather focus on the consequences of the results. Once energy and entropy changes are fundamentally linked to one another, the laws that on the one hand restrict in isolated systems average net energy changes to zero and on the other hand confine spontaneous net entropy changes to zero or more but not less, thus, condemn entropy to maximise over time, may become fundamentally linked as well. If our analysis on the thermodynamics of medium-sized objects, which can either be described by quantum physics or by classical physics, were generalizable to all scales, we were to conclude the following.

The First and Second Law of thermodynamics describe isolated multicomponent systems in the observable universe as objects that conserve their energy due to their very isolation and that spontaneously maximise their entropy over time. For the latter to be true, the objects' size must be sufficiently large for fully reversible changes, that is, exactly reversed changes in their microstates, to become too improbable to occur within their lifetime. Additionally, an isolated ensemble of similar objects in the same universe will spontaneously maximise its overall entropy over time in a way (at a rate) that reflects its overall energy and identity, thus, its compositional and structural characteristics that define it as an ensemble of similar objects. If the physical isolation of the ensemble confines its overall average energy changes to zero, the way (rate) of maximizing entropy can only change when the degree of similarity within the ensemble of objects changes as well. We conclude that, given a constant (accessed) overall volume of an ensemble, the higher the degree of similarity is among its objects the slower is their rate of spontaneous entropy maximization and the closer to maximum entropy they are. Hence, it seems as if the rate of maximizing overall entropy of an ensemble of objects were related to the similarity of what characterises the individual objects within the ensemble.

Here we present a statistical means of quantifying the degree of similarity, namely, through the linear regression coefficient obtained from the correlation of the difference with the ratio of two object characterizing parameters (energy U and entropy S) that both depend on one independent variable (absolute temperature T). We depict, using experimental numerical values, 3D projections of the 5D parameter space {U; S; T; UT·S; U/S}pV (at constant pressure and volume pV).


The vast majority of the primary data are experimental and about one third of those originate from differential scanning calorimetric experiments where both the energy change under constant pressure, i.e., the enthalpy change ΔH, and the position of thermodynamic equilibrium between two macroscopic states, i.e., the free enthalpy change ΔG (Gibbs free energy), are derived from equation 1. The measured heat capacity Cp (at constant pressure) is a function of temperature T within a T-range needed to observe both major macroscopic states (termed 'folded' and 'unfolded') in virtually quantitative abundance. Enthalpy changes in a system open to atmospheric pressure, ΔH = Hmacrostate 1Hmacrostate 2, and energy U in a closed system are linked through U = Hp·V. Likewise, the Gibbs free energy difference ΔG = Gmacrostate 1Gmacrostate 2 is a measure for the driving force towards macroscopic stasis under constant pressure, and free energy is linked through F = Gp·V. The corresponding change in entropy ΔS of the system is usually calculated from ΔG = ΔH ΔS (or ΔF = ΔU ΔS) rather than directly from equation 1.

C p = d H / d T = T · d S / d T = T · ( d 2 G / d T 2 )

Another definition of heat capacity is the mean squared fluctuation in energy scaled by kT 2, or the mean squared fluctuation in entropy scaled by k (the Boltzmann constant), as shown in equation 2 [1].

C p = < δ H 2 > / k T 2 = < δ S 2 > / k

The difference in specific heat capacity between both major macroscopic states is directly measured from ΔCp = Cp(T100% unfolded) – Cp(T100% folded); Δ always refers to the difference between two distinct macroscopic states. Both Cp(100% unfolded) and Cp(100% folded) are assumed to exert the same T-dependence, hence ∂ΔCp/∂T = 0, i.e. ΔCp ≈ const.

The other two thirds of experimental data originate from so-called van't Hoff experiments in which, instead of Cp, equilibrium constant K = (fraction macrostate 1)/(fraction macrostate 2) = exp[–ΔG/RT] (R = 1.9872 cal mol-1 K-1) is measured within an appropriate range of T or other parameter capable of completely shifting the thermodynamic equilibrium from one macroscopic state to another. For thermally induced macrostate changes the accompanying energy and entropy changes are elucidated from fitting the experimental data to equation 3:

R · ln K = Δ H / T + Δ S = Δ G / T

In the vast majority of published van't Hoff experiments heat capacity changes are ignored altogether: ΔCp ≈ 0. This approximation is justified by the usually observed linear relationship for lnK versus 1/T. In both kinds of experiments, calorimetric and van't Hoff, any true T-dependence of ΔCp may be neglected when compared to the one of ΔG = ΔH ΔS (or of ΔF = ΔU ΔS) over the measured T-range. In summary, classical thermodynamics provides us with equations 4 and 5 in the fundamental, most general case ΔCp = f(T) [2]. Equations 6 and 7 result from the 'calorimetric neglection' of the T-dependence of ΔCp. After a 'van't Hoff neglection' of ΔCp, ΔH and ΔS become constants with respect to T.

Δ H T = Δ H T ref + T ref T Δ C p ( T ) d T
Δ S T = Δ S T ref + T ref T Δ C p ( T ) T d T
Δ H T = Δ H T ref + Δ C p · ( T T ref )
Δ S T = Δ S T ref + Δ C p · ln ( T / T ref )


We extracted from the literature 1555 experimental datasets { Δ C p ; Δ H T ref ; Δ S T ref } on the thermal and non-thermal unfolding of proteins and nucleic acids. The vast majority of data was downloaded from the ProTherm database [3, 4] at and controlled in the original literature. For each dataset Tref = TΔH = ΔS= Tm. Tm is the so-called midpoint or equilibrium temperature, the temperature at which in a dynamic and fully reversible two-state equilibrium the fractions of both (two particularly stable and well observable) macrostates are equal, therefore Δ G T m = 0 (eqn. 3). We expanded the above datasets with an additional function each, the state function ΔGT = ΔHT ΔST, using equations 3 (right-hand side), 6 and 7. At that stage, no numerical values were attributed to T yet. Each dataset was now made up of five 'characterizing parameters' { Δ C p ; Δ H T m ; Δ S T m ; T m = Δ H T m / Δ S T m ; Δ G T = Δ H T T Δ S T } , all of which are dependent on one another through the fundamental thermodynamic equations 1 to 5, and of one 'independent variable' T. Note that all five parameters, despite being derived from Cp and T, bear distinct physical meanings (interpretations).

All 1555 datasets were then grouped into 154 families, according to the structural similarity of the members within each group (mostly 'single-chain link' variants, 'point mutants'). The datasets of each of the 154 groups were submitted to a group-specific correlation between the two combined (with respect to ΔH and ΔS) parameters ΔGT and Tm. An increasingly refined sampling of ΔGT on a representative part of the groups led to a complete correlation analysis Δ G T median vs. Tm of all groups at a group-specific T = Tmedian. Tmedian is the statistical median of all equilibrium temperatures Tm of a group.


The correlations between T = 273 and 373 K appeared visibly linear for the vast majority of the analysed groups, hence, a linear regression according to equation 8 was used to characterise every group.

Δ G T = h T T m · s T = h T ( Δ H T m / Δ S T m ) · S T

Detailed results are described in the additional files 1 and 2. Here it suffices to note that all members of the same group share the same 'group parameters' hT and sT which express nothing more than the average energy and, respectively, entropy of the group of similar objects. They are therefore only dependent on T and the choice of which individual members constitute 'a group'. The numerical values for the slope S T median are actually average values of all numerical Δ S T m values of each group member within one group. The numerical values for hT and all other sT depend on ΔCp(T), the more so the larger |TTmedian| is. According to equation 8 the T-dependence of hT and sT is the same as for ΔGT. For ΔCp = const. this T-dependence adopts the form f(T) = a + b·T + c·T· lnT, in which c is nil for ΔCp = 0 (eqns. 3, 6 and 7). We fitted this function to all experimental data, to obtain the 'group constants' (with respect to T) h0-2 and s0-2 for hT = h0 + h1·T + h2·T· lnT and sT = s0 + s1·T + s2·T· lnT. Note that h0-2 and s0-2 [see additional file 2] can all be derived from the Δ S T m , ΔCp and Tm values of a group with no additional information or assumptions (eqn. 42 [see additional file 1]).

The main result is that at Tmedian, at the temperature where the sum of ΔG of all group members within one group is closest to nil, the vast majority of experimental data produces a linearity of unexpected quality. The linearity as such remains visible but its quality, as expressed through the regression coefficient, degrades quite strongly and monotonously with increased |TTmedian| (Figures S14-S15 [see additional file 1]) and, in a non-trivial fashion, as we join evidently less similar objects into the analysed group (Figures S1, S5-S6, S10-S11 [see additional file 1]). The experimental group sizes vary between 4 and 68 (average 10). The regression coefficients r T median of all calorimetric groups lie between 0.90 and 0.999'999 with an abundance maximum between 0.999 and 0.9999 (Figures S12-S13 [see additional file 1]). The van't Hoff groups do not fall far behind (Figure S7 [see additional file 1]). In addition, the same correlation method was tested on the calculated thermodynamics of formation from the pure chemical elements in their standard state of a homologue series of PM3-calculated simple organic molecules, as well as of published ab initio-calculated water clusters [5], using statistical thermodynamics at 298 K. The somewhat lower correlation coefficients r298K as compared to the above experimental r T median values are due to the fact in part that at T = 298 K many calculated datapoints within one group do not center around ΔG = 0. The linearity of similar groups is nevertheless unambiguously apparent (Figures S37-S39 [see additional file 1]).


The mere fact that changes in energy and entropy are fundamentally correlated is not unexpected; after all, their temperature dependence is akin and dictated by the corresponding change in heat capacity (eqn. 1), i.e., their mean fluctuation (eqn. 2). A relationship between free energy and the temperature at which it vanishes is not astonishing either. Both ΔGT and Tm are commonly interpreted as a representation of 'thermodynamic stability', the former is expressed in energy units and depends on ΔCp(T), the latter lends its unit from the temperature scale and is untouched by any T-dependence of ΔCp. However, we were unable to find in the literature any systematic study that would demonstrate this particular linearity from experimental data, nor its strong dependence on the similarity of congeners, nor its highest quality at T = Tmedian. The distinct linear grouping of the theoretically calculated molecules (of chemically very different nature from that of proteins or nucleic acids) is at least inasmuch significant as their thermodynamic parameters are independently derived from partition functions rather than from experimental enthalpies or experimental equilibrium constants, and in spite of the not entirely exact nature of the calculation of S (due to the harmonic oscillation approximation).

Taken together, the similarity-dependent linearity of Δ G T median vs. Tm, quantified through the regression coefficient r T median , seems to be as general as the whole theory of thermodynamics is. It may thus be that this linearity's origin lies at least in part in the mathematical structure of thermodynamics, not entirely in the physics for which thermodynamics was designed to describe. Therefore we proceed with deriving general consequences, with respect to physics, such as the entanglement of the First and Second Laws for groups of similar objects as mentioned in the introduction. We continue with the mathematical and geometrical analysis of a function that was generated from the combination of equations 3, 8 (both right-hand side), 4 and 5 to give through the elimination of ΔGT equations 9 and 10, i.e., the fundamental energy-entropy relationship and mathematical basis for the 5D parameter space { Δ H T m ; Δ S T m ; T m = Δ H T m / Δ S T m ; Δ G T = Δ H T T Δ S T ; T } . Equation 9 is a simplified version for ΔCp = 0 (for clarity) of the general form as shown in equation 10. Both equations can be analytically solved for Δ S T m (eqn. 26 [see additional file 1]).

Δ H T m = T · Δ S T m · h T T + Δ S T m S T + Δ S T m
Δ H T m = T · Δ S T m · h T T m T Δ C p ( T ) d T + T · T m T ( Δ C p ( T ) T ) d T T + Δ S T m S T + Δ S T m

The above functions are variants of the well known quadric x = y·z of the shape of a hyperbolic paraboloid (where x = Δ H T m , y = Δ S T m and z = T), thus, of a single saddle point centered in the origin {x = 0; y = 0; z = 0} and the S4-symmetric function spreading from there with an all-negative Gaussian curvature (Figure 1). Any temperature dependence of ΔCp(T) is consistent with the hyperbolic paraboloid (eqn. 9) as shown in equation 10. For ΔCp = 0 (eqn. 9 with hT = h0 + h1·T and sT = s0 + s1·T from the van't Hoff datasets) the basic shape of the function does not change when compared to x = y·z, although the function area may be quite heavily 'distorted' (not shown). However, for ΔCp ≠ 0 = const. (eqn. 9 with hT = h0 + h1·T + h2·T· lnT and sT = s0 + s1·T + s2·T· lnT) the group constants h0-2 and s0-2 that were obtained from the experimental calorimetric datasets produced shapes of the eyebrow-rising kind. In Figure 2 four views of the same 3D-projection, ΔHT versus ΔST and T, of the thermodynamic 5D parameter space is shown for one particular but representative calorimetrically measured protein mutant group (mutants of Staphylococcal Nuclease). In Figure 3 one to two views of three different 3D-projections for the same mutant group are depicted. Both Figures 2 and 3 focus on the zone that contains the experimental data (yellow dots). The interested reader is welcome to copy any set of experimental group constants h0-2 and s0-2 [additional file 2], plot equation 9 at any scale (best solved for Δ S T m to suppress a maximum of asymptotic planes in certain 3D projections) and enjoy the shapes and wormholes created by the lnT terms. A more comprehensive study on the characteristics of this function shall be published elsewhere.

Figure 1

3D Projections at T = 298 Kelvin of a 5D hyperbolic paraboloid (eqns. 9, 10) where h T = T · s T and Δ C p = 0. Dimensions: ΔH, ΔS, ΔHT·ΔS (= ΔGT), ΔHS (= Tm) and T. The Gaussian curvature K of the hyperbolic paraboloid is negative everywhere: K = – (1 + ΔS2 + ΔH2)-2. Left: Tm = ΔHS. The function Tm = ΔHS is identical in shape as Tm = ΔGTS + T; in these projections the saddle point is at {ΔH = 0; ΔS = 0; Tm = 0}. Right: ΔHS = Tm = ΔH·298 K/(ΔH – ΔG298K). The function Tm = ΔH·T/(ΔH – ΔGT) is identical in shape as Tm = ΔGT·T/(ΔH – ΔGT) + T; in these projections the saddle point is at {ΔH = 0; ΔGT = 298K = 0; Tm = T = 298 K}. The vertical central lines mark ΔH = ΔS = ΔG298K = 0; the lower half of the plots have no physical meaning (quadrants where Tm ≤ 0 Kelvin).

Figure 2

3D-Projections Δ H T versus Δ S T and T of a 5D hyperbolic paraboloid using experimental T -independent Δ C p values. The function is specific of the protein mutant family Staphylococcal Nuclease at pH 7; the primary data were obtained from ProTherm entry numbers 107-120. (a) and (b): Two relatively narrow and orthogonally oriented wormholes in the central region of the quadric. (c) and (d): The smaller wormhole — also visible in (b) at higher temperatures — hosts the experimental data, cf. yellow dots in (d), from which the function was calculated using equations 8 and 9 where hT = h0 + h1·T + h2·T· lnT and sT = s0 + s1·T + s2·T· lnT [see additional file 2]. The narrowness of such wormholes is characteristic for a ubiquitous compensation of ΔHT against ΔST , as briefly mentioned in the introduction, and suggests why empirically good, albeit statistically questionable, 2D-linear relationships are found in a vast majority of experimental ΔHT vs. ΔST correlations in the literature.

Figure 3

3D-Projections other than Δ H T versus Δ S T of a 5D hyperbolic paraboloid using experimental T -independent Δ C p values. The function is specific of the same group as in Fig. 2. Dimensions (for ΔCp = const.): ΔHT, ΔGT, T and Tm. The yellow dots are the experimental datapoints {ΔHTm, Δ G T median } (b), {ΔHTm, Tm} (c) and {Tm, Δ G T median } (d). The yellow line is the linear regression in ΔGT vs. Tm at Tmedian = 320.2 K. With the exception of the projection (d), where wormholes are never found and the quality of the linear regression is best when the datapoints gather around an average zero free energy, the size of the wormholes that harbor datapoints in all other than ΔHT vs. ΔST projections, cf. (b) and (c), precludes any linearity through attempted empirical datapoint correlations. All plots generated by MATHEMATICA® (Wolfram Inc.) and edited in PHOTOSHOP® (Adobe).

The yellow line in Figure 3d, i.e. the experimental isotherm at T = Tmedian, lies in a 'valley' at Tmedian = 320.2 Kelvin created by the saddle of this particular hyperbolic paraboloid. It seems that this isotherm is the best defined of all T, therefore, producing the best linear regression coefficient r T median . Each straight line in ΔGT versus ( Δ H / Δ S ) T Δ G = 0 that represents a structurally similar group is, in geometric terms, a geodesic on the hyperbolic paraboloid. The corresponding group functions Δ H T m ( Δ S T m ) or Δ S T m ( Δ H T m ) , as expressed through equations 9 and 10 are therefore also geodesics. Geometric considerations indicate that the datapoints produce the best rT values in Δ G T median vs. ( Δ H / Δ S ) T Δ G = 0 when they are closest to the maximal negative curvature, thus, to the saddle point of the hyperbolic paraboloid (cf. Figure 3d). Flatter curvatures, thus, steeper surface areas of the hyperbolic paraboloid farther away from the saddle point (cf. Figure 1) allow for a higher dispersal of the datapoints owing to idiosyncratic ΔCp values, which leads to lower regression coefficients rT.

Independently of geometric considerations, we interpret this consistently observed linearity as a (physically) 'minimal expense' or (mathematically) 'minimal action' effect: The appearance or evolution of small structural changes within the same group, i.e., without touching essential framework structuring, can only result in constantly proportional, therefore, unevolving free energy changes being 'linear' with respect to their equilibrium temperature changes. A thermodynamic interpretation of this linear relationship would be that incremental irreversible changes within a group of reversibly dynamic similar but distincty different structures are just as reversible changes are: virtually uncoupled, therefore, additive and independent of the path taken in between, as is the prerequisite for obeying the Gibbs-Helmholtz equation and synonymous to ΔG and ΔF being state functions.

One might argue that the linearity of equation 8 is a simplified manifestation of the Taylor series expansion for any mathematical function f(x) = f(x0) + (df/dx)·(x – x0) + (d2f/dx2)·(x – x0)2 + (d3f/dx3)·(x - x0)3 +... which always becomes approximately linear for any slowly varying function f(x), Δ G T median in this case, sufficiently close to the reference point x0 (Tm or Δ H T M / Δ S T M in this case). In performing the linear correlations ΔGT versus Δ H T M / Δ S T M at T =Tmedian, we do not explicitly claim that the linear relation holds at all temperatures. We do claim, however, that a correlation between ΔGT and Tm at any temperature T using a polynomial of higher than first (linear) degree, as generalised in the above Taylor series expansion, will lead to an analytically solvable relationship for Δ H T m ( Δ S T m ) or Δ S T m ( Δ H T m ) . We did not prove the generality of this claim but solved ΔHT·ΔS = hT – [(ΔHSs1,T + (ΔHS)2·s2,T + (ΔHS)3·s3,T], which is a Taylor series-expanded version of equation 8 (where ΔCp = 0), for ΔH and ΔS, respectively. The expanded nonlinear variants with s3,T = 0 (quadratic) and s3,T ≠ 0 (cubic) did each result in at least one non-complex analytical solution for ΔHS) and ΔSH), albeit bearing a more complicated mathematical structure (not shown). In other words, we claim that a fundamental relationship between energy and entropy for a group of similar objects results from any analytically solvable relationship between ΔGT and Δ H T M / Δ S T M . We opt for the simplest, a linear solution: ΔGT and Δ H T M / Δ S T M are proportional over a reasonably large temperature range.

Most important for physics is the fact that group specific thermodynamic parameter spaces depict the only possible values that can be realised by a particular group of similar objects. The rest is void, terra incognita for the group members, unless an object changes its characteristics (structure, composition, etc.), unless it 'dissimilarises' off from 'its' group - most likely, to join some other one. The definition of a group, that is, how to determine whether a number of individuals belong to the same group or not, seems at first sight worrying or at least not clearly solved. However, when we think of individuals as being more or less similar to one another, we see that a clear distinction between different groups is not a fundamental issue. Similarity does exist; in the microscopic and macroscopic world it is often a matter of judgement according to some objective, statistically relevant technical signal (at highest available resolution) or at least a subjective physiological 'measurement' ("I know it when I see it", cf. Graphical Abstract). For microscopic objects such as molecules, one should never be tempted to define a group through a good linear regression coefficient only; independent knowledge and/or studies are mandatory. For instance, the advantage of studying mutant protein families not only means being able to analyse a large number of families and sometimes many congeners within one family. Most importantly, we are also certain that single or even multiple site mutants of the same protein do indeed belong to the same structural group, the mutants are undoubtedly similar to one another. Other molecular systems such as synthetic host-guest complexes or water clusters may be less evident to this respect. Still other objects might be even more readily grouped than mutant proteins (cf. Conclusion). The concept of similarity is intrinsically a not readily quantifyable one because intuitively it seems to be a not very objective 'measurement', at least down to Planckian scales: How similar and with respect to what exactly?

We are free to group similar objects essentially at will. For example, we can group one set of RNA hairpins into two families, the one that bears various all-Watson-Crick pairs and the one that contains various single-mismatched base pairs at different positions in the stem, the stem length and loop sequence being the same in both families [6]. We can overlook this subtle difference and treat those hairpins as one group that consist of the same loop sequence and stem length irrespective of single mismatches being present or absent in the stem. The outcome will be a slightly lower linear regression coefficient for this group. It can then be compared to another group of RNA hairpins showing, for example, the same stem length and stem sequence variations but a different loop sequence. We can treat protein mutant families with the same varied degrees of precision/resolution. We could define all known proteins as belonging to the same group and compare it to a more drastically different group of compounds (objects). Nothing prevents us from grouping objects at still lower resolution; the obvious trade-off will be increasingly lower linear regression coefficients. As a matter of fact, there is no a priori objection that we can think of to the grouping of the entire universe and comparing it to some other one, if it were observable. In principle, one would have to agree upon a set of observables (like energy, entropy and temperature), measure them on a statistically representative number of individual members of what we decide, through some hopefully objective criterium, to call a group, determine the corresponding group parameters and then gain easier access to more members of the same group but also, to obtain an objective means for the comparison of this group to another one. In practise, of course, as we embrace more and more dissimilar objects, we will probably evoke increasingly unacceptable linear regression coefficients. Where this limit of a meaningful group analysis lies remains to be seen.


In this study we introduce a geometrical parameter space description of thermodynamics and offer a general way of objectively quantifying similarity (to whatever resolution) of individual objects based on two well known abstract notions (not postulated 'empirical' physical parameters): the use of the knowledge of a group membership, and the mathematical relationship between difference and ratio being the results from the two most fundamental mathematical operations, substraction and, respectively, division. The latter notion opens access to a higher than three-dimensional (ΔH, ΔS, T) geometrical description of thermodynamics through expansion of the parameter space with ΔHT·ΔS and ΔHS. The combination of both notions indicates a group-related redundancy in the mathematical structure of thermodynamics; a redundancy which becomes evident when relating substraction and division for the characterisation of similar objects. This redundancy necessarily unravels a group-related fundamental relationship between energy and entropy for similar objects and, possibly, a general unified law of thermodynamics for structured matter. According to our findings, any group of similar objects may be characterised by precisely how the energy and entropy of each individual group member is related (coupled) to one another. We show that similar dynamic structures, for example molecules, 'minimise their action' on thermodynamic state changes such that, within a structural framework — within 'a group' as specified by the group parameters hT and sT using equations 8, 9 and 10 — the distinction between energy and entropy becomes a formal one.

The usually incomplete knowledge of all molecular properties of a thermodynamic system, such as differential solvation, salt, and bulk solvent effects in biomolecular systems, continues to confront us with the limitation of exactly calculating the free energy, the enthalpy, or the entropy from the fundamentals. However, having at hand reliable experimental or theoretical data of both ΔG and ΔH of as many group members of similar structures as possible, thus, of a statistically sufficient number of group members, we can predict from either ΔH or ΔG of more group members their respective ΔG or ΔH and concurrently ΔS. The relatively simple mathematical structure of group thermodynamics allows us to quantify through linear regressions the structural similarity imprinted into the thermodynamic behavior of, in principle, any structural framework. On a molecular scale, group thermodynamics may strongly simplify the elucidation of entropies of molecules that are known to belong to a group of similar compounds through a bypass of costly calculations of the vibrational components of idealised partition functions. With the knowledge of the group parameters hT and sT at hand, S can be calculated from U or H. In addition, it may be a possibly useful complement for cross-checking ΔG calculations that have been obtained from simulations using molecular dynamics techniques. Generally group thermodynamics may contribute to systematic analyses in biomolecular and chemical thermodynamics and, when applied to chemical reaction kinetics, in systems chemistry.

Theories from quite different domains such as, to name a few, probability theory [710], information theory and the emergence of complex systems [1118], quantum relativity/cosmology [1929] and string theory [30] operate with entropy and the Second Law of thermodynamics yet in conjunction with parameters different from the ones studied here. Urgent problems are being at least attacked, and possibly solved, through the insight into apparent and/or fundamental analogies between statistical thermodynamics and, for example (respectively), randomness of sequential irregularities ("algorithmic entropy", "approximate entropy"), computational compactness ("logical depth"), quality change of hereditary information (change in systemic "knowledge" through periodically discarded "Shannon entropy"), the dynamics of black holes ("Bekenstein-Hawking entropy"), and tracing back the microscopic origin of their area-entropy by counting the degeneracy of periodic and persistent topological defects (Bogomol'nyi-Prasad-Sommerfield soliton bound states) in certain kinds of supersymmetric branes that mimic the thermodynamics of idealised extremal, highly charged black holes. In all above cases the problem arises of how to reliably quantify or sample randomness, logical depth, knowledge, entropy, in order to understand their physical origins and perhaps their development over time. The energy-entropy relationship derived from thermodynamic group characteristics may help solve one or the other problem, in particular, when the to be analysed physical objects are not as potentially overwhelmingly dissimilar as chemical systems can be — in order to ease, for a start, the choice of groups.

Black holes, being the most immensely dense and, with respect to their composition, the perhaps most uniform objects known in physics, are all in a state of maximal entropy and are thought to differ from one another through, out of all known matter, the least of characterising parameters; only mass, angular momentum and, for some limited time period, electric charge makes them different: "black holes have no hair". In contrast, elementary particles may differ through a whole plethora of characteristics (according to the standard model) and the variability, thus, potential dissimilarity of objects that are composed of these elementary particles (of 'normal' nonrelativistic matter) multiplies, i.e., increases at a geometric rate with the number of involved particles. If micro black holes indeed existed and could be transiently generated in future Large Hadron Collider experiments, if different classes of such potentially highly similar objects could be observed and analysed, we would predict that the relationship between their gravitational energy and the surface area of their event horizon would correlate in a fashion that were characteristic for their kind: Energy (= mass) and entropy (= surface) would correlate, through equation 10, differently, i.e., with different group parameters for objects of a particular (range of) angular momentum and electric charge than for another. Distinct groups should appear and be best visible in free energy correlations as formulated in equation 8. A difficulty might arise from the fact that micro black holes are not expected to be formed in a thermodynamic equilibrium, but rather 'kinetically controlled'. How then to measure free energy? We imagine that a measure of free energy of micro black holes would be their abundance under given experimental conditions: Plot under maximum and constant total abundance ('steady state') conditions the logarithm of abundance (through counting) versus ratio of gravitational energy (mass) over surface (of the event horizon). The linearity should produce the best linear regression coefficients when, within a group of analysed micro black holes, the median mass is populated most.


  1. 1.

    Prabhu NV, Sharp K: Heat capacity in proteins. Annu Rev Phys Chem 2005, 56: 521–48. 10.1146/annurev.physchem.56.092503.141202

  2. 2.

    Benzinger TH: Thermodynamics, chemical reactions and molecular biology. Nature 1971, 229: 100–2. 10.1038/229100a0

  3. 3.

    Bava KA, Gromiha MM, Uedaira H, Kitajima K, Sarai A: ProTherm, version 4.0: thermodynamic database for proteins and mutants. Nucleic Acids Res 2004, 32: D120–21. 10.1093/nar/gkh082

  4. 4.

    Kumar MD, Bava KA, Gromiha MM, Prabakaran P, Kitajima K, Uedaira H, Sarai A: ProTherm and ProNIT: thermodynamic databases for proteins and protein-nucleic acid interactions. Nucleic Acids Res 2006, 34: D204–6. 10.1093/nar/gkj103

  5. 5.

    Dunn ME, Pokon EK, Shields GC: Thermodynamics of Forming Water Clusters at Various Temperatures and Pressures by Gaussian-2, Gaussian-3, Complete Basis Set-QB3, and Complete Basis Set-APNO Model Chemistries; Implications for Atmospheric Chemistry. J Am Chem Soc 2004, 26: 2647–53. 10.1021/ja038928p

  6. 6.

    Strazewski P: Thermodynamic Correlation Analysis: Hydration and Perturbation Sensitivity of RNA Secondary Structures. J Am Chem Soc 2002, 124: 3546–54. 10.1021/ja016131x

  7. 7.

    Chaitin GJ: Randomness in arithmetic. Sci Am 1988, 259: 80–5. 10.1038/scientificamerican0788-80

  8. 8.

    Pincus SM: Approximate entropy as a measure of system complexity. Proc Natl Acad Sci USA 1991, 88: 2297–301. 10.1073/pnas.88.6.2297

  9. 9.

    Pincus S, Singer BH: Randomness and degrees of irregularity. Proc Natl Acad Sci USA 1996, 93: 2083–88. 10.1073/pnas.93.5.2083

  10. 10.

    Pincus SM, Kalman RE: Irregularity, volatility, risk, and financial market time series. Proc Natl Acad Sci USA 1997, 101: 13709–14. 10.1073/pnas.0405168101

  11. 11.

    Kuhn H: Model Consideration for the Origin of Life. Naturwissenschaften 1976, 63: 68–80. 10.1007/BF00622405

  12. 12.

    Bennett CH: On the nature and origin of complexity in discrete, homogeneous, locally-interacting systems. Found Phys 1986, 16: 585–92. 10.1007/BF01886523

  13. 13.

    Bennett CH: Information, Dissipation, and the Definition of Organization. In Emerging Syntheses in Science. Edited by: Pines D. Addison-Wesley, Massachusetts; 1987:297.

  14. 14.

    Kuhn H: Origin of life and physics: Diversified microstructure - Inducement to form information-carrying and knowledge-accumulating systems. IBM J Res Devel 1988, 32: 37–46. 10.1147/rd.321.0037

  15. 15.

    Lloyd S, Pagels H: Complexity as Thermodynamic Depth. Ann Phys 1988, 188: 186–213. 10.1016/0003-4916(88)90094-2

  16. 16.

    Landauer R: A simple measure of complexity. Nature 1988, 336: 306–7. 10.1038/336306a0

  17. 17.

    Kuhn H: Origin of life - Symmetry breaking in the universe: Emergence of homochirality. Curr Op Colloid Interface Sci 2008, 13: 3–11. 10.1016/j.cocis.2007.08.008

  18. 18.

    Kuhn H: Is the transition from chemistry to biology a mystery? J Syst Chem 2010, 1: 3. 10.1186/1759-2208-1-3

  19. 19.

    Christodolou D: Reversible and irreversible transformations in black-hole physics. Phys Rev Lett 1970, 25: 1596–97. 10.1103/PhysRevLett.25.1596

  20. 20.

    Christodolou D, Ruffini R: Reversible transformations of a charged black hole. Phys Rev 1971, D4: 3552–55. 10.1103/PhysRevD.4.3552

  21. 21.

    Penrose R, Floyd R: Extraction of rotational energy from a black hole. Nature Phys Sci 1971, 229: 177–9.

  22. 22.

    Hawking SW: Gravitational radiation from colliding black holes. Phys Rev Lett 1971, 26: 1344–6. 10.1103/PhysRevLett.26.1344

  23. 23.

    Bekenstein JD: Black holes and the second law. Nuovo Cimento Lett 1972, 4: 737–40. 10.1007/BF02757029

  24. 24.

    Bekenstein JD: Black holes and entropy. Phys Rev 1973, D7: 2333–46. 10.1103/PhysRevD.7.2333

  25. 25.

    Bekenstein JD: Generalized second law of thermodynamics in black-hole physics. Phys Rev 1974, D9: 3292–300. 10.1103/PhysRevD.9.3292

  26. 26.

    Carter B: Rigidity of a black hole. Nature 1972, 238: 71–2. 10.1038/238098b0

  27. 27.

    Bardeen J, Carter B, Hawking S: The four laws of black hole mechanics. Comm Math Phys 1973, 31: 161–70. 10.1007/BF01645742

  28. 28.

    Hawking SW: Black hole explosions? Nature 1974, 248: 30–1. 10.1038/248030a0

  29. 29.

    Hawking SW: Particle creation by black holes. Comm Math Phys 1975, 43: 199–220. 10.1007/BF02345020

  30. 30.

    Strominger A, Vafa C: Microscopic origin of the Bekenstein-Hawking entropy. Phys Lett B 1996, 379: 99–104. [] 10.1016/0370-2693(96)00345-0

Download references


We thank Prof. Peter Schuster, Theoretische Chemie, Universität Wien, Prof. Emmerich Wilhelm, Physikalische Chemie, Universität Wien, and Prof. Irene Poli, Statistical Department, University Cà Foscari, Venezia, for critically reading an extended version of the manuscript, and Prof. Günter von Kiedrowski, Bioorganische Chemie, Ruhr-Universität Bochum, for critically reading many versions of the manuscript and important enlightening discussions about a Unified Law of Thermodynamics. We are indepted to Prof. Bertrand "BOP" Castro (ex Sanofi-Aventis, Gentilly), for calculating the formation thermodynamics of simple organic homologues, and to Prof. Hans-Christoph Im Hof, Mathematical Institute, University of Basel, for performing a differential geometry analysis on the Gaussian curvature and geodesics of x = y·z. A preliminary version of this manuscript was posted to on 15th June 2009. Last but not least we greatly acknowledge the European Cooperation in Science and Technology for their pioneering, ongoing and generous support of Systems Chemistry, in particular, through the COST Action CM0703, as well as the European Science Foundation for their support in divulging the contents of this recently constituted research community and

Author information

Correspondence to Peter Strazewski.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

PZ derived the Mathematical Appendix [see additional file 1] and contributed significantly to the correct description of the mathematical relationships (in particular, eqn. 10 and T-dependence of Cp, hT and sT) and much of the fundamental physics in the text. ST extracted all primary data from the ProTherm and ProNIT databases at, cross-checked the numerical values and analysed all error margins in the original literature, carried out [see additional file 2] and plotted all linear regressions and polynomial fittings (Figures S2, S3, S4, S8, S9 [see additional file 1]). PS derived equations 8, 9 and Δ S T m ( Δ H T m ) as shown in equation 26 [see additional file 1], conceived of the study and wrote the manuscript and both additional files. All authors read and approved the final manuscript and both additional files.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Zimak, P., Terenzi, S. & Strazewski, P. New concept for quantification of similarity relates entropy and energy of objects: First and Second Law entangled, group behavior of micro black holes expected. J Syst Chem 1, 2 (2010).

Download citation


  • Entropy
  • Black Hole
  • Water Cluster
  • Similar Object
  • Linear Regression Coefficient