Chemical Potential-Induced Wall State Transitions in Plant Cell Growth

The pH/T duality of acidic pH and temperature (T) action for the growth of grass shoots was examined in order to derive the phenomenological equation of wall properties for living plants. By considering non-meristematic growth as a dynamic series of state transitions (STs) in the extending primary wall, the critical exponents were identified, which exhibit a singular behaviour at a critical temperature, critical pH and critical chemical potential (μ) in the form of four power laws: fπτ∝τβ-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{\pi } \left( \tau \right) \propto \left| \tau \right|^{\beta - 1}$$\end{document}, fτ(π)∝π1-α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{\tau } (\pi ) \propto \left| \pi \right|^{1 - \alpha }$$\end{document}, gμ(τ)∝τ-2-α+2β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{\mu } (\tau ) \propto \left| \tau \right|^{ - 2 - \alpha + 2\beta }$$\end{document} and gτ(μ)∝μ2-α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{\tau } (\mu ) \propto \left| \mu \right|^{2 - \alpha }$$\end{document}. The indices α and β are constants, while π and τ represent a reduced pH and reduced temperature, respectively. The convexity relation α + β ≥ 2 for practical pH-based analysis and β ≡ 2 “mean-field” value in microscopic (μ) representation were derived. In this scenario, the magnitude that is decisive is the chemical potential of the H+ ions, which force subsequent STs and growth. Furthermore, observation that the growth rate is generally proportional to the product of the Euler beta functions of T and pH, allowed to determine the hidden content of the Lockhart constant Ф. It turned out that the pH-dependent time evolution equation explains either the monotonic growth or periodic extension that is usually observed—like the one detected in pollen tubes—in a unified account.


Introduction
The procedure for seeking explanations for the different modes of plant growth may include analytical models capable of describing the properties of a wall, cell time evolution, or both. 1 3 number of parameters can be misleading; this approach does not oversimplify the complexity of biological systems, but conversely takes all of them into account by building an outer scaffold that cannot be surpassed. One cannot dismiss the possibility that this is the role of phenomenology in making universal statements, in this case imposing a physical constraint (i.e. EWS) on the extending primary wall.
A commonly held view is that biological systems are complex (and molecular biologists are deluged with data), and that therefore modelling has to cope with the fact that the degrees of freedom are numerous, and correspondingly, the number of parameters should be high. The associations that exist between these parameters may in turn cause the emergent properties due to those relationships. Though physical systems on a microscopic scale are also complex and the number of degrees of freedom in real systems is abundant (e.g. Avogadro number, as in the solid-state physics), their behaviour can usually be described by a few parameters (e.g. Ginzburg 2004) within the framework of a phenomenological theory. A similar argument also concerns advanced quantum model(s) (e.g. Hubbard 1963), which has only two parameters: t-for the hopping integral and U-for the Coulomb interaction of electrons within narrow energy bands. In what follows, it is argued that, similar to physical systems, a low number of relevant parameters can also be sufficient to describe complex biological systems. To support this view, let me quote Portes et al. (2015): "This is not to invalidate the strive for detailed models, which are necessary to understand the role of specific components of interest, however, they require a compatible amount of data. In contrast, more general models help to unveil the fundamental principles driving a phenomenon of interest, helping to identify the key regulatory interactions, albeit lacking details". In what follows, the second route was chosen.
The growth of plant cells and organs (like extending cylindrical organs such as pollen tubes, root hairs or grass shoots) are affected by light and humidity, which along with pH and temperature are the most important factors that influence growth. It can be a rewarding enterprise to understand plant growth in terms of the time evolution of STs (e.g. from an ordered to a disordered state, from a stressed to a relaxed state, from a "covalent bonds state" to a "disruption of covalent bonds state", etc.) that take place in the peripheral cell wall. Changes in pH (Ortega et al. 2015;Rayle and Cleland 1992) have been shown to relax the wall stresses presumably by disrupting the load-bearing bonds. Such a transition ('state transition', ST, or equivalently 'state change', which should not be confused with the phase transition that is encountered in physics, since the system is not at an equilibrium steady state) occurs via the exchange of particles and energy-consuming metabolic processes, such as energyconserving ATP, with the inside of cell compartment-the vacuole and cytoplasm, which are treated as the reservoirs of molecules and heat (thermostat). Hence, living organisms were considered to be open systems that dissipate energy and exchange entropy (heat) and matter (molecules) with the environment (Barbacci et al. 2015), a thermostat. This process is accomplished by the presence of pH-and T-dependent Euler beta functions in model equations, which represent the modes of interaction (like heat and mass transfer) with the outside world of a growing plant. Even though I considered a single ST, a cascade (ratchet) of subsequent chemical potential-induced dynamic STs, which may constitute the extension of the cell wall and growth, can be imagined. The prevailing role of the chemical potential was demonstrated and it was proposed that the growth of plant cells and organs is regulated, at the molecular level, by the chemical potential (Matlak et al. 2004;van der Marel 2004) of H + ions. I also show, at a phenomenological level, that the ionic fluxes may generate power-law relations. The proposed model correctly accounted for the set of experimental results and predicted the growth of grass shoots from pH with a high fidelity. To my knowledge, this research is the first attempt to use a statistical physics methods in calculation of critical exponents to model the cell wall during expansive growth and I believe it can provide critical insights on cell wall dynamics at a molecular level (this behaviour could be for a specific type of plant).
Plant developmental systems have evolved within the universal limitations that are imposed by the plant cell wall (Lintilhac 2014). Modulation of the mechanical properties occurs through the control of the biochemical composition and the degree and nature of interlinking between cell wall polysaccharides (Bidhendi and Geitmann 2016). Plant cells encase themselves in a complex polysaccharide wall (Cosgrove 2005a) and characteristically obtain most of their energy from sunlight via the photosynthesis of the primary chloroplasts. The expansive growth of turgid cells, which is defined as an irreversible increase in cell volume, can be regarded as a physical process that is governed by the mechanical properties of the cell wall and the osmotic properties of the protoplast (Schopfer 2006). The precise biochemical mechanism that regulates the ability of the growth-limiting walls to extend irreversibly under the force of turgor pressure has not yet been identified (Kutschera 2000), but it is correlated with the loosening of cells walls, which increases the cell's susceptibility to expansion.
Growing plant cells characteristically exhibit acid growth (Rayle and Cleland 1970;Hager et al. 1971;Cosgrove 1989), which has been formulated in the form of the "acid growth hypothesis" originally independently proposed by Cleland and Hager (Cleland 1971;Hager et al. 1971;Taiz and Zeiger 2006). The acid growth hypothesis was proposed almost 50 years ago and is widely considered to play a role in expansive cell growth. The acid growth hypothesis postulates that both phytotoxin fusicoccin (FC) and the growth-promoting factor auxin (indole-3-acetic acid, IAA) cause wall loosening and produce the concomitant induction of growth (growth enhancement) through the rapid acidification of the extension-limiting cell wall (Cleland 1973;Kutschera 1994). These enigmatic "wall-loosening processes" are, in fact, minor changes within the polymer network of the extension-limiting walls, i.e. the incorporation of proteins, the enzymatic splitting of polymer backbones or covalent cross-links or the disruption of non-covalent interactions between wall polymers via expansin activity (McQueen-Mason and Cosgrove 1994;Dyson et al. 2012). Molecular studies revealed that expansins are a large family of proteins with 38 members in Arabidopsis thaliana, which can be divided into three subfamilies: α-, β-and γ-expansins (Li et al. 2002). It was shown that individual expansins might be tissue-or process-specific, and that a loss-of-function mutation resulted in a disorder in a specific aspect of plant growth or development, for example, in root hair development (reviewed by Marzec et al. 2015a).
Growth is accomplished through the enlargement of the cell volume owing to water uptake, which maintains the appropriate inside pressure in the vacuole as well as the irreversible extension of the pre-existing (primary) cell wall. Hereafter, an almost constant or slowly varying turgor pressure is assumed, which by definition is a force that is generated by water pushing outward on the plasma membrane and plant cell wall that results in rigidity in a plant cell and the plant as a whole.
Expansive growth is the result of the coupling effects (Barbacci et al. 2013) between mechanical (pressure), thermal (temperature) and chemical energy (pH). The thermal sensitivity of biochemical processes refers to the coupling effect between the thermal and chemical energies. In the search for a plant-specific EWS in growing biological cells or non-meristematic tissues, I considered two state variables, namely 'physiological' temperature T (from 0 to 45 degrees in the Celsius scale) and pH at a constant turgor. Note that pH is not sensu stricto a fundamental physical quantity, and therefore, cannot be treated as a usual intensive variable. Nonetheless, relying on the definition of pH and then considering the chemical potential (Baierlein 2001), the following model is proposed.
Evidence has accumulated that the final goal of auxin action (Steinacher et al. 2012;Lüthen 2015;Majda and Robert 2018) is to activate the plasma membrane (PM) H + -ATPase, a process that excretes H + ions into the cell wall. The auxin-enhanced H + -pumping lowers the cell wall pH, activates pH-sensitive enzymes and proteins such as expansins (Cosgrove 1993), xyloglucans (Fry et al. 1992 or yieldings (Okamoto-Nakazato et al. 2001) within the wall and initiates cell wall loosening and extension growth. It has also been observed that when auxin-depleted segments were submerged in an acid buffer solution, the segments started their elongation growth immediately ("acid growth"), whereas the auxin-induced growth ("auxin growth") began after a delay (lag phase), as is shown in Fig. 1 (curve 1) in Hager (2003); see also the remaining plots, which are of quite a different character, (curves 2-6) that correspond to "acid growth". From the many investigations that have been done for more than four decades, it has been deduced that protons, which are exerted into the wall compartment, are directly responsible for the wall-loosening processes through the hydrolysis of covalent bonds, transglycosylation or the disruption of non-covalent bonds (thus changing the state of the system). The growth effect is illustrated in Fig. 2 in Hager (2003), where the "Zuwachs" (increments in %) is plotted against pH.
However, pH can be expressed in terms of a more elementary quantity, namely the chemical potential μ H+ . Note that indirect measurements of the chemical potential by means of pH can be compared to previous measurements of the electromotive force (EMF) in order to detect phase transitions, which are localised by the kinks in the chemical potential in condensed matter physics ). Although pH is not a typical generalised coordinate, the microscopic state of the system can be expressed through it in a collective way.
Temperature is among the most important environmental factors that determine plant growth and cell wall yielding. Plants are responsive to temperature and some species can distinguish differences of 1 °C. In Arabidopsis, a warmer temperature accelerates flowering and increases elongation growth (Jung et al. 2016). In many cases, the process of showing the optimum conditions for plant cell/organ growth. Both coordinates, temperature (T) and acidity or basicity (alkaline medium) pH, are rescaled to unity. The best conditions for maximum growth (green) are found at the saddle point (a stationary point). The simulation parameters that were used are similar to those that were used for auxin-induced growth-see SI Table 1 (Color figure online) growth can be differentiated by the response of plants to temperature (Went 1953). Actually, only a few papers in which the temperature response is treated as a major problem can be mentioned (Yan and Hunt 1999 and papers cited therein), but some recent publications have highlighted the role of temperature as a factor that regulates plant development, i.e. via phytohormones (Dockter et al. 2014;Jung et al. 2016).
Inspired by the Ansatz using Euler beta function f (x, α, β) for growth empirical data, I will show that the final outcome of temperature and pH on plant growth is effectively similar, even though both of the triggers of these responses are apparently of a different nature.

