Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Numerical Investigation On Detonation Cell Evolution In A Channel With Area-Changing Cross Section

The two-dimensional cellular detonation propagating in a channel with area-changing cross section was numerically simulated with the dispersion-controlled dissipative scheme and a detailed chemical reaction model. Effects of the flow expansion and compression on the cellular detonation cell were investigated to illustrate the mechanism of the transverse wave development and the cellular detonation cell evolution. By examining gas composition variations behind the leading shock, the chemical reaction rate, the reaction zone length, and thermodynamic parameters, two kinds of the abnormal detonation waves were identified. To explore their development mechanism, chemical reactions, reflected shocks and rarefaction waves were discussed, which interact with each other and affect the cellular detonation in different ways.

Optimization of tried wave-function in quantum monte-carlo method

A method for optimizing tried wave functions in quantum Monte Carlo method has been found and used to calculate the energies of molecules, such as H-2, Li-2, H-3+, H-3 and H-4. Good results were obtained.

化学反应和球面聚心爆轰波的相互作用

考察球面爆轰波聚心传播过程中，波阵面附近压力和温度不断升高引起化学反应进程的改变；对比氢氧可燃气体与氮气的数值模拟结果，分析化学反应对波面温度和压力的影响，从而考察Zeldovich理论预测聚心爆轰波后参数的精确性．数值结果表明，爆轰波聚心传播初期，放热的燃烧反应对波后热力学参数起主导作用；传播后期，波阵面趋近于对称中心时，吸热的气体解离反应变得非常活跃，解离反应对后期的汇聚压力影响不大，但会在很大程度上限制汇聚温度的升高．

Control of Large-size SiC Growth Process Using Modeling tool

Large size bulk silicon carbide (SiC) crystals are commonly grown by the physical vapor transport (PVT) method. The PVT growth of SiC crystals involves sublimation and condensation, chemical reactions, stoichiometry, mass transport, induced thermal stress, as well as defect and micropipes generation and propagation. The quality and polytype of as-grown SiC crystals are related to the temperature distribution inside the growth chamber during the growth process, it is critical to predict the temperature distribution from the measured temperatures outside the crucible by pyrometers. A radio-frequency induction-heating furnace was used for the growth of large-size SiC crystals by the PVT method in the present study. Modeling and simulation have been used to develop the SiC growth process and to improve the SiC crystal quality. Parameters such as the temperature measured at the top of crucible, temperature measured at the bottom of the crucible, and inert gas pressure are used to control the SiC growth process. By measuring the temperatures at the top and bottom of the crucible, the temperatures inside the crucible were predicted with the help of modeling tool. SiC crystals of 6H polytype were obtained and characterized by the Raman scattering spectroscopy and SEM, and crystals of few millimeter size grown inside the crucible were found without micropipes. Expansion of the crystals were also performed with the help of modeling and simulation.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

碳氢燃料超音速燃烧的数值研究

本文对CH_4横向喷入超音速高温空气主流的二维流场进行了数值模拟，对CH_4-O_2的六方程反应模式进行了检验，借以对混合与燃烧过程的现象与机理加以研究，并与单方程反应模型及H_2横向喷射的情况进行了对比。得到了较为理想的结果。本文采用二维雷诺平均全N-S方程进行计算，采用热完全气体模型，用Baldwin-Lomax代数涡粘性湍流模型来模拟湍流效应。假定N_2不参加反应，CH_4-O_2反应机制选取六个基本反应，以及九个组元O、O_2、CH_3、CH_4、OH、CHO、CH_2O、CO和CO_2，应用空间二阶精度Harten-Yee隐式TVD格式，采用化学源项点隐的全隐方法数值求解。本文对多种超音速空气主流及边界条件的工部进行了数值模拟，并对其流场进行了分析。针对工作中出现的问题，提出了对下一步工作的展望。

Progress towards the total synthesis of a spermine-conjugated dynemicin analog

