Home Page Overview Site Map Index Appendix Illustration About Contact Update FAQ


Thermodynamics


Non-equilibrium Thermodynamics & Life (2020)

Equilibrium Thermodynamics +
       Entropy Dissipation in Irreversible Process
Local Equilibrium Thermodynamics
Steady State
Dissipative Processes, Structures and Life
Entropy and Information
Self-assembly and Self-organization + (2024 update)
Origin of Life

Equilibrium Thermodynamics + Irreversibility

For example, the variation of the internal energy U can be expressed as the difference between the initial (i) and finial (f) states :
U = TS - PV + N ----- (2)
where U = Uf - Ui, ... Variation of heat and work are defined by Q = TS and W = PV respectively. Note that the changes are induced by the extensive variables.

A more specific example is provided by the idealized Carnot Engine Cycle (Figure 02,a): See "Statistical Physics of Self-replication" for derivation of Eq.(2a).

[Top]


Local Equilibrium Thermodynamics

Local Equilibrium Thermodynamics Diffusion In Local Equilibrium Thermodynamics, the whole system is not in equilibrium; but the intensive variables such as T, P, ... can be defined within a small volume, i.e., they are homogeneous (well-mixed) inside that volume (Figure 04,a). In addition, if the system can be described by mathematical formulas, then these intensive variables have to be differentiable, i.e., the functions should be smooth, continuous, and no break, angle, or cusp so that the derivative exists across the

Figure 04 Local Equilibrium Thermodynamics

Figure 05 Diffusion [view large image]

whole system (Figure 04,b). Accordingly, the so-called Gibbs equation in the form of Eq.2 can be carried over with the difference "" replaced by the differential "d" (see Figure 04,b) : dU = TdS - PdV + dN ----- (3).
A simple example is provided by the process of diffusion, where dU = 0, and dV = 0. In terms of entropy per unit volume s, and number density n, the rate of change of s can be derived from Eq.(3) by considering the x direction only without loss of generality :

T(x)[ds(x)/dt] = - (n,x)[dn(x)/dt] ----- (4).

By the equation of continuity, the variation in time can be related to the flux J crossing a surface (see "Fluid Dynamics") :
n/t = - (J/x) ----- (5),
where "" is the partial derivative similar to "d" but is specific for a particular variable in case there are many.

