Topologically invariant reaction coordinates for simulating multistate chemical reactions.

Evaluating free energy profiles of chemical reactions in complex environments such as solvents and enzymes requires extensive sampling, which is usually performed by potential of mean force (PMF) techniques. The reliability of the sampling depends not only on the applied PMF method but also the reaction coordinate space within the dynamics is biased. In contrast to simple geometrical collective variables that depend only on the positions of the atomic coordinates of the reactants, the E(gap) reaction coordinate (the energy difference obtained by evaluating a suitable force field using reactant and product state topologies) has the unique property that it is able to take environmental effects into account leading to better convergence, a more faithful description of the transition state ensemble and therefore more accurate free energy profiles. However, E(gap) requires predefined topologies and is therefore inapplicable for multistate reactions, in which the barrier between the chemically equivalent topologies is comparable to the reaction activation barrier, because undesired "side reactions" occur. In this article, we introduce a new energy-based collective variable by generalizing the E(gap) reaction coordinate such that it becomes invariant to equivalent topologies and show that it yields more well behaved free energy profiles than simpler geometrical reaction coordinates.

Tests of an adaptive QM/MM calculation on free energy profiles of chemical reactions in solution.

We present reaction free energy calculations using the adaptive buffered force mixing quantum mechanics/molecular mechanics (bf-QM/MM) method. The bf-QM/MM method combines nonadaptive electrostatic embedding QM/MM calculations with extended and reduced QM regions to calculate accurate forces on all atoms, which can be used in free energy calculation methods that require only the forces and not the energy. We calculate the free energy profiles of two reactions in aqueous solution: the nucleophilic substitution reaction of methyl chloride with a chloride anion and the deprotonation reaction of the tyrosine side chain. We validate the bf-QM/MM method against a full QM simulation, and show that it correctly reproduces both geometrical properties and free energy profiles of the QM model, while the electrostatic embedding QM/MM method using a static QM region comprising only the solute is unable to do so. The bf-QM/MM method is not explicitly dependent on the details of the QM and MM methods, so long as it is possible to compute QM forces in a small region and MM forces in the rest of the system, as in a conventional QM/MM calculation. It is simple, with only a few parameters needed to control the QM calculation sizes, and allows (but does not require) a varying and adapting QM region which is necessary for simulating solutions.

Organic reactions in supercritical fluids: Continuous, selective and green

Supercritical fluids (SCFs) offer a wide range of opportunities as media for chemical reactions and supercritical CO2, ScCO2, is becoming increasingly important as a benign replacement for more toxic solvents.1 High pressure reactions, however, are more capital intensive than conventional low pressure processes. Therefore, supercritical fluids will only gain widespread acceptance in those areas where the fluids give real chemical advantages as well as environmental benefits. This lecture gives a brief account of the use of flow reactors for continuous reactions in supercritical fluids, particularly those of interest for the manufacture of fine chemicals.

Diffusion and chemical reaction in the porous structures of solid oxide fuel cells

The Rolls-Royce Integrated-Planar Solid Oxide Fuel Cell (IP-SOFC) consists of ceramic modules which have electrochemical cells printed on the outer surfaces. The cathodes are the outermost layer of each cell and are supplied with oxygen from air flowing over the outside of the module. The anodes are in direct contact with the ceramic structure and are supplied with fuel from internal gas channels. Natural gas is reformed into hydrogen for use by the fuel cells in a separate reformer module of similar design except that the fuel cells are replaced by a reforming catalyst layer. The performance of the modules is intrinsically linked to the behaviour of the gas flows within their porous structures. Because the porous layers are very thin, a one-dimensional flow model provides a good representation of the flow property variations between fuel channel and fuel cell or reforming catalyst. The multi-component convective-diffusive flows are simulated using a new theory of flow in porous material, the Cylindrical Pore Interpolation Model. The effects of the catalysed methane reforming and water-gas shift chemical reactions are also considered using appropriate kinetic models. It is found that the shift reaction, which is catalysed by the anode material, has certain beneficial effects on the fuel cell module performance. In the reformer module it was found that the flow resistance of the porous support structure makes it difficult to sustain a high methane conversion rate. Although the analysis is based on IP-SOFC geometry, the modelling approach and general conclusions are applicable to other types of SOFC.

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.

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.

