introduction introduction
play

Introduction Introduction Why to use a simulation Some examples - PowerPoint PPT Presentation

Introduction Introduction Why to use a simulation Some examples of questions we can address Molecular Simulations MD Molecular dynamics : solve equations of motion r 1 Monte Carlo : r 2 importance sampling r n Calculate


  1. Introduction

  2. Introduction Why to use a simulation • Some examples of questions we • can address

  3. Molecular Simulations MD • Molecular dynamics : solve equations of motion r 1 • Monte Carlo : r 2 importance sampling r n • Calculate thermodynamic MC and transport properties for a given intermolecular r 1 potential r 2 r n 3

  4. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system 4

  5. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system We assume the interactions between the particles are known! 4

  6. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system 4

  7. Uses of Molecular Simulations Exact= in the limit of infinitely long simulations the error bars can be made infinitely small The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system 4

  8. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system 4

  9. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system Pressure Heat capacity Heat of adsorption Structure …. 4

  10. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system 4

  11. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system Diffusion coefficient Viscosity … 4

  12. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system 4

  13. Uses of Molecular Simulations If one could envision an experimental system of these N particles that interact with the potential. The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system 4

  14. Uses of Molecular Simulations The idea for a given intermolecular potential “ exactly” compute the thermodynamic and transport properties of the system 4

  15. Why Molecular Simulations 5

  16. Why Molecular Simulations Paul Dirac, after completing his formalism of quantum mechanics: “ The rest is chemistry…”. 5

  17. Why Molecular Simulations Paul Dirac, after completing his formalism of quantum mechanics: “ The rest is chemistry…”. 5

  18. Why Molecular Simulations Paul Dirac, after completing his formalism of quantum mechanics: “ The rest is chemistry…”. This is a heavy burden the shoulders of “chemistry”: 5

  19. Why Molecular Simulations Paul Dirac, after completing his formalism of quantum mechanics: “ The rest is chemistry…”. This is a heavy burden the shoulders of “chemistry”: 5

  20. Why Molecular Simulations Paul Dirac, after completing his formalism of quantum mechanics: “ The rest is chemistry…”. This is a heavy burden the shoulders of “chemistry”: 5

  21. Intermolecular potential 6

  22. Intermolecular potential The intermolecular potential can: 6

  23. Intermolecular potential The intermolecular potential can: Mimic the experimental system as • accurate as possible: Replace experiments (dangerous, • impossible to measure, expensive, …) 6

  24. Intermolecular potential The intermolecular potential can: Mimic the experimental system as • accurate as possible: Replace experiments (dangerous, • impossible to measure, expensive, …) Make a model system: • Test theories that can not directly be • tested with experiment 6

  25. If we know/guess the “true” intermolecular potential

  26. Example 1: Mimic the “real world” Critical properties of long chain hydrocarbons 8

  27. Example 1: Mimic the “real world” Critical properties of long chain hydrocarbons To predict the thermodynamic properties (boiling points) of the hydrocarbon mixtures it is convenient (=Engineering models use them) to know the critical points of the hydrocarbons. 8

  28. Critical points of long chain hydrocarbons Pentane 9

  29. Critical points of long chain hydrocarbons Heptadecane Pentane 9

  30. Hydrocarbons: intermolecular potential United-atom model Fixed bond length • CH 2 CH 2 CH 3 Bond-bending • CH 2 CH 3 Torsion • Non-bonded: Lennard- • Jones 10

  31. OPLS (Jorgensen) Model

  32. Vapour-liquid 
 equilibria 12

  33. Vapour-liquid 
 equilibria Computational issues: 12

  34. Vapour-liquid 
 equilibria Computational issues: How to compute • vapour-liquid equilibrium? 12

  35. Vapour-liquid 
 equilibria Computational issues: How to compute • vapour-liquid equilibrium? How to deal • with long chain hydrocarbons? 12

  36. Vapour-liquid 
 Molecular dynamics: press enter and see … equilibria Computational issues: How to compute • vapour-liquid equilibrium? How to deal • with long chain hydrocarbons? 12

  37. Vapour-liquid 
 Molecular dynamics: press enter and see … equilibria Computational issues: How to compute • vapour-liquid equilibrium? How to deal • with long chain hydrocarbons? Molecular dynamics: press enter and see … 12

  38. Vapour-liquid 
 But my system is Molecular dynamics: press extremely small, is enter and see … equilibria the statistic reliable? Computational issues: How to compute • vapour-liquid equilibrium? How to deal • with long chain hydrocarbons? Molecular dynamics: press enter and see … 12

  39. Vapour-liquid 
 But my system is Molecular dynamics: press extremely small, is enter and see … equilibria the statistic reliable? Computational issues: How to compute • vapour-liquid equilibrium? How to deal • with long chain hydrocarbons? Molecular dynamics: But C48 moves much slower press enter and see … than methane (C1). Do I have enough CPU time 12

  40. Vapour-liquid 
 But my system is Molecular dynamics: press extremely small, is enter and see … equilibria the statistic reliable? Lectures on Free Computational issues: Energies and How to compute • Phase Equilibrium vapour-liquid equilibrium? How to deal • with long chain hydrocarbons? Molecular dynamics: But C48 moves much slower press enter and see … than methane (C1). Do I have enough CPU time 12

  41. Vapour-liquid 
 But my system is Molecular dynamics: press extremely small, is enter and see … equilibria the statistic reliable? Lectures on Free Computational issues: Energies and How to compute • Phase Equilibrium vapour-liquid Lectures on equilibrium? advanced Monte How to deal • Carlo with long chain hydrocarbons? Molecular dynamics: But C48 moves much slower press enter and see … than methane (C1). Do I have enough CPU time 12

  42. Critical Temperature and Density Nature 365 , 330 (1993). 13

  43. Example 2 Methane Storage

  44. Methane cars: the technological obstacle Gasoline, 1 liter CH 4 1 liter 0.036 MJ 34.2 MJ

  45. Methane versus gasoline LNG CNG Makal et a l. Chem. Soc. Rev. 2012 41.23, 7761-7779.

  46. 65 bar 5.8 bar Insufficient flow ~1 bar P H = 65 bar P L = 5.8 bar

  47. The deliverable capacity P L P H = 5.8 bar = 65 bar Methane adsorbed Methane adsorbed (v STP/v) (v STP/v) at tank charging at tank discharge pressure pressure ARPA-E (DOE) target: 315 m 3 STP methane/m 3 adsorbent

  48. An optimal heat of adsorption? Goal: maximize deliverable capacity

  49. An optimal heat of adsorption? Goal: maximize deliverable capacity ( ) opt = H 0 exp − q iso RT H CH 4

  50. An optimal heat of adsorption? Goal: maximize deliverable capacity ( ) opt = H 0 exp − q iso RT H CH 4

  51. An optimal heat of adsorption? Goal: maximize deliverable capacity “For methane, an optimal enthalpy change of [16.2] kJ/mol is found.” ( ) opt = H 0 exp − q iso RT H CH 4

  52. In silico screening of zeolites MFI expt’l data: Sun et al . (1998) J. Phys. Chem. B. 102(8), 1466-1473. Zhu et al. (2000) Phys. Chem. Chem. Phys. 2(9), 1989-1995. Force field: Dubbeldam et al . (2004) Phys. Rev. 93(8), 088302.

  53. In silico screening of zeolites C. Simon et al . (2014) Phys. Chem. Chem. Phys . 16 (12), 5499-5513

  54. Enthalpy vs. entropy Δ S not the same for all materials • Wide range of Δ H that yields optimal material •

  55. Can we find a material that meets the DOE target? Screening > 100,000 materials zeolites • Metal organic Frameworks, MOFs (Snurr and • co-workers) zeolitic imidazolate frameworks, ZIFs, • (Haranczyk) Polymer Porous Networks, PPNs (Haranczyk) •

  56. Insight from the model Empty tank

  57. Example 3: make a model system 26

  58. Example 3: make a model system Question: are attractive interactions needed to form a solid phase? 26

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend


More recommend