Pratical course Theoretical Chemistry CO on Palladium Dr. Carine MICHEL
|
|
- Sherilyn Gregory
- 5 years ago
- Views:
Transcription
1 Pratical course Theoretical Chemistry CO on Palladium adsorption energies vibrations diffusion path Dr. Carine MICHEL
2 Introduction This hands-on session aims at providing you a first experience in simulating a metallic catalyst (here Pd) and its interaction with a molecule (here CO). The calculations will be done at the Density Functional Theory (DFT) level using the VASP package ( In this code, the wave function is expanded on a plane wave basis set leading to 3D-periodic simulations. The visualization will be performed using VMD ( Editing and modifying a text file can be done using GEDIT. The operating system of the computers is linux-based. This text assumes that the participants are beginners and relies on a short manual, The Survivor Guide, to connect, get started and use VMD, GEDIT, etc. WARNING: This is a practical course. Parameters have been chosen in order to reduce as much as possible the computational cost, not to provide accurate results. 2
3 C. M ICHEL P RATICAL C OURSE Bulk Atomic positions The POSCAR file Pd crystalizes in the face centered cubic (fcc) lattice. The corresponding conventional cell contains 4 atoms and is represented in Figure 1. Figure 1: Conventional cell of the fcc lattice. This cube of side a contains 4 atoms: one at the corner, and three face-centered. The cell vectors can be written: A = ( a 0 0) B = (0 a 0) C = (0 0 a ) And using those three vectors as a basis, one can write the coordinates of the atoms: Pd1 = ( ) (1) Pd2 = ( ) (2) Pd3 = ( ) (3) Pd4 = ( ) (4) Those structural informations are collected in the POSCAR file. It is provided in the bulk directory. You can open it using GEDIT clicking on the file. You can also visualize the corresponding structure using VMD (see the Survivor Guide) clicking on poscar.vmd. 3
4 POSCAR Comments Pd bulk Title line 3.89 Cell parameter in Å(scaling factor applied to the cell vectors) A vector B vector C vector 4 Number of atoms in the cell Direct Type of coordinates (direct or cartesian) Pd Pd Pd Pd4 1.2 Atom kind The POTCAR file With a plane wave basis set, one has to use a pseudo-potential to replace the core electrons. The related information is in the POTCAR file. Open it and report here the following information: Atom name Atomic number Electronic configuration Minimal energy cutoff ENMIN Maximum energy cutoff ENMAX The DFT functional 4
5 1.3 Calculation parameters The INCAR file This file contains a list of keywords that controls the way the calculation is done (numerical precision, choice of the DFT functional, type of calculation, etc.) INCAR Pd bulk PREC = Normal ALGO = Fast ENCUT = 200 IBRION = -1 ISMEAR = -5 Comments Title line Global keyword that controls the numerical precision Type of algorithm to converge the electronic energy Energy cut off on the plane wave basis set Compute the energy of the system Smearing method 1.4 K-points grid The KPOINTS file The mesh of K-points is used to integrate over the 3D-Brillouin zone. It is defined by the KPOINTS file. KPOINTS Comments Automatic Mesh Title line 0 Automatic generation scheme Gamma Generate a Γ-centered mesh Monkhorst Pack Scheme used to generate the mesh Number of K-points in each direction Optional shift of the mesh 1.5 First calculation We are now ready to perform our first calculation. Click on vasp.j to run VASP. VASP generates many files by default. Some are deleted automatically at the end of the job. We want to know the energy of the system that can be found in several files. Open the OSZICAR file. Report here the energy: E0 = 5
6 1.6 Optimizing the cell parameter We want to find the cell parameter that minimizes the energy. Redo the previous calculation changing the cell parameter in the POSCAR file. You can share the calculations with other participants. Cell parameter a (Å) Energy (E0 = ) (ev) Table 1: Energy in function on the cell parameter for standard set up What is the cell parameter that minimizes the energy? What is the corresponding energy? How does this compare with the experimental cell parameter? 6
7 1.7 Influence of the energy cutoff Now, we want to see the importance of an appropriate energy cutoff on the plane wave basis set. We need to redo the same set of calculations for various values of the energy cut off ENCUT. The successive steps (changing the lattice parameter, click on vasp.j, report the energy found in the OSZICAR file) have been automatized in a script called bulk.j. Clicking on this file will generate a table in the file bulk.dat that collects the cell parameter and the corresponding energy. Modify the INCAR file to have the appropriate ENCUT value, run bulk.j, open bulk.dat and report the optimal cell parameter and the corresponding energy in Table 2. ENCUT (ev) Optimal cell parameter Energy (E0 =) (ev) ENMIN 200 ENMAX Table 2: Influence of the plane wave basis set (ENCUT) on the bulk calculation For the following, we will choose the ENCUT value to ensure an energy converged to 10 mev/atom: ENCUT =... 7
8 1.8 Influence of the k-point mesh In addition to the energy cutoff, an other important parameter is the quality of the K-points mesh. It is defined in the KPOINTS file. You should test this with a sufficient energy cutoff ENCUT (see the previous paragraph). What is the minimal K-point mesh to ensure an energy converged to 10 mev/atom? KPOINTS Optimal cell parameter Energy (E0 =) (ev) 7x7x7 9x9x9 11x11x11 13x13x13 Table 3: Influence of the K-points mesh on the bulk calculation at a given energy cutoff 1.9 The best set up To conclude this section, what is the best set up? and what is the cell parameter and the corresponding energy using it? 8
9 2 Surface Let s now get closer to a catalyst, cleaving the bulk along the (001) direction and the (111) and compare the two surfaces we obtain. The surfaces are modeled by a slab of finite thickness. Those two surfaces can be described using a small primitive cell (u, v) containing only atom in each layer. However, we want to study the adsorption of CO on each surface (see Paragraph 4 page 16). Then, we need larger cells to screen coverage lower than 1ML. For instance, a p(2x2) corresponds to a supercell (U, V) that is twice the primitive cell in the two directions (U = 2 u, V = 2 v). This supercell contains 4 Pd atoms per layer and the introduction of one CO in this cell leads to a coverage of 1/4ML. See Figure 2. (a) (b) Figure 2: Top view of the (a) Pd(001), p(2x2) cell (b) Pd(111), p(2x2) cell. Slabs are intrinsically 2D-system but VASP is a 3D-periodic code. So, our surface model consists in a slab of a given thickness, 2D periodic but repeated also periodically in the third direction (perpendicularly to the surface plane). The C vector has to be chosen large enough to avoid interactions between the periodic images of the slab in this direction. See Figure 3. 9
10 Figure 3: Side view of the Pd(111) slab. The three layers are represented in red/blue/white. The p(2x2) cell is repeated in U,V,z directions. 10
11 2.1 The (001) surface Energy calculation To start, we will compute the energy of a (001) surface as cut from the bulk. Here again, we need to use four input files. The POSCAR has been modified to provide the corresponding vectors and coordinates. We use the optimal Pd-Pd distance corresponding to the optimal cell parameter obtained in the previous section on bulk calculations. In addition, we will freeze the two bottom layers and optimize the positions of the two up layers using VASP. The POTCAR is unchanged. The INCAR has been adapted to slab calculation, changing the smearing method. Extra keywords are also added in this input file to perform a geometry optimization. More details can be found on the VASP website ( Last, in the direction perpendicular to the slab, we don t need k-points. Thus, the KPOINTS file is modified to reduce the k-points mesh at the Gamma point in this direction. It is also adapted to the size of the supercell p(2x2). See next page. 11
12 POSCAR Pd(001) 2.80 Pd-Pd distance in Å (scaling factor) U vector V vector C vector, 10 times the interlayer distance 4 Selective dynamics Only some coordinates will be optimized (T flags) Direct F F F Layer 1, x frozen, y frozen, z frozen F F F Layer 1, x frozen, y frozen, z frozen F F F Layer 1, x frozen, y frozen, z frozen F F F Layer 1, x frozen, y frozen, z frozen F F F Layer 2, x frozen, y frozen, z frozen F F F Layer 2, x frozen, y frozen, z frozen F F F Layer 2, x frozen, y frozen, z frozen F F F Layer 2, x frozen, y frozen, z frozen T T T Layer 3, x free, y free, z free T T T Layer 3, x free, y free, z free T T T Layer 3, x free, y free, z free T T T Layer 3, x free, y free, z free INCAR Geometry Optimisation Electronic minimization PREC = Normal ALGO = Fast ENCUT = 400 LREAL = Auto EDIFF = 1e-6 Ionic relaxation NSW = 50 IBRION = 2 POTIM = 0.5 EDIFFG = DOS related values ISMEAR = 2 Recommended for large systems Precision on the wavefunction Maximum number of optimization steps Geometry optimisation algorithm Step size Convergence criterium on the gradient Smearing method, different from the bulk one KPOINTS Automatic Mesh 0 Monkhorst Pack At the Gamma-point in the third direction 0 0 0
13 Run VASP clicking on vasp.j. Open the OSZICAR file. You can observe that you have performed several steps of optimization. Report the corresponding energies. E0 = E0 = E0 =... The first one corresponds to the initial structure, without any relaxation. The last one corresponds to the optimized structure. Other output files are of great interest. The XDATCAR contains all the geometries that VASP has tried to find the best possible one. The final geometry is reported in the CONTCAR file. You can visualize them with VMD and look at the Pd-Pd distance between the upper layers. You can also easily compute the interlayer distance using the CONTCAR file. How does it evolve during the geometry optimization? 2.2 The (111) surface In the slab111 directory, you will find the input files for VASP. Redo the same analysis than for the Pd(001) surface and compare the two facets. 13
14 3 CO molecule 3.1 Optimal distance and energy First, we need to know the equilibrium CO distance and the corresponding energy of CO in gas phase. We have to adapt the POSCAR, POTCAR and KPOINTS files. The INCAR is unchanged. Even if a CO molecule is aperiodic, VASP is a periodic code. Thus, we have to insert CO in a large box of 10Å 3 and do the computation at the Γ-point. The POTCAR concatenates the information concerning the C atom and then the O atom. Files are provided in the CO_gasphase directory. POSCAR CO molecule alone 1.0 scaling factor A vector B vector C vector carbon, 1 oxygen (the order is defined by the POTCAR) Cartesian Cartesian coordinates C atom at the center of the box O atom, CO along the x axis KPOINTS Automatic Mesh 0 Monkhorst Pack At the Gamma-point in the three directions Report here the energy and the CO distance in the optimal geometry. The energy will be used later to compute adsorption energies and the CO distance will be a good indicator of the back donation from the Pd d orbitals to the CO vacant π. 14
15 3.2 Vibration frequency We want to follow the evolution of the CO vibration frequency depending on the facet and the adsorption site. To start, we compute here the vibration frequency of CO vibration in vacuum in the CO_freq directory. We have to restart from the optimal geometry, provided in the POSCAR file. The INCAR is modified to do a frequency calculation. INCAR Freq PREC = Normal ALGO = Fast ENCUT = 400 IBRION = 5 POTIM = 0.02 NFREE = 2 NSW = 50 ISMEAR = 2 Frequency calculation Step used for the frequency calculations Number of points used for the frequency calculations Click on vasp.j. The frequency is provided in the freq.dat file: f = 15
16 4 4.1 CO adsorption For each fact (Pd(001) and Pd(111)), several sites have to be considered (see Figure 4) : top: the CO is sitting on a Pd atom bridge: the CO is bridging two Pd atoms hollow: the CO is coordinated to three (Pd(111)) or four (Pd(001)) Pd atoms. (a) (b) Figure 4: Top view of the (a) Pd(001): top (T), bridge (B) and hollow (h) sites (b) Pd(111): top (T), bridge (B) and hollow (fcc and hcp) sites We have to optimize the geometry of the CO adsorbed on each facet in each possible site. It is a bit long. Thus, the corresponding results are provided in the slab111_ CO_OPT and slab001_co_opt. The adsorption energy E ads is defined as the energy gained by the interaction between the surface and the molecule: E ads = E CO@slab E CO E slab (5) What is the most stable configuration on each facet? Can you comment the CO distance evolution? 16
17 E E CO E slab E ads (ev) CO distance (Å) Pd(001) top bridge hollow Pd(111) top bridge hcp fcc 17
18 4.2 CO vibration The frequencies corresponding to the optimal structures of CO adsorbed on Pd(111) and Pd(001) have been computed and collected in the slab111_co_freq and slab001_ CO_FREQ directories. They are listed in the freq.dat file. If the structure is a minimum of the potential energy surface (PES), all the frequencies are real. The highest wave number corresponds to the CO vibration. If the structure is a maximum of the PES, all the frequencies are imaginary. If the structure is a transition state, one and only one frequency is imaginary. Pd(001) top bridge hollow E ads (ev) Minimum? CO distance (Å ) σ(co) (cm 1 ) Pd(111) top bridge hcp fcc Could we distinguish the adsorption sites using infrared spectroscopies? 4.3 PDOS In this part, we will give a quick look at where are the electrons in CO@Pd(001), adsorbed in bridge position. The molecular orbitals of CO are given in the PRO- CAR file in the DOS directory. At each K-point, each band (=molecular orbital) is decomposed in term of atomic orbitals. C is the first atom, O the second one. Try 18
19 to complete the CO molecular orbitals diagram with the energies and the sketch of the MOs. E To know how those orbitals are modified upon adsorption, we consider now the density of states (DOS) of the CO@Pd(001) projected on those molecular orbitals. The projected DOS (PDOS) are given in the file PDOS-3 7.dat. The first column is the energy in ev. The second one is the PDOS on the MO3 of CO. The third one os the PDOS on the MO4 of CO, etc. The fermi level of the CO@Pd(001) is ev. You can plot the PDOS using Gnumeric clicking on the file PDOS-3 7.gnumeric. You could pay attention to: the dispersion of each band (the energy width ) the split of some bands the shift in energy the occupation of molecular orbitals that were previously unoccupied 19
20 4.4 CO diffusion on Pd(001) We will focus here on the CO diffusion on the Pd(001) surface as example of reactions paths. Three paths are possible (see Figure 5). You will find the results in the corresponding directories slab001_co_path. Path 3 Path 1 Path 2 Figure 5: The three possible paths for CO diffusion on Pd(001) NEB For each reaction path, we have optimized structures along the path (called images) using the Nudge Elastic Band method (NEB). The results are in the NEB subdirectory. The structures are reported in the XDATCAR and the energies in the energy profile.dat file. Plot the energy in function of geometrical parameters such as a Pd-C distance or the C-O distance. What is the structure closest to the transition state? 20
21 4.4.2 TS Starting from the structure the closest to the TS, we have optimize it. Looking at the results in TS_OPT sub-directory, report here the energy of this transition state. What is the energy barrier of this diffusion path? E0 = Activation energy = To check if the structure corresponds at a saddle point of order one, the frequencies have been computed (see TS_FREQ). Report the imaginary frequencies. Conclude. 21
22 4.4.3 Comparison Path Energy Profile TS Energy Activation Barrier Imaginary Frequencies 22
The VASP Scripter AddOn
The VASP Scripter AddOn Tutorial Version 11.8.1 The VASP Scripter AddOn: Tutorial Version 11.8.1 Copyright 2008 2011 QuantumWise A/S Atomistix ToolKit Copyright Notice All rights reserved. This publication
More informationGwangju Institute of Science and Technology (GIST) Super Computing & Collaboration Environment Technology Center (SCENT) Introduction to VASP 5.
Gwangju Institute of Science and Technology (GIST) Super Computing & Collaboration Environment Technology Center (SCENT) Introduction to VASP 5.2 Introduction to PLSI System Partnership & Leadership for
More informationA guide to the akmc script developed in Henkelman group
A guide to the akmc script developed in Henkelman group Lijun Xu and Graeme Henkelman Department of Chemistry and Biochemistry, University of Texas at Austin, Austin, Texas 78712-0165 (Dated: January 11,
More informationPhonons and Spectroscopy Tutorial
Phonons and Spectroscopy Tutorial The aims of this tutorial session is to introduce you to running larger-scale CASTEP jobs on supercomputer clusters, creation of input files for and setup of CASTEP jobs,
More informationThis file contains a quick introduction to the STM program: Hive (for VASP).
Dear user, This file contains a quick introduction to the STM program: Hive (for VASP). I ) Installation II) Condition of Usage III) Creating VASP input IV) Manual for Hive: How to generate pseudo-stm
More informationApril 26, Installation Requirements Getting started Help! The first run... 3
The cif2cell manual Torbjörn Björkman April 26, 2014 Contents 1 Installation 3 1.1 Requirements.................................. 3 2 Getting started 3 2.1 Help!.......................................
More informationOpenMX Installation and Basic Calculation Tutorial
OpenMX Installation and Basic Calculation Tutorial Mitchell Senger Oregon State University - Department of Physics Table of Contents 1 Preliminary Steps for Windows... 1 1.1 Installation of ssh software
More informationCastep Workshop 2007: Session 1
Castep Workshop 2007: Session 1 1 Login Before we start, you will have to log into the server that we are going to run Castep. Please follow the instructions given at the beginning of this practical session.
More informationUsing HyperChem for CH 463. Chun Park
Using HyperChem for CH 463 Chun Park Using HyperChem for CH 463 1 st Edition Written by Chun Park Oregon State University Table of Contents Introduction 1 Chapter 1 Getting Started 2 Starting HyperChem...
More informationAssignment #4 due 10/21/04
10.675 Assignment #4 due 10/21/04 In this problem set, you will use Car-Parrinello Molecular Dynamics (CPMD) to calculate the adsorption energy of oxygen atom on the Si(100) surface and compare it to the
More informationIntroduction to the CASTEP code. Stewart Clark University of Durham
CASTEP code Stewart Clark University of The CASTEP Developers Group (CDG) Stewart Clark () Keith Refson (STFC-RAL) Matt Probert (York) Chris Pickard (UCL) Phil Hasnip (York) Matt Segall (Industry/Cambridge)
More informationVirtual NanoLab. Tutorial. Version
Virtual NanoLab Tutorial Version 2014.2 Virtual NanoLab: Tutorial Version 2014.2 Copyright 2008 2014 QuantumWise A/S Atomistix ToolKit Copyright Notice All rights reserved. This publication may be freely
More informationOpenMX Tutorial. First trial On the manual Recommended trials OpenMX Viewer. Taisuke Ozaki (ISSP, Univ. of Tokyo)
OpenMX Tutorial First trial On the manual Recommended trials OpenMX Viewer Taisuke Ozaki (ISSP, Univ. of Tokyo) The Winter School on DFT: Theories and Practical Aspects, Dec. 19-23, CAS. First trial Login
More informationCHEM /12/01 INSTRUCTIONS FOR AB-INITIO CALCULATIONS Prepared by: Dr. Cynthia HARWOOD, Ahu AKIN, Prof. Tim KEIDERLING
CHEM 542 10/12/01 INSTRUCTIONS FOR AB-INITIO CALCULATIONS Prepared by: Dr. Cynthia HARWOOD, Ahu AKIN, Prof. Tim KEIDERLING SETTLING IN: - Type your login name : Chem542 and press Enter - Type your password
More informationNEAR-IR BROADBAND POLARIZER DESIGN BASED ON PHOTONIC CRYSTALS
U.P.B. Sci. Bull., Series A, Vol. 77, Iss. 3, 2015 ISSN 1223-7027 NEAR-IR BROADBAND POLARIZER DESIGN BASED ON PHOTONIC CRYSTALS Bogdan Stefaniţă CALIN 1, Liliana PREDA 2 We have successfully designed a
More informationAbinit Band Structure Maker MANUAL version 1.2
Abinit Band Structure Maker MANUAL version 1.2 written by Benjamin Tardif benjamin.tardif@umontreal.ca questions, bug report and suggestions are welcome! Purpose : The purpose of this program is to provide
More informationBasis Sets, Electronic Properties, and Visualization
Basis Sets, Electronic Properties, and Visualization Goals of the Exercise - Perform common quantum chemical calculations using the Gaussian package. - Practice preparation, monitoring, and analysis of
More informationTable of Contents. Table of Contents Coupling QuantumATK with Synopsys tools
Table of Contents Table of Contents Coupling QuantumATK with Synopsys tools Preparations Installing the addon New project Silicon crystal DFT model setup Running the calculation Visualizing the band structure
More informationTmoleX. A Graphical User Interface to the. Quantum Chemistry Program Package. Tutorial
TmoleX A Graphical User Interface to the TURBOMOLE Quantum Chemistry Program Package Tutorial COSMOlogic GmbH & Co. KG Burscheider Str. 515 D-51381 Leverkusen, Germany Phone +49-2171-363-668 Fax +49-2171-731-689
More informationU s e r' s G u i d e
U s e r' s G u i d e CALYPSO version 1.0 September 28, 2010 Written by Yanming Ma, Yanchao Wang, Jian Lv and Li zhu State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, China
More informationDynamic Photonic Crystal Superlattices
Dynamic Photonic Crystal Superlattices C. Neff, W. Park*, and C. J. Summers School of Materials Science & Engineering, Georgia Institute of Technology, Atlanta, GA *Department of Electrical & Computer
More informationAdvanced Graphics: NOMAD Summer. Interactive analysis and visualization of complex datasets
NOMAD Summer A hands-on course on tools for novel-materials discovery September 25-29, 2017, Berlin Advanced Graphics: Interactive analysis and visualization of complex datasets Michele Compostella Markus
More informationthink of the molecular atoms or ligands L
2 nd year Crystal and olecular Architecture course The Octahedral Point Group As an inorganic or materials chemist the octahedral point group is one of the most important, Figure o a vast number of T or
More informationHands-on Tutorial on Electronic Structure Computations: Transport with Wannier functions
Hands-on Tutorial on Electronic Structure Computations: Transport with Wannier functions tutor: Giovanni Borghi (giovanni.borghi@ep.ch) July 10, 2013 exercises by Giovanni Pizzi (2013) The les for this
More informationHands on Simulations of Interfacial Electron Transfer
Hands on Simulations of Interfacial Electron Transfer Snoeberger III, Robert C. Rudshteyn, Benjamin benjamin.rudshteyn@yale.edu Negre, Christian F.A. christianfannegre@gmail.com Ding, Wendu wendu.ding@yale.edu
More informationSingle-particle electron microscopy (cryo-electron microscopy) CS/CME/BioE/Biophys/BMI 279 Nov. 16 and 28, 2017 Ron Dror
Single-particle electron microscopy (cryo-electron microscopy) CS/CME/BioE/Biophys/BMI 279 Nov. 16 and 28, 2017 Ron Dror 1 Last month s Nobel Prize in Chemistry Awarded to Jacques Dubochet, Joachim Frank
More informationA Rough Guide to ELNES/EELS with Castep
A Rough Guide to ELNES/EELS with Castep Phil Hasnip March 27, 2008 1 Introduction Castep 4.3 is the first version of Castep to be capable of generating electron energy loss spectra (EELS); specifically,
More informationISOTROPY. Version 9.0. User s Manual. August 2007
ISOTROPY Version 9.0 User s Manual August 2007 Harold T. Stokes, Dorian M. Hatch, and Branton J. Campbell Department of Physics and Astronomy Brigham Young University ISOTROPY version 9.0 August 2007 isotropy
More informationpractical class: Step by step introduction to Yambo Excitons in LiF: comparison of theories & approximations H-chains: reasons for the ALDA failure
practical class: Step by step introduction to Yambo Excitons in LiF: comparison of theories & approximations H-chains: reasons for the ALDA failure Step by step introduction to Yambo * Technical introduction
More informationDFTB+ recipes. Release B. Aradi
DFTB+ recipes Release 11-09 B. Aradi November 20, 2014 Contents I Setting up DFTB+ 1 1 Compiling and installing 3 1.1 Prerequisites............................................. 3 1.2 Architecture dependent
More informationPorting CASTEP to GPGPUs. Adrian Jackson, Toni Collis, EPCC, University of Edinburgh Graeme Ackland University of Edinburgh
Porting CASTEP to GPGPUs Adrian Jackson, Toni Collis, EPCC, University of Edinburgh Graeme Ackland University of Edinburgh CASTEP Density Functional Theory Plane-wave basis set with pseudo potentials Heavy
More informationSIMULATION OF ARTIFICIAL SYSTEMS BEHAVIOR IN PARAMETRIC EIGHT-DIMENSIONAL SPACE
78 Proceedings of the 4 th International Conference on Informatics and Information Technology SIMULATION OF ARTIFICIAL SYSTEMS BEHAVIOR IN PARAMETRIC EIGHT-DIMENSIONAL SPACE D. Ulbikiene, J. Ulbikas, K.
More informationParallel Implementation of Nudged Elastic Band Method
Parallel Implementation of Nudged Elastic Band Method.33 Final Project Anubhav Sinha December, Background Methods in computational chemistry give potential energy surfaces (PES). Purpose: find the transition
More informationModeling Flow Through Porous Media
Tutorial 7. Modeling Flow Through Porous Media Introduction Many industrial applications involve the modeling of flow through porous media, such as filters, catalyst beds, and packing. This tutorial illustrates
More informationEngineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering
Engineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering Here SolidWorks stress simulation tutorials will be re-visited to show how they
More informationAPPENDIX 2. Density Functional Theory Calculation of Vibrational Frequencies with Gaussian 98W 1
APPENDIX 2 Density Functional Theory Calculation of Vibrational Frequencies with Gaussian 98W 1 This appendix describes the use of Gaussian 98W software for calculating vibrational frequencies for the
More informationOptimization of metallic biperiodic photonic crystals. Application to compact directive antennas
Optimization of metallic biperiodic photonic crystals Application to compact directive antennas Nicolas Guérin Computational Optic Groups (COG) IFH, ETHZ, http://alphard.ethz.ch Keywords: 3D modeling,
More informationBoth equations are solved using a finite differences (iterative relaxation) method, which takes some time to converge.
WEIGHTFIELD 2D Silicon Strip Detector Simulation V0.04 14 November 2011 markus.friedl@oeaw.ac.at Abstract WEIGHTFIELD is a program that allows simulating a silicon strip detector in two dimensions (crosssection).
More informationWhat is!?! Fortran95: math&sci. libxc etsf-io.
What is!?! Theory&approximations Many-Body Perturbation Theory Time-dependent density functional theory Bethe-Salpeter equation, GWA AHC electron-phonon coupling... Languages www.yambo-code.org Interfaces
More informationQuick Instruction. For EXAFS Data Analysis
Quick Instruction For EXAFS Data Analysis Using WinXAS, FEFF 7 and ATOMS Programs Farideh Jalilehvand Department of Chemistry University of Calgary e-mail: faridehj@ucalgary.ca Last updated: December 2002
More informationTEAMS National Competition High School Version Photometry Solution Manual 25 Questions
TEAMS National Competition High School Version Photometry Solution Manual 25 Questions Page 1 of 15 Photometry Questions 1. When an upright object is placed between the focal point of a lens and a converging
More informationAMS Manual. ADF Modeling Suite
AMS Manual ADF Modeling Suite 2018 www.scm.com Sep 18, 2018 CONTENTS 1 General 1 1.1 Overview................................................. 1 1.2 Motivation and progress.........................................
More informationOpen new file: File New Create MolGroup. Double Click at Ring Fragment Select Benzene Click at new window
Chem 4320 Working in GaussView: N N Open new file: File New Create MolGroup Double Click at Ring Fragment Select Benzene Click at new window Double Click Element Fragment on the Builder Click on nitrogen
More informationThe Cube and its Symmetries
The Cube and its Symmetries This workshop is designed to show Y8/9 children how symmetry operations interrelate. We take the cube as our example. Administration In your area is a visualiser, an OHP and
More informationFlair Geometry Editor Part I. Beginners FLUKA Course
Flair Geometry Editor Part I Beginners FLUKA Course Starting the Geometry Editor Click on icon or from Menu View Geometry Editor or with [F4] shortcut Either start flair with option -g 2 Geometry editor
More informationViewing Gaussian Structures with GaussView (GaussView Version: 3.09)
Viewing Gaussian Structures with GaussView (GaussView Version: 3.09) In GaussView, there are several ways to visualize structures from a Gaussian output file. The log file, checkpoint file, and formatted
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More information[ ] u 1. ME309 Homework #2 $ % & u = 1 s 2 " # u 2. s,u. A,E constant along length. 4Etc
ME09 Homework # OBJECTIVES: Introduction to convergence issues and modeling approaches Postprocessing procedures Element selection Experience with shape functions and stiffness terms NOTES: Problems 1and
More informationTEAMS National Competition Middle School Version Photometry Solution Manual 25 Questions
TEAMS National Competition Middle School Version Photometry Solution Manual 25 Questions Page 1 of 14 Photometry Questions 1. When an upright object is placed between the focal point of a lens and a converging
More informationInstallation Guide. Beginners guide to structural analysis
Installation Guide To install Abaqus, students at the School of Civil Engineering, Sohngaardsholmsvej 57, should log on to \\studserver, whereas the staff at the Department of Civil Engineering should
More informationRelease Notes: Thermo-Calc Software Package Version 2018a
Release Notes: Thermo-Calc Software Package Version 2018a 2018 Foundation of Computational Thermodynamics: Solna, Sweden Contents Release notes Thermo-Calc 2018a 2 Highlights 2 Software Development Kits
More informationNX Fixed Plane Additive Manufacturing Help
NX 11.0.2 Fixed Plane Additive Manufacturing Help Version #1 1 NX 11.0.2 Fixed Plane Additive Manufacturing Help June 2, 2017 Version #1 NX 11.0.2 Fixed Plane Additive Manufacturing Help Version #1 2 Contents
More informationSimple REFMAC tutorials
Simple REFMAC tutorials Prerequisites: To use this tutorial you need to have ccp4. For jelly body, automatic and local ncs restraints, occupancy refinement you need to have the latest version of ccp4-6.2
More informationSelective Space Structures Manual
Selective Space Structures Manual February 2017 CONTENTS 1 Contents 1 Overview and Concept 4 1.1 General Concept........................... 4 1.2 Modules................................ 6 2 The 3S Generator
More informationSIMULATION OF AN IMPLANTED PIFA FOR A CARDIAC PACEMAKER WITH EFIELD FDTD AND HYBRID FDTD-FEM
1 SIMULATION OF AN IMPLANTED PIFA FOR A CARDIAC PACEMAKER WITH EFIELD FDTD AND HYBRID FDTD- Introduction Medical Implanted Communication Service (MICS) has received a lot of attention recently. The MICS
More informationSupporting Information. Investigating the Nature of Palladium Chain-Walking in the Enantioselective Redox-Relay Heck Reaction of Alkenyl Alcohols
Supporting Information Investigating the Nature of Palladium Chain-Walking in the Enantioselective Redox-Relay eck Reaction of lkenyl lcohols Margaret J. ilton, 1 Li-Ping Xu, 2 Per-Ola Norrby, 3 Yun-ong
More informationAdvanced Visualization for Chemistry
Advanced Visualization for Chemistry Part 9 Tools in depth Mario Valle March 7 8, 2006 Plan 1. Introduce two tools 1. STM3 the tool used at CSCS (more developer oriented) 2. VMD widely used free tool (more
More informationApplied Reservoir Simulation. Course. Workshop Problems
Applied Reservoir Simulation Course Workshop Problems Reservoir Simulation Application Training Course and (Eclipse) Workshop GeoQuest Training and Development, And NExT Denver and Houston March 06 Applied
More informationOPSE FINAL EXAM Fall CLOSED BOOK. Two pages (front/back of both pages) of equations are allowed.
CLOSED BOOK. Two pages (front/back of both pages) of equations are allowed. YOU MUST SHOW YOUR WORK. ANSWERS THAT ARE NOT JUSTIFIED WILL BE GIVEN ZERO CREDIT. ALL NUMERICAL ANSERS MUST HAVE UNITS INDICATED.
More informationAssignment 4: Molecular Dynamics (50 points)
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard November 2, 2015 Assignment 4: Molecular Dynamics (50 points) In this assignment, you will use the Tinker molecular modeling software package to carry out
More informationStep by step introduction to Yambo
Step by step introduction to Yambo * Technical introduction to the code: # Flow of a calculation # Yambo command line options # Yambo I/O file START ~>a2y -F o.kss 1. Generate the core databases ~>yambo
More informationdiffraction patterns obtained with convergent electron beams yield more information than patterns obtained with parallel electron beams:
CBED-Patterns Principle of CBED diffraction patterns obtained with convergent electron beams yield more information than patterns obtained with parallel electron beams: specimen thickness more precise
More informationWorkshop 15. Single Pass Rolling of a Thick Plate
Introduction Workshop 15 Single Pass Rolling of a Thick Plate Rolling is a basic manufacturing technique used to transform preformed shapes into a form suitable for further processing. The rolling process
More informationInterference of Light
Interference of Light Young s Double-Slit Experiment If light is a wave, interference effects will be seen, where one part of wavefront can interact with another part. One way to study this is to do a
More informationManual for the mini-workshop
2014/01/07 00:40 1/11 Manual for the mini-workshop Manual for the mini-workshop Information SSH/STM4 softwares are installed in D:\ drive, you may find shortcut in D:\. To use STM4: double click on go.bat/stm4.bat
More informationHands-on introduction of OpenMX - Lecture -
1/31 Hands-on introduction of OpenMX - Lecture - Ozaki group Mitsuaki Kawamura Schedule: 14:30~15:30 Lecture 15:50~17:20 Practice 2018/7/2 ISS2018 1 st part 2/31 Outline Introduction : Basic OpenMX What
More informationnewfasant Periodical Structures User Guide
newfasant Periodical Structures User Guide Software Version: 6.2.10 Date: February 23, 2018 Index 1. FILE MENU 2. EDIT MENU 3. VIEW MENU 4. GEOMETRY MENU 5. MATERIALS MENU 6. CELL MENU 6.1. DEFINE CELL
More informationSupersonic Flow Over a Wedge
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 2 Supersonic Flow Over a Wedge Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification A uniform supersonic stream encounters a wedge
More informationWETTING PROPERTIES OF STRUCTURED INTERFACES COMPOSED OF SURFACE-ATTACHED SPHERICAL NANOPARTICLES
November 20, 2018 WETTING PROPERTIES OF STRUCTURED INTERFACES COMPOSED OF SURFACE-ATTACHED SPHERICAL NANOPARTICLES Bishal Bhattarai and Nikolai V. Priezjev Department of Mechanical and Materials Engineering
More information(Refer Slide Time: 00:10)
Fundamentals of optical and scanning electron microscopy Dr. S. Sankaran Department of Metallurgical and Materials Engineering Indian Institute of Technology, Madras Module 02 Unit-4 Phase contrast, Polarized
More informationCrystal Structure. A(r) = A(r + T), (1)
Crystal Structure In general, by solid we mean an equilibrium state with broken translational symmetry. That is a state for which there exist observables say, densities of particles with spatially dependent
More informationTiling Three-Dimensional Space with Simplices. Shankar Krishnan AT&T Labs - Research
Tiling Three-Dimensional Space with Simplices Shankar Krishnan AT&T Labs - Research What is a Tiling? Partition of an infinite space into pieces having a finite number of distinct shapes usually Euclidean
More informationspecular diffuse reflection.
Lesson 8 Light and Optics The Nature of Light Properties of Light: Reflection Refraction Interference Diffraction Polarization Dispersion and Prisms Total Internal Reflection Huygens s Principle The Nature
More informationc Fluent Inc. May 16,
Tutorial 1. Office Ventilation Introduction: This tutorial demonstrates how to model an office shared by two people working at computers, using Airpak. In this tutorial, you will learn how to: Open a new
More informationGetting Started with OpenMX Ver. 3.5
Getting Started with OpenMX Ver. 3.5 1. Getting the package and its installation 2. Running of test jobs 3. Input files 4. Output files 5. Basis sets and pseudopotentials 6. Functionals 7. Cutoff energy
More informationSupporting Information. High-Throughput, Algorithmic Determination of Nanoparticle Structure From Electron Microscopy Images
Supporting Information High-Throughput, Algorithmic Determination of Nanoparticle Structure From Electron Microscopy Images Christine R. Laramy, 1, Keith A. Brown, 2, Matthew N. O Brien, 2 and Chad. A.
More informationOverview of OpenMX: 3. Ongoing and planned works. 2. Implementation. 1. Introduction. 4. Summary. Implementation and practical aspects
Overview of OpenMX: Implementation and practical aspects 1. Introduction 2. Implementation Total energy, forces, numerical integrations, basis functions, getting SCF, geometry optimization, parallelization,
More informationInstallation. Chapter. Installing the Hardware Protection Key. Activating the Hardware. Changing Hardware or Versions of the MacOS
Chapter I Installation The application CrystalKit and its associated files are installed by double clicking on the installer package. After authorizing the installer with the administrator password, the
More informationTutorial 2. Modeling Periodic Flow and Heat Transfer
Tutorial 2. Modeling Periodic Flow and Heat Transfer Introduction: Many industrial applications, such as steam generation in a boiler or air cooling in the coil of an air conditioner, can be modeled as
More informationRadical Pd(III) C C Reductive Eliminations in Palladium Catalyzed Sequences
Electronic Supplementary Information for: Radical Pd(III) C C Reductive Eliminations in Palladium Catalyzed Sequences Giovanni Maestri, * Max Malacria, and Etienne Derat * Institut de Chimie des Substances
More informationPrincipal Roll Structure Design Using Non-Linear Implicit Optimisation in Radioss
Principal Roll Structure Design Using Non-Linear Implicit Optimisation in Radioss David Mylett, Dr. Simon Gardner Force India Formula One Team Ltd. Dadford Road, Silverstone, Northamptonshire, NN12 8TJ,
More informationPlane wave in free space Exercise no. 1
Plane wave in free space Exercise no. 1 The exercise is focused on numerical modeling of plane wave propagation in ANSYS HFSS. Following aims should be met: 1. A numerical model of a plane wave propagating
More informationDose Calculations: Where and How to Calculate Dose. Allen Holder Trinity University.
Dose Calculations: Where and How to Calculate Dose Trinity University www.trinity.edu/aholder R. Acosta, W. Brick, A. Hanna, D. Lara, G. McQuilen, D. Nevin, P. Uhlig and B. Slater Dose Calculations - Why
More informationSimulating the RF Shield for the VELO Upgrade
LHCb-PUB-- March 7, Simulating the RF Shield for the VELO Upgrade T. Head, T. Ketel, D. Vieira. Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil European Organization for Nuclear Research
More informationCrystal Fingerprinting. for USPEX output analysis
Crystal Fingerprinting and STM4 for USPEX output analysis Mario Valle USPEX 2011 Xi an workshop 04/08/2011 The CSCS Data Analysis and Visualization Group mvalle@cscs.ch 1 CSCS computer room CSCS usage
More informationPhysical Optics. You can observe a lot just by watching. Yogi Berra ( )
Physical Optics You can observe a lot just by watching. Yogi Berra (1925-2015) OBJECTIVES To observe some interference and diffraction phenomena with visible light. THEORY In a previous experiment you
More informationUsing ANSYS and CFX to Model Aluminum Reduction Cell since1984 and Beyond. Dr. Marc Dupuis
Using ANSYS and CFX to Model Aluminum Reduction Cell since1984 and Beyond Dr. Marc Dupuis 1980-84, 2D potroom ventilation model Physical model 1980-84, 2D potroom ventilation model Experimental results
More informationKey Features for OptiFDTD 14
14.0 New Features Created to address the needs of research scientists, photonic engineers, professors and students; OptiFDTD satisfies the demand of users who are searching for a powerful yet easy to use
More informationGeometry Optimization Made Simple with Translation and Rotation Coordinates
Geometry Optimization Made Simple with Translation and Rotation Coordinates Lee-Ping Wang 1 and Chenchen Song 2, 3 1) Department of Chemistry, University of California; 1 Shields Ave; Davis, CA 95616.
More informationMiroslav Kohout. STR-module for Avizo. User s Guide. Springer
Miroslav Kohout STR-module for Avizo User s Guide Springer Contents 1 Visualization with the STR module for Avizo...................... 1 1.1 Avizo module installation....................................
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature10934 Supplementary Methods Mathematical implementation of the EST method. The EST method begins with padding each projection with zeros (that is, embedding
More informationFLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016
FLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016 Note: These instructions are based on an older version of FLUENT, and some of the instructions
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 3 Common Elements of a Molecular Simulation David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Boundary Conditions Impractical
More informationChapter 36. Diffraction. Copyright 2014 John Wiley & Sons, Inc. All rights reserved.
Chapter 36 Diffraction Copyright 36-1 Single-Slit Diffraction Learning Objectives 36.01 Describe the diffraction of light waves by a narrow opening and an edge, and also describe the resulting interference
More informationAppendix Building Alkyl Semibenzoquinone Structures
WARNING NOTICE: The experiments described in these materials are potentially hazardous and require a high level of safety training, special facilities and equipment, and supervision by appropriate individuals.
More informationRevised Sheet Metal Simulation, J.E. Akin, Rice University
Revised Sheet Metal Simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering analysis.
More informationEnergy Minimization -Non-Derivative Methods -First Derivative Methods. Background Image Courtesy: 3dciencia.com visual life sciences
Energy Minimization -Non-Derivative Methods -First Derivative Methods Background Image Courtesy: 3dciencia.com visual life sciences Introduction Contents Criteria to start minimization Energy Minimization
More informationGoals of parallel computing
Goals of parallel computing Typical goals of (non-trivial) parallel computing in electronic-structure calculations: To speed up calculations that would take too much time on a single processor. A good
More informationGraphical User Interface
The NOMAD (Novel Materials Discovery) Laboratory a European Centre of Excellence Graphical User Interface Deliverable No: 2.4 Lead Beneficiary: Humboldt-Universitaet zu Berlin (HUB) Contributing Beneficiaries:
More informationEFFECT OF BAND STRUCTURE DISCRETIZATION ON THE PERFORMANCE OF FULL-BAND MONTE CARLO SIMULATION
EFFECT OF BAND STRUCTURE DISCRETIZATION ON THE PERFORMANCE OF FULL-BAND MONTE CARLO SIMULATION G. Karlowatz, W. Wessner, and H. Kosina, Institute for Microelectronics, TU Vienna, Austria Corresponding
More information