Time Evolution Equation of Plant Cells
Next, the question was posed of whether the cell volume extension or elongation growth of plants (represented by expanding volume V) can be adequately predicted by pH, temperature and pressure.
Modulation of the mechanical properties during expansive growth is still a 'hot topic' for the plant cell growth community, e.g. recently in Sridhar et al. (2018), Bidhendi and Geitmann (2016), Boudon et al. (2015), Barbacci et al. (2013) and Rojas et al. (2011), to name a few. Plant cell expansion, in turn, is controlled by the balance between the intracellular turgor pressure, cell wall synthesis and stress relaxation. Several mainstream explanations of the mechanical aspects of cell growth have been considered in publications, beginning with the cell wall model (Holdaway-Clarke and Hepler 2003) through the instability model (Wei and Lintilhac 2007) and the hydrodynamic model (Zonia and Munnik, 2007) to the statistical model of cell wall dynamics during expansive growth (Sridhar et al. 2018). I propose that these apparently conflicting visions of cell wall extension and growth may be consolidated in a unifying chemical potential-based explanation, which converges in a coherent picture of this somewhat elusive topic.
The last few years have witnessed significant progress in uncovering the molecular basis that underlies plant cell elongation and expansion. In many cases, the first step in cell growth is related to a restriction of the symplasmic communication between neighbouring cells, which permits a specialisation programme that is specific for the individual cell to begin (Marzec and Kurczynska 2014). This isolation enables cells from multicellular organisms to be considered to be individual cells. Afterwards, different signalling and developmental mechanisms are run in order to induce either the expansion of the entire cell or localised elongation. The roles of different protein classes have been described based on the analysis of the growth of pollen or root hair tubes, which are the models that are used in investigations of tip growth (Higashiyama and Takeuchi 2015;Marzec et al. 2015a). Among them, the function of expansins, the proteins that are involved in plant cell wall loosening, which is related to cell wall acidification, has been confirmed in both monocot and dicot species. Analysis of different mutants in the genes encoding expansins confirmed their role based on the analysis of plant phenotypes (Cosgrove 2015). Nevertheless, a mathematical model that permits the cell status  Table 1) of mutant and wild-type plants to be compared might prove useful for further investigations.
According to the acid growth hypothesis, auxin plays a pivotal role in the growth of plant cells (Hager 2003;Lüthen 2015;Majda and Robert 2018 for review). However, recent reports have indicated that a new class of plant hormones, strigolactones (SLs), cooperate with auxin in many developmental processes such as shoot and root branching under normal and stress conditions (Marzec et al. 2013b). Additionally, new functions of SLs, which indicate that those hormones may play a key role in the regulation of plant growth and development, have been proposed (Marzec and Muszyńska 2015). Until now, most data have indicated that SLs influence auxin transport via the distribution of PIN proteins (Shinohara et al. 2013;, but it has also been proposed that SLs influence secondary growth in plants (Agusti et al. 2011) and cell elongation (Hu et al. 2010;Koltai et al. 2010). Mathematical formulas that permit cell growth under different treatments to be measured and described may shed new light on the hormonal regulation of plant cell elongation.
From the physical point of view, the basic ingredients of an expanding plant cell are the cell wall, plasma membrane and cytoplasm. Together, they form a system that can be described as a two-compartment structure with the plasma membrane being treated as the interface between the cytoplasm and the wall. On the other hand, the intensive quantities can be attributed to the system of volume V, namely pressure, and the chemical potential, the latter of which is directly related to a measurable quantity-pH. The problem of the volumetric growth of plant cells, which lies at the intersection of physics, biochemistry and composite materials of a nanometric scale, is inherently difficult to capture in a single theoretical frame (Breidwood et al. 2014). Therefore, the goal was to deconstruct all of the above-mentioned approaches in order to assemble a model that is based on physical principles that are also able to admit other perspectives. It is believed that the EWS that was proposed above has the ability to do this, although its relevance to many of the growth details may be beyond the scope of this conjecture.

Chemical Wall-Loosening Approach
It is believed that the cell wall has active mechanisms for self regulation . The cell wall is formed through the activity of the cytoplasm and maintains a close relationship with it. Communication occurs between the cell wall and the cytoplasm and the plasma membrane occupies a pivotal position in transmitting particles and information. In this model, when applied to pollen tube growth, cell growth is ultimately dependent on the biochemical modification of the wall at the apical part where exocytosis is believed to occur (Holdaway-Clarke and Hepler 2003). Cell wall loosening then allows turgorinduced stretching, thereby pulling the flow of water into the cell (down its potential gradient). During oscillatory growth, enzymes that interfere with wall loosening are believed to be periodically activated or inhibited, which is supposed to drive oscillations in growth and the growth rate (as commented by Zonia 2010). In this approach, the passive role of turgor pressure is stressed.

Loss of Stability Model
A new look at the physics of cell wall behaviour during plant cell growth, which focuses on the purely mechanical aspects of extension growth such as wall stresses, strains and cell geometry, was investigated within the framework of a model that was based on the Eulerian concept of instability Lintilhac 2003, 2007). It was demonstrated that loss of stability (LOS) and consequent growth are the inevitable result of the increasing internal pressure in a cylindrical cell (such as the internodal cylindrical cells of intact Chara corallina plants) once a critical level of pressure-induced stress (P c ) is reached. This model apparently challenged the viscoelastic or creep-based models (Ray et al. 1972;Cosgrove 1993;Schopfer 2006), which were designated as chemical wall loosening and was questioned by Schopfer (2008), who asserted that stress relaxation can be attributed to chemorheological changes in the load-bearing polymer network that permit the plastic deformation of wall dimensions. In his words, the "growth of turgid cells is initiated and maintained by chemical modifications of the cell wall (wall loosening) followed by mechanical stress relaxation generating a driving force for osmotic water uptake". In response to Schopfer's letter, Wei and Lintilhac (2008) claimed that LOS did not challenge the principles of osmotic water relations or biochemically mediated wall loosening and emphasised that plant cells are both osmometers and pressure vessels, and therefore, both facts should be accommodated in any meaningful model. However, the problem of the water potential difference ΔΨ = Ψ o − Ψ i of a growing cell and a "mysterious pump" if pressures must rise to drive LOS behaviour still remained. Regardless of whether the new wall synthesis is homogeneous or patchy (the cell wall can always be conceptually shrunk to a patch), it will affect cell behaviour by modifying the turgor pressure in short time scale oscillations δP (Pietruszka and Haduch-Sendecka 2015a), which are specifically enhanced by the flow of the osmotic pressure through the entire cell volume (beyond the apex as well).

Hydrodynamic Model
Recent works have revealed an important role for osmotic pressure and the hydraulic features of the system in driving cell shape restructuring and growth (Proseus et al. 2000;Proseus and Boyer 2006;Zonia et al. 2006Zonia et al. , 2010Haduch-Sendecka et al. 2014). This concept was established by works that revealed an important role for hydrodynamics in pollen tube growth. The main message that hydrodynamics have the potential to integrate and synchronise the function of the broader signalling network (in pollen tubes) was expressed by Zonia (2010). Pascal's principle states that pressure exerted anywhere in a confined incompressible fluid (cell sap) is transmitted in all directions equally throughout the fluid in such a way that the pressure variations (initial differences) remain the same. Indeed, according to Pascal's principle, turgor pressure is a fast (propagation equal to the velocity of sound in the fluid-soap-medium) and isotropic messenger in the time-evolving cell. Moreover, hypotonic or hypertonic treatment of the periodically elongating cells of pollen tubes led to the recognition that the growth rate oscillation frequency was "doubled" or "halved" with respect to the unperturbed (isotonic) mode (Zonia 2010). The oscillation frequency in such a system has been determined (Pietruszka 2013) to obey the following relation: ω ~ √P, which means that the oscillation growth angular frequency is directly related to the pressure.
From these brief accounts, the question arises-What drives plant cell growth-cell wall loosening or osmotic pressure? Even though there is consent about the role of cell wall properties and turgor pressure in growth, a discrepancy about what drives the initial episode of cell extension exists (Zonia and Munnik 2007). One model proposes that cell wall loosening occurs first and permits the expansion of the cell due to decreased wall rigidity. There is considerable evidence to support this view (Cosgrove 2005b). However, cytomechanical studies have revealed that the viscoelastic properties of the pollen tube cell wall in the apical growth zone do not change considerably (Geitmann and Parre 2004) during growth oscillations. Hence, it was suggested that the cycles of cell wall loosening cannot be the driving force of oscillatory growth (Zonia and Munnik 2007). Moreover, it was shown that increased pressure induces the expansion of the cell wall and drives growth in Chara cells (Proseus et al. 2000;Proseus and Boyer 2006). An intermediate LOS model was put forward (Wei and Lintilhac 2003;Wei et al. 2006) that proposed a gradual increase of turgor pressure to a critical point at which the loss of stability initiates wall extension and growth.

