ReaxFF for Vanadium and Bismuth Oxides

1. ReaxFF for Vanadium and Bismuth Oxides Kim Chenoweth Force Field Sub-Group Meeting January 20, 2004

2. Overview • Significance of a Bi/V force field • ReaxFF: general principles • Force field optimization for V • Force field optimization for Bi • Future work

3. Cat Designing a Better Catalyst - I • 85% of industrial organic chemicals are currently produced by catalytic processes • 25% are produced by heterogeneous oxidation catalysis such as ammoxidation CH2=CHCH3 + NH3 + 3/2 O2 CH2=CHCN + 3 H2O • Bi-molybdates are currently used as the catalyst • Use of alkanes as a cheaper feedstock requires design of a selective catalyst • Promising catalysts are complex oxides containing Mo, V, Te, X, and O where X is at least one other element • Bismuth is one of the 19 elements listed in the Mitsubishi patent

4. Designing a Better Catalyst - II • Low-MW alkenes (i.e. ethene and propene) can be formed via non-oxidative dehydrogenation (ODH) of the corresponding alkane • Supported vanadia is the most active and selective simple metal oxide for alkane ODH1 • Due to its reducible nature, it leads to rapid redox cycles necessary for catalytic turnover • Local structure strongly influences ODH reaction rates and selectivity • Force field would allow for the study of large and complex systems with many atoms • Generate interesting structures for further study using QC methods • Optimize ratio of the various metals in the catalyst • Elucidate the purpose of the different metals 1Argyle et al, J. Catal. 2002, 208, 139

5. Atoms Molecular conformations years Design Electrons Bond formation FEA MESO Time Grids MD ReaxFF Grains QC Empirical force fields 10-15 ab initio, DFT, HF Ångstrom Kilometers Distance ReaxFFBridging the gap between QC and EFF • Empirical methods: • Study large system • Rigid connectivity • QC Methods: • Allow reactions • Expensive • ReaxFF: • Simulate bond formation in larger molecular systems

6. 3-body 4-body 2-body multi-body ReaxFF: Energy of the System • Similar to empirical non-reactive force fields • Divides the system energy into various partial energy contributions

7. Important Features in ReaxFF • A bond length/bond order relationship is used to obtain smooth transition from non-bonded to single, double, and triple bonded systems. • Bond orders are updated every iteration • Non-bonded interactions (van der Waals, coulomb) • Calculated between every atom pair • Excessive close-range non-bonded interactions are avoided by shielding • All connectivity-dependent interactions (i.e. valence and torsion angles) are made bond-order dependent • Ensures that their energy contributions disappear upon bond dissociation • ReaxFF uses a geometry-dependent charge calculation scheme that accounts for polarization effects

8. ReaxFF as a Transferable Potential General Rules: • No discontinuities in energy or forces even during reactions • No pre-defined reactive sites or reaction pathways • Should be able to automatically handle coordination changes associated with reactions • One force field atom type per element • Should be able to determine equilibrium bond lengths, valence angles, etc from chemical environment

9. Strategy for Parameterization of ReaxFF • Identify important interactions to be optimized for relevant systems • Build QC-training set for bond dissociation and angle bending cases for small clusters • Build QC-training set for condensed phases to obtain equation of state • Force field optimization using • Metal training set • Metal oxide clusters and condensed phases • Applications

10. Cluster Bonds -Normal, under-, and over-coordinated systems Angles O-V=O, V-O-V, O=V=O Condensed Phase Metal BCC, A15, FCC, SC, Diamond Metal Oxide VO (II) FCC V2O3 (III) Corundum VO2 (IV) Distorted rutile V2O5 (V) Layered octahedral Vanadium Training Set 1st row transition metal (4s23d3) • Successive bond dissociation of oxygen in V4O10

11. Bulk Metal - Vanadium QC ReaxFF • ReaxFF reproduces EOS and properly predicts instability of low-coordination phases (SC, Diamond)

12. Bond Dissociationin VO2OH

13. V=O Bond Dissociation in V4O10

14. Angle Distortion in V2O5 V-O-V Angle O-V=O Angle

15. Angle Distortion in VO2 O=V=O Angle

16. Angle Distortion in V2O6 V-O-O Angle

17. 4 3 1 2 2 1 3 3 6 4 2 5 1 7 Charge Analysis for VxOy Clusters in Training Set

18. Charge Analysis for VxOY Clusters in Literature(QC data taken from Calatayud et al, J. Phys. Chem. A 2001, 105, 9760.)

19. Bismuth Training Set Common oxidation states: 3, 5 • Cluster • Bonds • -Normal, under-, and over-coordinated systems • Angles • Bi-Bi=O, O=Bi-O • Condensed Phase • Metal • HCP, SC, BCC, A15, FCC, Diamond • Metal Oxide • BiO (II) • Trigonal • a-Bi2O3 (III) • Monoclinic • b-Bi2O3 (III) • Distorted cubic • Bi2O4 (BiIIIBiVO4) • Monoclinic • BiO2 (IV) • Cubic

20. Bulk Metal - Bismuth QC ReaxFF

21. Relative Stabilities of V and Bi Bulk Phases Bismuth Vanadium

22. 1900 K Application: Melting Point of Vanadium 55 molecules 900 K 1700 K 2500 K 1700 K 900 K • Melting point of Vanadium = 2163 K • Melting point obtained from simulation ~ 1900 K

23. Application: Melting Point of Vanadium 147 molecules 900 K 1700 K 2500 K 1700 K 900 K 2000 K • Melting point of Vanadium = 2163 K • Melting point obtained from simulation ~ 2000 K

24. Future Work • Bismuth oxide force field training set: • Optimization of Bi oxide force field • Add bond dissociation and bond angles for clusters • Add bismuth oxide condensed phases • Vanadium oxide force field training set: • Further optimization of vanadium oxide force field • Add successive V=O bond dissociation for V4O10 • Add vanadium oxide condensed phases Add to training set and continue optimizing force field Add to training set and continue optimizing force field