Received Date: October 21, 2016; Accepted Date: February 2, 2017; Published Date: February 8, 2017
Citation: Michaelian K, Simeonov A (2017) Thermodynamic Explanation for the Cosmic Ubiquity of Organic Pigments. Astrobiol Outreach 5:156. doi:10.4172/2332-2519.1000156
Copyright: © 2017 Michaelian K, et al. This is an open-access 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.
Visit for more related articles at Journal of Astrobiology & Outreach
There is solid evidence for the occurrence of large amounts of organic material in the cosmos, particularly in the form of aromatic compounds. These molecules can be found on the surface of Earth and Mars, in the atmospheres of the larger planets and on many of their satellites, on asteroids, comets, meteorites, the atmospheres of red giant stars, interstellar nebulae, and in the spiral arms of galaxies. Many of these environments are expected to be of low temperature and pressure, implying that the Gibb’s free energy for the formation of these complex molecules should be positive and large, suggesting that their existence could only be attributed to non-equilibrium thermodynamic processes. In this article we first review the evidence for the abundance of these molecules in the cosmos and then describe how the ubiquity can be explained from within the framework of non-equilibrium thermodynamics on the basis of the catalytic properties of these pigment molecules in dissipating photons of the ultraviolet and visible emission spectra of neighboring stars, leading to greater local entropy production. A relation between the maximum wavelength of absorption of these organic pigments and the corresponding stellar photon environment provides a guide to determining which aromatic compounds are most probable in a given stellar neighborhood, a postulate that can be verified on Earth. It is suggested that at least some of the baryonic dark matter may be associated with these molecules which emit in the extreme infrared with many, but weak, emission lines, thus so far escaping detection. This thermodynamic explanation for the ubiquity of these organic molecules also has relevance to the possibility of life, both as we know it, and as we may not know it, throughout the universe.
In 1937, the detection of the methylidyne (CH) radical in early type stars (O, B and A stars) [1,2] ushered in an era of astrochemical discoveries. It was simultaneously the first molecule and the first organic compound to be ever discovered in extrasolar space.
An upsurge of molecular discoveries followed during the sixties and seventies, with the discovery of ammonia , water , formaldehyde , methanol , and carbon monoxide . Of great importance to the question of the origin of life was the detection of N-containing organics like hydrogen cyanide , methylamine , and cyanamide , since it had been shown that HCN in the presence of water and a free energy source is sufficient to produce the DNA purines and pyrimidines [11-14]. It therefore came as no surprise when the nucleobases were discovered in carbonaceous meteorites, along with some unusual and terrestrially rare analogs (purine, 2,6- diaminopurine, and 6,8-diaminopurine), indicative of their extraterrestrial origin [15,16]. Another important discovery was that of interstellar acetylene  a precursor molecule to benzene and hence other aromatic and polyaromatic organic molecules .
Pivotal to the search and the accumulation of data has been the development of millimeter-wave and infrared spectroscopy which, through the detection of the stretching and bending vibrational modes of covalent bonds, as well as rotational transitions, has confirmed the existence of a vast array of these molecular types in space .
The large abundance of the element carbon (fourth most abundant element in the universe after H, He and O ), its unique covalent bonding flexibility, and the abovementioned pioneering discoveries of low molecular weight organics which can act as chemical precursors to other molecules, have spurred the search for larger, more complex organics in space.
In some unconventional and truly startling comparisons, Fred Hoyle and Chandra Wickramasinghe, the founders of the theory of panspermia, showed how the extinction spectra over the ultraviolet and visible of interstellar nebulae can be fitted well by the absorption spectra of freeze dried algae and bacteria . More recent higher resolution data refute the quality of the fit , but there can be little doubt that the nature and complexity of the absorbing material in the interstellar clouds must be somehow analogous to that of the biological macromolecules. The Extended Red Emission (ERE) feature, for example, which is a photoluminescence process powered by UV photons with a broad peak emission wavelength around 650-800 nm , most closely parallels the optical behavior of chlorophyll and other porphyrinic bio-pigments, that have absorbance in the near-UV/ blue and fluorescence in the red .
Moreover, the ERE is only one of a family of universal spectroscopic features, which, although not yet fully explained, are typically regarded as spectral “fingerprints” of complex organic compounds throughout the cosmos. These include: the Unidentified Infrared Emission (UIE) bands, the Diffuse Interstellar absorption Bands (DIBs), the 217 nm feature, and the 21 and 30 micron emission features .
In particular, the UIE bands, which are emission bands originating from numerous galactic and extragalactic astronomical environments, have traditionally been considered as the spectral signatures of polycyclic aromatic hydrocarbons (PAHs) [28,29], although newer spectroscopic analyses [30,31] revive earlier proposals  for even more complex chemical structure of the carriers of these bands; Mixed Aromatic/Aliphatic Organic Nanoparticles (MAONs), very reminiscent of the complex macromolecular structure of kerogen and coal on Earth, and the Insoluble Organic Matter (IOM) of carbonaceous chondrites .
There is thus unequivocal evidence for the formation of large and complex organic molecules throughout the universe. However, these complex molecules are prone to photo-ionization and destruction by charged particles, x-rays, and the same high-frequency UV photons that probably helped formed them. Therefore, their ubiquity, and thus the probable existence of life as we know it in the cosmos, requires an explanation that goes beyond near-equilibrium reaction pathways of formation. That these molecules are the result of non-equilibrium thermodynamic processes can be seen from the fact that their occurrence is dependent upon a continuous flow of UV photons from the neighboring star. However, this paper goes further to show that the ubiquity of organic pigments throughout the cosmos, wherever there exists UV light interacting with material, is a result of the non-linear, non-equilibrium nature of the photochemical process leading to the organic molecules. Specifically, it is shown that these molecules can be considered as microscopic self-organized dissipative structures which form and proliferate to dissipate the generalized thermodynamic force; the stellar photon potential formed between the high energy photon outflow of a neighboring star and the low energy (microwave) photon cosmic background of space. In other words, these pigment molecules owe their existence and proliferation to the entropy production that they perform by dissipating high energy photons into heat.
Evidence for the cosmic ubiquity of organic pigments
Infrared observations by Gillett et al.  and Russell et al. [35,36] revealed that the IR spectra of bright astronomical objects with associated gas and dust particles, such as HII regions and planetary nebulae, are dominated by strong, broad, and resolved emission bands. Since the carrier of the bands remained unidentified for almost a decade, they were dubbed the Unidentified Infrared Emission (UIE or UIR) bands. The fact that their intensity showed a direct correlation with carbon abundance within the source naturally implied a carbonbased carrier.
Knacke  and Duley and Williams  first suggested aromatic carriers, since the bands peak at wavelengths corresponding to the stretching and bending modes of various CH and CC bonds in aromatic hydrocarbons, so they also became known as the Aromatic Infrared Bands (AIB or AIR).
With the launch of the Infrared Space Observatory (ISO) mission in 1995 the UIE features were observed with better spectral resolution, covering a broader spectral range, and in a larger sample of sources. They were found to dominate the mid-IR emission of almost all astronomical objects except for asymptotic giant branch (AGB) stars and deeply embedded young stellar objects . It was estimated that around 20-30% of the infrared radiation of our galaxy is being emitted in these infrared emission bands , and up to 20% of the total radiation output of distant starburst galaxies [41,42], indicating that their carrier must be an extremely common and abundant cosmic material.
Solar System objects, such as Interplanetary Dust Particles (IDPs), carbonaceous meteorites, micrometeorites, Martian rocks , comets , Titan’s upper atmosphere  are also known to display these infrared features.
A wide variety of hydrocarbon and carbonaceous materials comprising aromatic units have been proposed as possible carriers of the AIB bands, including; Hydrogenated Amorphous Carbon (HAC), Quenched Carbonaceous Composites (QCC), Polycyclic Aromatic Hydrocarbons (PAHs), soot and carbon nanoparticles, fullerenes and fullerene-like particles, coal and kerogen, petroleum fractions etc. . Among all these models, the PAH hypothesis developed by Leger and Puget  and Allamandola et al. [28,29] has gained most attention, therefore the UIE bands are usually referred to in the literature simply as the PAH bands.
Polycyclic Aromatic Hydrocarbons (PAHs) are benzene rings of sp2- hybridized C-atoms linked to each other in a plane, with H atoms or other functional groups saturating the outer bonds of peripheral C-atoms. Electron delocalization over their entire carbon skeleton makes them exceptionally stable because photonic excitation energy quickly becomes dispersed over the entire molecule before it can break a single bond.
According to the most accepted hypothesis, the infrared bands are attributed to the vibrational relaxation of PAH molecules containing around 50 C-atoms on average, which are stochastically heated to high temperatures (~1000 K) by the absorption of individual high-energy UV photons [48,49]. In the low density environment of the interstellar medium, where collisional de-excitation is not possible, these UV-pumped, gas-phases PAH molecules undergo spontaneous de-excitation via infrared fluorescence which gives the AIB features.
While the hypothesis requires strong UV radiation field to excite the PAH molecules (PAHs absorb strongly only in the UV), Uchida et al. [50,51] observed UIE features in the colder radiation environments of reflection nebulae, implying that other types of molecules may be responsible at least for some part of the emissions.
In support of this supposition is the fact that the PAH model is unable to explain the appearance of some aliphatic emission features at 3.4, 6.9 and 7.2 μm and the superimposition of the AIBs upon broad and strong emission plateau features at 6–9, 10–15, and 15-20 μm, which were shown to be of aliphatic nature [52,53].
Kwok and Zhang [30,31] made spectral decomposition analyses of archival spectroscopic observations of the UIE bands and showed that the data are most consistent not with pure, free-flying PAH molecules, but with solid organic nanoparticles that have a mixed aromatic– aliphatic (sp2-sp3) structure. They called this new model MAON (Mixed Aromatic/Aliphatic Organic Nanoparticle). The proposed MAON emission carriers are quite different from the pure hydrocarbon, planar, gas-phase, relatively small molecular (~50 Catoms) PAHs, where the structure, no matter how large, is regular with repeatable patterns. MAONs are amorphous solids of hundreds or thousands of C-atoms and impurity elements such as O, N, and S, with a disorganized, three-dimensional macromolecular structure comprised of aromatic and aliphatic units, each with variable sizes and random orientations.
Since the late 1980’s Papoular et al. have argued that the observed spectral properties of the UIE bands most closely resemble those of coal and kerogen [32,54,55], which are actually organic, amorphous, macromolecular materials with randomly oriented clusters of aromatic rings linked by long aliphatic chains, with O, N, S functional groups and heterocycles in their structure. This complex macromolecular structure is also very similar to the structure of the Insoluble Organic Matter (IOM) found in the Murchison meteorite and identified in laboratory analyses by Derenne and Robert . More than 70% of the organic matter in carbonaceous chondrites, such as the Murchison, Murray, Tagish Lake, and Orgueil meteorites is in this form of insoluble amorphous solids [57,58]. The structure of the IOM in all these meteorites shows remarkable resemblance , while isotopic ratios point to their interstellar origin [60,61].
It is interesting that in the UV-poor environment of protoplanetary nebulae (PPN), the aliphatic chains of these UIE emitters constitute a significantly larger percentage of the entire macromolecular structure when compared to the UV-intense regions of Planetary Nebulae (PN) where the aromatic component is predominant [53,62]. For example, in the PPN phase of stellar evolution, the 3.4 and 6.9 μm aliphatic emission features are as strong as the 3.3, 6.2, 7.7, and 11.3 μm aromatic emission features, but as the star evolves to the PN stage (a process of less than few thousand years), the aromatic features become considerably stronger . This phenomenon has been explained as the result of photochemistry, where the onset of UV radiation modifies the aliphatic side groups through cyclizations and isomerizations, and transforms them into aromatic ring structures, making the macromolecule more aromatic [63,64].
Much of the diverse organic material found on the surface of Earth is in the form of pigments that absorb strongly across the UV and visible wavelength region and dissipate this energy rapidly into heat (vibrational energy) which is then further dispersed over the vibrational modes of the surrounding solvent water molecules [65-68]. This is true of DNA and RNA which absorb and dissipate strongly from 220-280 nm, for the aromatic amino acids tyrosine, tryptophan, and phenylalanine which absorb strongly from 220-300 nm, for the mycosporine-like amino acids and scytonemin which absorb over 300-500 nm, and for the other common pigments of phototropic organisms; chlorophylls, phycobilins, carotenoids, anthocyanins, which absorb from the visible to the infrared.
It is an intriguing fact that most of the fundamental molecules of life (those found in all three domains of life) are also pigments which absorb in the UV-C and have chemical affinity to RNA and DNA. Through Förster-type energy transfer, as well as through the coupling of vibronic modes , these molecules can also use RNA and DNA as acceptor quencher molecules to dissipate extremely rapidly (through a conical intersection) the absorbed photon excitation energy into heat [67,68].
Planet Earth has for a long time been considered as the only body in the solar system that contains organic matter, usually attributed to the activities of life. The conventional picture for the chemical composition of other solar system bodies (planets, satellites, asteroids, comets, and dust particles) is that they are entirely made up of metals, minerals, ices and gases, with either no or only traces of organics. However, with the development of more sophisticated technologies like sending spacecrafts to solar system bodies for close spectroscopic observation and direct sample collection, organic substances are increasingly being recognized as a major component of planetary chemistry.
One of the first hints of the presence of organic matter on the surface of asteroids came from their deep red colors and low albedos (0.01-0.15). Such intense colors and low albedos are difficult to explain considering only inorganic minerals and ices but can easily be explained by organic pigments with a complex kerogen-like macromolecular composition . Roush and Cruikshank  have noted that terrestrial dark pigmented organic solids like tar sands, asphaltite, pyrobitumen, and kerite, which all have complex aromatic/ aliphatic composition, have very low albedos and dark reddish-brown colors, reminiscent of asteroids.
Saturn’s largest moon Titan, the only known natural satellite to have a dense atmosphere, has attracted much attention since the Cassini– Huygens unmanned spacecraft, which arrived there in 2004, revealed that Titan’s atmosphere contains haze-like solid particles that are the result of condensation of organics. The observations support the hypothesis that methane and nitrogen molecules excited by UV photons react to form polymeric hydrogenated carbon-nitride compounds, called tholins, which give the distinctive thick layer of orange-brown haze in Titan’s lower stratosphere [72,73]. Cassini RADAR observations found that these organic nanoparticles condense on the surface of sand grains that are blown into longitudinal dark colored dunes by winds. Most interestingly, some of these organic nanoparticles end up dissolved in the numerous lakes and rivers of liquid methane and ethane in Titan’s polar regions, which exhibit active liquid-gas phase cycling, similar to the water cycle on Earth, although at much lower temperature .
The Cassini spacecraft has also found evidence of mixed aromatic/ aliphatic hydrocarbons on other Saturnian satellites: Iapetus, Phoebe and Hyperion [75-77]. Iapetus is the third largest satellite of Saturn, locked in synchronous rotation about the planet, with the leading hemisphere and sides of dark brown color and low albedo (0.03–0.05), and most of the trailing hemisphere and poles of bright color and high albedo (0.5–0.6) . This difference in coloring between the two hemispheres is striking. Temperatures on the dark region's surface reach about 129 K at the equator, while the bright surfaces reach only about 113 K due to less sunlight absorption . The low-albedo material exhibits mainly an aromatic hydrocarbon signature .
Trans-Neptunian Objects (TNO) are a group of minor bodies that orbit the Sun in the outer solar system beyond the planet Neptune. Photometric observations in the visible have found the colors of some of these objects to be intensely red , which can be an indication of the presence of organic material on their surface.
Exceptionally rich in organic content and organic diversity are comets . They represent agglomerates of rocky debris, frozen gases and ices, and are likely the most primitive bodies in the solar system . Dust samples from the coma of comet Wild 2 were collected by the Stardust spacecraft and returned to Earth in 2006. The samples were analyzed with various techniques, giving evidence of rich organic content, but mostly in the form of aromatic-aliphatic macromolecules similar to the IOM of meteorites but with higher N and O content. Of particular significance is that the aromatic compounds are not pure hydrocarbons but have a very high N content, where N is incorporated predominantly in the form of aromatic nitriles (R–C≡N); a fact of great astrobiological implication as aromatic nitriles are precursory to many organic pigments and other fundamental biological molecules, and comets are likely to have contributed to Earth’s prebiotic chemical inventory [44,84].
In the following section, we will show how, under an imposed UV-C light flux from nearby stars, non-equilibrium thermodynamic principles based on the general evolutionary criterion of Classical Irreversible Thermodynamics, leading to increases in entropy production through UV-photon dissipation, would indicate the formation, proliferation, and aromatization of these organic compounds and thus explain their ubiquity in the cosmos.
A non-equilibrium thermodynamic explanation for organic pigment ubiquity
It is thus now abundantly clear that a large class of pigment molecules exist throughout the cosmos in the form of tholins, polyaromatic hydrocarbons (PAHs) and/or Mixed Aromatic/Aliphatic Organic Nanoparticles (MAONs) and that these systems are absorbing light in the UV and visible spectrums and dissipating it into infrared emission bands. The dissipative process in which a part of the incident solar spectrum is absorbed and dissipated into heat by organic pigment molecules at the Earth’s surface is by far the most important entropy producing process occurring on Earth. Similarly, the interstellar reddening due to photon dissipation by the interstellar PAHs and/or MAONs is probably the most important entropy producing process occurring in the cosmos after nuclear fusion and the expansion of the universe.
It will now be shown that the formation of these pigment molecules can be viewed as microscopic self-organization of material (microscopic dissipative structuring) in response to the imposed UV photon potential from nearby stars and that pigment proliferation to well beyond expected equilibrium concentrations follows from the same non-linear, non-equilibrium thermodynamic principles. Our thermodynamic analysis will follow that of Prigogine  for a set of purely autocatalytic chemical reactions, the difference being that instead of dissipating chemical potentials, our cosmic pigment system will be dissipating photochemical potentials related to the difference in the photon pressure of the relevant stellar spectra and that of the cosmic microwave background spectrum.
Advantages for microscopic dissipative structuring utilizing photochemical reactions over normal thermally induced chemical reactions occurring in the electronic ground state are various ;
1. Reactions may occur that are very endothermic in the ground state since the absorbed photon donates its energy to the molecule. For example, the energy of a photon of 260 nm (UV-C) is 4.77 eV or 110 Kcal/mol.
2. In the electronic excited state antibonding orbitals are occupied and this may allow reactions which are not possible for electronic reasons in the ground state.
3. Photochemical reactions can involve singlet and triplet states. Thermal reactions usually only involve singlet states. Therefore, photochemical reaction intermediates may be formed which are not accessible under thermal conditions.
The coupled photochemical dissipative process involving UV-C light dissipation and cosmic pigment production, which drives pigment proliferation, is the transformation of a black-body spectrum at a high temperature, that of the surface of a nearby star (e.g., our sun at 5760 K) into an intermediate black-body spectrum at a lower temperature (that of the vibrational spectrum of the low molecular weight precursor molecules of the pigments) after absorbing UV-C photons, and finally the dissipation to a black-body spectrum at a temperature of the infrared emission bands. Besides this photon dissipation process, we also have a “parasitic” photochemical reaction producing pigments from the excited precursor molecules. These pigments then act as catalysts for the dissipation of the incident stellar spectrum into the infrared emission spectrum, thereby making the global dissipation process autocatalytic and this brings the non-linearity into the nonequilibrium thermodynamic process. It is this nonlinearity that leads to pigment proliferation just as it does to catalyst proliferation in a purely autocatalytic chemical reaction process.
We will show that in the thermodynamic stationary state, the concentration of the pigment molecules can grow many orders of magnitude greater than what would be expected if the system were closer to equilibrium or if the pigment were not an effective catalyst in the photon dissipation. This is the same reasoning which explains the concentration of a product catalyst in a normal autocatalytic chemical reaction growing many orders of magnitude larger than its expected equilibrium concentration due to its involvement in the dissipation of the imposed chemical potential .
Photons carry momentum, and thus for a collection of photons in equilibrium (a black-body spectrum), a unique equilibrium pressure can be defined which is found to go as the fourth power of the temperature. It is assumed here that the rate of energy conversion between spectra will be linearly proportional to the difference of the photon pressures of the different spectra. This is similar to assuming that the rate of a chemical reaction is proportional to the difference in the concentrations of the reactant and product (with equilibrium and rate constants set equal to one), or that the flow of material is proportional to the gas pressure difference, or that the flow of electrical current is proportional to the difference in electrical potential. These latter two linear relationships for transport processes are valid when the mean free paths (between scattering events) of the particles are small with respect to the size of the system, i.e. that the system attains local equilibrium, while for the former chemical reactions (or photon dissipation), it is only required that the reactants and products (or photons) retain a Maxwell-Boltzmann (or Boltzmann) distribution of the velocities (energies) of the particles (or photons) involved.
There are then four relevant photon pressures, 1) Ps corresponding to the stellar photon spectrum arriving at the site of the precursor molecules, an approximate black-body spectrum with a temperature equal to that of the surface of the local star, 2) PI corresponding to the vibrational, assumed black-body, temperature of the precursor molecules after absorbing the energy of a UV-C photon and then electronic de-excitation through a conical intersection, converting electronic energy into vibrational energy, and 3) PE the photon spectrum of the infrared emission bands, 4) Pp corresponding to the spectrum emitted by the produced pigment which will be dependent on the local concentration of the pigments and the energy diverted into forming their covalent bond structures.
Although using a black-body assumption for these spectra is in error since the incident radiation from the star has particular directionality properties and only approximates a black-body spectrum, the purpose of the calculation here is to show qualitatively how the proliferation of organic pigments in the cosmos can be explained from a non-linear irreversible thermodynamic analysis of the photon dissipation process. A more accurate derivation of this would employ Planck's formula for the entropy of an arbitrary beam of photons.
Employing the formalism of classical irreversible thermodynamics , the entropy production is a product of generalized forces, X, times generalized flows, J,
P= diS /dt = Σ k JkXk ≥ 0 (1)
Where the sum is over all the irreversible processes occurring within the system. As an example, for a chemical reaction, the generalized force X is the affinity over the temperature A/T, where the affinity A for the reaction is equal to the stoichiometric coefficients of the reactants multiplied by the chemical potentials minus the same of the products. The chemical potentials, in turn, depend on the concentrations of the reactants and products and the affinity can be derived from the Gibb’s free energy differences of the reactant and product states. The generalized flow J for the irreversible process of a chemical reaction is the rate of the reaction.
For our photon dissipation irreversible process to be considered here, the generalized forces are the photon spectrum affinities determined from the photon pressures, and the generalized flows are the flows of energy which are transformed from the incoming short wavelength spectrum to the emitted long wavelength spectra.
The change of the entropy production dP can be decomposed into two parts, one related to the change of forces and the other to the change of flows ,
dP = dXP + dJP = Σk Jkdx k+Σk xkdJk (2)
Where the sums are over all irreversible processes occurring within the system.
In the whole domain of the validity of thermodynamics of irreversible processes, and under the condition of constant external constraints, the contribution of the time change of the forces to the entropy production is negative or zero,
dXP ≤ 0 (3)
This is known as the general (or universal) evolutionary criterion, or the Glansdorff-Prigogine criterion, since it was first established by Glansdorff and Prigogine  (see also Prigogine ). The general evolutionary criterion is the most general result that has thus far been obtained from Classical Irreversible Thermodynamics. For systems with constant external constraints (for example, a constant flux of UV photons from a neighboring star), the system will eventually come to a thermodynamic stationary state in which case ,
dXP = 0 (4)
With this background in classical irreversible thermodynamics, we can now analyze our particular irreversible processes consisting of the conversion of energy through different photon spectra. First, the incident flow of energy of the photon spectrum coming from the surface of the star IS(λ), which can be approximated by a black-body spectrum for convenience, IS(TS), where TS is the temperature of the surface of the star (say 5760 K, for a G-type star like our sun), is converted, by absorption on precursor molecules followed by either a photochemical reaction leading to a pigment molecule, or by the dissipation directly into heat giving an intermediate spectrum of intermediate temperature (also assumed to be black-body), II(TI), with say TI=1000 K (corresponding roughly to the vibrational temperature of a molecule of 50 atoms or so after absorbing a single UV-C photon). This spectrum is then converted through interaction of the molecule with its gas or solid MAON surroundings into the red-shifted photon spectrum of the infrared emission bands IE(TE), with say TE=287 K. Note that in the very low density environment of space, aggregation of individual molecules into large structures like MAONs would be thermodynamically selected since this would allow distribution of the incident photon energy over more microscopic degrees of freedom, or, in other words, emission at lower temperatures, giving greater entropy production. Small isolated molecules in a vacuum would have to fluoresce at high frequencies in order to dissipate the excitation energy of a single UV-C photon, leading to less entropy production.
This process can be equivalently looked at as one in which, at the intermediate temperatures of the vibrationally excited precursor molecules, a photochemical reaction may take place in which some of the free energy of the vibrationally excited molecules is diverted into the production of organic pigments leading to an emitted photon spectrum of the newly formed pigments which is dependent on the concentration C of the pigments and how much free energy went into forming the relevant covalent bonding and is also assumed to be characterized by a black-body spectrum of temperature Tp.
These newly formed pigments themselves act as catalysts for the overall dissipation process of Is(Ts) → IE(TE), and can be even more catalytic by, for example, providing local heat of dissipation to increase the probability of the photochemical reaction leading to a pigment II(TI) → Ip(Tp). Photochemical reactions are temperature dependent, just like normal chemical reactions. These two effects make the dissipation process auto-catalytic. A schematic diagram for this auto-catalytic photochemical process can be given as follows;
Since the system is far from equilibrium, the backward rate constants for the spectra conversions can be considered as being essentially zero, which is also the case for the photochemical reaction leading to the pigments (reaction 3). Note, however, that if a high energy photon can produce a pigment, a similar photon could also destroy it. There is escape from this predicament if the pigment is fortuitously produced with a conical intersection allowing rapid nonradiative dissipation, thereby neutralizing the destructive potential of the UV-C photons. Pigments produced without a conical intersection, but having chemical affinity to a molecule that has one, so that the pair could operate in the donor-quencher mode, would also be spared from the destructive potential of the UV-C photons. Such a selection process would tend to build up pigments with conical intersections or pigments with chemical affinity to those that have one, forming complexes like the MAON structures.
Since the photon spectra are all considered to be black-body, they are completely characterized by one variable which is the black-body temperature T, or equivalently by the photon pressure P, which goes as T4. Thus, the conversion of the energy through the different spectra and the energy flow into the formation of the pigments can be characterized, in a first approximation, by the pressures corresponding to the assumed black-body spectra. We can therefore rewrite equation (5) alternatively in terms of pressures;
Where, as mentioned above, PP is the photon "pressure" related to the black-body temperature of the pigments after their production which is related to their concentration or free energy that went into their formation. It is also assumed that the rate of production of the organic pigments (reaction 3) is linearly related to the photon pressure difference PI-PP. This pressure difference is essentially related to the free energy available for the formation of the pigment. If the photon pressure PP associated with the newly formed pigments is low, then more photon energy will flow into the production process, pigments will be produced at a greater rate.
In terms of the affinities and flows of the different processes, we can write Prigogine's general evolutionary criterion, Eq. (4) once arriving at the stationary state, in the following form,
Where the sum is over the three photon dissipative processes, A1 is the affinity for the conversion of energy of the solar photon spectrum into energy of the vibrationally excited intermediate precursor molecule photon spectrum, A2 is the affinity for the conversion of energy of the excited intermediate precursor molecule photon spectrum into energy of the emitted photon spectrum in the infrared, and A3 is the affinity for the photochemical reaction producing the organic pigment.
For the case of equilibrium photon distributions (black-body spectra), the affinities will go as the logarithm of the ratio of the photon pressures ,
The ᶹK in Eq. (7) are the rates of the corresponding energy conversion (dissipation) processes which we assume are related to the differences of the photon pressures (see discussion at the beginning of this section) attributed to the different black-body spectra.
Since the organic pigments are assumed to act as catalysts for the conversion of energy from the stellar spectrum to the emitted spectrum (i.e., the local heat of photon dissipation by the pigment can be dissipated into the environment (reaction 1+2) or can act as a catalyst for the photochemical production of a new pigment (reaction 3)), and since Pp is related to the concentration of the pigments, then we could model the autocatalytic process by multiplying the first dissipation process, Ps → PI, by a factor (1+αPp) where α represents the effectiveness of the organic pigment to act as a catalyst for photon energy conversion to heat (i.e., α→∞ for an excellent catalyst, and α→0 for a completely ineffective catalyst). Therefore, the rates of conversion, assuming all constants of proportionality equal to one for convenience (again, we are only interested in showing qualitatively the non-equilibrium dynamics of pigment proliferation), are given by
Note the non-linear relation between the forces, Eq. (8), and flows, Eq. (9).
Using Eq. (7) for the steady state together with Eqs. (8) and (9), taking the Boltzmann constant KB=1 for convenience, and observing that the free forces can be characterized in terms of the two free pressures, PI and PP (since PS and PE are fixed and given by the local star surface temperature to the fourth power and the infrared emission temperature to the fourth power respectively) gives:
At the steady state, we have (Eq. (4)):
Which, since , the overall affinity is a constant, giving, gives, as the only solution to (12):
Which implies that at the steady sate no free energy is being diverted into further pigment production v3=0. Any mechanism of pigment dispersal, however, would imply a continual production.
Also, at the steady state, solving Eqs. (10) and (11) gives,
Where we define is, therefore, a measure of the "distance" from equilibrium of the system). Therefore, since is much greater than (the pressures go as the temperature to the fourth power for black-body spectra) equation (14) indicates that the photon pressure related to the organic pigments , or in other words the pigment concentration, or the free energy which has gone into pigment production, has increased due to its catalytic activity in forming new pigments and dissipating the stellar photon spectrum into the emitted spectrum.
The entropy production of the energy conversion processes, including catalytic activity of the organic pigment, is given by flows times forces,
It can also be shown that the entropy production at the stationary state shifts to larger values as a result of the catalytic activity (Prigogine , for the corresponding case of chemical reactions).
These results give the non-linear irreversible thermodynamic explanation for the proliferation of organic pigments in the cosmos. Their concentration can become much greater than that expected under near equilibrium conditions due to their catalytic nature in dissipating the stellar spectrum, depending on the ratio of which is a measure of the distance of the system from equilibrium.
In referring to purely chemical reactions, Prigogine , in fact, noticed that such a result sheds light on the problem of the occurrence of complicated biological molecules in steady state concentrations which are of orders of magnitude larger than the equilibrium concentrations. In his 1967 book "Thermodynamics of Irreversible Processes"  he states, "Thus, for systems sufficiently far from equilibrium, kinetic factors (like catalytic activity) may compensate for thermodynamic improbability and thus lead to an enormous amplification of the steady state concentrations of the catalyst. Note that this is a strictly non-equilibrium effect. Near equilibrium, catalytic action would not be able to shift in an appreciable way the position of the steady state."
The cosmos is replete with organic pigments which absorb at the UV wavelengths of prominence of the local star and dissipate these high energy photons into the infrared. The low probability of formation of these pigments under the low temperature and high vacuum conditions of space suggest that they are the result of non-equilibrium thermodynamic processes involving photochemical reactions. In this article, we have shown how these pigments should be considered as microscopic dissipative structures formed as a non-equilibrium thermodynamic response to dissipate the impressed stellar photon potential.
The proliferation of the pigments to concentrations many orders of magnitude greater than what could be expected under near equilibrium conditions can be explained within this non-equilibrium thermodynamic framework by invoking the autocatalytic nature of the pigments in dissipating the same photon potential that produced them.
Evidence in favor of our analysis is many fold:
1. Organic pigment ubiquity in interstellar space would not be expected under near equilibrium conditions.
2. Aromatic pigments are found to be more prevalent under environmental light conditions of higher energy UV photons, whereas aliphatic chain pigments are found under conditions of lower energy UV photons. That is, there is a direct relation between the UV flux from a local star and the energy of maximum absorption of the produced pigment (aromatic pigments absorb higher energy photons than aliphatic pigments).
3. The aggregation of the organic pigment molecules into MAON nano-sized structures can be explained thermodynamically on the greater dissipation provided by these structures with many more microscopic degrees of freedom for distributing the vibrational energy as compared to the isolated PAHs.
4. Proliferation of pigments by the same mechanism on Earth has led to the proliferation of the fundamental molecules of life, almost all of which dissipate in the UV-C range, exactly over a UV-C window which existed in Earth’s atmosphere during the Archean [66,68].
It may be that in regions of the cosmos, dissipative self-organization of material under stellar photon fluxes has evolved to the point where the emitted spectrum is so far towards the infrared such that the emission lines have become many and too weak to be detectable. It may thus be possible that some of this dissipative structured material is contributing to the baryonic dark matter. Dark matter has been determined to make up some 27% of the mass-energy density of the universe and baryonic dark matter makes up approximately 10% of this dark matter density if the universe is at the critical density for closure, and if the standard model for big bang nucleosynthesis is correct. Since the known baryonic matter makes up less than 5% of the mass-energy density , the amount of material structuring dissipatively into the type of baryonic dark matter suggested here could exceed 50% of all observable matter in the universe. Other explanations for the invisible baryonic dark matter such as supermassive black holes or MACHOS have difficulty accounting for all of the missing mass.
If such self-organization into microscopic photon dissipative structures is indeed occurring in the cosmos, then it could be expected that the fundamental molecules of life (those found in all three domains of life), which appear to have been self-organized dissipative structures in the UV-C , should be ubiquitous throughout the cosmos. This would imply that the origin of life, as we know it [65,90] based on these molecules, should also be very common throughout the universe. Other forms of biotic-abiotic biospheres, such as the organic molecules floating in the clouds of Jupiter, absorbing sunlight and channeling the heat of dissipation into the multiple vortices and storms of the Jovian atmosphere, or the organic molecules in the methane lakes of Titan fomenting the methane rain cycle, could be equally, or even more, common.
Finally, we end this paper with an extraordinary quote by Chandra Wickramasinghe on work he performed with Fred Hoyle in the 1970's and 1980's, from his book entitled "A Journey with Fred Hoyle" ;
"A discovery of a 3.28 μm emission feature in the diffuse radiation emitted by the Galaxy confirmed that aromatic molecules of some kind were exceedingly common on a galactic scale. We argued that the infrared emission, not only at 3.28 μm, but over discrete set of wavelengths: 3.28, 6.2, 7.7, 8.6, 11.3 μm must arise from the ultraviolet absorption of starlight by the same molecular system that degrades this energy into the infrared. We had shown much earlier that the 2175 Å extinction of starlight may be due to biological organic molecules, and it seems natural then to connect the two phenomena. Thus we developed a unified theory for the infrared emission and ultraviolet extinction by the same ensemble of aromatic molecules."
Wickramasinghe and Hoyle, therefore, very early on, characterized very well the fact that the same aromatic organic molecules were responsible for the UV absorption and infrared degraded emission spectra. Here we hope we have been clear in bringing out the thermodynamic significance of this discovery; that the "exceedingly common on a galactic scale" nature of these aromatic molecules is directly related to their ability to "degrade this energy into the infrared".
K.M. would like to thank DGAPA-UNAM project number IN-102316 for financial support.