An additional assumption of local equilibrium thermodynamics is that the flux (see Onsager reciprocal relations below)
J = L(/x) ----- (6)
where L is a phenomenological coefficient (in #density-cm2/erg-sec). Finally, we obtain the diffusion equation :

----- (7),
where D = L(/n) (in dimension cm2/sec).
A more general formalism considers the intensive variable as some sort of thermodynamic force Fi, which can be or T or P, ... It has been shown that in the regime where both the flux Ji are small and the thermodynamic forces vary slowly, the rate of creation of entropy is linearly related to the fluxes :

.

The second law of thermodynamics requires that the matrix L be positive definite. Statistical mechanics considerations involving microscopic reversibility of dynamics imply that L is also symmetric, i.e., Li,j = Lj,i. It is hence known as Onsager reciprocal relations.

Applying Eq.(8) to the example of diffusion yields : s/t = L(/x)2 = [L(/n)2](n/x)2 0 ----- (9),
while it looks different from Eq.(7), the dimension is correct at least.

[Top]


Steady State

Molecules and atoms in solid state can form structure of certain shape that can last for long time. The stability is sustained by force which holds them together. In liquid or gaseous state, energy infusion from external environment melts or evaporates the solid into form-less state (suitable for fluid-dynamics and/or thermodynamic treatment). However, in some circumstances, the energy or entropy flow in such a way that some sort of structure re-appeared. The Rayleigh-Benard convection provides a good example to demonstrate how may such structure be formed.

Rayleigh-Benard Convection A Rayleigh-Bénard (R-B) system is created by heating a thin film of viscous liquid (oil, for example) from the bottom. For certain threshold of heat infusion, convective cells of regular form appear in the system (Figure 06,a). The formation arises via the upwardly buoyant force against the opposing viscous resistance as shown mathematically by the equation of motion in fluid dynamics (see Figure 06,b for configuration and pattern of convective flow).

Figure 06 Rayleigh-Benard Convection [view large image]

    There are 2 ways to calculate the entropy production in the system :

  1. The variation of entropy S = Sp - SE, where SE is the total entropy infusion from the external environment, while Sp is the total entropy generating within the system. For steady state S = 0, and thus,
    Sp = SE = dQ/T ----- (14),
    where the heat flow to the system dQ is the difference between the inflow from the heating surface and the outflow to the inner surface of the container (in contact with the liquid). The integration is over the entire surface and duration.

    The conductive heat transfer is : dQ/dAdt = -k(Tq - T)/x, where A is the cross-section area, t the time, x the thickness of the conducting plate, Tq the temperature of the heat source, and k the conductivity in erg/cm-sec-Ko.

  2. Another method is to adopt the formula in Eq.(8) for local equilibrium thermodynamics.
Entropy Spatial Distribution
    It turns out that the 2 methods to calculate entropy production are essentially in agreement for wide range of Ra (Figure 09,b) except:

  • When the experiment (numerically in Figure 09) turns on initially, there is substantial difference until the system stabilizes to steady flow. The effect is attributed to the assumption of local equilibrium in computing Sg. Such assumption is not valid in a complex, unsteady state (see Figure 09,a).

Figure 08 Entropy Spatial Distribution, Ra = 20100 [view large image]

Figure 08 shows the entropy spatial distribution from two-vortex numerical simulation with Ra = 20100. The minima occur at the locations of the vertexes, i.e., there's more orderly formation at those spots.
Rayleigh Number
  • Figure 09,b shows increasing entropy production as the value of Ra goes up. The insert displays the entropy production difference (Sp - Sg), which increases with Ra. The effect is also attributed to the inadequate assumption of local equilibrium.
  • Figure 09 Benard Cell Entropy Production [view large image]

    See "The Character of Entropy Production in Rayleigh–Bénard Convection".

    [Top]


    Dissipative Processes, Structures and Life

    Table 01 lists the dissipative structures in steady state over the eons as the energy flow per unit mass increases (see Figure 11,a).

    System Source m Energy Flow (erg/s-g) Comments
    First Galaxy GM/R 1.5 M = 1012Msun, R = 10 kpc, contraction time ~ 300 My
    Sun (ongoing) Solar Constant 2 m = (Solar constant) x (4D2/Msun), D = 1 AU
    Earth GM/R 200 Accretion time of planetesimals ~ 100 yrs, see "Origin of Solar System"
    Small Plant (ongoing) Solar Constant 103 100 cm2 leaf area, mass ~ 100 gm; useful light ~ 30%, efficiency ~ 3%
    Human (ongoing) ~ 300 watts 5x104 Mass = normal weight per person ~ 62 kg
    Human Brain (ongoing) 20 watts 2x105 Human Brain weight ~ 1 kg
    2000-watt Society 2000-watts 3x105 Mass = normal weight per person ~ 62 kg
    Cyanobacteria (ongoing) Solar Constant 109 For density ~ 1 g/cm3, m ~ Flux (A/V) x 0.01(efficiency), r ~ 10-4 cm

    Table 01 Dissipative Structures

    Constants and unit conversions : Solar Constant = 1.36 kW/m2, watt = joule/sec = 107 erg/sec, 1 day = 86400 sec, 1 year ~ 3x107 sec, Gravitational constant G = 6.67x10-8 cm3/sec2-gm, Msun = 2x1033 gm, Rsun = 1.3x107 cm, 1 AU = 1.5x1013 cm, and 1 pc = 3x1018 cm.
    BTW, the ultimate source of energy from the Sun is the thermo-nuclear fusion at its core.

    Europa That brings up some ideas about life in the universe. It seems that unicellular organisms such as the cyanobacteria are the earliest life form that can adjust to environmental stress readily because of simpler body plan but enough bio-functions to respond. They can live in all sorts of extremely harsh condition even in other worlds such as Europa (Figure 13f). But they would never be able to send out radio signals to announce their presence. It is no wonder that the current search for Extra-Terrestrial Life by big radio telescopes have turned up nothing so far. The chance of detecting human-like creatures is very slim as we are the product of environmental changes under sequence of unique circumstances. Statistically, it can never be reproduced (see "Life on Mars").

    Figure 13f Europa
    [view large image]

    Table 02 below lists some major events about life on Earth to show the process and duration for its development with pictorial illustration in Figure 13g.

    Event Time (GYA) Interval (GYA) Comments
    Big Bang 13.8 ~ 0 Beginning of the Universe
    Solar System 4.5 9.3 Formation of the Earth
    LUCA, ancestor of prokaryotes 3.6 0.9 Hydrosphere + Chemical reactions + RNA world (see "Prebiotic")
    Cyanobacteria 2.4 1.2 Mutations + Proliferation + Photosynthesis Oxygenic atmosphere
    Eukaryotes 1.6 0.8 Novelty by mixing parts or whole from prokaryotes
    Cambrian Explosion 0.54 1.06 From uni- to multi-cellular, cause not clear (see Figure 13g)
    Present 0 0.54 Darwinian evolution with adoption to changing environment

    Table 02 Major Events of Life on Earth

                                        

    Figure 13g Ascending Complexity of Life on Earth [view large image]

    [Top]


    Entropy and Information

    When expressing thermodynamics in terms of statistical mechanics, the temperature T, pressure P, ... are considered to be macro-states. There could be many micro-states within each macro-state as long as they don't change the property of that macro-state. For exmaple in Figure 13h, there are 8 particles distinguishable by their color. The distribution of particle in each of the 2 partitions defines the macro-state A, B, C, D. The number of arrangements can be calculated according to the formula shown in the Figure. This example illustrates different significance conveyed by entropy and information because they are related to differently.

    In term of , entropy is defined as : S = kBlog2() ----- (17).
    Entropy and Information As shown in Figure 13h, the macro-state A has the lowest entropy with S = 0 for = 1. This is a very special state that need special arrangement to reach and is considered to be very orderly. While the macro-state D is the most commonly occurrence at equilibrium, there's a lot of micro-states in such configuration.

    Information is defined by : I = log2(1/) ----- (18).

    It signifies that one particular micro-state is selected out of so many others such as = 70 in macro-state D, thus conveying a very specific (a lot of) information. In other words, information creates

    Figure 13h Entropy and Information [view large image]

    order out of chaos. This is the kind of arrangement suitable for the formation of DNA molecules. It has to be a very specific arrangement out of the 5 nitrogen-bases for coding a very specific gene. Both processes occur in the life cycle as shown in Figure 02d.

    [Top]


    Self-assembly and Self-organization + (2024 update)

    As shown in Figure 14a,a, self-assembly extends the scope of chemistry to weaker interactions, thus encompassing more building blocks including those in nano- and other length scales. Self-assemble of parts into whole can be explained by minimizing the Gibbs free energy (Figure 14a,b). The dissipated energy is absorbed by the external environment. Water is an exceptionally good medium to remove such
    Information and Self-assembly energy and thus reducing the entropy inside. That's why many scenarios for the "origin of life" involves sea water either in tidal pool or hydrothermal mounds. Self-assembled structures can be thermodynamically very stable depending on the depth of the energy valley. The ossified bony structure in many fossils stands witness to its durability (Figure 14a,c).

    Figure 14a Minimum Energy, Self-assembly [view large image]

    See "OoLife, Self-assembly" and "virus" for further details.


    [2024 Update]

    The research article "Assembly theory explains and quantifies selection and evolution" published in October 2023 proposes a theory to link physics with biology. It is embodied in a single equation shown in Figure 14b. An example of assemble a protein is illustrated in Figure 14c. It shows 4 amino acids as building blocks (left); the actual object and its representation as a string (middle); and the assembly pathway to construct the proteins (right). Its relationship to the mathematical symbols :

    ni - it is the number of connecting lines (with an arrow head) from various boxes to form an object,
      for example :
    1. n1 = 4 to form HHH.
    2. n2 = 4 to form GGG.
    3. n3 = 2 to form HHHV.
    4. n4 = 2 to form GGGF.
    5. n5 = 2 to form HHHVG.
    6. n6 = 2 to form GGGFHHHVG.
    Assembly Equation Assembly Example Thus NT = 10, since NT = (ni - 1), where N = 6.
    The example assumes the assembly index ai = 8 for all i's. The positive ai implies that self-assembly is favorable. Self-assembly is un-favorable with ai < 0; for ai = 0, the assembly is neutral, i.e., A = 1.

    Figure 14b Assembly Equation
    [view large image]

    Figure 14c Assembly Example
    [view large image]

    See other examples in Figure 14a


    This new idea about self-assembly provides a mathematical frame work as shown in the Assembly Equation (Figure 14b) for the Effective Theory below :

    Some Reductionist's Systems A popular concept of the "Theory of Everything" considers such theory to consist the most basic ingredient upon which everything can be built or understood. The most often definition of the ingredient is on the size of the components. Unfortunately, this is not the case as the construction become very complicated very soon when it is used to describe the system in the next (larger) level.

    Figure 14d Reductionist's Systems [view large image]

    Size of Some Systems The reductionist's approach is to use different method/theory for each different level and merge the various levels together at the end. Usually the connection is to adopt the derived quantities or properties at the lower (smaller size) level as input parameters or established facts in the theory for the higher (larger size) level. Details in the lower level can be neglected by the theories for higher levels. Figure 14d and Table 03 illustrate some of the systems according to the size of their components. The classification is incomplete as there are many more systems below, above, and in between. The displays only try to emphasize the place of the world of nano-scale among the other systems. The image on the left shows finer size division for some systems and components in term of the viewing apparatus.

    See also the mathematical intricacy of the many body system either in classical or quantum theory.


    System System Size (cm) Component Size (cm) Interaction(s) Effective Theory
    Atom ~ 10-8 10-16(electron) - 10-12(nucleus) Electromagnetism Quantum Mechanics
    Cell 10-5 - 10-1 10-7(Water Molecule) - 10-6(Protein) Covalent and Hydrogen Bonds Molecular Biology
    Living World 109(Earth) 10-5(cell) - 3x103(blue whale) Light, Sound, Smell, Force Biology
    Solar System 5x1017 108(Pluto) - 4x1014(Sun) Gravity Newtonian Mechanics

    Table 03 A Reductionist's View of Some Systems

    Thus, according to the Effective Theory, the assembly paths in the Figure 14,c example can be interpreted as :
  • The symbols for H, G, V, F are the lower level components in the form of amino acids.
  • The connecting line represents the interacting force between the amino acids or the intermediate products.
  • The combinations are the proteins at the next size level.
  • The assembly index ai is the parameter linking the lower levels to the intermediate/final products. It is usually determined by experiment/measurement.

  • Here's an example which can calculate the assembly index (see Figure 14b) from theory in principle. It is the Quark-Proton transition occurred about 10-6 sec after the Big Bang at energy ~ 1 Gev (see "Cosmic Time Line"). The building blocks are 3 quarks with different color charges, i.e., n1 = 3 with N = 1 (see Figure 14e).
    Quark-Proton Transition The assembly index a1 ~ cross-section (1 Gev) / cross-section (threshold). The cross-sections for the reactions are influenced by the energy density, temperature, and other factors. Theoretical models, such as those based on quantum chromodynamics (QCD), are used to describe these interactions. However, obtaining precise and accurate cross-section values for the relevant processes is a challenging task, and it involves sophisticated mathematical and computational methods. Researchers typically use computer simulations and numerical techniques to model the early universe and the processes occurring during that time. These simulations take into account the fundamental forces and particles involved, incorporating the principles of quantum mechanics and relativistic physics.

    Figure 14e Quark-Proton Transition [view large image]

    [End of 2024 Update]

    Self-organization also assembles parts into whole; however, it is very different from self-assembly which ends up in an equilibrium state. Off-equilibrium processing is the main feature of self-organization. See "OoLife, Self-reorganization" and "The Fundament of Cell Biology" for further details.

    A very specific example so essential to life is the process of exciting the ground state of carbon to the SP3 state upon the infusion of 2 ev energy (Figure 15,b). It promotes the carbon atoms into a tetrahedral structure, which provides stable covalent bonds with other atoms. This property of carbon accounts for the large number of known compounds. At least 80 percent of the 5 million chemical compounds registered as of the early 1980s contain carbon. The affinity of carbon for the most diverse elements does not differ very greatly - so that even the most diverse derivatives need not varying very much in energy content.
    SP3 State of Carbon This ability allows the organic world to exist in a special form of thermodynamic stability for as long as the Sun keeps supplying the 2-ev energy. This is the same process that creates all the parts for organizing into a living cell as shown in Figure 15,c, and the basic reason for pumping energy into biological system to maintain metabolism and cellular structure.

    Figure 15 SP3 State of Carbon, Self-organization

    SP<sup>3</sup> State of Glucose The stand alone SP3 excited state has a half-life about 10-15 sec, the tetrahedra would break down and return to equilibrium (see Figure 02d) in absence of energy infusion in the form of Sun light to have it recycled (see "Photosynthesis - Respiration"). The stability of the SP3 state in organic molecules like carbohydrates (glucose) is maintained by the strength of the covalent bonds and is not subject to spontaneous de-excitation for the individual carbon atom (Figure 16). The energy stored in SP3 is released with the break up of the covalent bonds.

    Figure 16 SP3 State in Glucose

    The 2 ev energy released by de-excitation of SP3 is distributed in 6 ATP ~ 0.32 ev each.
    ChatGPT provides further details (in italic) :

    In general, during cellular respiration, glucose (or other organic molecules) is oxidized through a series of biochemical reactions, ultimately leading to the generation of ATP. The energy released during each step of these reactions is used to drive the synthesis of ATP from ADP (Adenosine Diphosphate) and inorganic phosphate (Pi).

    The energy released from the de-excitation of a carbon atom in an SP3 hybridized state, which you mentioned as 2 eV, contributes to the overall energy available for cellular respiration. However, it's important to understand that this energy is not directly used to "distribute" to ATP molecules. Instead, it's used indirectly through the electron transport chain and oxidative phosphorylation.

    In brief, during cellular respiration, the energy from the oxidation of glucose is used to create a proton gradient across the inner mitochondrial membrane. This proton gradient is then used to drive the synthesis of ATP from ADP and Pi by the enzyme ATP synthase. This process is called oxidative phosphorylation.

    So, while the energy from the de-excitation of the carbon atom contributes to the overall energy available for cellular respiration, it doesn't directly distribute to ATP molecules. Instead, it's used indirectly to generate the proton gradient that ultimately drives ATP synthesis.

    [Top]


    Origin of Life

    In the perspective of looking forward, Figure 21 portrays the supposed sequence of events from the formation of the Earth to LUCA.

                                        

    Figure 21 Transition to Life [view large image]

    LUCA Figure 22 shows a version of LUCA completed with all the essential parts and the life sustaining process from DNA transcription to protein production. It also indicates the flow of information and energy from/to the environment. This is just the process of entropy dissipation as mention earlier. The original article "Key Steps in the Early Evolution of Life from the Origin of Protein Synthesis to Modern Cellular Life" is a long long discourse on the origin of life starting from the RNA world to LUCA. It contains a lot more detail than the sketch presented here. In this article, it is the chemosynthesis which uses inorganic molecules (such as iron, hydrogen, sulfur, and methane)

    Figure 22 LUCA
    [view large image]


    as a source of energy and convert them into organic substances. The simple organism consumed the chemicals from the environment and employed enzymes to speed up chemical reactions, which released energy into phosphate bonds carried by ATP or other similar compounds. This scenario is more suitable for Hydrothermal Mounds origin of life as opposite to "The Hot Spring Hypothesis". (see "Bacterial Metabolism" for various ways).

    Go to Top of Page to Select
     or to Main Menu

    .