Dissipative Photochemical Abiogenesis of the Purines

We have proposed that the abiogenesis of life around the beginning of the Archean may have been an example of “spontaneous” microscopic dissipative structuring of UV-C pigments under the prevailing surface ultraviolet solar spectrum. The thermodynamic function of these Archean pigments (the “fundamental molecules of life”), as for the visible pigments of today, was to dissipate the incident solar light into heat. We have previously described the non-equilibrium thermodynamics and the photochemical mechanisms which may have been involved in the dissipative structuring of the purines adenine and hypoxanthine from the common precursor molecules of hydrogen cyanide and water under this UV light. In this article, we extend our analysis to include the production of the other two important purines, guanine and xanthine. The photochemical reactions are presumed to occur within a fatty acid vesicle floating on a hot (∼80 °C) neutral pH ocean surface exposed to the prevailing UV-C light. Reaction–diffusion equations are resolved under different environmental conditions. Significant amounts of adenine (∼10−5 M) and guanine (∼10−6 M) are obtained within 60 Archean days, starting from realistic concentrations of the precursors hydrogen cyanide and cyanogen (∼10−5 M).

From the perspective of the "Thermodynamic Dissipation Theory of the Origin of Life" [14][15][16][17][18][19][20][21][22][23][24][25][26][27][28], descriptions of molecular synthesis alone, however, are insufficient to capture life's vitality, which includes proliferation, dynamics, and evolution. The vitality of life requires instead a constant dissipation of an imposed thermodynamic potential. Of the above listed free energy sources employed in synthesis experiments, only UV light would have been sufficiently energy dense and continuously available throughout the Archean. Our theory therefore affirms that the "fundamental molecules of life" (those in the three domains of life) were, at their origin, microscopic photon dissipative structures (i.e., pigments) which arose "spontaneously" under the thermodynamic imperative to dissipate the prevailing UV-C (210-285 nm) solar photon spectrum ( Figure 1) [14,15,21,26]. Figure 1. The spectrum of UV-C light available at Earth's surface before the origin of life at approximately 3.9 Ga and until at least 2.9 Ga (curves black and red respectively) and perhaps throughout the entire Archean until 2.5 Ga. CO 2 and probably some H 2 S were responsible for absorption at wavelengths shorter than ∼205 nm, and atmospheric aldehydes (common photochemical products of CO 2 and water) absorbed between about 285 and 310 nm [29], approximately corresponding to the UV-B region. Around 2.2 Ga (green curve), UV-C light at Earth's surface was completely extinguished by oxygen and ozone resulting from organisms performing oxygenic photosynthesis. The yellow curve corresponds to the present surface spectrum. Energy fluxes are for the sun at the zenith. The names of the fundamental molecules of life are plotted at their wavelengths of maximum absorption; nucleic acids (black), amino acids (green), fatty acids (violet), sugars (brown), vitamins, co-enzymes and cofactors (blue), and pigments (red) (the font size of the letter roughly correspond to the relative size of their molar extinction coefficient). In this work, 125 W m −2 µm −1 is the "nominal" incident intensity at 260 nm, but we also consider lesser and greater intensities due to the uncertainties in the CO 2 atmospheric concentration at the origin of life [30]. Adapted from Michaelian and Simeonov [18].
In this paper we model the origin of life as a photochemical dissipative structuring process with reaction-diffusion occurring within a fatty acid vesicle floating at the ocean surface, permeable to the precursor molecules HCN, cyanogen and H 2 O, and under a UV-C photon spectrum arriving at the Earth's surface during the Archean (Figure 1). The concentrations of the relevant product molecules are determined as a function of ocean surface temperature, and as a function of incident light intensity, varying the amplitude but not the shape of the spectrum of Figure 1. Starting from a realistic environmental HCN concentration of 6 × 10 −5 M, we show that adenine and guanine concentrations of 2.4 × 10 −5 M and 7.6 × 10 −7 M respectively, could have been obtained in only 60 Archean days at neutral pH, 80 • C, and under a nominal UV-C flux of 125 W m −2 µm −1 at 260 nm. These final purine concentrations increase exponentially with temperature (roughly an order of magnitude for each 10 • C increase in surface temperature) and have a linear but complex dependence on light intensity.