Progress towards the synthesis of the spermine-conjugated Dynemicin analog 4 is described. The synthetic route starts with the Michael addition of menthyl acetoacetate to trans-ethyl crotonate followed by a Dieckman condensation to form the cyclohexanedione 14 which, through a series of chemical reactions, is transformed into the quinone imine 6. Key features in the route include the Suzuki coupling reaction of the aryl boronic acid 11 and the enol triflate 12, thermal deprotection/internal amidation of the biaryl 19, cis addition of the (Z)-enediyne 33 to the quinoline 25, intramolecular acetylide addition to a carbonyl within the ketone 29, and an addition/elimination of the cyanophthalide to the quinone imine 6 to form the anthraquinone 36 utilizing the Kraus and Sugimoto methodology.

Chelation-enforced metal-arene interactions : insights into substrate binding and catalysis by late transition metal complexes

Understanding and catalyzing chemical reactions requiring multiple electron transfers is an endeavor relevant to many outstanding challenges in the field of chemistry. To study multi-electron reactions, a terphenyl diphosphine framework was designed to support one or more metals in multiple redox states via stabilizing interactions with the central arene of the terphenyl backbone. A variety of unusual compounds and reactions and their relevance toward prominent research efforts in chemistry are the subject of this dissertation.

Chapter 2 introduces the para-terphenyl diphosphine framework and its coordination chemistry with group 10 transition metal centers. Both mononuclear and dinuclear compounds are characterized. In many cases, the metal center(s) are stabilized by the terphenyl central arene. These metal–arene interactions are characterized both statically, in the solid state, and fluxionally, in solution. As a proof-of-principle, a dinickel framework is shown to span multiple redox states, showing that multielectron chemistry can be supported by the coordinatively flexible terphenyl diphosphine.

Chapter 3 presents reactivity of the terphenyl diphosphine when bound to a metal center. Because of the dearomatizing effect of the metal center, the central arene of the ligand is susceptible to reactions that do not normally affect arenes. In particular, Ni-to-arene H-transfer and arene dihydrogenation reactions are presented. Additionally, evidence for reversibility of the Ni-to-arene H-transfer is discussed.

Chapter 4 expands beyond the chelated metal-arene interactions of the previous chapters. A dipalladium(I) terphenyl diphosphine framework is used to bind a variety of exogenous organic ligands including arenes, dienes, heteroarenes, thioethers, and anionic ligands. The compounds are structurally characterized, and many ligands exhibit unprecedented bindng modes across two metal centers. The relative binding affinities are evaluated spectroscopically, and equilibrium binding constants for the examined ligands are determined to span over 13 orders of magnitude. As an application of this framework, mild hydrogenation conditions of bound thiophene are presented.

Chapter 5 studies nickel-mediated C–O bond cleavage of aryl alkyl ethers, a transformation with emerging applications in fields such as lignin biofuels and organic methodology. Other group members have shown the mechanism of C–O bond cleavage of an aryl methyl ether incorporated into a meta-terphenyl diphosphine framework to proceed through β-H elimination of an alkoxide. First, the electronic selectivity of the model system is examined computationally and compared with catalytic systems. The lessons learned from the model system are then applied to isotopic labeling studies for catalytic aryl alkyl ether cleavage under dihydrogen. Results from selective deuteration experiments and mass spectrometry draw a clear analogy between the mechanisms of the model and catalytic systems that does not require dihydrogen for C–O bond cleavage, although dihydrogen is proposed to play a role in catalyst activation and catalytic turnover.

Appendix A presents initial efforts toward heterodinuclear complexes as models for CO dehydrogenase and Fischer Tropsch chemistry. A catechol-incorporating terphenyl diphosphine is reported, and metal complexes thereof are discussed.

Appendix B highlights some structurally characterized terphenyl diphosphine complexes that either do not thematically belong in the research chapters or proved to be difficult to reproduce. These compounds show unusual coordination modes of the terphenyl diphosphine from which other researchers may glean insights.

