Plasma Modeling and Prebiotic Chemistry: A Review of the State-of-the-Art and Perspectives

We review the recent progress in the modeling of plasmas or ionized gases, with compositions compatible with that of primordial atmospheres. The plasma kinetics involves elementary processes by which free electrons ultimately activate weakly reactive molecules, such as carbon dioxide or methane, thereby potentially starting prebiotic reaction chains. These processes include electron–molecule reactions and energy exchanges between molecules. They are basic processes, for example, in the famous Miller-Urey experiment, and become relevant in any prebiotic scenario where the primordial atmosphere is significantly ionized by electrical activity, photoionization or meteor phenomena. The kinetics of plasma displays remarkable complexity due to the non-equilibrium features of the energy distributions involved. In particular, we argue that two concepts developed by the plasma modeling community, the electron velocity distribution function and the vibrational distribution function, may unlock much new information and provide insight into prebiotic processes initiated by electron–molecule collisions.


Introduction
Earth's chemical history and the origin of life are two of the most crucial topics in astrochemistry and astrobiology [1][2][3][4]. Two main theories [5,6] have been proposed regarding the emergence of life on Earth: (i) exogeneous delivery and (ii) endogenous synthesis. According to the first theory, prebiotic molecules have reached the Earth's surface by means of comets, meteors and asteroids. The so-called delivery, that is, the actual transport of molecules from space to Earth, represents the key point of the Panspermia theory [7][8][9][10][11]. It was demonstrated that carbonate and sulfate micrometeoroids (especially calcite and anhydrite) may act as good prebiotic molecule carriers [12][13][14], in view of the well-known association between carbonates/sulfates and organic matter presence [15][16][17][18]. Recently, the role of the chemical composition of the primordial atmosphere was pointed out, in the context of the atmospheric entry of the aforementioned carbonate and sulfate micrometeoroids [19].
On the other hand, the endogenous synthesis approach, developed in the 1950s and 1960s, speculates that primordial atmospheres of terrestrial planets were filled up with methane, ammonia and hydrogen, and it theorizes that the primordial clouds provide the necessary chemical and energetic conditions to activate small molecules and link them into bigger molecules, in order to synthesize the building blocks of life. The topic of early Earth's atmosphere as a source of prebiotic matter comes from the Miller-Urey experiments [20][21][22], partially modeled on Oparin's hypothesis of reducing primordial atmosphere [23]. So, the endogenous theory states that the synthesis of organic molecules occurred directly on Earth starting from simple parent molecules. Actually, a third theory concerning the origins of life has been proposed, related to the production of organic molecules and even amino acids in hydrothermal vents [24].
Miller-Urey's experimental apparatus (left side of Figure 1) consisted of two glass flasks: the lower one was partially filled with liquid water-to mimic the primordial ocean-, and a gas mixture composed of methane, ammonia and hydrogen-reproducing the primordial atmosphere-filled the other balloon. In this primordial gas, Miller placed two electrodes to simulate a lightning-like electric discharge. Therefore, Miller and Urey recreated the disequilibrium chemistry conditions by applying an electrical discharge to the gases, which probably formed a highly reducing primordial atmosphere. The two flasks in the apparatus were connected in a hydrological cycle, which allowed it to mimic rain after the electrical discharge. The analyses demonstrated that, after the experiment, the primordial ocean was rich in carbon compounds, which contained significant amounts of amino acids. Miller and Urey demonstrated that organic molecules can form spontaneously from inorganic components, in such primordial reducing environments. It can be stated that the Miller-Urey experiment is the first abiotic synthesis of organic compounds, under simulated primordial Earth conditions, in the context of the studies on the origin of life.  Although recent progress in astrochemistry and astrobiology suggests a role of increasing importance of the activation process by UV radiation and shock waves [26,27], the following considerations must be taken into account: (i) The Miller-Urey experiment represents a milestone in prebiotic chemistry, from the conceptual, historical and pedagogical point of view; therefore, a better understanding of the experiment can add new interesting results in prebiotic studies; (ii) The chemical activation by shock waves, like that produced by meteors, proceeds, in the case of high Mach numbers, through the ionization of the gas. Fast energy electrons are still responsible for the chemical activation of the neutral gas. A better understanding of the Miller-Urey experience phenomenology can also help, therefore, in the scenarios of chemical activation by shock waves; (iii) The activation mechanism by electric discharge needs to be better understood in light of the growing knowledge of electrical phenomena in planetary atmospheres, such as those of Venus and Titan.
With respect to the studies carried out in the few years following Miller experiment, considerable progress has been made in terms of theoretical and methodological understanding (see [28,29] and references therein). Moreover, hypotheses regarding the composition of the primordial terrestrial atmosphere have changed considerably.
The composition and depth of the primordial atmosphere are the subject of a lively debate in the planetology community. Dating back to Oparin [23] and Haldane [30], most scientists believe that, at the very least, the primordial atmosphere was enriched in reducing low-molar-mass components (methane, ammonia and even hydrogen). Additionally, there is no reason to assume that an early Earth atmosphere (a mixture of components with very different molar mass, like nitrogen and hydrogen) should have uniform composition in its lower and upper regions: probably, the early atmosphere was enriched in low mass components at high altitudes.
Hypotheses regarding the chemical composition of the Earth's earliest atmosphere are numerous. Moreover, it is believed that there were many atmospheres over the vastness of Earth's geological time [31]. Geologically based argumentations, that treat the atmosphere as outgassed from solid Earth, indicate that the Earth's primitive atmosphere was composed mostly of H 2 O, CO 2 and N 2 , with traces of CO and H 2 [32][33][34].
Since the 1980s, geoscientists have been disputing Miller's hypothesis of highly reducing primordial atmosphere [35]. Many experiments were carried out using CO and CO 2 atmospheres [36,37], showing that a CO-dominant atmosphere might have produced organic compounds; however, the synthesis of organic compounds by the action of electric discharges in neutral gas mixtures is much less efficient, in this case. Results presented by Cleaves et al. [38] demonstrated that neutral atmospheres can produce amino acids in much higher amount than previously thought. A complete and exhaustive review about the origin and conditions of the Earth's earliest atmospheres can be found in the work by Zahnle et al. [39], with particular attention to the issue concerning the origin of life (∼4.6-4 Gys ago).
Under the hypotheses that the primordial Earth's atmosphere consisted predominantly of CO and CO 2 , the chemical channels previously considered for a reducing atmosphere (based on the preliminary production of the highly reactive HCN molecule) must be reformulated and subjected to new studies, bearing in mind that other components, even minority ones, can play a role in subsequent chemical reactions. Although theoretical studies on the CO and CO 2 system in the context of the Miller-Urey experiment are advanced [40][41][42], one aspect of these studies is still rather primitive, namely the role of electrons in the plasma (a word indicating an ionized gas introduced by I. Langmuir [43]) in producing the first excited and radical species and ions, which later enter several chemical channels to form prebiotic molecules. Studies of chemical kinetics performed in the case of gas discharges with chemical composition enormously simpler than the one used in Miller-type experiences have shown that ionization and excitation processes have an extremely complex and interesting kinetics. These studies show that the efficiency of radical formation can be much greater than the one observed in the much more phenomenological studies devoted to Miller's experience so far.
Several studies actually considered the chemical activation via plasma electron of planetary atmospheric compositions with a prebiotic perspective. For example, plasma chemistry has been studied in electric discharges fed with chemical compositions compatible with the present atmosphere of Titan [44][45][46][47][48] and for Venus and Jupiter atmospheres [49].
Titan resembles the terrestrial planets, and particularly Earth, as it is the only other body in our solar neighbourhood with liquid reservoirs on its surface. The familiar role of water in our atmosphere is replaced by methane, which survives as a liquid at the extremely low temperature conditions at the surface (−180 • C). The scheme in Figure 2 reports the chain of reactions of prebiotic chemistry in Titan atmosphere. The fact that plasma chemistry might be relevant to the interpretation of the Miller-Urey experience was mentioned in the 1978 Libby review [50]: in particular, this review discusses the capacity of the methyl cation CH + 3 , formed by the ionization of methane, to produce higher organic species by ion-neutral reactions. At that time, however, there was no possibility to model the kinetics of complex molecular plasma like ionized methane. This was developed more recently, mostly for the study of material science processes where CH 4 is ionized in order to ultimately fix carbon into diamonds [51,52], or carbon nanotubes [53].
All these considerations lead to the the idea of exploring prebiotic molecular mechanisms in the early stages of the production of new chemical species from a plasma kinetics point of view, highlighting therefore the importance of electrons and excited molecules in chemical channels.
The aim of this paper, in the framework of the endogenous synthesis, is to show how recent developments in the theoretical modeling of chemical reactions and transport phenomena, in weakly ionized gases and in thermodynamic non-equilibrium conditions, can shed a new light on the studies on the synthesis of prebiotic molecules in the primordial atmosphere. Actually, in our view, contributions from delivery mechanisms are not completely ruled out, since they can provide useful monomers which can react with active species produced from the native medium. Our primary goal is to fill the gap between plasma kinetics and prebiotic chemistry. We do believe that the original idea of bringing together two "distant" theories may represent a solid base for future important developments.
The paper is organized as follows: in Section 2, we introduce non-equilibrium plasma kinetics and, in particular, the role electrons can have in understanding primordial chemical processes. In Section 3, we present the kinetics of molecular activation by plasma electrons in plasma discharges, with the techniques, which have good chances to be suitable even for the complex molecules in prebiotic chemistry. Using a consistent set of accurate cross sections for many elementary processes, we show a realistic distribution of the kinetic energy of free electrons, assuming a realistic chemical composition for the gas for two primordial atmosphere compositions frequently discussed in previous studies. Section 4 reports an overview on ab initio calculations of cross sections. Final considerations and outlooks are presented in Section 5.