Statistical Model
Global biophysical models such as those of Lockhart and Ortega have provided crucial macroscopic understanding of the expansive growth process, but they lack the connection to molecular processes that trigger network rearrangements in the wall. Sridhar et al. (2018) attempted to fill in this gap by introducing the statistics of tethers connected to microfibrils constituting the mechanical scaffold of the extending primary wall. The tethers are governed by the Poisson point process that has the property that each point is stochastically independent to all the other points in the process. The main advantage of this approach is a consistency with the previous findings (expressed by the Ortega equation), while the disadvantage is the lack of temperature, decisive for the development of plants. The expansive growth, as presented by the authors is a passive process, forced by pressure. Nonetheless, the presence of viscoplastic and viscoelastic terms mapped onto a statistics of independently acting elastic springs (tethers) is very appealing. The detachment and reattachment events causing irreversible deformation of the expanding cell wall seem very convincing in the first order of approximation (independent events, tethers acting as a parallel assembly of springs). The tether properties are described by molecular scale parameters and, in principle, can be obtained from experiments. Interestingly, the calculated relaxation half-time t 1/2 is identical to the one obtained in Pietruszka (2011).
The statistical model can be compared with the concept of STs taking place in the extending primary wall. The latter, expressed in terms of this work, can occur at subsequent time slices, where the tethers statistically detach causing cell expansion. Indeed, the authors state that "wall length arise due to discrete events of tether detachment". If the statistics of this events was big enough causing a ST, a step in length increase can arise leading to pollen tube oscillations.
Without going further into the details of the ongoing discussion (e.g. Zonia and Munnik 2007;Winship et al. 2010Winship et al. , 2011Sridhar et al. 2018) and bringing together these different opinions, I propose a chemical potential based scenario. In what follows, I have intentionally omitted irrelevant details in order to make the main steps clear when introducing it and to build a story with a logical flow.

Plant Material
In order to check the equations, limited experiments were performed using 3-day-old seedlings of Zea mays L cv. Cosmo (Nasiona Kobierzyc. Sp. z o.o. Hodowla roślin rolniczych-Kobierzyce, Poland; courtesy of my PhD student Monika Olszewska). The maize seeds were soaked in tap water for 2 h, and then sown in moist lignin. The seedlings were grown in an incubator in darkness at 27 ± 0.5 °C for 3 days, because this temperature has been defined to grow maize the fastest. Ten mm segments were cut from the plants 5 mm below the tip and first leaf was removed.

Simultaneous Measurement of Growth and Acidification
The growth of intact seedlings of maize (Zea mays L.) and pH of the incubation medium was simultaneously investigated. The growth rates were measured by applying a non-invasive technique that records time-lapse images of the macroscopic elongation of the coleoptiles, while changes in the pH were monitored using a pH/Ion meter (see Olszewska et al. 2018, Supplementary Information, Fig. S1 for the experimental design). The experiment was performed for artificial pond water (APW: 0.1 mM NaCl, 0.1 mM CaCl 2 , 1 mM KCl; pH 6.5 established by HCl or NaOH) conditions. Thirty ten-mmlong coleoptile segments were prepared and then placed in two identical sets (stacks of 15 coleoptiles) of the elongationmeasuring apparatus in an aerated incubation medium (ibid.). The volume of the incubation medium in the apparatus was 5 ml per glass tube (0.3 ml of the incubation medium per segment). Measurements of the growth and pH were performed for 12.5 h and recorded every 15 min. The pH measurement was carried out using two CPI-501 pH-meters. The accuracy level was 0.002 pH, according to the manufacturer's information. Images of the segments were recorded using a CCD camera (Hama Webcam AC-150). Only one of the two probes in separate tubes was selected for analysis. The total length of coleoptiles in one column was 15 × (the length of a single coleoptile), corresponding to the number of repetitions equal to 15 times. The total length of the 15 coleoptile segments (each initially ten-mm-long) was measured in the ImageJ program (ImageJ program bundled with 64-bit Java 1.8.0) with the accuracy established at the ± 0.1 mm level. The relative elongation was calculated using the formula (l t − l 0 )/l 0 with l for "length" and l 0 = l(t = t 0 ) = 10 mm. Also, it was assumed that the coleoptile segment is a cylinder and therefore the crosssection area is a circle. Under an optical microscope (Motic Microscope RED233 at a × 20 magnification), the radius r of the circle for coleoptile cross-section was measured to be equal 0.344 ± 0.001 mm. Then, using the formula πr 2 , the surface of a circle was calculated to be 0.372 mm 2 . Then, the relative volume was obtained by multiplying the length of the coleoptile segment by the average (n = 15) cross-section area. The temperature during this auxiliary experiment was 25 ± 0.5 °C and was maintained using a water bath at a similar level. All manipulations were performed under dim green light.

Equation of Wall Properties for Plant Cells
Based on the "acid growth hypothesis" and relevant experimental data (Yan and Hunt 1999;Hager 2003;), the pH/T (or μ/T) duality of acidic pH or auxin-induced acidification and temperature for the growth of grass shoots was examined in order to determine the model equation (EWS) for the extending primary wall of living plants. In order to introduce pH and T scaling, let us start from basic definitions and data.

Construction of a Phenomenological Model
In the literature, pH is defined as the decimal logarithm of the reciprocal of the hydrogen ion activity a H+ in a solution pH = − log 10 a H+ = log 10 1/a H+ (Covington et al. 1985). The pH value is a logarithmic measure of H + -activity (the tendency of a solution to take H + ) in aqueous solutions and defines their acidity or alkalinity; a H+ denotes the activity of hydronium ions in units of mol/l. The logarithmic pH scale ranges from 0 to 14.
The temperatures at which most physiological processes normally occur in plants range from approximately 0 °C to 45 °C, which determines the physiological temperature scale (in Kelvin scale: [K] = [°C] + 273.15). The temperature responses of plants include all of the biological processes throughout the biochemical reactions (high Q 10 factor). This is clearly visible in the relative rates of all of the development or growth processes of maize as a function of temperature, which is illustrated in Fig. 3 in Yan and Hunt (1999) and ranges from 0 to 45 °C. Therefore, to be "on the safe side", I propose that the 'physiological' temperature scale for plants [0-45] °C, which corresponds to a [0, 1] interval after rescaling, be setup. Otherwise, diverse temperature scales can be introduced in different situations.
Which data have been adopted for the research and how those data have been handled are described below. For the analysis, I accepted the reliable published data, namely the pH plots that are presented in Hager (2003) and the temperature plots that are presented in Yan and Hunt (1999). These charts were first digitised by a GetData Graph Digitizer (version 2.26.0.20) for reconstruction and then rescaled (see SI Fig. 1). The original temperature range [0, 45] °C was converted to a [0, 1] interval by dividing by 45. The same procedure was conducted on the pH scale [0,14] in order to obtain the normalised [0, 1] interval (SI Fig. 1), which was indispensable for the further calculations. Very low or very high pH values are obviously not physiological at least for higher plants, however, such excessive conditions, which constitute unacceptable extremes for life to come to existence, are excluded by nature. Nonetheless, the [0, 14] pH scale should be mapped as a whole to a 0 to 1 range to establish the normalised scale, free form any arbitrary assignment.
The pH plots that are presented in Fig. 2 in Hager (2003) as well as the temperature plots that are presented in Fig. 3 in Yan and Hunt (1999) look similar when mirror symmetry is applied. Providing that a substitution x → (1 − x) is performed and a function f(x) of variable x and of its reflection (1 − x) is considered, one might expect that they will look approximately the same after the correct scaling. These properties can be adequately described by the Euler beta density distribution (or beta Euler, for short). In this approach, let us assume that x equals either pH or temperature (both rescaled to a [0, 1] interval), an assignment which makes sense since at a relatively low pH, which corresponds to a relatively high temperature in this unit-scaling, plant cells and organs grow (elongate) the fastest, which is clearly reproduced in SI Fig. 1.
The data from SI Fig. 1a were taken separately ("auxin" and "acid" growth; Hager, 2003) for the fitting procedure (SI Table 1 and SI Fig. 2), while the temperature data (Yan and Hunt 1999) from SI Fig. 1b were first normalised and then fitted by beta Euler. Furthermore, when comparing the data from different experiments (Yan and Hunt 1999;Hager 2003), it is observed that the x variable could, after rescaling, have the following meaning-either acidity (basicity) x ≡ pH or temperature x ≡ T (see SI Fig. 1) with a characteristic beta function cut-off at x = 0 (SI Fig. 1a and the inset) and at x = 1 (SI Fig. 1b); the normalised scale that was used is presented in SI Fig. 1b (inset).