Thermodynamics of Dissipative Structuring
We have identified the long wavelength part of the UV-C region (∼210-285 nm), plus part of the UV-B and UV-A regions (∼310-360 nm), of the solar spectrum as the thermodynamic potential which could have driven the molecular dissipative structuring [21,26], proliferation, and evolution relevant to the origin of life. This light prevailed at Earth's surface from the Hadean, before the origin of life (somewhere near the beginning of the Archean ∼3.9 Ga), and for perhaps as long as 1400 million years until the formation of an ozone layer at about 2.5 Ga [29,31,32] after natural oxygen sinks (e.g., volcanic reducing gases such as CH 4 , CO, H 2 S and Fe +2 ) became overwhelmed by organisms performing oxygenic photosynthesis, possibly fertilized by volcanic phosphorus [33].
The relevance of this particular region of the solar spectrum, corresponding to the Archean atmospheric window of transparency, to the dissipative structuring of the fundamental molecules is that longer wavelengths do not contain sufficient free energy to directly break double covalent carbon bonds, while shorter wavelengths contain enough free energy to destroy carbon based molecules through successive ionization or fragmentation.
A number of empirical evidences support our dissipative structuring conjecture for the fundamental molecules of life. First, the wavelength of maximum absorption of many of these molecules coincide with the predicted window in the Archean atmosphere ( Figure 1). Secondly, many of the fundamental molecules of life are endowed with peaked conical intersections [26] giving them broad band absorption and high quantum yield for internal conversion, i.e., extremely rapid (picosecond) dissipation of the photon-induced electronic excitation energy into vibrational energy of molecular atomic coordinates, and finally into the surrounding water solvent [21,34]. The most convincing evidence of all, however, is that many photochemical routes to the synthesis of nucleic acids [35], amino acids [36], fatty acids [23], sugars [37], and other pigments [18] from common and simple precursor molecules have been identified at these wavelengths and the rate of photon dissipation within the Archean window generally increases after each incremental step on route to synthesis, a behavior strongly suggestive of dissipative structuring in the non-linear non-equilibrium regime [21,23,26]. Finally, even small changes (e.g., tautomerizations) or additions to the fundamental molecules destroys these special optical properties.
The photochemistry involved in the non-equilibrium thermodynamics of microscopic dissipative structuring and proliferation of the fundamental molecules, and how this leads to the evolution of molecular concentration profiles with ever greater solar photon dissipation capacity (corresponding to increases in global entropy production) has been described in detail elsewhere (Michaelian [26]). Here we emphasize only those aspects most important to the dissipative structuring of the purines guanine and xanthine.
HCN is found throughout the cosmos [43,44]. Astronomical observational data indicates that HCN and water can be found together in unexpectedly high abundance in many hot star forming regions [45]. It may be, therefore, that HCN was incorporated into the bodies that formed the Earth and was gradually out-gassed to the surface, as is now suspected for Earth's water. Cyanide has, in fact, been found in abundance in carbonaceous chondrite meteorites [46]. Furthermore, if Earth's atmosphere during the Hadean and Archean had been weakly reducing [30] (see below), HCN production on Earth could have resulted from CH 4 photolysis by the solar Lyman alpha line (121.6 nm) in the upper atmosphere, leading to CH * radicals attacking N 2 , or, UV (145 nm) photolysis of N 2 and the resulting radicals attacking CH or CH 2 [47,48].
Cyanogen, NCCN or (CN) 2 , can be generated from HCN either photochemically [49] or thermally [50]. The Lyman-α line on atmospheric HCN will produce, through photolysis, the CN − radical with high quantum yield (reaction R177 from Table 1 of Zahnle [51]) and this radical gives cyanogen after interaction with a second HCN molecule (reaction R199 from Table 1 of Zahnle). Both HCN and cyanogen have been found in the atmosphere of Titan [48] where they probably derive from similar UV photochemical processes in the upper or middle atmosphere.
In Michaelian [26] we have addressed the issue of the expected low concentration of HCN in the Archean oceans leading to a low polymerization rate of HCN with respect to its rate of hydrolysis. Traditional approaches to the origin of life have "resolved" this issue by invoking eutectic freezing to increase the solute HCN concentration to values sufficient for biasing the reactions towards polymerization [52][53][54]. This low temperature approach, however, reduces significantly all thermal reaction rates. Instead, we have emphasized the existence of an ocean surface microlayer of high organic concentration, 10 4 or 10 5 times greater than that of bulk water (see Michaelian [26] for discussion and references). We also assume the existence of fatty acid vesicles of ∼100 µm diameter which would allow the incident UV-C light, as well as the small precursor molecules, HCN, cyanogen, and H 2 O, to permeate the vesicle bi-layer wall (Figure 2), while trapping inside the reaction products due to their larger sizes and larger dipole moments (Table 1). This allows the intermediate and product molecules, as well as the heat from UV-C photon dissipation, to accumulate within the vesicle.  (Table 1). Table 1. Nomenclature, chemical formula, abbreviation (used in the text and in kinetic equations), label in Figure 3, wavelength of maximum absorption λ max (Figure 1), molar extinction coefficient at that wavelength max , electric dipole moment µ, and the topological polar surface area (TPSA), of the molecules involved in the photochemical synthesis of the purines. Values marked with "*" are estimates, obtained by comparing to similar molecules, for cases where no published data have been found. The existence of amphipathic fatty acid hydrocarbon chains, which through Gibb's free energy minimization spontaneously form lipid vesicles at the ocean surface, is a common assumption in origin of life scenarios [66,67]. Deamer and Georgiou [68,69] have discussed in detail the many roles of lipid membranes in life's origin. In order to maintain vesicle integrity at the high surface temperatures considered here at the beginning of the Archean, ∼80 • C [70][71][72][73][74][75], these fatty acids would necessarily have been long (∼18 C atoms) and cross linked through UV-C light, which improves stability at high temperatures and over a wider range of pH values [23,76], and provides resistance to salt flocculation. In the anoxic Archean ocean surface environment, UV-C light would promote cross-linking between consecutive fatty acids rather than cause oxidative damage.
Archean fatty acids may have been aromatic, giving them improved stability at higher temperature for shorter lengths due to a van der Waals stacking interaction between the rings of the interior and exterior layers of the bilayer vesicle. There is, in fact, a predominance of 16 and 18 carbon atom fatty acids in the whole Precambrian fossil record [77,78]. Today, the aromatic fatty acid 10-phenyltridecanoic acid (16 carbon atoms) and other ω-phenyl fatty acids, such as 12-phenyldodecanoic acid (18 carbon atoms) and 14-phenyltetradecanoic acid (20 carbon atoms), are found in the highly salt resistant marine bacteria V. alginolyticus [79], and ω-cyclohexyl fatty acids are found today in some acidophilic-thermophilic bacteria isolated from hot springs [80].
From the perspective of our dissipative structuring scenario, aromatic fatty acids are also interesting because the conjugated aromatic rings have wide band UV absorption cross section within the 210-285 nm region and a conical intersection for rapid internal conversion to the ground state, giving them large photon dissipation efficacy and photochemical stability under this UV-C light. Indeed, we have suggested that fatty acids, just like the nucleic acids and aromatic amino acids, could have been dissipatively structured from HCN and CO 2 saturated water under this UV-C light at moderate temperatures on the ocean surface [23] because these optical characteristics are the hallmarks of dissipative structuring [21,26].

