Standard Thermodynamic Properties, Biosynthesis Rates, and the Driving Force of Growth of Five Agricultural Plants

Popovic, Marko; Minceva, Mirjana

doi:10.3389/fpls.2021.671868

ORIGINAL RESEARCH article

Front. Plant Sci., 31 May 2021

Sec. Plant Biophysics and Modeling

Volume 12 - 2021 | https://doi.org/10.3389/fpls.2021.671868

This article is part of the Research Topic Plant Thermodynamics View all 5 articles

Standard Thermodynamic Properties, Biosynthesis Rates, and the Driving Force of Growth of Five Agricultural Plants

$\r\nMarko Popovic*$ Marko Popovic^* $Mirjana Minceva\r\n$ Mirjana Minceva

Biothermodynamics, TUM School of Life Sciences, Technical University of Munich, Freising, Germany

Elemental composition of Gossypium hirsutum L. (cotton), Oryza sativa L. (Asian rice), Phaseolus vulgaris L. (common bean), Saccharum spp. L. (sugarcane), and Zea mays L. (corn) was used to calculate their empirical formulas (unit carbon formulas) and growth stoichiometry. The empirical formulas were used to find standard enthalpy of formation, standard molar entropy, standard Gibbs energy of formation, and standard molar heat capacity. A comparison was made between thermodynamic properties of live matter of the analyzed plants and other unicellular and multicellular organisms. Moreover, the growth process was analyzed through standard enthalpy, entropy, and Gibbs energy of biosynthesis. The average standard Gibbs energy of biosynthesis was found to be +463.0 kJ/C-mol. Thus, photosynthesis provides energy and carbon for plant growth. The average intercepted photosynthetic energy was found to be 15.5 MJ/C-mol for the analyzed plants. However, due to inefficiency, a great fraction of the intercepted photosynthetic energy cannot be used by plants. The average usable photosynthetic energy was found to be –2.3 MJ/C-mol. The average thermodynamic driving force for growth is –1.9 MJ/C-mol. Driving forces of growth of C3 and C4 plants were compared. It was found that C4 plants have a greater driving force of growth than C3 plants, which reflects the greater efficiency of C4 photosynthesis. The relationship between the driving force and growth rates was analyzed by determining phenomenological L coefficients. The determined phenomenological coefficients span two orders of magnitude, depending on plant species and environmental conditions. The L coefficient of P. vulgaris was found to be lower than that of other plants, due to additional energy requirements of nitrogen fixation.

GRAPHICAL ABSTRACT

