THE ROLE OF HYDROGEN CYANIDE IN CHEMICAL EVOLUTION by Rafael Navarro- Gonzalez Dissertation submitted to the Faculty of the Graduate School of The University of Maryland in partial fulfill~ent of the requirements for the degree of Doctor of Philosophy 1989 t . \ Advisory Committee: Professor Cyril Ponnamperuma, Chairman/Advisor Professor Norman J. Hansen Professor Joseph Sampugna Professor Raj K. Khanna I Professor Joseph Silve \r man Doctor Mitchell K. Hobish ? Copyright by Rafae l Navarro- Gonz.:1lez 1989 ABSTRACT Title of Dissertation: THE ROLE OF HYDROGEN CYANIDE IN CHEMICAL EVOLUTION Name of degree candidate: Rafael Navarro-Gonz~lez Degree and year: Doctor of Philosophy, 1989 Dissertation directed by: Cyril Ponnamperuma, Professor and Director Laboratory of Chemical Evolution, Department of Chemistry and Biochemistry Two maJor research areas are investigated: The electrosynthesis of hydrogen cyanide; and the role of cyanocomplexes in the free- radical oligomerization of hydrogen cyanide. The electric discharge production of hydrogen cyanide from a simulated primitive atmosphere composed of methane, nitrogen and water vapor was investigated. The radiation chemical yield {G) of formation of HCN was determined to be 0.26. A free radical mechanism was proposed to account for the observed chemical changes. Computer simulations of the reaction mechanism could effectively model the early stages of electrolysis of the gas mixture, and permitted the estimation of the rate of electrosynthesis of hydrogen cyanide under various atmospheric conditions . The possible role of cyanocomplexes of transition elements on the free- radical oligomerization of hydrogen cyanide was investigated. Aqueous, oxygen- free, dilute solutions of hydrogen cyanide and hexa- cyanoferrate(II) or (III) were submitted to various doses of gamma irradiation. The presence of either cyanocomplex led to a significant decrease in the rate of decomposition of hydrogen cyanide. The major products were ammonia and carbon dioxide . Computer simulations of these systems permitted the elucidation of the reaction mechanism and the derivation of rates of reactions of free- radicals with the cyanocomplexes. The results obtained provide an insight into the pos - sible role of cyanocomplexes of transition elements in chemical evolution. ii DEDICATION To my wife, Fabiola Aceves - Diaz, my son Rafael Navarro-Aceves, my mother, Carmen Gonz~lez de Navarro, and all members of our families. To the memory of my fath e r, Rafael Navarro- Contreras 1932- 1985 lll ACKNOWLEDGEMENT I am especially indebted to my advisor, Professor Cyril Ponnamperum~ for his constant support and encouragement throughout this study . He fostered creativity in my research by providing right degree of autonomy. I am thankful to Professor Joseph Silverman for a] lowing me have complete use of his ionizing radiation facilities, and to his graduate students Byron Lambert and Fuh - Wei Tang for their constant as - sisLance during the irradiation of samples. I appreciate the great interest and help of Dr. Yasuhiro Honda during the development of the dosimetric methods and preparation of samples for electric discharges. I am also very thankful to Dr . Mitchell K . Hobish for his great inte rest and help at severa] stages in this study. I am indebted to my Mexican advisor, Professor Alicia Negr6n- Mendoza, Chairman of the Chemistry Department of the Institute of Nuclear Sciences, U.N.A.M . , for her constant interest and support throughout this study; I am thankful to Professor lV Murcos Rosenbaum P., Director of this Institute, f o r his interest and advises during my academic training at Maryland. I conducted a portion of this work in their Institute , and am grateful for their consideration t o us e th e ionizing radiation sources and the laboratory fa c j] i1,jes durinc- my short stays at Mexico City. I want to acknowledge a fellowship from The National Autonomous University of Mexico (U.N.A.M.) during my graduate studies at Maryland. I wish to thank all members of the Laboratory of Chemical Evolution for providing a friendly and scien- tjfic atmosphere. Finu]ly to my wife, Fabiola Aceves - Diaz, and my sun Rafael Navarro- Aceves, for their constant encouragement and patient endurance, I am grateful in me asure beyond words. College Park, Md Rafael Navarro-Gonz~lez Spring, 1989 V TABLE OF CONTENTS Section List of Tables viii List of Figures ix Chapter 1 THE ABIOTIC CHEMISTRY OF HYDROGEN CYANIDE 1 1. 1 Introduction 1 1. 2 Abiotic synthesis of HCN 4 1. 3 Distribution of HCN 20 1. 4 Formation of biomolecules from HCN 26 1. 5 Research objectives 55 Chapter 2 EXPERIMENTAL PROCEDURES 58 2 .1 General 58 2.2 Preparation of electric discharge samples 59 2.3 Preparation of ll' - irradiated samples 62 2 . 4 Irradiation of samples 64 2 .5 Control Experiments 75 2.6 Analyses of samples 75 2.7 Radiation chemical yields: G 88 vi 2 . 8 Computer simulation of the chemical systems 89 2.9 Presentation of data 98 Chapter 3 ELECTROLYSIS OF A SIMULATED PRIMITIVE ATMOSPHERE. THE SYNTHESIS OF HYDROGEN CYANIDE 99 3.1 Introduction 99 3.2 Dosimetry of high voltage electric discharges 102 3.3 Electrolysis of gas mixture 122 3. 4 . Conclusions 152 Chapter 4 THE EFFECT OF CYANO COMPLEXES ON THE FREE- RADICAL OLIGOMERIZATION OF HYDROGEN CYANIDE 156 4. 1. Introduction 156 4 . 2 The ~-irradiation of aqueous hexacyanoferrate(II) 158 4 . 3 The b - irradiation of aqueous hexacyanoferrate(II)-HCN mixtures 174 4 . 4 The 'o' - irradiation of aqueous hexacyanoferrate(III)-BCN mixtures 196 4 . 5 The effect of cyanocomplexes in the free - radical oligomerization of HCN: Implications to chemical evolution 210 vii 4.6 Conclusions 217 Chapter 5 GENERAL CONCLUSIONS 219 5.1 Electrolysis of a simulated primitive atmosphere: The synthesis of hydrogen cyanide 221 5.2 The effect of cyanocomplexes on the free - radical oligomerization of hydrogen cyanide 224 REFERENCES 227 v i ii LIST OF TABLES Number Page 1.1. Estimated productions of HCN 21 1. 2 . Initial Radiation Chemical Yields 40 1. 3 . G 0 of formation of amino acids 46 1. 4 . Chemical Equations used in 53 1.5 . Percent Relative Rates of 55 2 . 1. Sensitivity factors for ions in 77 2 . 2 . Slopes (m) of calibration 87 2 . 3 . Radiation and molecular product 91 2 . 4 . Hydrogen atom reactions with 92 2 . 5 . Hydrated electron reactions with 93 2 . 6 . Hydroxyl radical reactions wi th 94 2 . 7 . Miscellaneous reactions with 94 2.8 . Relevant reactions in the fricke 96 3 . 1. Initial chemical yields {G 0 ) in 134 3 . 2 . Predicted radical and molecular 139 3.3 . Estimated productions of HCN in 151 4.1. Initial chemical yields {G 0 ) of 191 4 . 2. Initial chemical yields {G 0 ) of 206 4 . 3 . Initial radiation chemical yields 210 ix LIST OF FIGURES Numbe r 1.1. Formati on of DAMN : Initial Ste p s 27 1. 2 . Mechanism of Adenine Synthes i s 32 2. 1. Electric Discharge apparatus 60 2.2. Preparation of Oxygen-free 62 2. 3. Energy De pendance of the 67 2 . 4 . Schematic Representation of 68 2. 5. Experimental Layout for Measuring 71 2 . 6. Me asurement of Heat Generation by 72 2. 7 . Dose dependance of temperature in 74 2 . 8 . Calibration Curve for the 78 2. 9. Calibration Curve for the 78 2 . 10 . Computer Simulation of the Fricke 97 2 . C 1 o1 m. puter Simulations of the Effect 98 3 . 1. Time dependence of the 104 P 3 ea. k2 -. to- Peak and Zero- crossing 107 The Effect of Resistance 109 3 . 3 . Sch 3 . 4 . ematic Representation of 110 The Effect of Res 111 3 . 5 . istance of R:s on Dependence of the 3 V. 6 o. l tage (a) and 112 Pressure Depen 3 d. e7 n. ce of Voltage 113 3 T. e8 m. perature Increase as 1 a 1 6 Tem 3 p. e9 r. a ture Increase as 1 a 2 0 X 3. 10. Dose De p e ndance o f the 123 3.1 1 . Dose De p endance o f the 124 3. 12 . The Effect of Dose on the 125 3. 13. The Effect of Dose on the 127 3. 14. Gas Chromatogram of DNPH 128 3. 15. Mass Spectra of DNPH Derivatives 129 3. 16 . Dose De p e ndance of the Formation 131 3. 1'7. Ef f ect of Dos e on the Formation 132 3. 18 . Dose De pendance of the Chemical 133 3. 19 . Computer Simulations of the 145 3.20 . Computer Simulat ions of the 146 1. 1. Decomposition of 160 4 . 2 . Dose dependance of the 162 4 .3. Variation of pH with irradiati on 162 4 . 4 . Extinc tion coe ffi c ient of 163 4 . 5. Dose dependance of the 165 4 . 6 . Computed trends for the format ion 168 4 . 7 . Computed trends for the formation 170 4 .8. Dose dependance of the 173 4 . 9 . Variation of pH as a function of 175 4 . 10 . Decomposition of HCN during the 178 4 . 11. Formation of hexacyanoferrate 179 4 . 12 . Formation of molecular hydrogen 180 4 . 13 . Computer simulation for the 181 4 . 14 . Formation of CCb during the 182 4 . 15 . Formation of NH3 during the 182 ? Xl 4. 16. Gas chromatograms of DNPH 185 4. 17 . Electron impact (a) and chemical 186 4. 18. Formation of formaldehyde as a 187 4. 19 . Gas chromatograms of 188 4. 20. Dos e dependance of the formation 189 4. 21. El ectron impact and chemical 190 4.22. Dose d ependence for the formation 193 4 . ?.3. DPpendance of the decomposition 194 4.24. Ga s chromatograms of acids formed 195 4.25. Decomposi tion of HCN in solutions 197 4. 2G. Decomposition of hexacyanoferrate 199 4.27. Formation of molecular hydrogen 200 4.2D. Computer simulation fur 200 4 .20. Formation of carbon dioxide 201 4. 30. Effect of dose on the 202 4 . :H. G~s chromato~ram of DNPH 203 4.32. Electron impact and chemical 204 4. 3?.. Dose dependence for the formation 204 4.33. Gas chromatograms of acids methyl 205 4.34. Dose dependance of the formation 207 4. 3!J. Dose dependance of the formation 208 4.37. Dependance of the rate 214 4.38. Dependance of the rate 215 4.39 . Dependance of the rate 216 1 CHAPTER 1 CHEMISTRY OF HYDROGEN CY ANIDE THE ABIOTIC 1.1 Introduction scenarios of the origin of life on Earth Modern sed on ideas generated in dependently by Oparin are ba cording to these authors (1924) and Haldane (1929) . Ac pearance of the first liv ing cells was preceded the ap accumulation of organic c om- by abiotic syntheses and mplexity on the primitive Earth . pounds of increasing co nd These compounds were ess ential for the origin a lution of phase- separated systems that eventually evo ng to the emergence of the f irst populations of livi led 2; and Haldane, cells (Oparin, 1924, 1938 , 1957, 197 1938). g- The earliest published ex periment expressly desi of organic com- nated as demonstrating th e formation mosphere was pounds from a hypothetic al primitive at 2 that by Groth and Suess (1938). They exposed a C2 1 eV cm- 3 y- 1 can be derived from these values for the atmosphere at the surface of the Eart-h . An estimate of radiochemical prcx.iuction of HCN would therefore vary from 10-e to 10- 1 mole dm- 3 y- 1 bas e d on the fluctuation of G (HCN) determined under v~rious experimental conditions (0.01 to 0 . 9) by Zhdamirov et al. {1970) . 1.2.3 Electrolysis The mechanism of decomposition of methane is not well understocxl {Bossard et al., 1982). The major hydrocarbons {Ponnamperuma and Woeller, 1964; Ponnam- peruma and Pering, 1966; Ponnamperuma et al., 1969). Their relative yields depend on the pressure of the system and current of the electric discharge . At low 16 pressures ethane is the most abundant hydrocarbon (glow discharge ), but as the pressure is raised to atmos - pheric levels, a c etylene becomes the principal product (spark d ischarge ) (Fujio, 1930; We iner and Bu r ton , 1953 ; Sieck and Johnsen, 1963; Borisova and Eremin, 1968) . At highe r currents ( s ilent and corona d is- c harge s), ethan e is the most abundant hydrocarbon even at a~mosph e ri c pre ssures ( Li nd and G]ockl er , 1929 , 1930; Lind and Schltze, 1931; Ponnamperuma and Woeller, 1964; Ponnamperuma and Pering, 1966; Ponnamperuma et al . , 1969 ) . Th e fi rst report of electrosynthesis of HCN from a s imulated atmosphere of CH4 - Nz is that of Sanchez et al . ( 1966). Toupance et al. (1975) carried out a series of invest igations on the yield of formation of HCN as a function of N:z percentage in Cfk - N:z mixtures. They found that maximum production of HCN occurs in mixtures composed of 70% N:z and 30% Cfk (Toupance et al., 1975). Mo lecular hydrogen formed in situ (Raulin et al . , 1982) or added prior to electrolysis (Stribling and Mi lle r, 1987) decreases the rate of formation of HCN . Chemical yields (G) of products formed in electric d i scharge experiments are usually not reported because 17 of the difficulty in determining the energy introduced by electric discharges (Navarro-Gonzalez et al., 1986). In an attempt to measure the yields of HCN from CH4 - N2 mixtures subjected to electric discharges, Stribling and Miller (1987) estimai.,ed the energy introduced by comparing the heai., evolved in the experiment with that. of a known power source. The ability to reproduce i.,beir measurements was d etermined to be within a facto r of i.,wo (Stribling and Mill er, 1987). However, the authors did not take into account the energy converted into chemical form, light or that not absorbed by the gas. Gases have a low collisional mass stopping power when they are irradiated with low energy electrons (Spinks and Woods, 1976), as in the case of electric discharges. Therefore, their method underestimates the amount of energy supplied by electric discharges . Electric discharges are generally considered to be very efficient in synthesizing HCN based on the initial amount, of CH.q converted into HCN (Miller and Orgel, 1974; Raulin et al., 1982). However, further work needs to be done in order to understand the mechanism and efficiency of HCN synthesis induced by electric discharges in the primitive atmosphere . 18 1.2.4 Pyrolysis '!'he thermal decomposition of methane has been extensively studied under a variety of conditions (Skinner and Ruehrwein, 1959; Kevorkian et al . , 1960; Kozlov and Knorre, 1963; Palmer and Hirt, 1963; Kondratiev, 1965; Harting et al., 1971). The rate of decomposition of methane is independent of pressure , reac~ion conditions (flow, static, or shock wave pyro- lysis), and presence of an inert gas; the energy of activation calculated under various conditions is about 423 kJ mo le- 1 in the temperature range of 900-2,000?K (Skinner and Ruehrwein, 1959; Kevorkian et al., 1960; Kozlov and Knorre, 1963; Palmer and Hirt, 1963; Kondratiev, 1965; Harting et al., 1971; Chen et al . , 197 5). The dissociation reaction leading to the decom - position of methane is shown in reaction 1.17. Dimerization of CEh? or H? leads to the formation of the initial products: C2Ho or lb, respectively. Ethane is, however, rapidly consumed during the pyro- lysis leading to the formation of ethylene, the The latter further reacts to secondary product. produce acetylene, the tertiary product, as the extent of the reaction is increased (Chen et al., 1975; Roscoe 19 and Thompson, 1985). Heat (1.17) CH.,. CH3 ? + H? tmosphere The shock wave pyrolysis of a CH4-N:z a s of 5% CH.,. was studied by Rao et al. (1967). Mix~ure % N~ in argon were shock- heated to temperatures and 5 ide is formed between 1500 and 6000?K. Hydrogen cyan ' tarting from 2500?K and i ncreasin~ in yield up to s t 20% at 5000?K (Rao et a l., 1967). The second abou r rate of formation for HCN was determined as orde 1 1 n:r mole- s - 1. k (CH,, + N:z) = 9. 32X10- 3 - 5. 021 , in d - mechanism of formation w as not studied but a pos The ne ible scheme consistent wit h the mechanism of metha s fol - pyrolysis was suggested ( reactions 1.18 to 1.23, lowed by reactions 1.15 a nd 1.16). ( 1. 18) T22500?K N:z* N2 + H? ( 1. 19) CE-b: + N:z CH3 ? + N:z* + N:z + H? (1. 20) ? CH: CI--k: + N:z* ( 1. 21) ? CHN2* ? CH: + N:z* ( 1. 22) 2 HCN ? CHN:z* + ? CH: ( 1. 23) + . N: HCN ? CHN:z* 20 The chemical yield calculated for HCN was G 3.2 (Bar -Nun and Shaviv, 1975). Thi s value is high compared to that obtained by radiolysis (section 1.2 .3) and makes shock waves the most efficient energy source for producing HCN under thes e conditions. Bar- Nun and Taul>cr ( 1972) hav P- estimated thnt shock waves in the primitive atmosphere we re efficiently produced by lightine and met,eoritic impact, and an al>undance of 4.18 J cm- 2 y - 1 has Leen estimated. A d ose rate of 2.6X1019 eV cm- 3 y - 1 can be derived from this value. An estimate of the shock wav e production of HCN in the primitiv e atmosphere is lCJ" 2 mo l e dm- 3 y- 1 taking into a c count a of G=3.2 for HCN. 1.3 Distribution of HCN HCN is ubiqu i tous in nature. It has b een detected in the interstellar me dium {Snyder and Buhl, 1971), comets (Ulich and Conklin, 1974), and in the atmos - pheres of Jupiter {Tokunaga et al., 1981), and Titan {Hanel et al . , 1981, and Owen, 1982). Small amounts of HCN occur on the contemporary Earth; this is bio- synthesized from amino acids by several photosynthetic and non - photosynthetic bacteria (Nazly et al., 1981), algae {Vennesland et al., 1981), plants (Nartey, 1981), 21 and arthropods (Duffey, 1981). The probaLl e o c curre nce of HCN ou the primitive Earth is supported by its facile formation under a vari e ty of conditions . Table 1.1 summarizes the est ima t e d productions of HCN in the primitive atmos - ph c~ r e indueed by different energy sources. Table 1. 1. Estimated productions of HCN in the primi -- t i v e atmosphere by diffeninL energy sources . Energy source G 0 (HCN) HCN production (mole dm- 3 y- 1 ) UV light negligible 1 negligible Ionizing radiation: 0.01-0.97 Upper atmosphere 10- b - 10- .q Lowe r atmosphere 10- 5 - 10- 1 Electric Discharges undetermined Shoc k waves 3. 23 Total 1. Raulin et al., 1982; 2. Zhdamirov et al., 1970; 3. Bar- Nun and Shaviv, 1975. HCN synthesized in the primitive atmosphere would 22 h a ve been effici e ntly transported into the oceans by p reci pitation processes wh e re abiotic syntheses of compl e x organi c mo l ecul es may occur (Oro and KimLa ll, 196 1 ) . Fee ly et al. (1986} have d e rive d a HCN rainout r a Le of (3 - 14)X1011 mo le y 1 taking into account a He nry ' s law constant of 4Xl03 , a IJCN gloLa l mixing ra t io of 10- 6 , and a globa lly rainfall r a t e of 3.3X1(~ 6 It is difficult to estimate the concentration of HCN or any other compound in the primitive oceans. Miller and 0rgel (1974} have attempted to estimate the p oss ible uppe r limit of organic compounds dissolved in t h e prim i tive hydrosphere. They rationalized that jt i s v e ry unlikely that all of the surface carbon was eve r present as organic compounds, but 1 to 20% might have been. If the volume of the primitive oceans was as large as the present ones, this would have given a solution of 0.01 to 0 . 10 mole dm- 3 for a one- carbon- atom organic compound such as HCN (modified after Mi ller and 0rgel, 1974). This approximation places only an upper limit since not included in the calcula- tion is the loss of HCN by polymerization and/or hydro- lysis in the oceans. Such a high steady- state concen- tration does not seem likely since hydrolysis to formic a c id r e quires at most a short period of time ( 1()3 to ---- - 23 104 years) at moderate pH and temperature (Sanchez et al., 1966). The simplest mechanism that might have operate d on Ui e primitive Earth for concentrating organic compounds is evaporation {Miller and Orge l, 1974). Such a mechanism is unlikely to have contributed for the ac- cumulation of a volatile compound such as HCN. Sanchez et al. {1966) have suggested the eutectic freezing of HCN - water mixtures as a possible mechanism to increase the concentration of HCN. The eutectic in a HCN- water system occurs at -23.4?C and contains 74 . 5% w of HCN {Sanchez et al., 1966). This mechanism possibly ope rated on confined areas on the Earth such as polar regions, but was not dominant on a global perspective for the Earth since many of the physical conditions that prevailed in the early Earth were probably not significantly different from those that exist today (Ferris and Hagan, 1984); in particular, oceans and other bodies of water were subject to similar buffering processes, resulting in a pH that presumably was close to that of modern oceans (pH 8.0-8.5), or at least not far from neutrality (Ferris and Hagan, 1984). Navarro-Gonzalez (1983) has suggested the forma- tion of cyanocomplexes as a mechanism to protect HCN 24 from hydrolysis and to accumulate it in a nonvolatil e ?arm. Cyanocomplexes are readily formed when a cyanide s o lut,jon is in cont-act with rocks. Beck {1978) has shown tb aL when 1 R of powdered rock is suspended in 10 ml of 0.014 mo ] e d~ 3 aqueous solution of KCN, a varieLy of cyanocom11]exes of transition elements suc h ai; Mn , Fe, Co, Cu, and Mo are form ed . Equilibrium conditions are reac h ed within 36 hr, and the steady :::: t, a t e concentration of cyano complexe s varies from 10- 5 to 10- 4 mol e dm- 3 depending on the type of roc k; bexa- cyanoferrate ( II) and/or (III) were in most cases the principa l cyanocomp]ex formed (Beck, 1978). Kobayashi and Ponnamperuma (1985b) found that the yic Jd of HCN in the electric discharge of C[l3 to 104 years) at moderate pH and temperature (Sanchez et al ., 1966 )_ The source of transition metals for the pr imordial hydrosphere was the lithosphere (Beck, 197 8). Among the different transition elements present in the lithosphere and hydrosphere iron is by far the most abundant (Weast et al., 1985). The concentrati on of ely iron in the contemporary hydrosphere is ex trem small, e.g., 3.5X10- 0 mole dm- 3 for sea water (Brewer, 1975) and l.2Xl0-0 mole dm- 3 for river water l concentration 1S a (Mackenzie, 1975). This sma l consequence of the low solubility of iron( III) and the absence of complex- forming ligands. hydroxide, During the early stages of the primitive h ydrosphere, iron was probably present in its more solu ble form, iron(II). Hydrogen cyanide may have played an important role in increasing the dissolutio n of iron and other metal ions from the lithosphere (Navarro- onzlez et al., 1989). Cyanocomplexes of ir on are G 4 characterized by their high solubilities in water, and ole dm-3 4 Maximum solubility varies from 0.06 to 0.22 m for potassium or sodium bexacyanoferrate(I I) at 0 and l00C, respectively (Weast, et al., 1985). 26 h i gh s t ability constantse . 1.4 Formation of biomolecules from BCN I t is gene rally accepted that HCN synthesized in the prim i tive atmosphere would have bee n efficiently trans porte d into the oceans where it played a signifi - can t r o l e in the process that preceded the emergence o f life {Ore and Kimball, 1961; Sanchez et al . , 1966; Ferris and Hagan, 1984). A number of studies has been devoted not only to the identification of prcx:lucts r esulting from the oligomerization of HCN but also to t he elucidation of the mechanisms for their formation . This sect ion reviews the current status o f the studies of the ol i g omerization of HCN in aqueous solution initiated by ionic (section 1.4.1) or radical processes (section 1 . 4 . 2) . It also examines two possibly c ompeting processes in the primitive oceans: Hydrolysis vs oligomerization of HCN (section 1.4.3) . e The stability constant varies from 2.5Xl of rea c tion ( 1. 29) or llU? wj Lh HCN w :i nr; com?e ?Ljtion kine ?Lics 1n i r-rHd j H Led solu L j on:?: flCN -1 110? (1 . 29) 1 mq)1 1 :..; ' for r-Pac :Lion 1. 2~' using pulse radiolysis coup 1 e el tu fast, k inr~L i r.::; i; p r?c troscopy, and optimizing Tb c struc i, ure o f UH ' radi cal produced by this reaction wa: ; j nfcrrcd from i -t1(' s Lull I I products ( Drc:1Ganic et a 1 . , 1873) and also from specLro:;copy (Buchler et al., 197G; Several studies haVf ' lH!en don e t ,o determine the n1te of reaction (1. 30 ) o f 110? with - cN usine steady- state radiolysis (Kralji c and Trumbore, 1965) and pulse r1111 the hydroxyl radical dissociates (pK- 11 . 9) producincr tbe hydroxyl radical anion (?O- ) (Dragan i c and Drag ani c , 1971). This radica] readily r eact s witb cyanide ion (reac tion 1.31), k =2.6X10'3 dm~ mo ] e- 1 s - 1 : - CN -1 ? o- (1.31) Tb ,1 trans i e nt radical produced by this reaction has h P-en dFd,ected :in pulse radiolysis studies coupled to fa s t spectros c opy (Behar, 1974; Buchler et al., 1976). The short- lived intermediates produced in reac- tion ~ 1 . 26 t u 1.31 react among themselves and/or with HCN/CN- leading to the decomposition of HCN and to the formation of a large variety of compounds. Table 1.2 summarizes the initial radiation chemical yields (GD} 01 decomposition of HCN at various pH's. GD ( -HCN) is maximum at pH 6 (Draganic and Draganic, 1980). The variation of GD ( - HCN) with pH is primarily due to the 40 Tah l e 1.2. Init i al Radiation Chemical Yields (G 0 ) of Decomposi t,ion of HCN and - cN at d i ff e r ent pH 1 ? Compound pil Go HCN 2.4 5 . 8 6 . 0 13.8 NH..qCN 9.2 6 .5 NaCN (or KCN) 11. 3 5.2 1. Draganic and Draganic, 1980. diffe r e nces of rates of reactions of H?, HO?, and e-. q Cons iderable progress has been made in the charact,erization of some of the stable products produced in the free- radical oligomerization of dilut,e0 aque ous solutions of HCN (Draganic and Draganic, 1980) . About 79% of the nitrogen and carbon atoms from HCN mo l ecules decomposed are built into larger molecules, mainly various oligomers. A summary of the products that have been characterized is given below . Section 1.4 . 2.1 examines the mechanism of formation and 0 Unless otherwise stated, the data given is for 0.1 mole dm- 3 HCN aqueous solutions. 41 chemical y ie l ds 7 (G) of low mo l ecular weigh t, p roduct,s wh:i le t;ect, i on 1. 4. 2. 2 d eals wi tb large r mol ecul es and / or small mo l eL:ule~~ released upon thei r hyd r olysis. 1.4.2.1 Formation of low molecular weight prcxiucts in d ecomposi t ion of wa t e r mo l ecules {reaclion 1 . 25). Hydro1,;cn, G 0 =0. 4, a ccumulat,es dur i n ~ irradiation of t,be Fl CN so lution ( Draf; a ni c et al., 1973 ) . Hydroge n p e roxide , G 0 =0.3 for 0.006 mol e dm- 3 HCN, d isappe ars 1n re,H""! i , i ons wiU1 rad i olyt,i c produc t, s (Drartm1ic et al., 19'/3; Ogu n~, 1968 ) . .C.11r b cln_ d i92(:jdE;., G 0 =0.8, JS formed with a signi ficant yield, but G rapidly dec reases wiU1 dos e in the kGy? ranee to about, 0 . 5. It,s mechani s m o f forma - tion i s shown by reactions 1.32 and 1.33 . 7 Unl ess olb erwi se stat,ed, the chemical yi e lds are the initial value s (G 0 ): Those determined at the lowes t pos s ible doses. 0 Gray, Gy, i s the commonly used unit for radiation dose and its equivalent in SI units is l.8X10- 7 J mole- 1 ? 42 2Il0- CH -= N? llC( =NU)OH + HO- CN ( 1. 32) HO- CN + H20 C~ + NH~ ( 1. 33) Carbon dioxide accurnulat..es in solution due to its bir;h soluhi li ?Ly, and then re1,:H.:ts w.i.Lh hydrated e l ecL rons a s it., is indicated by r eaction 1.34 (Drat~m1ic et al ., 1973; Dragani c and Drar;anic, 1980 ) . C a nd 1. 3G. At h igh d oses o f irrad i a - tion, CH-. CllO, CII-. COCHO, and Cl:-b ( C0h-Clb are f ormed, anrl could o rig inate by short-ch a in olir:o rn eri:,:;at ion of the rad j c1-1 l ~~ produced by r eac-L ions 1. 2G t ,o 1 . ~i l ( Nc\v ar r o - Gonz~] e z , 1983}. Carh_gx~]_i~ Ac~-!Q!J ~ Fo rm ic ac i d, G 0 -=- 0. 1 f o r 0. ()()(") mo l e d m?- 3 HCN, may b e forme d by rea c t.ion;-; 1. 34 and J . 37 A large numbe r o f mon o - , di - , and tri - Crir boxyli c a c id s we r e identified in HCN (pH=6) and NB-,CN ( p f1 =9. 2) : a c onitic, adipic, 1,2,4- butene- tricar- boxy li c, butyri c , carboxysucc inic , citric, c itraconic, 1,2- dime thyl - succ ini c , fumaric, glutaric, itaconic, mal c ic , mal i c, ma lonic, me thyl - maloni c , 2 - methyl - tri - ca rl m llyli c, oxalie, pime lic, succinic, and tricarbal - lyl ic (N a varro - Gonzal ez , 1983; Negron - Me ndoza et al., 1 982, J 983 ) . Isolated, short- chain oligomerization mi ght, explain th e f o rmation of carboxylic acids or thei r amide s (Negr6n- Me ndoza et al . , 1983). Qre~ is an important product in the aque ous solu- tion s of a mmon i um c yanide , pH 9.2, (G =0 . 65, 180 kGy), 44 but is formed in low yield in aqueous solutions o f hydrogen cyanide, pH 6 (G =0.03, 130kGy). It may b e formed by reactions 1.39 to 1.42 (Dragani c and Draganic, 1980; Niketic et al., 1983). ? C ( OH ) =N- + fb 0 OH- + ? C (OH) =NH ( 1. 39) ? C(OH)=NH NfkC=O ( 1. 40) Nfk C=O ? Nfb ? NfkC( =O)Nf-k ( 1. 41) HO - CN + NH3 Nf-k C ( =O ) Nfb ( 1. 42) Carbamyl glycinamide, IbNC(=O)NHCfbC( =O)NEk, is formed in low yield in HCN at pH=6 (G =0 . 02, 130 KGy). It has not yet been detected in NfkCN solutions, pH 9.2. Carbamyl glycinonitrile, H:zNC( =O)NHCE-kCN, forms in low yield (G =0.05) in HCN solutions, pH 6 (130 KGy) and NI:-kCN solutions, pH 9.2 (180 KGy). These prcxlucts could originate by short- chain oligomerization of the radicals prcxluced by reactions 1.26 to 1.30, 1.38 and 1 . 42 (Niketic et al., 1983). Amino acids. Trace amounts of amino acids were found in some cases of irradiated solutions. Their origin was attributed to partial hydrolytic breakdown of oligomeric material (Draganic et al., 1976a). ~ - --------------------------- 45 Heteroc clic com ounds. A large number of products is formed; however, the identification of most of them still remains to be done. They represent 2 - 5% of the total mass o:f all non- volatile products. Adenine, uracil, thymine, and cytosine are formed. There is evidence of the presence of triazines, but the y have not yet been identified (Ponnamperuma, 1965; Negron - Mendoza and Draganic, 1984). Heterocyclic products are most likely fragments of short, side- chains of ol igomers which are formed as a resu lt of free radical processes {Negron- Mendoza and Dragan ic, 1984). D~MN is formed with a low radiation chemical yield (G 0 =10-- 3 ) in dilute aqueous solutions (Ogura and Kondo, 1967; Draganic and Draganic, 1980); however, as the concentration of HCN is increased {0.5 mole dm- 3 ), it becomes an important product (G 0 =0.45) in the irradia- tion {Ogura and Kondo, 1967). 1.4.2.2 Formation of larger molecules Between 50% (NfkCN, pH=9.2) and 79% (HCN, pH=6) of nitrogen and carbon atoms from decomposed starting 46 material are Luilt inLo Janre r rad:i olyti c prcxlucts. The oligorm?:rjc maLPrirc11i:anic nnd Dr?,11~anic:, lUU(: ). Ac id bydrolys is of Uw nonvolatile products re] eases several protu.i ni c und nonproteinic amin o acid~:; Th e compu~;j t, jun of amino acids in the Tab] e 1.3. G 0 of forroat,i()fl or amino H<: id~~l. - ---------?--? ------- - -- - ---- System Aru.in o a cid Il CN NH ... CN NaCN pll 6 pII =9 . 2 pH =ll.3 -- ---- ----------------- ---- -- Ala O.OO!i 0.004 0 . 022 Asp 0.0 12 0 . 010 0.008 Gly 0.14 0 0 . 600 0.240 Glu <0.001 <0 . 001 0.002 Hi s 0.002 Ser 0 . 002 0 . 005 0 . 018 Thr 0.018 <0.001 <0.001 Tota l 0. J f)() 0.620 0.290 1. Dragan:ic et al., 1976a; and Draganic and Draganic, 1977. 47 mixture does not depend strongly on the pH of system , but th~ir yi e lds do vary considerably depending on the initial pH of the solution (Table 1.3) (Draganic et al _, 1976a, 1977a; Draganic and Draganic, 1977). FracLionation of the nonvolatile material yield s basic (17% (HCN) and 38% (NfJ,,.CN) of total weight). acid ic 3% (HCN and NH.,.CN) of total weight) and neutral (69% (HCN) and 41% (NH4CN) of total weight) fractions. Th e analysis of the acidic fraction indi cates the l'I"f!'.;c?nce of dicarboxyl ic acids in the unhydrolyzed sample (Draganic et al., 1980; Niketic et al., 1983). Hydrolysis of these fractions releases several amino 1:wids in diffe rent percentages: 83.9% in the neutral fraction; 13.8% in the basic fraction; and 2.3% in the acid fraction (Draganic et al . , 1980). The oligomers from ammonium cyanide (pH=9.2) and hydrogen cyanide (pH=6) solutions are very similar (Dracranic et al., 1980; Niketic et al - , 1982, 1983). At present, the following oligomers have been detected: Oli_gg_tmin~, -(CH:z NH)"-, is the main constituent of the neutral fraction, and represents 16.4% of the total It releases trace mixture of radiolytic products. amounts of amino acids upon hydrolysis, but it is not 48 c ] e ave d by proteolytic enzymes (Draganic et al . , 1980; Niket j c et a I . , 1982) . lS present in the neutral fra c tion, and represents 8.2% of t.b e toi...al mass of 1.,be mixture of produc Ls . After hydro lys i s , glyc ine is released, constituting 55% of the mass of the oljgome r. The oligomer lacks peptidic bonds, and upon hydrolysis ammonia is released (Dra~cmic et al . ,, 1980). Other oligome rs witb urea- aldehyde fragments have p e ptidi c bonds and are isolated in different fractions because of various side- chain modificat.ions present in They represent t h ue p to oli g 4o 0%m e or fs . the non- The length of the v o pl ea pt -i -le radiolytic products. tidic fragments varies, but it is estimated to be up to 30% of the oligomer mass (M.W . 5,000) in some neutral fract,ious (Draganic and Draganic, 1980). Olj,gomers wi th~p_tid_ic bonds""? The presence o:f these oligome rs was established based upon :four tests: Amino acid content on acid hydrolysis a. (Draganic et al., 1976a, 1977a; Draganic and Draganic, 1977}. 49 b. Positive biuret reaction (Draganic et al., 1976a; Draganic and Draganic, 1977}. c. Infrared analysis (Draganic et al., 1977b) . d. Enzymal.ic c leavage wi l.b aminopepl.idas e M and pronase (Draganic et al., 1980; Nikctic et al., 1982). The oligomers are isolated from the bas ic fract.ion and rc"!present 5. 3% (NH"'CN) and 14-30% (HCN) of the toLal weight of the mixture of radiolyt.ic prcxlucts. Their moleeular weights are in the range 1,500- 20,000. Enzymati c treatment only cleaves 30% of the peptidic bonds (Draganic and Draganic, 1980; Niketic et al., 1982). Some of thes e oligomers show catalytic activity as demonstrated by the hydrolysis of p - nitrophenylace - tate ( pNPA ). The rate of hydrolysis9 of pNPA varies from 1. 0 to 7. 8 depending of the oligomer type, and is o:f the same order of magnitude as that of chymotryp- sin1 0 (Niketic, 1984b}. 9 Xl0- 8 mole (pNPA) s - 1 g- 1 (substance) . 10 k(pNPA ? chymotrypsin} =3. 7X10- a mole s - 1 g- 1 ? 50 In spite of the prog ress made in the elucidation of products from t,b e :free - rad i cal oligomerization of UCN, furLh e r work i s s1, i l J needed to chara cterize the comp]ex mjxt,urc of r adjolyLic products. 1 . 4 . 3 Hydrolysis vs oligomerization of HCN Sanchez et a 1 . ( 196 7) invest.igljt.ec.l Lhe extent of hydrolysis and i oni c o liaomerization of HCN in a basi c medium . They found that hydrolysis is the major reac- Lion channel at conce ntration s ~ 0.01 mole dm- 3 HCN . Ionic oligomerization i s important at concentrations 2 0. 1 mole dm- 3 HCN; t.h e ro 1.n: i mum conversion is centered at pH 9. 2, but, sharply f a ll c: at 8. 5 ~pH29. 8. Mixed kinetics occur in the range 0.01 - 0.1 mole dm- 3 HCN (Sanchez et al ., 1967) . Ogura et al . (1972) st,udied the fre e - radical olj - gomerization of HCN in the p1 range from 1 to 10, and in the concentration range from 0.0005 to 0.1 mole dm- 3 HCN. They found that the radiation chemical yield of d ecomposi-Lion of HCN is maxi mum (G 0 =6) in the pH range from 6 to 8 at concentrations of 0.006 mole dm- 3 HCN, and decreases to about half at 3~pH21 0. The authors also found that G 0 ( - HCN) i s not linearly r e lated with 5l t,b e initial concentra t , j o n of HCN, and varies from about, 3 -L o 9 . 5 i n t ,h e concentra1-ion range from 0. 0005 to 0. 1 mo l e drn 5 . J,'ormj c a c:i d, a hydrolyzed produc t of HCN, is f ormc~d 1t1 low yie]cl (Ogura, 1968 ); howe v e r, its o r ig?jn .i s by fr ee rad i(;al r e act,:i ons (1.29, 1.34, 1.37), m,d ri uL by b a~;r ?d o r a c :id cata ]y~ed r e a c tion ~-; (1.43). Th e l iaGf-: and /or a c id C<:} L1-1 lyzed hydro] ys is of HCN :is j n :..; ir~n ifir, ,:rnL und e r the author's condit,ions since the r- a t .e of prod u c Lion of radicals in very high, : 5X10- s HO- ~' HC( =O)Nib ( 1. 43 a ) HCN ILO ff' _,_ _ .;. HC ( =O) NH::.: ( 1 . 4 3 l,) Tht : ra L e of hydrolys is of HCN by base or acid catalysis i s low; howev e r, 1,he la1,ter :i s more si~riifi - c ant. The Arrhenius parame ters for the s e reactions h n v c b een r ccen1, ly repor1-ed by S-Lri bl ing and Mill er ( 1 :1 8 7): for k =k o exp( - E,,. R- 1 1 1 ), whe re, k ,.., = 4 . 8X101 1 dm~- mole- 1 s - 1 , and .Ec. = (85. 0) kJ mole- 1 , in the base catalyzed reaction (1.43a); and, ko= l. 8X1Cf' dm"3 mole- 1 s - 1 , and .Ec. = (90. 4) kJ mole- 1 , 1n the acid catalyzed reaction (1.43b). 52 With the data given in the previous and present sections, we are in a position to evaluate the possible extent of hydrolysis and oligomerization of HCN in the primitive hydrosphere. Such an estimate has not yet been available in the literature; however, the kinetic p a rameters for the major processes are known. There- f o r e , we can simplify such processes by a small set of simultaneous chemical equations {Table 1.4), and numerically integrate them with the use of a computer. The details of the computer program {ACUCHEM) are d e scribed detail in Chapter 2, section 2.7. Equations 1 and 2 in Table 1.4 are the equilibrium reactions for the dissociation of liquid water and HCN, respectively. The losses of HCN by base and acid hydrolysis are given by equations 3 and 4, successively. Equation 5 is the reaction leading to the ionic oligomerization of HCN. The production of radicals (R?) such as H-, HO-, and e..q- formed by interaction of various energy sources on the primitive hydrosphere 1s difficult to determine. If we only consider those produced. by radioactivity, and use a moderate energy flux of 20 J cm- 2 y- 1 (pages 14- 15, this study), the dose rate for the hydrosphere would be about 5Xl0-- 7 kGy s- 1 ? The total yield of radicals from 53 Chemical Equations used in Simulation. No. Equation Constant 1. J . Ib O fl()- + IJ-+- K=l. ox10- 1 "' 2 . IICN -- cN + H. .. K=6. ox10- 1 0 EhO 3 . HCN -I fm? -+ IJCON[b -+ on- k =6. OXl 0-"' ,. H::- 0 c:: 4 . HCN -I R? ? HCONib + w? r? k =2. 6Xl0- 1 1 rr,:?? r? 5 . BCN + -- cN Ion--oligomer k =2. 5X10- 7 :1 G. TL -0 R? k =5.2X10- 11 .=:. r 7 . HCN/CN-- + R? FR--oligomc r k = l. OX10'3 C , , , 8. 2R? R2 k = l.OX1010 r 1 . Equilibrium or rate constants are distinguished by t.b e symbol Kor k, respectively. The units are s -- 1 , and dm3 moJe-- 1 s-- 1 for first and second order rat.e constants, consecutively. The values are cRlculated for 25?C from equations given in this chapter. radiolysis of liquid water 1s G=6. 4; however, if we conside r a low production, G=2, a rate of radical formation of 5. 2x10-- 11 s -- 1 is derived (equation 6). This calculation underestimates G (R?) from radiolysis of liquid water, and does not include the cont.ribution of other energy sources, which may have contributed 54 significantly {page 35). The rates of reactions of these radicals with HCN/CN- are all known {page 36 to 39) , and we use a representative value of 108 dm3 mo le- 1 s - 1 {equation 7). Dimerization of this generic radical (equation 8) has to be taken into account as a like ly competition for equation 7. The rates of radical - radical reactions from transients of water are all known (Buxton et al., 1988), and are typically ( C around 1010 dm3 mole- 1 s - 1 I ; therefore. this value was I I I used in the computer simulation. I The pH used in the simulation is 8, and this value is generally accepted for the primitive hydrosphere (Fox and Dose, 1977). As pointed out on page 22, we are not yet able to estimate the concentration of HCN in the primitive hydrosphere. A concentration range of 0 . 01 - 0.1 mole dm- 3 was used as the upper limit in the computer simulation. Table 1 . 5 summarizes the net rates of hydrolysis and oligomerization via ions and free-radicals. The integration time was set for a relatively short period of time in geological terms, 2 months in order to avoid over-decomposition of HCN. As it is evident from these results, the decomposition of HCN is primarily via free - radical oligomerization reac- tions, 288%, at concentrations ~10- 3 mole dm- 3 ? Hydro- lysis is unimportant (~12%), but ionic oligomerization 55 Table 1.5 . Percent Relative Rates of Hydrolysis and Ol igomeriza~ion o f HCN for the Primitive Hydrosphe r e . HCN Hydrolysis Oligomerization Ioni c Free- Radical 1. ox10- 4 0 . 5 <0.1 99.5 1 . ox10-::s 11. 5 0.1 88 . 4 1. ox10- 2 54.1 2.8 43 . 1 1.ox10- 1 60.0 31. 5 8.5 is negligible {~0.1). Oligomerization via ions could b e come significant only at high concentrations {~ 0.1 mole dm- 3 ), but under this conditions hydrolysis also enhance s . 1.5 Research objectives Great progress has been achieved in the under- standing of mechanisms and efficiencies of HCN produc- tion triggered by ultraviolet light, ionizing radia- tion, and pyrolysis. The contribution of electric discharges to the synthesis of HCN in the primitive atmosphere is generally considered significant (Miller 56 and Orgel, 1974). Nonetheless, no accurate evaluation of such a contribution has been presented. The major obstacle has been the inability to measure the amount of energy transferred from the high-voltage, high- fre- quency electric discharges. Our research goals were centered in obtaining a better insight into two major aspects of role of hydrogen cyanide in chemical evolution . In particular, one of these was to accurately measure the chemical yield of formation of HCN in the electrolysis of a simulated primitive atmosphere and to elucidate its me chanism of synthesis. This work is described in chapter 3 . The other maJor research goal was to investigate the effect of cyanocomplexes of transition elements on the free- radical oligomerization of HCN. Specifically, we were interested in hexacyanoferrate(II) and (III), since they probably accumulated in the primitive oceans (section 1.3, this study). There is great in- terest involved in this investigation since only a limited number of studies in abiotic chemistry has been devoted to examine the effects of metal ions in catalyzing reactions of significance to the origins of life (Kobayashi and Ponnamperuma, 1985a). The limited 57 reports indicate that metal ions could significantly increase the yields of prebiologically and/or bio- logically important compounds from a variety of abiotic type reactions (Beck, 1978; Kobayashi and Ponnamperuma, 1985b). The systems HCN/Fe were ~-irradiated to trigger free - radical reactions in aqueous solutions, and various mixing ratios for iron were investigated ranging from 100.00 to 0.01%. The results and implications of this work are described in chapter 4. A general discussion and conclusions of all the data is given is chapter 5. 58 CHAPTER 2 EXPERIMENTAL PROCEDURES 2.1 General Al l chemicals were of the highest purity comme r - cial ly available, and were used without further purifi - cation un l ess otherwise stated . Their sources were: Ald r i c h Chem. Co . , Inc.; Allied Chemical; Eastman Kodak, Co; Air Products & Chemicals, Inc . ; EM Industries, Inc . ; Fisher Scientific Co.; J.T. Baker Chem . Co.; Mallinckrodt, Inc.; Matheson Gas Products; Sigma Chem . Co. Methanol and benzene were purified according to the procedures of Perrin ?t al . (1966). Deionized-distilled water was used for cleaning of glassware and for analysis of irradiated (electric discharge or ~-irradiation) samples. It was purified by flowing tap water through a Millipore system (Mil- lipore, Co.) consisting of charcoal activated and 59 anion- cation exchanger cylinders, and finally distilled into a glass container . Tr i ple distilled water was used for preparation and irradiation of samples . This was obtained by succ essive distillations from an alkaline { pff:, 10) solution of 0.01 mole dm- 3 potassium permanganate, an acid {p[?:F=2 ) solution of 0. 01 mole dm- 3 potassium dichromate, and finally distilled without reagents {Draganic and Draganic, 1971). In all cases, only f res hly prepared water was used. All glassware was made of borosilicate glass {"Pyrex"), and was cleaned according to standard procedures used 1n radiation chemistry {Swallow, 1960) . These consisted of immersing the glassware in concen- trated sulfuric acid and warm 7 mole dm- 3 nitric acid, followed by their abundant washing with tap water, deionized- distilled water, and triple distilled water, and finally dried overnight at 2200?C in an oven . 2.2 Preparation of electric discharge samples Three gaseous systems were investigated. A mix- ture of C~ - Nz - IkO was used to simulate an early primi - 60 tive atmosphere. This system has been extensively used in this laboratory to obtain basic information on abiotic synthesis of amino acids {Kobayashi and Ponnam- peruma, 1985b). The other systems consisted of gaseous CH4, and air {N:z - 0.Z); these were used to study the properties of electric discharges. C C r 2.2.1 Clk-N:z-EkO system rr 'r El e ctric discharge experiments of a simulated primitive atmosphere were performed using a 4.3 dm3 dumbbe ll shaped flask {Figure 2 . 1). 0 . 1 dm3 of 0. 05 mole dm- 3 ammonium chloride buffer at pH 8.7 was Figure 2.1. Electric Discharge apparatus. AC: Alter- nating current power; Cl: Coil; Cr: Condenser; E: Elec- trodes; G: Ground; H: Heating tape; HM: Heating mantle; I: Vacuum and gas inlet; MS: Magnetic stirrer; T: Transformer. 61 filtered through a Millipore membrane (pore size 0.22 u m, type GS) and then introduced into the flask. The syste m wa s air-evacuated by a diffusion vacuum pump. The gases dissolved in the ammonium chloride buffer solution were removed by 4 freeze - thaw cycles while the va c uum pressure was maintained at 30 mTorr. 200 Torr of me thane (Matheson Gas Products) and 200 Torr of nitrogen (Air Products & Chemicals, Inc.) of purity 299.99%, were introduced into the flask to yield a gaseous solution of 0.0108 mole dm- 3 of methane, 0.0108 mole dm- 3 of nitrogen, and 0.00782 mole dm- 3 of water vapor after allowing the system to equilibrate at 60?C prior to electrolysis . 2.2.2 e&.. and air (N:z/0:z=3.7) systems These systems were prepared in a 4.3 dm3 dumbbell shaped flask (Figure 2.1) or in a 2.1 dm3 round flask. Air taken directly from the atmosphere was introduced into the flask and the pressure was adjusted from 1 to 760 Torr with a pressure gauge. Methane (Matheson Gas Products) systems were prepared by removing the air from the flasks with a vacuum diffusion pump to about 50 mTorr, followed by filling with methane gas of purity 299.99%,, and its re-evacuation. The flasks 62 were filled with methane at 200 or 760 Torr. 2.3 Preparation of ~-irradiated samples Aqueous , oxygen- free, 0.1 mole dm- 3 HCN solutions containing {10- 3 - 10-e mole dm- 3 ) potassium hexacyano- ferrate{II) or {III) at pH~ 6 can be conveniently prepared using a simple setup {Figure 2 . 2). 0.3 dm3 of triple distilled water or an aqueous solution of (10-e _ 9 1 2 3 LI 5 Figure 2.2. Preparation of Oxygen- free Aqueous Solu- tions of HCN. 1. Argon cylinder; 2. Tank containing triple distilled water; 3. Reactor vessel containing aqueous sulfuric acid (1/1 v); 4. Syringe filled with saturated aqueous KCN or NaCN solution; 5. Syringe for injecting or removing sample solution; 6. Vessel con- taining HCN sample; 7. NaOH or KOH pellets; 8. Satura- ted NaOH or KOH solution to trap undissolved HCN; 9. To hooo extractor. 63 10- 3 mole dm- 3 ) potassium hexacyanoferrate(II) or (III) (Fisher Scientific Co. or Allied Chemical, res - pectively) is placed into a vessel (6 in figure 2.2), and the solution is degassed at room temperature by bubbling with 99.998% argon (Air Products & Chemicals, Inc .) for at least 30 min. HCN is generated in the reactor vessel (3) by slowly intrcxlucing 0.01 dm3 of a degassed saturated NaCN or KCN (98.5% purity from J.T. Baker Chem. Co.) solution (4) into 0.02 dm~ of degassed sulfuric acid (J.T. Baker Chem. Co.) solution (1/1, v/v). Gaseous HCN is intrcxluced into the sample flask (6), and the concentration of HCN is monitored by withdrawing aliquots with a syringe (5). Undissolved HCN is trapped in containers (8) with saturated KOH or NaOH. Aqueous, oxygen- free, 0.02 mole dm- 3 potassium hexacyanoferrate(II) (Fisher Scientific Co.) solutions (pH 6.7) were prepared using the same setup (Figure 2.2), except that HCN was not generated. 0.01 - 0.10 dm3 deoxygenated solutions of HCN and/or hexacyanoferrate{II) or {III) were prepared just prior to irradiation and were kept in special borosilicate glass syringes with caps during irradiation. When larger amounts of sample were needed for analysis, the 64 solution was introduced into special Pyrex vessels (0.45- 0 . 75 dm3). These vessels permitted gaseous radiolytic products to escape from the system prevent- ing a build- up of an over-pressure above the liquid phase, and precluding air from entering into the solu- tion. 2.4 Irradiation of samples 2.4.1 Electric discharge experiments The electric discharge across two electrcxies separated by a 0.1 dm gap was produced using high fre- quency generators (Tesla coils) from Electro-Technic Products, Inc. Three different Tesla coils were used: Two were model BD-50 (A and B) and one was model BD- 50E. The electrodes were made from 0.13 mm diameter tungsten rods from Laboratory Supply & Equipment Co. The intensity of the discharge was controlled by the output adjustment knob of the Tesla coil. In order to determine the amount of energy introduced by an electric discharge, we investigated the properties of the electric discharge at various spark gaps (0-0.1 dm) 65 and pressures {1-760 Torr) in three gas systems: air, CH,. - ~ - fk O, and CH. ... Three different methods were developed to deter- mine the amount of energy introduced into the systems by the electric discharge. The first was an indirect me thod, the che111ical dosi?eter {Navarro-Gan~ lez et a 1 . , 1986). This consisted of measuring the extent of decomposition of methane produced by the electric dis- charge and correlating it with a study where both the energy input and extent of decomposition of methane were measured. The other two were direct methods and consisted of measuring the voltage and current of the electric discharge with an oscilloscope and by calori - metry . The energy input per mole {E~) was calculated as E~ - IVtn- 1 ; where I and V are the discharge current and voltage, respectively, tis the reaction time, and n is the number of moles of the gas present. In order to ascertain the effectiveness of the calorimetric method reported by Stribling and Miller {1987), we also measured the energy input of the electric discharge according their procedure. Energy measurements using methane as chemical dosi- meter. In a chemical dosimeter the energy input is determined from the chemical change produced in a sub- 66 stra1-e. Calculation of the energy input requires a knowledge of the chemical yield (G) for the reaction or product, which is found by comparing the chemical sys - tern with some form of absolute dosimeter. Prior to 1986, absolute dosimetric methods of high voltage high frequency electric discharges were not available. Absolute dosimetry for low voltage electric discharges has long been recognized (Wiener and Burton, 1953), in which the energy input per mole (E~) was calculated from the discharge current and voltage. Wiener and Burton (1953) studied the electrolysis of methane at atmospheric pressure using a low voltage (~500 volts) direct current electric discharge. They measured the yields (G) of decomposition of methane and formation of acetylene, the principal product, as a function of energy input. Since this system is very simple to reproduce in the laboratory, we proposed its use as an indirect method to determine the energy in the high voltage electric discharge experiments (Navarro--Gon- zalez et al., 1986). Figure 2.3 shows graphically the raw data obtained by Wiener and Burton (1953). The extrapolation of G 0 was obtained using different type of polynomial regressions, and the best estimate was 67 2.0 ?-- ......, i 8 1.5 ....... .-.1.? t.O Q 3 1.0 2 .0 s.o 4 .0 0.6 b ~ 0 .. ( 0 ~ .......... D.'1 0 ,0-t--- -------~---i 0.0 0.6 1.0 1.6 2.0 Figure 2.3. Energy Dependance of the Decomposition of Me thane (a) and its Chemical Yield (b) in the Atmos - pheric Electric Discharge of Methane. Raw data was taken from Wiener and Burton, 1953. G 0 ( - CH~) =8.34, with an error of 20%. The power (P) of the high voltage electric discharges can be deter- mined by measuring the rate of disappearance of methane, k (-C~), and is calculated according to the following expression: 68 p - rges introduced (E1) b y electric discha is The energy - 1 ? computed by E1 = pt n ents with an oscillo scope. measurem fultage and curre nt ge and current of the volta The circuit for m easuring shown in figure 2. 4. is Vo ltage the electric disch arge R1 TC VM CM DM Rz AC =G c Cir- chematic Represen tation of Electri . Figure 2.4. S re nt ring the Discharge Voltage and Cur its for Measu urr~nt measuremen t; cu rrent power; CM: C AC: Alternating c u ital_mult1meter; G : DC: Discharge cham ber; DM: Dig coil; VM: Voltage : Tesla Ground; R: Resist ors; TC ~.R:z22.0 k Q; ~easurement. 20.0< Ri20.5 MD; 10.0 - lOOQ~}hL 50 D.. 69 ark discharge were ma de with a measurements of the sp tenfold attenuator Tektronix 515A oscill oscope using a lectricity supplied t o The e Probe (PG105 Tektron ix}. was Tesla coils was 120 V and 0.1-0.3 A, and the 020 B digital multim eter from John measured with a 8 resistors (R} used w ere of 10 F'luke MGL Co . , Inc. The %. W Power rating with a tol erance of 2-5 discharge were con-tric Resistor 1 and the e lec allel to the voltage source and the same nected in par ting path. The high conduc Voltage was applied to each and the output was attenu ated by resistor 1 Voltage of tor 2. sultant voltage was m easured across resis the re ies, the e resistors 1 and 2 w ere connected in ser Sinc e in this path was th e same. The current everywher was, output voltage of the Tesla coil (Vo} difference measured across here, Vm is the potentia l W sistance of R1 and 1k was varied resistor 2. The re .0 k fl, respective- from 0.5 to 20.0 M ~ . and 2.0 to 10 ly_ 70 t (I~) supplied by the Tesla coil The output curren arge and the path Was split between the spark disch nts in both containing R1 and Rz. The sum of curre is - l paths was equal to I ~. Since the spark d Paralle current harge was connected in series with Ih, the c e same at every point in this path. The voltage was th measured in order to d etermine the cur- across 1h was this path. The curre nt supplied by the rent through Tesla coil was, oss Vm' s the potential diff erence measured acr Where, 1 to . The resistance of R :s was varied from 50 resistor 3 1000 Q orimetry. fultage and current m easurements by cal ltage 2.5 shows the set up u sed to measure the vo Figure the electric discharge by calorimetry. and current of nitored by the The heat dissipated (Q1 ) from R1 was mo r1?metric box made of temp erat lo ure ri? se 1? n a ca temperature was meas ured with a Polyurethane. The nce thermometer from moctel 9540 digital pla tinum resista 71 T TC DM p AC C xperimental Layout fo r Measuring the Figure 2.5. E by Calorimetry. AC: Discharges Energy Input of Elec tric current power; C: Ca lorimetric box; DC: A~ternating multimeter; G: Ground ? Discharge Chamber; D M: Digital al therm~- hermometer probe; R: Resistors; T: Digit P : T 5 k D. eter; TC: Tesla Coil. 500~R1220 k D; R::z= m of h an accuracy of 0.00 3?C. The voltage Guildline wit ccording to the follo wing the discharge was ca lculated a equation: e was obtained from t he The current of the di scharg from both R1 and R::z, and was heat dissipated (Q 1 ? 2 ) following expression: Calculated according to the 72 0 I _ [ _Q_-_<_vz_Rz_l_R_i_) J ? ~ The res istance of R1 and Rz were varied from 20 to 500 k D ; and 5 k ~ . respectively. Energy determination by calorimetry. This method was suggested by Stribling and Miller {1986, 1987). It consists of measuring the rise of temperature in the ealorime t-ric box produced by the sparking of the electric discharge sample {Figure 2.6). The amount of TC - G p C Figure 2.6. Measurement of Heat Generation by the Electric Discharge according to Stribling and Miller {1987). C: Calorimetric box; DF: Discharge Flask; E: Electrode; G: Ground; P: Thermometer probe; T: Ther- mometer; TC: Tesla Coil. 73 energy introduced by this methcxl is not accurate since several aspects not taken into account are the energy converted into chemical form, or into light or not absorbed by the gas. Gases have a low collisional mass stopping power when they are irradiated with low energy electrons {Spinks and Woods, 1976), as in the case of electric discharges. Therefore, the methcxl of Miller and Stribling underestimates the amount of energy supplied by electric discharges; however, we also measured the energy of our experiments with this methcxl in order to compare its effectiveness. 2.4.2 ~-Irradiations experiments ~ - Irradiations were carried out using three radio- active 6 ?Co units located at the Institute of Nuclear Sciences, Universidad Nacional Autonoma de ~xico {U . N.A.M.) at Mexico City (dose rates of 0.7 and 20 kGy h- 1 ) and the Department of Chemical and Nuclear Engineering, University of Maryland (U.M.) at College Park (dose rate of 36.6 kGy b- 1 ). The irradiation dose was varied from 0.014 to 600 kGy. 74 The irradiations were initiated at room tempera- ture. The temperature of the samples did not increase during its irradiation at low dose rate (0.7 kGy h- 1 , U. N. A. M. ); however, it increased up to 45?C or 55?C when the samples received a total dose of 80 kGy using dose rates of 20 kGy h- 1 or 36.6 kGy h- 1 , respectively. Figure 2.7 shows the variation of temperature as a function of dose produced by the 0 ?Co source located at U. M. The increase in temperature during the irradia- tion of the aqueous solutions does not lead to a sig- nificant effect on the radiation chemical yield of decomposition of liquid water, G 0 (-EkO); G 0 (-Ek0) 60 0 0 ._,, Q) L :::, +' 40 0 L t> Q_ F t> f- 20 0 50 100 15C Dose ( kGy) Figure 2.7. Dose dependance of temperature in the ~- irradiation of aqueous dilute solutions. 0 ?Co source: Department of Chemical and Nuclear Engineering, Univer- sity of Maryland at College Park. 75 ly from 4.41 at 23?C to 4.54 at 65?c and varies slight .07% per degree (Dragan ic and Draganic, increases 0 rate 1971). The effect of temperatur e on reaction ts is also not large, in creasing only by a few constan nd Draganic, tenths of a percent per degree (Draganic a 1971). 2.5 Control experiments samples were prepared a nd analyzed simul - Control ey were taneously to irradiated samples except that th ubmitted to electrolysis or radiolysis. The not s absence of analysis of these sampl es corroborated the d the analysis contamination during the preparation an of the samples. 2.6 Analyses of samples eir The samples submitted to irradiation and th diately after irradiation controls were analyzed imme below. according to procedures described 76 2 6 Determination of gas products - -1 gas from the electric 3 An aliquot of 0 . 05 dm of harge ct ? arge sample was transfe rred from the disc lsch iner and flask to a high- vacuu m sealed gas conta on- a model 21-620A mass spectrometer from C analyzed by ics Co. The gases we re iden- solidated Electrodyna m porting ied by their fragmen tation pattern. Sup tif eir identity was obta ined by gas chroma- evidence of th reported previously ample, and was tography of the gas s gas obayashi and Ponnampe ruma (1985b). The by K were quantified by th e intensity of the Products the inant fragment ion wh en molecular ion or a p redom d by a dif- ormer overlapped wit h a fragment produce f of nt Table 2.1 summari zes the sensitivities fere gas. spectrometry. Ethyle ne th he mass e different ions in t -, antified because its molecular ion (C:z&.. Was not qu t (C2H3+, 27) overlap ped 2 ) and the predominan t fragmen 8 +, 28). ith those of ethane a nd also nitrogen (N2 W olution and Dissolved gases were extracted from s - a mass spectrometer (Consolidate Electro injected into s chromatograph dynamic s c o. ) 00 ga or 1? nt o a Varian 14 ivity detector (Varia n equipped with a therm al conduct ied Toepler pump acco rding to a Associates) by a mod if Torres et al., 1982; Procedure already de scribed ( 77 Tab le 2 . 1 . Sens itivity fa c tors fo r ions in the mass spectrome ~ ry ana lys is . Parent Molecule Frag me n t Se n s it i vity (MW) Factor1 CH.,._ Cfl.,. + (16) 2.1 C2Ek 40 - 6 o- + HCk - HO- + . Ck - 4. OX10'3 41 - 0- + EkCk EkO + -(k- ~5 . OX10'3 95 Table 2 . 7. Miscellaneous reactions with transients from wat er (Continue). No . Reaction k1 4 2 . o- + Cb . Cb- 3. 6X109 . o- EkO 4 3 + ? Q;_, - 44 . o- Cb + 2HO- 6. 0X108 + -(h- (k + (k2- 8. 0X108 45 -0:z- + w- Heb? 1.ox1011 46 . 0::z - + I-bO Heb? + HO- 1. 0X109 H::zO 47 . 0:z - + H 20,000 volt) at high frequencies (> 0.5 MHz}. Voltage measurements in spark discharges ~ere previously reported from this laboratory using a high voltage attenuator probe coupled to an In a different o as pc -illoscope (Park et al., 1975). Proach, the energy input was estimated by measuring the input voltage and current supplied to a Tesla coil (Shimoyama et al., 1978). The energy determined by this last method includes, however, the energy intro- 101 duced into the system not only as spark discharge but also the electricity converted into heat in the trans- former. A different method has been recently reported in which the energy introduced into the system was estima- ted by comparing the heat evolved during the sparking of an insulated discharge system with that produced by a known power input across a resistor (Stribling and Miller, 1986, 1987). However, the authors in that study did not take into account the energy converted into chemical form, light, or that not absorbed by the gaseous system. Navarro-Gonzalez et al. (1986) have calculated that up to 24% of the energy introduced by the electric discharge can be converted into chemical form. In addition, gases have a low collisional mass stopping power when they are irradiated with low energy electrons (Spinks and Woods, 1976) such as are produced by electric discharges. As a consequence the method of Stribling and Miller underestimates the amount of ener- gy intrcxluced by electric discharges. The purpose of this study is to develop dosimetric methcxls for high voltage, high frequency electric dis- determ ch inar eg e ths e raa dn iatio o na ately , d t ccur chemical yield (G) of HCN in the electrolysis of a 4.2 102 dm3 gaseous mixture at 60?C and 600 Torr composed of 0.0108 mole dm3 of methane, 0.0108 mole dm- 3 of nitroge n, and 0.00782 mole dm- 3 of water vapor in equi - librium with a liquid phase of 0.1 dm3 of 0.05 mole dm- 3 aqueous ammonium chloride buffer at pH 8.7. Three different methods were successfully applied to determine the dose supplied by electric discharges. The methodology is described in Chapter 2, section 2.4.1, and the results are reported here, in section 3.2 . Once an accurate knowledge of the dose introduced inLo the electric discharges experiments was concluded, we could measure the initial chemical yields (G 0 ) of HCN and other early products from the electrolysis of a CH,.q.-tb - fkO mixture. The results and a discussion are given in section 3.3.1. A mechanism to explain the electrolysis of the CH....-tb-EbO mixture, and the forma - tion of HCN is given in section 3.3.2. The implica- tions of this study to chemical evolution are given in section 3.3.3 and the conclusions in section 3.4. 3.2. Dosimetry of high voltage electric discharges The dosimetry of electric discharges was performed by chemical dosimetry, and electric and calorimetric 103 ts of discharge voltage and current. The measuremen easurements are given below . results from these m 3.2.1 Chemical dosimetry ethod The chemical dosimeter is an indirect m r and dose rate supplied developed to estimate the powe ctric discharges using the by alternating current ele chemical system investigat ed by Wiener and Burton ethane at 25?C and 760 Torr . (1953) : gaseous m alysis of the data reported by Wiener and Burton An ) on the direct current el ectric discharge of (1953 thane at atmospheric press ure indicated that the me ulated Power (P) of the electric discharge could be calc on 2.4.1, according to the following equation (see Secti Pages 65 - 68): 1.6022X10- 1 ' k(-Cfk) p - G O (-CU..) ction of methane in ~here, k(-CH ) is the rate of destru 4 104 0 e electric discharge; and G ( - CH..i.) is the initial th of methane. chem ica l yield of decomposition of Calculating P required an accu rate knowledge O O G - CH4} . G (-CH-i.) was not determined by Wiener and ( alculated from their Burton (1953), and therefore, w as c data for this study (see pages 65- 68, this disserta- O .,.) with tion}; a value of 8. 34 was deri ved for G (-CH an error of 20% (see Figure 2.3 ). ane was determined The rate of destruction of meth in the alternating current elec tric discharge of ) methane at atmospheric pressure (Figure 3.1); k(-C~ N N I 1.00 .0.. . X .---. 0 .76 l'l I) "a 0.60 ..&"..I a0 ...., 0.20 ..,. ::r:: 0.00 u 0 3000 600 0 9000 I Time (sec) ime dependence of the destructi~ n of Figure 3.1. T ane in the alternating current electric discharge meth of methane at 25?C and 760 Tor r. 105 was calculated from the slope of the curve shown figure 3 . 1 by linear regression analysis, and was determined to be : 8.73X1017 molecules sec- 1 with a correlation coefficient of 0.986. The standard error of the data was calculated to be within 6%. The power of the electric discharge was, therefore determined to be 1.68 W, with an estimated error of 20%. The error is mainly due to the deviation in G O ( - CH ... ) . The dose rate (Gy s- 1 ) of the electric discharge experiments was calculated from the equation, Dose rate P I m where, Pis the discharge power; and mis the number of kilograms of gas in the electric discharge system, and is 2.65Xl0- 3 kgr for our CH...-~-IkO system. The dose rate is, therefore, 634 Gy s - 1 , or 2.28 MGy hr- 1 ? This method was the first dosimetric technique that was developed. Initially, the values of power and dose rate of the electric discharges were considered as estimates since we did not know if these parameters depended on the composition of the gas phase. Subse- 106 quently, two additional methods were developed (Sec- tions 2.2.2 and 2.2.3), and it was found that these parameters do not depend on the chemical composition of the system. Furthermore, the dose rates determined independently by the different methods agreed to within a 15% deviation (Section 3.2.4). The chemical dosimetric method described in this section was simple and reproducible. The determination of power and dose rate was obtained in a matter of hours for a single assay: this includes sample prepara- tion, irradiation, and analysis. Several samples must be analyzed at various time intervals in order to obtain an accurate determination. The only precaution that must be taken into account is to measure the amount of methane decomposed when its rate of destruc- tion is linearly related with irradiation time; this is usually at methane conversions ~50% (Figure 3.1). 3.2.2 Electric measurements of voltage and current Asymmetrical sinusoidal signals with a frequency of 0.3 MHz were produced for both discharge voltage and current (Figure 3.2). This frequency was independent of the chemical composition and conditions of the 107 60 -> ~ 50 Ill .Isl l ~ 0 wIll j - 30 CJ !! A - 60 0 5 10 15 mo ....... 1.0 .: ........ b ... O.t'> f ~ 0.0 Ill ~"? -0.6 .cl () I? iS -1.0 0 6 10 16 190 Time (?sec) Figure 3.2. Peak- to- Peak and Zero-cro ta ssk ie nn g 1 Vro am lu et sh e Oscilloscope were plott c eu dr v be y f si mtt oi on thg to obtain the Voltage (a O ) s ac ni dll a Cti uo rrn es n tp r (o bd )u ced by an Electric Disc a ht a rgro eo m in te Am irp erature and atmospheric pr g ea sp s) u. r e (1 cm system. The discharge voltage and current collapsed in about 10 wsec, and re-generated in about 190 ?sec. A cluster of about 14 electric discharges were produced every 7.2 msec; the average frequency was calculated to 108 be 514 wsec per electric discharge. The total voltage and current {Y) supplied by the electric discharge were calculated according to the following expression, n A { ::::: X.,,.i ) i =l y - F whe r e , X.:C, is the absolute peak value of each positive and negative components during a transient discharge occurring in an interval of 10 wsec; A= 0.637, and is a conversion factor to calculate the area of a sinusoi - dal signal from its peak values{~) (Malmstadt et al., 1981); and F = 514 wsec, and is the frequency of elec- tric discharge. The total voltage was found to depend under cer- tain circumstances on the resistors (R) used to measure this quantity. Figure 3.3 shows the dependance of voltage on the total resistance of R1 and R:z. At low --------------------------------------------- 109 15()0 ,,...... -> ( Rt + R2) 1000 -QI ..I.I. f ~ 5()0 0 5 10 15 20 Resistance (MO) Figure 3.3. The Effect of Resistance of Ri + 1k on the Voltage of a 1 cm gap Electric Discharge in Air. total resistance for Ri and lk, the voltage increased up to about 3M ~ , where it reached a steady value of about 470 V (Figure 3.3). This effect was attributed to the facility by which electricity could flow more freely through the pathway containing Ri and Rz as the resistance on this pathway decreases {see Figure 3.4 for the electric circuit): Wben the resistance on such a pathway was zero, electric discharges were no longer produced. 110 TC CM VM DM AC =G Figure 3 . 4. Schematic Representation of Electric Cir- cui~s for Measuring the Discharge Voltage and Current. AC: Alternating current power; CM: Current me asureme D nC t: ; Discharge chamber; DM: Digital multimeter; G: Ground; R : Resistors; TC: Tesla coil; VM : Voltage measurement. 20.0~RiL0.5 M u ; 10.0~lkL2.0 k Q; 1000~R"32. 50 ~1. The voltage measurements for dosimetry of the electric discharges were taken at high resistance for R, and Rz, usually above 5 Mu , in order to minimiz e any interference on the electric discharge. The current of the electric discharge was not altered by the resistance used for lh. This is demons - trated in figure 3.5 where the current is plotted against a wide resistance range from 50 to 1000 u for 111 0 .009 ........ R3 ct ........ 0 .006 ...... 0 (.) ~ 0 0 '9 .. C.I.I 0 0 .. :: 0, .003 0 u 0.000 0 250 600 750 1000 Ri::iehilll!lnce ( 0) F igure 3.5 . Th e Effect of Resistance r oen ft ~ o of n a th1 e cm C ug ra -p Electric Discharge in Air. Electric discharge experiments were performed at a 1 c m dis c harge gap; however, in order to gain insight on -Lb e effe ct of discharge gap size on voltage and current, we studied the electric quantities {e.g., vol - tage , curre n t ) as a function of discharge gap size. The total voltage increased linearly with the length of the discharge gap. Figure 3 . 6a shows this effect from 0 to 1 cm discharge gap; at high discharge gaps {23 cm) , electric discharges were no longer produced. The total current of the electric discharge also increased linearly with the length of the discharge gap . This effect is shown in Figure 3.6b . 112 600 a ,..,. .>..., . .t,i 260 ~ 0 > ? ' 0.4 O.B 1.2 B.O ~.... b I>< ~ ......... 4 .0 .... A f ~ u 0.0 0.0 DA 0.8 1.2 Di?charge Gap (cm) Figure 3.6. Dependence of ( thb) e w Voit lh ta gth e e ( aE )l e ac nt dr i Cc uD rri es nch t p ao rgse ed g ao pf iA ni r a a st y sA tet mm o cs op mhe -ric Pressure. The pressure dependence of voltage and current was studied in different gaseous systems: air, CH4, and The electric discharges were produced with different Tesla coils and using different outputs. The experiments were performed at various discharge gaps from 0 . 1 to 1 cm. In order to normalize the data, the voltage and current per gap length were plotted 113 v ersus pressu r e (Figure 3 . 7) . The voltage per g ap length linearly i n c r e a sed with the tot.a l p ressu re of 1~ ,----------------, ?. ? , A X 10- 3 A. c:m- 1 b ._ X 10-? 10 t -? - ?- ?- - ? .- ?- ?-?- __. . _ -?-?-?- ?-?- ?- ?- ?- ?- ? ..... 6 ? a ? ? j - -i - -~ - - - - _. . _ - "- - -- - - -,-? 0 .___ __ ,__ __, ___ __ .__ __ _, o 200 .oo eoo eoo Pre111ru re ( 'l'orT' ) F igure 3. 7. Pressure Dependence of V C ou lr tr ae gn e t (a(L ) ) anp de r Discharge Gap Lengt c ho . i l N::z/0:z =m 3( o .A d 7) , e l TB eD s- l5 a0 , Output 8, disc e hx ac re gp et ga at p =7 06 .0 7 T co mr , r, where it varied { fr? o m) ; 0 .C 1l- k t, o D 1i .s 0c h ca mr ge Gap=0. 7 cm, T B eD s- la5 0 c, oo ilu t (p Bu )t m4 o. d5 e l( I! ); N::z/CH.,.=1, D Te is scla h arc go ei l g am po =d l el c mB ,D - 50E, output 5.6 co { i &l ) and{A) m o Td ee sl l aB D- 50, output 8 { L ). 114 the system and did not depend neither on the h . c em1ca1 c omposition of the system nor on the instrum t en used produce the electric discha tr og e (Figure 3 . 7a). The current per gap did not depend strongly on the total pressure of the system but did depend on b th 0 the Te s 1 a coil and the output used (Figure 3.7b). Tl ue minimum pressure plotted in Figure 3.7 corresponds to 1 T orr; however at lower pressures (~50 mTorr), electric dis - charges were no longer produced. The measurement of discharge currents and voltages was carried out in a matter of minutes and was Per- formed routinely for each experiment. During low vol - tage and current measurements a noisy baseline was obtained on the oscilloscope due to reception of radio waves (around 5 MHz) produced by Tesla coils. These radio waves were visible on the oscilloscope even with- out direct connection to the electric circuit. In order to eliminate the effect of noisy baselines, the determinations of voltage and current were obtained with a number of resistances for Resistors 1 to 3. R 1 was usually varied from 3 to 20 M Q ; 1k was changed from 2 to 10 k Q; and~ was modified from 0.05 to 1 k Q _ The standard deviation from all of these dif- ferent sets of measurements was determined to be 12%. 115 The power of the electric discharge was calculated by the equation, P =VI; where V and I are the total voltaae and current of the electric discharge, respec- tively. Its value depended on the discharge gap length and the total pressure of the system and ranged from 0.01 to 2.29 W. The power of the electric discharge did not depend on the chemical composition of the sys - t e m. From the above studies, we conclude that the power of the electric discharge delivered in the electrolysis of the CB.. . - N:z - fkO mixture was 2.29 W with a deviation of 0.27 Wat 60?C, 600- 160 Torr, and at 1 cm discharge gap length. The dose rate was calculated to be 864 G s - 1 or 3.11 MGy hr- 1 ? The standard deviation was deLermined to be within 12%. 3.2.3 Calorimetric measurements of voltage and current The values of voltage and current were derived from the measurement of the heat dissipated from resis - tors 1 and 2 in an insulator box during the sparking of the electric discharge experiment (see Figure 2.5 for the experimental diagram). Figure 3.8 shows a typical experiment in which the rise of temperature was 116 - 28 u 0 ........ Ql ~ -+-> Z6 ~ "Q"l s0- . II) E-< u 0 2000 4-000 6000 6000 Time ( sec ) Figure 3.8 . Temperature Increase as a Function of Time during the Calorimetric measurement of Voltage ( 0 ), I l' and Voltage and Current { 6 ) of the Electr ,,i c Discharge of Air at 760 Torr and 1 cm discharge gap. Ri=0.5 M Q ; lk =5 k D ; { 0) Qi=1.16X10- 4 ?C sec- 1 ; ( L. ) Q1 .... 2 =8.87X10- 4 ?C sec- 1 ; C:,=1594.4 W sec 0 C- 1 . measured as a function of time for voltage and both voltage and current of the electric discharge. The heat (Q) dissipated from Ri and R1+lk was calculated from the slopes (m) of the curves (Figure 3 . 8) and from the heat capacity (Cp) of the calorimeter: Q = Cp m. The beat released by R1 was unmeasurable at resistances of 21 M 2 , and only resistances of ~0.5 M 2 could be used to accurately determine the beat dis - sipated . The voltage of the electric discharge was 117 determined using resistors from 0.5 M Q to 0.02 M 2 without affecting its value. In the previous section, it wa s sh own that at low resistances for R1 and Ik, the voltage of the electric discharge is diminished {Figure 3.3), and consequently less current flows through the e l ectric discharge pathway and more through that containing R1 and lk. In spite of this alteration, the total amount of electricity that flows in both pathways was constant {see Figure 2.5 for the electric diagram). This conclusion was reached after determining that the t ota l consumption o-f e]ectricity by the Tesla coil was unaffected by connecting R1 and lk to the electric ,. discharge. Be cause of this I effect, the voltage and '., current we re simultaneously determined by measuring the heat dissipated by R1 and lk in order to quantify the total amount of electricity that flow in both pathways. Independently, the voltage was also measured to be abl e to solve for the value of current from the previous determination . The resistance of 1k was adjusted to 5 k Q and was not varied since the current of the elec- tric discharge was not interfered by the value used for Ik {see Figure 3.5). The power of the electric discharge experiments was determined to be 2.07 Wat 60?C, 600-760 Torr, and an 1 cm discharge gap length. It was found to indepen- 118 d ent of the chemical composition of the system use d. The dos e rate was d e termined to be 781 Gy s - 1 or 2.81 Th e standard d e v i ation was calculated to be 26% . Th e d os ime Lry of an electric di sct1arge experiment by t h i s method was simple and could be accomplished in a relatively short period of time, about 6 hrs. The only inconvenience was to determine the appropriate resistance of R1 in order to obtain a measurable beat dissipat ion rate; this was about ~0.5 M ~ at 1 cm dis - charge gap and 760 Torr of total pressure for the sys - ,? tem. :,I ' 3.2.4 Comparison of the different methods The three methods described above (sections 3.2.1 to 3.2.3) were effective in determining the dose rates of the electric discharge experiments. The dose rate vary from 2.28 MGy hr- 1 (using the chemical dosimeter) to 2 . 81 MGy hr- 1 (calorimetric measurements of V and I) or 3.11 MGy hr- 1 (electric measurements of V and I), and in all cases was independent of the chemical com- position of the system. The average dose rate was determined to be 2.73 MGy hr- 1 ? The standard deviation 119 among the three different methods was reasonable, 15% . This value is to b e contrasted to a 50% deviation ob- taine d by Stribling and Miller (1987) using a different approac.;h to determine the energy introduced by electric discharges. The method suggested by Stribling and Miller (1987 ) consists of determining the power of the elec- tric discharge by measuring the heat evolved during the s parking of the electric discharge sample (Stribling and Miller, 1987) . We have already discussed that this approach underestimates the amount of energy supplied by the electric discharge since it does not taken into account the energy converted into chemical form, light and that not absorbed (see section 1.2.3). In order to experimentally demonstrate the weakness of such a meth od, we also determined the power in our electric discharge experiments according to their method (see page 72 for experimental conditions). Figure 3.9 shows the temperature increase as a function of time during the sparking of an electric discharge experiment. The power (P) of the electric discharge was calculated from the slope of the curve (m) and the heat capacity of the calorimeter (Cp): P = Cp m. The average value of P was 0.37 W with a 120 26 26 24 ----~----,c-----------j 0 2000 4000 Figure 3.9. Temperature Increase as a Function of Time during the Measurement of Power of a 1 cm Electric Discharge in air at 760 Torr. m=4.89Xl0- 4 ?C sec- 1 ; C.C. =1594 . 4 w sec 0 c- 1 ? standard deviation of 10%. The dose rate was deter- mine d to be 140 Gy s - 1 or 0.50 MGy hr- 1 ? Therefore, the method of Stribling and Miller underestimates the dose rate by about a factor of 5. 3.3.4 Evaluation of the dosimetry of electric discharges Electric discharges have been extensively used to induce abiotic syntheses in the last forty years (see for instance, Fox and Dose, 1977). The measurement of the amount of energy supplied from the electric dis - 121 charges ha s b een usually neglected mainly bec aus e of the difficulty of d etermini ng such a parameter , and only r ecently h a v e there b een atte mpts to estimate this quani ,ity by Navarro- Gonz.:i lez et al . (1986 ) and Stri - bling and Mi ller (1986, 1987). Howeve r, this is the rirsL study 1n whi ch the dose rate of abioti c discharge expe riments is comprehensively investigated . Th e average dose rate of our electric discharge experiments was determined to b e 2.73 MGy hr- 1 at a 1 cm gap electric discharge and 600 Torr of pressure . This value depends on the discharge gap lengLh and the total pressure but does not depend on the chemical com- position of the system. Knowledge of the dose rate of electric discharge experiments is important because it allows to determine the dos e absorbed in the experiments and to calculate the chemical yields of the products. Furthermore, it p e rmits direct comparison with other energy sources in order to contrast the relative contributions that they may have played in the syntheses of organic compounds in the primitive Earth. The next sections on this chapter demonstrate the accomplishments that are ob- tained by knowing the dose rate of electric discharge experiments. 122 3.3 Electrolysis of gas mixture The electrolysis of a 4 . 2 dm3 gaseous mixture at 60?C and 600 Torr composed of 0.0108 mole dm3 of methane , 0.0108 mole dm- 3 of nitrogen, and 0.00782 mole dm- 3 of water vapor in equilibrium with a liquid phase of 0.1 dm3 of 0.5 mole dm- 3 aqueous ammonium chloride buffe r at pH 8.7, was investigated in the dose range from 0.38 to 152.30 MGy. 3.3.1 General Qecomposition of initial......ru!ses. Methane and nitrogen were identified and quantified in the electrolyzed samples and blanks by mass spectrometry. Figure 3.10 shows the decomposition of methane as a function of dose. The concentration1 of methane linearly decreased with dose up to about 20 MGy and then decreased steadily up to about 150 MGy; less than 6% of the 1n1- tial methane remained in the dose range from 80 to 150 1 The concentrations of solutes and/or prcxlucts are given in terms of moles per dm- 3 of gas phase. 123 - 0.012 "f' E '"O 0.008 ? 0 E -...., 0.00-4 ? ,......., -.t :r: ,_l_), ? 0.000 D -40 80 120 160 Dose ( MGy) Fig ure 3 . 10 . Dose De pendance of the Dec ompos ition of Me t h ane in the Hig h Vo ltage Electric Discharge o f CH4 - N2 - fk O . MGy. The initial che mical yield2 (G 0 ) of the d ecom- posit ion of methane was determined to be: G O ( - Cll4 ) =6 . 5 4 . Th e d ecomposition of nitrogen as a function of dos e i s shown in Figure 3.11. The concentration of nitrog e n linearly decreased with incre asing dose up to about 40 MGy, where it reached a steady value of about 0.0081 mole dm- 3 in the dose range from 40 to 150 MGy . Th e chemical yield of decomposition of nitrogen wa s 2 The chemical yield (G 0 ), number of molecules formed or destroyed per 100 eV absorbed, is calculated from the slope of a dose versus dose curve, e.g. as Figure 3 . 10. 124 r;") 0 .0 12 I 8 "O ..C...l). s0 ,......., N :7. 0 .006 +-- --~--- ~--~----a '--' 0 20 40 60 80 Dase ( MGy) Fi~ure 3 . 11 . Dose De pendance of th e Decomposition of Nitro{;en in the High Voltage Electric Discharge of en,. - N2 - fh. o. d e t e rmj ned t o b e: G O ( - ~) 1.26. Nitro~en was quite stable to electrolysis; its decompos ition occurred only in the dos e range where a significant amount of the original methane remained in the system down to about 20% of CH4. This suggests that the channel of decomposition of nitrogen depends on methane and/or a transient che mical species formed in the decomposition of methane during the bombardment of electrons into the system. Formation of hydrogen cyanide . HCN was identified and quantified by its titration with silver nitrate and 125 wi t h a s e lective cyanide electrode . Its concentration increased linearly at l ow dose and then steadily above 30 MGy (Figure 3 . 12); its chemical yi e ld of formation wa s d e termined to b e : G 0 (HCN) =0.26. About 10% of the initial me Lhane or nitro~t!O wa s converted into HCN at 150 KGy. t') 0 -' 1.2 X ........ 0 (") I 0.8 ~ Q) 0s 0.4 .__., ,....., zu 0.0 .:.:.r.:. . 0 60 100 150 DoBe ( UGy ) Figure 3.12. The Effect of Dose on the Synthesis of Hydrogen Cyanide in the Electric Discharge of CH,. - N:2 - fbO . In order to get a better insight into the chemical mechanism of formation of HCN during the electrolysis of the simulated atmosphere, other initial products 126 such as hydrocarbons and aldehydes were also - investi - gate d . 1''ormation of ___bygroca~bons. (b fb and (b Ho were identified and quantified in the experiments by mass spectrometry; however, (bfk could not be quantified because its molecular ion (C2H4+, 28) and the Pre- dominant fragment (C2H3+, 27) formed in the mass spec- trometry analysis overlapped with those of ethane and also nitrogen (ttz+, 28); the latter were extremely abundant. Figure 3.13 shows the formation of (bfb and ?tk as a function of dose. At low doses, their concentra- tions linearly increased with dose up to a maximum located at about 30 MGy; these maxima also coincided with the maximum decomposition of methane (Figure 3. 10). C2fb was the most abundant hydrocarbon, and its chemical yield of formation was determined to be: G D ( C2 H:z ) =2 . 14 . C:zHo was formed with a smaller chemi- cal yield: G 0 (C2Ho)= 0.57. The concentrations of C:zlk and C2Ho decreased rapidly with increasing dose start- ing after about 50 MGy. 127 ~ 2.0 2.0 1 I"; I'll .., I 1.6 n 1.6 I ~ ~ ? 1.0 'o 1.0 .. ,; ..f.l. , a 0.6 ? Oi 'i .l..)... 0.0 60 .. ~..... 0 100 160 ? Doee ( WGy} Figure 3 . 13 . The Effect of Dose on the Formation of C2H:z and C:zI:k, in the High Voltage Electric Discharge of CH,. - N2 - Ik O. Formation ol aldehydes. Figures 3.14 and 3.15 show a gas-chromatogram and mass spectra, respectively, of aldehyde products that were derivatized via 2,4- dini - trophenyl-hydrazine (DNP) into 2,4-dinitrophenyl- hydrazones (DNPH) derivatives. Acetaldehyde, propionaldehyde and butyraldehyde were identified by gas chromatography based on their retention times and coinjections with standard DNPH derivatives; additional evidence was provided by detection of the mass ions and the mass ions plus one in the electron impact and chemical ionization mass spectrometric analyses, respectively (see Figure 3.15). 128 1 Q) rtJ A 0 p, (I) Cl,) n:: 2 3 0 5 10 Time (min) Figure 3.14. Gas Chromatogram of DNPH Derivatives of Aldehydes Formed at 6.35 MGy. 1. Solvent (ClbCb); 2. DNP; 3. Formaldehyde; 4. Acetaldehyde; 5. Propional - dehyde; 6. butyraldehyde; 7 and 8. Unknowns. 4' and 6' are syn isomers. A peak (3) in the gas chromatographic analysis of DNPH derivatives at high doses, 21.59 MGy, (Figure 3.14) was resolved with a retention time corresponding to that of formaldehyde; coinjection of the electric discharge sample with a standard solution of formal- dehyde- DNPH resulted in an increment of the area under peak 3. The electron impact mass spectrum of this fraction (Figure 3.15a) demonstrated the presence of a low abundance, ~ 1%, mass fragment, 210 m/e, that could originate from formaldehyde-DNPH (M.W. 210); however, the mass ion plus one was not produced by chemical ionization mass spectrometry under similar conditions 129 100 122 5 0 152 1EIO 262 131 ~ lt'i w 0.. 0 50 100 150 200 250 300 Q) r.a ~ lil ..... b 225 0 100 ~ 50 122 l&i 206 0 60 100 160 200 250 300 m/e Figure 3.15. Mass Spectra of DNPH Derivatives Formed at 1.59 MGy . Electron impact mass spectrum (a). Mass ions: Acetaldehyde (224); propionaldehyde (238), and hutyraldehyde (252). Chemi ca l ionization mass spectrum (b). Mass ions plus one: Acetaldehyde (225), propion- aldehy< 9.0 ........ . t ai' .) X 10 I ~ 6.0 DX5 ..C...l.) 0 -s s.o ,....., 0 :u:c 0.0 ~ 0.0 . 2__ ., 0 4.0 6.0 Dose ( M:Gy) Figure 3.16. Dose Dependance of the Formation of Alde- hydes in the Electric Discharge of Cfk - ~-EkO: Acetal- dehyde ( 0 ); Propionaldehyde ( L ); and butyraldehyde ( D ) . Formation of oligomeric material. A water-insoluble material deposited on the walls of the flask during the elec trolysis of the gas sample. The material was analyzed for H, C, and N content in order to evaluate its yield of formation and to perform a material balance of major products from the electric discharge experiments. Figure 3.17 shows the formation of such a material as a function of dose. Its abundance linearly increased with dose; its rate of formation was calculated to be 2.12x10- 2 gr dm- 3 MGy 1 , and was 132 M I 0.06 El 'O a.. 1111 ....... 0.04 0 0 r...-..., <11 -~ cu 0.02 0 +J iU !:I "O :..c:i a .OD 0 tf.l 0 50 100 L60 '--' Dose ( MGy) Figure 3.17. Effect of Dose on the Formation of Solid Ma t e rial in the Electric Discharge of CH.:i-~-I:-kO . calculated from the slope of the curve given in Figure Figure 3 . 18 shows the chemical yields of this water - insoluble material in terms of carbon, hydrogen, nitrogen, oxygen, and chloride as calculated from eleme ntal analysis obtained at various doses. The chemical yields for this material were calculated from extrapolation to infinitesimal conversion and were 133 ,oo 0 1 ft 10- V G V ------,,.- - 2 D 10 ,a D - 3 - ... 10 ... J 0 50 100 150 Dose ( MGy) Figure 3.18. Dose Dependance of the Chemical Yield of Soljd Mater ia l in terms Carbon ( 0 ), Hydrogen ( ~ ) , Ni t rogen ( D ) , Oxygen ( v ) and Chloride ( i, ) . d e termined t o be: G 0 (H, C) =0.25; G 0 (N) =0 . 042; G 0 (0) =0 . 054 ; and G 0 (Cl) =0 . 0019 . Th e fun c tional groupings of the oligome ric mate - rial were not investigated. The empirical formula was c a l c ulated to be : C12eH12eCkaN:zCl . The main atomic constituents were carbon and hydrogen at one to one ratio . The primary interest in the quantitative analysis of this material was exclusively to obtain a balanc e for the prcxlucts as completely a s pos sible. 134 Material balance . Table 3.1 summarizes the initial chemical yi e lds {G O ) for the analyzed products in the electrjc discharge expe riments. Ac :et.ylene and ethylene were the most abundant products, followed by the oligo- me ric materia l , acetaldehyde, hydrogen cyanide, butyr- aldehyde, and finally propionaldehyde. In order to maini,a in a mat,erial b a lance, tbe chemica l yield(s) of decomposition of solute(s) must equal to the chemical yi e ld s of identified and unknown products. For the Tabl e 3.1. Initial chemical yields (GD) in the elec- trolys is of gaseous CH..,. - N2 - fbO. Compound GD - CH..,. 6.54 - N2 1. 26 C2fk 2. 14 C2Ho 0.57 HCN 0.26 CfhCHO 0.13 CH3 CfkCHO 0.011 CH3 ( Cfb ) z CHO 0.016 Solid Material: C 0 . 25 H 0.25 0 0 . 054 N 0.0042 Cl 0.0019 135 case of carbon, the following expression is derived: G O O O O ( - CH. . ) - 2G ( C2 Eh ) -t 2G ( C2 H.,,, ) + G ( HC N ) + 2G O ( CH:,s CHO ) + 3G O ( CH"' CI-b CHO ) + 4G 0 (Cl-b(Cfb) 2 CHO) + G 0 (Solid mater1aJ) + G 0 (Unknown products); where, 6. 54 :_::, 6. 29 + G O ( Unknown products); G O ( Unknown carbon containing products) ::c? 0. 25. The identified products represent approximately 96% of the initial carbon-contain ing products in the electrolysis of the CH.,. - N~ - Ik O mixture. For case of nitrogen, HCN and the oligomeric material were the only nitrogen-products identified, and the following expression is obtained: 2G O ( - N:z) = G O (HCN) -t G O ( Solid material) + G 0 (Unknown products); where, 2 . 52 = 0.26 + G 0 (Unknown products); therefore, G 0 (Unknown products)= 2.26 . These calculations indicate that about 96% of the total methane and 10% of the total nitrogen were con- verted into the identified products in the early stages of electrolysis of the gas mixture. A number of other 136 producLs are also formed in the experiments such as amino acids (Kobnyashi and Ponnamperuma, 1985); however , they are formed through secondary reacLions invo]ving the primary products, e.g., ~Eb, C2 fb , CH ~CHO, HCN, eLe. Th e resu]Ls presented in -Lhis sec-Lion ennl>le us to infer possible reaction schemes that may operate in the c lec-Lrolysi s of a CEL. - N2-fbO mixture. The next section deRJs with this topic . 3.3.2. Reaction mechanism A pos s i bl f! reaction mechard sm was developed to accour1L for the observed chemical changes described in the previous section. The reaction mechanism was tested by computer simulation using ACUCHEM (see Chap- ter 2, section 2.8) to evaluate its potential to ex- p] ain the early stages of electrolysis of a CI:-L. - N::e - EkO mixture. During the irradiation of the gas mixture, the electrons from the electric discharge interact with gas molecules producing excited gas molecules (Reactions 3.1 to 3.3) (Bossard et al., 1983). The probabi]ities 137 of interaction depend on the concentration 0 f each gas in the mixture, and are 37.5% for either methane and nitrogen, and 25% for water vapor (deno~ed by H~O) _ e CB,. CH,. ... ( 3 . l ) e NZ ~ .. (3 . 2) e EbO - ---t Ek 0-- ( 3. 3) Th e exc jte d molecules produced in reac tions 3 _ 1 to 3. 3 may undergo dissociatjon producine fref-! r1:1djc:a]s_ Exci Le d water molecules (EkO*) dissoc iate producing hydrogen atoms (H? } and hydroxyl radicals (HO? ) (Reac - tion 3 . 4) (Spinks and Woods, 1976; Maurey et al., 1901) . Excited methane could dissociate producin11 Cfb?, : Cfb and :CI1? radicals, reaction 3.5 (Slager and Blac k, 1982; Bossard et al . , 1983). Computer fittinrr of the experimental data suggests that channel 3.5a is the major dissociation reaction accounting for 87.5%; whereas channel 3.5b and c were less important. H2 0- H? + HO? (3 . 4) 138 CH ~. - + H. (87 . 5%) (3 . 5 a ) I cu .. - I : : Cfb + 2H? (11. ox) (3 . 5b) : CH? + 3H? (l.5%) {3 . 5 c ) Nitrogen mo l ecules we r e found t o be stable to e l cc: t .rolysis in a CII,. - N::z - f-h- 0 mixture , and only a small p e r centane of the nitrogen decompos ed in the early s t .a ges o f e lectrolysi s , reac tion 3.6a . Computer simu - lation suggests that channels 3.6a and b represent 0. 2 5% a nd 99.75%, respectively. ?N: + ?N: (0.25%) (3.6a) N::z + hv {99.75%) {3.6h) Compute r fitting of the reaction mechanism (reac - tions 3 . 1 to 3.30) with the experimental results of decomposition of CH,. and formation of C2H:z, C:zHo, HCN, and C::z H.q, indica-Le that the rate of production of elec- trons {e ) in the electric discharge was 7.95X10- 3 elec - tro n s MGy- 1 and its yield was calculated to be G 0 , (e) = 121. 5. Table 3.2 summarizes the predicted radical and mo lecular yields of primary species obtained by 139 Table 3.2. Predicted radical (R) and molecular (M) yield s (GR,,M} in the electrolysis of CH.. - N:z, - fko. Spec ies Gf"II./M Cfh? :Clk 2 6.42 : CH ? 3.32 H? 0.45 7 H 0O .? 76 . N : 36.34 0. - 2C 7H ... - N~ 30.20 0.14 -l kO 36.34 simulation of the electrolysis of CH ... -N2-EkO. Cfh?, H?, and HO? are the major radicals, and :Cfk, :CH? and ? N : are produced in low yield. Supporting evidence for the formation of CEbs?, :CEk, :CH?, Ek and H? is provided by their spectroscopic identification in the electric discharge of gaseous methane (Peters and Wag- ner, 1931; Harkins and Jackson, 1933; Willey, 1934). Water vapor was constantly produced during the ex- periments by maintaining the aqueous phase at 60?C. The vapor pressure of water was 149.38 Torr at 60?C (Weast et al., 1985). The equilibrium reaction of 140 liquid water {denoted by fkO (1)) going into water vapor is given by reaction 3.7; the equilibri?um con- stant was calculated to be l.3Xl0-- 4 , taking into ac- count the concentration of water vapor and liquid wate r. Ek O ( 1) fbO (3 . 7) An additional source of water vapor was its refor- mation from the reactions of hydroxyl radicals with hydrogen atoms {reactions 3.8 and 3.9). The rate con- stant for reaction 3.8 is 2.2X1016 1 2 dm6 mole-2 ? s-1 (Tsang and Hampson, 1986), while that for reaction 3.9 is 1. 5Xl019 1 2 6 ? dnf' mole- 2 s - 1 (Wilson, 1972). H? + HO? + tb fbO + tb (3.8) H? + HO? + fbO fbO + fbO (3.9) In addition to the formation of water molecules, hydroxyl radicals could dimerize producing hydrogen peroxide (reaction 3.10); the rate constant for such a reaction is 2.4X1011 (T/300)-0 7 ? dnf' mole-2 s - 1 (Baulch, et al., 1982). 141 HO? + HO? (3. 10) Molecular hydrogen was not examined in thi s work, but supporti ng evidence for its formation was provided by the studies of Wi ene r and Burton (1953), Toupance et al. (1 975), and Bossard et al. (1983) on the electro - lysis of methane and/or me thane - nitrogen mixtures . Formation of molecular hydrogen is predicted by di - me ri zat ion of hydrogen atoms, reactions 3.11 to 3.12; the rate constants are 5.4X1012 1 1 3 dm6 mo1~- 2 ? s- 1 ( Cohen and Westberg, 1983; Tsang and Hampson, 1986), and 1 . 0X10 1 3 1 1 dm6 mole- 2 s - 1 (Cohen and Westberg, 1983), respectively . H? + II? + N:z 1-k + N:z (3.11) H? + H? + 1-kO 1-k + E-bO (3.12) Methane molecules were re- formed in the electric discharge according to reaction 3.13; the rate constant for such a reaction is l.2X1012 1?? 4 dm3 mole- 1 s - 1 (Tsang and Hampson, 1986). + H? (3.13) Methyl radicals may also dimerize leading to the formation of ethane, reaction 3.14; the rate constant 142 f o r this reaction is 1.0X101 2 1?? 04 dm3 mol e- 1 s - 1 (Tsang and Hamps on, 1986). A minor los s of me thyl radi ca l s may b e attr i buted to its reaction with me thyl - ene radicals which could lead to the formation of ethylene, reaction 3 . 15; the rate constant of such a r eac tion i s 4 . 2X10 1 0 dm3 mol e- 1 s - 1 (Tsang and Hampson , 198 6) . Cih ? + c~- (3.14) CH-.., ? + :CE-k (3.15) Acetyl e ne may be formed by collision of two CEb radicals, reaction 3.16 (Weiner and Burton, 1953); this r eaction bas a rate constant of 3. 2X1010 dm3 mole- 1 s - 1 (Tsang and Hampson, 1986). Acetylene may also be forme d by dimerization of two CH radicals, r e action 3. 17 (Braun et al., 1967). Tbe rate constant for suc h a r eaction i s 1 . 2X1011 dm3 mole- 1 s - 1 (Braun et al., 1967) . An additional source for CH radicals could b e a t tribute d to the reaction of hydrogen atoms with methylene radicals, reaction 3.18; the rate constant for this reaction is 1 . 6X1011 drrr mole- 1 s - 1 (Tsang and Hampson, 1986). 143 :Clb + :CI:-k Ckib + I:-b (3.16) :CH? -f : CH? C:z I:-b (3.17) :CH::2 + H? : CH? + I:-k (3.18) Formation of hydrogen cyanide may be explained by t,b p, reaction of CH"" and Cfk radicals with atomic nitrogen, reactions 3.19 and 3 . 20, respectively (Bos - sard et al., 1983) . Tbe rate const,ants for these reac- tions have been estimated by Yung et al. (1984) to be 3 . 0X1010 exp( - 250/T) dm"" mole- 1 s - 1 (13). CH-:!- ? + . N: HCN + Ek (3.19) :Cfk + . N : HCN + H? (3 . 20) Nitrogen atoms may also dimerize to form molecular r1itrogen, reaction 3.21; the rate constant for such a reaction has not yet been determined; however, it can bP. estimated by analogy to the dimerization of atomic oxyeen (reaction 3.22), for which a rate constant has been de1.,ermined to be: 1.9X107 exp(900/T) dm6 mole- 2 s - 1 {Tsang and Hampson, 1986). 2 ? N: N::z N::z {3.21) 2 ? O? + N::z ~ (3 . 22) 144 He ac ti o n s 3 . 1 to 3.22 were used to mode l the elec- tro lys i s of a Cfl,. N2 - ~ 0 mixture. Fi~ure 3 . 19 s hows t h e comr1u-Ce. .. z .._, z 6 u g 0.12 u (.') 0.20 8 J: (_') 50 J: 0 ~ 0 "0' 0 .00 D 0 .0 0 ().0 0 0 .2~ o.-,o 0.7!) 1.00 0.00 0.25 0.50 0.75 1.oi nNz / ( nCH~ + nN2 ) nH2 / ( nCH4 + nN2 + nH2 ) Preigure (TDrT) 200 -4-00 600 BOD C d ,-.. .z.- ,, 0.30 z u 0 ~J 0.15 :r ., '-..I 0 0.00 0.00 0.25 0.50 0.75 1.00 50 7';) 100 Temperoture ("C) Figure 3.20. Computer Simulations of the Effect of G 0 (HCN) on the molar ratio of Nitrogen (a), Hydrogen (b), Wai.er vapor (c), and also as a function of Pres- sure and Temperature (d) during the Electrolysis of a CII" - N:2-IkO mixture. Lines were obtained by ACUCHE/'1 u s ing reactions 3.1 to 3.22 (a, c, d) or to 3.23 (b) and were plotted with published values for comparison: ( 0) Raulin et al. (1982); and (?)Stribling and Miller (1987). The total concentration of gas phase molecules was 2.879X10- 2 mole dm- 3 except in (d) where it was changed accordingly to the pressure change. The concentration of water vapor was kept constant in all cases (7.19Xl0-3 mole dm- 3 ) except in {c) where it was varied from Oto 2.879Xl0- 2 mole dm- 3 . CI:k and N:z were present in one- to-one ratio, except in (a) where it varied from Oto 1. 147 Raulin et al. (1982) and Stribling and Miller (1987) have investigated the effect of added hydrogen on the synthes is of HCN. Computer simulations can al so pre d ict the trend of G 0 HCN a s a function of the molar ratio for f:b, but i t is n eces s a ry to include rea ction 3.23 t o a ccoun t for the electric discharge dissociation of mo l e cular hydrogen into atomic hydroge n. Figure 3.20b shows the computer trend (solid line) and expe r i - mental behavior (symbols) as a function of the molar rati o for hydrog en. G 0 (HCN) is maximum in the absence of hydrogen ; the yield decreases line arly with increas - ing hydroge n content, d own to zero in the absence of eithe r CH4 and N:z. A similar behavior as for hydrogen , is al so predi c ted as a fun c tion of water vapor in the mixt ure {see Figure 3 . 20c ). e Eb H? + H? (3.23) The total pressure and/or the temperature of the gas phase does not have any effect on the yield of HCN as predicted by computer simulation (see Figure 3.20d). Therefore , G 0 (HCN) depends only on the composition of the gas mixture . 148 The r esults presented in this sec tion ma y be u sed to p redict the p ossible produc tion of HCN by el ectri c discharges i n t h e pr i mi tive atmosphere ; this is con- sid ered in the n e xt s ection. 3.3 . 3 . HCN synthesis in primitive atmospheres At mosphe ric elec tricity is a poorly understood s ubje ct , and no estimates of the energy available from this s ource on the primitive Earth have been made {Fox a nd Dose , 1977) . Cloud electrification is a common phenomenon in the atmosphere {Weast, et al., 1986) . The p roces s es by wh i ch the charges in a cloud are fo rmed and separated are still controversial (Bar- Nun, 1981 ) . The most common forms of neutralizat ion of the c h a rge s are by intracloud and cloud- to- ground dischar- ges ; the former is more abundant by a factor of about 5 than the latter {Bar- Nun, 1981) . It has been calculated that the flux of lightning ene r gy i n the tropics is about l . 5X1019 J cm- 2 y- 1 {Bar- Nun, 1981) . In mid- latitudes, the number of lightning strokes is much smaller and in high lat i tudes t hunde rstorms are rare {Bar- Nun, 1981). The globally ,------- 149 averaged energy input by electrical discharges projected per square- centimeter surface area of the Earth is about 16.74 J cm- 2 y - 1 (Fox and Dose , 1977). This value includes 3.77 J cm- 2 y- 1 due to lightning and about 12 . 55 J cm- 2 y- 1 due to corona discharges from pointed objects (Miller and Orgel, 1974). ? The quantity of electric discharges that occurred on the primitive Earth is unknown (Fox and Dose, 1977); however, it is generally accepted that electric dis - charges were as abundant in the primitive atmosphere as are today; and therefore, an energy flux of 16.74 J cm- 2 y - 1 has been extensively used to estimate the electrosynthesis of organic compounds in the primitive Earth (Miller and Orgel, 1974; Stribling and Miller, 1987) . A dose rate of 1.04X102? eV cm- 3 y- 1 or 13.9 MGy y- 1 is derived for electric discharges. The composition of the primitive atmosphere is not known precisely (Fox and Dose, 1977); although it is usually accepted that an atmosphere composed mainly of methane, nitrogen and water could be reasonable for the initial stages of development of the Earth (Ferris and Chen, 1975; Bossard et al., 1982); however, there is 150 no estimate on the relative proportions of CH N 4 , ~ and EkO in the primitive atmosphere. In this study we carried out an extensive 1?n ves- tigation to determine the chemical yield of formation of HCN from the electrolysis of a gas mixture composed of C~ (36.7%), N:z (36.7%), and EkO {26.6%); and a G0 {HCN) =0.26 was experimentally determined. Computer simulations predict that 0 G (HCN) could vary from a factor of 2 to 4 depending on the composition of water vapor or nitrogen in the gas mixture, respectively_ If we consider these effects, G0 (HCN) could have varied from 0.1 to 0.4 depending on the initial nitrogen and water vapor content of the primitive atmosphere. If we take into account the dose rate of electric discharges in the atmosphere, and consider that G 0 (HCN) could have varied from 0.1 to 0.4 depending on the atmos - pheric conditions, an estimate for the rate of electro- synthesis of HCN in the primitive atmosphere would vary Table 3.3 summarizes the chemical yields of HCN synthesis and estimated HCN production in the primitive atmosphere triggered by various energy sources. Shock waves are the most efficient energy source for the synthesis of HCN, followed by ionizing radiations and 151 Table 3. 3. Est imated productions of HCN in the pr1m1.- ti v e a-Lmosphere by diffe r ent energy sources. En e rgy source G O ( HCN) HCN production (mole dm- 3 y - 1 ) UV light negligible 1 n eg ligibl e Ionizing radiation: 0.01 - 0 . 9 2 Upper atmosphere 10- 0 - 10- 4 Lower atmosphere 10- '5-10- 1 Electric Disc harges 0.1 - 0 . 43 Shock waves 3. 2 4 Tota l 1 . Raulin et al., 1982; 2. Zhdamirov et al., 1970; 3. This study; 4. Bar- Nun and Shaviv, 1975. electric discharges. Ultraviolet light has been found to contribute negligibly to the HCN synthesis (Raulin et al. , 1982). The estimated total production of HCN in the primitive atmosphere would b e about 1010 - 1019 Fegly et al., (1986) have derived a HCN rainout rate of (3 - 14)X1011 mole y- 1 ; this indi - cates that HCN synthesized in the primitive atmosphere would be efficiently transported into the oceans by precipitation processes. 152 The physical conditions of the primitive Earth are not, known but it is believed that they did not change significantly during the history of the Earth (Ferris and Hagan, 1984). If we assume that the volumes of the at.mosphere and hydrosphere have not changed, and if all UCN synt.hesized in th e atmosphe re would b e dissolved in the hydrosphe re without any loss due to HCN hydrolysis and / or oligomerization, the concentration of HCN in tbe hydrosphere would vary from 10- 3 to 10-- 2 mole dm- 3 y- 1 ? These figures represent only the upper limits since HCN would b e hyd rolyzed and/or oligomerized in hydrosphere; howe v e r, as it was discussed in Chapter 1 (section 1.3) , these values could be increased by concentrating mec hanisms suc h as the eutectic freezing of water volumes (Sanchez et al., 1967) and by the formation of stable cyanocomplexes as ferro and/or ferricyanide (Na- varro- Gonz.alez , 1983; Navarro- Gonzalez et al ., 1989). 3.4. Conclusions Thi s study describes and characterizes three new me U1od s 1,o . OCQoo. -X Cb .. ..... .. ..... '?? .... ' 0 ?????'. t0 1 .9 I C'l E "D 1 .8 (1) 0.0 2 .0 4.0 6 .0 8.0 0 E ,.__,.. 2 . b ,---, I "< 0 .6 o? ..--,.., I") 0 I 0 .3 0 E o ?? -0 n., ~ 0 .0 Cj) 0.0 4.0 8.0 12.0 0 E --..__; 9.0 b r--, I I"') 6 .0 ....<...O.... z (.) --..__; Q) L.. '---' 0 . 0 200 400 600 BOO Dose ( kGy) Figure 4.5. Dose dependance of the concentration of hexacyanoferrate {III) at low (a) and high dose of ~- irradiation of a 0.02 mole dm- 3 K4Fe{CN)6 aqueous solu- t ion. Dotted line was obtained by computer simulation of the reaction mechanism. Hexacyanoferrate{II) may also react with hydrated electrons leading to the formation of hexacyanoferrate {I), reaction 4.7. The rate constant for this reaction has been determined to be ~7X104 dm3 mole- 1 s- 1 in pulse radiolysis experiments carried out by Zehavi and Rabani {1974). 166 (4.7) h hexacyanofe rrate(II) ct wit Hydrogen atom s rea the triple car - rogen in ing to the ad dition of hyd lead ction 4.8); th e d (rea -nitrogen bond of the ligan bon ined by Zehavi and Rabani te constant w as determ ra and was found e radiolysis e xperiments, (l9 74) in pul s d Rabani, 197 4). 7 drrr mole- 1 s- 1 (Zehavi an to be 3.9X10 )aFr- (4.8 ) Fe(CN H + ses, by two consec utive proces followed Reaction 4.8 is (reaction 4.9 ), e first one is first order of Which th ediate s a reaction o f the interm d one i While the sec on avi and (reaction 4.10 ) (Zeh ith hexacyano ferrate(III) ~ Rabani, 1974) . (4.9) N)aH4 + Fe(C N)a3 - Fe(C - CN (4.10) CN)a 4 - + H Fe(CN)~3- + Fe( 9 was found to te constant o f reaction 4. The ra 03 s- 1 at pH 0-4 (Zeh avi from 1CJ5 to 1 decrease with pH used in and the lowe st value was ), and Rabani, 1974 of eriments, whil e the rate computer simu lation exp 167 1 1 reaction 4.10 was set to 2x10~ dm 3 mole- s - (Zehavi and Rabani, 1974). Zehavi and Rabani (1974) have suggested that pen- tacyanoferrate(I) may react with hexacyanoferrate(II) and (III), reactions 4 . 11 and 4.12, respectively. The authors found that the rate constants of these reac- tions are pH dependent, and a value of 1X108 dm3 mole- 1 s - 1 {Zehavi and Rabani, 1974) was used in computer modeling . N)6 3 - --+ Fe(CN)o3 4 Fe(CN)o 4 + Fe{C - + Fe{CN)6 - {4.11) - Fe(CN)o4 4- + Fe(CN)6 ---+ Fe(CN)e 3 - + Fe(CN)6e- {4.12) Pentacyanoferrate{II) formed in reactions 4.10 to 4.12, may react with hydrogen peroxide {reaction 4.13) or its anion (reaction 4.14) leading to the formation of pentacyanoferrate(III). The rate constants for these reactions were determined by Davies and Garafalo (1976) to be: 107 and 3 dm3 mole- 1 s - 1 , respectively. Fe(CN) o 3 - + Ek< ,...._, n I 4.0 E "O ,,,,.,-P Q) 0 2.0 ..E_,, ,....., N 0. I L...I 0 .0 0.4 0.8 1.2 Dose ( kGy) F igure 4 . 8. Dose dependance of the concentration of mo lec ular hydrogen in the irradiation of aqueous solu- tions of hexacyanoferrate(II). Dotted line was obt a i n e d by computer simulation . The results presented. in this section clearly d e monstrate that hexacyanoferrate(II) is very reactive toward free radical reactions. However, the cyano groups coordinated to the metal ion do not participate in the formation of organic compounds as does when it is not coordinated.. Therefore, it seems likely that HCN might have not played a significant role on the synthesis of organic compounds in the primitive Earth if it was mainly present as a ligand coordinated to me tal ions such as iron . 174 The effects of hexacyanoferrate(II) in aqueous mixtures containing free HCN and hexacyanoferrate(II) are discussed in the next section. 4.3 The ~-irradiation of aqueous hexacyano- ferrate(II)-HCN mixtures Aqueous, oxygen-free solutions (pH :.:-f 6) of 0.1 mole dm- 3 HCN and 0.001 mole dm-3 K,.Fe(CN)o were irradiated with 0 ?Co-ti' - rays in a wide dose range up to about 100 kGy. Hydrogen cyanide is initially present in solution in its non-dissociated form; however during the irradiation of the aqueous solution, the pH increases from 6 to about 7.5 (Figure 4.9). Under these condi - tions, hydrogen cyanide is reversibly converted into cyanide ion according to reaction 4.23; the rate con- stants of both forward and reverse reactions are known, and are: 5.2X104 dm3 mole- 1 s- 1 , and 3.7X1

Q.llfllF-----,--------.------.-----' .z.... ... 0.0 5.0 10.0 15.0 DoGe ( kGy) Figure 4.15 . Formation of N~ during the irradiation of a solution of 0.1 mole dm- 3 HCN and 0.001 mole dm- 3 K .... Fe( CN )6. 183 HO - CH =N? or - o - CH=N? leads to the formation of cyanic acid or cyanate with rate constants of l.1Xl09 dm3 mole- 1 s - 1 (reaction 4.35) (Behar, 1974; Buchler et al., 1976) or 6.5Xl0'3 dm3 mole- 1 s - 1 (reaction 4.36) (Behar, 1974), respectively. Cyanic acid or cyanate undergo hydrolysis and/or ammonolysis leading to the formation of CO:z, NH3 and urea (reactions 4.37 to 4. 43). The second order rate constants (k) given in these reactions have units of dm3 mole- 1 s- 1 and were taken from Wilson, 1972a. In addition, the acid-base equilibrium reactions for ammonia, reactions 4.44 and 4.45 (Castellan, 1983), and for carbon dioxide, reac- tions 4 . 46 and 4.47 (Weast et al., 1985) were taken into account in computer simulations. 2HO- CH =N? HO-CH=NH + HO-CN {4.35) 2- 0 - CH=N? HCN + HO- + -o-CN (4.36) HO- CN - o - CN + Er" pKa=3.46 (4.37) fkO HO- CN + g+ N~ ... + 00:z (k=6. ox10- 2 > (4.38) HO-CN + EkO Nib + 00:z (k=7. 8Xl0- 4 ) (4.39) EkO -o- CN + 7EkO Nib + He~- (k=8.3Xl0 ) (4.40) 184 fh-, O - o - CN + H. ... Nib + CCb ( k =43 } (4_ 41) HO-CN + NR"3 Nf-bCONI-b (k=17 . 8} (4 . 42) - o - CN + NH. .. .... NI-b CONH:z (k=9. 6X1 0-- 0 } (4 . 43} NB:,, + f-b O NH. .. .... + HO- ( k =6. OXl O'.l C 4 ' 2 .. 6 0 (1 ? ., 4 $') 2 () !> a:: 6 12 0 6 12 0 6 12 0 Time (min) grams of DNPH deriv atives of Figure 4.16 . Gas chromato (b) and 95 (c) kGy in aldehydes formed a t 30 (a), 61 HCN and 0.001 mole of 0.1 mole dm- 3 al- :g~;ous solutions DPH; 3. Form K4Fe(CN)6. Legend : 1. Solvent; 2 . m ehyde; 5. Propional dehyde; 6 . Butyr- dehyde; 4. Acetald 8ldehyde. identified by gas chromatography bu~Yraldehyde were shows the dose and mass spectrome try. Figure 4.18 e formation of form aldehyde. The dependence of th ulated emical yields of fo rmation were calc ~adiation ch 2 ; 0 3 (CihCH0)=2.70X10- 0 to be : G (HCH0)=9 .10Xl0- ; G 26X10-3 0 CH0)=7. ? G? fhCfkCH0)=2. 09X10- 3 ; G (Cfk (Cfk )2 (C 186 91 122 100 8 ~2 162 163 224 210 50 106 198 131 281 ~ td 266 Q) 0... 0 60 100 160 200 250 300 Q) rn aj m ....., 199 0 100 b ~ 226 263 50 211 104 165 267 0 50 100 150 200 250 300 m/e Figure 4 . 17. Electron impact (a) and chemical ioniza- tion (b) mass spectra of the DNPH fraction formed at 95 kGy . Molecular weights of DNPH derivatives: Formal - d e hyde 210; acetaldehyde: 224; propionaldehyde: 238; butyraldehyde:252; valeraldehyde : 266 . 187 LO 0 12.0 X ,......, 0 ,..., I E 8.0 'O ., 0 E 4.0 ,..., 0 I u 0.0 I 0 25 50 75 L.....J Dose ( kGy) Figure 4.18. Formation of formaldehyde as a function 01 dose in aqueous solutions of 0.1 mole dm- 3 HCN and 0 . 001 mole dm- 3 K4Fe(CN)6. A variety of polycarboxylic acids are formed in irradiated HCN solutions. Figure 4.19 shows gas chromatograms obtained at different doses. Except for oxalic acid (peak 1 in Figure 4.19), their concentra- tions exponentially increased with dose. Figure 4.20 shows the variation of concentration of oxalic and malonic acids with dose. The identifications of poly- carboxylic acids were confirmed by mass spectrometry using electron impact and chemical ionization techniques {Figure 4.21) . Table 4.1 summarizes the identification type for each carboxylic acid as well as their initial radiation chemical yield. 188 a ? 20 40 b e ,, 0 20 40 ~ 7 C 0 4 ? a ,, 0 20 40 Time ( min ) Figure 4.19. Gas chromatograms of polycarboxylic acids methyl esters from HCN - f<...Fe(CN)o mixtures at 30 (a), 46 (b) and 95 (c) kGy. 1. Oxalic; 3. Malanie; 4. Fumaric; 5. Succinic; 6. 1,2 Dimethyl succinic; 7. Maleic; 8. Glutaric; 10. Malic; 11. Carboxysuccinic. 189 ~ ~ ~ ~ X X 4.0 ,... ~ 4.0 n n I I E E 'O '0 ., ., 0 0 E 2.0 2.0 E ....,, ,...... ,....., "O "O ?.:; u 0 0 .',,! 0 C: ~ 0. 0.0 0 0 75 100 c ........ ::E D ~ Dose ( kGy) Figure 4.20. Dose dependance of the formation of oxalic and malonic acids in HCN- ~Fe(CN)o mixtures. 190 Oxalio methyl e&1ter 59 119 .:,t 100 100 ..>t. 0 0 .ID. Cl) Q.. 0..,. ~ 50 50 Ill C-D m 0 0 ~ 118 105 100 150 200 50 100 150 Molonlc methyl eeter 59 100 ~ 100 .:1 0 Q Cl) Cl. 101 Cl) a.. Cl) ., ~ 50 50 to 0 C Q'.l ...D.. 0 't ~ 1.32 101 ~ 100 150 200 50 100 150 Moleic methyl ester 102 ~ 100 172 100 _,., g 59 69 i a.. 145 0.. a,; 50 102 50 I 0 C CD ID 113 134 I 100 150 200 50 100 150 m / e m/e Figure 4.21. Electron impact (left) and chemical ionization (right) mass spectra of polycarboxylic acids formed at 95 kGy. 191 Table 4 . 1 . In i tial chemical yields (GD) of formation of p o lycarboxylic acids . Acid Identificat i on GD type* Oxa lic ,6 D ... 2 ... ?? 2.34X10- Malanie 4 D ? Fumar i c ~ D ...... 2.24Xl0- Succ in ic D ?? 3 . 39Xl0-15 D. . 61Xl0-15 1 , 2 - Dime thyl succinic 6 D ... 6 0.2 :J "U 'in (J) I,... -- -- -CT c -- -- -B--- -EJ- -- -- --0. 0 0 25 50 75 100 Dose ( kGy) Figure 4 . 22. Dose dependance for the formation of oligomeric material measured as total dry residue {a) and as specific elemental composition in the residue (b). Legend for b: Carbon ( 0 ); nitrogen ( ~ ); hydroge n ( D ); and oxygen ( 7 ). Figure 4.23 shows the effect of the concentration of hexacyanoferrate(II) on the free radical oligomeri - zation of HCN. The amount of HCN decomposed is only 194 ""' 0.06...----------------~ n I E "O Cl 0 E ......... .--, z u .I._ . I 0.02-+--+-t' - --~--------r---------j 0 Figure 4.23. Dependance of the decomposition of HCN as a function of the concentration of ~Fe(CN)o in aqueous solutions of 0.1 mole dm- 3 HCN that have been expos e d to a dose of 61 kGy. markedly decreased at high concentrations of Fe(CN)o 4 - Figure 4.24 shows gas chromatograms of a variety of polycarboxylic acids that are synthesized. Their yields are only markedly decreased at high concentrations (10-3 mole dm- 3 ) of Fe(CN)o 4 - but are not significantly affected at con- centrations 210-4 mole dm- 3 (Figure 4 . 24). 195 0 1 J 7 C G) II) C. 0 a 1/) II (I) a::: 11 10 ,2 0 20 4-0 0 20 40 I J 7 b I J 7 11 d G.) Ul C 6 0 Q_ Ill (U a::: 0 20 40 0 20 40 Time ( min ) Time (min) Figure 4.24. Gas chromatograms of polycarboxylic acids formed in the b' - irradiation of an aqueous solution of 0 . 1 mole dm- 3 HCN and K4Fe{CN)6 that have been exposed 61 kGy . The concentration of K4Fe {CN)6 was set to {mole dm- 3 ): 1X10-3 {a); 1X10-4 (b); 1Xl0- 6 (c); and no additive {d). Based on the previous discussion it is evident that the effect of hexacyanoferrate{II) on the free - radical oligomerization of HCN was only its inhibition. No new synthetic pathways nor enhancement in prcxlucts yields was influenced by the reactions of hexacyanofer- rate{II) with the free- radicals produced from BCN. On the contrary, bexacyanoferrate{II) only led to the re- formation of BCN by its reaction with the radicals: 196 HO- CH =N? and -o- CH=N? (reactions 4.32 and 4.34). Therefore, in contradiction to the general ideas of the possible catalytic role of cyanocomplexes in chemical evolution (Beck, 1978), it seems reasonable to conclude that hexacyanoferrate(II) did not play a role in the oligomerization of HCN via free-radicals. The next section examines the effect of hexacyano- ferrate(III) on the free-radical oligomerization of HCN. 4.4 The ~-irradiation of aqueous hexacyano- ferrate(III)-BCN mixtures Aqueous, oxygen-free solutions (pff.::i 6) o:f 0.1 mole dm?- :s HCN and 0. 001 mole dm- 3 K:sFe(CN)6 were irradiated with 6 ?Co-d-rays in a wide dose range up to about 100 kGy. Figure 4.25 shows the decomposition of HCN as function of dose. The initial radiation chemical yield of decomposition of HCN was determined to be 2.00. This value is considerably smaller than that obtained 197 0.12~------ --------- I") I ([) ??? ?? ?Q .... . Q E ?? ???? .. Q ..... .P .. . "O O.OB ?? ? ?? I) 0 E 0 .04 .--, z u I '--' 0.00+---------------~ 0 25 50 75 100 Do:1e ( kGy) Figure 4.25. Decomposition of HCN in solutions of HCN and K~Fe(CN}6. Dotted line was obtained by com- puter simulation {see text for details}. in pure aqueous HCN solutions {G 0 =7.60}. The predicted rate decomposition of HCN {dotted line in figure 4 . 25} is consistent with experimental results. Modeling of this system was carried out taken into account reactions 1 to 54 given in Tables 2.3 to 2.7 {Chapter 2, section 2.7) as well as reactions 4.2 to 4 . 54 considered in sections 4.2 and 4.3 in this chapter. The decomposition of HCN was attributed to reac- tions 4.24 to 4.31. A possible reaction {4.56) of 198 hexacyanoferrate(III) with EkC=N? radicals could ac- count for the marked decreased of G 0 {-HCN) in this system: EbC=N? + Fe(CN)o 3 - Fe(CN)o 4 - + H~ + HCN {4.56) This reaction was included in computer modeling in order to derive its rate constant; however, computer simulations indicated that such a reaction does not participate in the overall reaction mechanism of decom- position of HCN, and therefore, was ruled out. The marked decrease in G 0 ( - HCN) is due to a scavenging effect of the species, H? by hexacyanofer- rate{III): About 64% of hydrogen atoms react with hexa- cyanoferrate(II) according to reaction 4.18, and 36% react with hydrogen cyanide according to reaction 4.24. Figure 4.26 shows the decomposition of hexacyano- ferrate(III) as a function of dose. The radiation chemical yield of decomposition of hexacyano- ferrate(III) was determined to be 0.18. Computed and experimental trends agree quite well at low doses {<5kGy). Hexacyanoferrate(III) mainly disappears by 199 I"') 1.01ar-----------------, 0 ,--, X I .......... ~ .....c..:.o.. I") z I 0.5 u E '-" -0 t.::i L . ..J a.> 0 E: .__, 0.0 ..-----.-------.-----.-- - ----j 0 20 40 60 80 Dose ( kGy ) Figure 4.26. Decomposition of hexacyanoferrate(III) in solutions of HCN and K"5Fe(CN)o. Dotted line was obtained by computer simulation (see text for details. reactions 4.18 and 4.19. At higher doses, the computed and experimental trends deviate. This might be due to reactions of pentacyanoferrate(II) with the free radi - cals from water radiolysis. Figure 4.27 and 4.28 show the experimental dose dependence for the formation of hydrogen (G 0 =0.45) and predicted formation of hydrogen peroxide (G 0 =0.57), respectively. Their mechanism of formation was given in detail in Chapter 2, section 2.7, and is summarized by reaction 4.1. 200 LO 0 6.D X 0 2.0 .D ? ..E..., 0 'rN 0.0 L...J 0.0 0.4 0.8 1.2 Dose ( l< ,..._ n I E -0 2.0 ~ 0 E ....... ,...., N ~ 0.0 ? r.: 0.0 2.0 4.0 6.0 1..-1 Dose ( kGy) Figure 4.28. Computer simulation for the formation of hydrogen peroxide during the ~-irradiation of an 3 oxygen- free, aqueous solution of 0.1 mole dm- HCN and 0.001 mole dm- 3 ~Fe(CN}o. 201 Carbon dioxide (Figure 4.29) and ammonia (Figure 4.30) are the major products formed from HCN- f6Fe(CN) o mixtures. Their initial radiation chemical yields of format ion are 1 . 00 and 1.10, respectively. Their mec hanism of formation may be explained by reactions 4.35 and 4 . 36, followed by reactions 4.37 to 4 . 43. In addition, ammonia may also formed by hydrolysis of IkC=Nfl (reaction 4 . 48) and HO- CH=NH (reaction 4.49) . ..,. 0 X 1.5 ......,_ r<') I E , .o "'O I) 0 E 0.5 ..._, ,......, <'4 0 o. .(....).. 0.0 0.4 0.8 1.2 Dose ( kGy) Figure 4.29. Formation of carbon dioxide during the 'ti- irradiation of an oxygen-free, aqueous solution of 0.1 mole dm- 3 HCN and 0.001 mole dm- 3 lbFe(CN)o. Dotted line was obtained by computer simulation. 202 ,,., 0 1.5 X ...-.. I') I 1.0 E "O ~ 0 E 0 .5 ......., ~ : 0 . zr 0.0 5.0 10.0 15.0 L-1 Dose ( kGy} Figure 4.30 . Effect of dose on the concentration of ammonia formed radiolytically in 0.1 mole dm- 3 HCN and 0 . 001 mo l e dm- 3 K3 Fe(CN)o aqueous solutions. Figure 4.31 shows gas chromatograms of aldehydes formed in the radiolysis of HCN- I6Fe(CN)o mixtures. Formaldehyde (G 0 =9.10X10- 3 ), acetaldehyde (G 0 =3 . 41X10-2 ); propionaldehyde (G 0 =7.19X10-3 ), and butyraldehyde {G 0 =1.75X10- 2 ) were all identified by their retention times and coinjection of standards in gas chromatography, and also by their fragmentation pattern by electron impact and chemical ionization mass spectrometry (Figure 4.32). Their radiation chemical yields were calculated from the slope of concentration vs dose plots as that shown in Figure 4.33. The con- centration of formaldehyde linearly increase with 203 a b C 2 Q} 4 (/) C 46 0 6 D... (/) 2 3 Q} 5 Cl'.: 5 2 0 6 12 0 6 12 0 6 12 Time (min) Figure 4.31. Gas chromatogram of 2,4-dinitrophenyl - hydrazone derivatives of aldehydes formed at 27 (a), 61 (b) and 95 (b) kGy of b' - irradiation of an oxygen- free, aque ous solution of 0.1 mole dm- 3 HCN and 0.001 mole dm- 3 K3 Fe(CN)6. See legend in Figure 4.16 for peak assignment . irradiation dose. Its mechanism of formation may be e x p lained by reaction 4.48. 204 a 11 1112 JOO ... ao lOtl .!?II P"-. 0 lid1tt 60 lOO t6D 200 2110 300 Al ?IQ .p ..:.i un 0 100 b a,a a( 111 1DO - atw 0 M 100 l!W 800 2ffO 300 m/e Figure 4.32. Electron impact {a) and chemical ioniza- tion {b) mass spectra of the 2,4- dinitrophenyl - hydra - zone fract ion formed at 95 kGy in the HCN- K"" Fe(CN)o sy s t e m. IJ') 0 >< 12.0 I') 0 I E 8 .0 "U CV ci E '4-.0 ....., ,---, 0 I u 0.0 :r: 0 25 50 75 100 L.......l Dose ( kGy) Figure 4 . 33. Dose dependence for the formation of for - maldehyde in the HCN - K:sFe(CN)o system. 205 Oxalic acid was identified as the major polycar- boxylic acid formed in irradiated HCN-KsFe(CN).,, solu- tions (F igure 4.34). A variety of other acids were , 0 e .. ,o,, 0 20 40 b ~ C CJ 0 a_ 7 In CV a::: ? 11 0 20 4-0 3 7 C 0 Time (min) Figure 4 . 34. Gas chromatograms of polycarboxylic acids methyl esters formed at 27 {a), 46 {b) and 95 {c) kGy in the HCN - K3 Fe{CN).,, system. See Figure 4.19 for peak assignment . 206 also formed but in low yield. Their identification was based on retention times and coinjections with stan- dards in gas chromatography. Supporting evidence was based on the GC-MS analysis carried out for the HCN- K...Fe(CN)6 system (Section 4.3). Table 4.2 summarizes the initial radiation chemi - cal yields of formation of polycarboxylic acids. 0 Table 4.2. Initial chemical yields (G ) of formation of polycarboxylic acids. Compound 2 Oxalic l.60X10- 4 Malonic 1.1sx10- Fumaric 3.90X10-!5 Succinic 2.57X1Q-!5 1,2- Dimethyl succinic 9.95Xl0- 6 Maleic l.14Xl0- 4 Glutaric 2.83X10-!5 Adipic 5. 12Xl0- 6 Malic 5.87Xl0- 6 207 The con c entrat ion of polycarboxylic acids exponentially i n c r e ases wi th dose with the exce ption of oxalic acid . Figure 4 . 3 5 s h ows these b e h a v i or f or oxa l ic and ma leic acid s. ?0 ~ r X X ~.o 4.0 ~ n I I [ E ~ ~ C w 0 0 E 2.0 E ~ u ~ ~ u 0 0 ~ ~ 0 0. 0 .0 ~ .~.... . 25 50 75 100 ! D 4 ~~ (~) Fig ure 4.35. Dose dependance of the formation of oxali c and maleic acids in the HCN - KsFe(CN)6 system. Oligomeric material was formed in low yield in HCN- K3 Fe{CN)6 mixtures (Figure 4.36). At low doses its is very low but exponentially increases with irradia- tion dos e . The initial radiation chemical yields of 208 formation for this materia l we r e determined to b e: G 0 {C) =0 . 87; G 0 (N) =0.82; G 0 (H) =0.88; and G 0 {0) =0.09. ,......_ I") 2 .0 I E. a u Q\ '-./ 1.0 _____,-9- (IJ ::, u ' in D~ Q) e? ~ ~ c 0 0 .0 0 25 50 75 100 ,--. I") I I r 0.4 b 'O ' c,, ......, Cl) ::, "C 'in .Q..J. ?:- 0 0 2b 50 Dose ( kGy) Figure 4.36. Dose dependance of the formation of oligomeric material measured as total dry residue {a) and as specific elemental composition in the residue (b): Carbon { 0 ); nitrogen ( 6 ); hydrogen ( D ); and oxygen ( 7 ) . 209 The results given in this section indicate that h exacyanoferrate(III) does not lead to an enhancement o f p r oducts yields nor contribute to new synthet ic pathwa y s . Hexacyanoferrate(III) markedly inhibits the d ecomposition and oligomerization of hydrogen cyanide . The mechanism responsible for such an effect is due to the s c avenging of H- and e..q by hexacyanoferrate(III), r eactions 4 . 18 and 4.19. Therefore, it seems r easonable to conclude that hexacyanoferrate(II) did not play a role in the oligomerization of HCN via free - radicals . A comparison of the results obtained in the radio- lysis of hydrogen cyanide in the presence or absence of h exac yanoferrate(II) or (III) is given in section 4.4 . 210 4.5 The effect of cyanocomplexes in the free-radical oligomerization of BCN: Implications to chemical evolution Tabl e 4 . 3 summarizes the initi a l radiation chem i - c a l yields of major produc ts obtained in the radiolys is of aqueous solutions of 0.1 mole dm- 3 HCN in the Table 4.3. Initial radiation chemical yields (G 0 ). Compound Go HCN No additive K. .. Fe( CN )b ~Fe(CN)6 -HCN 7.60 2.00 2.00 F'e( CN ) 6 3 - 4. 10 - 0.18 Ek 0. 3S- 0.45 0.45 CO:z 0. 80- 0.60 1.00 Nlb 2. 4?>- 1.01 1.10 HCHO 0. 18 0.009 0.009 ( COzH) 0.22 0.022 0.016 Oligomeric material: ( - HCN) 5.80 1. 22 0.87 -nraganic et al., 1973. 211 absence or presence of 0.001 mole dm- 3 of ~Fe(CN) 0 or f6 Fe(CN)o. The decomposition of hydrogen cyanide is markedly dec reas e from 7.60 to 2.00 in the presence of hexa- cyanoferrates. Carbon dioxide and ammonia were the major products in the presence of cyanocomplexes. The yi elds o f othe r products such as formaldehyde and oxalic acid were decreased at least by an order of mag nitude . About 76% of decomposed HCN molecules are bu i lt up into larger molecules in pure HCN solutions. The yield of this oligomeric material decreases to about 61% or 44% in the presence of hexacyanofer- rate (II) or (III), respectively. The effect of hexacyanoferrate(II) or (III} on the radiolysis of aqueous solutions of HCN was the marked decrease in the yield of decomposition of HCN. In the case of hexacyanoferrate{II), this effect was attributed to the re-formation of HCN by the reactions of Fe{CN)o 4 - with HO- CH=N? and - o - CH=N? , reactions 4.32 and 4.34: 212 HO- CU =N? + Fe(CN)6 4 - --+ Fe{CN)6 3 - + HCN + HO- {4.32) - o - CH =N? + Fe{CN)o 4 - --+ Fe{CN)o 3 - +HO-+ - cN (4 . 34) Whereas in the case of hexacyanoferrate(III), the ma rke d decreas e in G 0 ( - HCN) was attributed to a scav e ngin~ effect of H? by hexacyanoferrate(III): 64% of hydrogen atoms react with hexacyanoferrate{II) ac- cording to reaction 4.18, and 36% react with hydrogen cyanide according to reaction 4.24. 4 ---t_ _: Fe(CN)6 - + w? {4.191:1) H? F e (CN) a 3 - Fe(CNh, 3 - + HCN {4.19b) H? + HCN I-bC=N? {4.24) In spite of the differences in mechanisms between hexacyanoferrate(II) and (III), the overall effect is the same, a marked decrease in G 0 ( - HCN). Furthermore, the initial yields of certain products such as H:.z, NH3 and HCHO are almost similar, and at higher doses of ir- radiations the differences in product yields becomes less apparent between the two systems (see for instance GC profiles for carboxylic acids, Figures 4.19 and 4. 34). The lack of a significant difference between the two systems, particularly at high doses, is due to 213 the fast interconversion of hexacyanoferrate(II) into (III) and viceversa. Therefore, after a few Gy of radiation dose, both hexacyanoferrate(II) and {III) are participating in the reaction mechanism. The effect of other cyanocomplexes of transition elements on the free- radical oligomerization of HCN is expected to be similar to that found for hexacyanofer- rate(II) and (III). This is because the rates and modes of reactions of cyanocomplexes of transition elements with hydroxyl radicals (Figure 4.37), hydrogen atoms (Figure 4.38), and hydrated electrons (Figure 4.39) are very similar. 214 ,,......._, 11 10 -r-- I (/) Pi(II) 9 tvlo( IV) Os(II) Fe(II) .- 10 g 0 0 I 0 0 (I) Ru(II) Ni(II) Au(!) 0 E 5 10 t'1 E 'U 3 ....._.,. 10 ~ 1 2 3 4 5 6 7 8 9 10 Number of d electrons Figure 4.37. Dependance of the rate of reaction {k) of hydroxyl radicals with cyanocomplexes of transition eleme nts. Values were taken from Ross and Ross, 1977, and Buxton et al . , 1988 . 215 , 1 ---- 10 ,-- I rt) Ni( II ) 9 Fe(III) 0 ,-- 10 0 I -Cl) 0 5 0 Pt(II) Cd(II) E 10 Co(III) n E -0 3 ....__ 10 _y 1 2 3 4 5 6 7 8 9 10 Number of d electrons Figure 4.38. Dependance of the rate of reaction (k) of hydrogen atoms with cyanocomplexes of transition ele- ments. Values were taken from Anbar and Ross, 1975, and Buxton et al., 1988. 216 11 ..-.... 10 r- fl) Mn(II) Co(ll) g Hg(l i) Mo(IV) Ni(II) Au(I) ,--- 10 0 Fe(III) 0 0 I Cr(l11) 0 0 -Q) Cr(II) Co(III) 0~ 0 Pd(II) 5 0 Cd(II) E 10 Pt(II) Cu(I) 0 Zn(II) r<) [ TI J ..,___, 10 ..:::L. 1 2 3 4 5 6 7 8 9 10 Number of d electrons Figure 4.39. Dependance of the rate of reaction (k) of hydra~~d electrons with cyanocomplexes of transition elements. Values were taken from Anbar et al., 1973, Ross, 1975, and Buxton et al., 1988 . 217 4.6 Conclusions A comprehensive investigation of the effect of hexacyanoferrate(II) and (III) on the free-radical oligomerization of HCN was carried out. The t - radio- lysjs of aqueous solutions of 0.02 mole dm- 3 K4Fe(CN) 6 , 0.1 mole dm- 3 HCN in the presence of (10-~-10-3 mole dm- 3 K4Fe(CN)o or 10-3 ) K~Fe(CN)o were investigated in a wide dose range. The radiation chemical yields for a variety of products were calculated, these included among others: fb, C.N~ Ch~m- 63 , 793. Fegly, B.Jr., Prinn, R.G., Hartman, H. and Watkins, G.H.: 1986, Nature 319, 305. Ferris, J .P., Donner, D.B. and Lotz, W.: 1972, _.J_,__Am__ Ch~ m. Soc. 94, 6968. Fe rris, J.P. and Ryan, T.J.: 1973, J. Org. Che~ . 38, 3302. Ferris, J.P., Sanchez, R.A. amd Mancuso, R.W.: 1973, Organic Syntheses 5, 32. Ferris, J.P. and Nicodem, D.E.: 1974 , in Dose, K. , Fox, S.W., Deborin, G.A. and Pavlvskaya, T.E. (eds.), Th~ Origin of Life _and Evolutionary Biochemistry, Plenum Press, New York, 107 . Ferris, J.P., Williams, E.A . , Nicodem, D.E., Hubbard, J.S. and Voecks, G.E .: 1974a, Nature 249, 437. Ferris, J.P., Wos, J.D . and Lobo, A.P.: 1974b, J. Mol . _Evol. 3, 311. Ferris, J.P. and Chen, C.T.: 1975, Natur~ 258, 587. Ferris, J.P . and Joshi, P.C.: 1978, Science 201, 361. Ferris, J.P. and . Edelson, E.H.: 1978, J. Org. Chem . 43, 3989. Ferris, J.P., Joshi, P.C., Edelson, E.H. and Lawless, J.G.: 1978, J. Mol. Evol. 11, 293. Ferris, J.P., Alwis, K.W., Edelson, E.H., Mount, N. and Hagan, W.J.Jr .: 1981a, in Wolman, Y. (ed.), Origin of Life, D. Reidel Publ. Co., 125. Ferris, J.P., Edelson, E.H., Auyeung, J.M . and Joshi, P .C.: 1981b, J. Mol. Evol. 11, 69 . 232 Fe rris, J . P . , Hagan , W.J.Jr . , Alwi s , K. W. and Mc Crea , J_ : 1982 , J . Mo ]. Ev o ] . 18 , 304. Ferris, J.P. and Hagan, W.J.Jr: 1984, Te trah ~d ~Qn 40 , 1093. Fi e l d, F.H., Franklin, J.L. and Munson, M.S.D.: 1963, J: . . Am. _ Che m. So c . 85, 3 5 7 5 . Fox, S.W., and Dos e, K.: 1977, Mo lecular Evolution and tb e Oxia in _of _L_ife , Re vis ed Ed . , Ma r ce l Dekk e r Inc., ? Ne w York, 370 p_ Fripiat, , J . J. and Cruz- Cumplido, M. I.: 1974, hPQ .,_ _ R~v . .E;g rth ...?nd Plane t . Sc i_ 2, 239. FrohP.n F . W.: 1974, Be r. _ Bunse n - Gesellschaft _Phys. Gb e m. 78 , 18 4 . Fug i o, C.: 1930, Bu] 1. Chem ~ Soc. (Japan) 5, 249. G~rriso n, W.M., Morrison , D.C., Hamilton, J . G., Be n s on, A.A. and Calvin, M.: 1951, Sc ience 114, 416. Go rd o n , R. and Ausloos, P . : 1967, J . Ph s. Chem . 46, 4 823 . Grot h , W. and Suess, H.: 1938, Naturw i sser!schaften 26, 77. Ha is s insky, M., Koulkes, A.M. and Masri, E.: 19G6, L Chim. Phys. 63, 1129. Haldane, J.B.S.: 1929, The Origin of Life, The Rat jonalist Annual. {in J.D.: 1967, The Origin of Life, We idenfeld and Nicolson, London, 242). Hanel, R.A . , Conrath, B., Flaser, F.M., Kunde, V., Maguire, W., Pearl, J., Pirraglia, J., Samuelson, R. , Herath, L., Allison, M., Cruikshank, D., Gautier, D., Gierasch, P., Horn, L . , Koppany, R. and Ponnamperuma, C.: 1981, Science 18, 157. Harada, K. and Fox, S.W.: 1964, Nature 201, 335. Harada, K. and Iwasaki, T.: 1975, Chem. Letters 185. 233 Ha rada, K. , Suzuki, S. and Ishida, H. : 1978, Expe ~ j~Q - tifl 34 , 17. Harada, K., Has e gawa, N. and Shimoyama, A. 1984, Pri_g ins of Life 14 , 11 5 . Ha rkins , W. and ,Jac k s on , J . M.: 193~3 , J. Ch e m. Phys. 1, 37. Uart eck, P . and Dond es , S. : 1957, J. Ch e m. Phys. , 27, 5 4 G. Harte ck, P . and Dondes, S.: 1958a, J...,__ C.hem . Phy$. , 28 , 9 7 5 . Ha r1-eck , P. and Dondes , S .: 1958b, J:~ y s . Ch e m., 63, 95G. Harting, R . , Troe , J. and Wagner, H.G . G.: 1971, Utb furn )O_~ ium on Combu s t i on, Combus ti on Inst itute, P it,t s burg b , 147. Ha rting , K . J . and Getoff, N.: 1980, J._f.b_Qtochem . 13, 207 . Hauser , W. P. : 1964, J. Phys . Chem . 68 , 1576. Hoing , E.R. and Sheppard, C. W. : 1946, J . Ph s. Che ID . 50, 119. Hug , G.L.: 1981, Optical Spectra of _Nonme1,a l l .ic Ing r_-:: gani Q... j'ran1;i ients S ecies in A ueous So ]JJt ion, NSRDS- NBS-69, Nat i onal Bureau of Standards, Washington , D. C. Jame s, A.D. and Murray, R.S.: 1975, J. Chem , Soc . Dal _ ton 1530. Kammaludin, Singh, M. and Deopujari, S . W. : 1986, Origins of life 17, 59. Kevorkian, V., Heath, C. E . and Boudart, M. : 1960, J:.,__ Phys. Chem. 64 , 964. Kliss, R. M. and Matthews, C. N.: 1962, Proc. Nat. Ac ad. Sci. USA 48 , 1130. 234 Kobayashi, K. and Ponnamperuma, C.: 1985a, Ori_g_ins of 1-,i f.~ 16, 41. Kobayashi, K . and Ponnamperuma, C.: 1985b, OrifLins of J, ife 16, 57. Kobay ashi, K. and Ponnamperuma, C.: 1986, Viva Origino 14 , 41. Kohe rstein, E.: 1973, Ind. En_gng_. Chem. Process. Des. J)~ v . 12, 444. Kondratiev, V.N.: 1965, 10th Sm osium on Combustion, Combustion Institute, Pittsburgh, 319. Koz lov, G. I. and Knorre, V.G .: 1963, RtJ.~? J. Phys. _Chem . 9, 1128. Kraljic, I. and Trumbore, C.N.: 1965, L_hw.,.___i;bem. Soc . 87 , 2547. Kuhn, W.R. and Atreytl, S . K.: 1979, lcaru~ 37, 207. Lampe , F.W.: 1957, J. Am. Chem. Soc . 79, 1055 . Laufer, A. and Bass, A.: 1978, Combustion and Flame 32, 215. Lasana, A.C., Holland, H.D. and Dwyer, M.J . : 1971, Science 174, 53. Levine, J.S.: 1982, J ,_ Mol. Evol. 18, 161. Levine, J.L., Augustsson, T.R. and Natarajan, M.: 1982, Qrigins of Life 12, 245. Lind, S.C. and Bardwell, D.C. : 1926, J. Am. Chem. Soc. 48, 2335. Lind, S.C. and Glockler, G.: 1929, J. Am. Chem. Soc. 51, 2811. Lind, S.C. and Glockler, G.: 1930, J. Am. Chem. Soc. 52, 4450. Lind, S.C. and Schltze, G.R.: 1931, J. Am . Chem. Soc . 53, 3355. 235 Loew, G.H. and Chang, S .: 1971, TetrahedroQ 27, 298 . Loew, G.H.: 1971, J . Th e or . Biol. 33, 121. Lone, D.A., George, W. O. and Williams, A. E.: 1960, l:_,:-oc. _Qhero__,_ Soc_,_ , 285. Lowe, C.U . , Rees, M. W. and Markham, F.R .S .: 1963, N_<1 = :t;ure 199, 219. Mackenzie, F.T. : 1975, in Riley J.P. and Skirrow, G. ( eds. ), ~hemical Oceanography, Academic Press, London, 1. Magee, E.M .: 1963, J. Chem . Phy~. 39, 855. Mahan, B . H. and Mandal, R.: 1962, J. Chem. Phys _ 37, 207. MHlmstarlt, H.V., Enke, C.G. and Crouc h, S.R.: 1981, EJ-~ctr_gnics and Instrumentation for ScieQtists, The Benjamin/Cummings Publishing Co., Inc., Men lo Park, p. 543. Masri, E. and Haissinsky, M.: 1963, J ~ hj m. Phys. 60, 397. Matthews, C.N . and Moser, R.E.: 1966, roe. Nat. Acad. ~p_i. 56, 1087 . Maurin, J.: 1962, J. Chim. Phys . 59, 15. Mc Govern, W.E.: 1969, J. Atmos. Sc i. 26, 623. Miller, S.L.: 1953, Sc ience 117, 528. Mi ller, S.L.: 1957a, Biochim. et Biophys. Ac ta 23, 480 . Miller, S.L.: 1957b, Ann. N.Y. Acad. Sci. 69, 260. Miller, S.L.: 1974, in Oro, J., Miller, S.L., Ponnamperuma, C. and Young, R.S. (eds.), osmochemical Evolution and the Origins of Life, I.D. Reide l Pub . Co., Holland, p. 139. 236 Mi ll e r, S.L., and Orge l L.E .: (1974) Th e _O rigins _of J!i f ~ on ....: P1~ E,u-.th, Prentice- Hall, Inc., Ne w J ersey, 229 P- Moffat , J . B. : 197 5, .J ._ Chem . Soc. Ch em. Commun . 888. Maurey, D., Raulin, F. and Toupance , G.: 1981, in Wo l - man , Y. ( ed. ), Qrjgj_ n o:f_J,j fc;, D. Reide l Puhlishincr Co ., Holland, 73. Narty, F.: 1981, in Vennesland, R., Conn, E.E., Know- less, C.J., West ley, J . and Wissing, F . (ed s ), .Qyanide in...B_igJ_~, Academic Press, London, 115. N~varro- Gonzalez, R.: 1983, ldentiLication _of Carboxylic Acids in Irradiated Samples of Sim le Nit,rj]es and C_.Yanides, Bachelor Science Thesis, Univer- sidad Nacional Autonoma de ~ xico, C.U. ~xico. Nava rro-Gon zalez, R. and Neeron-Mendoza, A. : 1983a, An,J.lY;;;_j_s of Ca.f':hony L Comp_Q_unds. Part_L_ Pape_r Chroma - J.,onraph_y ..Qf .... 8ld ~h_yq ~? ~nQ ....K~t....9n~:c;, Technical Report -Q- 05-83, Center of Nuclear Studies, U. N.A . M., C.U., Mexico D.F. 08500. Na varro- Gonzalez, R. and Negron- Mendoza, A. : 1983b, Anal_ysis of Carbonyl Com___Rounds. Part II. Ga s Chroma - toi;;raphy of Aldeh_ys}es and Ke-Lones, Technical Report Q- 06-83, Center of Nu c lear Studies, U.N .A. M. , C.U., Mexico D. F. 08500. Navarro- Gorn!:a Jez , R., Negron- Mendoza, A. and Ponnam- peruma, C.: 1986, Origins of Life 16, 301. Navarro- Gonzalez, R., Negron- Mendoza, A., Aguirre - Calder6n , M.E. and Ponnamperuma, C.: 1989, Advances 1n s.r>ace Research in press. Nazly, N. , Collins, P.A. and Knowles, C.J.: 1981, in Vennesland, B., Conn, E.E., Knowless, C.J., Westley, J . and Wissing, F. (eds), Cyanide in Biology, Academic Press, London, 289. Negron- Mendoza, A. and Navarro- Gonzalez, R.: 1982, Determination of Carbox lie Acids b Gas Chromato- g:r~@...Y., Technical Report Q-12- 82, Center of Nuclear Studies, U.N . A.M ., C. U., ~xico D. F . 08500 . Negr6n-MP-ndoza , A., Navarro-Gan~ lez, R., Draganic, z. and Draganic, I.: 1982, ResumeQes d el IV Simposio sabre ~\ct!!JlL!:~~ _lh!_r~l tJ~L. Rad 10-9!:l ::i m::i Qg __ y _Q.y_:j_mj _Q_Q _d e Rad i acignes, U.N.A.M . , Mexico , 34. Negr6n- Mendoza , A . , Dragard c, Z. D., Navarro- Gonzalez, R . and Dragan ic, I.G.: 1983, R~d - Res. 95, 248. Negr6 n - Mendoza, A. and Draganic, Z.D.: 1984, Adv. Spac~ Res. 4, 12L Negr6n- Me ndoza , A., Na v arro - Gonza l ez , R . and Azamar, J . A. : 1984, Analysis of Aldehydes. Keton es , Alcoholes, .Bnd _t,ioI)oc~rboxyJ i Q. __ AQjQ~ b:y____ _fo~~--~o l id Chrom_ntography, Technical Report Q-04-84, Center of Nuclear Studies, U. N. A. M. , C.U., Mi;>xico D.F. 08500. Neta, P. and Fessenden, R.W.: 1970, J_ Phys_ Chem_ 74, 3362. Njkc-t, jc , V ., Draganic, Z., Neskovic, S. and Dragani c, I.: 1982, J. Mo l. Evol. 18, 130. Njkctic, V., Draganic, Z . D., Neskovic, S. and Draganic, I.G.: 1983, J. Mo l. Eva]. 19, 184. Niketic, V.: 1984a, Origins of Life 14, 251. Niketic, V.: 1984b, Adv. Space Res. 4, 107. Ogura, H. and Kondo, M.: 1967, Bull. Chem. Soc. Japan 40, 2448. Ogura, H.: 1968, Bull . Ch~W- S09. Jap~n 41, 2871. Ogura, H., Fujimura, T., Murozono, S. Hirano, K. and Kond o, M.: 1972, J. Nuclear Sci. Technology 9, 339. Oparin, A.I.: 1924, Proiskhozhdenie zhizny, Moscow . Izd. Moskovshii Rabochii. {Translation in Bernal, J.D. 1967, Jhe Origin of Life, Weidenfeld and Nicolson, London, 199). Oparin, A. I.: 1938, The Origin of _l,_ife, Dover, New York. 238 Opa rin, A.I. : 1957, The Origin of Life on th~ _E~rth, 3rd Ed . OJ iver & Boyd, Lond on . Oparin, A.I.: 1972, in Ponnampe ruma, C., (ed.), Com~ Qarative Plane tolggy, Academj c Press , Ne w York, 1. Or6, J. 1961, Nf.!t,_~re 191, 1193. Oro , ,T. a nd Kiml ia ll, A.P .: 1961, Arc b. Bioc h e 111_. ~iQ1>h y ~ . 94 , 217. Or6, J . and Kam f:l t , S . S.: 1961, N~ture 190, 442. Oro, J., Miller, S.L., and Urey, IJ.C.: 1977, inBuvet,, R. , ed., Living Systems as Ene r Converter~, Elsevier/North- Holland Biomedical press, Amste rdam , 7 . Owe n, T.: 1982, J_. Mol. Evol . 18, 150. Palm, C. and Calvin, M.: 1962, J. Am. Che m. SoQ. 84, 2115. Palm~ r , H.B . and Hirt, T,J.: 1963, J. Ph s . Cb~W- 67, 709. Park, W.K., Hoc hst,im, A.R . and Ponnamperuma , C. : 1975, Origins of Life 6, 99 . Perrin, D.E., Amnre~o , W.L. and Perrin, D.R. : 1966, Purification of Laboratory Chemicals, Pergamon Press, Oxford . Pet,ers, K. and Wagner, O.H.: 1931, Z. Phys_ Chem. A153, 161. Ponnamperuma, C.: 1965, in Fox, S.W. (ed.), The Origj_n Qf _ Prebiolo ? stems and of Their Molecular Mat,rices, Academic Press, New York, 221. Ponnamperuma, C., Lemmon, R.M., Mariner, R. and Calvin, M.: 1963, Proc. Nat. Acad. Sci. USA 49, 737. Ponnamperuma, C. and Woeller, F.: 1964, Nature 203, 272. 239 Ponnamperuma, C. and Pering, K.: 1966, NaJ~11r~ 209, 979 . Ponnamperuma, C., Woeller, F., Flores, J., Rom ie~ , M. and Allen, W.: 1969, Adv. Chem Ser. 80, 280. Rabani, J . and Ma theson, M. S. : 1966, ~J.,_ Eh:t.s. _{;h e m. 70, 761 . Ilabani, J. and Me y e r s t,e i n , D. 1968 : J . .l 'hys. Che D1. 72 , 1599 . Rao, V.V., Mackay, D. and Trass, 0 . .Ch Q.m _,___EQfI . 4 5 , 6 1 Rasmussen, 0. L . and Bjergbakke , E. : 1984, f:HE~SIMlU.c.. I-) P rogram Packa e for Numerical s_i~ulatlon of Chg mic~l E~action Syst~w~. Ris?- R- 395 , Ris~ National Laboratory, DK- 4000 Roskilde, Denmark. Raulin, F. and Toupance , G.: 1975a, Q_..r-:ir: :i n s of I.:ife 6, 91. Raulin, J<'. and Toupancc , G. : 1975b, Origins of Life 6 , 507 . Hau ] in, F., BossRrd, A. and Toupance, G. 38, 358. Raulin, F., Maurey, D. and Toupance, G. 1982, Oriain.? _gf_ Life 12, 267. ReLbe rt, R.E. and Ausloos, P. 1972, J__. Pbotochem . 1, 171. Rebbc r-L, R.E. and Ausloos, P . 1973, J _.. ___Re ~. N.B.S. 77A , 109. Roscoe, J.M . and Thompson, M.J.: 1985, Int. J. Chem. Kipetics 17 , 967. Ross, A.B.: 1975, Se]ected__fu)ecific Rates o~_Reactions _g_f_Jransients From Water in Aqueous Solution. 1. H dra - ted_E_l~ctron.,__Qypfilemental D~t9 , NSRDS- NBS- 43 Supple- me nt, National Bureau of Standards, Washington, D.C. 2 40 Ross, F. and Ross , A.B.: 19 77, .Se l ected Speci f ic Rates of React.jon s _ of Transients l'.:_rom _Wa ter _jn Ag?eou ,; Sol u - .ti.Qn. J I L liYdroxyl__R~_di_QgJ_ _ansLPerhy(jroxy] Radi.9al a nd lht:! ir Ha d i~:R l Ions, NSRDS-NBS- 59 , Nati on a l Bureau o f Sta ndard s, Wash i n~ ton, D.C . Ross , A . B. and Ne 1,a, P: 1979, Ra 1,e_Con:.t,an t s for _Re a Q..-::- 1,i q_n s of Ino r gaQ_ic Radicals in AIBJe ou s Solu tion, NSRDS- NBS- 65 , Na t.. i ona l Bu rea u of Stand ards, Wa s hington, D. C. Ross , A. IL a nd Neta , P: 1982 , Rcrt_~ Co n stants f o r R~ar, = tion s o f Alip h at i c _Carb on- Cente r e d Radica l s in ~ 91!...~ou~ S g l t!t;i9n, NSRDS- NBS-70 , Na t iona l Bu reau o f St a nda rds, Was h i nat.on, D. C. Sanche z, R .A . , Ferris, J . P . a nd Orge l, L.E.: 1966, Sc j_ ~ _ri cQ 154 , 7 84. San c h ez , R. A., Ferr i s , J .P. a n d Orge l, L . E. 1967, s.L_ _ Mol _. Bj o l. 30, 22 3 . Sch aefer, F . C. : 197 0 , in Zuirap popo r t, S. ( ed. ), Tb e ~ h e mj~tJ:.y_of the Cyano Group, Interscience , N. Y. , 240. Sharpe , A . G. : 1976, Th e Chemistry of Qy~p__Q_ C,q_1fil~.l ~:x~x-~f .tbe Transit i on Me tals , Ac ade mic Pres s , London, 113. Shimo y a ma, A., Blair, N. and Ponnampe ruma , C.: 1978, 1n Noda , H. ( ed), Origin of Life, Jap an Scienti f ic So c ieties Press, To kyo , p . 95. Sid l e, A. B. : 1967, Nature 216, 408. Sieck, L.W. and Johnsen, R.H . : 1963, J. Phys. Chem . 6 7, 2281.. Ski nn e r, G. B . and Ruehrwe in, R. A.: 1959, J . Phys. Che m. 63 , 1736. S l a nger, T . G. and Black, G. 1982, J. Che m. Phys_ 77, 2 432. Smith, R . M. and Ma rte ll, A.E.: 1976, Crit ical Stability Constants, Volume 4 : Inorganic Complexes, Plenum Press , Ne w York, 26 . 241 Snyder, L.E . and Buhl, D.: 1971, Astrophys. J. Le tt . 163, IA?. Spinks J. W. T. and Woods, R.J.: 1976, An Introductc;i.Qn t o Radjat,ion _Cbe mistry, J ohn Wil ey & Sons, Ne w York, 503 p. St,anley, J.B., Brown, D. F., Senn, V.J. and Dollean, F.G.: 1975, J . Food Sci. 40 , 11 34. Striblinc-, R. a nd Miller, S.L.: 1986, Qri_g_;i.ns Qf__ Life 16, 279. Striblin(r, R. and Miller, S.L .: 1987, Ori_gjns of ... Life 17 , 26 1. Sun, H., and Weiss l er, G.L.: 1955, ~- Ch em . Phy$. 23 , 1160. Swallow, A.: 1960, Radiation Ch e mistry __ of_ Or_gaqj_g CQ!!I --: ounds, Pergamon Press, Ne w York, 45, 243. Tain- Jen, Y., Sutherland, R. G. and Verra ll, R. E .: 1980, Can J . Cbem~ 58 , 1909 . Terasawa, J . and Harada , K. : 1980, Ch ~m- Letters 73. Thomas, J.K. and Hart , E . J .: 1962, Rad. Res . 17, 408. Tokunaga, A. T., Beck, S.C . , Geballe, T.R . , Lacey, L .H. and Serabyn, E.: 1981, lea ~~~ 48 , 283. Torres, J.L., Negron - Me ndoza, A. and Castillo, S. : 1982, CaUhration of Carbq.n_Monox1d~ _ ..M~thane Ethane and Carbon Dioxid e b_y Gas Chromatography, Technical Report Q-11-82 , Ce nter of Nuclear Studies , U. N.A. M. , C.U., Mexico D.F . 08500 . Toupance, G., Sebban, G. and Buvet, R.: 1970, J . Ch i m. Phys. Physicochim. Biol. 67, 1870 . Toupance, G. , Raulin, F. and Buvet, R.: 1975, Origins of Life 6, 83 . Toupance, G., Bossard, A . and Raulin, F. : 1977, Q;rigins of Life 8, 259. 242 Tsang, W. and Hamp~on, R.F.: 1986, J. Phys_ Che~. Re f. D_fl_tn 15 . 1 087. Ulich, B. L. and Conklin, B. K. : 1974, Nat.l)_:r~ 248 , 121. Ve nnesland, B .? Pistorius , E.K. and Gewitz, H.S.: 1981, in Ve nnesland, B., Conn, E . E., Knowless, C. J., Westley, J. and Wissing, F. (eds), Cyanide in Biology. Academic l'resf; , London, 349 . Voe-L, A.B. and Schwart;z, A.W.: 1982, Orig:ip s of Life 12. 45 . Von-L, A.B . and Schwartz, A.W. 1983, Bioor_g_ Che{D. 12, 8. Volker , T.: 1957, Ang~~-- Chem. 69, 728. Volker, T.: 1960, Angew. Chem. 72, 379. Weast , R.C., Astle, M.J. and Beyer, W. H. (eds.): 1985, ~RC Handbook of Chemistry and Physics, CRC Press, F] orida. Weissler, A.: 1953, J . Acoustical Soc. Ame r . 25, 651. Wjener, H. and Burton, M.: 1953, J. Am. Chem. Soc. 75, 5815. Wiley, E.J.: 1934, Trans. Farady Soc. 30, 230. Wilson, W. E.: 1972, J. Phys. Chem. Ref. Data 1, 535. Wilson, I.R.: 1972a, in Bamford, C.H. and Tipper, C.F.H ( eds ), .Qh~I!!!_Q?!l Kinetics, Elsevier Publishing Company, Amsterdam, 217. Yang, W., Minard, R. and Matthews, C.: 1976. J. Theor. Biol. 56, 111. Yuasa, S. and Ishigami, M.: 1975, Origins of Life 6, 75. Yung, Y.L., Allen, M. and Pinto, J.P.: 1984, Astro- phisical J. Supplement Series 55, 465. 243 Ze h a vi, D. and Raban i, J . : 1972 , J. Phys _ Che m. 76 , 3703. Zehav i , D. and Ralla ni, J.: 1974, J ,__:PbYf. . __che I!). 78, 1368 . Zhdamj rov, G. G. , Kornienko, V. N. , Bo rosov, E . A. , Rumyantsev, Y.M. and Dzanti e v, B. G.: 1971, High En e rgy Ch e mi ?: t,ry 5, 475. CURRICOLUM VITAE Name? Rafae l Navarro - Gon z ~lez Permanent address: Inst i t uto d e Ci e n c ias Nucleares , Uni v ersidad Nc1cional Aut,onoma d e Mex ico, Ciudad Universitaria, Mex ico D.F. 04.510 MEXICO lJl:/; ret: aTJd d at,e Lo h e conferred: Ph.D., 1989 DuLc? uf hiri ,h: April 25, 195 9 PJace of birth : Me xico, D. F., Mexico Col lcr;iai.,e lnstituL.ior is attended: D/\TE OF DEGREE Onjvcr:.,; jt,y of' Maryland , Collet;e Park 1984- 1989 Ph.D . 1989 Biuc bemisLry tJniVf )Y-;,jdad Nac iona] Autonoma de Me xico, Me xico City 1978 - 1981 B. Sci. 1983 Bi o logy llonor~;: 1 . 1980- 1982. Research Fellowship from Th e Nationa] Autonomo us University of Mexico. 2. 1984- 1989. Researc h F e llowship from Th e National Autonomous University of Me xi c0. Me mbe rsb.ips: 1. International Soc iety for the Study of Origins of life . Associate Me mbe r, 1981 - 2. Mex ican Chemical Society. Member, 1986- 3. Radiation Researc h Society. Student Member, 1986- Puli l i cat ions : 1 . Ne;Gr6n - M<.in