Modeling of H2 and H2/CH4 Moderate-Pressure Microwave Plasma Used for Diamond Deposition

One-dimensional transport models of moderate-pressure H ₂ and H ₂ /CH ₄ plasmas obtained in a diamond deposition microwave reactor are presented. These models describe the plasma as a thermochemically nonequilibrium flow with three different energy modes. The solution of the one-dimensional plasma transport equations enabled the estimation of plasma species concentrations and temperatures on the axis of the reactor. As far as pure H ₂ plasmas are concerned, results showed that the model predictions of gas and vibration temperatures are in good agreement with experimental measurements. The model also yields a relatively good qualitative prediction of the variations of H-atom mole fraction with the power density absorbed by the plasma. The results obtained for H ₂ /CH ₄ discharges showed that the model prediction on the variations of H-atom mole fraction with methane percentage in the discharge is in good qualitative agreement with experimental results. They also showed that methane is rapidly converted to acetylene before reaching the discharge zone. The concentrations of neutral hydrocarbon species in the reactor are mainly governed by thermal chemistry. The addition of methane strongly affects the ionization kinetics of the plasma. Three major ions are generally obtained in H ₂ /CH ₄ plasmas: C ₂ H ₂ ⁺ , C ₂ H ₃ ⁺ , and C ₂ H ₅ ⁺ . The relative predominance of these ions depends on the considered plasma region and on the discharge conditions. The ionic species concentrations are also mainly governed by chemistry, except very near the substrate surface. Finally the use of this transport model along with the surface chemistry model of Goodwin ⁽¹⁾ enabled us to estimate the diamond growth rate for several discharge conditions.     

Experimental and Theoretical Studies of a Pulsed Microwave Excited Ar/CF4 Plasma

The present work deals with a pulsed microwave discharge in an Ar/CF ₄ gas mixture under a low pressure (1–10 mbar). The discharge chamber developed has a cylindrical geometry with a coupling window alternatively made of quartz or alumina. The setup allows one to investigate the plasma–wall interactions (here etching of the quartz window) and the ignition process of the pulsed microwave plasma. Microwave pulses with a duration of 50–200 s and repetition rate between 1 and 10 kHz are typical for the experiments. The space-time behavior of the fluorine number density in the discharge has been investigated experimentally by optical actinometry. The discharge kinetics is modeled using electron-transport parameters and rate coefficients derived from solutions of the Boltzmann equation. Together with the solution of the continuity and electron balance equations and the rate equations describing the production of CF _x (x=2, 3, 4) radicals and F atoms, a good agreement between experimental and theoretical data can be achieved.     

Classification of auxin plant hormones by interaction property similarity indices

Although auxins were the first type of plant hormone to be identified, little is known about the molecular mechanism of this important class of plant hormones. We present a classification of a set of about 50 compounds with measured auxin activities, according to their interaction properties. Four classes of compounds were defined: strongly active, weakly active with weak antiauxin behaviour, inactive and inhibitory. All compounds were modeled in two low-energy conformations, P and T, so as to obtain the best match to the planar and tilted conformations, respectively, of indole 3-acetic acid. Each set of conformers was superimposed separately using several different alignment schemes. Molecular interaction energy fields were computed for each molecule with five different chemical probes and then compared by computing similarity indices. Similarity analysis showed that the classes are on average distinguishable, with better differentiation achieved for the T conformers than the P conformers. This indicates that the T conformation might be the active one. Further, a screening was developed which could distinguish compounds with auxin activity from inactive compounds and most antiauxins using the T conformers. The classifications rationalize ambiguities in activity data found in the literature and should be of value in predicting the activities of new plant growth substances and herbicides.     

Reparametrization of models of non-stationary chemical kinetics as a method of investigation of the inverse problem

A new approach to solve the inverse problem of the non-stationary chemical kinetics by a reparametrization of the initial model is suggested. This approach permits to obtain all solutions of IP in the case of a finite number or some parametric functions in the case of the infinite number of solutions.     

Prediction of 13C NMR chemical shifts in substituted naphthalenes