Graphical Abstract. Intercepted photosynthetic energy q is partly converted into usable photosynthetic energy, Δ_psG, with an efficiency μ. A part of Δ_psG is used for biosynthesis Δ_bsG, which is represented by the growth reaction in the bottom, where “Bio” denotes newly formed live matter. The remaining energy is the driving force of growth, Δ_DFG, used to make growth at a desired rate. The background for the figure was taken from Pixabay (https://www.pixabay.com/), image ID number: 4779065, under Pixabay license.

Introduction

Thermodynamics has a great potential in life sciences (Von Stockar and Liu, 1999; Von Stockar, 2010, 2013a; Lucia, 2015). Thermodynamics has been used in systems biology to analyze metabolic pathways and analyze metabolic strategies evolved by microorganisms (Von Stockar and Liu, 1999; Von Stockar, 2013b). Thermodynamics is used more and more often in biology for analysis of microorganism cells and communities (Hellingwerf et al., 1982; McInerney and Beaty, 1988; Del Giorgio and Cole, 1998; Kleerebezem and Van Loosdrecht, 2010; Soh and Hatzimanikatis, 2010), as well as plants (Dragicevic and Sredojevic, 2011; Keller, 2013). Morowitz (1968) was the first to roughly estimate thermodynamic properties of ubiquitously present microorganisms using statistical thermodynamics. Afterward, the interest in biothermodynamics increased, which motivated researchers to study the subject in more detail. Wang et al. (1976) and Duboc et al. (1999) reported experimental enthalpies of combustion of microorganisms. The next great advance in the field was made by Battley et al. (1997) by experimentally determining the entropy of Saccharomyces cerevisiae. Based on this data, Battley (1998) has made a very detailed analysis of S. cerevisiae growth, determining the change in specific enthalpy, entropy, and Gibbs free energy for the process in both aerobic and anaerobic modes and on several substrates. Barros and Feijóo (2003) made a quantitative analysis of soil microorganism growth, describing it with a growth reaction. Multicellular and unicellular organisms represent a highly organized amount of substance, containing nucleic acid, proteins, lipids, and carbohydrates, characterized by specific elemental composition (empirical formula) and thermodynamic properties (enthalpy, entropy, and Gibbs energy) (Von Stockar, 2013a; Ozilgen and Sorguven, 2016). The conversion of energy by plants must be performed in full accordance with the general laws of thermodynamics (Ksenzhek and Volkov, 1998). Plant growth rate is proportional to the product of the metabolic rate and the metabolic efficiency for production of anabolic products (Criddle et al., 2005).

Thermodynamics has been used to study plant growth and photosynthesis. Hansen et al. (1994) related metabolic heat and CO₂ rates with the plant growth rate and carbon-conversion efficiency, by developing a model based on physiological parameters. Criddle et al. (2005) found that plants have evolved to have maximum growth rates at temperatures most common in their environment. Moreover, analysis of the energetics of respiration shows that measurements of metabolic heat and CO₂ rates by calorespirometry, combined with estimates of substrate and biomass composition, are sufficient to calculate substrate carbon-conversion efficiencies, anabolic rates or rates of growth and development, and relative activities of metabolic paths (Macfarlane et al., 2002; Hansen et al., 2005). Dragicevic and Sredojevic (2011) used thermodynamic principles to analyze plant physiological processes, including seed germination and plant growth. Keller (2013) analyzed the energetics of photosynthesis. Enquist et al. (1998, 1999) found that energetic constraints on metabolism are reflected in population density, and other ecological and evolutionary phenomena, including resource allocation among species in ecosystems. Attorre et al. (2019) developed a model of plant growth based on non-equilibrium thermodynamics, which was able to predict CO₂ fixation by plants, as well as their response to changes in their environment. Colombi et al. (2019) used isothermal calorimetry to study energy used for root growth and penetration through the soil. Herrmann and Colombi (2019) used calorespirometry to study energy use efficiency of root growth. Ksenzhek and Volkov (1998) analyzed the applications of thermodynamic principles to plant growth and metabolism.

Growth of organisms requires a thermodynamic driving force, in the form of negative Gibbs energy (Von Stockar, 2013b). While the value of the driving force of higher plants remains unknown, it has been determined for heterotrophic microorganisms (Von Stockar and Liu, 1999; Von Stockar, 2013b) and microalgae. Von Stockar et al. (2011) used calorimetry to measure thermodynamic properties of photosynthetic growth of the microalgae Chlorella vulgaris, strain CCAP 211/11B. They found that the culture absorbed –5,000 kJ of light per C-mol of biomass formed (Von Stockar et al., 2011). This energy is partly used to cover the cost of biosynthesis of new biomass, which was measured to be +458 kJ/C-mol (Von Stockar et al., 2011). The remaining –4,500 kJ/C-mol represents the driving force for growth (Von Stockar et al., 2011). This excess driving force is required to keep the reaction sufficiently far from equilibrium for it to occur at a desired rate, since the reaction rate is proportional to distance from equilibrium (Hellingwerf et al., 1982; Westerhoff et al., 1982; Balmer, 2011; Von Stockar, 2013b,c; Demirel, 2014; Popovic and Minceva, 2020a,b).

Solar radiation is the ultimate source of energy for most life, with the exceptions of hydrothermal vents in the deep ocean. It is harvested by autotrophs and then distributed through food webs in the form of chemical energy. The role of photosynthesis in plant growth has been extensively analyzed in the literature (Pearcy et al., 1981; Allen et al., 1988; Kirschbaum, 2011; Holding and Streich, 2013; Weraduwage et al., 2016). However, the absorbed light energy goes through intermediate steps before it can be used to drive metabolic processes. It is first converted into chemical energy in the form of photosynthate, which is then partly oxidized in catabolism to produce ATP for metabolic processes (e.g., maintaining homeostasis and growth). A detailed discussion of thermodynamics of plant catabolism and respiration is given in Hansen et al. (1994, 2004, 2021), Gary et al. (1995), and Ellingson et al. (2003).

The goal of this research is to make a thermodynamic analysis of growth of five widely cultivated plants: Gossypium hirsutum L. (cotton), Oryza sativa L. (Asian rice), Phaseolus vulgaris L. (common bean), Saccharum spp. L. (sugarcane), and Zea mays L. (corn). Elemental compositions of the analyzed plants will be used to calculate empirical formulas, thermodynamic properties of plant tissues, and their growth, as well as to formulate growth reactions. After that, the calculated thermodynamic properties will be compared with those of viruses, bacteria, fungi, archaea, and human tissues. This will show whether there are significant differences in the driving force of growth of phototrophic and chemotrophic organisms, as well as between prokaryotes and eukaryotes. Moreover, the driving force will be related to the growth rate using nonequilibrium thermodynamics.

Materials and Methods

Like all other organisms, plants have characteristic empirical formulas and thermodynamic properties of live matter. Their determination, starting from elemental composition data, is discussed in the Properties of Live Matter section. The empirical formulas are then used to find stoichiometry of growth and thermodynamic properties of biosynthesis, as discussed in the Properties of Growth and Biosynthesis section. The determined properties have an uncertainty, as discussed in the Uncertainties section.

Properties of Live Matter

The starting point for this research is elemental composition of plant dry matter reported in the literature. Elemental composition of live matter is measured as mass fractions, which are then converted into unit carbon formulas, also known as C-mol or empirical formulas (Battley, 2013). The general unit carbon formula of the analyzed plants is C_nCH_nHO_nON_nNP_nPS_nSNa_nNaK_nKMg_nMgCa_nCaCl_nCl, where n_J is the atomic coefficient of element J, which can be calculated using the equation

n_{J} = \frac{w_{J}}{w_{C}} \frac{M_{C}}{M_{J}} (1)

where w_J and w_C are the mass fractions of element J and carbon in the biomass, respectively, while M_J and M_C are the molar masses of element J and C, respectively (Duboc et al., 1999). Empirical formulas are reported for organism dry matter (excluding water) (Duboc et al., 1999). Molar masses of empirical formulas, M_bio, can be found through the equation (Duboc et al., 1999)

M_{b i o} = \frac{1}{w_{C}} M_{C} (2)

Thermodynamic properties of live matter were calculated as described in, Patel and Erickson (1981); Hurst and Harrison (1992), Battley (1998, 1999), and Popovic (2019). Elemental composition of animate matter can be used to determine its enthalpy of formation through the Patel–Erickson equation and classical reaction thermochemistry. The Patel–Erickson equation relates heat released during combustion to number of electrons transferred to oxygen (Patel and Erickson, 1981; Battley, 1998):

△_{C} H^{0} (b i o) = - 111.14 \frac{kJ}{mol} \cdot E (3)

where Δ_CH⁰ is standard enthalpy of combustion of live matter and E is the number of electrons transferred to oxygen during combustion to CO₂(g), H₂O(l), N₂(g), P₄O₁₀(s), and SO₃(g) (for a discussion on other conventions concerning SO₃ please see Popovic, 2019). During combustion, a C atom gives its four valence electrons to O, H gives one, N gives none since it is converted to N₂, P gives five, and S gives six. Inorganic ions, like Na⁺ and Mg²⁺, are not included, since they are already in their highest oxidation state and cannot transfer any electrons to oxygen (Battley, 1998). Thus, E is calculated through the equation:

E = 4 n_{C} + n_{H} - 2 n_{O} - 0 n_{N} + 5 n_{P} + 6 n_{S} (4)

where n_C, n_H, n_O, n_N, n_P, and n_S are, respectively, the number of C, H, O, N, P, and S atoms in the biomass empirical formula (Patel and Erickson, 1981; Battley, 1998). If any of these atoms are not present in the empirical formula, they are just neglected during the calculation (Battley, 1998). Standard enthalpy of combustion of live matter Δ_CH⁰(bio) is the standard enthalpy change of the reaction.

\begin{array}{l} C_{n C} H_{n H} O_{n O} N_{n N} P_{n P} S_{n S} K_{n K} M g_{n M g} C a_{n C a} A l_{n A l} S i_{n S i} M n_{n M n} F e_{n F e} \\ C l_{n C l} + (n_{C} + ¼ n_{H} + 1¼ n_{P} + 1½ n_{S} + ¼ n_{K} + ½ n_{M g} \\ + ½ n_{C a} + ¾ n_{A l} + n_{S i} + n_{M n} + ¾ n_{F e} - ½ n_{O} - ¼ n_{C l}) O_{2} \to \\ n_{C} C O_{2} + ½ n_{H} H_{2} O + ½ n_{N} N_{2} + ¼ n_{P} P_{4} O_{10} + n_{S} S O_{3} + ½ n_{K} \\ K_{2} O + n_{M g} M g O + n_{C a} C a O + ½ n_{A l} A l_{2} O_{3} + n_{S i} S i O_{2} + n_{M n} \\ M n O_{2} + ½ n_{F e} F e_{2} O_{3} + n_{C l} H C l (5) \end{array}

Thus, Δ_CH⁰ can be used to find standard enthalpy of formation, Δ_fH⁰(bio), of live matter, through the Hess law (Battley, 1998; Popovic, 2019).

\begin{array}{l} △_{f} H^{0} (b i o) = n_{C} △_{f} H^{0} (C O_{2}) + \frac{1}{2} n_{H} △_{f} H^{0} (H_{2} O) + \frac{1}{4} n_{P} △_{f} H^{0} \\ (P_{4} O_{10}) + n_{S} △_{f} H^{0} (S O_{3}) + \frac{1}{2} n_{K} △_{f} H^{0} (K_{2} O) + n_{M g} △_{f} H^{0} (M g O) \\ + n_{C a} △_{f} H^{0} (C a O) + \frac{1}{2} n_{A l} △_{f} H^{0} (A l_{2} O_{3}) + n_{S i} △_{f} H^{0} (S i O_{2}) \\ + n_{M n} △_{f} H^{0} (M n O_{2}) + \frac{1}{2} n_{F e} △_{f} H^{0} (F e_{2} O_{3}) + n_{C l} △_{f} H^{0} (H C l) \\ - △_{C} H^{0} (b i o) (6) \end{array}

Elemental composition can also be used to determine standard molar entropy of live matter, S⁰_m(bio), through the Battley equation (Battley, 1999)

S_{m}^{o} (b i o) = 0.187 \sum_{J} \frac{S_{m}^{o} (J)}{a_{J}} n_{J} (7)

where n_J is the number of atoms of element J in the empirical formula of the biomass, S⁰_m(J) is standard molar entropy of element J, and a_J is the number of atoms per molecule of element J in its standard state elemental form. For example, the standard state elemental form of carbon is graphite, which is simply written as C, which makes a_C = 1. On the other hand, hydrogen, oxygen, and nitrogen are in their standard state elemental forms of all diatomic gasses H₂, O₂, and N₂, respectively, which implies that a_H = a_O = a_N = 2. The summation is over all elements constituting the dry live matter. The Battley equation can also be used to determine standard entropy of formation of live matter Δ_fS⁰(bio). In this case, it takes the form (Battley, 1999)

△_{f} S^{0} (b i o) = - 0.813 \sum_{J} \frac{S_{m}^{o} (J)}{a_{J}} n_{J} (8)

Finally, by combining standard enthalpy and entropy of formation, it is possible to calculate standard Gibbs energy of formation of live matter Δ_fG⁰(bio), as

△_{f} G^{0} (b i o) = △_{f} H^{0} (b i o) - T △_{f} S^{0} (b i o) (9)

where T is the temperature.

Standard molar heat capacities of the analyzed plants were determined using the Hurst–Harrison equation. The Hurst–Harrison equation is a group contribution model, describing standard molar heat capacity of a substance, C⁰_p,_m, as a sum of contributions from its constituent elements

C_{p, m}^{o} = \sum_{J} n_{J} c_{J} (10)

where n_J is the number of atoms of element J in the substance formula, while c_J is an empirical parameter that describes the contribution of element J to heat capacity (Hurst and Harrison, 1992; Ozilgen and Sorguven, 2016). The values of c_J can be found in Hurst and Harrison (1992).

Properties of Growth and Biosynthesis

Growth is a chemical process in which an organism converts nutrients into new live matter, a process that can be represented by chemical reactions, known as growth reactions (Raven, 1985; Battley, 1998, 2013; Von Stockar, 2013b; Pineda-Insuasti et al., 2014; Popovic and Minceva, 2020a). Growth reactions describe the stoichiometry with which nutrients are converted into live matter.

The non-metals C, H, N, O, P, and S, also known as CHNOPS elements, are essential for all living systems, forming the majority of cell structures (Wackett et al., 2004). Carbon enters photosynthetic tissues of plants as carbon dioxide. Hydrogen and oxygen enter plants mostly through the root as water. Nitrogen enters plants mainly through roots in the form of ammonium or nitrate salts. Phosphorus enters plants predominantly from soil in the form of phosphate salts. Sulfur is taken up by plants mostly from soil in the sulfate form (Linzon et al., 1979). Except through soil, plants can take many elements through leaves from foliar fertilizers, which can reduce the total amount of fertilizer applied and achieve high fertilizer efficiency (Niu et al., 2021).

Except for the CHNOPS elements, plants use a great number of metallic ions, which are taken predominantly from the soil through the roots. Potassium greatly influences plant growth and is taken as K⁺ (Ragel et al., 2019; Xu et al., 2020). Magnesium is important in protein synthesis and is imported by plants in the form of Mg²⁺ (Verbruggen and Hermans, 2013; Guo et al., 2016; Hauer-Jákli and Tränkner, 2019). Calcium is taken as Ca²⁺ (White et al., 2002; White and Broadley, 2003) and has a dual function as a structural component of cell walls (Hepler, 2005; Thor, 2019) and membranes and as intracellular second messenger (Knight et al., 1997; Blume et al., 2000; Michard et al., 2011; Monshausen et al., 2011; Wilkins et al., 2016; Ortiz-Ramírez et al., 2017; Tang and Luan, 2017; Kudla et al., 2018; Thor, 2019). Silicon is taken in the form of silicic acid Si(OH)₄ (Ma and Yamaji, 2006) and has beneficial effects on growth, development, yield, and disease resistance of many plant species (Ma, 2004; Ma and Yamaji, 2006; Luyckx et al., 2017). Aluminum is present in plants as Al³⁺ ion (Bojórquez-Quintal et al., 2017). In small amounts, aluminum stimulates root growth, increases nutrient uptake, and improves enzyme activity but inhibits root growth at higher concentrations (Bojórquez-Quintal et al., 2017). Manganese is taken as Mn²⁺ and is an important cofactor for catalyzing the water-splitting reaction in photosystem II (Alejandro et al., 2020). Iron is present in the soil mostly as Fe³⁺ ions but can also be present as Fe²⁺ in smaller amounts. Iron is imported by plants using various strategies (Morrissey and Guerinot, 2009; Jacoby et al., 2017; Mhlongo et al., 2018) and is used in many processes, including electron transfer chains and photosynthesis (Schmidt et al., 2020). Chlorine is taken as Cl^– and is important in osmoregulation, stomata opening, and photosynthesis (Critchley, 1985; Coleman et al., 1987; Chen et al., 2010; Marschner, 2011; Colmenero-Flores et al., 2019).

Other elements, such as Cu and Zn, are also used by plants. However, the literature did not contain sufficient data to include these elements into the calculations. Since their amount is relatively small compared with the CHNOPS elements (Wackett et al., 2004), they do not influence thermodynamic properties of plants. However, it would be interesting to consider their influence on plant nutrition in the future from the perspective of biothermodynamics and growth reactions.

Based on the discussion above, the general growth reaction proposed for the analyzed plants is

\begin{array}{l} C O_{2} (g) + H_{2} O (l) + N H_{4}^{+} (a q) + H_{2} P O_{4}^{-} (a q) + S O_{4}^{2 -} \\ (a q) + K^{+} (a q) + M g^{2 +} (a q) + C a^{2 +} (a q) + A l^{3 +} (a q) \\ + S i (O H)_{4} (S) + M n^{2 +} (a q) + F e^{3 +} (a q) + C l^{-} (a q) \to \\ B i o + O_{2} (g) + H^{+} (a q) (11) \end{array}

where “Bio” denotes live matter (excluding water). Stoichiometric coefficients of nutrients and metabolic products are specific for each plant species and are discussed in the Live Matter Elemental Composition and Stoichiometry of Growth section.

Growth reaction thermodynamic parameters were calculated using classical thermochemistry

△_{b s} H^{0} = \sum_{p r o d u c t s} ν △_{f} H^{0} - \sum_{r e a c t a n t s} ν △_{f} H^{0} (12)

△_{b s} S^{0} = \sum_{p r o d u c t s} ν S_{m}^{o} - \sum_{r e a c t a n t s} ν S_{m}^{o} (13)

△_{b s} G^{0} = \sum_{p r o d u c t s} ν △_{f} G^{0} - \sum_{r e a c t a n t s} ν △_{f} G^{0} (14)

where ν’s are stoichiometric coefficients of species participating in the reaction, while Δ_bsH⁰, Δ_bsS⁰, and Δ_bsG⁰ are standard enthalpy, entropy, and Gibbs energy of biosynthesis, respectively (Atkins and de Paula, 2011; Atkins et al., 2017).

Uncertainties

Thermodynamic properties of the studied plants were determined from elemental composition using empirical relations and thus have some uncertainty. Δ_CH⁰ was found using the Patel–Erickson equation, the uncertainty of which is 5.36% (Popovic, 2019). The determined Δ_CH⁰ values were then subtracted from standard enthalpies of formation of oxides (Equation 6) to find Δ_fH⁰(bio). Since standard enthalpies of formation of oxides were precisely determined by experiment (more details in Chase, 1998), they have a negligible error compared with that in Δ_CH⁰. Thus, the uncertainty in standard enthalpy of formation of live matter, δ(Δ_fH⁰(bio)), is equal to the error in Δ_CH⁰.

δ (△_{f} H^{0} (b i o)) = 0.0536 \cdot [111.14 \frac{k J}{m o l} (4 n_{C} + n_{H} - 2 n_{O} - 0 n_{N} + 5 n_{P} + 6 n_{S})] (15)

The uncertainty in estimation of entropy using the Battley equation is 2% for dry matter and 19.7% for hydrated matter (Battley, 1999). The greater of the two values was taken, and hence the uncertainty in standard molar entropy of live matter, δ(S⁰_m(bio)), is

δ (S_{m}^{0} (b i o)) = 0.197 \cdot S_{m}^{0} (b i o) (16)

Δ_fS⁰(bio) is defined as the standard entropy change of the reaction of formation of live matter from its constituent elements in their standard forms.

\begin{array}{l} n_{c} C (s, g r a p h i t e) + ½ n_{H} H_{2} (g) + ½ n_{O} O_{2} (g) + ½ n_{N} N_{2} (g) \\ + n_{p} P (s, w h i t e) + n_{s} S (s, r h o m b i c) + n_{N a} N a (s) + n_{K} K (s) \\ + n_{M g} M g (s) + n_{C a} C a (s) + ½ n_{C l} C l_{2} (g) \to C_{n C} H_{n H} O_{n O} \\ N_{n N} P_{n P} S_{n S} N a_{n N a} K_{n K} M g_{n M g} C a_{n C a} C l_{n C l} (17) \end{array}

and is defined as the difference in S⁰_m(bio) and standard molar entropies of the elements, which have been determined with great accuracy by experiment (Chase, 1998). Thus, the uncertainty in Δ_fS⁰(bio) is equal to that in S⁰_m(bio) (Popovic, 2019).

Δ_fH⁰(bio) and Δ_fS⁰(bio) are used to find Δ_fG⁰(bio). Therefore, the uncertainty in the standard Gibbs energy of formation of live matter, δ(Δ_fG⁰(bio)), is

δ (△_{f} G^{0} (b i o)) = \sqrt{{[δ (△_{f} H^{0} (b i o))]}^{2} + {[T \cdot δ (S_{m}^{0} (b i o))]}^{2}} (18)

Finally, the uncertainty in Δ_fG⁰(bio) is equal to that in Δ_bsG⁰, since it is the greatest source of uncertainty in its determination. Δ_bsG⁰ is determined using Equation (14), as the difference of Δ_fG⁰ values of reactants and products of reaction (Equation 11). The Δ_fG⁰ values of all reaction participants, except for live matter, have been determined with great accuracy by experiment (Chase, 1998). Thus, uncertainty in Gibbs energy of biosynthesis, δ(Δ_bsG⁰), is equal to δ(Δ_fG⁰(bio)). Similarly, δ(Δ_bsH⁰) and δ(Δ_bsS⁰) are equal to δ(Δ_fH⁰(bio)) and S⁰_m(bio), respectively.

Results and Discussion

Elemental composition and thermodynamic properties of plants and their growth are discussed in the following sections.

Live Matter Elemental Composition and Stoichiometry of Growth

Elemental compositions of the five plant species were taken from the literature (Latshaw and Miller, 1924; Medic et al., 2012; Titiloye et al., 2013; Pineda-Insuasti et al., 2014; Nam and Capareda, 2015; Fialho et al., 2019). Information about the plant samples is summarized in Table 1, including geographic origin, plant parts, and whether experimental enthalpy of combustion was reported. Calorimetrically determined enthalpies of combustion were reported for Gossypium hirsutum L., Oryza sativa L., Saccharum spp. L., and Zea mays L. stalk and cobs from Ghana (Titiloye et al., 2013; Nam and Capareda, 2015). Since for these samples Δ_CH⁰(bio) was determined directly by experiment instead of Equation (4), their standard enthalpies of formation, Δ_fH⁰(bio), are more accurate. Elemental compositions of all the analyzed plants have been reported for dry matter.

TABLE 1

Table 1. Information about the plant samples considered in this research.

Based on the elemental composition of plants, unit carbon formulas were calculated, as described in the Properties of Live Matter section. The results are shown in Table 2. It is interesting to compare the empirical formulas in Table 2 with average formulas of some of the major classes of organisms: bacteria (CH_1.₇O_0.₄N_0.2), fungi (CH_1.₇O_0.₅N_0.1), and algae (CH_1.₇O_0.₅N_0.1) (Popovic, 2019). The C:H ratio of the analyzed plants varies from 1.54 to 1.85, with an average of 1.63, which is close to that of the other classes of organisms. However, the analyzed plants have a C:O ratio from 0.71 to 0.90, with an average of 0.80, which is a significantly higher than the C:O ratio of the other classes of organisms. On the other hand, the C:N ratio of the analyzed plants varies from 7 × 10^–4 to 0.15, with an average of 0.02, which is lower than that of the other classes of organisms. The great span of C:N ratios in plants originates from its dependence on growth conditions, including soil, meteorological factors, agricultural practices, and application of nitrogen fertilizers, as well as availability of other nutrients involved in N accumulation and cycling in crops. Thus, the average oxidation state of carbon in the analyzed plants is 0.00 (from –0.24 to 0.33), which is slightly higher than the averages for bacteria (–0.3), fungi (–0.4), and algae (–0.4). Moreover, the oxidation state of carbon in carbohydrates is also 0.00, which is in accordance with the fact that carbohydrates are a major constituent of plants.

TABLE 2

Table 2. Elemental composition of the plants considered in this work.

Based on the elemental composition, growth reactions were made and stoichiometric coefficients calculated, as described in the Properties of Growth and Biosynthesis section, which are summarized in Table 3 and Supplementary Material 3. Growth of plants is autotrophic, using CO₂ as a source of carbon, which is reduced using energy from light. Moreover, data in Supplementary Material 3 show that the greatest amount of O₂ is produced by growth of Z. mays grain (1.06 O₂-mol/C-mol), while the lowest is that of Phaseolus vulgaris waste (0.92 O₂-mol/C-mol), with the average being 1.00 O₂-mol/C-mol. This is in good agreement with the average oxidation state of carbon of 0.00, since while going from +4 in CO₂ to 0.00 in biomass, 1 mol of carbon takes four electrons, which is exactly the amount released when 1 mol of O₂ is produced from H₂O (half reaction: 2H₂O → O₂ + 4H⁺ + 4e^–).

TABLE 3

Table 3. Formulas giving stoichiometric coefficients for the plant growth reactions.

Thermodynamic Properties of Live Matter and Biosynthesis

Based on elemental composition, standard thermodynamic properties of the analyzed plants were determined, as described in the Properties of Live Matter section. Table 4 gives standard thermodynamic properties of the analyzed plants and plant organs, including standard enthalpy of formation, standard molar entropy, standard Gibbs energy of formation, and standard molar heat capacity at constant pressure. All the standard enthalpies of formation are negative, implying that the biomass has a lower energy content than its constituent elements, due to partial oxidation of all other elements by oxygen and nitrogen (Erickson et al., 1978; Popovic, 2019). All the standard molar entropies are positive, in accordance with the third law of thermodynamics (Atkins and de Paula, 2011; Atkins et al., 2017). Standard molar entropies, S⁰_m, per C-mol of live matter are around 37 J/C-mol K, laying between that of graphite, 5.740 J/mol K, and carbon in gaseous state, 158.10 J/mol K (Atkins et al., 2017). This indicates that the mobility of C atoms in live matter is greater than in graphite but lower than in the gaseous state. All the standard Gibbs energies of formation are negative, with the average value being –126 kJ/C-mol. The standard molar heat capacities at constant pressure are around 35 J/C-mol K.

TABLE 4

Table 4. Thermodynamic properties of live matter of the analyzed plants.

Thermodynamic properties of the analyzed plants are compared with those of other organisms in Figure 1. The data for other classes of organisms are summarized in Supplementary Material 2. Data in Figure 1 represent the average of values for 14 bacteria, 11 fungi, 11 algae, and 18 virus species, as well as 30 human tissues. Figure 1 shows that standard enthalpy and Gibbs energy of formation of the analyzed plants is among the most negative measured, but within the span of algae, which are also photosynthetic organisms. Plants have a more negative standard enthalpy and Gibbs energy of formation, indicating a greater energy content and lower oxidation state of carbon in the biomass, compared with other classes of organisms. The reason for this might be that plants gain energy directly from the sun, while heterotrophic organisms have to obtain energy from food. This is supported by the fact that algae, which are also photosynthetic organisms, also have highly negative enthalpy and Gibbs energy of formation. Standard molar entropies of all the organisms in Figure 1 are very similar. Furthermore, plants use a part of the stored energy in catabolic processes to maintain homeostasis, as well as for growth.

FIGURE 1

Figure 1. Comparison of thermodynamic properties of formation of live matter: (A) standard enthalpy of formation, Δ_fH⁰; (B) standard molar entropy, S_m⁰; (C) standard Gibbs energy of formation Δ_fG⁰. The columns represent the average thermodynamic parameters for the organism groups, while the error bars show the spread in thermodynamic properties within the groups. Data for bacteria, fungi, and algae taken from Popovic (2019). Data for human tissues taken from Popovic and Minceva (2020c). Data for viruses are taken from (Popovic and Minceva, 2020a,b).

Thermodynamic properties of biosynthesis were calculated for the analyzed plants and plant organs, as described in the Properties of Growth and Biosynthesis section. The results are given in Table 5. The average standard enthalpy and Gibbs energy of growth are +454.4 and +463.0 kJ/C-mol, respectively, while the average standard entropy of growth is –29.1 kJ/C-mol. To compare, Gibbs energies of biosynthesis of several organisms have been determined for aerobic growth on a mixture of amino acids and glucose: Escherichia coli –45.25 kJ/C-mol, Pseudomonas (strain C12B) –18.67 kJ/C-mol, Bacillus subtilis –31.75 kJ/C-mol, Cyanobacteria synechocystis (strain PCC 6803) –13.74 kJ/C-mol, and Saccharomyces cerevisiae –15.90 kJ/C-mol (Popovic and Minceva, 2020a). For S. cerevisiae, Gibbs energy of biosynthesis was determined for several environments: anaerobic on glucose –11.05 kJ/C-mol, aerobic on glucose –35.76 kJ/C-mol, aerobic on ethanol +28.86 kJ/C-mol, and aerobic on acetic acid +24.10 kJ/C-mol (Battley, 1998). Chlorella vulgaris was found to have a Gibbs energy of biosynthesis of +458 kJ/C-mol, during phototrophic growth (Von Stockar et al., 2011). Therefore, Gibbs energy of biosynthesis of heterotrophic organisms is slightly negative or sometimes slightly positive. The reason for this is that heterotrophic organisms use substrates of biological origin, with a Gibbs energy similar to their own live matter. On the other hand, the analyzed plants and C. vulgaris are phototrophic organisms and have highly positive Gibbs energies of biosynthesis of around +500 kJ/C-mol. This is due to the high energy requirement for reduction of carbon from CO₂. However, phototrophic growth is still thermodynamically feasible, due to photosynthesis (Von Stockar, 2013b).

TABLE 5

Table 5. Thermodynamic properties of biosynthesis of the analyzed plants: standard enthalpy (Δ_bsH⁰), entropy (Δ_bsS⁰), and Gibbs energy (Δ_bsG⁰) of biosynthesis.

Standard Gibbs energy of formation is by definition the Gibbs energy change of reaction (17), the hypothetic formation of live matter from elements in their standard states (Atkins and de Paula, 2011). The elements represent merely a reference state, relative to which Gibbs energy is measured (Atkins and de Paula, 2011). This reference state was chosen because it is not possible to determine absolute Gibbs energy. On the other hand, Gibbs energy of biosynthesis is the Gibbs energy change of reaction (11), representing formation of live matter from nutrients in its environment (Von Stockar, 2013b). Depending on the nutrients and environment, Gibbs energy of biosynthesis can vary even for a single organism (Battley, 1998).

Photosynthetic Energy and the Driving Force of Growth

Energetics of photosynthesis is summarized in Figure 2. The relationship between the biomass accumulation and light energy intercepted by plants is known in agriculture as radiation use efficiency, RUE, an important quantifier of crop production, combining photosynthesis and the efficiency of a plant in producing live matter (Monteith and Moss, 1977; Hatfield and Dold, 2019). RUE is defined as the slope of the line expressing the relationship between plant biomass, m_bio, and intercepted solar radiation

m_{b i o} = R U E \cdot \sum_{e m e r g e n c e}^{h a r v e s t} f_{i} R_{s p} (19)

FIGURE 2

Figure 2. Distribution of energy in photosynthesis. (A) Total energy of the sun is partly not in the photosynthetic 400–700 nm region and is partly reflected, absorbed, or transmitted by leaves. The remaining energy is used in photosynthesis and is called intercepted photosynthetic energy (q, shown in orange). (B) Intercepted photosynthetic energy is not the actual energy used by the plant. A part is lost due to quantum efficiency requirements, in conversion to glucose, and in dark and photorespiration. The remaining energy is usable photosynthetic energy (Δ_psG, shown in orange). (C) Usable photosynthetic energy is used partly to provide energy for biosynthesis (Δ_bsG) and partly dissipated to make growth occur at a desired rate (Δ_DFG).

where f_i is the fraction of intercepted photosynthetically active radiation at stage i of plant’s development and R_sp is the incoming photosynthetically active radiation (Monteith and Moss, 1977; Sadras et al., 2016; Hatfield and Dold, 2019). Photosynthetically active radiation denotes light in the 400–700 nm range (Murchie et al., 2019). The value of f_i can be found through the equation f_i = K_C,_i – 0.3, where K_C,_i is the crop coefficient at stage i of plant’s development (Sadras et al., 2016). Crop coefficients are calculated, according to the Food and Agriculture Organization (FAO) method (Doorenbos and Pruitt, 1977), by dividing the crop cycle into four stages: initial, rapid growth, maximum, and declining (Allen et al., 1998; Sadras et al., 2016). For each stage, K_C,i values have been tabulated for many plants (Allen et al., 1998). Similarly, RUE has been determined for many plants and depends on environmental conditions, such as availability of nitrogen and other nutrients, temperature, soil water, atmospheric humidity, and wind (Sadras et al., 2016; Hatfield and Dold, 2019; Murchie et al., 2019). More information about RUE can be found in Murchie et al. (2019).

RUE, defined as the biomass grown per unit light energy, can be used to calculate the intercepted photosynthetic energy, q, needed by plants to incorporate 1 mol of carbon into new biomass

q = \frac{M_{b i o}}{R U E} (20)

where M_bio is the molar mass of live matter empirical formula [mass per 1 C-mol, calculated using Equation (2)]. Equation (20) is Equation (19) written for 1 C-mol of plant biomass, with q representing the sum term. The values of q for the analyzed plants were calculated using literature RUE values and are given in Table 6.

TABLE 6

Table 6. Photosynthetic energy and driving force of growth.

It would be interesting to compare the q-values from Table 6 with the result of an alternative calculation. Thus, we will estimate q in another way. Approximately 5.7 × 10²⁴ J of solar energy is irradiated to the Earth’s surface on an annual basis (Miyamoto, 1997). Of this energy, 98–99% is reflected or absorbed by objects, leaving (1.5 ± 0.5)% to be captured by photosynthesis (Thompson et al., 2018). Thus, the total energy absorbed by photosynthetic organisms is (8.6 ± 2.5) × 10²² J/year. Photosynthetic organisms utilize this solar energy to fix 2 × 10¹¹ tons of CO₂ per year (Miyamoto, 1997). Since 1 mol of CO₂ contains 1 mol of carbon and the molar mass of CO₂ is 44 g/mol, the total amount of carbon fixed by plants is found to be 4.5 × 10¹⁵ C-mol/year. Therefore, by dividing the total energy (8.6 ± 2.5) × 10²² J/year by the amount of carbon fixed 4.5 × 10¹⁵ C-mol/year, the intercepted photosynthetic energy is found to be q = (18.8 ± 6.3) MJ/C-mol. This value is in good agreement with the q-values of the analyzed plants from Table 6.

The intercepted photosynthetic energy represents all light energy received by a plant and seems to be very large. However, photosynthesis has several sources of inefficiency: (1) 47% is lost since it is not in the 400–700 nm region used for photosynthesis; (2) 30% is lost in reflection, absorption, and transmission by leaves; (3) 24% for quantum efficiency requirements; (4) 68% in conversion to glucose; and (5) 35–40% on the processes of dark and photorespiration (Sudhakar and Mamat, 2019). Inefficiency sources (1) and (2) are taken into account in Equation (19) by R_sp and f_i, respectively. The first source of inefficiency is taken into account when calculating RUE, since Equation (19) contains only R_sp which represents only photosynthetically active 400–700 nm light. Moreover, f_i in Equation (19) takes into account the correction for reflection, absorption, and transmission by leaves. Thus, these two sources of inefficiency are taken into account when calculating q using Equation (20), which is derived from Equation (19). However, the remaining three sources of inefficiency have not been taken into account. Thus, the photosynthetic efficiency, μ, in converting the intercepted photosynthetic energy q into energy usable to plants is

μ = [100 % - (24 % quantum efficiency loss)] \cdot [100 % - (68 % conversion to glucose loss)] \cdot [100 % - (40 % photorespiration loss)] = 15 % (21)

Thus, the usable photosynthetic energy, Δ_psG, is the photosynthetic energy actually available for plants to drive their photosynthetic metabolism and is lower than the intercepted photosynthetic energy, q

△_{p s} G = - μ \cdot q (22)

The negative sign was added because plants receive energy from the sun as q and then spend it on metabolic processes as Δ_psG. A more formal derivation of Equation (22) is given in Supplementary Material 1.

The sum of the usable photosynthetic energy and the Gibbs energy required for biosynthesis, Δ_bsG⁰, represents the driving force of growth, Δ_DFG (Von Stockar et al., 2011; Von Stockar, 2013b).

△_{D F} G = △_{p s} G + △_{b s} G^{0} (23)

The values of q, Δ_psG, and Δ_DFG for the analyzed plants are given in Table 6. The Δ_DFG values greatly vary, even within a species. For example, the driving force for growth of Saccharum spp. ranges between –1.0 and –3.8 MJ/C-mol. The reason for this is great variation in reported RUE, which is very different for plants grown with irrigation and in dry farming (rain-fed) (De Silva and De Costa, 2012). Moreover, RUE of Z. mays depends on its growth location (Hatfield and Dold, 2019). However, all the driving force values, despite their differences, have a general trend: they are all negative, making growth of plants a thermodynamically favorable process.

The results for vascular plants obtained in this work can be compared with those for unicellular algae C. vulgaris. C. vulgaris was found to have an intercepted photosynthetic energy of q = 5,000 kJ/C-mol and a Gibbs energy of biosynthesis of Δ_bsG = +458 kJ/C-mol (Von Stockar et al., 2011). The average values for the plants analyzed in this research are q = 15,537 kJ/C-mol and Δ_bsG = +469 kJ/C-mol. To compare the q-values, the photosynthetic efficiency, μ, has to be taken into account through Equation (22). Thus, usable photosynthetic energy, Δ_psG, is –750 kJ/C-mol for C. vulgaris and –2,330 kJ/C-mol for the plants analyzed in this research. These values are then added to the corresponding Δ_bsG (Equation 23) to obtain the driving forces of growth, Δ_DFG, which are –292 kJ/C-mol for C. vulgaris and –1,861 kJ/C-mol for the plants analyzed in this research. The difference between these two values most likely originates from additional processes present in multicellular organisms, such as nutrient uptake by the root (Sondergaard et al., 2004) and great energy used by the root to penetrate through the soil during growth (Colombi et al., 2019; Herrmann and Colombi, 2019). In nutrient limited conditions, plants will spend energy to acquire the limiting nutrient, in accordance with the law of the minimum. Moreover, photosynthesis does not directly drive growth but generates sugars that are oxidized in aerobic respiration to make ATP, which is hydrolyzed to release energy required for plant growth (Hansen et al., 1994, 2004, 2021; Gary et al., 1995; Ellingson et al., 2003). Each step in this process dissipates energy gained during photosynthesis.

The comparison can be carried further by including heterotrophic organisms. Heterotrophic organisms have a driving force of growth of approximately Δ_DFG = –500 kJ/C-mol (Heijnen and Van Dijken, 1992; Von Stockar and Liu, 1999; Liu et al., 2007; Von Stockar, 2013b). The corresponding Δ_DFG values for photosynthetic growth are –292 kJ/C-mol for C. vulgaris and –1,861 kJ/C-mol for the plants analyzed in this research. The value for C. vulgaris is on the same order of magnitude as that for heterotrophic organisms. However, the plants analyzed in this research have a significantly greater driving force, which is most likely due to greater energy demands required by multicellular organisms, as discussed above.

Data in Table 6 can be used to compare energetics of C3 and C4 plants. Saccharum spp. L. have on average the greatest usable photosynthetic energy Δ_psG and driving force Δ_DFG. The reason is that Saccharum spp. L. is a C4 plant and is thus more efficient than C3 plants. Representatives of C3 plants, G. hirsutum L., and Phaseolus vulgaris L. have lower photosynthetic energy Δ_psG and driving force Δ_DFG. The lowest values are those of O. sativa L., which is also a C3 plant. Z. mays L. has the greatest span of Δ_psG and Δ_DFG, due to many reported values. Many values have been reported because Z. mays is grown throughout the world, under a wide range of conditions. Under optimal conditions, Z. mays has a driving force of Δ_DFG = –(1.3 + 1.5) = –2.8 MJ/C-mol, which is the one of the highest values for the analyzed plant species. This is in good agreement with it being a C4 plant.

Driving Force and Growth Rate

The driving force of growth, in the form of negative Gibbs energy, exists to insure that growth occurs at a desired rate. If a process occurs near equilibrium, its Gibbs energy is very close to zero, making it thermodynamically very efficient (Atkins and de Paula, 2011; Atkins et al., 2017). However, processes near equilibrium occur at an infinitely slow rate (Atkins and de Paula, 2011; Demirel, 2014; Atkins et al., 2017), which would not be practical for growth of organisms. Thus, organisms have to make a compromise, between energetic efficiency and growth rate (Von Stockar and Liu, 1999; Von Stockar, 2013c). In other words, they have to waste a certain amount of energy to make growth occur at a desired rate (Von Stockar and Liu, 1999; Von Stockar, 2013c). This wasted energy is the driving force for growth (Von Stockar and Liu, 1999; Von Stockar, 2013c). The driving force, Δ_DFG, is hence related to the growth rate of the plant, r, through the equation

r = - \frac{L}{T} △_{D F} G (24)

where L is a phenomenological coefficient and T is temperature (Hellingwerf et al., 1982; Westerhoff et al., 1982; Von Stockar, 2013c; Demirel, 2014; Popovic and Minceva, 2020a,b). The phenomenological coefficient is a proportionality constant that takes into account the kinetic factors (Demirel, 2014) and can be calculated from Equation (24) if the driving force and growth rate are known.

Plant growth data were taken from the literature for G. hirsutum L. (Shi et al., 2013), P. vulgaris L. (Shi et al., 2013), and Z. mays L. (Shi et al., 2013; Koca and Erekul, 2016). The growth data were analyzed using the three-phase linear model (Buchanan et al., 1997), which distinguishes three phases in growth: lag, exponential, and stationary (Buchanan et al., 1997; Pérez-Rodríguez and Valero, 2013; Pla et al., 2015). These three phases are analogous to the initial, rapid growth, maximum crop cycle stages, from the FAO method for calculating crop coefficients (Doorenbos and Pruitt, 1977). The original three-phase linear model was modified to describe plant growth

m = m_{0} for t \leq t_{l a g} (25a)

m = m_{0} + r \cdot (t - t_{l a g}) for t_{l a g} < t < t_{m a x} (25b)

m = m_{m a x} for t \geq t_{m a x} (25c)

r = \frac{m_{m a x} - m_{0}}{t_{m a x} - t_{l a g}} (25d)

where m is plant mass at time t, m₀ is plant mass before the intensive growth phase, t_lag is the duration of the lag phase (time from planting to the intensive growth phase), r is the growth rate, and t_max is the moment the steady state is reached and growth ceases. The model has four adjustable parameters: m₀, m_max, t_lag, and t_max. The three-phase linear model is simple, making it a good primary model for obtaining input data for studying more complicated questions, such as the environmental influence on the growth rate (Garthright and Buchanan, 1997). In our case, the three-phase linear model was chosen since it provides a simple and impartial way of determining plant growth rates from growth data. Fitting was made using least squares regression. The fitting results are presented in Figure 3. Growth rate r was calculated as the slope of the fitted line using Equation (25d).

FIGURE 3

Figure 3. Plant growth curves. The blue circles ( www.frontiersin.org ) represent experimental data, while the orange lines () represent a fit made with the three-phase linear model (Equation 25). Experimental data sources: (A–C) from Shi et al. (2013) and (D) from Koca and Erekul (2016).

The growth rates obtained from fitted functions were combined with Δ_DFG (Table 6) to find the phenomenological coefficients using Equation (24). The results are given in Table 7. The obtained phenomenological coefficients vary greatly, depending on plant species and growth conditions. Phenomenological coefficients are specific for every plant species and environmental conditions, including soil, insolation, and water. Thus, the values of the phenomenological coefficients can vary even within one species, if it is cultivated under various conditions. For example, Z. mays cultivated in Turkey and China have different phenomenological coefficients, which are influenced by use of different seeds (hybrids) and growth on different soil, as well as temperature, water, and nutrients (Niu and Masabni, 2018).

TABLE 7

Table 7. Growth rates, r, and phenomenological L coefficients of the analyzed plants.

Data in Table 7 can be interpreted in terms of metabolic processes performed by the analyzed plants. P. vulgaris L. has the lowest phenomenological L coefficient. The reason is that it is a legume. Legumes can fix nitrogen and thus have a high nitrogen content (P. vulgaris L. Waste 1 has the greatest nitrogen content in among all the samples from Table 2). However, nitrogen fixation requires a lot of energy, making less energy available to drive growth (Sinclair and Muchow, 1999). Moreover, by comparing Z. mays L. grown in China and Turkey, it can be seen that various sorts of the same species grown under different conditions can have very different phenomenological L coefficients and growth rates.

Conclusion

Empirical formulas with both macro and micro elements were calculated for Gossypium hirsutum L. (cotton), Oryza sativa L. (Asian rice), Phaseolus vulgaris L. (common bean), Saccharum spp. L. (sugarcane), and Zea mays L. (corn). Moreover, stoichiometry of growth was determined and expressed in the form of growth reactions. Plant live matter and growth process were analyzed by determining their standard thermodynamic properties of formation and biosynthesis, respectively. The average standard Gibbs energy of formation of the analyzed plants was found to be –126 kJ/C-mol. The analyzed plants were found to have among the lowest enthalpies and Gibbs energies of formation relative to the elements, along with algae, compared with bacteria, fungi, viruses, and human tissues. The highly negative enthalpies and Gibbs energies of formation of plants and algae indicate a high energy content, originating from their ability to obtain energy directly from the sun. Standard Gibbs energies of biosynthesis are positive for all the analyzed plants, since photosynthetic growth is based on reducing carbon in CO₂. The average standard Gibbs energy of biosynthesis of the analyzed plants is +463.0 kJ/C-mol.

Photosynthesis was analyzed through three parameters: intercepted photosynthetic energy, usable photosynthetic energy, and the driving force of growth. Energetics of photosynthesis was analyzed starting from RUE, biomass grown per unit light energy, reported in the literature. The average RUE for the analyzed plants was found to be 1.8 g/MJ. RUE was used to find the intercepted photosynthetic energy, light energy in the photosynthetic 400–700 nm range received by plants. The average intercepted photosynthetic energy was found to be 15.5 MJ/C-mol for the analyzed plants. However, intercepted photosynthetic energy is not completely usable to plants. Thus, it was corrected for photosynthetic efficiency of 15%, to find the usable photosynthetic energy, photosynthetic energy actually available for plants to drive their metabolism. The average usable photosynthetic energy was found to be –2.3 MJ/C-mol for the analyzed plants. The usable photosynthetic energy was added to the standard Gibbs energy of biosynthesis, to find the driving force of growth—energy spent on keeping growth sufficiently far from equilibrium for it to occur at a desired rate. The average driving force of growth of the analyzed plants was found to be –1.9 MJ/C-mol. This value is greater than that of heterotrophic and phototrophic unicellular organisms reported in the literature, most likely due to additional energy requirements by multicellular organisms. Driving forces of growth of C3 and C4 plants were compared. It was found that C4 plants have a greater driving force of growth than C3 plants, which reflects the greater efficiency of C4 photosynthesis.

The driving force of growth is related to the growth rate, the proportionality constant being phenomenological L coefficients. Phenomenological L coefficients were determined for the first time for the analyzed plants. The phenomenological L coefficients of the analyzed plants span two orders of magnitude, depending on plant species and growth conditions. The L coefficient of P. vulgaris was found to be lower than that of other plants, due to additional energy requirements of nitrogen fixation.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

MP conceived and designed the research, collected the data, analyzed and interpreted the data, and wrote the manuscript. MM analyzed and interpreted the data, and wrote the manuscript. Both authors contributed to the article and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.671868/full#supplementary-material

References

Alejandro, S., Höller, S., Meier, B., and Peiter, E. (2020). Manganese in plants: from acquisition to subcellular allocation. Front. Plant Sci. 11:300. doi: 10.3389/fpls.2020.00300

PubMed Abstract | CrossRef Full Text | Google Scholar

Allen, R. G., Pereira, L. S., Raes, D., and Smith, M. (1998). “ETc-single crop coefficient (Kc),” in Crop evapotranspiration, (FAO). Available online at: http://www.fao.org/3/x0490e/x0490e00.htm#Contents (accessed May 16, 2021).

Google Scholar

Allen, S. G., Idso, S. B., Kimball, B. A., and Anderson, M. G. (1988). Relationship between growth rate and net photosynthesis of Azolla in ambient and elevated CO2 concentrations. Agric. Ecosyst. Environ. 20, 137–141. doi: 10.1016/0167-8809(88)90153-3