Supporting Information
|
|
- Junior Greer
- 5 years ago
- Views:
Transcription
1 Prebiotic NH 3 formation: Insights from simulations András Stirling 1, Tamás Rozgonyi 2, Matthias Krack 3, Marco Bernasconi 4 1 Institute of Organic Chemistry, Research Centre for Natural Sciences of the Hungarian Academy of Sciences, POB 286, Budapest, 1519, Hungary 2 Institute of Materials and Environmental Chemistry, Research Centre for Natural Sciences of the Hungarian Academy of Sciences, POB 286, Budapest, 1519, Hungary 3 Laboratory for Reactor Physics and Systems Behaviour, Paul Scherrer Institute, 5232 Villigen PSI, Switzerland 4 Department of Materials Science, University of Milano-Bicocca, Via R. Cozzi 55, Milano, Italy Supporting Information S1
2 Contents 1 Optimization of NO 3 and NO 2 on pyrite surfaces S3 2 NEB calculations S3 3 Maximally localized Wannier functions S6 4 Metadynamics simulations S7 5 Error Sources S8 5.1 Issues of the quantum chemical method S8 5.2 Issues of the free energy samplings S8 6 Representative snapshots of transition states S11 7 Free energy surfaces S12 References S14 S2
3 1 Optimization of NO 3 and NO 2 on pyrite surfaces We have identified a large number of possible minima on the potential energy surfaces of the adsorption of NO 3 and NO 2 on perfect and defective pyrite surfaces without water. For the optimization we have used the same quantum chemical setup as in the MD simulations. The main conclusions of the optimizations are the following: NO 3 : on the perfect and defective surfaces the optimal molecular adsorption is an O,O bidentate coordination mode to two neighbour iron sites. Monodentate modes are not found. Interestingly, the reverse is true in the presence of hot-pressurized water where an O-monodentate mode became the preferred coordination. This can be explained by the favorable interaction of the ion with the H 2 O molecules of the water layer. NO 2 : on the perfect surface the most stable adsorption form is the O,O bidentate coordination mode to two neighbour iron sites. However, the N-monodentate coordination is only 2.5 kcal/mol less stable. The presence of water has an important effect: beside these two coordination forms, the O-monodentate form also becomes accessible, because this provides the easiest contact with the water layer. On the defective surface the most stable form is the N,O-bidentate coordination on the defective iron sites. In the presence of water both this coordination form and the N-monodentate form could be observed during the unbiased simulations. 2 NEB calculations Climbing image nudged elastic band (CI-NEB) calculations 1 were performed to investigate NO x dissociation on dry pyrite surfaces using the same quantum chemical method and parameters and the same periodic pyrite slab as in the MD simulations but without water. The spring parameter was 0.6 au for all the calculations. The convergence criteria on the band were the following: maximum value of displacement was au, maximum value of forces was au, root mean square (rms) of displacement was au and rms of forces was au. In case of multistep dissociation processes each step was optimized separately. For each individual step (i.e., featuring only one transition state (TS)) 8 or 12 images were used and the two end-point configurations for each steps were freely optimized. Figures S1 and S2 depict the initial, the transition and the final states of the dissociation reactions for the perfect and the defective surfaces, respectively. The dissociation of NO 3 turned out to be a two-step process on the perfect surface and a threestep process on the defective surface. The first step is the breaking of one of the N-O bonds followed by the barrierless transfer of this O atom to a neighboring sulfur site. The second step of the process is the coordination of the dangling O atom to the surface iron site which has a much lower energy barrier. The dissociation of NO 3 on the defective surface is similar but the O atom fills the vacancy site. The reaction energy is considerably larger and the overall barrier is considerably smaller for the NO 3 dissociation on the defective than on the perfect surface (see Tables 1 and 2 ). The dissociation of NO 2 on both perfect and defective surfaces is a single step process. The S3
4 Figure S1: Initial, transition state and final geometries of the oxygen transfer reactions occuring on the perfect FeS 2 surface as obtained from NEB calculations. First row: NO 3 NO 2 + O, the transition state with the highest barrier is shown. ; second row: NO 2 NO + O ; third row: NO N + O. For color code see Fig. 1 of the article. Table 1: Reaction energies and barrier heights for the dissociation of NO x surface. Energies are in kcal/mol. species on perfect pyrite process reaction energy barrier height NO 3 NO 2 + O NO 2 NO + O NO N + O dissociation of NO 2 on perfect surface starts with the breaking of one of the N-O bonds and then the O atom is transferred to a neighboring sulfur site, while the product NO reorients to an Fe-NO configuration without barrier resulting in an overall slightly endotherm process. The dissociation of NO 2 from its most stable N,O bidentate coordination on the defective surface is slightly exotherm whereas a significant barrier arises for the N-O bond breaking which is followed by the movement of the dissociated O atom into the vacancy. The dissociation of the NO on perfect surface is a single step process. The TS configuration is an NO coordinated with the N and O to surface Fe and S atoms, respectively. The stable S4
5 Figure S2: Initial, transition state and final geometries of the oxygen transfer reactions occuring on the defective FeS2 surface as obtained from NEB calculations. First row: NO 3 NO2 + O, the transition state with the highest barrier is shown; second row: NO 2 NO + O ; third row: NO N + O, occurring in the presence of a vacancy filled by an O atom. For color code see Fig. 1 of the article. Table 2: Reaction energies and barrier heights for the dissociation of NO x species on defective pyrite surface. Energies are in kcal/mol. process NO 3 NO 2 NO 2 NO NO N reaction energy barrier height +O O O dissociated form is then a surface iron-nitride species and a neighbor monoxidized surface sulfur site. Similarly, on the defective surface the dissociating O atom forms a bond with a neighbor sulfur atom. The dissociation of NO2 and NO on perfect FeS2 surface has been evaluated also by Sacchi et al.2 using Linear Synchronous Transit and Quadratic Synchronous Transit methods.3 They obtained 46.1 and 61.3 kcal/mol activation barriers for the O abstraction from NO2 or from NO, respectively, in reasonable agreement with our results. S5
6 3 Maximally localized Wannier functions In order to have a chemically meaningful picture about the valance electron distributions we have localized the Kohn-Sham orbitals according to the maximal localization criteria by Marzari and Vanderbilt 4 and analyzed the centers of these Wannier orbitals. These orbitals are analogous of the localized molecular orbitals of clusters in the traditional quantum chemistry. Figure S3 illustrates this analysis for four selected examples of the adsorbed species. Note that our methodology treats explicitely only the outer electrons whereas the core electrons and the nuclei are replaced by pseudopotentials. In particular, we use an iron pseudopotential which corresponds to an extended valence state: 3s 2 3p 6 4s 2 3d 6. For the other atoms the number of valence electrons are: 6 for O and S, 5 for N and 1 for H. Panel A shows an NO 3 adsorbed on a surface iron site. The iron is in its 2+ state therefore it has 14 valence and semicore electrons, as shown by the 7 WCs. The NO 3 unit features 12 WCs, which clearly indicate a nitrate anion (a neutral NO 3 would have 23 valence electrons, whereas the 12 WCs represent 24 electrons). On panel B an NO anion adsorbed on an iron site is depicted. Again the iron is in its Fe(II) state with seven WCs. The NO unit features 6 WCs which indicate a NO charge state (neutral NO has 11 valence electrons). Panel C shows an Fe(VI)-nitrido species. 5 WCs surround the iron ion, 3 WCs constitute the triple Fe-N bond and a Wannier center rests on the nitrogen. This is the bonding mechanism given for this high-valent iron state. 5, 6 The Fe-N unit has a positive charge (in neutral state it would have 14+5=19 electrons, whereas the 9 WCs represent 18e charge). On panel D an adsorbed NH 2 O species is seen. The iron ion is surrounded by 7 WCs indicating Fe(II) state. Another 7 WCs can be found around the NH 2 O species, which indicates a negatively charged state (neutral NH 2 O would have 13 valence electrons). A B C D Figure S3: Centers (green balls) of the maximally localized Wannier functions. For the further colors see Fig. 1 of the article. A: NO 3 adsorbed on an iron site. B: NO anion adsorbed on an iron site. C: iron(vi)-nitrido state. D: NH 2 O adsorbed on iron. S6
7 4 Metadynamics simulations We have selected various functions of coordination numbers 7 (CNs) as CVs. Table 3 summarizes these CN-s and their functions as CVs. Table 3: Coordination numbers applied in the present study and their action. Distances are in Å. A B r c p q role of the CV CN1 N O NOx N-O bond breaking CN2 N O NOx N-O bond breaking CN3 Fe surf, S surf O NOx O interaction with the surface CN4 Fe surf, S surf O NOx O interaction with the surface CN5 Fe surf, S surf O NOx, N NO 2 adsorption, O transfer CN6 N H water proton interaction with N CN7 Fe surf, S surf O water, water adsorption on surface CN8 O water H water O-H bond breaking functions of CNs CV9 CN4-CN2 O transfer from NO x to the surface CV10 CN8-CN7 O-H bond breaking upon adsorption Most of the CVs involves many degrees of freedom which implies that they are truly collective in nature and the corresponding biasing potential makes a large variety of reaction channels accessible. We have used metadynamics in its direct approach without introducing auxiliary variables for the CVs. In exploring a reaction route we have first performed an exploratory run using Gaussian biasing potentials with a height of k B T and width of 0.1 (in coordination number units). Then in subsequent refinements we started the simulations from a situation where the first free energy well was filled roughly half-way and continued the metadynamics exploration with lower Gaussian hills. A similar strategy has been successfully employed in previous studies on the simulation of the peptide synthesis at the FeS 2 surface. 8, 9 With this strategy the expected error of the metadynamics sampling is of the order of k B T (see further discussion below). Inspection of the temperature of the individual atoms has shown that thermal equilibration during the metadynamics runs was sufficient in most of the cases and that in less efficient cases refinement of the simulations enhanced the temperature homogeneity. In the present study, we have focused on the forward paths. This implies that a full free energy profile cannot be constructed. For the barrier of a given step several runs have been performed and the average barriers are presented. We did not seek to obtain relative stabilities of the intermediates and backward free energy barriers. As shown in the main article the thermodynamic force is very large toward the NH 3 formation which indicates that the overall exergonic reaction can be analyzed and interpreted by the barriers obtained from metadynamics. S7
8 5 Error Sources 5.1 Issues of the quantum chemical method The methods employed in the present study have their own limitations. In the article we have already discussed the performance of the PBE functional in reproducing the properties of bulk FeS 2. The same functional has been shown to perform well in modeling reaction occurring at the pyrite-water interface 8, 9, 10 and also in bulk water as well. 8, 9, 11 Nevertheless it is important to keep in mind that semilocal functionals yield somewhat underestimated barriers because they tend to overstabilize delocalized states (transition states usually feature more delocalized electron distribution than intermediates) due to the uncomplete correction of the self-interaction in the exchange-correlation energy. 12 An important issue is the inclusion of empirical dispersion corrections. 13 Such corrections are highly important when long range interactions play an important role. In our systems the reactive events are localized, therefore we do not expect large improvements from disperion corrections. In addition we have found that including dispersion correction has negative effect on the pyrite cell constant. In the light of the above discussion we have not used empirical dispersion corrections in the calculations. Additional errors may arise from the finite basis sets used to expand the orbitals and densities. To this end, test calculations have been performed on small molecules and on bulk pyrite. On the basis of these results we selected the MOLOPT basis set of double-zeta quality and augmented it with polarization functions for the S, O, N, and H atoms. This basis set yielded nice agreement between theory and either experiment or other calculations for structural properties of the surface and individual molecules and for bulk properties (cell parameter, bulk modulus, band-gap). In the simulations we have considered reaction paths featuring closed shell configurations. Certainly the extreme conditions and the presence of the iron may initiate open shell (radical) routes. However, their exploration requires far more computational efforts. In this regard it is necessary to note that the barriers found in the present simulations represent upper bounds for the NO x NH 3 reaction because possible more favorable radical routes will necessarily have lower barriers. 5.2 Issues of the free energy samplings The statistical error of the free energies obtained from metadynamics depends on the various parameters of the metadynamics simulations. For ideal cases Laio et al. derived suitable formulas which can be, however, extremely time-consuming to evaluate in real cases. 14, 15, 16 For example in a domain of the CVs with cubic shape the error can be estimated by the following equation: ( ) π2 k 2 δs 2 2S 2 ɛ 2 ( = S2 w βdτ g δs ) d 2π exp S k 0 π 2 k 2 (1) S8
9 where S is the size of the free energy well, δs is the hill width (assuming spherical hills), w is the hill height, τ g is the hill deposition frequency, β = 1/k B T, D is the diffusion coefficient of the CVs and d is the dimensionality of the CV space. k is a d dimensional vector with positive or null integers and k 2 is its square norm. The diffusion coefficient D can be calculated from the velocity-velocity autocorrelation function of the CVs. A more practical approach is to calculate the standard deviation of the free energy values obtained from repeated simulations with different initial conditions. In this study we have followed this strategy. For each barrier at least 3 independent estimates have been obtained. Table 4 summarizes the error estimates for the calculated barriers obtained in this way. Table 4: Free energy barriers and their standard deviations. Values are in kcal/mol. process activation free energy error perfect surface NO 3 NO NO 2 Fe-NO defective surface NO 3 NO NO 2 Fe-NO Fe-NO Fe N Fe-NO Fe-NHO Fe N Fe-NH Fe-NHO Fe-NH Fe-NHO Fe-NH 2 O Fe-NH Fe-NH Fe-NH 2 O Fe-NH Fe-NH 2 Fe-NH The finite size of the model system limits the accessible concentrations of the solute in water. In a real experiment the concentration of the reactants and products may be significantly different from those we could model in our simulation cell. The effect of this uncertainty on the barrier height can be estimated by F = F (c 2 ) F (c 1 ) = k B T ln c 1 c 2 (2) Even at a 100-fold concentration variation this translates to only 4.6 k B T (7.4 kcal/mol) stabilization, ie. the barrier heights may increase by this amount of energy. This increase should be taken into account only for the reactions of NO 3 and NO 2 surface. as the other intermediates prefer to stay at the Additional conceptual issues arise from the complexity of the real system that we wish to model. S9
10 Concentration range of the reactants, ph, fluctuations in temperature and pressure, various surface imperfections, different roles of possible radicals are not taken into account in the present models. However, the very large number of combination of these issues is presently not tractable with the available computational technology. Instead we focused on well-defined model reactions to address the NH 3 formation at a molecular scale. S10
11 6 Representative snapshots of transition states A B Figure S4: Representative configurations of the transition state ensembles for reactions occurring on the perfect FeS 2 surface. A: NO 3 NO 2. B: NO 2 NO. For color code see Fig. 1 of the article. A B C D E F G H I J Figure S5: Representative configurations of the transition state ensembles for reactions occurring on the defective FeS 2 surface. A: NO 3 NO 2. B: NO 2 NO. C: NO Fe N. D: Fe-NO Fe-NHO. E: Fe N Fe-NH. F: Fe-NHO Fe-NH. G: Fe-NHO Fe-NH 2 O. H: Fe-NH Fe-NH 2. I: Fe-NH 2 O Fe-NH 2. J: Fe-NH 2 Fe-NH 3. For color code see Fig. 1 of the article. S11
12 7 Free energy surfaces We present here the free energy surfaces (FES) for the most important routes which have not been shown in the main article. The CVs of the simulations are associated with the CVs of Table 3. Figure S6: FES of the NO 3 NO2 transformation occurring on a defective FeS2 surface. The label indicates the reaction stages according to the article. CV1 = CN1; CV2 = CN4. The energy values are in kcal/mol. Figure S7: FES of the NO transformation occurring on a defective FeS2 surface. The 2 NO label indicates the reaction stages according to the article. CV1 = CN1; CV2 = CN4. The energy values are in kcal/mol. S12
13 Figure S8: FES of the NO Fe N Fe-NH Fe-NH2 NH3 transformation occurring on a defective FeS2 surface. The labels indicate the reaction intermediates according to the article. CV1 = CV9; CV2 = CN6. The energy values are in kcal/mol. The double minima at the nitride intermadiate region correspond to the built-in of the dissociated O atom into a neighbour Fe and S atom in the uppermost pyrite layers. Figure S9: FES of the NO NHO NH2 O NH2 NH3 transformation occurring on a defective FeS2 surface. The labels indicate the reaction stages according to the article. CV1 = CV9; CV2 = CN6. The energy values are in kcal/mol. The double minima at the NH2 region correspond to the built-in of the dissociated O atom into a neighbour Fe and S atom in the uppermost pyrite layers. S13
14 References (1) Henkelman, G.; Uberuaga, B. P.; Jonsson, H. J. Chem. Phys., 2000, 113, (2) Sacchi, M.; Galbraith, M. C. E.; Jenkins, S. J. PCCP, 2012, 14, (3) Halgren, T. A.; Lipscomb, W. N.; Chem. Phys. Lett, 1977, 49, 225. (4) Marzari, N.; Vanderbilt, D. Phys. Rev. B 1997, 56, (5) Berry, J. F.; Bill, E.; Bothe, E.; George, S. D.; Mienert, B.; Neese, F.; Wieghardt, K. Science 2006, 312, (6) Hohenberger, J.; Ray, K.; Meyer, K. Nature Comm. 2012, 3, 720. (7) Sprik, M. Faraday Discuss. 1998, 110, 437. (8) Schreiner, E.; Nair, N. N.; Marx, D. J. Am. Chem. Soc., 2008, 130, (9) Schreiner, E.; Nair, N.N.; Wittekindt, C.; Marx, D. J. Am. Chem. Soc., 2011, 133, (10) Nair, N. N.; Schreiner, E.; Marx, D. J. Am. Chem. Soc., 2006, 128, Pollet, R.; Boehme, C.; Marx, D. Orig. Life Evol. Bios., 2006, 36, (11) Schreiner, E.; Nair, N. N.; Marx, D. J. Am. Chem. Soc., , (12) Parr, R. G.; Yang, W. Density functional Theory of Atoms and Molecules, Oxford University Press, (13) Grimme, S. WIREs Comput. Mol. Sci. 2011, 1, (14) Bussi, G.; Laio, A.; Parrinello, M. Phys. Rev. Lett., 2006, 96, ; (15) Iannuzzi, M.; Laio, A.; Parrinello, M. Phys. Rev. Lett., 2003, 90, (16) Laio, A.; Gervasio, F. L. Rep. Prog. Phys., 2008, 71, S14
Supporting Information for. Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite. Edge Water Interface
Supporting Information for Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite Edge Water Interface Zhen Jiang, Konstantin Klyukin, and Vitaly Alexandrov,, Department
More informationCurrent address: Department of Chemistry, Hong Kong Baptist University, Kowloon Tong, Hong Kong,
Hydrolysis of Cisplatin - A Metadynamics Study Supporting Information Justin Kai-Chi Lau a and Bernd Ensing* b Department of Chemistry and Applied Bioscience, ETH Zurich, USI Campus, Computational Science,
More informationAu-C Au-Au. g(r) r/a. Supplementary Figures
g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationSupporting Information for
Supporting Information for Pb-activated Amine-assisted Photocatalytic Hydrogen Evolution Reaction on Organic-Inorganic Perovskites Lu Wang *,,, Hai Xiao, Tao Cheng, Youyong Li *,, William A. Goddard III
More informationYuan Ping 1,2,3*, Robert J. Nielsen 1,2, William A. Goddard III 1,2*
Supporting Information for the Reaction Mechanism with Free Energy Barriers at Constant Potentials for the Oxygen Evolution Reaction at the IrO2 (110) Surface Yuan Ping 1,2,3*, Robert J. Nielsen 1,2, William
More informationSupplementary Information
Pt II as proton shuttle during C H bond activation in the Shilov process Pietro Vidossich,* a Gregori Ujaque, a Agustí Lledós a a Departament de Química, Universitat Autònoma de Barcelona, 08193 Cerdanyola
More informationDesign of Efficient Catalysts with Double Transition Metal. Atoms on C 2 N Layer
Supporting Information Design of Efficient Catalysts with Double Transition Metal Atoms on C 2 N Layer Xiyu Li, 1, Wenhui Zhong, 2, Peng Cui, 1 Jun Li, 1 Jun Jiang 1, * 1 Hefei National Laboratory for
More informationSupporting Information for. Dynamics Study"
Supporting Information for "CO 2 Adsorption and Reactivity on Rutile TiO 2 (110) in Water: An Ab Initio Molecular Dynamics Study" Konstantin Klyukin and Vitaly Alexandrov,, Department of Chemical and Biomolecular
More informationSelectivity in the initial C-H bond cleavage of n-butane on PdO(101)
Supporting Information for Selectivity in the initial C-H bond cleavage of n-butane on PdO(101) Can Hakanoglu (a), Feng Zhang (a), Abbin Antony (a), Aravind Asthagiri (b) and Jason F. Weaver (a) * (a)
More informationLi ion migration in Li 3 PO 4 electrolytes: Effects of O vacancies and N substitutions. Winston-Salem, North Carolina 27106, USA
75 Downloaded 22 Dec 28 to 52.7.52.46. Redistribution subject to ECS license or copyright; see http://www.ecsdl.org/terms_use.jsp ECS Transactions, 3 (26) 75-82 (28).49/.35379 The Electrochemical Society
More informationChemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase
Supplementary Information Chemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase Jia Chen, Ye-Fei Li, Patrick Sit, and Annabella Selloni Department of Chemistry,
More informationSupporting Information
Supporting Information Indirect Four-Electron Oxygen Reduction Reaction on Carbon Materials Catalysts in Acidic Solutions Guo-Liang Chai* 1, Mauro Boero 2, Zhufeng Hou 3, Kiyoyuki Terakura 3,4 and Wendan
More informationFinding minimum energy paths using Gaussian process regression
Finding minimum energy paths using Gaussian process regression Olli-Pekka Koistinen Helsinki Institute for Information Technology HIIT Department of Computer Science, Aalto University Department of Applied
More informationJoint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012
2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.
More informationSupporting information for
Supporting information for What is the role of pyridinium in pyridine-catalyzed CO 2 reduction on p-gap photocathodes? Martina Lessio a and Emily A. Carter* b Contents: 1) Cluster calculations: general
More informationTECHNIQUES TO LOCATE A TRANSITION STATE
32 III. TECHNIQUES TO LOCATE A TRANSITION STATE In addition to the location of minima, the calculation of transition states and reaction pathways is an interesting task for Quantum Chemistry (QC). The
More information22 Path Optimisation Methods
22 Path Optimisation Methods 204 22 Path Optimisation Methods Many interesting chemical and physical processes involve transitions from one state to another. Typical examples are migration paths for defects
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationSupporting Information for
Supporting Information for Role of Four-Fold Coordinated Titanium and Quantum Confinement in CO 2 Reduction at Titania Surface by Donghwa Lee and Yosuke Kanai Computational Details We employ Density Functional
More informationSupport Information. For. Theoretical study of water adsorption and dissociation on Ta 3 N 5 (100) surfaces
Support Information For Theoretical study of water adsorption and dissociation on Ta 3 N 5 (100) surfaces Submitted to Physical Chemistry Chemical Physics by Jiajia Wang a, Wenjun Luo a, Jianyong Feng
More informationAtomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationPrerequisites for reliable modeling with first-principles methods. P. Kratzer Fritz-Haber-Institut der MPG D Berlin-Dahlem, Germany
Prerequisites for reliable modeling with first-principles methods P. Kratzer Fritz-Haber-Institut der MPG D-14195 Berlin-Dahlem, Germany Prerequisites for modeling (I) Issues to consider when applying
More informationInteraction of O 2 with Gold Clusters: Molecular and Dissociative Adsorption
4066 J. Phys. Chem. A 2003, 107, 4066-4071 Interaction of O 2 with Gold Clusters: Molecular and Dissociative Adsorption Bokwon Yoon, Hannu Ha1kkinen,* and Uzi Landman School of Physics, Georgia Institute
More informationLinking electronic and molecular structure: Insight into aqueous chloride solvation. Supplementary Information
Linking electronic and molecular structure: Insight into aqueous chloride solvation Ling Ge, Leonardo Bernasconi, and Patricia Hunt Department of Chemistry, Imperial College London, London SW7 2AZ, United
More informationCrystallographic Dependence of CO Activation on Cobalt Catalysts: HCP versus FCC
Crystallographic Dependence of CO Activation on Cobalt Catalysts: HCP versus FCC Jin-Xun Liu, Hai-Yan Su, Da-Peng Sun, Bing-Yan Zhang, and Wei-Xue Li* State Key Laboratory of Catalysis, Dalian Institute
More informationSupplementary Information
Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 2015 Supplementary Information Insights into the Synergistic Role of Metal-Lattice
More informationThe electric field as a novel switch for uptake/release of hydrogen storage in nitrogen. doped graphene
The electric field as a novel switch for uptake/release of hydrogen storage in nitrogen doped graphene Z. M. Ao, 1,* A. D. Hernández-Nieves, 2,3 F. M. Peeters 3 and S. Li 1 1 School of Materials Science
More informationCatalysis Lectures W.H. Green 5.68J/10.652J Spring Handouts: Norskov et al., J. Catalysis Imbihl and Ertl, Chem. Rev. (partial) Homework
Catalysis Lectures W.H. Green 5.68J/10.652J Spring 2003 Handouts: Norskov et al., J. Catalysis Imbihl and Ertl, Chem. Rev. (partial) Homework Major points: 1) Why reactions have barriers, and how catalysts
More informationElementary Steps of the Catalytic NO x Reduction with NH 3 : Cluster Studies on Reactant Adsorption at Vanadium Oxide Substrate
Elementary Steps of the Catalytic NO x Reduction with NH 3 : Cluster Studies on Reactant Adsorption at Vanadium Oxide Substrate M. Gruber and K. Hermann Inorg. Chem. Dept., Fritz-Haber-Institut der Max-Planck-Gesellschaft,
More informationDehydrogenation of cyclohexaphenylene on Cu(111): possible reaction pathways from ab initio simulations
2nd CP2K Tutorial, ETHZ 7-11 February 2011 Materials Science & Technology Nanostructures and adsorption on metallic surfaces: Dehydrogenation of cyclohexaphenylene on Cu(111): possible reaction pathways
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationEnduring Understandings & Essential Knowledge for AP Chemistry
Enduring Understandings & Essential Knowledge for AP Chemistry Big Idea 1: The chemical elements are fundamental building materials of matter, and all matter can be understood in terms of arrangements
More informationNucleophile. Reaction Intermediate. Introduction to Reaction mechanisms. Definitions 2/25/2012
Introduction to Reaction mechanisms Definition A reaction mechanism is the step by step sequence of elementary reactions by which overall chemical change occurs. It is also a detailed description of the
More informationSupporting Information: Surface Polarons Reducing Overpotentials in. the Oxygen Evolution Reaction
Supporting Information: Surface Polarons Reducing Overpotentials in the Oxygen Evolution Reaction Patrick Gono Julia Wiktor Francesco Ambrosio and Alfredo Pasquarello Chaire de Simulation à l Echelle Atomique
More informationPhotoinduced Water Oxidation at the Aqueous. GaN Interface: Deprotonation Kinetics of. the First Proton-Coupled Electron-Transfer Step
Supporting Information Photoinduced Water Oxidation at the Aqueous Interface: Deprotonation Kinetics of the First Proton-Coupled Electron-Transfer Step Mehmed Z. Ertem,,,* eerav Kharche,,* Victor S. Batista,
More informationSupporting Information
Inhibition at Perimeter Sites of Au/TiO2 Oxidation Catalyst by Reactant Oxygen Isabel Xiaoye Green, Wenjie Tang, Monica McEntee, Mattew Neurock, and John T. Yates, Jr. Supporting Information Table of Contents:
More informationFORMAL CHARGE AND OXIDATION NUMBER
FORMAL CHARGE AND OXIDATION NUMBER Although the total number of valence electrons in a molecule is easily calculated, there is not aways a simple and unambiguous way of determining how many reside in a
More informationImproved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 22 8 DECEMBER 2000 Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points Graeme Henkelman a) and
More informationPerforming Metadynamics Simulations Using NAMD
Performing Metadynamics Simulations Using NAMD Author: Zhaleh Ghaemi Contents 1 Introduction 2 1.1 Goals of this tutorial.......................... 2 2 Standard metadynamics simulations 3 2.1 Free energy
More informationThe Study of Chemical Reactions. Mechanism: The complete, step by step description of exactly which bonds are broken, formed, and in which order.
The Study of Chemical Reactions Mechanism: The complete, step by step description of exactly which bonds are broken, formed, and in which order. Thermodynamics: The study of the energy changes that accompany
More informationSupporting Information
alladium monophosphine d(h 3 ): does it really exist in solution? ietro Vidossich*, Gregori Ujaque,* and Agustí Lledós* Supporting Information Contents: 1. Schemes S1 and S2 2. Computational details 3.
More informationSupplementary Information
Supplementary Information Supplementary Figure 1. X-ray diffraction patterns of (a) pure LDH, (b) AuCl 4 ion-exchanged LDH and (c) the Au/LDH hybrid catalyst. The refined cell parameters for pure, ion-exchanged,
More informationStructure and dynamics of the diarsenic complex in crystalline silicon
Structure and dynamics of the diarsenic complex in crystalline silicon Scott A. Harrison, Thomas F. Edgar, and Gyeong S. Hwang* Department of Chemical Engineering, University of Texas, Austin, Texas 78713,
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationCHEM J-5 June 2014
CHEM1101 2014-J-5 June 2014 The molecular orbital energy level diagrams for H 2, H 2 +, H 2 and O 2 are shown below. Fill in the valence electrons for each species in its ground state and label the types
More informationSupporting Information
Supporting Information The effect of strong acid functional groups on electrode rise potential in capacitive mixing by double layer expansion Marta C. Hatzell, a Muralikrishna Raju, b Valerie J. Watson,
More informationExam 1 Chemistry 142, Spring 2005 March 2, 2005
Exam Chemistry 4, Spring 005 March, 005 Part. Answer 7 of the following 8 multiple choice questions. If you answer more than 7 cross out the one you wish not to be graded, otherwise only the first 7 will
More informationPre-yield non-affine fluctuations and a hidden critical point in strained crystals
Supplementary Information for: Pre-yield non-affine fluctuations and a hidden critical point in strained crystals Tamoghna Das, a,b Saswati Ganguly, b Surajit Sengupta c and Madan Rao d a Collective Interactions
More informationMatter. Atoms and Elements
Matter Atoms and Elements Chemistry: The study of matter and its changes So, what is matter and what isn t? Matter isn t: Energy in the form of waves-- Light, Electricity, Heat, Magnetic Fields, Sound,
More informationKohn Sham density functional theory [1 3] is. Role of the Exchange Correlation Energy: Nature s Glue STEFAN KURTH, JOHN P. PERDEW.
Role of the Exchange Correlation Energy: Nature s Glue STEFAN KURTH, JOHN P. PERDEW Department of Physics and Quantum Theory Group, Tulane University, New Orleans, Louisiana 70118 Received 11 March 1999;
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/315/5819/1692/dc1 Supporting Online Material for Enhanced Bonding of Gold Nanoparticles on Oxidized TiO 2 (110) D. Matthey, J. G. Wang, S. Wendt, J. Matthiesen, R. Schaub,
More informationChapter Seven. Chemical Reactions: Energy, Rates, and Equilibrium
Chapter Seven Chemical Reactions: Energy, Rates, and Equilibrium Endothermic vs. Exothermic 2 Endothermic: A process or reaction that absorbs heat and has a positive ΔH. Exothermic: A process or reaction
More informationProperties of Individual Nanoparticles
TIGP Introduction technology (I) October 15, 2007 Properties of Individual Nanoparticles Clusters 1. Very small -- difficult to image individual nanoparticles. 2. New physical and/or chemical properties
More informationUniversity of Zurich. Proton transfer in imidazole-based molecular crystals. Zurich Open Repository and Archive. Iannuzzi, M.
University of Zurich Zurich Open Repository and Archive Winterthurerstr. 190 CH-8057 Zurich http://www.zora.uzh.ch Year: 2006 Proton transfer in imidazole-based molecular crystals Iannuzzi, M Iannuzzi,
More informationLow pressure CO 2 hydrogenation to methanol over gold nanoparticles activated on a CeO x /TiO 2 interface
Low pressure CO 2 hydrogenation to methanol over gold nanoparticles activated on a CeO x /TiO 2 interface 1 Xiaofang Yang, 1 Shyam Kattel, 1 Sanjaya D. Senanayake, 2 J. Anibal Boscoboinik, 3 Xiaowa Nie,
More informationDefects in TiO 2 Crystals
, March 13-15, 2013, Hong Kong Defects in TiO 2 Crystals Richard Rivera, Arvids Stashans 1 Abstract-TiO 2 crystals, anatase and rutile, have been studied using Density Functional Theory (DFT) and the Generalized
More information10/26/2010. An Example of a Polar Reaction: Addition of H 2 O to Ethylene. to Ethylene
6.5 An Example of a Polar Reaction: Addition of H 2 O to Ethylene Addition of water to ethylene Typical polar process Acid catalyzed addition reaction (Electophilic addition reaction) Polar Reaction All
More informationThe Low Temperature Conversion of Methane to Methanol on CeO x /Cu 2 O catalysts: Water Controlled Activation of the C H Bond
The Low Temperature Conversion of Methane to Methanol on CeO x /Cu 2 O catalysts: Water Controlled Activation of the C H Bond Zhijun Zuo, a Pedro J. Ramírez, b Sanjaya Senanayake, a Ping Liu c,* and José
More informationElectronic structure simulations of water solid interfaces
Electronic structure simulations of water solid interfaces Angelos Michaelides London Centre for Nanotechnology & Department of Chemistry, University College London www.chem.ucl.ac.uk/ice Main co-workers:
More information1 Molecular collisions
1 Molecular collisions The present exercise starts with the basics of molecular collisions as presented in Chapter 4 of the lecture notes. After that, particular attention is devoted to several specific
More informationSection 8.1 The Covalent Bond
Section 8.1 The Covalent Bond Apply the octet rule to atoms that form covalent bonds. Describe the formation of single, double, and triple covalent bonds. Contrast sigma and pi bonds. Relate the strength
More informationCHEM 2010 Symmetry, Electronic Structure and Bonding Winter Numbering of Chapters and Assigned Problems
CHEM 2010 Symmetry, Electronic Structure and Bonding Winter 2011 Numbering of Chapters and Assigned Problems The following table shows the correspondence between the chapter numbers in the full book (Physical
More informationAb-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University
Ab-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University Hands-on Tutorial Workshop on Ab-Initio Molecular Simulations, Fritz- Haber-Institut, Berlin, July 12-21,
More informationIonic Bonds. H He: ... Li Be B C :N :O :F: :Ne:
Ionic Bonds Valence electrons - the electrons in the highest occupied energy level - always electrons in the s and p orbitals - maximum of 8 valence electrons - elements in the same group have the same
More informationSTRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY S. Hara, T. Kumagai, S. Izumi and S. Sakai Department of mechanical engineering, University of
More informationSupporting Information
1 Supporting Information 2 3 4 Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds 5 6 7 Fabio Pietrucci 1, José C. Aponte 2,3,*, Richard Starr 2,4, Andrea Pérez-Villa
More informationSupporting Online Materials: Nature of proton. transport in a water-filled carbon nanotube and in. liquid water
Supporting Online Materials: Nature of proton transport in a water-filled carbon nanotube and in liquid water Ji Chen, Xin-Zheng Li,,, Qianfan Zhang, Angelos Michaelides, and Enge Wang, ICQM and School
More informationTheoretical studies on the bifunctionality of chiral thiourea-based organocatalysts: Competing routes to C C bond formation
Supporting Information Theoretical studies on the bifunctionality of chiral thiourea-based organocatalysts: Competing routes to C C bond formation Andrea Hamza, a Gábor Schubert, a Tibor Soós b and Imre
More informationSupplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water
Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Thomas D. Kühne 1, 2 and Rustam Z. Khaliullin 1, 1 Institute of Physical
More informationDescribe how the inter-conversion of solids, liquids and gases are achieved and recall names used for these inter-conversions
Understand the arrangements, movements and energy of the particle in each of the 3 states of matter : solid, liquid and gas Describe how the inter-conversion of solids, liquids and gases are achieved and
More informationChem 101 General Chemistry Practice Final Exam
Name h = 6.626 x 10-34 J s (Planck s Constant) c = 3.00 x 10 8 m/s (speed of light) R H = 1.097 x 10-7 m -1 (Rydberg Constant) Chem 101 General Chemistry Practice Final Exam Multiple Choice (5 points each)
More informationChapter 9 Ionic and Covalent Bonding
Chem 1045 Prof George W.J. Kenney, Jr General Chemistry by Ebbing and Gammon, 8th Edition Last Update: 06-April-2009 Chapter 9 Ionic and Covalent Bonding These Notes are to SUPPLIMENT the Text, They do
More informationName: Hr: 8 Basic Concepts of Chemical Bonding
8.1-8.2 8.3-8.5 8.5-8.7 8.8 Name: Hr: 8 Basic Concepts of Chemical Bonding 8.1 Chemical Bonds, Lewis Symbols, and the Octet Rule State the type of bond (ionic, covalent, or metallic) formed between any
More information4 Diatomic molecules
s manual for Burrows et.al. Chemistry 3 Third edition 4 Diatomic molecules Answers to worked examples WE 4.1 The Lewis model (on p. 174 in Chemistry 3 ) Use the Lewis model to describe the bonding in (a)
More informationImproved initial guess for minimum energy path calculations. Søren Smidstrup, 1, 2 Andreas Pedersen, 2, 3 Kurt Stokbro, 1 and Hannes Jónsson 2, 4
Improved initial guess for minimum energy path calculations Søren Smidstrup, 1, 2 Andreas Pedersen, 2, 3 Kurt Stokbro, 1 and Hannes Jónsson 2, 4 1) QuantumWise A/S, Lersø Parkallé 107, DK-2100 Copenhagen,
More informationCrossing the barriers - simulations of activated processes
Crossing the barriers - simulations of activated processes Mgr. Ján Hreha for 6 th Student Colloquium and School on Mathematical Physics Faculty of Mathematics, Physics and Informatics Comenius University
More informationCampbell Biology, 11e (Urry) Chapter 2 The Chemical Context of Life. 2.1 Multiple-Choice Questions
Campbell Biology, 11e (Urry) Chapter 2 The Chemical Context of Life 2.1 Multiple-Choice Questions 1) About 25 of the 92 natural elements are known to be essential to life. Which 4 of these 25 elements
More informationThiourea Derivatives as Brønsted Acid Organocatalysts
Supporting Information Thiourea Derivatives as Brønsted Acid Organocatalysts Ádám Madarász, Zsolt Dósa, Szilárd Varga, * Tibor Soós, Antal Csámpai, Imre Pápai * Institute of Organic Chemistry, Research
More informationUnderstanding Chemical Reactions through Computer Modeling. Tyler R. Josephson University of Delaware 4/21/16
Understanding Chemical Reactions through Computer Modeling Tyler R. Josephson University of Delaware 4/21/16 A little about me B.S. in Chem E from U of M, 2011 Currently, Ph.D. student at University of
More informationBasis Set for Molecular Orbital Theory
Basis Set for Molecular Orbital Theory! Different Types of Basis Functions! Different Types of Atom Center Basis Functions! Classifications of Gaussian Basis Sets! Pseudopotentials! Molecular Properties
More informationCHAPTER 2. Structure and Reactivity: Acids and Bases, Polar and Nonpolar Molecules
CHAPTER 2 Structure and Reactivity: Acids and Bases, Polar and Nonpolar Molecules 2-1 Kinetics and Thermodynamics of Simple Chemical Processes Chemical thermodynamics: Is concerned with the extent that
More informationWater structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials
Bull. Mater. Sci., Vol. 31, No. 3, June 2008, pp. 525 532. Indian Academy of Sciences. Water structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials
More informationStudy Guide Exam 1 BIO 301L Chinnery Spring 2013
Study Guide Exam 1 BIO 301L Chinnery Spring 2013 Lecture 1: Biology and Life What are the big picture messages from this lecture? How can you associate the contents of this lecture with those of the other
More informationXiang-Kui Gu,, Botao Qiao,,, Chuan-Qi Huang, Wu-Chen Ding, Keju Sun, Ensheng Zhan,, Tao Zhang, Jingyue Liu*,,, and Wei-Xue Li*,
Supported Single Pt 1 /Au 1 Atoms for Methanol Steam Reforming Xiang-Kui Gu,, Botao Qiao,,, Chuan-Qi Huang, Wu-Chen Ding, Keju Sun, Ensheng Zhan,, Tao Zhang, Jingyue Liu*,,, and Wei-Xue Li*, State Key
More informationDioxygen molecule adsorption and oxygen atom. oxygen atom diffusion on clean and defective Aluminum (111) surface using first principles calculations.
Dioxygen molecule adsorption and oxygen atom diffusion on clean and defective Aluminum (111) surface using first principles calculations Mathilde Guiltat, Marie Brut, Sébastien Vizzini, Anne Hémeryck To
More information(Ox) 6I - (aq) + BrO 3 - (aq) + 6H + (aq) 3I 2 (aq) + Br - (aq) + 3H 2 O(l)
Creating an Energy Profile For the Aqueous xidation of odide by omate in Acid The net ionic equation for the process we are investigating (x) is depicted below. (x) 6 - (aq) + - (aq) + 6 + (aq) 2 (aq)
More informationChemistry Final Review 2017
Chemistry Final Review 2017 Atomic/Molecular Structure and Periodic Trends 1. What is the atomic number trend on the periodic table? 2. On the following periodic table label metals, nonmetals, Alkali metals,
More informationChapter 6 Chemical Reactivity and Mechanisms
Chapter 6 Chemical Reactivity and Mechanisms 6.1 Enthalpy Enthalpy (ΔH or q) is the heat energy exchange between the reaction and its surroundings at constant pressure Breaking a bond requires the system
More informationAtoms, electrons and Solids
Atoms, electrons and Solids Shell model of an atom negative electron orbiting a positive nucleus QM tells that to minimize total energy the electrons fill up shells. Each orbit in a shell has a specific
More informationCHEM 150. Time: 90 Mins ATTEMPT ALL THE QUESTIONS
CHEM 150 Section 01, Q2 2016 Midterm 1 Student name... Student number... Time: 90 Mins ATTEMPT ALL THE QUESTIONS 1 Formulae and constants pv = nrt P 1 V 1 T 1 = P 2V 2 T 2 Ptotal = p1 + p2 + p3 +... U
More informationTitle of file for HTML: Supplementary Information Description: Supplementary Figures, Supplementary Tables and Supplementary References
Title of file for HTML: Supplementary Information Description: Supplementary Figures, Supplementary Tables and Supplementary References Title of file for HTML: Supplementary Movie 1 Description: This movie
More informationBreakdown of cation vacancies into anion vacancy-antisite complexes on III-V semiconductor surfaces
Breakdown of cation vacancies into anion vacancy-antisite complexes on III-V semiconductor surfaces A. Höglund and S. Mirbt Department of Physics, Uppsala University, Box 530, SE-75121 Uppsala, Sweden
More informationAb- ini&o quantum simula+ons of complex photocataly+c reac+ons
Ab- ini&o quantum simula+ons of complex photocataly+c reac+ons With: Dmitry Vinichenko, Grigory Kolesov, Georgios Tritsaris Department of Physics, Department of Chemistry and Chemical Biology, and School
More informationarxiv: v2 [physics.chem-ph] 12 Dec 2013
Car-Parrinello Simulation of the Reaction of arxiv:18.4969v [physics.chem-ph] 1 Dec 13 Aluminium with Oxygen Marius Schulte and Irmgard Frank, Ludwig-Maximilians-Universität München, Department Chemie,
More information9.1.1 CHEMICAL EQUATIONS AND STOICHIOMETRY
9.1.1 CHEMICAL EQUATIONS AND STOICHIOMETRY Work directly from Zumdahl (Chapter 3). Work through exercises as required, then summarise the essentials of the section when complete. A chemical equation is
More information