Tools from Plasma Kinetics
Plasma kinetics is the study of chemical reaction rates in ionized gases [54]. It is a very active field of study with many applications in microelectronics, aerospace, medicine and astrochemistry, just to name a few. Thanks to many previous studies [55][56][57], we have knowledge of the physical conditions of a plasma associated with natural phenomena, such as for example lightning discharges.
Plasma kinetics was born with the intent to understand more quantitatively the physics of systems with a considerable degree of non-equilibrium, in particular ionosphere chemistry. Over time, it has developed its own calculation and prediction methods, principally coming from statistical mechanics, which are very specialized with respect to those used in chemical kinetics of gases and in chemical kinetics in general.
In such non-equilibrium (sometimes called "low-temperature", or "cold") plasmas, electrons have a much greater kinetic energy (from few eV to few tens of eV) than that of atomic and molecular neutral species, and distributions of speed and energy are not trivial. These peculiarities are produced by energy flows and energy exchanges within the system, whose main origin is the acceleration of electrons (and, to a lesser extent, of ions) by the electric field; the electric field is applied to the plasma by means of two electrodes at different electrical potential. Electrons maintain a high average energy because of the slow processes of dissipation of electronic kinetic energy, due to the unfavorable mass-ratio between electrons and other species (which slows down energy transfer in elastic collisions), and the small cross sections for inelastic processes (compared to the corresponding elastic ones). The energy distribution of free electrons is therefore a complex function that cannot be reduced to the Maxwell-Boltzmann distribution, as normally assumed in the chemical kinetics of gases. This is of paramount importance, because the violation of the Maxwell-Boltzmann distribution rules out the renowned Arrhenius law, used in kinetic studies of prebiotic processes.
In view of the knowledge produced in plasma science, there is little doubt that the earlier studies on the kinetics of the Miller-Urey experience must be reconsidered. The primary chemical activation of the gas is associated with chemical dissociation and ionization of small monomers molecules present in the original Miller-Urey mixture and modern variants: these molecules are CO 2 , CO, H 2 , NH 3 , CH 4 , H 2 O, NO, and a few others. While a huge number of studies has been produced concerning chemical reactions of activated monomers to produce amino acids, aldehydes and sugars, very little has been studied in the direction of a better understanding of the primary chemical activation. In this process, the presence of high energy free electrons (whose share is to be determined), fast ions, quantum mechanics and statistical mechanics of activation are new fields where progress occurs on a month-by-month basis.
Very effective deterministic, stochastic and hybrid methods have been developed over the years to determine the electron energy distribution, typically as a function of the reduced electric field E /N (E is the magnitude of the electric field, N is the gas number density) and the neutral gas composition [58][59][60][61][62][63][64][65][66]. The strong gap between the average energy of the electronic component and that of the neutral component leads in turn to a chemical kinetics that requires adequate models, so far not very much used outside plasma electronics, aerospace and a few studies on the chemistry of the interstellar medium. These models root in the determination of reaction channels following the flow of energy provided by free electrons, which constantly produces reactive species such as ions, free radicals and excited states of atoms and molecules. Such a state-of-the-art entails concrete possibilities that the know-how implemented and validated in the field of plasma kinetics can lead to remarkable results if applied to the study of those systems in astrobiology where similar circumstances are met. Taking into account specifically the chemical peculiarities of Miller-Urey's experience, particularly in the new formulations that are emerging, these aspects will be illustrated in the following sections by means of two case studies, which, however, only aim to represent a small example of the great variety of useful research that can be performed on the same basis.