Derivation of the EWS for Plants
The probability density function of the beta distribution also called the Euler integral of the first kind (Polyanin and Chernoutsan 2011) for 0 ≤ x ≤ 1 and shape parameters α and β > 0 is a power function of the variable x where B = B(α, β) is the normalisation constant. Based on the experimental observations as discussed in the previous section and presented in SI Fig. 1, the left side of Eq. (1), rewritten twice, first for x = pH, then x = T, can be merged into a single expression with a condition that T and pH are normalised and dimensionless; c pH and c T are constants normalising the maxima of the beta functions to 1, and "pH" is hereafter treated as a non-separable variable name. The lower indices "T" and "pH" in Eq. (2) should be read "at constant temperature" and "at constant pH", respectively, which is valid for the preparation of strict experimental conditions to achieve a stable equilibrium (Lüthen et al. 1990;Hager 2003). The appropriate (i.e. symmetric with respect to x = 0.5 line) values of T (1) and pH should be taken into account in Eq. 2). The equality sign and the equalisation 1 = 2 ≡ , 1 = 2 ≡ as well as c T = c pH in Eq.
(2) seem to be justified for biological systems (see SI Fig. 1.1) in the first order of approximation. Such a modified Eq.
(2) will be considered further as a constraint on the two independent variables T and pH. By defining a new function f at constant turgor pressure P (isobaric conditions: δP ≈ 0) according to that is an analytic function of pH and T on the physically accessible interval [0, 1) for each variable (here function f constructs a bridge between the physical and chemical features of the living system), and assuming adequate water uptake to fulfil this requirement, we determine from Eq. (2) After rearrangement, the constitutive dimensionless EWS for living plant cells reads which has the elegant structure of a double power law. T in Eq. 5 is the temperature in Celsius rescaled to a [0, 1] interval as pH also is. Here, the constants α and β are the shape exponents (SI Table 1 and SI Fig. 2). The saddle-point type of Eq. (3) that is presented in Fig. 1 represents a concave surface in the form of a contour plot, which is roughly accurate at low pH and moderate temperatures. Equation (5) becomes increasingly inaccurate (divergent) at very low or very high pH values (close to zero or one in the normalised scale, which is used throughout of the article), though such extreme conditions are not admitted by nature and can be ignored. The isotherms of Eq. (3) are presented in Fig. 2a, while the lines of constant pH are presented in Fig. 2b.
Even though the derivation was apparently launched from the "acid growth theory", a closer look at the model assumptions revealed that, in fact, it was based on the typical raw experimental data of pH-and T-dependent growth (SI Fig. 1). For that reason, Eq. (5) is independent of whether the acid growth theory applies or not, and therefore, seems universal, at least for tip-growing grass shoots. Note, Eq. (5) is not an evolution (growth) equation, but a system (cell wall) property equation. The solutions of Eq. (5) for α and β in the form of contour plots are presented in Fig. 3.
The essential relations are property relations and for that reason they are independent of the type of process. In other words, they are valid for any substance (here: the primary wall of a given plant species) that goes through any process or mode of extension. Yet, in the present case, we should be aware of the empirical origin of Eq. (5) and remember the fact that a fundamental variable of the system is the chemical potential μ at a constant pressure and temperature (here: μ = μ H+ (T) for H + ions). By taking the total derivative, where the implicit μ H+ (T) dependence is taken into account, we get where In approximation, since pH is usually measured in isothermal conditions (SI Fig. 1a), let us temporarily suspend the pH(T) dependence in Eq. (6a). Equations (3) and (6) describe the mathematical relations between pH, T and μ H+ variables.
Equation (5) is the EWS that describes the plant cell wall properties, where the SFs are the scaled temperature and pH of the expanding cell volume. Here, pH, which is intimately connected with the acid growth hypothesis, introduces a direct link between the physical variables that represent the state of the system (growing cell or tissue) and the biological response (via pH altering-H + -ATPase), which results in growth. Clearly, Eq. (5) can be further verified by the diverse experiments that have been conducted for many species in order to determine the characteristic triads (α, β, Ω S ) that belong to different species or families (classes) taxonomically, or that are simply controlled by growth factors (growth stimulators such as auxin, fusicoccin or inhibitors such as CdCl 2 ). However, the real test of the utility of Eq. (5) would include the predictive outputs for any perturbations that could then be subject to experimental validation. From the dyads (α, β), Eq. (5) should also be able to connect with either the underlying physiological processes or with their molecular drivers. In particular, whether the underlying molecular mechanism is identical or different should be reflected in the critical exponents (α and β pair should be identical for the same molecular mechanism and dissimilar for diverse mechanism). Despite the fact that f is not a generalised homogeneous function, it can also be used to check for a kind of "universality hypothesis" (Stanley 1971) and to validate the evolutionary paradigm with genuine numbers.

Calculating Critical Exponents
Analytical evaluation must lead to predictions that can be used as a Litmus test in order to cut through the myriad of conceptual models that have been published in the biophysical community. In plant science, given that the calculation of critical exponents can be connected with the optimum growth, it can be helpful to identify the molecular mechanism that is responsible for growth (in principle, similar exponents should be connected with the same mechanism). It could also be a viable method when growth is stimulated (laser, UV, electromagnetic field, etc.). In physics, phase transitions occur in the thermodynamic limit at a certain temperature, which is called the critical temperature, T c , where the entire system is correlated (the radius of coherence ξ becomes infinite at T = T c ). Here, I want to describe the behaviour of the function f(pH, T), which is expressed in terms of a double power law by Eq. (3) that is close to the critical (optimum) temperature and, specifically in the considered case for practical reasons, also about the critical (optimum) pH = pH c . Critical exponents describe the behaviour of physical quantities near continuous phase transitions. It is believed, though not proven, that they are universal, i.e. they do not depend on the details of the system. Per analogia, since growth (such as cell wall extension in pollen tubes or the elongation growth of the shoots of grasses, coleoptiles or hypocotyls) can be imagined as advancing over the course of time, a series of non-equilibrium STs (see "Discussion"), the critical exponents can be calculated.
Lets introduce the dimensionless control parameter (τ) for a reduced temperature and similarly (π), for a reduced pH which are both zero at the transition, and calculate the critical exponents. It is important to keep in mind that (in physics) critical exponents represent the asymptotic behaviour at a phase transition, thus offering unique information about how the system (here: the polymer built cell wall) approaches a critical point. A critical point is defined as a point at which ξ = ∞ so in this sense T = T c and pH = pH c are critical points of Eq. (5). This property can help to determine the peculiar microscopic mechanism(s) that allows for wall extension (mode of extension) and growth in future research. The information that is gained from this asymptotic behaviour may advance our present knowledge of these processes and their mechanisms and help to establish the prevailing one. By substituting Eqs. (7) and (8) into (3), one can calculate the critical exponents (when T → T c and pH → pH c ) from the definition (Stanley 1971), which will give us and here λ 1 = β − 1 and λ 2 = 1 − α are the critical exponents (λ 3 = 0). The above equations result in two power relations that are valid in the immediate vicinity of the critical points where lower indices denote constant magnitudes. Equations (12) and (13) represent the asymptotic behaviour of the function f(τ) as → 0 or f(π) as → 0 . In fact, we can observe a singular (divergent) behaviour at T = T c (Fig. 4a) and pH = pH c (Fig. 4b). At this bi-critical point, the following relation holds: α + β = 2. Thus, the exponents are related through scaling relations, in analogy with the static critical phenomena.
Retaining τ for a reduced temperature, let us introduce another-microscopic-control parameter, μ, for the chemical potential Substituting τ and Eq. (14) into Eq. (6b) for g(T, μ) and assuming the reference potential μ 0 = 0, one can calculate the critical exponents (when T → T c and H + → c ) to get The above limits result in two more power relations, which are valid in the immediate vicinity of the critical points Equations (18) and (19) represent the asymptotic behaviour of the function g(τ, μ) as → 0 or g(τ, μ) as → 0 . At the Fig. 4 Critical behaviour at the STs that take place in the cell wall (pH, T-variables). a A "kink" (discontinuity in the first derivative) at the critical temperature T = T c . b a "log-divergent" ("λ"-type) solution at the critical pH = pH c . Control parameters: reduced temperature (τ) and reduced pH (π). Simulation parameters: a β = 1.93 and b α = 1.7 (SI Table 1). This could be representative of distinct wallloosening mechanisms microscopic level, we observe a singular behaviour at T = T c (Fig. 5a) and μ = μ c (Fig. 5b, c).
For the bi-critical point, the exact relation β = 2 holds (compare with SI Table 1 and SI Fig. 2), thus leaving us with the only free parameter (α) of the model. Note the broad peak at the critical temperature (indicating that optimal growth occurs within a certain temperature range) and a sharp, more pointed peak for the critical chemical potential. In this representation, growth can be treated as a series of subsequent STs that take place in the cell wall. It looks as though a cardinal quantity that is decisive for growth in this particular approach is the chemical potential.
Complex systems such as the cell wall of a growing plant are solid and liquid-like at the same time. At the molecular level, they are both ordered and disordered (allowing for the transition between these states). There are obviously complex interactions between the model parameters and inputs from many other kinds of inputs. However, the EWS describes the properties of the wall and unknown complex interactions in a general way. Whether the wall structure is ordered or disordered, universal features can be identified using simple scaling laws.

Time Evolution Equation of Plant Cells
If one imagines that the growth rate, as a function of pH, has a (Euler) beta function form for any fixed temperature (T) and a similar dependence as a function of T for a fixed pH (see SI Table 1 and SI Fig. 1 and SI Fig. 1.1), one can conclude that the volumetric growth rate is proportional to the product of the beta functions of T and pH (and to the effective turgor pressure P − P c ). This assertion can be the starting point for the heuristic derivation, in the first approximation in which one can neglect the nested pH(μ H+ (T)) dependence, of the equation below.

