Simulations of Slow Condensed Phases
|
|
- Eugene McCarthy
- 5 years ago
- Views:
Transcription
1 Simulations of Slow Condensed Phases Part 1: Metallic Glasses, Alloying & Vibrating Nanoparticles J. Daniel Gezelter Department of Chemistry & Biochemistry University of Notre Dame The role of the experimentalist is to perform crude analog simulations of theoretical predictions. - attributed to Frank Stillinger
2 Overview Two very simple models for diffusion Simulating metallic systems Results from simple glass simulations Nanoparticle dynamics 1. Alloying dynamics of core-shell particles 2. Vibrational and Elastic properties of Au nanoparticles
3 Non-Arrhenius diffusion and non-exponential decay Empirically, self-diffusion in fragile glasses exhibits a crossover from Arrhenius behavior to super-exponential ( e (T 0/T ) 2 ), power law ((T T c ) α ), or Vogel-Fulcher-Tamman ( e B/(T T 0) ) temperature dependence when approaching T g. At these temperatures, time correlation functions (van Hove, cage correlation, etc.) also exhibit stretching with 0 < β < 1: C(t) = e (γt)β Are these observations connected? Will simple models for diffusion capture the phenomena that lead to these behaviors?
4 Zwanzig s interrupted oscillation model Starting with the Green-Kubo formula, D = v(t) v(0) dt Expand v(t) v(0) in terms of the harmonic oscillations of the normal modes of the inherent structures in the liquid: v(t) v(0) = ρ(ω) cos(ωt)e t/τ dω The normal modes are damped by a characteristic lifetime, e t/τ Swap the order of integration, and the resulting expression is: D = k BT m ρ(ω) τ 1 + ω 2 τ 2 dω R. Zwanzig, J. Chem. Phys. 79, 4507 (1983) Stillinger & Weber, Phys. Rev. A 25, 978 (1982) Stillinger & Weber, Phys. Rev. A 28, 2408 (1983) Weber & Stillinger, J. Chem. Phys. 80, 2742 (1984) Stillinger & Weber, J. Chem. Phys. 83, 4767 (1985)
5 Continuous Time Random Walk (CTRW) model Random walks take place on a regular lattice with spacing σ 0 between lattice points, and with a distribution of waiting times, ψ(t) = 1 τ e t/τ γ = 1 ( t τ ) 1 γ 0 < γ < 1 In the γ = 1 limit, this model returns the expected diffusive behavior, D = σ2 0 6τ When γ < 1, however, it is also possible to derive a long-time approximate solution for the mean square displacement, r 2 (t) = 4σ2 0 3 π ( t τ ) γ ( 2 γ γ ) γ 1 (2 γ) 2 Γ(7 2 γ) The CTRW model predicts sub-diffusive transport when the waiting time distribution is pathological. Blumen, Klafter & Zumofen, Phys. Rev. B 27, 3429 (1983) Klafter & Zumofen, J. Phys. Chem. 98, 7366 (1994) Shlesinger, Klafter & Zumofen, Am. J. Phys. 67, 1253 (1999)
6 Common Features & Questions Both transport models require what looks like the same time constant (τ) which measures the characteristic time between diffusive hops. Which model is a better approximation to the transport dynamics in a super-cooled liquid? Does sub-diffusive behavior even make sense on a longer time scale? Figure from: Walter Kob, J. Phys. Condens. Matter 11, R85 (1999) r 2 (t) (a) T=5.0 T=0.466 A particles time t
7 Why metallic systems? They are good experimentally-known glass-formers (e.g. Ag 6 Cu 4 ). Alloys of different sized atoms have built-in lattice frustration. They form glassy systems relatively quickly. There are no orientational degrees of freedom. The potentials are simple interesting.
8 Simulating Transition Metals Spherically symmetric electron density + Positively charged core Potential for Embedding an atom at this location V i = F [ ρ i ] + j i φ ij (r ij ) Embedding Functional Pairwise terms (mainly repulsive) There are three empirical functions that are required for this approach: Electronic Density Function: ρ i = j i ρ j (r ij ) Effective Charge Distribution Function: Density or Embedding Functional: φ ij (r ij ) = Z i(r ij )Z j (r ij ) r ij F [ρ]
9 Other Embedded Atom Methods Finnis-Sinclair: Empirical (but simple) functions for densities, pair potentials and functionals Parametrized for: Cu, Ag, Ni half as fast as LJ M.W. Finnis and J.E. Sinclair, Philos. Mag. A. 50, (1984). EAM: ( ) Densities from Roothan-Hartree-Fock ρ i ρi F i [ ρ i ] = A i E 0 Rose s Universal Embedding Functional ρ 0 ln i ρ 0 i Parameterized for: Cu, Ag, Au, Ni, Pd, Pt and their alloys half as fast as LJ S.M. Foiles, M.I. Baskes, and M.S. Daw, Phys. Rev. B 33, (1986). MEAM: Densities include non-spherical symmetry (s, p, d, and f terms) Analytic Density Functional Screening Function one tenth the speed of LJ Parameterized for: Most transition metals + C, Si, Ge, H, N, O, and compounds M.I. Baskes, Phys. Rev. B 46, (1992). All three methods share the same basic potential structure: V = F [ ρ i ] + i i j>i φ ij (r ij )
10 Results: Structural Features K Below 540K, the second peak in g(r) is split K g(r) 600 K 550 K K K 300 K 0 10 R(Å) 0.3 This coincides with the temperature predicted by the (archaic) Wendt- Abraham measure for Tg. (i.e. when the system crosses an ad hoc value of for the ratio between the first minimum and first peak height in g(r)...) R WA The liquid still has appreciable transport at that temperature, however. T g WA = 547K Temperature (K)
11 Dynamics!% This temperature also appears to mark the onset of fragile behavior in the alloy ::.;9<72!5!4!3!"!!!#!"!!$#!"!!%# $&'()*(+,-.+(/012 ",# ", (768953:4 D )"++ ;<94=3>(?(0:<67(106E57FA()! " )"++ G'HI(?(G9>6!E: :4.()/! (01!& + The CTRW model has the correct temperature dependence of the selfdiffusion constant. Two separate versions of the Zwanzig model show agreement at very low temperatures, but incorrect high-temperature limiting behavior: ",! ",& "!"" #"" $"" %"" &""" &!"" &#"" '()*+
12 Deeper into the CTRW model If we want exact behavior from the CTRW model, we can fix τ and vary σ 0 to get our first unphysical result: the average jump distance increases with decreasing temperature! 4 3! 0 (Å) T (K)
13 Can we extract meaningful waiting time information? Fix σ 0 = Å and γ = 1, use <r 2 (t)> to obtain hopping times: Relax linearity in <r 2 (t)> fits to yield γ for CTRW Cage-Correlation Function CTRW fits to <r 2 (t)> ! hop (fs) ! from Cage Correlation function " from CTRW expression for <r 2 (t)> T (K) T (K)
14 Summary Even with the pathological distribution of hopping times, the CTRW model appears to do a better job than Zwanzig s model at capturing the temperature dependence of the diffusion constant. Structural measures like the Wendt-Abraham estimate for Tg may tell us more about where to look for the mode coupling temperature. They don t appear to be useful at finding Tg. Open Questions: Is there a theory for the hopping distance, σ 0? Does the waiting time distribution require a long-time pathology to result in stretched-exponential correlation functions? If we observe anomalous dynamics at intermediate times, will these become linear at longer times? Are other (more reasonable) waiting time distributions amenable to analytical treatment?
15 Before we started calling them nanoparticles...
16 Preparation of Core-Shell nanoparticles T. Shibata, B.A. Bunker, Z. Zhang, D. Meisel, C.F. Vardeman II, and JDG, J. Am. Chem. Soc. 124, (2002). Spherical Metallic nanoparticles were prepared in sizes ranging from nm in diameter via the reduction of metal salts (Au(CN) 2 and AgClO 4 ) using an organic reducing agent (citrate), or by generating the reducing organic radical radiolytically by γ-irradiation A. Henglein and D. Meisel, Langmuir 14, 7392 (1998).
17 Alloying Dynamics in Core-Shell Nanoparticles Room-temperature alloying with a time scale of approximately 3 days was observed using X-ray absorption fine structure (XAFS) data for the Au(core)-Ag(shell) 12.5 Å-19.8Å particles. Larger particles alloyed on a substantially slower time scale. x Ag Au (A) (B) (D) (C) (E) (F) (H) (G) 2.5nm 3.5nm 4.6nm 8.3nm 20 nm Total Nanoparticle Size (nm) (I) (J) (K) (L) (M)
18 What s the mechanism? Solid-state diffusion constants are 7-9 orders of magnitude too small to explain alloying on this timescale. (A) (B) (C) Diffusion of atoms in a liquid droplet should alloy the metals in a matter of microseconds! (E) (F) (G) What can simulation tell us about the mechanism for alloying? (H) 5nm (J)
19 Initial Nanoparticle Geometries r core 2Å Interface r shell
20 Analysis of Mixing Times In the infinite time limit, the mean square displacement of particles confined to a spherical volume of radius R asymptotically approaches a value of 6R 2 /5. Similarly, the mean square displacement of the radial coordinate (relative to the center of the spherical volume) asymptotically approaches 3R 2 /40. The solutions to the diffusion equation for particles with reflecting boundary conditions at R are: sin(nπr/r) ρ(r, t) = a n r n=1 Many diffusional modes contribute to the particle dynamics. e n2 π 2 Dt/R The KWW function (which is often used to approximate sums of exponentials) and the infinite time limit are helpful in predicting a mixing time: r(t) r(0) 2 (1 ) 6R2 e ( t/τ)β 5 r radial (t) r radial (0) 2 (1 ) 3R2 e ( t/τ)β 40 In these equations, β is not the result of anomalous dynamics. It simply reflects the contributions of many diffusional modes to the overall shape of the mean square displacement. r(t) θ(0,t) r(0)
21 Projected Dynamics at 300K /.!6*7189:3;<7;0*=;<;1<794.6*7189:3;<7;0*=;<;1<794 $!6*7189:3;<7;0*=;<;1<794 /! 01+!*2345- %. Projected Alloying Times at 300K: 0% vacancies: ~ years 5% vacancies: 2.9 days 10% vacancies: 1.5 hours %!!"!!!#!"!!$!"!!$%!"!!$&!"!!$'!"!!$#!"!!% $()*+,!$ -
22 Structural Features!(r)/! 800K, structures sampled from last 1 ns of a 24 ns trajectory 3 0% interfacial vacancies 2 Ag Au % interfacial vacancies!(r)/! r (Å)
23 Summary of Interdiffusional Dynamics Vacancy migration seems to be required to explain the fast alloying dynamics. We expect that a few percent of interfacial vacancies can lead to the dramatic speedup observed via XAFS. Dissimilar metals in core-shell particles may be even more prone to lowtemperature interdiffusion due to lattice mismatch.
24 Elastic Properties of Au Nanoparticles P. Conforti, M. Sprague, C. F. Vardeman II and JDG, in press in J. Phys. Chem. A (2003). Hartland and Hu have shown that the time scale for laser heating (determined by the e-ph coupling constant) is faster than a single period of the breathing mode for spherical nanoparticles. Hot-electron pressure and direct lattice heating contribute to the thermal excitation of the atomic degrees of freedom, and both mechanisms are rapid enough to coherently excite the breating mode of the spherical particles. Intensity (a.u.) 0 Au 15 nm Diameter!"!"!#$!# e-ph coupling T e = T l ph-ph coupling Delay Time (ps) Is the dephasing seen following coherent excitation due to inhomogeneous size distributions in the sample? Or is it due to softening of the breathing mode vibrational frequency following a melting transition? Period (ps) Is the flattening in the period vs. laser intensity due to melting? Anharmonicity? Something else? Intensity (µj) 8 10
25 Simulating Laser Heating We have looked at the initial (40 ps) response to instantaneous heating in a range of spherical nanoparticles: 20 Å (1926 atoms), 25 Å (3884 atoms), 30 Å (6602 atoms), 35 Å (10606 atoms) hν The initial heating is easy: simply resample the velocities from a Maxwell- Boltzmann distribution at twice the target temperature.
26 Extracting the Dynamics of the Breathing Mode How does one measure the volume of a collection of atoms (without a periodic box)? We use the volume of the convex hull as an approximate measure of particle volume Convex hulls overestimate the volume (and can fail when the particles split). But they are better estimates than the minimum enclosing sphere!
27 Breathing Dynamics 2.0e Å particle Hulll Volume (Å 3 ) 1.9e e e e e e e Å particle time (ps)
28 Beat frequencies in the volume fluctuations 3 6")7)+89%&:;( ρ V (ω) = V (t) V (0) e iωt dt / #6!! 0!1*%-!1*!-2. $! /. $! $!)7)+89%&:;( Beat pattern is striking at 800K #$!! ##!! #!!! 5!! 3!! 4!! #6!! #$!! ##!! #!!! 5!! 3!! 4!!!$! " #! #" $! %&'()*+,-
29 Dependence on Temperature and Radius Å Radius 25 Å Radius 30 Å Radius 35 Å Radius Period (ps) Period (ps) Intensity (µj) Experiments show that increasing the laser intensity increases the period of the oscillation up to a threshold T (K) Simulations Show: Period is increasing with increasing temperature even beyond the melting transition. Beat Frequencies.
30 Three approaches to calculating the Bulk Modulus The Thermodynamic Approach: K = 1 κ = V ( P V ) The Linear Response Approach: The Extended System Approach: lim δρ( k) 2 = k BT ρ k 0 V K 2 δρ( k) = V e i k r (ρ( r, t) ρ ) dv H NV T = U + fk B T ext ( τ 2 T χ(t)2 2 + t 0 χ(s)ds ) This is, to within a constant, the Helmholtz free energy, A. Since P K = V K V P = T ( S V ( ( T ( = ( A V ) T ) T T ( U V ) T ) ( 2 S V 2 U 2 V T )T 2 2 U V 2 )T ) K = V ( 2 H NV T V 2 )T H NVT is nearly (but not perfectly) conserved in NVT simulations.
31 The bulk modulus and the surface-melting transition (## '"#,*+789- '## "# 0 1 *+2*345!' *,!' - # $#!# "# /#.#."*6.#*6 ("*6 (#*6 (#!"# $"# %"# &"# '#"# ''"# '("# )*+,-
32 The Lamb Model for Elastic Spheres H. Lamb. Proc. London Math. Soc. 13, (1882). G.V. Hartland, M. Hu and J.E. Sader, J. Phys. Chem. (2003). τ t = 2πR θc t where tan θ = 3θ 3 θ 2 τ l = 2πR ηc l where tan η = 4η 4 η 2 c 2 l /c2 t In a liquid with surface stress, γ, the transcendental equation that must be solved is: τ = 2R ηc l (l) where η cot η = 1 η2 c 2 l 4c 2 l 2γ ρr In the 35 Å particle, the Lamb theory (and liquid extension) give a reasonable explanation for the splittings at low and high temperatures: Period of Toroidal Mode (at 300K) Period of Longitudinal (breathing) mode (at 300K) Period of Liquid droplet breathing mode (at Tm) 2.46 ps 2.56 ps 2.73 ps
33 Conclusions The initial drop-off in the bulk modulus can be attributed to surface melting The peak in the bulk modulus can be explained by lattice reformation upon compression. The volume fluctuation power spectra are picking up both longitudinal and toroidal modes at low temperatures. At higher temperatures, we see the melting take place during the duration of the simulation as evidenced by the two frequency components present in the volume fluctuation power spectra. Since core-shell melting will only drop the frequency of the volume fluctuations (and won t give two frequency components), we can safely say that we see both the crystalline and liquid droplet forms of the nanoparticles in the first 40 ps.
34 Acknowledgments Support From: DOE / NDRL NSF CAREER (CHE ) NSF MRI (DMR ) Dreyfus New Faculty Award Alfred P. Sloan Foundation Megan Sprague Computational Geometry & Vacancy Migration Charles F. Vardeman II All of the Above Patrick Conforti Elastic Nanoparticles
T. 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 informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationATOMISTIC MODELING OF DIFFUSION IN ALUMINUM
ATOMISTIC MODELING OF DIFFUSION IN ALUMINUM S. GRABOWSKI, K. KADAU and P. ENTEL Theoretische Physik, Gerhard-Mercator-Universität Duisburg, 47048 Duisburg, Germany (Received...) Abstract We present molecular-dynamics
More informationA Review of Liquid-Glass Transitions
A Review of Liquid-Glass Transitions Anne C. Hanna December 14, 2006 Abstract Supercooling of almost any liquid can induce a transition to an amorphous solid phase. This does not appear to be a phase transition
More informationGlass-Transition and Side-Chain Dynamics in Thin Films: Explaining. Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate)
Supporting Information for Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate) David D. Hsu, Wenjie Xia, Jake
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationDynamics of Supercooled Liquids The Generic Phase Diagram for Glasses
Dynamics of Supercooled Liquids The Generic Phase Diagram for Glasses A normal liquid will crystallize at a melting temperature T m as it is cooled via a first-order phase transition (see figure above).
More informationFundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres
Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Basic Principles Equations of Hydrostatic Equilibrium and Mass Conservation Central Pressure, Virial
More informationTransient lattice dynamics in fs-laser-excited semiconductors probed by ultrafast x-ray diffraction
Transient lattice dynamics in fs-laser-excited semiconductors probed by ultrafast x-ray diffraction K. Sokolowski-Tinten, M. Horn von Hoegen, D. von der Linde Inst. for Laser- and Plasmaphysics, University
More informationFlow of Glasses. Peter Schall University of Amsterdam
Flow of Glasses Peter Schall University of Amsterdam Liquid or Solid? Liquid or Solid? Example: Pitch Solid! 1 day 1 year Menkind 10-2 10 0 10 2 10 4 10 6 10 8 10 10 10 12 10 14 sec Time scale Liquid!
More informationMicroscopic Picture of Aging in SiO 2 : A Computer Simulation
Microscopic Picture of Aging in SiO 2 : A Computer Simulation Katharina Vollmayr-Lee, Robin Bjorkquist, Landon M. Chambers Bucknell University & Göttingen 7 td tb r n (t) 6 R 5 4 3 2 ti t [ns] waiting
More informationPhysics of disordered materials. Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University
Physics of disordered materials Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University Course plan Familiarity with the basic description of disordered structures
More informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
More informationInteraction X-rays - Matter
Interaction X-rays - Matter Pair production hν > M ev Photoelectric absorption hν MATTER hν Transmission X-rays hν' < hν Scattering hν Decay processes hν f Compton Thomson Fluorescence Auger electrons
More informationAcoustic study of nano-crystal embedded PbO P 2 O 5 glass
Bull. Mater. Sci., Vol. 9, No. 4, August 6, pp. 357 363. Indian Academy of Sciences. Acoustic study of nano-crystal embedded PbO P O 5 glass SUDIP K BATABYAL, A PAUL, P ROYCHOUDHURY and C BASU* Department
More informationStructural Signatures of Mobility in Jammed and Glassy Systems
Lisa Manning Sam Schoenholz Ekin Dogus Cubuk Brad Malone Tim Kaxiras Joerg Rottler Rob Riggleman Jennifer Rieser Doug Durian Daniel Sussman Carl Goodrich Sid Nagel Structural Signatures of Mobility in
More informationhydrate systems Gránásy Research Institute for Solid State Physics & Optics H-1525 Budapest, POB 49, Hungary László
Phase field theory of crystal nucleation: Application to the hard-sphere and CO 2 Bjørn Kvamme Department of Physics, University of Bergen Allégaten 55, N-5007 N Bergen, Norway László Gránásy Research
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 informationSub -T g Relaxation in Thin Glass
Sub -T g Relaxation in Thin Glass Prabhat Gupta The Ohio State University ( Go Bucks! ) Kyoto (January 7, 2008) 2008/01/07 PK Gupta(Kyoto) 1 Outline 1. Phenomenology (Review). A. Liquid to Glass Transition
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 8 Oct 1996
December 21, 2013 arxiv:cond-mat/9610066v1 [cond-mat.stat-mech] 8 Oct 1996 Some Finite Size Effects in Simulations of Glass Dynamics Jürgen Horbach, Walter Kob, Kurt Binder Institut für Physik, Johannes
More informationWhat 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 informationChapter 4. The Effect of Elastic Softening and Cooperativity on the Fragility of
Chapter 4 The Effect of Elastic Softening and Cooperativity on the Fragility of Glass-Forming Metallic Liquids Key words: Amorphous metals, Shear transformation zones, Ultrasonic measurement, Compression
More informationEnergetics in Ice VI
Energetics in Ice VI Nishanth Sasankan 2011 NSF / REU PROJECT Physics Department University of Notre Dame Advisor: Dr. Kathie E. Newman Abstract There are many different phases of Ice, which exist in different
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
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 informationA Toy Model. Viscosity
A Toy Model for Sponsoring Viscosity Jean BELLISSARD Westfälische Wilhelms-Universität, Münster Department of Mathematics SFB 878, Münster, Germany Georgia Institute of Technology, Atlanta School of Mathematics
More informationThe glass transition as a spin glass problem
The glass transition as a spin glass problem Mike Moore School of Physics and Astronomy, University of Manchester UBC 2007 Co-Authors: Joonhyun Yeo, Konkuk University Marco Tarzia, Saclay Mike Moore (Manchester)
More informationCase study: molecular dynamics of solvent diffusion in polymers
Course MP3 Lecture 11 29/11/2006 Case study: molecular dynamics of solvent diffusion in polymers A real-life research example to illustrate the use of molecular dynamics Dr James Elliott 11.1 Research
More informationStatistical Mechanics of Jamming
Statistical Mechanics of Jamming Lecture 1: Timescales and Lengthscales, jamming vs thermal critical points Lecture 2: Statistical ensembles: inherent structures and blocked states Lecture 3: Example of
More informationModeling of Micro-Fluidics by a Dissipative Particle Dynamics Method. Justyna Czerwinska
Modeling of Micro-Fluidics by a Dissipative Particle Dynamics Method Justyna Czerwinska Scales and Physical Models years Time hours Engineering Design Limit Process Design minutes Continious Mechanics
More informationRelevance of jamming to the mechanical properties of solids Sidney Nagel University of Chicago Capri; September 12, 2014
Relevance of jamming to the mechanical properties of solids Sidney Nagel University of Chicago Capri; September 1, 014 What is role of (dis)order for mechanical behavior? Andrea J. Liu Carl Goodrich Justin
More informationMonte Carlo simulations of alloy segregation in PtAg octahedral nanoparticles
Monte Carlo simulations of alloy segregation in PtAg octahedral nanoparticles Louis C. Jones 6/8/12 Abstract Simulations were carried out to investigate phase segregation of insoluble alloy nanoparticles
More informationVelocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores
PHYSICAL REVIEW E VOLUME 62, NUMBER 1 JULY 2000 Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores A. Verdaguer and J. A. Padró Departament de Física
More informationSupporting Information
Supporting Information Improved Working Model for Interpreting the Excitation Wavelength- and Fluence-Dependent Response in Pulsed aser-induced Size Reduction of Aqueous Gold Nanoparticles Daniel Werner
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 informationQENS in the Energy Domain: Backscattering and Time-of
QENS in the Energy Domain: Backscattering and Time-of of-flight Alexei Sokolov Department of Polymer Science, The University of Akron Outline Soft Matter and Neutron Spectroscopy Using elastic scattering
More information12. Spectral diffusion
1. Spectral diffusion 1.1. Spectral diffusion, Two-Level Systems Until now, we have supposed that the optical transition frequency of each single molecule is a constant (except when we considered its variation
More informationExploring the anomalous behavior of metal nanocatalysts with finite temperature AIMD and x-ray spectra
Exploring the anomalous behavior of metal nanocatalysts with finite temperature AIMD and x-ray spectra F.D. Vila DOE grant DE-FG02-03ER15476 With computer support from DOE - NERSC. Importance of Theoretical
More informationPlan of the lectures
Plan of the lectures 1. Introductory remarks on metallic nanostructures Relevant quantities and typical physical parameters Applications. Linear electron response: Mie theory and generalizations 3. Nonlinear
More informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
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 informationSupporting Information for Effect of. Polymer/Solid and Polymer/Vapor. Instantaneous Interfaces on the Interfacial
Supporting Information for Effect of Polymer/Solid and Polymer/Vapor Instantaneous Interfaces on the Interfacial Structure and Dynamics of Polymer Melt Systems Selemon Bekele and Mesfin Tsige Department
More informationModeling aluminum using molecular dynamics simulation
Journal of Materials and Environmental Sciences ISSN : 2028-2508 CODEN : JMESCN Copyright 2017, University of Mohammed Premier Oujda Morocco J. Mater. Environ. Sci., 2018, Volume 9, Issue 1, Page 93-99
More informationHeat Transport in Glass-Forming Liquids
Heat Transport in Glass-Forming Liquids by VINAY VAIBHAV The Institute of Mathematical Sciences CIT Campus, Taramani, Chennai 600 113, India. A thesis submitted in partial fulfillment of requirements for
More informationLocalization I: General considerations, one-parameter scaling
PHYS598PTD A.J.Leggett 2013 Lecture 4 Localization I: General considerations 1 Localization I: General considerations, one-parameter scaling Traditionally, two mechanisms for localization of electron states
More informationP. W. Anderson [Science 1995, 267, 1615]
The deepest and most interesting unsolved problem in solid state theory is probably the theory of the nature of glass and the glass transition. This could be the next breakthrough in the coming decade.
More informationAn EAM potential for the dynamical simulation of Ni-Al alloys
J. At. Mol. Sci. doi: 10.4208/jams.022310.031210a Vol. 1, No. 3, pp. 253-261 August 2010 An EAM potential for the dynamical simulation of Ni-Al alloys Jian-Hua Zhang, Shun-Qing Wu, Yu-Hua Wen, and Zi-Zhong
More informationCoarse-grained Models for Oligomer-grafted Silica Nanoparticles
Modeling So+ Ma-er: Linking Mul3ple Length and Time Scales KITP Conference, Santa Barbara, June 4-8, 2012 Coarse-grained Models for Oligomer-grafted Silica Nanoparticles Bingbing Hong Alexandros Chremos
More informationProducing a Sound Wave. Chapter 14. Using a Tuning Fork to Produce a Sound Wave. Using a Tuning Fork, cont.
Producing a Sound Wave Chapter 14 Sound Sound waves are longitudinal waves traveling through a medium A tuning fork can be used as an example of producing a sound wave Using a Tuning Fork to Produce a
More informationMSE 383, Unit 3-3. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.
Dynamic Mechanical Behavior MSE 383, Unit 3-3 Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept. Scope Why DMA & TTS? DMA Dynamic Mechanical Behavior (DMA) Superposition Principles
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 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 informationSpeeding up path integral simulations
Speeding up path integral simulations Thomas Markland and David Manolopoulos Department of Chemistry University of Oxford Funded by the US Office of Naval Research and the UK EPSRC Outline 1. Ring polymer
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
More informationKinetics of the scavenger reaction
J. Phys. A: Math. Gen. 17 (1984) L451-L455. Printed in Great Britain LETTER TO THE EDITOR Kinetics of the scavenger reaction S Redner and K Kang Center for Polymer Studiest and Department of Physics, Boston
More informationGlass Formation and Thermodynamics of Supercooled Monatomic Liquids
pubs.acs.org/jpcb Glass Formation and Thermodynamics of Supercooled Monatomic Liquids Vo Van Hoang* Department of Physics, Institute of Technology, National University of HochiMinh City, 268 Ly Thuong
More informationG. R. Strobl, Chapter 5 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). J. Ferry, "Viscoelastic Behavior of Polymers"
G. R. Strobl, Chapter 5 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). J. Ferry, "Viscoelastic Behavior of Polymers" Chapter 3: Specific Relaxations There are many types of relaxation processes
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 informationIntermediate valence in Yb Intermetallic compounds
Intermediate valence in Yb Intermetallic compounds Jon Lawrence University of California, Irvine This talk concerns rare earth intermediate valence (IV) metals, with a primary focus on certain Yb-based
More informationMolecular dynamics simulations of EXAFS in germanium
Cent. Eur. J. Phys. 93 2011 710-715 DOI: 10.2478/s11534-010-0074-0 Central European Journal of Physics Molecular dynamics simulations of EXAFS in germanium Research Article Janis Timoshenko Alexei Kuzmin
More informationCHAPTER 9 Statistical Physics
CHAPTER 9 Statistical Physics 9.1 9.2 9.3 9.4 9.5 9.6 9.7 Historical Overview Maxwell Velocity Distribution Equipartition Theorem Maxwell Speed Distribution Classical and Quantum Statistics Fermi-Dirac
More informationMuon Spin Relaxation Functions
Muon Spin Relaxation Functions Bob Cywinski Department of Physics and Astronomy University of eeds eeds S 9JT Muon Training Course, February 005 Introduction Positive muon spin relaxation (µsr) is a point-like
More informationSTRONG CONFIGURATIONAL DEPENDENCE OF ELASTIC PROPERTIES OF A CU-ZR BINARY MODEL METALLIC GLASS
Chapter 3 STRONG CONFIGURATIONAL DEPENDENCE OF ELASTIC PROPERTIES OF A CU-ZR BINARY MODEL METALLIC GLASS We report the strong dependence of elastic properties on configurational changes in a Cu-Zr binary
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 informationSupplementary Materials
Supplementary Materials Sample characterization The presence of Si-QDs is established by Transmission Electron Microscopy (TEM), by which the average QD diameter of d QD 2.2 ± 0.5 nm has been determined
More informationModified Maxwell-Garnett equation for the effective transport coefficients in inhomogeneous media
J. Phys. A: Math. Gen. 3 (998) 7227 7234. Printed in the UK PII: S0305-4470(98)93976-2 Modified Maxwell-Garnett equation for the effective transport coefficients in inhomogeneous media Juris Robert Kalnin
More informationThermodynamic aspects of
Thermodynamic aspects of nanomaterials Advanced nanomaterials H.HofmannHofmann EPFL-LTP 2011/2012 ÉCOLE POLYTECHNIQUE FÉDÉRALE DE LAUSANNE Thermodynamic properties p of nanosized materials 100000 120 Total
More informationSupporting Materials
Supporting Materials Figure S1 Experimental Setup Page Figure S (a) (b) (c) Feynman Diagrams Page 3-6 Figure S3 D IR Spectra Page 7 Figure S4 Kinetic Model Page 8 Figure S5 Van t Hoff Plots Page 9 1 k
More informationInteratomic Potentials. The electronic-structure problem
Interatomic Potentials Before we can start a simulation, we need the model! Interaction between atoms and molecules is determined by quantum mechanics: Schrödinger Equation + Born-Oppenheimer approximation
More informationExperimental Physics EP3 Atoms and Molecules Spectroscopy X-rays, lasers
Experiental Physics EP3 Atos and Molecules Spectroscopy X-rays, lasers http://research.uni-leipzig.de/valiu/ Experiental Physics III X-ray and laser spectroscopy 1 Bresstrahlung Tungsten ev E ax 1 ax h
More informationGravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces
Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry Planetary Surfaces Gravity & Rotation Polar flattening caused by rotation is the largest deviation from a sphere for a planet sized object (as
More informationAnharmonicity and Quantum Effect in Thermal Expansion of an Invar Alloy
Anharmonicity and Quantum Effect in Thermal Expansion of an Invar Alloy 1 Institute for Molecular Science (IMS), Okazaki The Graduate University for Advanced Studies (Sokendai) Toshihiko YOKOYAMA 1,, Keitaro
More informationAs in my QM3 and QM4, it is useful to replace the phase by its gradient: In the present, nonlinear case the hydrodynamic equations take a form:
1 Lecture 3 BEC What happens if we subject the condensate to an action of a time dependent potential? In the mean field approximation we may expect that the single particle wave function will satisfy the
More informationHow DLS Works: Interference of Light
Static light scattering vs. Dynamic light scattering Static light scattering measures time-average intensities (mean square fluctuations) molecular weight radius of gyration second virial coefficient Dynamic
More informationLight in Matter (Hecht Ch. 3)
Phys 531 Lecture 3 9 September 2004 Light in Matter (Hecht Ch. 3) Last time, talked about light in vacuum: Maxwell equations wave equation Light = EM wave 1 Today: What happens inside material? typical
More informationX-ray Photon Correlation Spectroscopy (XPCS) at Synchrotron and FEL sources
X-ray Photon Correlation Spectroscopy (XPCS) at Synchrotron and FEL sources Christian Gutt Department of Physics, University ofsiegen, Germany gutt@physik.uni-siegen.de Outline How to measure dynamics
More informationICONE COMPUTER SIMULATION OF DIFFUSION OF PB-BI EUTECTIC IN LIQUID SODIUM BY MOLECULAR DYNAMICS METHOD
Proceedings of ICONE TH International Conference on Nuclear Engineering Arlington, VA, USA, April -, ICONE-3 COMPUTER SIMULATION OF DIFFUSION OF PB-BI EUTECTIC IN LIQUID SODIUM BY MOLECULAR DYNAMICS METHOD
More informationRabi oscillations within TDDFT: the example of the 2 site Hubbard model
Rabi oscillations within TDDFT: the example of the 2 site Hubbard model Johanna I. Fuks, H. Appel, Mehdi,I.V. Tokatly, A. Rubio Donostia- San Sebastian University of Basque Country (UPV) Outline Rabi oscillations
More informationFrequency- and Time-Domain Spectroscopy
Frequency- and Time-Domain Spectroscopy We just showed that you could characterize a system by taking an absorption spectrum. We select a frequency component using a grating or prism, irradiate the sample,
More informationThermal and Mechanical Properties of Pt-Rh Alloys
arxiv:cond-mat/9611241v1 [cond-mat.mtrl-sci] 28 Nov 1996 Thermal and Mechanical Properties of Pt-Rh Alloys G. Dereli, T. Çağın, M. Uludoğan, M. Tomak Department of Physics, Middle East Technical University,
More informationCalculating the hopping rate for diffusion in molecular liquids: CS 2
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 7 15 FEBRUARY 1999 Calculating the hopping rate for diffusion in molecular liquids: CS 2 J. Daniel Gezelter, Eran Rabani, and B. J. Berne a) Department of
More informationKinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys
Kinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys, Scott Dunham Department of Electrical Engineering Multiscale Modeling Hierarchy Configuration energies and transition
More informationSUPPLEMENTARY NOTE 1: ADDITIONAL CHARACTERIZATION OF NANODIAMOND SOLUTIONS AND THE OVERHAUSER EFFECT
1 SUPPLEMENTARY NOTE 1: ADDITIONAL CHARACTERIZATION OF NANODIAMOND SOLUTIONS AND THE OVERHAUSER EFFECT Nanodiamond (ND) solutions were prepared using high power probe sonication and analyzed by dynamic
More informationNonequilibrium transitions in glassy flows. Peter Schall University of Amsterdam
Nonequilibrium transitions in glassy flows Peter Schall University of Amsterdam Liquid or Solid? Liquid or Solid? Example: Pitch Solid! 1 day 1 year Menkind 10-2 10 0 10 2 10 4 10 6 10 8 10 10 10 12 10
More informationStructural characterization. Part 1
Structural characterization Part 1 Experimental methods X-ray diffraction Electron diffraction Neutron diffraction Light diffraction EXAFS-Extended X- ray absorption fine structure XANES-X-ray absorption
More informationarxiv: v1 [cond-mat.soft] 21 Aug 2015
Particle dynamics in fluids with random interactions Lenin S. Shagolsem and Yitzhak Rabin Department of Physics, and Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat Gan 59,
More informationGlass Transition as the Rheological Inverse of Gelation
NNF Summer reading group, July 18 th 2017 Glass Transition as the Rheological Inverse of Gelation ACS Macromolecules 46, 2425-2432 (2013) H Henning Winter Department of Chemical Engineering and Department
More informationCooling dynamics of glass-embedded metal nanoparticles
Cooling dynamics of glass-embedded metal nanoparticles Vincent Juvé, Paolo Maioli, Aurélien Crut, Francesco Banfi,, Damiano Nardi, Claudio Giannetti, Stefano Dal Conte, Natalia Del Fatti and Fabrice Vallée
More informationHow Different is the Dynamics of a Lennard-Jones Binary Fluid from One-Component Lennard-Jones Fluid? 1
How Different is the Dynamics of a Lennard-Jones Binary Fluid from One-Component Lennard-Jones Fluid? Takayuki NARUMI and Michio TOKUYAMA Summary We investigate the dynamics of liquids and supercooled
More informationLecture 06. Fundamentals of Lidar Remote Sensing (4) Physical Processes in Lidar
Lecture 06. Fundamentals of Lidar Remote Sensing (4) Physical Processes in Lidar Light interactions with objects (continued) Resonance fluorescence Laser induced fluorescence Doppler effect (Doppler shift
More informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationKinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationSupporting Information. Influence of Vapor Deposition on Structural. and Charge Transport Properties of. Ethylbenzene Films
Supporting Information Influence of Vapor Deposition on Structural and Charge Transport Properties of Ethylbenzene Films Lucas W. Antony, Nicholas E. Jackson,, Ivan Lyubimov, Venkatram Vishwanath, Mark
More informationMonte Carlo strategies for first-principles simulations of elemental systems
Monte Carlo strategies for first-principles simulations of elemental systems Lev Gelb Department of Materials Science and Engineering, University of Texas at Dallas XSEDE12 Lev Gelb (UT Dallas) Monte Carlo
More informationIntroduction To Materials Science FOR ENGINEERS, Ch. 5. Diffusion. MSE 201 Callister Chapter 5
Diffusion MSE 201 Callister Chapter 5 1 Goals: Diffusion - how do atoms move through solids? Fundamental concepts and language Diffusion mechanisms Vacancy diffusion Interstitial diffusion Impurities Diffusion
More informationOrigins of Mechanical and Rheological Properties of Polymer Nanocomposites. Venkat Ganesan
Department of Chemical Engineering University of Texas@Austin Origins of Mechanical and Rheological Properties of Polymer Nanocomposites Venkat Ganesan $$$: NSF DMR, Welch Foundation Megha Surve, Victor
More informationMolecular Dynamics Study on the Binary Collision of Nanometer-Sized Droplets of Liquid Argon
Molecular Dynamics Study on the Binary Collision Bull. Korean Chem. Soc. 2011, Vol. 32, No. 6 2027 DOI 10.5012/bkcs.2011.32.6.2027 Molecular Dynamics Study on the Binary Collision of Nanometer-Sized Droplets
More informationMolecular Dynamics Simulations of Glass Formation and Crystallization in Binary Liquid Metals
Citation & Copyright (to be inserted by the publisher ) Molecular Dynamics Simulations of Glass Formation and Crystallization in Binary Liquid Metals Hyon-Jee Lee 1,2, Tahir Cagin 2, William A. Goddard
More informationMechanical properties of polymers: an overview. Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore
Mechanical properties of polymers: an overview Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore UGC-NRCM Summer School on Mechanical Property Characterization- June 2012 Overview of polymer
More informationEffective Temperatures in Driven Systems near Jamming
Effective Temperatures in Driven Systems near Jamming Andrea J. Liu Department of Physics & Astronomy University of Pennsylvania Tom Haxton Yair Shokef Tal Danino Ian Ono Corey S. O Hern Douglas Durian
More information