Kinetics of Molecular Activation by Plasma Electrons
The theoretical study of reaction rates involving molecules has progressed remarkably since the first formulation of Miller-Urey's experience. A particularly important technique, developed in parallel in plasma kinetics and in astrophysics, is the so-called state-to-state (STS) approach. In this formulation of chemical kinetics, the individual excited levels of molecules are treated as individual species, with their own characteristic reaction data. The STS is able to comply with the fact that the specific reaction rates of the different excited states of most molecules may differ by orders of magnitude. Furthermore, the distribution of populations of the different levels in a system far from equilibrium, such as a cold plasma, can differ considerably from the results obtained in statistical mechanics applied to equilibrium systems, that is, the Boltzmann distribution.
As a result, the mole fraction of tens or hundreds, sometimes thousands, of individual "components", χ i , must be calculated by solving extended systems of ordinary differential equations (ODE's) of the form d dt addressed in this context by the term Master Equation (ME). R i is the reaction rate for the i-th STS component, which is a function potentially of all of the χ's: it includes many terms describing elementary reactions between STS components. The ME is tricky to solve accurately, since the different contributions to its r.h.s. may differ by orders of magnitude, making the system of ODEs so-called "stiff" in computational terms. Nevertheless, many fruitful and concrete applications of this approach have been demonstrated for reactive fluids, lasers and electric discharges [67][68][69][70][71][72]. An example is provided by the chemical activation of CO 2 , which is a very stable molecule, often abundant in the atmospheres of rocky planets. It appears that some processes, under the conditions of a primordial atmosphere, can activate chemically this molecule with the help of physical energy in the form of, for example, kinetic energy of the free electrons. In this respect and in recent years, great progress has been made in the theoretical understanding of the chemical plasma activation of CO 2 , due to the importance that this process plays in the CO 2 recycling study for the mitigation of the climatic consequences of CO 2 accumulation in the atmosphere [73,74].
Since the most effective dissociation path of CO 2 is through the asymmetric stretching vibrational mode, it is necessary to distinguish the 22 vibrational levels of this mode [73]. Electron-impact vibrational excitation processes, the so-called "eV-processes", however, can only pump energy into the lowest vibrational levels, say the first 4 or 5, because cross sections of the direct excitation of higher levels are very low: Subsequently, CO 2 molecules reach the dissociation limit with the help of a different process, named VV 1 , where a vibrationally excited CO 2 molecule interacts with a more excited one pushing it to higher excitation: The mechanism in Equation (3), proposed by Russian scientists about 40 years ago [75], pushes molecules up the internal energy ladder until they reach the dissociation threshold, i.e.: The process described by the reactions in Equations (2)-(4) is extremely effective in increasing the efficiency of activation, since the direct dissociation channel with electronic states as intermediates has an energetic cost in the order of 10 eV, whereas the energy cost for vibrational excitation is less than 1 eV.
In the last few years, many researchers, including two authors of the present paper, aimed to achieve a better understanding of this mechanism using computer simulations. For example, a recent scheme is based on the characterization of the upwards flux in energy space, J [76][77][78]. J is a functional of the function f ( ), which is the distribution of molecules according to their internal energy . The theoretical basis of this approach is the relative smallness of energy difference between consecutive vibrational quantum states compared with the binding energy for the majority of molecules; that is, in many molecules the vibrational levels for a given electronic curve are numerous. For example, in the case of the nitrogen molecule, they are around 50 for the electronic ground state, depending on the potential energy curve assumed; a lower number, 15, holds for the hydrogen molecule.
Most of the chemical processes between excited vibrational levels involve the exchange of no more than one quantum. We are therefore in the transition/continuum regime as schematically represented in Figure 3. Consequently, it is possible to formulate the problem by means of a diffusion equation, or more generally of a drift and diffusion equation, that is the Fokker-Planck equation: It has been shown [76][77][78] that the two transport coefficients in Equation (5), namely the two jump moments a 1 ( ) and a 2 ( ), can be calculated starting from the rate coefficients of the chemical reactions involving the different vibrational states. The calculation is not simple and the use of formulas requires the evaluation of numerous data; nevertheless the calculation is possible and produces results that are reliable and comparable with those based on the explicit discretization of the internal energy. However, discretization requires very long computational times and can hardly be applied to future theoretical formulations of Miller's experience, especially considering that in this experience complex species formed starting from monomers with numerous levels of internal energy are also involved. On the left, the arrow represents a one-quantum collisional transition between two consecutive levels. On the right: the level manifold is approximated by a continuum of energy, the kinetics is sampled by Brownian motions, the energy transport is drift-diffusive. Drawing adapted from [76].
The mathematics of diffusion, namely the continuity equation of the flux through the internal energy space, demonstrates that, at steady-state, the aforementioned flux functional J, whose expression is: where the internal energy drift b( ) and diffusion D( ) coefficients are: is a constant, that is, it does not depend on the internal energy . Furthermore, it can be shown that J[ f ] is the dissociation rate. f ( ) can be easily calculated by solving the Fokker-Planck equation based on rate coefficients of vibrational processes and the known value of J [78]. But, since J needs f to be calculated, the process is iterated until J is self-consistency attained, and the problem is solved. The behavior of f as a function of J can provide additional insight into the process of molecular dissociation. This is shown in the plots in Figure 4, which are based on calculations of the described type: in all cases, on the x-axis we set the internal energy, which goes from the value 0, corresponding to the zero-point energy of the quantum oscillator that describes the asymmetric stretching motion, to the dissociation energy. The plotted functions are precisely the internal energy distributions. The two plots differ in the value of T v , the excitation temperature of the 0 → 1 transition.  Comparing the plot on the left in Figure 4 (lower temperature) with the one on the right (higher temperature) and observing the values of J corresponding to the different curves shown in each of the two figures, two completely different behaviors can be observed. While in the case of the lower temperature, the function can only be modified with considerable variations of the functional J, in the case of high temperature, completely different energy distributions are obtained for very small variations of the functional J. On the one hand, the low-temperature case exemplifies that the energy distribution determines the dissociation rate, especially in the high-energy region, which is closer to where the dissociation takes place. On the other hand, in the high temperature case, a regime is realized in which the functional J, and therefore the dissociation rate of the molecules, changes very little even in the hypothesis of large variations of the high energy part of the distribution.
This last result is very important in a reconsideration of the Miller-Urey experiment with an atmosphere rich in CO 2 . In such an atmosphere, locally produced plasmas may produce radicals from CO 2 with a very high energetic efficiency, which was unforeseeable at Miller's times.
A similar synergistic mechanism could speed up considerably the dissociation of other molecules relevant to modern studies of planetary atmospheres in a prebiotic perspective: one of those is methane (CH 4 ). A continuum model for the internal energy flux in methane is completely lacking at present, but it is within reach of the present knowledge and methods and could advance the understanding of prebiotic chemistry, again, with the help of methods from plasma kinetics.
As regards electron kinetics, the characterization of the distribution of free electrons is fundamental to describe the thermodynamic state and its departure from equilibrium [65]. From the distribution function, it is also possible to calculate macroscopic quantities, such as electron swarm parameters [79], and electron impact rate coefficients. To calculate the electron distribution function, it is necessary to solve the Electron Boltzmann Equation (EBE). Several codes are available as open-source for calculation of the Electron Velocity Distribution Function (EVDF). Most of them are based on the two-term expansion of the EVDF in Legendre polynomials, such as BOLSIG+ [64] and LoKI-B [80]. At the expenses of computational time, other codes employ different numerical methods based on multiterm expansion of the EVDF in Legendre polynomials (e.g., the MultiBolt code [81]) and Monte Carlo (MC) simulations (e.g., the METHES code [82]), in order to obtain more accurate calculations.
Recently, a code based on the Monte Carlo Flux (MCF) method [59] has been implemented and benchmarked by some of the present authors against two-term and multi-term Boltzmann solvers [66,83]. The MCF method provides detailed calculations of the EVDF through a variance reduction technique, that provides a highly computationally efficient solution of the electron transport problem, compared with conventional MC simulations. In MCF, the phase space of electrons is partitioned into cells, and the problem of electron transport is reformulated in the following form [65]: where n i (t) is the number of electrons in the i-th cell at time t, ∆t is the MC time step for calculation of the transition probabilities q ij (∆t). Hence, a set of MC simulations is performed to calculate the coefficients q ij by placing systematically sample electrons in the i-th cell and collecting the number of electrons in any j-cell after a time ∆t [65] ( Figure 5). Then, the time evolution of the electron distribution is obtained by an iterative solution of Equation (8). As an alternative, the steady-state solution can also be calculated by means of an eigenvalue method [59]. The advantage of MCF is that it can calculate EVDFs with uniform statistical accuracy in low density regions that cannot be treated by the MC method. Moreover, in [83], the MCF method has also been extended to take into account the thermal velocity distribution function of gas molecules, together with an accurate description of rotational and vibrational excited states for a pure CO 2 gas.
In the present work, MCF has been employed for simulations of electrons in mixtures of gases representative of the composition of the primordial Earth atmosphere [84]. In particular, two different compositions have been analyzed:  [35]. However, for simplicity, this process has been neglected in the present study.
Here, the cross sections sets that are used as input for the simulations are illustrated for each gas. In this work, we include also superelastic collisions from vibrational states, i.e., collisions of electrons with excited states leading to de-excitation. The cross sections for these processes have been calculated with the formula of Klein-Rosseland [86].  [89]. The elastic momentum transfer cross sections from the Biagi database is corrected to take into account the population of vibrational bending mode levels and quadrupole rotational collisions, according to Vialetto and co-authors [83]. • Carbon monoxide: Electron impact cross sections for CO are the ones from the Magboltz source code v 11.9 [87]. In [90], this set has been used for calculations of electron transport parameters in CO, with particular focus on the treatment of dipole rotational collisions. A Morse anharmonic oscillator is used for calculations of the energy of vibrational states, whereas a rigid rotator model is assumed for rotational levels [91]. • Hydrogen: Electron impact cross sections for H 2 are taken from the IST-Lisbon database of LXCat [92]. • Nitrogen: Electron impact cross sections for N 2 are taken from the IST-Lisbon database of LXCat [92].
MCF simulations have been performed for different values of the reduced electric field, namely 0.1, 10, 20 and 50 Td (1 Td = 10 −21 V m 2 ). The population of rotational and vibrational levels is calculated according to a Boltzmann distribution at 300 K. Zeroth ( f 0 ) and first order ( f 1 ) Legendre polynomials coefficients calculated from MCF are shown in Figure 6 for the highly reduced and mildly reduced mixture of gases. In both conditions, the simulated EEDFs show strong deviations from a Boltzmann equilibrium distribution, due to the several electron-molecules interactions that are taken into account. In particular, in Figure 6a, an energetic tail of the EEDF is obtained for energies above 1 eV. This is due to the presence of a deep Ramsauer-Townsend minimum in the elastic momentum transfer cross section of CH 4 , that is the dominant gas of this mixture. For energies lower than 1 eV, electron impact rotational collisions with water molecules are dominant and lead to a peak in the body of the EEDFs. In Figure 6b, a more structured EEDF is obtained with the presence of plateaux and drops in the distributions. This feature is typical of molecular gases, such as N 2 and CO 2 , that present resonant vibrational processes. In fact, those inelastic processes act as a sharp barrier that electrons have to overcome to reach higher energies [91]. Moreover, in this case, the magnitude of the first Legendre polynomial coefficient of the EVDF expansion ( f 1 ) is comparable with the zeroth order one ( f 0 ). This indicates that the popular two-term approximation that is employed in many codes (e.g., in [64,80]) would not be accurate for calculations of steady-state EVDFs.
In both cases, the calculated distribution of kinetic energy differs by orders of magnitude, in specific regions, from the Maxwell Boltzmann distribution. We therefore already have the certainty, in extreme cases of composition, that the calculation of the specific reaction rates performed with traditional methods of the primary processes of Miller experience is totally inadequate.
To summarize, calculations show that: (i) a strong deviation from the equilibrium distributions is obtained considering both a Miller-Urey-type composition of gases and a more realistic primordial Earth atmosphere; (ii) computationally efficient calculations can be performed using MCF, where great details are included in the molecular physics of vibrational and rotational states; (iii) consistent and complete sets of cross sections can be gathered from the literature to study the electron kinetics in gas mixtures relevant for the interpretation of Miller-Urey's experiment as well as its modern variants.
Future work should be focused on coupling the MCF method with other codes that calculate the non-equilibrium vibrational distribution function of the molecules, such as the Fokker-Planck method that has been introduced before.

Quantum Calculations for Elementary Processes
As emerged in Sections 2 and 3, a crucial point that may facilitate the communication between the two disciplines -prebiotic chemistry and plasma kinetics -relies on the abundance of molecular data and the knowledge of reaction rates. As we have pointed out, kinetic studies related to Miller's experience have so far largely underestimated the fact that the primary events in this kind of experiments are induced by impacts of molecules with electrons and other high energy secondary particles produced ultimately by the presence of a strong electric field. On the other hand, these aspects have always played a central role in the study of plasma reactors.
In the plasma simulation community, the problem of availability of cross sections and specific reaction rates of elementary processes emerged and has been constantly considered for decades. This has emerged also very clearly from the discussion of the cross sections used for EVDF calculations in the previous section. A constant dialogue with the elementary process community is producing new sets of data optimized for the use with specific plasma modeling methods.
As a consequence of this virtuous long-term collaboration, in the last few years extremely effective and rigorous calculation methods have been developed-based on quantum mechanical foundations of the processes-that allow one to calculate cross sections of excitation and ionization processes for chemical species of any reasonable complexity in the context of plasma sources. Currently, numerous databases relating to these processes have been developed-see for example the KIDA database for astrochemistry reactions [93], the HITRAN database for molecular spectroscopy [94], the LXCat database for plasma modeling [95] and ExoMol database for exoplanet atmospheres [96]-and constantly updated. They are managed by many active research groups in mutual communication.
It is not possible to describe in details here the different methods that have proven to be the most effective. In general, they range from relatively elementary and ready-touse techniques, based for example on improvements of the Born approximation, to very advanced techniques that solve fine details of quantum aspects of the interaction between the particles involved. In particular, these advanced methods are promising when applied to reviewing historical experiences on the origin of chemical species of biological interest; as to physical rigor, they are characterized by an appropriate versatility, that is useful not only for the simpler aspects related to primary events occurring in plasma (excitation and ionization of small molecules), but also for the analysis of relatively simpler chemical processes among those which involve species of biological interest.
Among the advanced techniques, the R-matrix method, developed in the field of astrochemistry by the Tennyson's group [97], with which one of the authors of this work has recently collaborated, is particularly interesting. In the R-matrix approach, the space around the molecule is divided into an inner and an outer region. The methodology is based on the self-consistency of two solutions of the wave equations with many particles in these two regions: (i) in the inner part, related to interactions occurring at short distances, full ab-initio quantum chemistry calculations are used, whereas (ii) in the outer region, only effective long range interactions are taken into account.
However, methods like the R-matrix can give information only for a given frozen configuration of the molecule, normally at equilibrium geometry. But this is not enough in cases where vibrations and rotations of molecules must be taken into account. This is the case of non-equilibrium systems. In order to take in account internal degrees of freedom, other techniques have to be considered, such as the Local-complex-potential (LCP) [98] and the Multichannel-Quantum-Defect-theory (MQDT) [99].
As an example of a relatively simple application, the LCP method can be used for the vibrational excitation of the diatomic CO-type (relevant for the kinetics of Millertype experiences); there are no special difficulties in applying the same technique to different chemical processes and to more complex molecules such as water, ammonia, and even simple amino acids. Figure 7 reports the potential energy curves for CO and the corresponding resonant state CO − as a function of the internuclear distance R and schematically represents the vibrational excitation reaction: where the CO molecule is supposed to be in the initial vibrational state v i = 5 and the incoming electron to have an energy of 7 eV. The system then proceeds with a sequence of two steps: (i) the electron is captured into the target, the system then jumps from the initial wave function towards the resonant state; (ii) the resonant state is metastable for some time and decays into the final vibrationally excited CO and ejects an electron. If the potential energy curves and the couplings between the target CO molecule and the resonant CO − state were known, it would be possible to determine the cross sections of the process. In recent years, motivated by the great importance of CO, CO 2 , CO + and NO molecules in astrophysics, aerospace and laser technologies, great improvements, theoretically and experimentally, in the quantitative understanding of the reactivity of these species in a weakly-ionized plasma have been achieved. Figure 8 summarizes the results for cross sections by electron impact for CO [100,101] and NO [102,103] molecules obtained with the LCP method, for CO + [104] molecular ion obtained in the MQDT framework and for experimental data for the CO 2 molecule [105].  An important aspect of these studies, and generally of all modeling studies on plasma kinetics, is the validation. This is a complex question because many details of the cross sections and kinetics revealed by these models are difficult to verify by direct experimental investigation. As far as cross sections are concerned, the most common procedure is verifying the set of cross sections as a whole, by comparing the transport quantities calculated by means of the type of calculations described in previous section. Since these calculations depend on parameters, such as the reduced electric field, which are easily modifiable in a wide range both in the model and in the experiment, the possibility of validation is concrete [106,107]. Several techniques in development, including the possibility of employing Machine Learning methods, promise to make progress in this type of procedure [108].

Conclusions and Perspectives
In the present paper, it has been shown that the methods used in the plasma modeling community reveal crucial aspects of the chemical processes that have never been taken into account in previous studies of prebiotic processes. Actually, the energy and population distribution deviate considerably from the prediction of traditional statistical mechanics, and the reaction rates are strongly dependent on such distributions. Thereby, we propose to apply the theoretical methodologies that have been developed for the study of ionized gases to the Miller-Urey experience. In this regard, in the last decades, many so-called self-consistent models were developed in the field of plasma kinetics, which can reliably calculate the kinetics of the species produced in the plasma starting from physical conditions in the reactor. Therefore, in the light of real progress in the field of prebiotic chemistry, the first step consists of determining the conditions of the spark region from the point of view of the plasma physics and chemistry: the ionization degree of the gas, the energy distribution of free electrons and the formation of the simplest reactive species that underlie the complex chemical network.
In conclusion, in this paper, we have presented a new perspective on future research on the chemical kinetics of the first active species in the primordial atmosphere and their subsequent role in the formation of the first prebiotic species. The framework of this study is the Miller-Urey experiment, its improvement and understanding, and the clues it provides into the possibility that life machinery arose as a result of an abiotic process. Two new elements of knowledge may produce a synergistic push towards further progress: the first is the awareness that the primordial atmosphere was probably not the strongly reducing mixture with essentially solar-nebula composition it was believed to be in Miller's time. The second is the development of new methods in the context of the computer modeling of the kinetics of plasmas, motivated by ecological, astrophysical and aerospace problems. The communication between the two communities of plasma kinetics and prebiotic chemistry can therefore help, in the future, to attain a better understanding and new insight into the chemical kinetics of a historical experiment, which has changed our ideas on the genesis of prebiotic molecules on the primordial Earth.