Budding E^{*}, Gündüz G, Ozel ME and Demircan O  
Astrobiology Group, Çanakkale Onsekiz Mart University Çanakkale, TR 17020, Turkey  
Corresponding Author :  Budding E Astrobiology Group Çanakkale Onsekiz Mart University Çanakkale TR 17020, Turkey Tel: 0064 4 232 6388 Email: [email protected] 

Received January 30, 2015; Accepted February 20, 2015; Published February 28, 2015  
Citation: Budding E, Gündüz G, Ozel ME, Demircan O (2015) Biospheric Energization. Astrobiol Outreach 3:128. doi:10.4172/23322519.1000128  
Copyright: © 2015 Budding E, et al. This is an openaccess article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.  
Related article at Pubmed Scholar Google 
Visit for more related articles at Journal of Astrobiology & Outreach
We discuss a physical model relating to the origin and development of lifelike processes, paying particular attention to the continual growth of systemic energization and complexity from very simple beginnings. To approach recognizable life forms in the terrestrial case, we refer to a molecular (‘ABC’) process consisting of the following elements: A  a molecular structure that can temporarily stores a local energy excess; B  a catalyst that, making use of the local energy excess, catalyses also itself, and comes in a small number of slightly different forms involving different reaction steps; C  a disposable tool molecule that assists and tunes the response of the A and B molecules to the (timedependent) local source of energy. This local source of energy is external to the ABC system, and we take it to be solar photons (or a comparable flux from a solarlike star).
We show that the system can evolve from very simple beginnings to a progressively more highly tuned, energized and complex responding biosphere, that grows exponentially; albeit with a very low net growth factor. We note the inherent instability suggested by this can be related to the “low L” solution of the Drake equation. We consider possibilities for similar processes in the outer regions of other planets of the solar system or elsewhere. We outline some possible observational checks for the principles of this model.
Keywords  
Biosphere; Energization; Star; Boltzmann; Solar flux  
Definition of Life  
NASA has proposed that life be taken as a selfsustained chemical system capable of undergoing Darwinian evolution. This definition can be found consistent with general physical principles. The life chemistry encountered on Earth should accord with these principles, but as one particular case. This viewpoint allows us to seek for life in more generalized environments, although in this article we confine attention to planets and possible living systems having understandable comparisons to the terrestrial example. Further discussion of this and related matters was given by Budding et al. [1].  
Basic Physics and Life  
Budding et al. [1] contrasted the energization of locomotion in classical mechanics with that of a living system, noting the distinguishing role of metabolism, which implies a continuous supply of power directed specifically to maintain living [2]. Metabolism’s highly ordered energization system essentially involves that field of chemistry called biochemistry, and it is within the spirit of Darwinism to regard the latter as having evolved from the former. We may thus regard life as a selfsustained and evolving chemical system, but reducing this selfsustenance and evolutionary character to physical terms calls for close consideration.  
The least action principle of classical mechanics has parallels with the probability maximization underlying Boltzmann’s distribution of the energies of an assembly of particles in statistical mechanics in that they both involve mathematical extreme. Boltzmann’s distribution refers to a temperature T_{E}, which, in turn, implies a surrounding environment and associated radiation field. The argument is direct so long as T_{E} can be specified, even if the assembly requires some response (or ‘relaxation’) time τ_{R} for the distribution to be attained. It becomes less direct if the environment’s temperature is allowed to vary, on a timescale τ_{E}, say. If τ_{E}<<τ_{R} Ethen we can expect a Boltzmann distribution for any given TE. We could also expect a comparable distribution when τ_{R}>>τ_{E}, if there is some average <T_{E}> that applies to the long timescale τ_{R}. It becomes more complex if, for some system, is of the same order as τ_{E}. A quasiresonant configuration can be anticipated by moving to the frequency domain ν (=1 /τ ), with . In particular, if the possibility of driven or parametric resonance occurs: the responsive system becoming progressively more energized from the cyclic driving variation affecting the environment. This idea strongly influences our approach. We visualize a steady progression from a simple, low Qfactor (energy stored/energy lost per cycle), fairly uniformly dispersed, and marginally ‘disequilibrium’ chemical activity to an ever more stringently defined, highly differentiated, complex and hierarchically structured biochemical matrix.  
This possibility is not overt in Boltzmann’s [3] view of a predominating universal tendency towards random disorder (‘heat’), in line with the second law of thermodynamics. While sensing the disparity between this tendency and the existence of living systems that grow, evolve and become more complex, Boltzmann and others have assumed that the solar flux can drive the biosphere, rather in the way that the potential energy gradient from an inland water head drives a (not necessarily simple or static) river system. Boltzmann then regarded Darwinian evolution as a competition for entropy within living systems. Entropy, in the sense of number of available states, or information required to cover the essential properties of a system, may well increase in evolving lifeforms. However, such forms are not consistent with the specific probabilities for the arrangements that would be found in thermal equilibrium (e.g. a Boltzmannian distribution of particle energies). A satellite picture of the Earth at night, for instance, readily reveals the considerable excess of emitted highenergy photons above what would correspond to the nighttime temperature of the planet.  
Schrödinger [4] probed further into this issue, paying particular attention to the problem of “order from disorder”. Schrödinger interpreted this as a local absence of equilibrium (i.e. disequilibrium) of the living subsystem: the measure of disorder could be reduced in this subsystem at the expense of an increase in that of the surroundings overall. Within an enclosed system the continuation of such disequilibrium contravenes thermodynamical understanding. However, if the containing environment is unbounded (‘open’) the overall entropy may still increase, while a subsystem becomes more ordered. While the maintenance of thermal disequilibrium in the biosphere thus necessitates an open system, lack of closure of the overall system has still not been shown sufficient to generate and progress living forms to the degree observed and in the time available.  
More recently, a distinction has been drawn between increasing disorder and entropy. When a disequilibrium substructure occurs, ordered processes can be set up that work against the energy gradients producing the disequilibrium more effectively than in a randomly disorganized configuration. Such processes have been termed “dissipative structures” [5], and their action consistent with a maximal “energy degradation” [6]. This approach anticipates lifelike processes developing in any driven disequilibrium in a general way, while not yet specifying the driving mechanisms or any cycling (periodic) interactions that may be found in particular identifications of life.  
Living systems, while generally composed of a large assembly of molecular particles^{1}, are usually recognized, biologically, as individuals, having a high degree of ordered selfcontainment. Biologists often confine attention, regarding what is agreed to be fully alive, to such ordered structures, either monocellular or groups of cells with sufficiently strong interactions to constitute a multicellular organism. But, although life on Earth has evolved to this highly differentiated condition, in which cellbased life commands center stage, there are reasons to take an integrated or ‘biospheric’ view.  
A natural place to look for the origin of the power driving the terrestrial biosphere is the Sun. The solar flux amounts to around 1.36 kW for each projected square metre above the Earth’s atmosphere. The total incident power over the whole Earth is then around 2 × 10^{14} kW or ~7 × 10^{21} J per day. There may be other energy sources, but the location of living matter about the Earth and the central role of photosynthesis within biology [7] push solar flux as the most likely main driver. A reasonable estimate for the total energy driving the entire biosphere is ~10^{12} kW [8]. The solar flux thus has enough power to drive this ‘engine’: the problem is about how it works. Life is not simply a dissipative mechanical system like the weather: it evolves to greater complexity while exhibiting Darwinian selfsustenance, whose nature invites clarification.  
Chemical Potential Energy of the Biosphere  
We can consider the interplay of potential and kinetic energies also at the microscopic level of the arrangements of charged particles within atoms or molecules. ‘Chemical potential energy’ (CPE) then refers to the electrostatic potential associated with the relative excitation levels of valence electrons. Chemical interactions are mostly accounted for by the changing configurations, or energy states, of such electrons. If, as a result of some anisotropy, exothermic reactions occur, particulate potential energy tends to diffuse irreversibly into electromagnetic radiation. An endothermic process implies an atypical initial configuration within a given environmental temperature, which becomes regularized into a more probable arrangement. In either case, physical results are comparable to the degradation of macroscopic motions [1]. In normal conditions, naturally occurring processes tend toward configurations whose typical valence electron energies are comparable to those obtaining randomly at the ambient temperature of the local environment, i.e. a Boltzmann distribution.  
But this is not the case in the surface regions of planet Earth! There is a large residual CPE among the reagents of the biosphere, whose existence can be seen as the key physical problem posed by life. This energy can be estimated to be of the order of 10^{23} J in total. This begins to compare with the total thermal energy of the Earth’s atmosphere of ~10^{24} J; that is ~1.5 x 10^{5} J kg^{1}, while the average energy density of biotic matter is now some 3 x 10^{7} J kg^{1} [9,10]. The physical significance of the ratio of these two energy densities (~200) was pointed out by [1].  
A similar point arises from the existence of 3.71 x 10^{19} moles of free oxygen in the Earth’s atmosphere, that can be compared with a more feasible initial planetary atmosphere made up of lower chemical potential (more tightly bound) molecular arrangements, such as carbon dioxide, water, perhaps methane, or other constituents. Methaneoxygen combustion yields some 572.2 J per mole [11], so that a reverse process of oxygen liberation can be interpreted as a net potential gain of ~2 x 10^{23} J [1]. The very significant approximate equality of this datum to the potential energy yield of all organic matter around the Earth’s surface was pointed out already a long time ago by Haldane [12], and it was reiterated in ref [13] account of the biosphere’s energization.  
Typical chemical reactions that might occur naturally near the temperature and pressures about the Earth’s surface involve such short interaction times that local excess molecular energies should be swiftly radiated away or lost through kinetic collisions with other molecules, until approximate equilibrium with the thermal energy of the surroundings is achieved [14]. Life on earth is thus generally recognized to be distant from the thermal equilibrium of standard thermodynamics [15]. But we can reasonably assume that the excess CPE of this disequilibrium came from the solar influx converted and retained by living forms: if not, then from where?  
Homeostasis, Restoration and Darwinian Selfsustenance  
A characteristic property of a living system is its homeostasis: a restorative behaviour ‘attracted to’ continued living [16,17]. This attraction is strongly related to its Darwinian character or ‘aim to survive’. Most explanations of biological behaviour refer either directly or indirectly to this property; through its role in cell metabolism or implications for evolution. We focus now on the essential nature of the restorative mechanism underlying homeostasis.  
Budding et al. [1] background this point in more detail. They consider, as a general case of restoration, the SturmLiouville type secondorder differential equation  
(1)  
(time derivatives denoted here by D{ }); known as the parametric oscillator equation. The left members of this equation include first the accelerative term in D^{2}, followed by the dissipative term λ (t)D. The oscillatory coefficient , with mean value ω_{0}, includes also a driving component with mean frequency , that can cause a steadily increasing energy of motion for suitable values of period and amplitude.  
The responsive oscillation y, satisfying Eqn (1), has an amplitude y_{0} that depends exponentially on the time t, i.e.  
(2)  
That grows or declines according as the stimulation coefficient α exceeds that of the dissipation μ. The exponential character of this energization should be noted, even if the growth rate may be very small. The stimulation α depends on the basic driving frequency. The latter can be taken to be ~1 day^{−1} for the insolation, but in our scheme there are actually two processes separated by ~ 0.5d, so parametric driving could occur if the system’s natural mean response time was about a day, and having an initial energy amplitude of the same order as the external driving function f_{0}. The dissipation μ depends on the integrated effect of the damping coefficient (λ in Equation 1), which in turn relates to the available behaviour of the system, i.e. at what rate it loses potential energy.  
If we suppose that the CPE of a chemically responsive system, that can be regarded as a parametric oscillator subject to a daily solar energy input f(t), scales with the amplitude y_{0}(t) as above, then the current average CPE excess of biotic matter (~200, see Section 3) implies that 1012 (α −μ ) ≈ ln 200 , or the mean value of (α −μ ) ≈ 5 x 10^{12} day^{−1}.  
Consider now the growth of oxygen in the Earth’s atmosphere [18,19] from the Earth’s formation through to the present. The average slope for the logarithmic growth can be worked out as ~ (4.30/4.8By), which when translated into the exponential growth rate per day amounts to ~ 5 x 10^{12}: essentially the same as the foregoing estimate. The same agreement is reflected in the comparability of the CPE of the biosphere regarded (a) as its presently combustible total energy, and (b) the energy required to liberate the mass of free oxygen in the Earth’s present atmosphere from likely initial constituents noted in the preceding section. The liberation of bonded oxygen, with its implication of CPE enhancement, gives direct support for our resonating oscillator model of biosphere formation and development.  
The selective quality built into the solution (2) can also be seen. Different chemical pathways can be associated with different values of the net energization coefficient (α −μ ) . If such differences correspond to alternative possibilities for the same components in the same location and subject to the same dissipative effects, then the pathway that locally maximizes (α −μ ) will predominate there, due to its energy advantage. This interpretation physically parallels the Darwinian view of an advantaged mode having an enhanced chance of survival (‘natural selection’)^{2}.  
Molecular Mechanism  
Budding et al. [1] sketch a molecular model that can behave as a parametric oscillator according to the foregoing, with references to relevant earlylife studies that may substantiate their scheme experimentally. Here, we summarize the main points. The biochemical model for a circadian bacterial oscillator of Pavlidis and Kauzmann ([20] hereafter PK) can be compared with this scheme, although ours is directed in the first instance to a prebiotic ensemble of relatively simple molecules, from which a continually developing and recognizably Earthlike, biosphere evolves. Other evidence on the fundamental role of rhythmicity in biological systems was cited by Budding et al. [1].  
Numerous alternative pathways and branchings may occur in relevant biochemical reactions. This issue was discussed by King [21]. with a schematization similar to the one of Budding et al. [1]. King gave a necessary condition on the ratio of sidereactions to promoting reactions that should be satisfied in the milieu in which an autocatalytic sequence could be maintained. While such schemes tend to be directed initially towards originating a sugarbase combination comparable to a simple nucleoside, they could apply at different levels. With this in mind, the ‘ABC’ formulation which follows should be regarded as having a wider application than might hold for specific molecular combinations. The underlying issue is one of suitably timed molecular bindingenergy decreases (potential energy increases) that make available work used to support or build structures that enhance stability of the autocatalysis, with feedback in the interactions of mixed chemical systems a key element. The first stage in Budding et al. scheme is an endothermic reaction sequence, induced by the solar energy flux f, in the presence of a catalyst (B). In the terrestrial case, it should have originated from some very simple protocarbohydrate, such as formaldehyde (CH_{2}O), that would engender more complex subsequent reactions [2224]. The main outcome is a molecule deriving from an oxygen reducing process with enhanced CPE that is stable over the timescale of at least a halfday. This end product ‘A’ is then a photosynthetically produced sugar: probably, but not necessarily, along the lines considered by Garrison et al. [25,26]. This ‘Asubprocess’ we write as:  
(2)  
Where n is an appropriate small integer. The (organic) molecule A can be regarded as a temporary energy source and driver of subsequent molecular interactions. Pertinent to the origin of modern life forms was the discovery by Powner et al. [27] of a pathway to the nucleoside ribocytosine involving, as a key reactant, 2 aminooxazole. Figure 1 of Powner et al. [27] shows this latter molecule in an A type subprocess, in which glyceraldehyde (C_{3}H_{6}O_{3}) is catalysed to arabinose (C_{5}H_{10}O_{5}). Nucleosidic sugar + base assembly then follow, even if not by the more direct pathway supposed by early experimenters. We indicate in (I) the key catalytic role of the second conceptual basic molecule by B: the brackets around it are meant to suggest that other molecules may affect its action. Continuity of the AB cycle is achieved during the hours of darkness, when the A type molecule is involved in a mildly exothermic (energy release φ is less than f absorbed in the A process) reproduction of the B type molecule in the `Bsubprocess’.  
Evolved B type molecules on the Earth could be expected to include ring structures (as A), and, keeping in mind the amino basis of enzymes, involve nitrogen [2830]. PK denote this catalyst by E, with an excited state as E+, identified as adenosine diphosphate (ADP). The PK modelling underlines the intimate relationship of the two components A and B, condensing to the compound AB at the center of the sequence and controlling its outcome.  
Using the formalism of the previous section, we schematize the A subprocess in the form  
(3)  
where D{A} denotes the rate of sugar molecule formation, in dependence on the local density of B type catalysts, but with dependence on the variable solar insolation, mainly through the environmental function (e.g. local temperature) embodied in the term f_{1}(t)^{3}. In practice, the coefficient k1 would also be, to some extent, timedependent. Within the present highly simplified context, we expect other environmental factors to remain essentially constant.  
Powner et al (2009) note the primary role of urea as a source of the B type component in their nucleoside generation model. With this in mind, we can write for a Bsubprocess:  
(II)  
In this very simple approach, B is a derivative of the sugaramine condensation in an ammonia solution. In a similar way to Eqn (3) we can write for this B subprocess  
Ω (4)  
Where A indicates the amount of intermediate AB compound (=EX in PK), for which we write  
(4a)  
If we differentiate the first relation again and substitute from the one just given we find  
(5)  
With ω^{2} =k_{1}(k_{3}k_{4}  k_{2}), λ =k_{3}k_{4} and f_{2}(t) is the oscillatory function f_{1}’ + k_{3}k_{4}f_{1}, together with other small contributions arising from time derivatives of the coefficients k i. Equation (5) has the parametric oscillator form, yielding a consequent exponential growth in the appearance of sugar type molecules, if α > μ.  
In a similar way we can show that B satisfies the similar form  
(6)  
With f_{3}=(k_{2} −k_{3}k_{4}) f_{1}. Since we expect oxygen to be liberated as a byproduct of the AB process it too shows a similar steady exponential growth, and, by implication, the chemical potential energy of the process as a whole.  
A key point is that the natural timescales (02π /ω_{0}) for both A and B subprocesses should be of about a day in typical terrestrial conditions. Such timescales, demonstrated by Weber [30] for feasible real molecules, are relevant to the susceptibility to diurnal resonance for the whole scheme. The intermediate AB compound has an important regulatory role. The AB process overall is controlled, multiplicative and CPEenhancing, thus showing elementary lifelike properties.  
The simple reactions sketched as (I) and (II) above point to biochemistry evolving from a simpler chemistry in a terrestrial environment, with formative roles for the quasiorganic (i.e. naturally occurring) molecules formaldehyde and urea.  
Budding et al. [1]. took the original Asubprocess to be natural photosynthesis, following the idea of Garrison et al. [25]. Adopting some orderofmagnitude parameter values, they derived a natural oceanic daily oxygen yield of ~10^{9} gm, which would act as a seed for the subsequent gradual atmospheric buildup. The present atmospheric proportion of ~21% points to an advanced stage of the terrestrial biosphere.  
LotkaVolterra Formulation  
The oscillatory behaviour coming from Eqns (5) and (6) can be reconciled with a LotkaVolterra (LV)like formulation of the system, albeit with some small modification. Let us write  
(7)  
(8)  
Equation (7) accounts for a numerical proportion α of free A type molecules in a local environment, as in the Asubprocess, regarding this as dependent on the proportion of catalytic B type molecules β, as well as some (oscillatory) external power source (driving function) f_{1}, while the AB compound has sequestered free A types that would otherwise be available for such condensation, to an extent that depends on the product αβ. In a similar way, Equation (8), parallels the B subprocess, where B type molecules are liberated from their ABtype source, but may then be removed from the milieu by association with ambient free A types. Equations (7) and (8) are complementary to the standard LV form, normally written (without the driving term) as  
(7a)  
(8a)  
But it can be shown that the “collaborative” forms (7) and (8) share similar oscillatory behaviour to the “competitive” relations (7a) and (8a).  
In Figure 1 we show the results of a numerical integration of the system (7)(8), against linearly scaled (but by arbitrary amounts) axes. The lowest curve follows the changes of α (A type molecules), and the middle curve the complementary changes of β (B types). In terms of the original [1], scenario, the peaks of Atype production could be expected around noon, or in the early afternoon. B type maximization has occurred by the preceding early morning hours. The highest curve follows the more constrained variation of AB type molecule development, having twice the frequency of its progenitor molecules.  
Predatorprey interactions have been presented more in terms of stochastic encounters by Maynard Smith [31], and the corresponding AB synthesis problem can be so formulated.  
Figure 2 shows an example of such interactions where the diurnal input power term appears as the large amplitude blue sine wave at the bottom. The upper red curve corresponds to Atype molecule concentration, with Btypes appearing as the green curve. The AB molecular density shows up in purple in this diagram. The longerterm exponent indicates the steady growth in the numbers of such molecules.  
Growth, Evolution and Complexity  
The AB cycle thus acts as a parametric oscillator in a nearresonant condition with the diurnal rotation and molecular condensation timescales. During the day, B catalyses the production of more A molecules from initial constituents. Towards night, Atype molecules condense with the B type and at lower nocturnal temperatures this leads, via the AB intermediary, to the formation of more B molecules [20]. Sometimes a somewhat different molecule B’ may be involved. The possibility of rival catalysts B and B’, or control compounds AB or AB’, institutes molecular evolution. This is on the basis of preferential selection in a given environment depending on efficiency of the resonance (B or B’) for CPE generation, even at this fundamental level [32]. This simple AB or AB’ dichotomy has an analogue in the ‘black’ or ‘white’ conceptualization of Daisyworld [33]. Early biospheres on planets with moderately uniform topology could thus be expected to show a zonal structure. Zonal propagation of the energization may also be approachable via the Telegrapher Equation in spherical polar coordinates, but we consider here other issues arising with the developing energization.  
If driving exceeds dissipation, the process must grow, but not simply in linear amplitude of excitation or numbers of A and B molecules. The available modes of carbonbased molecular chemistry allow for the formation of more complicated, higher CPE, structures. With ongoing driving, the catalytic function becomes not just from single molecules, but groups of more complex, finely tuned, molecular sequences, genes, or macrostructures, cooperating towards the resonancedriven CPE enhancement. There should then be not just one, but a spectrum of AB cycles, filling in the available bandwidth for which α >μ.  
But how does controlled organization develop? It is the AB compound that promotes survival of the surrounding protobiotic molecular ensemble through its negative feedback. We indicated above that the B subprocess is associated with the release of energy φ, but this is at a rate determined symbolically by the value of the quantity k4. At this point, we encounter a third molecule − the C type − that chemically relates to the AB compound. In effect, reactions involving, or ‘targetting’ such C types by the AB (or AB’) types influence the dissociation rate k4 in (4a). In turn, this affects the natural frequency ω and subsequent rate of multiplication and evolution, depending on the separation of from Ω or possibly submultiples of Ω. C type molecules, coding optimal responses for the relevant local environmental conditions, can be regarded as proteins in the terrestrial case. They are outside the catalytic ABAB’ protonucleotide sequence, but would be maintained through subsequent multiplication if promoting the underlying energizing resonance. C type molecular structures will naturally evolve an increasing role in regulating local conditions to promote internal homeostasis of surviving forms.  
In due course, and with the steady increase of biospheric CPE, AB type molecules themselves compound to more complex ABAB’ forms, made possible by the versatility of suitably energized carbon bonding. This development can be expected around nodes in the environment, where advantageous circumstances for the growth of one mode (AB, say) give way to another (AB’). The appearance of ABAB’ forms at a node can recapitulate the same modified LV equations (78) as for the original AB combination, with α now standing for AB and β for AB’. Note then the early and essential appearance of recapitulation as a characteristic of living processes.  
Maintenance of the driving resonance implies that the mean daily B molecule production must be maintained by the overall ensemble, implying, at a fundamental level, the continued demise of individual AB types, i.e. the death of particular or individual forms. Genetic ‘immortality’ [34] might perhaps be countenanced in that a given driving function in a given environment will correspondingly reproduce, by the most advantaged pathway, essentially similar advantaged molecules.  
ABAB’ (etc.) molecules retain more information about environmental variation, assisting multiplication and CPE increase, but in a more specifically defined local environment; the different links in an ABAB’ molecule having complementary regulating actions. The slow buildup of more complex ABAB’ structures and corresponding C or C’ molecules, also deposits a proportion of the net CPE increase into fewer molecules, giving a proportionally large increase in the energization per molecule of longer protonucleotide chains. Surviving development, involving these more complex sequences, manifests as more organized behaviour, with highly energized long molecules imparting hierarchical structure into the control function.  
Local conditions harmonizing with this increasing specification and control may tend to concentrate advantaged lifelike processes into lacunae containing local center’s of organization. Rock or tidal pools, as well as clay molecular binding templates, have been surmised to provide an early natural realization of this scenario [3537]. This relates to the feedback affecting the distribution of A, B, AB, ABAB’ (etc., as well as corresponding C types) within a lacuna. This control distribution can be expected to develop a pyramidical or ‘sandpile’ form [38], with characteristic slope γ coming from the stability requirement that δU / δE ≤ γ, where δU represents a mean fluctuation of local energy density (e.g. local temperature changes) and δE corresponds to the typical interlink binding energy of an ABAB′ type molecule. Numbers of corresponding AB, ABAB’, ABAB’AB” (etc.) molecules would then tend to decrease geometrically in proportion to γ, while the overall height of the pyramid (i.e. length of the longest molecule in the pile) increases in proportion to the general growth (αμ)t.  
On the basis of the observed number n of AB types, in a molecular chain containing N such pairs that may become indefinitely large, we find the number of permutations of the AB complexions . The complexity measure , regarded as a function of time, using observed fossil data, can be seen to show a logarithmic character in its general growth [39]. Using the Stirling formula for log N!, we can confirm, since the argument of the complexity measure is itself a logarithm, that the general growth of the number of AB type molecules must itself be exponential, in order that the outer logarithm’s argument be approximately linear with time. Budding et al. [1] argued, in this way, that the exponent of complexity growth derived from fossil evidence is generally consistent with the value inferred from the current energization of the biosphere derived in Section 4. More detailed features in the growth of complexity require further consideration of the particular circumstances of development.  
We argued above that the restorative development in parametrically driven molecular ensembles implies that neighbouring particle combinations (in space and phase) other than those initially propelled would tend to take up the same local energization, provided they have the same, or very similar, properties. An interesting question is then this: at what stage does this ‘bottom up’ restorative development of the biosphere start to look more like ‘top down’ regulated replication of its infrastructures? Growing complexity of the latter seems likely increasingly to mask the underlying biospheric driver (cyclic insolation) as evolution proceeds: however, the fundamental role of biospherically generated global reagents, such as O_{2} or CO_{2}, shows that lacunal life can never be entirely selforganized.  
Do the organized molecular lacunae become cells? Cells would be a type of lacuna, in the way we have presented this organized energization scenario; but it is more useful to keep the lacuna concept that can pertain to a variety of levels separate from that of a cell, which has a much more clearly defined identity and function within biology. The rate of development of biospheric complexity indicates that cells have probably grown and evolved from the more primitive, but relatively large, viral proportion of the terrestrial biomass. In effect, these are energized ABAB’ structures having particular C or C’ type molecules regulating their function sequence and natural timescale. The formation of collaborative substructures within a distribution of energized individuals to optimize net performance is a wellknown issue in organization theory. It is reasonable to seek the origin of cells from viruslike precursor combinations as an example of this.  
One single resonating AB process can be conjectured to have become the predominant driver of early life like development on Earth. This is suggested by the single A type molecule (ribose) and the prevailing chirality of evolved biotic molecules, as well as a preference for the isotope C_{12} in living structures [40]. This same basic AB process may still be present within the ongoing biosphere, but the original molecular behaviour would by now have been considerably developed by secular evolution in a parallel way to the changes of atmospheric composition.  
Recapitulation may help gain insights of this. The requirement for homeostasis within a surviving multicellular organism recapitulates that of individual cells. Stability of the biosphere recapitulates that within its lacunae. Surviving lacunae should thus interact collaboratively in general, though competitive episodes are clearly visible in nature. The exponential development of the biosphere is also echoed in Haeckel’s conjecture that ‘ontogeny recapitulates philogeny’. The rate of phylogeny maps the logarithmic measure of biospheric complexity, so that time intervals for given increases in log (complexity) of species within the biosphere are proportional to their corresponding dates of first appearance. The complexity of an individual soma of an advanced species, meanwhile, should proceed linearly with time along its genetic (molecular) map during its own growth. Haeckel’s conjecture implies that the early stages of this growth in complexity are telescoped into a relatively very short timescale.  
Observational Points  
A corollary to our growth model for biospheric energization is the implied longterm instability coming from the small, but positive, exponent in the solution of the underlying parametric oscillator equation. Although not identified as an inherent feature of the biosphere, the likelihood of its ultimate catastrophe has been referred to from various points of view [41]. The situation in which the lifetime L of an advanced technological civilization was essentially short (i.e. L<~100 y) was discussed by Bracewell [42]. In reviewing alternative scenarios of interstellar communications corresponding to different solutions of Drake’s (1961) equation [43]. This alternative is the one most easily reconciled with Fermi’s famous “where are they?” question. We do not detect “them”, nor they us, because, cosmically, they (and ourselves) are not around long enough for any chance of meaningful dialogue.  
Other observational implications of this energization model were presented by Budding et al. [1], who referred also to other planets in the solar system. New, highprecision exoplanet photometry may offer interesting opportunities in this respect. Close monitoring of the Jovian atmosphere shows that the Great Red Spot (GRS) must be dissimilar to general atmospheric cyclones that follow the formation and dissipation timescales of expectable meteorological effects (i.e. a few rotation periods). The GRS, on the other hand, actively propels itself against dissipative turbulence in its surrounding belt. The longterm survival of the GRS against this atmospheric turbulence can be estimated, from simple aerodynamic comparisons, to involve a power injection of order 10^{16} W. The mean insolation over the GRS could be reckoned as of this order or somewhat greater: so that the Sun is able to drive the process energetically. The question of how this power might be transferred into organized restoration of this largescale dissipative structure remains to be clarified, but the GRS is definitely a candidate for protobiotic energization. Excitingly, surface inhomogeneities comparable to the GRS but on the recently discovered ‘Jupiters’ of other stars may become photometrically detectable in the foreseeable future.  
The discussion presented in this article is indicative and exploratory only. But we believe it provides a logical framework for consideration of the energization and progressive evolution of biospheres in a cosmic context, as well as an objective basis for experimental testing. It links together various fields of science and implies a multidisciplinary approach. In this pursuit, we would encourage, and have hopes for, further broad international collaborations to help drive and coordinate functions.  
Acknowledgements  
The authors acknowledge support and helpful discussions from members of the COMUbased Astrobiology Group ([email protected]).  
^{1}The term molecule in this article denotes a very small particle, i.e. a somewhat looser usage than would normally apply to a quantified chemical reaction.  
^{2}A physicsbased approach to the originoflife problem published by England (2013) recapitulates the point made about growth in the foregoing model (cf. Budding et al., 2012). Referring to England's (2013) statistical physics of a selfreplicating population of individual life forms, the population is set as n at time t. The population being given, the probability p(n)t is here 1. At time t + Δt, when the population is n + Δn the probability that the number is n has now changed by –Δq from 1, since Δn is generally nonzero. The probability p must always decline from certainty, in this context, (Δp < 0) but the change of population can be either positive or negative (Δn > or < 0). England's (2013) formula (9) recognizes that this probability depends on the number n itself, and introduces two other terms g and δ, corresponding to given or adopted rates of growth and decay, respectively, to handle the variation of Δn. England's formula gives, in effect, (1a) where Δq_{b} and Δq_{f} correspond to backward and forward displacements from the situation at time t, to a time when the population number has changed by Δn. When n becomes large, the increments small and the behaviour of the population changes smooth, we can reasonably consider the changes −Δq_{b} and −Δq_{f} to be numerically the same as the change Δp over the increment Δt. We can then cancel out the equal factors, and write (2a) or, by integration, (3a) where n_{0} is the size of the population at an initial time. But (3a) is the same as Equation (2) above, supporting the relevance of the driven parametric oscillator model of biospheric energization to the growth of biospheric life forms, within the thermodynamic argument of England [22]. 

^{3}Note here that f_{1} contains a suitable coefficient ensuring that Eqn (3)'s units – i.e. those of the local concentration rate – balance.  
References  

Figure 1  Figure 2 