1 Bulk Simulations in a HCP lattice
|
|
- Berniece McBride
- 5 years ago
- Views:
Transcription
1 1 Bulk Simulations in a HCP lattice 1.1 Introduction This project is a continuation of two previous projects that studied the mechanism of melting in a FCC and BCC lattice. The current project studies the melting process in a HCP lattice which is different than the previous two in some aspects. The method in which is employed in this study is molecular dynamics, this demands a search for a stable potential which can describe good physical behavior and results compared with experimental results and physical behavior. A potential was found for magnesium which has a HCP structure at no external press, the potential was constructed especially for temperatures close to the melting temperature. The potential is part of a family potential called an embedded atom potential, the potential has two parts, one is a pair potential like a Lennard Jones potential and the second is an embedded atom function. The latter holds the effective potential that an atom feels as a result from its surroundings. The first step in a melting simulation is to test the potential, it was chosen to check the thermal expansion and the elastic constants. It is expected that the thermal expansion 1 dl α = LdT will rise slightly with temperature and will have the same order of magnitude like the experimental thermal expansion. This formula is the thermal expansion along one side of the simulated box. From the elastic constants it is expect the decrease with temperature and of course to have the same order of magnitude as the experimental ones. The simulation is run in the NtT ensemble isothermal isobaric ensemble. In this ensemble the number of atoms is held fix, the external temperature and external pressure. This ensemble is chosen so that with the desired temperature and pressure held fixed the simulated box can expand from ground state (zero temperature) to its relaxed state with respect to the external variables. When the system is relaxed the lattice constants can be calculated, they are needed later for the surface simulations where it requires working in a NVT ensemble, V stands for a constant volume that is taken from the NtT simulation by measuring the lattice constants. 1.2 HCP lattice This project is different from the previous two because of the geometric structure of magnesium and the different symmetry compared with FCC and BCC. BCC and FCC both have a cubic symmetry, they belong to a class of structures called Bravais lattice which means that a crystal can be described by using only 3 vectors and their products with an integer. A HCP lattice is not a Bravais lattice and it does not have a cubic
2 symmetry. HCP stand for Hexagonal Closed Packed which is a Bravais lattice of hexagonal with a basis where the base composed of two atoms. Hexagonal lattice have 6-fold symmetry in the XY plane and a 2-fold symmetry in the Z axis. One other difference is that FCC and BCC lattice are defined by only one lattice constant while HCP needs two lattice constants. In the XY plane the atoms minimal distance is the lattice constant ' a ', the XY planes are separated by a lattice constant ' c 2'. All of those differences complicate the understanding of the system as well as the ability to model it correctly and do simulations. 1.3 NtT ensemble The algorithm that produces the NtT ensemble keeps the system coupled to two external baths, the heat bath and the pressure bath, by introducing to the equation of motion, of a particle, imaginary non physical forces. Each bath has it one imaginary force which after a few iterations is periodic in time, this periodic depends on the bath's mass. The baths mass is an arbitrary constant which should be chosen wisely, in the end it doesn't have any effect on the average values but to get the correct averaged values many periods should be performed which sums at the end in the total running time of one simulation. In this algorithm the simulated box is describes with 3 time dependent vectors abc,, defined from the origin of the coordinates. These 3 vectors form a 3 3 matrix h det h = a b c as a result h components are regarded as dynamical such that ( ) variables describing changes in the volume and shape of the simulated box. This procedure scales the atoms position into a simulated unit box where each atom coordinate r i is scaled down into a scaled coordinate s i via a realationship r i = hs. i The Lagrangian for the TtN ensemble T ( ) ( ) 1 2 T 1 L = mf i s i Gs i ϕ ri rj + WTrace h h 2 i i< j VTr 0 ( σε ) + Qf ( 3N+ 1) kbt0ln f 2 where QW, are the masses of the heat and pressure baths respectively, m is the atom's mass, σ is the external stress tensor, k B is Boltzmann constant, T 0 is the heat T bath's temperature, G = h h V = det h is the initial is called the metric tensor, ( ) T volume of the simulated box, N is the number of atom, ε ( h ) strain tensor. The equations of motion can be derived from the Lagrangian: 1 1 ( Gh 1) 0 0 = 2 is the
3 where p int ( ij ) 1 ϕ r s f T 1 T h = pintv ( h ) V0hh0 σ ( h 0 ) W 1 T kt B 0 f = mi fs i Gs i ( 3N + 1) Q i f ij 1 si = 2 G G+ si mf i j rij rij f is the microscopic internal stress tensor ( r ) ϕ p m f hs hs r r α, β 1 2 α β ij α β int = i i i ij ij V i i< j rij 1.4 Integrator The simulation uses a predictor-corrector integrator to solve the equation of motion. 1.5 Stages in the simulation The simulation start by initializing all of the atoms in their ground state, then it gives them random initial velocities depending on the temperature and Maxwell's distribution. The simulation than creates for each atom an array of 12 for its neighbors, this list is constant throughout the entire simulation, later the list is used for the lattice constants calculations. The simulation every 100 time steps, each time step is femtosecond, updates a neighbor list for each atom with respect to some cutoff radius, it then used in the calculation of the forces acting on that atom. The order of the simulation: 1. Building the neighbor list for each atom. 2. Computation of the forces on each atom. 3. Computation of the value of f. 4. Computation the element matrix of h. 5. Prediction of the next position and velocities of each atom. 6. Prediction of the next values of the f and f. 7. Prediction of the next values of the matrices h and h. 8. Computation of the forces on each atom calculated from the potential. 9. Computation of the value of f. 10. Computation the element matrix of h. 11. Adding to the forces acting on each atom the correction terms (the second and third term in s ). 12. Correction of the next position and velocities of each atom. 13. Correction of the next values of the f and f.
4 14. Correction of the next values of the matrices h and h. 15. Repetition of steps 5 to 12 until the total number of step is executed. Every 100 time steps the simulation calculates and writes the following in files: 1. The diagonal elements of the matrix h. 2. The off-diagonal elements of the matrix h. 3. Kinetic energy, potential energy, current temperature and pressure. 4. Lattice constants and volume of the simulated box. 5. Order parameters, one describing the XY plane and the second the Z axis. 1.6 Calculation of lattice constants and thermal expansion The lattice constants are calculated from the diagonal elements of the averaged h matrix over a long period of time. The spacing between XY planes is c = h 33 avg 2 ( N ) where N Atom is defined at the beginning of the simulation, it is the number of atoms in each direction of the simulation box. Atom The minimal distance between atoms in the XY plane is 1 h11 2 h22 avg avg a = + 2 ( NAtom 2) 3 ( NAtom 2) The factor 2 3 is a result from the vectors that composite the XY plane: a 3 a1 = [ a 0 0 ] ; a2 = a The calculated lattice constants enable the calculation of their thermal expansion: 1 dl 2 ΔL 2 1 N Atom N Atom α = = = a( T +ΔT) a( T) L T dt N a T ΔT N a T ΔT 2 2 ( ) ( ) ( ) Atom 1 1 α = +Δ a T ΔT ( ) Atom ( a( T T) a( T) )
5 2 Results All results are averaged starting from time step 4000 and ending at time step Figure 1: fluctuation of the matrix element h 11 as a function of time steps in 5 different temperatures
6 Figure 2: fluctuation of the matrix element h 22 as a function of time steps in 5 different temperatures Figure 3: fluctuation of the matrix element h 33 as a function of time steps in 5 different temperatures Figure 4: lattice constants as a function of temperature
7 Figure 5: thermal expansion of the lattice constants Figure 6: lattice constants as a function of temperature, each point is an average over 1000 time steps starting from step number The solid line connects between the averages of all the point in the same temperature
8 Figure 7: thermal expansion of the lattice constants from figure 5 vs. thermal expansion as a result of the segment division. Figure 8: fluctuation of the matrix element h 11 as a function of time steps, 5 different runs in the same temperature T = K
9 Figure 9: fluctuation of the matrix element h 22 as a function of time steps, 5 different runs in the same temperature T = K Figure 10: fluctuation of the matrix element h 33 as a function of time steps, 5 different runs in the same temperature T = K
10 Figure 11: average of each lattice constant for every one of the runs, the title contents the average of the averages and its standard error deviation. Figure 12: figure 4 vs. lattice constants as a function of temperature with new value at T = K
11 Figure 13: figure 5 vs. thermal expansion of the lattice constants with new value at T = K a(t) [Angstrom] c(t) [Angstrom] temperature [K] temperature [K] Figure 14: lattice constants as a function of temperature, each point is an average over 1000 time steps starting from step number The solid line connects between the averages of all the point in the same temperature. This figure includes all the runs for T = K
12 2.48 x 104 Volume expansion as a function of temperature Volume(T) [Angstrom 3 ] Volume(T) [Angstrom 3 ] temperature [K] 2.46 x Volume at T=K as a function of run number run # Figure 15:
What is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationCrystal Structure. Dr Bindu Krishnan
Solid State Physics-1 Crystal Structure Dr Bindu Krishnan CRYSTAL LATTICE What is crystal (space) lattice? In crystallography, only the geometrical properties of the crystal are of interest, therefore
More informationPhys 412 Solid State Physics. Lecturer: Réka Albert
Phys 412 Solid State Physics Lecturer: Réka Albert What is a solid? A material that keeps its shape Can be deformed by stress Returns to original shape if it is not strained too much Solid structure
More informationIntroduction to Solid State Physics or the study of physical properties of matter in a solid phase
Introduction to Solid State Physics or the study of physical properties of matter in a solid phase Prof. Germar Hoffmann 1. Crystal Structures 2. Reciprocal Lattice 3. Crystal Binding and Elastic Constants
More informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
More informationFrom Atoms to Materials: Predictive Theory and Simulations
From Atoms to Materials: Predictive Theory and Simulations Week 3 Lecture 4 Potentials for metals and semiconductors Ale Strachan strachan@purdue.edu School of Materials Engineering & Birck anotechnology
More informationMSE 201A Thermodynamics and Phase Transformations Fall, 2008 Problem Set No. 7
MSE 21A Thermodynamics and Phase Transformations Fall, 28 Problem Set No. 7 Problem 1: (a) Show that if the point group of a material contains 2 perpendicular 2-fold axes then a second-order tensor property
More informationCrystallographic structure Physical vs Chemical bonding in solids
Crystallographic structure Physical vs Chemical bonding in solids Inert gas and molecular crystals: Van der Waals forces (physics) Water and organic chemistry H bonds (physics) Quartz crystal SiO 2 : covalent
More informationLecture contents. Stress and strain Deformation potential. NNSE 618 Lecture #23
1 Lecture contents Stress and strain Deformation potential Few concepts from linear elasticity theory : Stress and Strain 6 independent components 2 Stress = force/area ( 3x3 symmetric tensor! ) ij ji
More informationHomework 4 Due 25 October 2018 The numbers following each question give the approximate percentage of marks allocated to that question.
Name: Homework 4 Due 25 October 218 The numbers following each question give the approximate percentage of marks allocated to that question. 1. Use the reciprocal metric tensor again to calculate the angle
More informationBulk Structures of Crystals
Bulk Structures of Crystals 7 crystal systems can be further subdivided into 32 crystal classes... see Simon Garrett, "Introduction to Surface Analysis CEM924": http://www.cem.msu.edu/~cem924sg/lecturenotes.html
More informationChapter 1. Crystal structure. 1.1 Crystal lattices
Chapter 1 Crystal structure 1.1 Crystal lattices We will concentrate as stated in the introduction, on perfect crystals, i.e. on arrays of atoms, where a given arrangement is repeated forming a periodic
More informationChem 728 Introduction to Solid Surfaces
Chem 728 Introduction to Solid Surfaces Solids: hard; fracture; not compressible; molecules close to each other Liquids: molecules mobile, but quite close to each other Gases: molecules very mobile; compressible
More informationProblem with Kohn-Sham equations
Problem with Kohn-Sham equations (So much time consuming) H s Ψ = E el Ψ ( T + V [ n] + V [ n] + V [ n]) ϕ = Eϕ i = 1, 2,.., N s e e ext XC i i N nr ( ) = ϕi i= 1 2 The one-particle Kohn-Sham equations
More informationPre-yield non-affine fluctuations and a hidden critical point in strained crystals
Supplementary Information for: Pre-yield non-affine fluctuations and a hidden critical point in strained crystals Tamoghna Das, a,b Saswati Ganguly, b Surajit Sengupta c and Madan Rao d a Collective Interactions
More informationPractice Problems Set II
P1. For the HCP crystal structure, (a) show that the ideal c/a ratio is 1.633; (b) show that the atomic packing factor for HCP is 0.74. (a) A sketch of one-third of an HCP unit cell is shown below. Consider
More information2. Surface geometric and electronic structure: a primer
2. Surface geometric and electronic structure: a primer 2.1 Surface crystallography 2.1.1. Crystal structures - A crystal structure is made up of two basic elements: lattice + basis Basis: Lattice: simplest
More informationCRYSTAL STRUCTURES WITH CUBIC UNIT CELLS
CRYSTAL STRUCTURES WITH CUBIC UNIT CELLS Crystalline solids are a three dimensional collection of individual atoms, ions, or whole molecules organized in repeating patterns. These atoms, ions, or molecules
More information3-D Crystal Lattice Images
3-D Crystal Lattice Images All of the following images are crossed-stereo pairs. To view them, cross your eyes and focus. Author's note this material has been expanded and updated, and can be found at
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 informationUNIT I SOLID STATE PHYSICS
UNIT I SOLID STATE PHYSICS CHAPTER 1 CRYSTAL STRUCTURE 1.1 INTRODUCTION When two atoms are brought together, two kinds of forces: attraction and repulsion come into play. The force of attraction increases
More informationChapter 2. Atomic Packing
Chapter 2. Atomic Packing Contents 2-1. Packing of directional bonding atoms 2-2. Packing of indirectional bonding in same size atoms 2-3. Packing of indirectional bonding in different size atoms 2-4.
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More information... 3, , = a (1) 3 3 a 2 = a (2) The reciprocal lattice vectors are defined by the condition a b = 2πδ ij, which gives
PHZ646: Fall 013 Problem set # 4: Crystal Structure due Monday, 10/14 at the time of the class Instructor: D. L. Maslov maslov@phys.ufl.edu 39-0513 Rm. 114 Office hours: TR 3 pm-4 pm Please help your instructor
More informationAtomic Arrangement. Primer Materials For Science Teaching Spring
Atomic Arrangement Primer Materials For Science Teaching Spring 2016 31.3.2015 Levels of atomic arrangements No order In gases, for example the atoms have no order, they are randomly distributed filling
More informationSurface and Interface Science Physics 627; Chemistry 542. Lectures 3 Feb 1, Surface Structure
Surface and Interface Science Physics 67; Chemistry 54 Lectures 3 Fe, 03 Surface Structure References: ) Zangwill, Pp. 8-3 ) Woodruff & Delchar, Chapter 3) Masel, Chapter 4) Ertl & Kuppers, 0-07 5) Luth,
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 informationMonte Carlo Study of the Crystalline and Amorphous NaK Alloy
Available online at www.sciencedirect.com ScienceDirect Procedia Computer Science 08C (207) 25 22 International Conference on Computational Science, ICCS 207, 2-4 June 207, Zurich, Switzerland Monte Carlo
More informationRustam Z. Khaliullin University of Zürich
Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation
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 informationDepartment of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an
The statistical characteristics of static friction J. Wang, G. F. Wang*, and W. K. Yuan Department of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an 710049, China * E-mail: wanggf@mail.xjtu.edu.cn
More informationn-dimensional, infinite, periodic array of points, each of which has identical surroundings.
crystallography ll Lattice n-dimensional, infinite, periodic array of points, each of which has identical surroundings. use this as test for lattice points A2 ("bcc") structure lattice points Lattice n-dimensional,
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationClass 29: Reciprocal Space 3: Ewald sphere, Simple Cubic, FCC and BCC in Reciprocal Space
Class 29: Reciprocal Space 3: Ewald sphere, Simple Cubic, FCC and BCC in Reciprocal Space We have seen that diffraction occurs when, in reciprocal space, Let us now plot this information. Let us designate
More informationAtomic Arrangement. Primer in Materials Spring
Atomic Arrangement Primer in Materials Spring 2017 30.4.2017 1 Levels of atomic arrangements No order In gases, for example the atoms have no order, they are randomly distributed filling the volume to
More informationAtoms, electrons and Solids
Atoms, electrons and Solids Shell model of an atom negative electron orbiting a positive nucleus QM tells that to minimize total energy the electrons fill up shells. Each orbit in a shell has a specific
More informationT. Egami. Model System of Dense Random Packing (DRP)
Introduction to Metallic Glasses: How they are different/similar to other glasses T. Egami Model System of Dense Random Packing (DRP) Hard Sphere vs. Soft Sphere Glass transition Universal behavior History:
More informationMolecular Dynamics Simulation of Argon
Molecular Dynamics Simulation of Argon The fundamental work on this problem was done by A. Rahman, Phys. Rev. 136, A405 (1964). It was extended in many important ways by L. Verlet, Phys. Rev. 159, 98 (1967),
More informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
More informationPhys 460 Describing and Classifying Crystal Lattices
Phys 460 Describing and Classifying Crystal Lattices What is a material? ^ crystalline Regular lattice of atoms Each atom has a positively charged nucleus surrounded by negative electrons Electrons are
More informationProject 4/5 - Molecular dynamics part II: advanced study of Lennard-Jones fluids, deadline December 1st (noon)
Format for delivery of report and programs The format of the project is that of a printed file or hand-written report. The programs should also be included with the report. Write only your candidate number
More informationLab 1: Handout GULP: an Empirical energy code
3.320/SMA 5.107/ Atomistic Modeling of Materials Spring 2003 1 Lab 1: Handout GULP: an Empirical energy code We will be using the GULP code as our energy code. GULP is a program for performing a variety
More informationIndex. A Ab-initio, 8 Accomodation, 7 Accuracy, 37 hardware, 37 software, 37 Ackland, 36, 55, 157 Andersen, 44, 47 Austenite, 3
Index A Ab-initio, 8 Accomodation, 7 Accuracy, 37 hardware, 37 software, 37 Ackland, 36, 55, 157 Andersen, 44, 47 Austenite, 3 B B2, 151 154, 156 Bain, 9, 10, 72, 73 parameters, 9 transformation, 9, 72
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 informationSupplementary Materials
Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials
More informationWe need to be able to describe planes and directions.
We need to be able to describe planes and directions. Miller Indices & XRD 1 2 Determining crystal structure and identifying materials (B) Plastic deformation Plastic deformation and mechanical properties
More informationEGN 3365 Review on Bonding & Crystal Structures by Zhe Cheng
EGN 3365 Review on Bonding & Crystal Structures 2017 by Zhe Cheng Expectations on Chapter 1 Chapter 1 Understand materials can be classified in different ways by composition, property, application, or
More informationComprehensive Nuclear Materials Modelling Section (R. Stoller, Editor) Chapter 128 Molecular Dynamics.
draft 5/10/10 revised SY 12/5/10 Comprehensive Nuclear Materials Modelling Section (R. Stoller, Editor) Chapter 128 Molecular Dynamics Wei Cai 1, Ju Li 2, Sidney Yip 3 1 Department of Mechanical Engineering,
More informationCrystals Statics. Structural Properties. Geometry of lattices. Aug 23, 2018
Crystals Statics. Structural Properties. Geometry of lattices Aug 23, 2018 Crystals Why (among all condensed phases - liquids, gases) look at crystals? We can take advantage of the translational symmetry,
More information5 Symmetries and point group in a nut shell
30 Phys520.nb 5 Symmetries and point group in a nut shell 5.1. Basic ideas: 5.1.1. Symmetry operations Symmetry: A system remains invariant under certain operation. These operations are called symmetry
More informationStructures of Solids. Unit Cells - Not(?) Chapter 4 Ionic and Other Inorganic Solids. CHEM 462 Wednesday, September 22 T.
Chapter 4 Ionic and Other Inorganic Solids CHEM 462 Wednesday, September 22 T. Hughbanks Structures of Solids Many dense solids are described in terms of packing of atoms or ions. Although these geometric
More informationPartitioning 3D space for parallel many-particle simulations
Partitioning 3D space for parallel many-particle simulations M. A. Stijnman and R. H. Bisseling Mathematical Institute, Utrecht University, PO Box 80010, 3508 TA Utrecht, The Netherlands G. T. Barkema
More informationChapter 3. The structure of crystalline solids 3.1. Crystal structures
Chapter 3. The structure of crystalline solids 3.1. Crystal structures 3.1.1. Fundamental concepts 3.1.2. Unit cells 3.1.3. Metallic crystal structures 3.1.4. Ceramic crystal structures 3.1.5. Silicate
More informationLecture 2. Unit Cells and Miller Indexes. Reading: (Cont d) Anderson 2 1.8,
Lecture 2 Unit Cells and Miller Indexes Reading: (Cont d) Anderson 2 1.8, 2.1-2.7 Unit Cell Concept The crystal lattice consists of a periodic array of atoms. Unit Cell Concept A building block that can
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationMSE 102, Fall 2014 Midterm #2. Write your name here [10 points]:
MSE 102, Fall 2014 Midterm #2 Write your name here [10 points]: Instructions: Answer all questions to the best of your abilities. Be sure to write legibly and state your answers clearly. The point values
More informationActivity 5&6: Metals and Hexagonal Close-Packing
Chemistry 150 Name(s): Activity 5&6: Metals and Hexagonal Close-Packing Metals are chemicals characterized by high thermal and electrical conductivity, malleability and ductility. Atoms are the smallest
More information2. As gas P increases and/or T is lowered, intermolecular forces become significant, and deviations from ideal gas laws occur (van der Waal equation).
A. Introduction. (Section 11.1) CHAPTER 11: STATES OF MATTER, LIQUIDS AND SOLIDS 1. Gases are easily treated mathematically because molecules behave independently. 2. As gas P increases and/or T is lowered,
More informationCondensed Matter Physics April, 8, LUMS School of Science and Engineering
Condensed Matter Physics April, 8, 0 LUMS School of Science and Engineering PH-33 Solution of assignment 5 April, 8, 0 Interplanar separation Answer: To prove that the reciprocal lattice vector G = h b
More informationGeometry of Crystal Lattice
0 Geometry of Crystal Lattice 0.1 Translational Symmetry The crystalline state of substances is different from other states (gaseous, liquid, amorphous) in that the atoms are in an ordered and symmetrical
More informationIntroduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou
Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties
More informationPhysics 211B : Problem Set #0
Physics 211B : Problem Set #0 These problems provide a cross section of the sort of exercises I would have assigned had I taught 211A. Please take a look at all the problems, and turn in problems 1, 4,
More informationElastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals
PHYSICAL REVIEW B 70, 214113 (2004) Elastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals Sang Kyu Kwak and David A. Kofke Department of Chemical and Biological
More informationMP464: Solid State Physics Problem Sheet
MP464: Solid State Physics Problem Sheet 1) Write down primitive lattice vectors for the -dimensional rectangular lattice, with sides a and b in the x and y-directions respectively, and a face-centred
More informationME 254 MATERIALS ENGINEERING 1 st Semester 1430/ st Med-Term Exam (1.5 hrs)
ME 254 MATERIALS ENGINEERING 1 st Semester 1430/1431 1 st Med-Term Exam (1.5 hrs) قسم الهندسة الميكانيكية Question 1 a) Classify the materials based on their properties and performance, give some examples.
More informationThe Periodic Table and Chemical Reactivity
The and Chemical Reactivity Noble gases Less electronegative elements More electronegative elements Then what is electronegativity? The tendency of an atom to attract an electron (or electron density)
More informationSection 10 Metals: Electron Dynamics and Fermi Surfaces
Electron dynamics Section 10 Metals: Electron Dynamics and Fermi Surfaces The next important subject we address is electron dynamics in metals. Our consideration will be based on a semiclassical model.
More information3.3 Show for the body-centered cubic crystal structure that the unit cell edge length a and the atomic radius R are related through a =4R/ 3.
EGR 1 Materials Science, HW, CH :, 6, 7, 11, 0,, 0a-c,. Show for the body-centered cubic crystal structure that the unit cell edge length a and the atomic radius R are related through a =4R/. Consider
More informationMolecular Dynamics Simulation of Nanometric Machining Under Realistic Cutting Conditions Using LAMMPS
Molecular Dynamics Simulation of Nanometric Machining Under Realistic Cutting Conditions Using LAMMPS Rapeepan Promyoo Thesis Presentation Advisor: Dr. Hazim El-Mounayri Department of Mechanical Engineering
More informationRam Seshadri MRL 2031, x6129, These notes complement chapter 6 of Anderson, Leaver, Leevers and Rawlings
Crystals, packings etc. Ram Seshadri MRL 2031, x6129, seshadri@mrl.ucsb.edu These notes complement chapter 6 of Anderson, Leaver, Leevers and Rawlings The unit cell and its propagation Materials usually
More informationContents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21
Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic
More informationSolid State Physics 460- Lecture 5 Diffraction and the Reciprocal Lattice Continued (Kittel Ch. 2)
Solid State Physics 460- Lecture 5 Diffraction and the Reciprocal Lattice Continued (Kittel Ch. 2) Ewald Construction 2θ k out k in G Physics 460 F 2006 Lect 5 1 Recall from previous lectures Definition
More informationHW# 5 CHEM 281 Louisiana Tech University, POGIL(Process Oriented Guided Inquiry Learning) Exercise on Chapter 3. Structures of Ionic Solids. Why?
HW# 5 CHEM 281 Louisiana Tech University, POGIL(Process Oriented Guided Inquiry Learning) Exercise on Chapter 3. Structures of Ionic Solids. Why? Many ionic structures may be described as close-packed
More informationLab 1: Handout GULP: an Empirical energy code
Lab 1: Handout GULP: an Empirical energy code We will be using the GULP code as our energy code. GULP is a program for performing a variety of types of simulations on 3D periodic solids, gas phase clusters,
More informationMP464: Solid State Physics Problem Sheet
MP464: Solid State Physics Problem Sheet 1 Write down primitive lattice vectors for the -dimensional rectangular lattice, with sides a and b in the x and y-directions respectively, and a face-centred rectangular
More informationBravais Lattice + Basis = Crystal Structure
Bravais Lattice + Basis = Crystal Structure A crystal structure is obtained when identical copies of a basis are located at all of the points of a Bravais lattice. Consider the structure of Cr, a I-cubic
More informationSupplementary Figures
Supplementary Figures 8 6 Energy (ev 4 2 2 4 Γ M K Γ Supplementary Figure : Energy bands of antimonene along a high-symmetry path in the Brillouin zone, including spin-orbit coupling effects. Empty circles
More informationThe Generalized Interpolation Material Point Method
Compaction of a foam microstructure The Generalized Interpolation Material Point Method Tungsten Particle Impacting sandstone The Material Point Method (MPM) 1. Lagrangian material points carry all state
More informationLab 4: Handout Molecular dynamics.
1 Lab 4: Handout Molecular dynamics. In this lab, we will be using MOLDY as our molecular dynamics (MD) code. The MOLDY manual is online. It is a well-written manual, and is an excellent resource for learning
More information1.4 Crystal structure
1.4 Crystal structure (a) crystalline vs. (b) amorphous configurations short and long range order only short range order Abbildungen: S. Hunklinger, Festkörperphysik, Oldenbourg Verlag represenatives of
More informationLab 1: Empirical Energy Methods Due: 2/14/18
Lab 1: Empirical Energy Methods Due: 2/14/18 General remarks on scientific scripting Scientific scripting for managing the input and output data is an important component of modern materials computations,
More informationStudy of Phase Transition in Pure Zirconium using Monte Carlo Simulation
Study of Phase Transition in Pure Zirconium using Monte Carlo Simulation Wathid Assawasunthonnet, Abhinav Jain Department of Physics University of Illinois assawas1@illinois.edu Urbana, IL Abstract NPT
More informationECE 474: Principles of Electronic Devices. Prof. Virginia Ayres Electrical & Computer Engineering Michigan State University
ECE 474: Principles of Electronic Devices Prof. Virginia Ayres Electrical & Computer Engineering Michigan State University ayresv@msu.edu Lecture 06: Completed: Chapter 01: quantify physical structures
More informationQuantum Condensed Matter Physics Lecture 4
Quantum Condensed Matter Physics Lecture 4 David Ritchie QCMP Lent/Easter 2019 http://www.sp.phy.cam.ac.uk/drp2/home 4.1 Quantum Condensed Matter Physics 1. Classical and Semi-classical models for electrons
More informationCrystal Structure and Chemistry
Crystal Structure and Chemistry Controls on Crystal Structure Metallic bonding closest packing Covalent bonding depends on orbital overlap and geometry Ionic bonding Pauling s Rules Coordination Principle
More informationLecture 2 Symmetry in the solid state -
Lecture 2 Symmetry in the solid state - Part II: Crystallographic coordinates and Space Groups. 1 Coordinate systems in crystallography and the mathematical form of the symmetry operators 1.1 Introduction
More information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More informationMachine learning the Born-Oppenheimer potential energy surface: from molecules to materials. Gábor Csányi Engineering Laboratory
Machine learning the Born-Oppenheimer potential energy surface: from molecules to materials Gábor Csányi Engineering Laboratory Interatomic potentials for molecular dynamics Transferability biomolecular
More information2. As gas P increases and/or T is lowered, intermolecular forces become significant, and deviations from ideal gas laws occur (van der Waal equation).
A. Introduction. (Section 11.1) CHAPTER 11: STATES OF MATTER, LIQUIDS AND SOLIDS 1. Gases are easily treated mathematically because molecules behave independently. 2. As gas P increases and/or T is lowered,
More informationSolid State Physics Lecture 3 Diffraction and the Reciprocal Lattice (Kittel Ch. 2)
Solid State Physics 460 - Lecture 3 Diffraction and the Reciprocal Lattice (Kittel Ch. 2) Diffraction (Bragg Scattering) from a powder of crystallites - real example of image at right from http://www.uni-wuerzburg.de/mineralogie/crystal/teaching/pow.html
More informationIonic Bonding. Example: Atomic Radius: Na (r = 0.192nm) Cl (r = 0.099nm) Ionic Radius : Na (r = 0.095nm) Cl (r = 0.181nm)
Ionic Bonding Ion: an atom or molecule that gains or loses electrons (acquires an electrical charge). Atoms form cations (+charge), when they lose electrons, or anions (- charge), when they gain electrons.
More informationSurface Structure and Morphology 2D Crystallography
Surface Structure and Morphology 2D Crystallography Selvage (or selvedge (it. cimosa)): Region in the solid in the vicinity of the mathematical surface Surface = Substrate (3D periodicity) + Selvage (few
More informationLecture Note on Crystal structures Masatsugu Sei Suzuki and Itsuko S. Suzuki Department of Physics, SUNY at Binghamton (Date: February 03, 2012)
Lecture Note on Crystal structures Masatsugu Sei Suzuki and Itsuko S. Suzuki Department of Physics, SUNY at Binghamton (Date: February 03, 2012) This is a part of lecture note on solid state physics (Phys.472/572)
More informationIntermolecular Forces and Monte-Carlo Integration 열역학특수연구
Intermolecular Forces and Monte-Carlo Integration 열역학특수연구 2003.3.28 Source of the lecture note. J.M.Prausnitz and others, Molecular Thermodynamics of Fluid Phase Equiliria Atkins, Physical Chemistry Lecture
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More information3d Crystallization. FCC structures. Florian Theil. Macroscopic Modeling of Materials with Fine Structure May Warwick University
FCC structures 1 Mathematics Institute Warwick University Macroscopic Modeling of Materials with Fine Structure May 2011 Outline 1 2 Outline 1 2 Extremum principles Crystalline structure is the minimizer
More informationInterface Resistance and Thermal Transport in Nano-Confined Liquids
1 Interface Resistance and Thermal Transport in Nano-Confined Liquids Murat Barisik and Ali Beskok CONTENTS 1.1 Introduction...1 1.2 Onset of Continuum Behavior...2 1.3 Boundary Treatment Effects on Interface
More informationA web based crystallographic tool for the construction of nanoparticles
A web based crystallographic tool for the construction of nanoparticles Alexios Chatzigoulas 16/5/2018 + = 1 Outline Introduction Motivation Crystallography theory Creation of a web based crystallographic
More informationChm October Molecular Orbitals. Instructions. Always show your work for full credit.
Inorganic Exam 2 Chm 451 29 October 2009 Name: Instructions. Always show your work for full credit. Molecular Orbitals 1. (4 pts) An unusual structure, and one for which there are probably no real examples
More informationMaterials for Civil and Construction Engineers CHAPTER 2. Nature of Materials
Materials for Civil and Construction Engineers CHAPTER 2 Nature of Materials Bonds 1. Primary Bond: forms when atoms interchange or share electrons in order to fill the outer (valence) shells like noble
More information