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.

The influence of substrate on the properties of Er2O3 films grown by magnetron sputtering

The structural properties and the room temperature luminescence of Er2O3 thin films deposited by RF magnetron sputtering have been studied. Films characterized by good morphological properties have been obtained by using a SiO2 interlayer between the film and the Si substrate. The evolution of the properties of the Er2O3 films due to rapid thermal annealing processes in O2 ambient performed at temperatures in the range 800-1200 °C has been investigated in details. The existence of well-defined annealing conditions (temperature of 1100 °C or higher) allowing to avoid the occurrence of extensive chemical reactions with the oxidized substrate has been demonstrated and an increase of the photoluminescence (PL) intensity by about a factor of 40 with respect to the as deposited material has been observed. The enhanced efficiency of the photon emission process has been correlated with the longer lifetime of the PL signal. The same annealing processes are less effective when Er2O3 is deposited on Si. In this latter case interfacial reactions and pit formation occur, leading to a material characterized by stronger non-radiative phenomena that limit the PL efficiency. © 2006 Elsevier B.V. All rights reserved.

Optical and structural properties of Er2O3 films grown by magnetron sputtering

The structural properties and the room temperature luminescence of Er 2O3 thin films deposited by magnetron sputtering have been studied. In spite of the well-known high reactivity of rare earth oxides towards silicon, films characterized by good morphological properties have been obtained by using a SiO2 interlayer between the film and the silicon substrate. The evolution of the properties of the Er2O3 films due to thermal annealing processes in oxygen ambient performed at temperatures in the range of 800-1200°C has been investigated in detail. The existence of well defined annealing conditions (rapid treatments at a temperature of 1100°C or higher) allowing to avoid the occurrence of extensive chemical reactions with the oxidized substrate has been demonstrated; under these conditions, the thermal process has a beneficial effect on both structural and optical properties of the film, and an increase of the photoluminescence (PL) intensity by about a factor of 40 with respect to the as-deposited material has been observed. The enhanced efficiency of the photon emission process has been correlated with the longer lifetime of the PL signal. Finally, the conditions leading to a reaction of Er2O3 with the substrate have been also identified, and evidences about the formation of silicate-like phases have been collected. © 2006 American Institute of Physics.

Single event molecular signalling for estimation and control

Cell biology is characterised by low molecule numbers and coupled stochastic chemical reactions with intrinsic noise permeating and dominating the interactions between molecules. Recent work [9] has shown that in such environments there are hard limits on the accuracy with which molecular populations can be controlled and estimated. These limits are predicated on a continuous diffusion approximation of the target molecule (although the remainder of the system is non-linear and discrete). The principal result of [9] assumes that the birth rate of the signalling species is linearly dependent on the target molecule population size. In this paper, we investigate the situation when the entire system is kept discrete, and arbitrary non-linear coupling is allowed between the target molecule and downstream signalling molecules. In this case it is possible, by relying solely on the event triggered nature of control and signalling reactions, to define non-linear reaction rate modulation schemes that achieve improved performance in certain parameter regimes. These schemes would not appear to be biologically relevant, raising the question of what are an appropriate set of assumptions for obtaining biologically meaningful results. © 2013 EUCA.

Reactive Many-Body Expansion for a Protonated Water Cluster.

We generalize the standard many-body expansion technique that is used to approximate the total energy of a molecular system to enable the treatment of chemical reactions by quantum chemical techniques. By considering all possible assignments of atoms to monomer units of the many-body expansion and associating suitable weights with each, we construct a potential energy surface that is a smooth function of the nuclear positions. We derive expressions for this reactive many-body expansion energy and describe an algorithm for its evaluation, which scales polynomially with system size, and therefore will make the method feasible for future condensed phase simulations. We demonstrate the accuracy and smoothness of the resulting potential energy surface on a molecular dynamics trajectory of the protonated water hexamer, using the Hartree-Fock method for the many-body term and Møller-Plesset theory for the low order terms of the many-body expansion.

Compatibility relationship among reaction equilibria, equivalence principle of reaction approaches, and silicon contamination in semiconductors

This paper discovers some shortcomings in the algorithm for the incorporation of Si into GaAs in the GaAs VPE process. These faults arise from neglecting a link, the compatibility relationship, in chemical thermodynamics. The meaning of said relationship is as follows: In an equilibrium complex system, each species can only contribute one and the same quantity (its equilibrium quantity) to the different equilibria of the various reactions involving it; yet even under this restriction, every equilibrium constant is satisfied, and all the reaction equilibria coexist compatibly in the system. Only by adding the relationship can the equilibrium theory for the complex system be complete. This paper also tells its position in chemical thermodynamics. Such a compatibility concept directly leads to an equivalence principle: In a complex system, a certain species can usually be simultaneously formed by many chemical reactions; when the system has reached equilibrium under fixed environmental conditions, the equilibrium quantity of said species calculated according to each chemical equation of these reactions will be equal and the various reaction approaches will be equivalent, provided that for all the reactants and all the other products of these reactions their equilibrium quantities in the system are respectively taken as corresponding knowns for the calculations, which is extremely useful for seeking a functional relation among the species' equilibrium quantities in a system (Si contamination is one of the examples). Under the guidance of those arguments, the various schools' algorithms for the Si contamination can be uniformized and simplified, and the contamination quantity relation between Si and O, two very important impurities, is found.