The Photochemical Route to the Purines
Our proposed photochemical route to the purines from the precursors H 2 O, HCN and cyanogen is based on the experimental data of Ferris and Orgel [35] and Sanchez et al. [6], as well as the quantum mechanical calculations of Boulanger et al. [60], and is given in Figure 3.  (Ferris and Orgel (1966) [35,60]), and a cyanogen molecule for guanine (Sanchez et al. (1968) [6]). Four molecules of HCN are transformed into the smallest stable oligomer (tetramer) of HCN, known as cis-2,3-diaminomaleonitrile (cis-DAMN) (2), which, under a UV-C photon flux isomerizes into trans-DAMN (3) (diaminofumaronitrile, DAFN) which by absorbing two more UV-C photons transforms into an imidazole intermediate, 4-amino-1H-imidazole-5-carbonitrile (AICN) (7). Hot ground state thermal reactions with another HCN molecule, or its hydrolysis product formamide (or ammonium formate), leads to the purine adenine (8). Hydrolysis of adenine (8) leads to hypoxanthine (9). Hydrolysis of AICN (7) leads to AICA (10) and a hot ground state thermal reaction with a cyanogen molecule leads to guanine (11). Hydrolysis of guanine (11) leads to xanthine (12). This is a microscopic dissipative structuring process under UV-C light [21,26] which ends in the purines; adenine, hypoxanthine, guanine and xanthine, all UV-C pigments with large molar extinction coefficients, at 260, 250, 252 and 268 nm respectively, and peaked conical intersections providing rapid dissipation of photons around the wavelength of maximum intensity of the Archean surface solar UV-C spectrum ( Figure 1).
The route to guanine results from the hydrolysis of 4-aminoimidazole-5-carbonitrile (AICN) to form 2-amino-3-iminoacrylimidoyl cyanide (AICA) (step (7) to (10) of Figure 3) which then combines through a thermal chemical reaction with cyanogen (NCCN) to form guanine plus a HCN molecule [6]. The reactions on route to adenine and guanine are in competition with hydrolysis and UV lysis, and these relative rates are dependent on concentrations, temperature, pH, the presence of metal ion-and organic intermediate product catalysts, as well as the wavelength dependent intensity of the incident UV spectrum (Section 4). Some of these complexities have been studied by Sanchez et al. [5,6].
In the following subsection we present our simplified out-of-equilibrium kinetic model for 5HCN → adenine and 4HCN + H 2 O + NCCN → guanine + HCN (Figure 3) photochemical reaction system occurring within a fatty acid vesicle floating within the high organic concentration microlayer of a hot (∼80 • C [70][71][72][73][74][75]) Archean ocean surface under the UV surface spectrum of Figure 1. Favorable comparisons of our simulated photochemical production rates to experimental data have been given in Michaelian [26] (this Special Issue).