Chemical sputtering of ta-C: Implications for the deposition of carbon nitride

The majority of attempts to synthesize the theoretically predicted superhard phase β-C3N4 have been driven towards the use of techniques which maximize both the carbon sp3 levels and the amount of nitrogen incorporated within the film. However, as yet no attempt has been made to understand the mechanism behind the resultant chemical sputter process and its obvious effect upon film growth. In this work, however, the chemical sputtering process has been investigated through the use of an as-deposited tetrahedrally bonded amorphous carbon film with a high density nitrogen plasma produced using an rf-based electron cyclotron wave resonance source. The results obtained suggested the presence of two distinct ion energy dependent regimes. The first, below 100 eV, involves the chemical sputtering of carbon from the surface, whereas the second at ion energies in excess of 100 eV exhibits a drop in sputter rate associated with the subplantation of nitrogen within the carbon matrix. Furthermore, as the sample temperature is increased there is a concomitant decrease in sputter rate suggesting that the rate is controlled by the adsorption and desorption of additional precursor species rather than the thermal desorption of CN. A simple empirical model has been developed in order to elucidate some of the primary reactions involved in the sputter process. Through the incorporation of various previously determined experimental parameters including electron temperature, ion current density, and nitrogen partial pressure the results indicated that molecular nitrogen physisorbed at the ta-C surface was the dominant precursor involved in the chemical sputter process. However, as the physisorption enthalpy of molecular nitrogen is low this suggests that activation of this molecular species takes place only through ion impact at the surface. The obtained results therefore provide important information for the modeling and growth of high density carbon nitride. © 2001 American Institute of Physics.

SU8 bio-chemical sensor microarrays

Microfabricated cantilevers have recently attracted considerable attention as novel label-free chemical and biological biosensors which translate surface reactions into nanomechanical bending motion. However these studies have primarily focused on commercially available silicon cantilevers and relatively little work has been performed on cantilevers fabricated from other materials. Polymeric materials, offer significant advantages over silicon by virtue of the low Young's modulus, ease of microfabrication and reduced cost. In this paper, we report a non-vacuum fabrication process to produce arrays of SU8 cantilevers and demonstrate their application as chemical sensors using in situ reference cantilevers. © 2006 Elsevier B.V. All rights reserved.

Cold-gas chemical vapor deposition to identify the key precursor for rapidly growing vertically-aligned single-wall and few-wall carbon nanotubes from pyrolyzed ethanol

Vertically-aligned carbon nanotubes (VA-CNTs) were rapidly grown from ethanol and their chemistry has been studied using a "cold-gas" chemical vapor deposition (CVD) method. Ethanol vapor was preheated in a furnace, cooled down and then flowed over cobalt catalysts upon ribbon-shaped substrates at 800 °C, while keeping the gas unheated. CNTs were obtained from ethanol on a sub-micrometer scale without preheating, but on a millimeter scale with preheating at 1000 °C. Acetylene was predicted to be the direct precursor by gas chromatography and gas-phase kinetic simulation, and actually led to millimeter-tall VA-CNTs without preheating when fed with hydrogen and water. There was, however a difference in CNT structure, i.e. mainly few-wall tubes from pyrolyzed ethanol and mainly single-wall tubes for unheated acetylene, and the by-products from ethanol pyrolysis possibly caused this difference. The "cold-gas" CVD, in which the gas-phase and catalytic reactions are separately controlled, allowed us to further understand CNT growth. © 2012 Elsevier Ltd. All rights reserved.