Secondary Organic Aerosol Composition Studies Using Mass Spectrometry

Trace volatile organic compounds emitted by biogenic and anthropogenic sources into the atmosphere can undergo extensive photooxidation to form species with lower volatility. By equilibrium partitioning or reactive uptake, these compounds can nucleate into new aerosol particles or deposit onto already-existing particles to form secondary organic aerosol (SOA). SOA and other atmospheric particulate matter have measurable effects on global climate and public health, making understanding SOA formation a needed field of scientific inquiry. SOA formation can be done in a laboratory setting, using an environmental chamber; under these controlled conditions it is possible to generate SOA from a single parent compound and study the chemical composition of the gas and particle phases. By studying the SOA composition, it is possible to gain understanding of the chemical reactions that occur in the gas phase and particle phase, and identify potential heterogeneous processes that occur at the surface of SOA particles. In this thesis, mass spectrometric methods are used to identify qualitatively and qualitatively the chemical components of SOA derived from the photooxidation of important anthropogenic volatile organic compounds that are associated with gasoline and diesel fuels and industrial activity (C12 alkanes, toluene, and o-, m-, and p-cresols). The conditions under which SOA was generated in each system were varied to explore the effect of NO_x and inorganic seed composition on SOA chemical composition. The structure of the parent alkane was varied to investigate the effect on the functionalization and fragmentation of the resulting oxidation products. Relative humidity was varied in the alkane system as well to measure the effect of increased particle-phase water on condensed-phase reactions. In all systems, oligomeric species, resulting potentially from particle-phase and heterogeneous processes, were identified. Imines produced by reactions between (NH₄)₂SO₄ seed and carbonyl compounds were identified in all systems. Multigenerational photochemistry producing low- and extremely low-volatility organic compounds (LVOC and ELVOC) was reflected strongly in the particle-phase composition as well.

Forced convection film boiling of subcooled water around a sphere

An experimental investigation was made of forced convection film boiling of subcooled water around a sphere at atmospheric pressure. The water was sufficiently cool that the vapor condensed before leaving the film with the result that no vapor bubbles left the film. The experimental runs were made using inductively heated spheres at temperatures above 740°C. and using inlet water temperatures between 15°C. and 27°C. The spheres used had diameters of 1/2 inch, 9/16 inch, and 3/8 inch and were supported by the liquid flow. Reynolds numbers between 60 and 700 were used.

Analysis of the collected non-condensables indicated that oxygen and nitrogen dissolved in the water accumulated within the vapor film and that hetrogeneous chemical reactions occurred at the sphere surface. An iron-steam reaction resulted in more than 20% by volume hydrogen in the film at wall temperatures above 900°C. At temperatures near 1100°C. more than 80% by volume of the film was composed of hydrogen. It was found that gold plating of the sphere could eliminate this reaction.

Material and energy balances were used to derive equations which may be used to predict the overall average heat transfer coefficients for subcooled film boiling around a sphere. These equations include the effect of dissolved gases in the water. Equations also were derived which may be used to predict the composition of the film for cases in which an equilibrium exists between the dissolved gases and the gases in the film.

The derived equations were compared to the experimental results. It was found that a correlation existed between the Nusselt number for heat transfer from the vapor-liquid interface into the liquid and the Reynolds number, liquid Prandtl number product. In addition, it was found that the percentage of dissolved oxygen removed during the film boiling could be predicted to within 10%.

1. Ultrasonic studies of binary liquid structure in the critical region. Theory and experiment for the 2,6-lutidine/water system. 2. Hartree-Fock calculations of electric polarizabilities of some simple atoms and molecules, and their practicality. 3. Calculation of vibrational transition probabilities in collinear atom-diatom and diatom-diatom collisions with Lennard-Jones interaction