The prediction of the ¹³C NMR signals for derivatives of naphthalene has been investigated using mathematical modeling techniques. Two empirical multiple regression models which utilize the field, resonance, and Charton's steric parameters together with molar refractivity were developed, one for α- and the other for β-substituted naphthalene derivatives. In the α case the model had a correlation coefficient of observed versus predicted line positions of r = 0.973 with a standard deviation of 2.2 ppm while in the β case r = 0.979 with the standard deviation being 2.3 ppm. The database consisted of 3152 signals from 394 naphthalene derivatives. We also report the use of the Taft steric parameter in place of the Charton steric parameter in the above- mentioned prediction equations. Received: 19 June 1998 / Accepted: 20 October 1998 / Published online: 16 March 1999     

Finite difference modeling of the gas transport process in glassy polymers

Finite difference modeling has been used to predict the results of gas transport experiments for a concentration-dependent diffusion coefficient. Experiments on the transport of CO₂ in poly(ethylene terephthalate) and poly(ethylene naphthalate) had previously shown a difference between the effective diffusion coefficients for absorption and desorption runs of a double-sided experiment, but this effect had not been seen for single-sided experiments. The finite difference calculations show that such results are to be expected, and the parameters included in the models that attempt to describe the diffusion process in glassy polymers, such as the dual-mode model, and which lead to concentration-dependent diffusion coefficients, can be found by fitting the experimental data for the double-sided experiment using finite difference modeling. The dependence of the effective diffusion coefficient on pressure for the single-sided experiment can be correctly predicted using results from the double-sided experiment for an identical sample. © 1996 John Wiley & Sons, Inc.     

A Three-Dimensional Two-Phase Model for Thermal Plasma Chemical Vapor Deposition with Liquid Feedstock Injection

In this paper, a comprehensive model for thermal plasma chemical vapor deposition (TPCVD) with liquid feedstock injection is documented. The gas flow is assumed to be steady, of a single temperature. Radiation and charged species contributions are excluded, but extensive homogeneous and heterogeneous chemistry is included. The liquid phase is traced by considering individual droplets. Discussion on the model's application to diamond production from acetone in a hydrogen–argon plasma is included. The major conclusions are: (1) Liquid injection possesses a capability to deliver the hydrocarbon precursor directly onto the deposition target. (2) For the case of complete evaporation of the droplet before reaching the substrate, the deposition rate is similar to that obtained with gaseous precursors. (3) The computational results compare well with experimental data. The modeling results can be used to optimize the injection parameters with regard to the deposition rate.     

The computer program LUDI: A new method for the de novo design of enzyme inhibitors

Summary A new computer program is described, which positions small molecules into clefts of protein structures (e.g. an active site of an enzyme) in such a way that hydrogen bonds can be formed with the enzyme and hydrophobic pockets are filled with hydrophobic groups. The program works in three steps. First it calculates interaction sites, which are discrete positions in space suitable to form hydrogen bonds or to fill a hydrophobic pocket. The interaction sites are derived from distributions of nonbonded contacts generated by a search through the Cambridge Structural Database. An alternative route to generate the interaction sites is the use of rules. The second step is the fit of molecular fragments onto the interaction sites. Currently we use a library of 600 fragments for the fitting. The final step in the present program is the connection of some or all of the fitted fragments to a single molecule. This is done by bridge fragments. Applications are presented for the crystal packing of benzoic acid and the enzymes dihydrofolate reductase and trypsin.     

Mechanism of the gas-phase oxidation of carbon disulfide at elevated temperatures (the C-S-O system)

Machine modeling was used to investigate the kinetic mechanism of the gas-phase oxidation of carbon disulfide at elevated temperatures. The scheme of the process examined includes 70 elementary reactions. On the basis of the kinetic data on these reactions available in the literature the rates of the oxidation process, the concentrations of end products and intermediate particles were calculated as a function of the time under various conditions (composition, pressure, temperature). Comparison of the results obtained with the experimental data showed satisfactory agreement and also revealed some important elementary reactions.N. N. Semenov Institute of Chemical Physics, Russian Academy of Sciences, Moscow 117977. Translated from Izvestiya Akademii Nauk, Seriya Khimicheskaya, No. 4, pp. 812–822, April, 1992.