Mystery of the Lockhart Coefficient Explained?
Simply speaking, the above statement for the relative growth rate can be expressed in a way that is analogous to the Lockhart (1965) equation as an Ansatz ("educated guess") and can be postulated in the general form Fig. 5 Critical behaviour at the STs that take place in the cell wall (for μ, T-variables). a A "log-divergent" ("λ"-type) solution at T = T c . b-c Solution at the critical value of the chemical potential μ = μ c ; b "acid growth", c "auxin growth" (SI Table 1). Control parameters: reduced temperature (τ) and chemical potential (μ). Simulation parameters: a α = 3.52, β = 1.93; b α = 1.87, β = 3.17; c α = 1.7 and β = 3.52 (SI Table 1) at a constant turgor pressure P and yield threshold P c (both in [MPa]). Here, c is a scaling constant that is proportional to the reciprocate of the product of the beta function normalisation constants, see Eq.
(1) and can be incorporated to Ф 0 that was introduced for dimensionality reasons: [Ф 0 ] = MPa −1 × s −1 , as in the Lockhart equation. In Eq. (20) I assume concomitant water uptake; note that the P − P c term can be substituted with a more appropriate turgor plus osmotic pressure term in order to allow for water transport (via aquaporins) due to the changes in the osmotic and water potential (Obermayer, 2017 in Pollen Tip Growth)). V = V(t) is the increasing volume of a growing cell over the course of time. Here, "pH" is treated, as before, as a non-separable variable name. In isothermal conditions, the integration of Eq. (20) yields the solution for the enlarging volume in the form where 0 and t denote the initial and final time. Note that, unlike the Lockhart equation, the 'acidification' that is present in Eq. (21), which is represented by the value of dynamic pH, prevents a plant cell from unlimited exponential growth (burst). In contrary, the solution of the Lockhart equation V(t) = V(0)exp(Ф 0 (P − P c ))t means that plants grow with no limitation, which is an apparent contradiction to real-life phenomena. Equation (21) correctly delivers not only the initial phase of slow, and next of rapid growth, but also the deceleration and saturation phases (Fig. 6). Equation (21) holds for dynamic pH-in Fig. 7 changes in the volume and saturation effect depend on the simulation parameters. The model suggests that the yielding behaviour encapsulated in the classical Lockhart equation can be explained by the strongly non-linear dependence of V on pH and T. By recalling the canonical form of the Lockhart equation with the extensibility Ф controlling the rate of growth and comparing Eqs. (20) and (22), one obtains the pHdependent estimate for the Lockhart constant: where the exponents α and β are numbers that are independent of pH (or μ) and T, whereas pH = pH(t). Note that Eq. (20) turns into its canonical-Lockhart form, Eq. (23), for α = β = 1 and c = 1. Hence, the Lockhart coefficient Ф, which is responsible for viscoplastic extensibility, is an analytic function of pH and T on the physically accessible interval from 0 to 1 for each variable. However, the augmented form of the Lockhart equation, which is expressed by Eq. (20), provides the possibility to not only describe monotonic cell expansion correctly (Figs. 6, 7), but also the oscillatory mode of the extension of the pollen tubes. Apparently, especially in the case of oscillatory or just changing turgor pressure, the additional (viscoelastic) term 1/ε × dP/ dt can be introduced (added) into Eq. (20) on the right side, to become the Ortega-like equation (Ortega, 1985), for completeness. However, another proposal may include a variability of P in agreement with pH changes (see the inset in Fig. 9) in a self-consistent solution of Eq. (20). This, in turn, can lead to the growth rate oscillations that are observed in pollen tubes. Additionally, the remarks concerning the possible changeability of P c (= Y, a yield threshold) during Fig. 6 a Growth rate plot obtained for Eq. (20) where pH was assumed to be inversely proportional to time (Lüthen et al. 1990, Fig. 6). Simulation data: α = 3, β = 2; arbitrary units. b Typical growth curve obtained from Eq. (21)-all growth phases are clearly visible: initial phase of slow growth, accelerated growth, quasi-linear phase, deceleration and saturation (Fogg, 1975). See also Fig. 7 for parameterisation of Eq. (21) with respect to pH (7a), temperature (7b) and turgor pressure (7c) growth (Pietruszka 2012) can be relevant in the higher order of approximation. Moreover, the periodically elevated turgor pressure P by a small δP(t) amount (Benkert et al. 1997) can reflect the viscoelastic/viscoplastic dynamic transitions caused by minute (but many) pH changes similar to the one presented in the inset to Fig. 9-leading to macroscopic pollen tube oscillations. Note that Eq. (20) guarantee that the maximum growth rate is for the maximum of beta f(α, β, T) function. It is also evident that oscillations of pH ≈ sin(Ω·t) will produce oscillations of the exponent function, however, shifted in phase like cos(Ω·t). This, in turn, can be responsible for-as yet unexplained-about 90° phase shift observed by many authors (e.g. Messerli et al. 1999) in pollen tubes.
Calculation of the α and β shape coefficients for the actual empirical data is commented below.

Determination of the α and β Shape Coefficients for Experimental Data
Potential applications of the presented formalism lie in the α and β parameters and the direct use of pH(t) data in calculating the volume increase dV. However, in Eq. (21), one  Lüthen et al. (1990). Squares and triangles-converted (recalculated) separately (point by point) using Eq. (21) from the digitised data presented in Fig. 6 by Lüthen et al. (1990); dashed curves-sigmoid curve (Boltzmann) fit (Microcal Origin). Conversion data: α = 3 and β = 2. Simulation parameters: T = 0.62 (~ 26 °C) and P -P c = 0.3 MPa obtains exp(f(α, β, c) × integral (α, β)) dependence-such problems have several local extremes. Therefore, the solution of Eq. (20) belongs to the class of highly non-linear optimisation problems and will be presented in a separate paper.
I suppose that the biological meaning of the α and β shape exponents that is not discussed here may acquire connotation for more abundant sets of data. Nonetheless, it seems that Eq. (21) can be treated as the transformation pH(t), P, T , ⟶ V(t) for the increasing volume of a plant cell (corresponding to the pH-driven cell enlargement).

Monotonic Cell Elongation
In order to check the Ansatz, which is expressed by Eq. (20), I also examined the decreases of pH that were induced by auxin (IAA) and fusicoccin (FC) in Lüthen et al. (1990), which is presented in Fig. 6 therein. There, the additions of IAA and FC resulted in rapid drops of pH. The authors asserted that these lower pH values are responsible for the rapid increase in the growth rate. I digitised the pH data that were presented by Lüthen et al. (1990) and inserted these (point by point) into Eq. (21) separately. As a result, I received two fragments of growth curves, which are presented in Fig. 8. Note that to determine the cell volume evolution over time, we only needed the values of the following dynamic variables: pH, temperature and turgor pressure; the initial volume V 0 was needed for scaling; T = const. and P ≈ const. though the latter may also change in time. Interestingly, the simulation data for the α and β parameters, which represent the properties of the wall, were adapted by combining the EWS with the evolution Eq. (21). These data, for which Eq. (21) is not divergent (unlike the Lockhart equation), closed the "system of equations" for growth. By the "system", I mean (a) the EWS, which describes the properties of the cell wall irrespective of the extension mechanism and (b) the equation of time evolution, Eq. (21), which describes the dynamics of the growing system at temperature T (see also Pietruszka, 2012). Moreover, the direct use of Eq. (21) to our own experimental data is presented in Fig. 9, where pH(t), P and T transform directly into the elongating height (volume V) of a plant tissue (note the correctly determined initial value of 1).
In this paragraph, I demonstrated that Eq. (21) correctly describes monotonic growth in the context of the "acid growth hypothesis" (note that pH is explicitly included in this equation; see also Figs. 8 and 9 and Eq. (9) in Pietruszka and Haduch-Sendecka (2016) and the conclusions therein), which validates its applicability not only for fusicoccin, but also for auxin (Hager 2003). The high accuracy of the theoretical fits (solid lines) with those that originated experimentally (Fig. 6 in Lüthen et al. 1990) or especially to the data, which is presented in Figs. 8 and 9, respectively, strongly support the presented approach. However, the question about the potential applicability of Eq. (21) to oscillatory growth still remains.

Periodic Cell Elongation
The pH value is a logarithmic measure of the tendency of a solution to take H + in aqueous solutions. Yet, pH can be directly expressed in terms of the chemical potential (μ H+ ) through the relation where R denotes the gas constant and 0 H + is the reference potential. pH is not a typical generalised coordinate, although the microscopic state of the system can be expressed through it in a collective manner. Just as temperature determines the 1 3 flow of energy, the chemical potential determines the diffusion of particles (Baierlein 2001).
After inserting Eqs. (24) into (21), one can calculate the influence of the chemical potential oscillations of H + ions, μ H+ = μ H+ (t), on the volumetric expansion. A simplistic numerical simulation reveals that the volumetric expansion represents a typical growth (sigmoid) curve (like the one presented in Fig. 6b). However, in this case, when superimposed onto the simulated oscillations of the chemical potential, they produce the growth rate oscillations in a straightforward manner (Pietruszka and Haduch-Sendecka 2015a). This type of oscillations was observed in the rapid expansion of pollen tubes (Plyushch et al. 1995;Hepler et al. 2001;Feijo et al. 2001), and the source of this periodic motion has been the subject of a long-lasting debate among plant biologists. Portes et al. (2015) noted that "Since the ultradian rhythms present in pollen tubes are not temporally coincident with protein synthesis and degradation, a simpler mechanism might be triggering the overall oscillation". Further, the authors claimed that "the strongest candidates seem to be the membrane trafficking machinery and ion dynamics, since they are the only candidate mechanisms that may oscillate independently of tube extension (Parton et al. 2003)". Indeed, in what follows a chemical potentialinduced (ultimately connected with ion dynamics since it measures the tendency of particles to diffuse) recurrent model, which triggers the overall oscillation, is proposed to explain this oscillatory mode of a single cell extension.