Part 1. Many interesting visual and mechanical phenomena occur in the critical region of fluids, both for the gas-liquid and liquid-liquid transitions. The precise thermodynamic and transport behavior here has some broad consequences for the molecular theory of liquids. Previous studies in this laboratory on a liquid-liquid critical mixture via ultrasonics supported a basically classical analysis of fluid behavior by M. Fixman (e. g., the free energy is assumed analytic in intensive variables in the thermodynamics)--at least when the fluid is not too close to critical. A breakdown in classical concepts is evidenced close to critical, in some well-defined ways. We have studied herein a liquid-liquid critical system of complementary nature (possessing a lower critical mixing or consolute temperature) to all previous mixtures, to look for new qualitative critical behavior. We did not find such new behavior in the ultrasonic absorption ascribable to the critical fluctuations, but we did find extra absorption due to chemical processes (yet these are related to the mixing behavior generating the lower consolute point). We rederived, corrected, and extended Fixman's analysis to interpret our experimental results in these more complex circumstances. The entire account of theory and experiment is prefaced by an extensive introduction recounting the general status of liquid state theory. The introduction provides a context for our present work, and also points out problems deserving attention. Interest in these problems was stimulated by this work but also by work in Part 3.

Part 2. Among variational theories of electronic structure, the Hartree-Fock theory has proved particularly valuable for a practical understanding of such properties as chemical binding, electric multipole moments, and X-ray scattering intensity. It also provides the most tractable method of calculating first-order properties under external or internal one-electron perturbations, either developed explicitly in orders of perturbation theory or in the fully self-consistent method. The accuracy and consistency of first-order properties are poorer than those of zero-order properties, but this is most often due to the use of explicit approximations in solving the perturbed equations, or to inadequacy of the variational basis in size or composition. We have calculated the electric polarizabilities of H₂, He, Li, Be, LiH, and N₂ by Hartree-Fock theory, using exact perturbation theory or the fully self-consistent method, as dictated by convenience. By careful studies on total basis set composition, we obtained good approximations to limiting Hartree-Fock values of polarizabilities with bases of reasonable size. The values for all species, and for each direction in the molecular cases, are within 8% of experiment, or of best theoretical values in the absence of the former. Our results support the use of unadorned Hartree-Pock theory for static polarizabilities needed in interpreting electron-molecule scattering data, collision-induced light scattering experiments, and other phenomena involving experimentally inaccessible polarizabilities.

Part 3. Numerical integration of the close-coupled scattering equations has been carried out to obtain vibrational transition probabilities for some models of the electronically adiabatic H₂-H₂ collision. All the models use a Lennard-Jones interaction potential between nearest atoms in the collision partners. We have analyzed the results for some insight into the vibrational excitation process in its dependence on the energy of collision, the nature of the vibrational binding potential, and other factors. We conclude also that replacement of earlier, simpler models of the interaction potential by the Lennard-Jones form adds very little realism for all the complication it introduces. A brief introduction precedes the presentation of our work and places it in the context of attempts to understand the collisional activation process in chemical reactions as well as some other chemical dynamics.

Evaluating boundary dependent errors in QM/MM simulations.

Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations provide a powerful tool for studying chemical reactions, especially in complex biochemical systems. In most works to date, the quantum region is kept fixed throughout the simulation and is defined in an ad hoc way based on chemical intuition and available computational resources. The simulation errors associated with a given choice of the quantum region are, however, rarely assessed in a systematic manner. Here we study the dependence of two relevant quantities on the QM region size: the force error at the center of the QM region and the free energy of a proton transfer reaction. Taking lysozyme as our model system, we find that in an apolar region the average force error rapidly decreases with increasing QM region size. In contrast, the average force error at the polar active site is considerably higher, exhibits large oscillations and decreases more slowly, and may not fall below acceptable limits even for a quantum region radius of 9.0 A. Although computation of free energies could only be afforded until 6.0 A, results were found to change considerably within these limits. These errors demonstrate that the results of QM/MM calculations are heavily affected by the definition of the QM region (not only its size), and a convergence test is proposed to be a part of setting up QM/MM simulations.

