Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals
|
|
- Cleopatra O’Connor’
- 5 years ago
- Views:
Transcription
1 Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals Kevin Wang (1), Pilar Ariza (2), and Michael Ortiz (3) (1) Department of Aerospace and Ocean Engineering Virginia Polytechnic Institute and State University (2) Escuela Técnica Superior de Ingeniería Universidad de Sevilla, Spain (3) Division of Engineering and Applied Science California Institute of Technology
2 MOTIVATION Hydrogen energy in aerospace and ocean engineering Type 212 submarine (Germany) UAVs submarines and UUVs - attractive attributes long endurance near-silent operation ease of control air-independent propulsion (for submarines and UUVs) - hydrogen fuel cell vehicles would increase mission endurance and stealth, while reducing logistical burdens and costs -- office of naval research
3 HYDROGEN STORAGE Metal hydride-based hydrogen storage (MgH 2, PdH, LiH, etc.) - advantages: high volume density, safety - topic of interest: design and fabrication of nanomaterials - challenge: slow deformation-diffusion coupled process at atomic scale initial Barnoush, Hydrogen Embrittlement, 2011 after 10 cycles Li et al, JACS (2007) Mg nanowires collapse after 10 cycles of charging/discharging; length scale: 1 nm; time scale: seconds to hours!
4 MULTISCALE FRAMEWORK Time-scale gap: from femtoseconds to hours experiments simulations molecular dynamics, APPLIED OPTICS (2012) A. Lakatos et al. Vacuum (2010) ab initio MD J. Molecular Liquids (2012) - atomic level rate-limiting processes: thermal vibration, individual hops of atoms (fs) - macroscopic processes of interest: hydrogen diffusion, embrittlement, etc. (s h) - multiscale model for long-term atomistic simulation fine scale: non-equilibrium statistical thermodynamics coarse scale: discrete kinetic laws
5 NON-EQUILIBRIUM THERMODYNAMICS Statistical thermodynamics for multi-species crystals - N particles, M species occupancy function 1, particle i is from species k for i = 1,, N, k = 1,, M, n ik = 0, otherwise - assumption: separation of time scales fine scale: thermal vibration (~fs) modeled p1 p2 coarse scale: global relaxation (>>fs) solved for probability: q, p, {n} ~ p q, p, {n}, - maximum entropy principle: max p with subject to: S p = k B log p 1 f = h 3N n O Γ NM (information entropy) f q, p, {n} p q, p, {n} dqdp H h i = Ee i i = 1,, N (atomic internal energy) n N ik = n ik x ik = Ni k = 1,, N (atomic molar fraction) i=1 pa qi pi p3 Position Momentum
6 NON-EQUILIBRIUM THERMODYNAMICS Thermodynamic potentials - solve the constrained optimization problem using Lagrange multipliers: L p; β, γ = S p k B β T h + k B γ T { n } ( constrained entropy ) p q, p, n ; β, {γ} = 1 Ξ exp β T h + γ T n S = k B log ρ = k B log Ξ + β T e + γ T {x} Φ β, γ = sup L ρ; β, γ = k B log Ξ( β, γ ) ρ β i = 1 S particle temperature: T k B e i = 1 i k B β i γ ik = 1 k B S x ik (total entropy) (grand-canonical free entropy) particle chemical potential: μ ik = γ ik β i = k B T i γ ik - however, p, Ξ, S, and Φ all depend on Hamiltonian h i = h i ( p, q, n ). in general, the closed-form formulation p, Ξ, S, and Φ are intractable
7 NON-EQUILIBRIUM THERMODYNAMICS Meanfield approximations - we define trial Hamiltionians of the form h 0i q, p ; q, p, m, ω = 1 p 2 m i p 2 i + m iω i i 2 p 0 = {π} p 0 ( q, p, {n}) = 1 Ξ 0 exp argmax L[p 0, β, {γ} ] p i, q i, ω i,{ m i } uncoupled harmonic oscillators N i=1 1 β i p 2 m i p 2 i + m 2 iω i i 2 Gaussian distribution meanfield parameters: q i (mean atomic position), p i (mean atomic momentum), ω i (frequency of oscillation), m i (mean atomic mass) 2 q i q i 2 given atomic temperature {β} molar fraction {x}, the local equilibrium states { q i, p i, ω i, m i } can be uniquely determined q i q i 2
8 DISCRETE KINETIC LAWS Detailed balance of energy - balance of energy at a particle u i = w i + μ i T x i + q i v h ( w i = i α=1 A A α is the external power) α requires kinetic models for x i - assume x i can be divided into particle-to-particle fluxes of the form with - assume linear kinetics J ij = D ij (γ i γ j ) particle-to-particle mass fluxes D ij : pairwise diffusivity coefficient - enforce consistency with Fick s law for continuum D ij = D 0 2a 2 (x i + x j ) (for fcc and bcc lattices) Venturini, Wang, Romero, Ariza, Ortiz, JMPS (accepted)
9 APPLICATION Long-term atomistic simulation of mass and heat transport nanovoid growth in Cu Ariza, et al. IJ Fracture (2013) twin boundary migration in Mg Wang, Ramabathiran, Ortiz (in prep) heat conduction in SiNW Venturini, Wang, et al. JMPS (accepted)
10 APPLICATION Hydrogen diffusion in Pd nanofilms - experiment (Y. Li and Y.-T. Cheng, 1996) specimen: Pd (111) thin films with thickness L = 22 nm, 46 nm, 135 nm method: electrochemical stripping step 1: introduce H to Pd by cathodic polarization (set ΔV = 0.83V) Pd + H 2 O + e Pd H ads + OH step 2: remove H from Pd-H by anodic stripping (switch ΔV to 0.5V) Pd H ads + OH Pd + H 2 O + e l electrolyte solution Pd film substrate: Au-coated Ni measurement: time history of discharge current at the surface of Pd film C l, t I t = FS J H (l, t) = FSD H = FSD HZ x l, t l l=0 N A V cell l l=0
11 Computational model APPLICATION FCC Pd unit cell - 1D, constant temperature (300K) - meanfield positions of Pd atoms are fixed - linearization of Pd-H interatomic potential local equilibrium condition γ i = log x i 1 x i + j I H,j i C ij k B T x j + B i k B T - diffusion of hydrogen: discrete kinetic law x i = j Nei(i) D ij (γ j γ i ) x i 1D computational domain two phases!
12 MODEL VALIDATION Validation for H diffusion in α-phase PdH electolyte Dirichlet b.c. x = Pd film: x 0 = substrate Neumann b.c. flux = 0 - total simulation time: 1 sec. - time-step: ~10-6 sec. - fully atomistic Wang, Ortiz, Ariza, IJ Hydrogen Energy (submitted)
13 H molar fraction H molar fraction PHASE SEPARATION Propagation of α/β phase boundary 1 electolyte Dirichlet b.c. x = 0.99 classical model (Fick s 2 nd law) Pd film: x 0 = substrate Neumann b.c. flux = 0 new model Pd film (thickness: 46 nm) Pd film (thickness: 46 nm) Wang, Ortiz, Ariza, IJ Hydrogen Energy (submitted)
14 Velocity of phase boundary PHASE SEPARATION - velocity of phase boundary: ~100 nm/s - in many materials, phase boundaries move slowly. The prediction of phase boundary velocity is an important yet challenging task in multiscale analysis Wang, Ortiz, Ariza, IJ Hydrogen Energy (submitted)
15 SUMMARY AND FUTURE WORK Long-term atomistic simulation of mass transport - non-equilibrium statistical thermodynamics model maximum entropy principle meanfield approximation - mass transport: discrete kinetic law - coupling: a partitioned procedure - validations and applications Barnoush, (2011) Future work - large-scale, 3D simulations of hydrogen diffusion in metals - 3D validations A. Lakatos et al. Vacuum (2010) Acknowledgement: ARL, VT new faculty support (K.W.)
Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals
Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals K.G. Wang 1, M. and M. Ortiz 2 Universidad de Sevilla 1 Virginia Polytechnic Institute and State University 2 California Institute of Technology
More informationof Long-Term Transport
Atomistic Modeling and Simulation of Long-Term Transport Phenomena in Nanomaterials K.G. Wang, M. and M. Ortiz Virginia Polytechnic Institute and State University Universidad de Sevilla California Institute
More informationLong-Term Atomistic Simulation of Heat and Mass Transport
Long-Term Atomistic Simulation of Heat and Mass Transport Kevin G. Wang 1, Mauricio Ponga 2, Michael Ortiz 2 1 Department of Aerospace and Ocean Engineering, Virginia Polytechnic Institute and State University
More informationPALLADIUM HYDROGEN System
The PALLADIUM HYDROGEN System F. A. LEWIS Department of Chemistry The Queen's University of Belfast Northern Ireland ACADEMIC PBESS 1967 LONDON NEW YORK CONTENTS Preface.......... v Chapter 1. General
More informationDiffusion in multicomponent solids. Anton Van der Ven Department of Materials Science and Engineering University of Michigan Ann Arbor, MI
Diffusion in multicomponent solids nton Van der Ven Department of Materials Science and Engineering University of Michigan nn rbor, MI Coarse graining time Diffusion in a crystal Two levels of time coarse
More informationI. Collective Behavior, From Particles to Fields
I. Collective Behavior, From Particles to Fields I.A Introduction The object of the first part of this course was to introduce the principles of statistical mechanics which provide a bridge between the
More informationBridging to the Continuum Scale for Ferroelectric Applications
Bridging to the Continuum Scale for Ferroelectric Applications Shanfu Zheng and Alberto Cuitiño Mechanical and Aerospace Engineering, Rutgers University Alejandro Strachan Materials Engineering, Purdue
More informationSimilarities and differences:
How does the system reach equilibrium? I./9 Chemical equilibrium I./ Equilibrium electrochemistry III./ Molecules in motion physical processes, non-reactive systems III./5-7 Reaction rate, mechanism, molecular
More informationCollective behavior, from particles to fields
978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the
More informationATOMISTIC MODELING AND SIMULATION OF LONG-TERM TRANSPORT PHENOMENA IN NANOMATERIALS
IS Atomistic - Evolving modeling Microstructures and simulation in Solids of long-term transport phenomena in nanomaterials XIII International Conference on Computational Plasticity. Fundamentals and Applications
More informationfiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES
Content-Thermodynamics & Statistical Mechanics 1. Kinetic theory of gases..(1-13) 1.1 Basic assumption of kinetic theory 1.1.1 Pressure exerted by a gas 1.2 Gas Law for Ideal gases: 1.2.1 Boyle s Law 1.2.2
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
More informationIdeal Gas Behavior. NC State University
Chemistry 331 Lecture 6 Ideal Gas Behavior NC State University Macroscopic variables P, T Pressure is a force per unit area (P= F/A) The force arises from the change in momentum as particles hit an object
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationIn-class exercises. Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 8 Exercises due Mon March 19 Last correction at March 5, 2018, 8:48 am c 2017, James Sethna,
More informationEntropy as an Indication of Damage in Engineering Materials
Entropy as an Indication of Damage in Engineering Materials Mohammad Modarres Presented at Entropy 2018: From Physics to Information Sciences and Geometry Barcelona, Spain, May 14-16, 2018 16 May 2018
More informationConcepts in Surface Physics
M.-C. Desjonqueres D. Spanjaard Concepts in Surface Physics Second Edition With 257 Figures Springer 1. Introduction................................. 1 2. Thermodynamical and Statistical Properties of
More informationApplications to simulations: Monte-Carlo
Applications to simulations: Monte-Carlo A.C. Maggs ESPCI, Paris June 2016 Summary Algorithms Faster/simpler codes Thermodynamics of Electric fields Partition function of electric field Thermal Casimir/Lifshitz
More informationUniversity of Ljubljana. Faculty of Mathematics and Physics. Department of Physics. High-entropy alloys. Author: Darja Gačnik
University of Ljubljana Faculty of Mathematics and Physics Department of Physics Seminar I 1 st year of Master study programme High-entropy alloys Author: Darja Gačnik Mentor: prof. Janez Dolinšek Ljubljana,
More information(a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron, but neglecting spin-orbit interactions.
1. Quantum Mechanics (Spring 2007) Consider a hydrogen atom in a weak uniform magnetic field B = Bê z. (a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron,
More informationDislocations in graphene
Dislocations in graphene M. Ortiz California Institute of Technology In collaboration with: M.P. Ariza, Universidad de Sevilla Symposium on Multiscale Dislocation Dynamics UCSD, La Jolla, January 16-17,
More informationEntropic methods to study the evolution of damage and degradation of materials
Entropic methods to study the evolution of damage and degradation of materials Mohammad Modarres Presented at 14th International Conference on Fracture Rhodes, Greece, June 18-23, 2017 20 JUNE 2017 Department
More informationENGN2210 CONTINUUM MECHANICS
School of Engineering Brown University ENGN2210 CONTINUUM MECHANICS Allan Bower Fall 2016 What is continuum mechanics? Boeing Dreamliner Chevy Volt What is continuum mechanics? A continuum idealized form
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
More informationThermodynamics of Damage: A Reliability Engineering Perspective
Thermodynamics of Damage: A Reliability Engineering Perspective Mohammad Modarres Presented at Department of Materials Science and Engineering Seminar UCLA, Los Angeles CA February 17, 217 Department of
More informationI. BASICS OF STATISTICAL MECHANICS AND QUANTUM MECHANICS
I. BASICS OF STATISTICAL MECHANICS AND QUANTUM MECHANICS Markus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus/cours.html
More informationDensity Functional Modeling of Nanocrystalline Materials
Density Functional Modeling of Nanocrystalline Materials A new approach for modeling atomic scale properties in materials Peter Stefanovic Supervisor: Nikolas Provatas 70 / Part 1-7 February 007 Density
More information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More informationThe Microcanonical Approach. (a) The volume of accessible phase space for a given total energy is proportional to. dq 1 dq 2 dq N dp 1 dp 2 dp N,
8333: Statistical Mechanics I Problem Set # 6 Solutions Fall 003 Classical Harmonic Oscillators: The Microcanonical Approach a The volume of accessible phase space for a given total energy is proportional
More informationDIFFUSION IN SOLIDS. IE-114 Materials Science and General Chemistry Lecture-5
DIFFUSION IN SOLIDS IE-114 Materials Science and General Chemistry Lecture-5 Diffusion The mechanism by which matter is transported through matter. It is related to internal atomic movement. Atomic movement;
More informationSteady-state diffusion is diffusion in which the concentration of the diffusing atoms at
Chapter 7 What is steady state diffusion? Steady-state diffusion is diffusion in which the concentration of the diffusing atoms at any point, x, and hence the concentration gradient at x, in the solid,
More informationNanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons
Nanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons Gang Chen Massachusetts Institute of Technology OXFORD UNIVERSITY PRESS 2005 Contents Foreword,
More informationStatistical Mechanics. Atomistic view of Materials
Statistical Mechanics Atomistic view of Materials What is statistical mechanics? Microscopic (atoms, electrons, etc.) Statistical mechanics Macroscopic (Thermodynamics) Sample with constrains Fixed thermodynamics
More informationElectrochemical System
Electrochemical System Topic Outcomes Week Topic Topic Outcomes 8-10 Electrochemical systems It is expected that students are able to: Electrochemical system and its thermodynamics Chemical reactions in
More informationContents. 2. Fluids. 1. Introduction
Contents 1. Introduction 2. Fluids 3. Physics of Microfluidic Systems 4. Microfabrication Technologies 5. Flow Control 6. Micropumps 7. Sensors 8. Ink-Jet Technology 9. Liquid Handling 10.Microarrays 11.Microreactors
More informationDiffusion in Dilute Alloys
Chapter 3 Diffusion in Dilute Alloys Our discussion of the atomistic mechanisms for diffusion has been confined to the situations where the diffusing species is chemically identical to the host atom. We
More informationMolecular Driving Forces
Molecular Driving Forces Statistical Thermodynamics in Chemistry and Biology SUBGfittingen 7 At 216 513 073 / / Ken A. Dill Sarina Bromberg With the assistance of Dirk Stigter on the Electrostatics chapters
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationSolutions for Assignment-6
Solutions for Assignment-6 Q1. What is the aim of thin film deposition? [1] (a) To maintain surface uniformity (b) To reduce the amount (or mass) of light absorbing materials (c) To decrease the weight
More informationCombined atomistic-continuum modeling of short-pulse laser melting and disintegration of metal films
Combined atomistic-continuum modeling of short-pulse laser melting and disintegration of metal films Dmitriy S. Ivanov and Leonid V. Zhigilei* Department of Materials Science and Engineering, University
More informationAssigning Oxidation Numbers
RULES Assigning Oxidation Numbers Examples 1. Each Uncombined Element has an 2Na + Cl 2 2NaCl oxidation number = 0 Na = 0 or written Na 0 Cl 2 = 0 or written Cl 0 Monatomic ions have an oxidation number
More informationMolecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao gaofeium@umich.edu Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationModeling the next battery generation: Lithium-sulfur and lithium-air cells
Modeling the next battery generation: Lithium-sulfur and lithium-air cells D. N. Fronczek, T. Danner, B. Horstmann, Wolfgang G. Bessler German Aerospace Center (DLR) University Stuttgart (ITW) Helmholtz
More informationFast Adiabatic Control of a Harmonic Oscillator s Frequency
Fast Adiabatic Control of a Harmonic Oscillator s Frequency Peter Salamon Dept. of Math & Stats San Diego State University KITP; June, 2009 Ensemble of Independent Harmonic Oscillators Sharing a Controlled
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature17653 Supplementary Methods Electronic transport mechanism in H-SNO In pristine RNO, pronounced electron-phonon interaction results in polaron formation that dominates the electronic
More informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationMultiscale modeling and simulation of shock wave propagation
University of Iowa Iowa Research Online Theses and Dissertations Spring 2013 Multiscale modeling and simulation of shock wave propagation Cory Nelsen University of Iowa Copyright 2013 Cory Nelsen This
More informationGear methods I + 1/18
Gear methods I + 1/18 Predictor-corrector type: knowledge of history is used to predict an approximate solution, which is made more accurate in the following step we do not want (otherwise good) methods
More informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationReactive potentials and applications
1.021, 3.021, 10.333, 22.00 Introduction to Modeling and Simulation Spring 2011 Part I Continuum and particle methods Reactive potentials and applications Lecture 8 Markus J. Buehler Laboratory for Atomistic
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationElectronic Squeezing by Optically Pumped Phonons: Transient Superconductivity in K 3 C 60. With: Eli Wilner Dante Kennes Andrew Millis
Electronic Squeezing by Optically Pumped Phonons: Transient Superconductivity in K 3 C 60 With: Eli Wilner Dante Kennes Andrew Millis Background: Mean-Field Theory of Simple Superconductivity If the effective
More informationIAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.
IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm
More informationi i ne. (1) i The potential difference, which is always defined to be the potential of the electrode minus the potential of the electrolyte, is ln( a
We re going to calculate the open circuit voltage of two types of electrochemical system: polymer electrolyte membrane (PEM) fuel cells and lead-acid batteries. To do this, we re going to make use of two
More informationELECTROCHEMICAL CELLS NAME ROW PD
4-26-12 NAME ROW PD (1) Which statement describes the redox reaction that occurs when an object is electroplated? The diagram below shows the electrolysis of fused KCl. A) It is spontaneous and requires
More informationBasic overall reaction for hydrogen powering
Fuel Cell Basics Basic overall reaction for hydrogen powering 2H 2 + O 2 2H 2 O Hydrogen produces electrons, protons, heat and water PEMFC Anode reaction: H 2 2H + + 2e Cathode reaction: (½)O 2 + 2H +
More informationCalculation of thermal expansion coefficient of Fe3Al with the addition of transition metal elements Abstract Introduction Methodology
Calculation of thermal expansion coefficient of Fe Al with the addition of transition metal elements Tatiana Seletskaia, Leonid Muratov, Bernard Cooper, West Virginia University, Physics Departement, Hodges
More information6. Computational Design of Energy-related Materials
6. Computational Design of Energy-related Materials Contents 6.1 Atomistic Simulation Methods for Energy Materials 6.2 ab initio design of photovoltaic materials 6.3 Solid Ion Conductors for Fuel Cells
More informationAn Atomistic-based Cohesive Zone Model for Quasi-continua
An Atomistic-based Cohesive Zone Model for Quasi-continua By Xiaowei Zeng and Shaofan Li Department of Civil and Environmental Engineering, University of California, Berkeley, CA94720, USA Extended Abstract
More informationMesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai
Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationPhonon wavefunctions and electron phonon interactions in semiconductors
Phonon wavefunctions and electron phonon interactions in semiconductors Bartomeu Monserrat bm418@cam.ac.uk University of Cambridge Quantum Monte Carlo in the Apuan Alps VII QMC in the Apuan Alps VII Bartomeu
More informationJavier Junquera. Statistical mechanics
Javier Junquera Statistical mechanics From the microscopic to the macroscopic level: the realm of statistical mechanics Computer simulations Thermodynamic state Generates information at the microscopic
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationBatteries (Electrochemical Power Sources)
Batteries (Electrochemical Power Sources) 1. Primary (single-discharge) batteries. => finite quantity of the reactants 2. Secondary or rechargeable batteries => regeneration of the original reactants by
More informationPolarons. University of Ljubljana Faculty of Mathematics and Physics. Seminar - Ib, 1. year II. cycle degrees
University of Ljubljana Faculty of Mathematics and Physics Seminar - Ib, 1. year II. cycle degrees Polarons Author: Jaka Vodeb Advisor: prof. dr. Viktor Kabanov Kranj, 017 Abstract After the discovery
More informationElectrochem 1 Electrochemistry Some Key Topics Conduction metallic electrolytic Electrolysis effect and stoichiometry Galvanic cell Electrolytic cell Electromotive Force (potential in volts) Electrode
More informationDynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics
for resubmission Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics Aram Davtyan, 1 Gregory A. Voth, 1 2, a) and Hans
More informationMOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM
MOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM MOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM By YUNFEI BAI, B. Eng.
More informationNANO/MICROSCALE HEAT TRANSFER
NANO/MICROSCALE HEAT TRANSFER Zhuomin M. Zhang Georgia Institute of Technology Atlanta, Georgia New York Chicago San Francisco Lisbon London Madrid Mexico City Milan New Delhi San Juan Seoul Singapore
More informationLectures on: Introduction to and fundamentals of discrete dislocations and dislocation dynamics. Theoretical concepts and computational methods
Lectures on: Introduction to and fundamentals of discrete dislocations and dislocation dynamics. Theoretical concepts and computational methods Hussein M. Zbib School of Mechanical and Materials Engineering
More informationQuantum Spin-Metals in Weak Mott Insulators
Quantum Spin-Metals in Weak Mott Insulators MPA Fisher (with O. Motrunich, Donna Sheng, Simon Trebst) Quantum Critical Phenomena conference Toronto 9/27/08 Quantum Spin-metals - spin liquids with Bose
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationElectrochemical Cell - Basics
Electrochemical Cell - Basics The electrochemical cell e - (a) Load (b) Load e - M + M + Negative electrode Positive electrode Negative electrode Positive electrode Cathode Anode Anode Cathode Anode Anode
More informationName: Regents Chemistry Date:
Name: Date: 1. The reaction CuO + CO CO 2 + Cu is an example of (A) reduction, only (B) oxidation, only (C) both oxidation and reduction (D) neither oxidation nor reduction 6. In which compound does chlorine
More informationWe can use chemistry to generate electricity... this is termed a Voltaic (or sometimes) Galvanic Cell
Unit 6 Electrochemistry Chemistry 020, R. R. Martin Electrochemistry Electrochemistry is the study of the interconversion of electrical and chemical energy. We can use chemistry to generate electricity...
More informationInterpolation. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Interpolation 1 / 34
Interpolation Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Interpolation 1 / 34 Outline 1 Introduction 2 Lagrange interpolation 3 Piecewise polynomial
More informationII Relationship of Classical Theory to Quantum Theory A Quantum mean occupation number
Appendix B Some Unifying Concepts Version 04.AppB.11.1K [including mostly Chapters 1 through 11] by Kip [This appendix is in the very early stages of development] I Physics as Geometry A Newtonian Physics
More informationFemtosecond laser microfabrication in. Prof. Dr. Cleber R. Mendonca
Femtosecond laser microfabrication in polymers Prof. Dr. Cleber R. Mendonca laser microfabrication focus laser beam on material s surface laser microfabrication laser microfabrication laser microfabrication
More informationLecture 2+3: Simulations of Soft Matter. 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications
Lecture 2+3: Simulations of Soft Matter 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications D. Frenkel, Boulder, July 6, 2006 What distinguishes Colloids from atoms or
More informationTransport (kinetic) phenomena: diffusion, electric conductivity, viscosity, heat conduction...
Transport phenomena 1/16 Transport (kinetic) phenomena: diffusion, electric conductivity, viscosity, heat conduction... Flux of mass, charge, momentum, heat,...... J = amount (of quantity) transported
More informationModèle de liaisons fortes au 4ème moment pour traiter l ordre-désordre dans les alliages
Modèle de liaisons fortes au 4ème moment pour traiter l ordre-désordre dans les alliages Jan Los, Christine Mottet, Guy Tréglia CINaM, Marseille Christine Goyhenex IPCMS, Strasbourg Outline Context Tight
More informationStatistical Mechanics
Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2
More informationElectrochemical Cells Homework Unit 11 - Topic 4
Electrochemical Cell Vocabulary Electrochemical Cells Homework Unit 11 - Topic 4 Electrode Anode Voltaic Cell Oxidation Electrolytic Cell Cathode Salt Bridge Reduction Half Reaction Refer to Table J: For
More informationYiping Zhao and Yuping He
Yiping Zhao and Yuping He Department of Physics and Astronomy Nanoscale Science and Engineering Center The University of Georgia, Athens, GA 30602 April 29, 2010 Outline: Introduction Hydrogenation behaviors
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More informationCoarse-graining molecular dynamics in crystalline solids
in crystalline solids Xiantao Li Department of Mathematics Penn State University. July 19, 2009. General framework Outline 1. Motivation and Problem setup 2. General methodology 3. The derivation of the
More informationMODELING EXCESS HEAT IN THE FLEISCHMANN-PONS EXPERIMENT
Theory (I) Session 2 O_1 MODELING EXCESS HEAT IN THE FLEISCHMANN-PONS EXPERIMENT Peter L. Hagelstein 1 and Irfan U. Chaudhary 2 1 Massachusetts Institute of Technology, Cambridge, MA, USA 2 University
More informationIntroduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012
Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 K. Kremer Max Planck Institute for Polymer Research, Mainz Overview Simulations, general considerations
More informationStructure of Surfaces
Structure of Surfaces C Stepped surface Interference of two waves Bragg s law Path difference = AB+BC =2dsin ( =glancing angle) If, n =2dsin, constructive interference Ex) in a cubic lattice of unit cell
More informationPREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures
PREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures Kevin G. Wang Virginia Tech Patrick Lea, Alex Main, Charbel Farhat Stanford University
More information(a) What are the probabilities associated with finding the different allowed values of the z-component of the spin after time T?
1. Quantum Mechanics (Fall 2002) A Stern-Gerlach apparatus is adjusted so that the z-component of the spin of an electron (spin-1/2) transmitted through it is /2. A uniform magnetic field in the x-direction
More informationThermal and Statistical Physics Department Exam Last updated November 4, L π
Thermal and Statistical Physics Department Exam Last updated November 4, 013 1. a. Define the chemical potential µ. Show that two systems are in diffusive equilibrium if µ 1 =µ. You may start with F =
More informationChemical reactions as network of rare events: Kinetic MonteCarlo
Chemical reactions as network of rare events: Kinetic MonteCarlo Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum
More informationChapter 18. Electrochemistry
Chapter 18 Electrochemistry Section 17.1 Spontaneous Processes and Entropy Section 17.1 http://www.bozemanscience.com/ap-chemistry/ Spontaneous Processes and Entropy Section 17.1 Spontaneous Processes
More informationJournal of the Mechanics and Physics of Solids
Accepted Manuscript Diffusive molecular dynamics simulations of lithiation of silicon nanopillars J.P. Méndez, M. Ponga, M. Ortiz PII: S0022-5096(17)30880-3 DOI: 10.1016/j.jmps.2018.03.008 Reference: MPS
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More information9/19/2018. Corrosion Thermodynamics 2-3. Course Outline. Guiding Principles. Why study thermodynamics? Guiding Principles
Kwame Nkrumah University of Science & Technology, Kumasi, Ghana Week 1 Course Outline Topic Introduction: Reactivity types, corrosion definition, atmospheric corrosion, classification, effects, costs,
More information