The Kinetic Equations
The kinetic equations for the chemical and photochemical reactions are resolved numerically. For this perturbed non-linear reaction-diffusion system it will be shown in Section 4 that various stationary state solutions exist with spatial symmetry breaking resulting in the highest concentration of purines at the center of the vesicle. This would then facilitate a subsequent UV-C polymerization of nucleobases into short oligos (of course assuming the possibility of UV-C-assisted synthesis of ribose from similar precursor molecules [42,81], a UV-C induced nucleoside formation from a base and ribose [82], and a temperature [83] or formamide-catalyzed [84] phosphorylation of the nucleosides [85,86], but these reactions will not be considered here). This stationary state coupling of reactions to diffusion, leading to symmetry breaking (i.e., particular regions of high concentration of the products), was shown to occur for purely thermal reactions with different activator and inhibitor diffusion rates by Turing [87] and studied more generally as dissipative structures under the framework of Classical Irreversible Thermodynamic (CIT) theory by Glansdorff and Prigogine [88].
Nomenclature, chemical formula, and abbreviations used throughout the text for the concentrations of the participating chemical species involved in the photochemical reactions leading to the purines shown in Figure 3 are given in Table 1. Their photon extinction coefficients, polar surface area and dipole moments (related to permeability of the vesicle wall), are also given.