Carboxilesterase na adaptação do metabolismo de lipídios do pacu à hipoxia

Os peixes são vertebrados que vivem em vários habitats. Adaptações de sua fisiologia e de sua bioquímica são estudadas para entender como conseguem sobreviver aos desafios presentes em cada ambiente. A diminuição da concentração de oxigênio dissolvido na água é um fenômeno natural cíclico em águas do Pantanal e da Amazônia. A crescente poluição antrópica dos rios da Amazônia e do Pantanal pode, somada à hipoxia, oferecer ameaça à permanência de muitas espécies de peixes. Ao estudarmos a carboxilesterase (CarbE), enzima importante para a biotransformação de xenobióticos, observamos que sua atividade diminuía em plasma de pacu, um peixe típico do Pantanal, mantido em hipoxia durante 42 horas. As carboxilesterases participam de inúmeras reações químicas no organismo, incluindo uma transesterificação capaz de produzir ésteres etílicos tóxicos de ácido graxo (FAEE, fatty acid ethyl esters) a partir de etanol e ésteres de ácidos graxos. Já que a diminuição da atividade de CarbE poderia ser uma vantagem, pois que o etanol (um produto da glicólise em peixes sob hipoxia) seria menos esterificado, resolvemos saber mais sobre a bioquímica da CarbE do plasma e do fígado de pacus. Os pacus foram submetidos à hipoxia por diminuição da concentração até 0,5 mg O2/L por meio de borbulhamento de nitrogênio na água. Os animais ficaram nessas condições por 42 horas, quando então coletamos sangue e retiramos seus fígados. A atividade de CarbE ensaiada foi 50% menor no soro e 25% menor nos microssomos de fígado se comparada com a de peixes sob 6 mg O2/L. A CarbE isolada do soro dos pacus em normoxia possui massa molecular relativa de 56.000. A eletroforese em gel desnaturante com a fração purificada rendeu três bandas, mas o gel nativo apresentou só duas bandas com atividades sobre α-naftil acetato. Inferimos que mais do que uma isoforma da enzima está no plasma dos pacus. A CarbE isolada do soro não possui atividade de lipase sobre o Tween 20, mas os microssomos dos fígados dos animais em normoxia e hipoxia possuem. No entanto, a CarbE possui atividade sobre acil-CoA assim como o microssomo de fígado. A enzima pura apresenta Vmáx três vezes maior e uma KM quatro vezes maior do que a atividade presente nos microssomos. Além disso, os microssomos do fígado dos pacus em normoxia podem hidrolisar acil-CoA com o dobro da velocidade daqueles em hipoxia. A atividade clássica de CarbE do soro foi inibida por 4-HNE a 4 mM. A atividade de acil-CoA dos microssomos de fígado também foi sensível a 4-HNE a 4 mM, enquanto 2 mM de 4-HNE inibiu metade desta atividade do soro.

Thermal evolution of Er silicate thin films grown by rf magnetron sputtering

Stoichiometric Er silicate thin films, monosilicate (Er2SiO 5) and disilicate (Er2Si2O7), have been grown on c-Si substrates by rf magnetron sputtering. The influence of annealing temperature in the range 1000-1200 °C in oxidizing ambient (O 2) on the structural and optical properties has been studied. In spite of the known reactivity of rare earth silicates towards silicon, Rutherford backscattering spectrometry shows that undesired chemical reactions between the film and the substrate can be strongly limited by using rapid thermal treatments. Monosilicate and disilicate films crystallize at 1100 and 1200 °C, respectively, as shown by x-ray diffraction analysis; the crystalline structures have been identified in both cases. Moreover, photoluminescence (PL) measurements have demonstrated that the highest PL intensity is obtained for Er2Si2O7 film annealed at 1200 °C. In fact, this treatment allows us to reduce the defect density in the film, in particular by saturating oxygen vacancies, as also confirmed by the increase of the lifetime of the PL signal. © 2008 IOP Publishing Ltd.