Discussion
The physicist Erwin Schrödinger (1944) may have been the first to consider the thermodynamic constraints within which life evolves, thereby raising fundamental questions about organism evolution and development. Here, it was attempted to present a more modest description in terms of the EWS and methods such as the calculations of critical exponents. More precisely, this work seeks to model the relationship between plant growth and temperature (or pH), thereby describing growth optima as functions of the state variables and critical exponents. Apparently, there is already a simple and compelling explanation for the growth optima of plant growth from the biological point of view. Growth is mediated both directly and indirectly by enzymes and increases broadly in line with the frequently observed effect of temperature on enzyme activities (Berg et al. 2002). However, once temperatures are sufficient to denature proteins, enzyme activities rapidly decrease and plasma membranes and some other biological structures also experience damage (compare SI Fig. 1b at a high temperature end).

Sensing Structural Transitions
A magnitude that infinitesimal changes result in symmetry change always exists in physics in continuous phase transitions (Landau and Lifshitz 1980). However, structural changes such as the ones considered here usually take place in non-equilibrium conditions (Racz 2002). Seemingly, such (order/disorder) symmetry changes (De Gennes 1991) in biological systems may be connected with a mechanism (Rojas et al. 2011) by which chemically mediated deposition causes the turnover of cell wall cross-links, thereby facilitating mechanical deformation. It may also reflect the pectate structure and distortion that was suggested by Boyer (2009) (Fig. 5c, d) in Chara cell walls by placing wall polymers in tension and making the load-bearing bonds susceptible to calcium loss and allowing polymer slippage that irreversibly deforms the wall. Proseus and Boyer (2007) suggested that a ladderlike structure would be susceptible to distortion and that the distortion would increase the distance between adjacent galacturonic residues (Fig. 5d in Boyer 2009), and therefore, the bonds may lengthen and thus weaken and decrease their affinity for Ca 2+ . Dissociation may then occur, thus allowing an irreversible turgor-dependent expansion (see also Fig. 4 in Pietruszka 2013; this is also applicable for non-isochronous growth in pollen tubes). The direction of the maximal expansion rate is usually regulated by the direction of the net alignment among cellulose microfibrils, which overcomes the prevailing stress anisotropy (Baskin 2005). However, the organisation of cellulose microfibrils depends on the organisation of the cortical microtubules (Paredez et al. 2006), which is modified by many different factors such as mechanical stimuli, light or hormones (Fishel and Dixit 2013). These modifications are regulated via DELLA proteins (Locascio et al. 2013) and/or arabinogalactan proteins (Willats and Knox 1996;Marzec et al. 2015b). The transient changes in the wall composition and the deformation of bonds may be a hallmark of symmetry change and a kind of ST (orchestrated instability) that is taking place in the cell wall.
The cell wall is supplied with the precursors of polysaccharides and also some proteins through exocytosis. This process is a form of active transport in which a cell transports molecules out of the cytoplasm into the extracellular matrix or intracellular space by expelling them in an energy-using process. Endocytosis is its counterpart. The loaded vesicles are transported by kinezine through microtubules into the plus end while by dinezine towards the minus ends, and ATP is utilised in both binding and unbinding processes of the motors heads. Besides microtubules, the actin cytoskeleton represents the most demanding energy-using process. In fact, in plant cells, both exocytosis and endocytosis are F-actin dependent processes that are related to cell wall assembly and maintenance. Moreover, there is strong evidence that plant cell elongation itself also depends on motors actin cytoskeleton rearrangement, see a series of studies by Baluška et al. (1997Baluška et al. ( ), (2001Baluška et al. ( ), (2002 and Wojtaszek et al. (2007). Indeed, these energy-consuming processes are apparently not included in the EWS.
A situation that was already pointed out in the Introduction in which the protons that are excreted into the wall compartment are directly responsible for the wall-loosening processes through the hydrolysis of covalent bonds or the disruption of non-covalent bonds may also be a signature for the symmetry change and STs that occur in the wall compartment of the growing cell. Though the underlying microscopic mechanisms are not well recognised as yet (and as a result the order parameters are difficult to identify-tethers 'melting' from the attached to detached state against pH or T can be a first candidate, Sridhar et al. 2018), a pH-driven "lambda" ST may be attributed to the maximum activity of PM H + -ATPase, while temperature-driven STs in the cell wall can be directly related to the maximum elevation of the effective diffusion rate (k 2 coefficient in Pietruszka 2012)-it was considered primary (diffusive) growth throughout this article. As an aside, from the calculation of cross-correlations, we may draw further conclusions that are connected with the biochemical picture of the acid growth theory. The result presented in SI Fig. 3 is a clear characteristic that temperature-induced growth and auxin-induced acidification (PM H + -ATPase) growth correlate strikingly well (there is good experimental evidence that a higher temperature induces auxin biosynthesis in some plants). However, the convolution of acidic pH growth and temperature growth is less shifted away from zero (lag) and is even more pronounced; see also recent outcomes in Pietruszka and Haduch-Sendecka (2016) as well as Olszewska et al. (2017Olszewska et al. ( , 2018 for cell wall pH (proton efflux rate) and growth rate, which are directly co-regulated in growing shoot tissue, thus strongly supporting the empirical foundations of EWS. Even more supportive for the EWS-model is the result obtained by Damineli et al. (2017), where the extracellular H + flux and growth rate curves practically overlap (see Fig. 2d therein); the periods (47.7 ± 3.2 s and 48.7 ± 4.3 s, respectively) for both processes are almost identical within experimental error. This strict quantitative result may contribute to the acid growth hypothesis that was developed by plant physiologists and make a contribution to resolving the problems that occur in the long-lasting discussion. However, a problem still arises. It is unclear to what degree Eq. (5) applies universally as the published raw data that constitute the experimental basis are only for grass shoot elongation growth. The long-standing acid growth theory postulates that auxin triggers apoplast acidification, thereby activating the cell wall-loosening enzymes that enable cell expansion in shoots. This model remains greatly debated in roots.
However, recent results (Barbez et al. 2017) in investigating apoplastic pH at a cellular resolution show the potential applicability of Eq. (5) for roots.

Migration of Plants
I now show that my findings may be embedded in the evolutionary context that is connected with the migration of plants away from the Equator (with changes in latitudes) as the climate changed as well as their adaptation to the spatial distribution of pH in the soil as a substitute for a high temperature. As an example that exhibits the capabilities of Eq. (5), let us consider the following. The dominant factors that control pH on the European scale are geology (crystalline bedrock) in combination with climate (temperature and precipitation) as was summarised in the GEMAS project account (Fabian et al. 2014). The GEMAS pH maps primarily reflect the natural site conditions on the European scale, while any anthropogenic impact is hardly detectable. The authors state that the results provide a unique set of homogeneous and spatially representative soil pH data for the European continent (note Fig. 5, ibid.). In this context, the EWS that is expressed by Eq. (5) may also have an evolutionary dimension when considering the spatial distribution of pH data as affected by latitude (angular distance from the Equator) elevation on the globe, which is associated with climate reflection (ibid.). pH is strongly influenced by climate and substrate-the pH of the agricultural soils in southern Europe is one pH unit higher than that in northern Europe. See also Fig. 4 in Fabian et al. (2014) in which the pH (CaCl 2 ) of soil samples that are grouped by European climate zones are measured (descending from pH 7 (Mediterranean), pH 6 (temperate) through pH 5.25 (boreal) to pH 4.80 (sub-polar)). This fact may be reinterpreted in terms of the pH-temperature duality that is considered in this work. When the migration of plants away from the Equator occurred, lower pH values further to the North might have acted in a similar way on the growth processes as high temperatures in the tropic zones. This situation could have led to energetically favourable mechanisms (adaptation) such as the amplified ATP-powered H + extrusion into the extending wall of plant cells and the secondary acidification of the environment (released protons decrease pH in the incubation medium), irrespective of the initial pH level, which is the main foundation for the acid growth hypothesis. Apparently, the evolutionary aspect of Eq. (5), which suggests a possible link with self-consistent adaptation processes, as one of its first applications, seems difficult to underestimate, and the potentially predicative power of EWS looks remarkable considering the modest input that was made by deriving Eq. (5). In this perspective, the pH-response could potentially provide a plant with an adaptive advantage under unfavourable climate conditions. Together with the discovery of the function of phytochromes as thermosensors (Jung et al. 2016), the EWS could be useful in breeding crops that are resilient to thermal stress and climate change.

Classes of Universality
The EWS for a biological system has not been reported until this study. To further verify the EWS and its implications, it was noted that based on the results that were obtained (the calculated critical exponents), a kind of universality hypothesis (Stanley 1971) can be tested experimentally in order to extract "classes of universality" for plant species (Ω S ), thereby substantiating the taxonomic divisions in the plant kingdom via quantitative measures (Table 1). This issue, which may be called "extended plant taxonomy" (similar to "dimensionless numbers" proposed by Ortega 2018 for Plant Biology), demands further experimentation and is beyond the scope of this theoretical work.

Further Perspectives
It seems the EWS may explain a number of phenomena and is powerful in its simplicity. It positions itself at the brink of biology and nano-materials. It may be helpful in identifying the processes of the self-assembly of wall polymers. In agricultural implementations, any departures from the optimum growth may simply be corrected by appropriately adjusting pH or T.
The finding of the EWS for plants adds an significant dimension in the biophysical search for a better explanation of growth-related phenomena in a coherent way through its applicability to the viscoelastic (or viscoplastic) monotonically ascending and asymptotically saturated (Pietruszka 2012) cell wall extension that is usually observed, as well as to pollen tube oscillatory growth. The underlying biochemical foundation, which is expressed in the form of the acid growth hypothesis, could also help to understand the growth conditions that are ubiquitous in biological systems through coupling to EWS. In this respect, the presented method is accompanied by its partial validation-its application to the important biological question concerning the acid growth hypothesis. In the above context, it is not surprising that the growth of the cell wall in plants can be thought of as a timedriven series of STs (via the breaking of polymer bonds or other mechanisms that were mentioned earlier), although they are ultimately connected with the physical conditions and chemistry by the Eqs. (5) or (6b). Since non-equilibrium stationary states, which always carry some flux (of energy, particle, momentum, information, etc.), are achieved after discrete intervals of time for single cells such as pollen tubes, such a (chemical potential ratchet) mechanism can lead to a kind of leaps in the macroscopically (of μm scale) observed growth oscillations (Geitmann and Cresti 1998;Hepler et al. 2001;Zonia and Munnik 2007). Similar mechanisms were observed in multicellular models for plant cell differentiation such as rhizodermis (Marzec et al. 2013a) where the growth of root hair tubes occurs by leaps and bounds (Vazquez et al. 2014). This phenomenon may be treated as the final result of the subsequent STs that take place in the cell wall.
On the other hand, spontaneous bond polarisation/breaking may be synchronised through Pascal's principle by pressure P or pressure fluctuations δP (Pietruszka and Haduch-Sendecka 2015a), thereby acting as a long-range (of the order of ξ) messenger at the state change (the radius of coherence → ∞ at a Γ-interface in Pietruszka 2013). When the entire system is correlated and the stress-strain relations are fulfilled (ibid.), the simultaneous extension of the cell wall at the sub-apical region may correspond to pollen tube oscillation(s). Hence, an uncorrelated or weakly correlated extension would result in the bending that is observed in the directional extension of the growing tube. In multi-cell systems (tissues) that have a higher organisation (coleoptiles or hypocotyls in non-meristematic zones), the EWS for plants will manifest itself by the emergent action of an acidic incubation medium or will be induced by endogenous auxin acidification as is presumed in the acid growth theory. In this aspect, the EWS may serve as a new tool for the further investigation and verification of claims about the acid growth hypothesis (Barbez et al. 2017). This should greatly facilitate the analysis of auxin-mediated cell elongation as well as provide insight into the environmental regulation of the auxin metabolism.
It is important to note that "Cellular growth in plants is constrained by cell walls; hence, loosening these structures Table 1 Preliminary data for α and β exponents for maize (Zea mays L.) and wheat (Triticum vulgare Vill.), which belong to the same family Poaceae but to different subfamilies: Panicoideae and Pooideae, respectively; the calculated values for α and β are similar The common bean (Phaseolus vulgaris L.) and cucurbita (Cucurbita pepo L.) belong to the same clade Eudicots but to different orders: Fabales and Cucurbitales, respectively; the calculated values for α and β differ significantly. The calculation that is based on Fig. 1 in Lewicka and Pietruszka (2008) and SI Eq. (1). is required for growth. The longstanding acid growth theory links auxin signalling, apoplastic pH homeostasis, and cellular expansion, providing a conceptual framework for cell expansion in plant shoots. Intriguingly, this model remains heavily debated for roots" (Barbez et al. 2017). The authors' findings permitted the complex involvement of auxin in root apoplastic pH homeostasis to be elucidated, which is important for root cell expansion, and shed light on the acid growth mechanism in roots. This is, in turn, the outcome of the model equations, which are apparently applicable for elongating roots (root hair), shoots or even a single plant cell (pollen tube). Moreover, by using the EWS, it will be possible to answer the question about role of the crosstalk between auxin and other hormones such as strigolactones (Marzec and Muszynska 2015), in growth of plant cells. It also seemingly delivers a narrative that provides a biophysical context for understanding the evolution of the apoplast, thereby "uncovering hidden treasures" in the as yet unscripted biophysical control systems in plants.
Summarising, I took a physical approach to plant growth to simplify the universal features that may explain the observations. I have fitted a simple algebraic function (a beta Euler function) to some growth rate data, thus capturing the variations with respect to temperature and pH, which show a surprisingly conjugate nature. Cross-disciplinary research at the interface between the physical and life sciences was accomplished. I first considered the Euler beta function enigma, which is the duality of acidic pH (or auxin-induced acidification) and temperature, in order to obtain the EWS for a living system. I started from the striking similarity (mirror symmetry) between the elongation growth of the grass shoots that were incubated in different pH environments or auxin-induced elongation in water (or a neutral buffer) with the respective growth at different temperatures. These observations were based on the hypothesis that the effect of temperature on elongation growth is effectively equivalent for the relative growth increments such as those that are caused by a change of pH or endogenous auxininduced acidification in the incubation medium. In order to resolve this ambiguity, I first used the beta function for both dependences, which were normalised prior to the comparison, in order to obtain the beta function shape parameters α and β. It turned out that even without referring to the biochemical underpinning, one might have concluded that the acidic conditions of the incubation medium or auxininduced acidification and environmental temperature can act interchangeably, at least when they are considered at a phenomenological (mechanistic) level. The numerically verified high accuracy of these complementary representations allowed us, among others, to treat this dual approach as a new apparatus for predicting the outcomes in different growth conditions, which is especially useful in changing climate environments. I presumed that by applying a beta distribution and continuously changing its character with the values of the shape parameters, our findings might be related to the evolutionary context that is connected with the migration of plants away from the Equator as the climate changed and their adaptation to the spatial distribution of pH in the soil as a substitute for high temperatures. The EWS may be also helpful in delivering analytical solutions (the respective (α, β, Ω S )-triads) for the assisted migration of plant species (Vitt et al. 2010) that are at the risk of extinction in the face of rapid climate change. In this application, the EWS tool can constitute a complementary method for genetic modifications in the assisted migration processes by permitting the best species to be selected for such an introduction from those having a similar (to within experimental error) triad.

Application of Critical Exponents
In physics, it is believed (Fisher 1966;Griffiths 1970) that the critical exponents are "universal", i.e. independent of the details of the Hamiltonian (energy function) that describes a system. The implications of these, however, are far reaching. One could take a realistic and complicated Hamiltonian, 'strip' it to a highly idealised Hamiltonian and still obtain precisely the same critical exponents. For instance, on these grounds, it is believed that carbon dioxide, xenon and the three-dimensional Ising model should all have the same critical exponents (Baxter 1989) and this appears to be the case (Hocken and Moldover 1976). In this context, since growth phenomenology (non-equilibrium orderings) and the STs approach have converged to a form an EWS for plants, the potential role of the critical exponents in discriminating the different modes of the wall extension (Breidwood et al. 2014 for review) of growing cells or tissues seems promising in opening new avenues of research.

Chemical Potential-Induced Ratchet-Like Growth in Plants
It is known that the group of oscillating variables during growth includes at least ion fluxes (Sukhov et al. 2011) and internal free concentrations of Ca 2+ , Cl − , H + and K + ions (Holdaway-Clarke et al. 1997;Holdaway-Clarke and Hepler 2003), the cytoskeleton (Geitmann 2010), membrane trafficking (Feijó 1999;Feijo et al. 2001;Zonia and Munnik 2008) and cell wall synthesis (see Portes et al. 2015 andHemelryck et al. 2017 for review). However, despite the progress made in this domain, the central core-controlling mechanism is missing . Moreover, how it produces a macroscopic outcome in the form of apical growth also remains uncertain. A possible link between this oscillatory component at the molecular level and the structurally (spatial arrangement) and temporally organised apical (periodic) growth is described below.
Indirect measurements of the chemical potential (μ H+ ) using pH can be compared to our previous measurements of the electromotive force (EMF) in order to detect phase transitions, which are localised by the 'kinks' in the chemical potential in condensed matter physics . Such kinks in the chemical potential always localised the occurrences of different types of phase transitions (where the change of state of a macroscopic system occurs) such as ferromagnetic, anti-ferromagnetic, superconducting, re-entrant or structural phase transitions . Based on these observations and considerations on the EWS for plants, one can combine these outcomes together in a single coherent account.
Launching reasoning from the EWS and putting the differing approaches that were outlined in the Introduction together, a plant cell expansion recurring model can be assembled. The latter, though unsophisticated for obvious reasons, not only describes one cell growth cycle in physical (mechanistic) terms, but also demonstrates how subsequent growth cycles can arise and even how the entire process may eventually decay (compare with Rounds et al. 2010).
The hypothetical steps in the transport of cations (such as H + ) against its chemical gradient by an electrogenic pump and its consequences are presented in Fig. 10 and SI Figs. 4, 5, 6. These can be described as follows: A. The initial state of a growing cell (I assume the chemical potential μ = μ 0 for reference only). B. The outward, active proton transport across the plasma membrane interface creates gradients of pH or equivalently gradients of the chemical potential μ = μ H+ , which drives the transport of ions and uncharged solutes. The plasma membrane (marked by a vertical line) depicts the boundary between the cell wall and the cytoplasm. C. H + -ATPase activity regulates cytoplasmatic pH (and controls cell turgor pressure). At the microscopic level, we can observe a singular behaviour at a certain critical value of the chemical potential μ = μ c (see also Fig. 5c). This chemical potential activity can be interpreted as a signature for an ST that is occurring in the cell wallduring an ST of a given medium, certain properties of the medium change, often discontinuously, as a result of a change in some condition, here: μ. D. The increased involvement of H + ions in cell wall building and exocytosis (Holdaway-Clarke and Hepler 2003) by δV raises the turgor pressure inside the cell by δP, thus producing force F, which may cause cell wall loosening and elongation by δl (see Wei and Lintilhac 2007 for a loss of stability approach and Sridhar et al. 2018 for a statistical model of expansive growth). The physical quantity that may be responsible for signal transduction here is pressure (e.g. Zonia 2010), which is a scalar quantity. E. The increased pressure at the apical plasma membrane by δP induces the stress relaxation of the cell wall, Fig. 10 Plant cell expansion recurring model-a proposal. a Initial state: the chemical potential μ (μ 0 -reference potential) equilibrated in both compartments; turgor pressure equals P; cell volume equals V. b Energy-consuming proton pump driven by the hydrolysis of ATP transports protons (ions) against the electrical gradient. c Chemical potential difference is established between the compartments-the system is far from equilibrium; turgor pressure equals P; cell volume V. d Synthesis of cellulose at the plasma membrane and of the pectin and hemicelluloses components with a Golgi apparatus depositing a layer on the inside of the existing cell wall, thus causing an increase of the cytoplasm pressure by the infinitesimal value of δP: P + δP(t), though the cell volume still remains V. d New wall synthesis. Elevated turgor pressure P by a δP(t) amount (Benkert et al. 1997) produces force F, which acts on the wall (energy stored equals PδV). The chemical potential starts to drop due to the overall particle loss in the system. e Increasing turgor pressure releases energy to extend the existing cell wall chamber by δl; the chemical potential further decreases before the next cycle (arrow). After the completed growth cycle, turgor pressure diminishes by the δP(t) value due to leakage or the other mechanisms illustrated in the text. Another growth cycle begins (arrow), however, with a new initial value of P equal to P -δP(t), V equal to V + δV and μ equal to μ -δμ (exaggerated). Legend Left compartment: plasma membrane (PM) and cell wall. Right compartment: cytoplasm. Grey-PM, blue-H + -"Fermi sea", greennew cell wall (Color figure online) relaxation-driven cell expansion by δV and triggers the start of the next growth cycle with a slightly modified turgor pressure ~ P − δP (see the linearly descending pressure plotted in Fig. 1 in Benkert et al. 1997). The system restores the initial configuration and allows a new growth cycle to begin (explicit in pollen tube oscillations). The latter means it is reset at the end of the cycle (Hemelryck Van et al. 2018).
As a result, the problem of causality, which was mentioned in the Introduction, is proposed to be resolved in a self-consistent manner, thus potentially leading to a paradigm change that can be expressed by a self-consistently evolving recurring model. Obviously, the latter should be further examined by numerical calculations (self-consistent solutions at subsequent time intervals) and confronted with the existing representations (Baskin 2015, see Fig. 1.1 therein for a complementary Maxwell-type model). The possible mechanisms of the participation of pH changes on photosynthesis (Grams et al. 2009;Sukhov et al. 2014;Sherstneva et al. 2015;Sukhov et al. 2016) can also be included in the form of a phase-locked loop, thus closing the 'system of equations'. The existence of oscillating behaviours in the investigated system, however, implies the existence of a control by feed forward or feedback loops, which are inserted into a self-regulating mechanism (Portes et al. 2015). The latter statement can be translated in mathematical terms into a self-consistent time evolution equation that corresponds with growth.

The Recurrent Model
Let us briefly discuss the major novelty of the recurrent model as compared to the earlier work in the context of the previous models. The main approaches that were outlined in the Introduction emphasised a single attribute that is responsible for expansion (or periodic) growth. However, it seems that none of these factors alone (such as oscillatory wall synthesis, excessive turgor pressure or turgor pressure oscillations) can explain all of the observable modes of cell expansion and growth. Moreover, because of the intricate dependence between the processes that are involved in tip growth, causal relations cannot be easily inferred solely from phase relations (Pietruszka and Haduch-Sendecka 2015b) as simultaneous information about the role of the participating ions and the actual growth rate amplitude are required.

Pollen Tube Growth and Oscillations
Furthermore, intracellular ion dynamics has been shown to play an essential role in pollen tube growth (Portes et al. 2015), while calcium, protons and chloride are also considered to be of major importance. Because the [Ca 2+ ] i pool at the apex oscillates, it has also been considered as supporting growth (Messerli and Robinson 1997). Additionally, it has been experimentally verified that pollen tubes stop growing when the [Ca 2+ ] i gradient is dissipated (Pierson et al. 1994). An interesting fact regarding chloride dynamics is that in contrast to Ca 2+ and H + , the ionic fluxes occur in the reverse direction at the apex, thereby elevating the chemical potential gradients even further.
Protons have received much attention as they play a crucial role in pollen tube growth. Since water, which ionises spontaneously, is the major component of living cells, H + ions are known to be involved in enzymatic activity, endo/ exocytosis and many other processes (Portes et al. 2015 for review). Therefore, pH inside cells has to be regulated in harmony with the growth phases (Messerli and Robinson 1998). It was found that pollen tubes showed an acidic domain that was located close to the membrane apex-the latter was present only during growth. In addition, the localisation of pollen-specific H + -ATPase provides evidence that indicates correlations with the proton extracellular fluxes and the acidic domain at the apex (Certal et al. 2008). It looks as though proton dynamics is involved in the mechanisms for maintaining the polarity that ensures pollen tube growth (Portes et al. 2015). Moreover, it has been established (Agudelo et al. 2016) that the pollen tube response is correlated with the conductivity of the growth medium under different AC frequencies, which is consistent with the notion that the effect of the field on pollen tube growth may be mediated via its effect on the motion of ions. This is, however, in accord with the presented chemical potential ratchet scenario and the representative time series of the oscillatory traces of incoming cations (H + and Ca 2+ ) or outgoing anions (Cl − ) at the growing tip. Hence, from the point of view of the recurrent model, the equilibration of the chemical potential for all of these species should also be the major mechanism for the cessation of growth. The ratchet scenario seems also compatible with tethers attachment/detachment/ reattachment dynamics in the molecular model proposed by Sridhar et al. (2018).
Ion homeostasis and signalling are crucial in regulating pollen tube growth and morphogenesis and affect upstream membrane transporters and downstream targets (Michard et al. 2016). Pollen tube growth is strictly dependent on ion dynamics. Ion fluxes and cytosolic gradients of concentration have been mechanistically associated with the action of specific transporters, especially protons. In the recurrent model I put forward, though extremely simplistic, H + chemical potential dynamics may serve as a "still missing" (Portes et al. 2015) central core-controlling mechanism or "pacemaker" that is able to produce a macroscopic outcome, i.e. structurally and temporally organised apical growth, especially since H + usually oscillates with the period similar to growth rate-see Fig. 2d in Damineli et al. (2017).
On the other hand, the short period oscillations of pH, which induce changes in wall properties, can be translated into the volume pulsatile extension by Eq. (20). Provided that cross-correlation (vs. time) of pH and P exists, turgor pressure should vary slightly ( Fig. 1 in Benkert et al. 1997) in accord with minute pH drops (note the pressure linear decrease); Eq. (20) allows for such analysis. This should also explain the observed phase delays in growth rate and ionic fluxes in pollen tubes in mechanistic terms (e.g. Messerli and Robinson 1998). This is, however, beyond the scope of this work and should be analysed in another study.

Final Comments
The aim of this work was to develop a unifying model to simulate plant cell expansion and growth. In the chemical potential representation, growth can be treated as a dynamic series of STs that can take place in the primary plant cell walls during subsequent time intervals. I suppose that such changes would proceed, e.g. from an ordered to disordered state, from a stressed to relaxed state, from a "covalent bonds state" to "disrupted covalent bonds state" or even through more subtle pathways in the wall polymer network, which themselves constitute an intriguing research task. Furthermore, the singular causality, whereby only biochemical wall loosening can lead to stress relaxation and turgor loss, is replaced by a self-consistent recurring model, in which all of the prevailing variables constitute the macroscopic evolution of the plant cell. A common denominator appears to emerge that ties all of the growth factors together, namely the chemical potential of intervening ions. [Oscillatory ion movement in the electrolyte solution ) accompanies the uni-directional increase in pollen tube, and changes the chemical potential (Fermi redox energy) of particular ions]. Further investigations can go along the lines clearly illustrated in Smithers et al. (2019).
At the beginning I asked, whether the volumetric growth of plants (in the dark) could be calculated from pH, T and P. It turned out that the answer is "yes", provided that it is delivered a posteriori, i.e. the pH growth data are already known. This reasoning, however, can be reversed-by knowing pH(t), T and P, one can precisely predict the cell volume V expansion over time using model Eq. (21).
As we found out, from the physical point of view, cell volume expansion in plants follows the pH(t), P, T , ⟶ V(t) analytic transformation through the coupling between the mechanical (pressure), thermal (temperature) and chemical (pH or μ) degrees of freedom. In this process, the chemical potential of the intervening ions can play an important role-distant (non-local) correlations are no longer controlled externally, but by the system itself, which is similar to the self-organised critical states Bak et al. 1988;Hesse and Gross 2014;Pietruszka and Olszewska 2019) and remains in agreement with the statistical approach (viscoplastic and viscoelastic terms mapped onto a statistics of independently acting tethers, Sridhar et al. 2018) or recent opinion delivered by Mangano et al. 2018 concerning pH oscillations in pollens. Indeed, it has been suggested that a growing cell proceeds along a critical trajectory over 1/f ridge (here f-fluctuation frequency of ionic channels), which is typical for regular growth (Pietruszka and Olszewska 2019;also Haduch-Sendecka et al. 2014 for the minimal energy loss in the natural growth conditions-a biological illustration of the principle of least action of Feynman and Hibbs 1965). The latter observation can be confirmed by the cross-correlations that are maximum for the control (normal) growth, see Figs. 4, 5 in Olszewska et al., 2018. Because of the pH(μ) dependence, the signalling cascade during the growth of a cell seems to be triggered and controlled by the chemical potential of vital ions.

Conclusions
By resolving the duality of low pH or auxin action (producing acidic pH) against temperature, not only it was introduced an EWS for the realm of plants, but also, by considering wall extension growth as a dynamic cascade of chemical potential driven STs, the critical indices (exponents) for this phenomenon were identified, which exhibit a singular behaviour at a critical temperature and critical pH in the form of power laws. This is concluded that the cell growth evolves along a path (trajectory) of the least action ('best-path'), thereby optimising growth under any physiological conditions. The inherent pH dynamics in close-to-natural conditions appears to essentially be responsible for this extreme trajectory, thereby providing a highly non-linear transformation pH(t), P, T , ⟶ V(t) , where V and P stands for the volume and internal pressure of a cell, respectively. A common unifying principle that can be applied for either condensed matter (studied in physics) or living matter (a subject of biology) is proposed for all of the chemical potential-driven phase or state transitions, respectively.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creat iveco mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.