# Reaction Reaction Constants
The following is a description of the reactions given in Table 2 by reaction number. A detailed description is given only for those reactions relevant to the production of guanine and xanthine. For a detailed description of all other reactions, see Michaelian [26].  Table 6 of Sanchez et al. [5]). The catalytic effect of trans-DAMN on the tetramization of HCN was incorporated into the model by reducing the energy of the activation barrier such as to give the same amplification factor of 12 due to the catalytic effect of the inclusion of 0.01 M trans-DAMN in the HCN solution observed in the experiments of Sanchez et al. [5] at a temperature of 20 • C.
The catalytic effect of cyanogen (Cg) is taken to be the same as that of trans-DAMN (T) since Sanchez et al. determined both of these to be strong catalysts [5]. As mentioned above, cyanogen is a precursor needed for guanine obtained from the Lyman-α line (121.6 nm) on atmospheric HCN giving, through photolysis, the CN − radical with high quantum yield (reaction R177 Table 1 of Zahnle [51]), and this can form cyanogen (CN) 2 by interacting with a second HCN molecule (reaction R199, Table 1 of Zahnle); 8. Trans-DAMN also acts as an auto-catalyst for its own thermal production from 4HCN.
Cyanogen is also an effective catalyst for this reaction (Table 6  These reactions, with large quantum efficiencies, represent the bulk of the flow of energy from the incident UV-C spectrum to the emitted outgoing infrared ocean surface spectrum and therefore contribute most to photon dissipation, or entropy production. To simplify the kinetic equations for the photochemical reactions listed in Table 2, we assume that the molecules only absorb within a region ±10 nm of their maximum absorption wavelength λ max and that this absorption is at their maximum molar extinction coefficient max (Table 1), and finally that these wavelength regions do not overlap for calculating the shadowing effect of molecular concentrations at higher positions in the vesicle.
We assume that the vesicle is at the ocean surface and the depth coordinate is divided into i = 20 bins of width ∆x = 5 µm and the time step for the recursion calculation for the concentrations is 10 ms. The recursion relation for the factor of light intensity L λ (i, C) for a concentration C of the molecule, at a depth x(i) = i · ∆x below the ocean surface will be, where α λ is the absorption coefficient of water at wavelength λ and λ is the molar extinction coefficient of the particular absorbing substance which has concentration C(i) at x(i).
The kinetic equations giving the increments in concentration after each time step ∆t ≡ dt, for use in a discrete recursion relation, at a depth x below the ocean surface are determined from the reactions listed in Table 2 to be the following: where the differentials are calculated discretely (e.g., dH/dt ≡ ∆H/∆t) and all concentrations are calculated at discrete time steps of ∆t = 10 ms and the calculated value of the change (e.g., ∆H(j)/∆t) for time step j is summed to the previous value (e.g., H(j − 1)). The day/night factor d is equal to 1 during the day and 0 at night. I 220 , I 298 , I 313 , I 275 and I 250 are the intensities of the photon fluxes at 220, 298, 313, 275 and 250 nm respectively (Figure 1).
λ are the coefficients of molar extinction for the relevant molecule at the corresponding photon wavelengths λ.

Vesicle Permeability and Internal Diffusion
The permeability of the fatty acid vesicle will decrease both with increasing size of the molecule and with increasing size of its electric dipole moment, and increase with temperature. The final and intermediate product molecules have large sizes and large dipole moments compared with the precursors, implying a tendency towards entrapment within the vesicle. We assume that the vesicle is completely permeable to H 2 O, HCN (H) formimidic acid (Fa) and cyanogen (Cg) but impermeable to all the other intermediate and final products. Note that ammonium formate (produced within the vesicle in reaction #4) would be in its ionic form and therefore also unable to cross the fatty acid membrane (permeability across lipid membranes are reduced by orders of magnitude if the molecules are polar or charged [101]).
The diffusion constant D Y for the molecule Y within the inner aqueous region of the vesicle will depend on the interior solution viscosity, which is dependent on the concentration of organic material. Studies of inter-cellular diffusion of nucleotides indicate three factors influencing diffusion rates, besides temperature, at high solute densities; the viscosity of the solvent, collision rate dependent on solute concentration, the size of the molecules, and the binding interactions between molecules [102]. The diffusion constant of adenine in pure water has been determined to be D A = 7.2 × 10 −6 cm 2 s −1 [103] while the measured diffusion rates in the cytoplasm of a cell can be as low as 1.36 × 10 −6 cm 2 s −1 [102].
Films of organics and trace metals, with a high density of lipids and other hydrocarbons, produced, for example, by the ultraviolet spectrum of Figure 1 on CO 2 saturated water [23], would be expected on the ocean surface during the Archean. Diffusion constants within this sea surface microlayer would then be significantly smaller than for bulk water. Diffusion rates inside the vesicle will depend on the amount of organic material captured at the air/water interface during the formation of the vesicle (which would have considerable spatial variability) and on the amount of ongoing organic synthesis within the vesicle.
The Stokes-Einstein diffusion constant of uncharged spherical molecules is given by: where T is the absolute temperature, η the dynamic viscosity of the solvent, and r the radius of the particle (assumed spherical). For polar solvents such as water, molecules with a charge or dipole moment experience further drag due to the charge-dipole or dipole-dipole interaction [104]. All diffusion constants can, therefore, be approximated relative to that measured for adenine at 300 K through the formula: where A A is the polar surface area of adenine, and µ A and µ Y are the dipole momenst of adenine and the molecule Y respectively ( Table 1). The factor of 2 in the equation weighs the importance of including the dipole moment in the diffusion calculation (Equation (17)), which is complex but relatively unimportant to final results. Our analysis is performed with the diffusion constant for adenine found in present day cytoplasm [102] (i.e., D A (300) ≈ 1 × 10 −6 cm 2 s −1 ). Using the values given in Table 1 for the molecular dipole moments and their polar areas, and Equation (18), in Table 3 we list the results for the diffusion constants of all relevant molecules at 80 • C (353.15 K) with respect to that for adenine at 300 K, D A (300). Table 3. Diffusion constants at 80 • C (353.15 K) relative to that of adenine at 300 K (i.e., d Y = D Y (353.15)/D A (300) obtained from Equation (18). In this paper, D A (300) is taken to be 1 × 10 −6 cm 2 s −1 . Cyclical boundary conditions are assumed for diffusion, except for HCN (H), formimidic acid (Fa) and cyanogen (Cg) which, because of their small size and small dipole moment, can permeate the vesicle wall and, therefore, at the wall they are given their fixed value specified in the initial conditions of the environment outside the vesicle (see below). The second order derivatives for calculating the diffusion were obtained using the second order finite difference method with double precision variables.

Initial Conditions
Miyakawa, Cleaves and Miller [91] estimated the steady state bulk ocean concentration of HCN at the origin of life assuming production, through electric discharge on a neutral atmosphere containing some methane, of radicals which attack N 2 , leading to an input rate to the oceans of 100 nmole cm −2 y −1 , and loss of HCN due to hydrolysis and destruction at submarine vents with a 10 million year cycling time of all ocean water for an ocean of 3 Km average depth. For an ocean of pH 6.5 and temperature of 80 • C, they obtained a value of [HCN] = 1.0 × 10 −10 M [91].
Not included in the Miyakawa et al. calculation, however, is, as mentioned above, the fact that HCN can also be produced through the solar Lyman alpha line (121.6 nm) photo-lysing N 2 in the upper atmosphere giving atomic nitrogen which then combines with CH and CH 2 to give HCN, or through 145 nm photolysis of CH 4 leading to a CH * radical which attacks N 2 to give HCN [47]. Including this UV production would increase the input of HCN concentration to the oceans by a factor of at least 6 [51,105,106]. More importantly, the first ∼100 µm of the ocean surface is now known to be a unique region (the hydrodynamic boundary layer) where surface tension leads to enriched organics with densities up to 10 4 or 10 5 times that of organic material in the water column slightly below [107]. Langmuir circulation, Eddy currents, and the scavenging action of bubbles, tend to concentrate organic material into this surface film. Trace metal enhancement in this microlayer can be one to three orders of magnitude greater than in the bulk [107,108].
This high density of organic material trapped through hydrophobic and ionic interactions at the ocean surface leads to significantly lower diffusion rates within the surface microlayer as compared to the ocean bulk [107]. Little diffusion and turbulence imply little mixing. The ocean microlayer is, therefore, a very stable layer which would not be recycled through ocean vents and thus allowed to accumulate. Finally, although HCN is very soluble in bulk water, molecular dynamic simulations have shown that it concentrates to about an order of magnitude larger at the air-water interface due to lateral HCN dipole-dipole interactions, and that HCN evaporates at lower rates than water [109]. Therefore, rather than assuming the low HCN bulk concentrations of Miyakawa et al., we instead consider a higher initial surface concentration of HCN (H) and cyanogen (Cg) of 6 × 10 −5 and for formamide (F) and formimidic acid (Fa) of 1 × 10 −5 , the latter resulting from a photochemical tautomerization of formamide, the hydrolysis product of HCN (reactions #1 and #2 of Table 2). We also allow for perturbation of the system by considering the probable existence of small and sparse patches of much higher concentration of the precursors, justified by the characteristics of the ocean microlayer and the dipole-dipole interaction between HCN molecules. The system is perturbed by a single patch of relatively high concentration (0.1 M) of HCN, cyanogen (Cg), and formimidic acid (Fa) (see Michaelian [26] for a description) into which our vesicle is assumed to drift (at day 10.4) only once for 2 min during 60 Archean days. The initial concentrations of all other reactants inside the vesicle (assumed impermeable to these) are taken to be 1.0 × 10 −10 M.
Our system is under a diurnal 8 h constant flux of radiation followed by an 8 h period of darkness during which thermal reactions proceed, but not photochemical reactions. Different surface UV-C light intensities will be considered due to the uncertainty in the Archean atmospheric gas concentration profile. The "nominal" intensity corresponds to a weakly reducing and self-shielding (for λ < 220 nm) atmosphere as proposed by Sagan and Chyba [30], which assumes that most of the greenhouse warming (required to solve the Faint Young Sun paradox [30]) is provided by CH 4 and a small amount of ammonia (mixing ratio 1 × 10 −5 ) which is shielded by the products resulting from CO 2 and CH 4 photolysis (e.g., HCN and aldehydes). The corresponding integrated (210-285 nm) UV-C light intensity for this "nominal" scenario ( Fig. 1) is 4 W m −2 . On the other hand, Cnossen et al. [32] suggest a neutral atmosphere in which the partial pressure of CO 2 has to be very large, 0.24 bar, in order to achieve the required global warming with the faint young sun. This results in a much reduced integrated (210-285 nm) UV-C surface flux of 0.10 W m −2 . We therefore perform our simulations for a number of different incident photon intensities covering both atmospheric profile scenarios.

Evolution of the Concentration Profile
In Michaelian [26] we have validated our photochemical reaction-diffusion model by comparing results to experimental data for the time dependence of a number of different intermediate products.
In Figures 4-7 we present the concentrations of the relevant molecular products within the vesicle as a function of time in Archean days (16 h) in the dissipative photochemical synthesis of adenine and guanine, obtained by solving simultaneously the differential kinetic Equations (2)- (16), for the initial conditions and diffusion constants listed in the figure captions.
The concentration profiles of the molecules evolve over time because of accumulation of photo-products within the vesicle and due to a deliberate external perturbation of the non-linear system at 10.4 days (0.1 M of HCN, cyanogen, and formimidic acid) which leads it to a new stationary state in which conversion of the environmental precursor molecules HCN and cyanogen into adenine and guanine occurs at greater rates. This leads to greater photon dissipative efficacy of the system, i.e., to a concentration profile of the product molecules which dissipates more efficiently the incident UV-C spectrum.
In Figure 6 we plot the concentration profiles of the photochemical products obtained also at 80 • C but with 1/100 of the incident light flux (i.e., UV-C flux 1.25 W m −2 µm −1 at 260 nm, integrated flux 210-285 nm = 0.040 W m −2 ). This would correspond approximately to the illumination on Earth's surface for a very neutral atmosphere if the only greenhouse gas to resolve the Faint Young Sun paradox (and increase the surface temperature to about 80 • C) was CO 2 (see Section 3.2), assuming the standard solar model. It is interesting to note that the adenine concentration after 60 days actually increases slightly with this lower surface UV-C light flux, while all other purine concentrations decrease. This is because less UV-C light implies more formamide (F) being available (less of photoreaction #2) and thus more ammonium formate (Af) (reaction #4). Formamide (F) is a catalyst in the most important thermal reaction #13 which has a low activation energy barrier for converting AICN (I) and ammonium formate (Af) into adenine (A).
Our resulting adenine concentrations are between about one and two orders of magnitude greater than those of guanine, which is consistent with experimental data of Levy et al. [10] using ammonium cyanide as a precursor under long time eutectic freezing conditions. This greater difficulty of producing guanine abiotically through both processes, eutectic freezing and photochemical dissipation, may explain why the major energy storage molecule of life is adenine triphosphate (ATP) rather than guanine triphosphate (GTP).
Given the fixed concentrations of HCN (H), cyanogen (Cg) and formimidic acid (Fa) in the environment, to which the vesicle is permeable, photochemical reactions occur during daylight hours (marked as violet colored sections of the horizontal dashed line labeled as D/N). This gives rise to the observable diurnal oscillations in the concentrations of trans-DAMN (T) and AICN (I) since these are direct products of photochemical reactions.  . This sudden impulse in HCN, cyanogen, and formimidic acid concentration gives rise to rapid increases in all concentrations within the vesicle, in particular formamide (F), the hydrolysis product of HCN, which is an important catalyst for reaction #13 which produces adenine (A) from AICN (I) (see Table 2) and this reaction has a low activation barrier energy. Ammonium formate (Af) is used up in this reaction so its concentration decreases rapidly after the perturbation. More importantly, however, immediately after the perturbation there is a greater production of trans-DAMN (T) in the vesicle and since T acts as a catalyst for the polymerization of HCN (H) (reactions #7 and #8), this produces greater metabolism of HCN into DAMN within the vesicle, leading to a stronger diffusion of HCN into the vesicle from the outside environment. In other words, the reason that a short impulse of HCN, cyanogen, and formimidic acid give rise to an important increase in the rate of production of adenine and guanine is that the vesicle's semi-permeable wall, together with the set of equations describing the photochemical and chemical reactions, Equations (2)- (16), form a non-linear (auto-and cross-catalytic) system with multiple solutions at any given time. The perturbation at 10.4 days causes the system to leave the attraction basin of one solution determined by its initial conditions and evolve towards a different, and more probable, stationary state with a higher rate of production of adenine and guanine (slope of the black and blue traces respectively of Figure 5). The new stationary state is more stable than the initial state because its concentration profile is more dissipative, i.e., with more molecules having conical intersections dissipating the absorbed photon energy rapidly into heat, and, therefore, there is less probability of photochemical reactions reverting the concentration profile to the initial state. There is thus a thermodynamic driving force for evolution to the new stationary state after the perturbation due to the system obtaining a greater photon dissipative capacity [26]. For systems where local thermodynamic equilibrium (in space and time) can be defined (see Section 3 of Michaelian [26]), this implies that the entropy production of the system increases (see below). This example of the dissipative structuring a UV-C pigment concentration profile, elucidates the physics and chemistry behind a thermodynamic selection relevant to biological evolution at its earliest stages [26]. The temperature dependence of the concentrations of the product molecules obtained after 60 Archean days is given in Figure 8. Ammonium formate (Af) is produced by the hydrolysis of first HCN (H) to formamide (F) (reaction #1) and then hydrolysis of formamide to Af (reaction #4). Both of these reactions have high activation energies and this results in Af only being produced in important quantities at temperatures greater than about 80 • C. Most of the adenine (A) production occurs through reaction #13 which consumes Af and therefore high temperatures are important to the production of adenine. The hydrolysis of AICN (I) to AICA (L) (step (7) to (10) of Figure 3) also has a high activation energy barrier (reaction #12), also implying high temperatures necessary for the production of hypoxanthine (Hy), guanine (G) and xanthine (Xa), as can be seen by the slopes of their traces on Figure 8. In Figure 9 we plot the concentrations of the product molecules obtained at a temperature of 80 • C after 60 Archean days as a function of the incident light intensity in the 210-285 nm region, by scaling the spectrum given in Figure 1. Maximum adenine concentration is obtained for fluxes of about 1.25 W m −2 µm −1 at 260 nm, however, at this flux, the concentration of guanine is 3 orders of magnitude smaller. At the nominal intensity of 125 W m −2 µm −1 at 260 nm (Figure 1), the guanine concentration is about 30 times lower than that of adenine which remains high at 2.4 × 10 −5 M. In Figure 10 we plot the entropy production as a function of time in Archean days due solely to the absorption and dissipation of photons into heat through the conical intersections of the product molecules associated with reactions 19 to 28 of Table 2 for the particular concentration profile at the given time. For a description of the entropy production and its evolution through time in relation to the non-equilibrium thermodynamic Glansdorff-Prigogine universal evolutionary criterion, see reference [26]. Figure 10. The entropy production as a function of Archean days for our UV-C photochemical dissipative structuring system within a vesicle floating at a sea surface at a temperature of 80 • C leading to the purines. Entropy production increases monotonically as photochemical reactions convert HCN into the photon dissipative product molecules, including adenine and guanine. During the day, entropy production is due to the dissipation of the UV-C spectrum into heat by the corresponding product concentration profile. During the night, entropy production is due to thermal chemical reactions, however, this entropy production is small and not included in the figure. At 10.4 Archean days, the system is perturbed and moves into a new stationary state with the entropy production increasing by more than 4 orders of magnitude and continues to increase as the concentration profile evolves (see Figure 5).

Discussion
As Figure 8 illustrates, high temperatures are very important to a building up high concentrations of purines within the vesicle. High temperatures also foment phosphorylation with phosphate salts using formamide as a catalyst, favoring the formation of acyclonucleosides and the phosphorylation and trans-phosphorylation of nucleosides which only occurs efficiently at temperatures above 70 • C [83,84]. High temperatures are similarly important for the abiotic production of amino acids and fatty acids. Such high temperatures are, in fact, consistent with the geochemical data of the early Archean [70][71][72][73][74] and with the genome comparison of strongly conserved ribosomal RNA [75].
Important concentrations of all purines buildup within the vesicle within only about 60 Archean days under the nominal environmental conditions considered here. There is no need to begin with unrealistically large initial concentrations of HCN by invoking eutectic freezing, and there is no need for highly alkaline conditions in order to favor HCN polymerization over hydrolysis.
It is instructive to compare our overall non-equilibrium results obtained with the model of UV-C dissipative structuring of adenine and guanine from HCN within a lipid vesicle at 80 • C with the quasi-equilibrium experiments of Ferris et al. [39] using high HCN concentration, alkaline pH, and room temperature. Starting with a concentration of HCN (0.1 M) in water (pH 9.2), and allowing this solution to polymerize in the dark at room temperature for a 7 months, and then subjecting these polymers to hydrolysis at 110  (Figure 4).
As detailed in Michaelian [26], besides entrapment of product molecules inside the vesicle, another concentration mechanism for these molecules that may arise is the coupling between reaction and diffusion in the non-linear regime which leads to the breaking of spatial symmetry. For low diffusion rates, the homogeneous stationary state is no longer stable with respect to space dependent perturbations and intermediate and final products may become preferentially concentrated and be consumed in the center of the vesicle (see Michaelian [26]).
The dependence of the final concentrations on the light intensity is complex (Figure 9) but important values of both adenine and guanine can be obtained for surface intensities corresponding to the nominal value (for a weakly reducing atmosphere) of 125 W m −2 µm −1 at 260 nm ( Figure 1) and for intensities either one order of magnitude smaller or greater. This result may have a dependence on the particular shape of the incident spectrum, but this was not studied here.
Inorganic catalysts have not been included in our reaction scheme. These can increase significantly the rate of purine production. For example Cu +2 ions increase the rate constant for the conversion of HCN (H) to cis-DAMN (C) [5] (reaction # 5). Metal ions would have been in high abundance at the ocean surface microlayer [107,110].
The dissipative structuring of the pyrimidines under the same Archean UV spectrum ( Figure 1) will be considered in a future article.

Conclusions
We have presented a simple model of UV-C photochemical dissipative structuring of the purines within a lipid vesicle permeable to the common precursor molecules H 2 O, HCN, cyanogen, and formimidic acid, but impermeable to the reaction products, floating on a hot Archean ocean surface. The chemical and photochemical kinetics are based on published experimental and quantum mechanical ab initio data. The surface physical conditions are consistent with biochemical and geochemical fossil evidence from the early Archean.
High temperatures are important to obtaining high purine final concentrations which increase by roughly an order of magnitude for each 10 • C increase in surface temperature. Important concentrations of all purines buildup within the vesicle within only about 60 Archean days. There is no need to assume unrealistically high concentrations of HCN by invoking eutectic freezing or highly alkaline conditions. Our resulting adenine concentrations are between 1 and 2 orders of magnitude greater than those of guanine.
The final product concentrations have a complex dependence on light intensity. The largest concentrations of adenine occur for light intensities of between 1 and 125 W m −2 µm −1 at 260 nm, with guanine concentrations reaching about 1/30 of adenine concentration at 125 W m −2 µm −1 at 260 nm ( Figure 9). This adenine to guanine ratio is consistent with experimental data of Levy et al. [10] employing long time eutectic freezing conditions. This greater difficulty of producing guanine abiotically, through both eutectic freezing and photochemical dissipation, may explain the fact that the major energy storage molecule of life is adenine triphosphate (ATP) rather than guanine triphosphate (GTP).
Perturbations, caused by the vesicle floating into isolated patches of higher concentration precursors floating on the ocean surface, can provoke the non-linear auto-catalytic system into new states of higher product productivity with a greater "metabolism" of both precursor molecules (HCN and cyanogen) from the environment and the incident photons, leading to a concentration profile of much greater photon dissipative efficacy. This evolution towards stationary states of usually higher photon dissipative efficacy is the essence of biological evolution which continues today, but now involving the whole solar spectrum and the entire biosphere.
Evolution during the Archean was towards concentration profiles of product molecules with an absorption maximum closer to the peak intensity of the incident UV-C spectrum at 260 nm ( Figure 1) and with peaked conical intersections to internal conversion, both increasing the overall efficacy of the dissipation of the incident UV-C solar flux.
For low diffusion rates, there can be significant coupling of reactions with diffusion, leading to non-homogeneous distributions of some of the intermediate products, in particular, with greater concentration of the products at the center of the vesicle (see Michaelian [26]). Such spatial symmetry breaking could facilitate further dissipative structuring, such as polymerization of the nucleobases into nucleic acid.
Dissipative structuring, dissipative proliferation, and dissipative selection (as explained in detail in Michaelian [26]) are the necessary and sufficient components for an explanation in physical-chemical terms of the synthesis, proliferation, and evolution of organic molecules on planets, comets, asteroids, and interstellar space [18] and, in particular, for explaining the origin and evolution of the fundamental molecules of life on Earth.