Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
|
|
- Logan Kelley
- 5 years ago
- Views:
Transcription
1 Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao
2 Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of forces Classical nuclei
3 Computer Simulation Multi-Scale Ø Time scale for radiation damage evolution and the corresponding simulation methods C + -Si<100> Defect properties BCA Binary Collision Approximation MD Molecular Dynamics KMC Kinetic Monte Carlo
4 Computer Simulation Multi-Scale Ø Time and length scales are always challenging for materials science
5 Can t Just Simulate Everything at Atomic Level: The Human Simulator? Person length scales Size 1 m mols atoms Person time scales Lifetime 100 years 10 9 s Time scale for atomic motions is s, so time step is s Total number of time steps Person simulation At 10 calculations per each atom per each time step we get 10 x x = calculations Present computing capability is teraflop FLOPS (floating point operations per second) Person simulation takes s years.
6 Limitations of MD Rare Events, System Size Rare Events Rare events occur very infrequently, so do not happen on MD time scales. E.g., diffusion (D cm 2 /s) has atom hops only every 10-3 s. Cannot see them with normal MD simulations (can use temperature accelerated MD, but that is an advanced topic). System size Time for a simulation scales with number of atoms, i.e, O(N). Programs can be parallelized to work by spatial decomposition of different regions on different processors Can treat up to atoms, and increasing, but still limited.
7 Need for Accelerated Molecular Dynamics Methods Many important problems lie beyond the timescales accessible to conventional MD, due to the infrequent surface or bulk diffusion mechanisms involved, including: Vapor-deposited thin film growth (metals and semiconductors) especially heteroepitaxial layers Stress-assisted diffusion of point defects, defect clusters, and dislocations Ion implantation or radiation damage annealing
8 Infrequent Event Systems
9 Transition State Theory (TST)
10 Three Recent Methods for Infrequent Events
11 Hyperdynamics: Basic Approach
12 Parallel Replica Dynamics: Basic Concept
13 Temperature Accelerated Dynamics
14 Temperature Accelerated Dynamics Key issue: As TAD relies upon harmonic TST for validity, any anharmonicity error at T high will lead to a corruption of the dynamics.
15 Example of parallel replica hyperdynamics Ag 10 /Ag(111) diffusion Cu/Cu(100) epitaxial growth A. F. Voter, Phys. Rev. Lett. 78, 3908 (1997). 10 monolayers (ML) per second
16 Example of parallel replica hyperdynamics Competing Kinetics and He Bubble Morphology in W W vacancies W interstitial Additional vacancies q Fast and slow growth regimes are defined relative to typical diffusion hopping times of W interstitials around the He bubble L. Sandoval et al. Phys. Rev. Lett. 114, (2015).
17 TAD Example: 1/4 ML Cu/Cu(111) at 150 K M.R. Sørensen and A.F. Voter, J. Chem. Phys. 112, 9599 (2000).
18 Comparison of Different Acceleration Methods
19 Accelerated Molecular Dynamics q Accelerated Dynamics v Hyperdynamics boost potential ΔV(r) v Parallel replica method for dynamics v Temperature accelerated dynamics v Long-time dynamics based on the dimer method v ABC slowly add penalty functions + static relaxation q Issues " Solve the rough potential surface problem for dynamics simulation " Coupling motion of fast and slow dynamics Energy Landscape - Rough Initial Final
20 Steady-State Accelerated MD q SSAMD method (rough potential surface and slow dynamics) Ø Add boost potentials slowly S 6 S 5 A B S 1 S2 S 1 S 2 S 4 S 1 S 3 A B Ø How to determine boost potential during MD simulation Ø In normal MD, increase T - the total vibration magnitude of atoms, driving the system to change state smoothly from A to B, and overcoming the energy barrier to state B A B
21 Basic Ideas q Vibration magnitude of atoms at a given temperature A i =! R i t! R i 0 λ = n i=1 A i λ (nm) λ of 2000 atoms at 400 K in Fe ª The system reaches an equilibrium state λ becomes a constant
22 Basic Ideas q Vibration magnitude of atoms at a given temperature A i =! R i t! R i 0 λ = n i=1 A i λ (nm) 2000 atoms ª The system occupies the different steady states at different temperatures. ª The states approach a continuum function as the increase in temperature is infinitesimal (temperatureinduced states).
23 Basic Ideas q Vibration magnitude of atoms at a given temperature A i =! R i t! R i 0 λ = n i=1 A i λ (nm) 2000 atoms ª The system occupies the different steady states at different temperatures. ª The states approach a continuum function as the increase in temperature is infinitesimal (temperatureinduced states).
24 Ø Form boost potential ΔV = F(λ S ) Basic Ideas λ S = A i,s 1) A i, S - the vibration magnitude of ith atom in S state 2) Activated volume (only boost atoms in the activated volume) AV
25 Boost Potential Ø Two conditions: 1. both the original and biased systems obey transition state theory ΔV Ab (r) = 0 k TST A = v A δ A (r)e βδv b (r) e βδv b (r) A b A b The bias potential vanishes at any dividing surface ΔV(r) Ab = 0 where δ Ab (r) 0 2. Concave boost potential!!"#$ (!!,,!!!" ;!) =!! (!) 1!!!,..,!!!"!!!!(!(!)!!!,..,!!!" )! E b is constant parameter related to state s and we slowly increase E b (t) and q(t)
26 General Picture Ø System evolves from one equilibrium state to another: λ (Å) ª The boost potential allows the system to jump from S 0 S 1 S 2. ª When the system reaches state S 1, the vibration λ reaches to q 1. ª The system can move along the corresponding constant energy contour. q k = q k 1 + dq E b,k = E b,k 1 + de
27 Transition State Ø The time processing involving the boost potential at step k " t k = t MD exp$ ΔVk # k B T % ' & t MD : normal MD time step Ø Determine initial q!!"#$ (!!,,!!!" ;!) =!! (!) 1!!!,..,!!!"!!!!(!(!)!!!,..,!!!" )! At each state, there exists an average value of total vibration magnitude of n atoms: q s0 -> initial λ at a given temperature q S1 -> state S 1 q- k = q k 1 + dq
28 Validate SSAMD q Migration of a mono vacancy in Fe and Surface Diffusion of Clusters in Cu Ø Ackland s potential for Fe and Cu, and Fe-He potential by PNNL Ø Time step 1 fs, Temperature 200 ~ 500 K. Ø MD simulation in the canonical ensemble of NVT is performed to get the initial equilibrium state and determine the initial q (initial λ at a given temperature). Ø Set appropriate dq and de, and add boost potential to system. At each state, the system is relaxed with boost potential and then evolved to next state (dq ~ A and de ~ 0.3 ev). Ø When a transition occurs, the simulation returns to normal MD which will bring the system to another equilibrium state. Ø Repeat the above steps, and explore the potential surface
29 Vacancy Migration One of the 1 st nearest neighbor jumps to vacancy site Fe atom vacancy site Ø At 300 K, a jump ~ 2.8 x 10-5 s Ø 100 jumps ~ 2.8 x 10-3 s E m ~0.6 ev E m (NEB) ~ 0.63 ev (Ackland JPCM)
30 Surface Diffusion of Cu Atom Ø Surface diffusion atoms go underground Nature Materials 2(2003)211
31 Surface Diffusion of Cu Atom Ø Surface diffusion hop-exchange of a dimer T = 300 K Simulation time ~ 5.6 x10-3 s (5.6 ms)
32 Material Issues in Fission and Fusion Reactors Radiation damage defects He accumulation alters microstructures and properties q Annual defect production in fission and fusion reactor environments Defects Fission (LWR) Fusion (3.5 MW/m 2 ) Displacement per atom/year 5 40 Helium (appm/year) Hydrogen (appm/year) Ø Helium has significant effects in fusion structure materials, affecting the global evolution of microstructures and degrading the mechanical properties of materials. Ø One of the most important issues is helium bubble formation and growth
33 HeV 2 Cluster at 400 K q Two paths Path I Fe atom He atom vacancy Ø Both paths are consistent with ab initio calculations (CC Fu, PRB) Ø Vacancy driven migration of the H n V m clusters (m>n) Path II
34 Migration of He n V m Cluster (n>m) q He 6 V 3 (He 9 V 6 ) (a) (b) v2 v1 v3 v4 v1 Vacancy He atom v2 v4 v4 v1 v3 v1 v2 v2 (c) (d) (1) Creation of SIA (2) Orientation change of the dumbbell interstitial to crowdion (3)The crowdion diffuses along the cluster surface Inters88al v3 (4) Combination with original vacancy: V +SIA->0 Ø SIA driven migration of He n V m (n>m) clusters (a two-step mechanism - creation of an SIA and recombination of the SIA with a vacancy)
35 SSAMD Migration of He-V clusters Total diffusion time is around 0.22s. He 6 V 4 T = 500 K He atom Vacancy Interstitial Mass center of He atoms Original mass center of He atoms Ø An self-interstitial atom driven migration
36 Migration of He n V m Cluster (n>m) q Migration energies of He n V m clusters (200, 300, 400 and 500 K) D = D 0 exp( E a k B T ) Ø D 0 and E a for the He 6 V 4 are e+14 and 0.88 ev, respectively. Ø D 0 and E a for the He 9 V 6 are e+13 and 0.87 ev ª The migration energy is close. ª The mechanisms are same for these two He-V clusters.
37 Growth of A He n V m cluster from Small Clusters q He 6 V 4 and He 9 V 6 He 13 V 10 (600K) He atom Vacancy Interstitial (a) He 9 V 6 He 6 V 4 (b) He 9 V 5 (c) He 10 V 7 He 6 V 5 (d) (e) He 13 V 10 (f) He 5 V µs He 2 The temporal vacancy and interstitial are formed: He 6 V 4 - >He 6 V 5 +SIA Mass transfer occurs between them.
38 SSAMD Growth of He-V clusters Total time is around 6.79µs. Interaction between He 9 V 7 and He 6 V 3 T = 600 K He atom Vacancy Interstitial Ø The small He-V disappears and transfers to the larger one Ø Ostwald Ripening (mass transportation) an important mechanism He 13 V 10 38
39 Electronic Effects: Two Temperature Model 39
40 Electronic Effects: Two Temperature Model 40
41 41
42
43
44
45
46
47
Multiscale modelling of D trapping in W
CMS Multiscale modelling of D trapping in W Kalle Heinola, Tommy Ahlgren and Kai Nordlund Department of Physics and Helsinki Institute of Physics University of Helsinki, Finland Contents Background Plasma-wall
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 informationLocal Hyperdynamics. Arthur F. Voter Theoretical Division Los Alamos National Laboratory Los Alamos, NM, USA
Local Hyperdynamics Arthur F. Voter Theoretical Division National Laboratory, NM, USA Summer School in Monte Carlo Methods for Rare Events Division of Applied Mathematics Brown University Providence, RI
More informationJoint ICTP-IAEA Workshop on Physics of Radiation Effect and its Simulation for Non-Metallic Condensed Matter.
2359-16 Joint ICTP-IAEA Workshop on Physics of Radiation Effect and its Simulation for Non-Metallic Condensed Matter 13-24 August 2012 Introduction to atomistic long time scale methods Roger Smith Loughborough
More informationKinetic Monte Carlo: from transition probabilities to transition rates
Kinetic Monte Carlo: from transition probabilities to transition rates With MD we can only reproduce the dynamics of the system for 100 ns. Slow thermallyactivated processes, such as diffusion, cannot
More informationMolecular Dynamics and Accelerated Molecular Dynamics
Molecular Dynamics and Accelerated Molecular Dynamics Arthur F. Voter Theoretical Division National Laboratory Lecture 3 Tutorial Lecture Series Institute for Pure and Applied Mathematics (IPAM) UCLA September
More informationTime accelerated Atomic Kinetic Monte Carlo for radiation damage modelling
PERFORM 60 FP7 Project Time accelerated Atomic Kinetic Monte Carlo for radiation damage modelling C. Domain, C.S. Becquart, R. Ngayam-Happy EDF R&D Dpt Matériaux & Mécanique des Composants Les Renardieres,
More informationFirst-Passage Kinetic Monte Carlo Algorithm for Complex Reaction-Diffusion Systems
First-Passage Kinetic Monte Carlo Algorithm for Complex Reaction-Diffusion Systems Aleksandar Donev 1 Lawrence Postdoctoral Fellow Lawrence Livermore National Laboratory In collaboration with: Vasily V.
More informationModeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques
Loughborough University Institutional Repository Modeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques This item was submitted to Loughborough University's
More informationInteraction of ion beams with matter
Interaction of ion beams with matter Introduction Nuclear and electronic energy loss Radiation damage process Displacements by nuclear stopping Defects by electronic energy loss Defect-free irradiation
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 informationDiffusion mechanisms in Cu grain boundaries
PHYSICAL REVIEW B VOLUME 62, NUMBER 6 Diffusion mechanisms in Cu grain boundaries 1 AUGUST 2000-II Mads R. So rensen Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 Yuri
More informationComparisons of DFT-MD, TB- MD and classical MD calculations of radiation damage and plasmawallinteractions
CMS Comparisons of DFT-MD, TB- MD and classical MD calculations of radiation damage and plasmawallinteractions Kai Nordlund Department of Physics and Helsinki Institute of Physics University of Helsinki,
More informationJ. Boisse 1,2, A. De Backer 1,3, C. Domain 4,5, C.S. Becquart 1,4
MODELLING SELF TRAPPING AND TRAP MUTATION IN TUNGSTEN USING DFT AND MOLECULAR DYNAMICS WITH AN EMPIRICAL POTENTIAL BASED ON DFT J. Boisse 1,2, A. De Backer 1,3, C. Domain 4,5, C.S. Becquart 1,4 1 Unité
More informationSCIENCE CHINA Physics, Mechanics & Astronomy
SCIENCE CHINA Physics, Mechanics & Astronomy Article April 2012 Vol.55 No.4: 614 618 doi: 10.1007/s11433-012-4679-8 Stability and diffusion properties of self-interstitial atoms in tungsten: a first-principles
More informationMolecular Dynamics Saddle Search Adaptive Kinetic Monte Carlo
Molecular Dynamics Saddle Search Adaptive Kinetic Monte Carlo Samuel T. Chill and Graeme Henkelman Department of Chemistry and the Institute for Computational Engineering and Sciences, The University of
More informationCherry-Pit Structures in Binary Immiscible Alloy Under Ion Irradiation
Cherry-Pit Structures in Binary Immiscible Alloy Under Ion Irradiation Shipeng Shu, Kenneth Tussey May 8, 2011 Abstract We study an special microstructure (matrix atom riched small clusters inside the
More informationKinetic Monte Carlo modelling of semiconductor growth
Kinetic Monte Carlo modelling of semiconductor growth Peter Kratzer Faculty of Physics, University Duisburg-Essen, Germany Time and length scales morphology Ga As 2D islands surface reconstruction Methods
More informationUniform-acceptance force-biased Monte Carlo: A cheap way to boost MD
Dresden Talk, March 2012 Uniform-acceptance force-biased Monte Carlo: A cheap way to boost MD Barend Thijsse Department of Materials Science and Engineering Delft University of Technology, The Netherlands
More informationRadiation damage I. Steve Fitzgerald.
Radiation damage I Steve Fitzgerald http://defects.materials.ox.ac.uk/ Firstly an apology Radiation damage is a vast area of research I cannot hope to cover much in any detail I will try and introduce
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 informationExperience with Moving from Dpa to Changes in Materials Properties
Experience with Moving from Dpa to Changes in Materials Properties Meimei Li, Argonne National Laboratory N. V. Mokhov, Fermilab 46 th ICFA Advanced Beam Dynamics Workshop Sept. 27 Oct. 1, 2010 Morschach,
More informationComputing free energy: Replica exchange
Computing free energy: Replica exchange 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 ls Macroscopic i a t e regime
More informationLecture 10 Thin Film Growth
Lecture 10 Thin Film Growth 1/76 Announcements Homework: Homework Number 2 is returned today, please pick it up from me at the end of the class. Solutions are online. Homework 3 will be set Thursday (2
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 informationTwo simple lattice models of the equilibrium shape and the surface morphology of supported 3D crystallites
Bull. Nov. Comp. Center, Comp. Science, 27 (2008), 63 69 c 2008 NCC Publisher Two simple lattice models of the equilibrium shape and the surface morphology of supported 3D crystallites Michael P. Krasilnikov
More informationDerivation of TTT diagrams with kinetic Monte-Carlo simulations
Derivation of TTT diagrams with kinetic Monte-Carlo simulations A. Gupta, B. Dutta, T. Hickel, J. Neugebauer Department of Computational Materials Design Düsseldorf, Germany 15. September, 2016 BRITS Workshop,
More informationKinetic Monte Carlo (KMC) Kinetic 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 informationModelling of radiation damage in tungsten including He production
Modelling of radiation damage in tungsten including He production C.S. Becquart 1, C. Domain 2 A. De Backer 1 M.F. Barthe 3 M. Hou 4, C. Ortiz 5 1 Unité Matériaux Et Techniques, UMET, UMR 8207, Villeneuve
More informationComputer Simulations and Nanotechnology
Computer Simulations and Nanotechnology Intro: - Making use of high speed modern computers, and therefore nanotechnology - Contributing to the development of nanotechnology Tools: A. Starting from basic
More informationAtomic-Scale Simulations in the Nanoscience of Interfaces
Atomic-Scale Simulations in the Nanoscience of Interfaces Kristen A. Fichthorn Departments of Chemical Engineering and Physics The Pennsylvania State University University Park, PA 16802 USA We Care About
More informationMolecular Dynamics & Adaptive Kinetic Monte Carlo Growth Study of Hydrocarbon Flakes
Molecular Dynamics & Adaptive Kinetic Monte Carlo Growth Study of Hydrocarbon Flakes Amit R. Sharma Department of Physics Wright State University Dayton, Ohio, USA amit.sharma@wright.edu Hydrocarbon co-deposits
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationReconstruction and intermixing in thin Ge layers on Si 001
Reconstruction and intermixing in thin Ge layers on Si 001 L. Nurminen, 1 F. Tavazza, 2 D. P. Landau, 1,2 A. Kuronen, 1 and K. Kaski 1 1 Laboratory of Computational Engineering, Helsinki University of
More informationQUASI-EQUILIBRIUM MONTE-CARLO: OFF-LATTICE KINETIC MONTE CARLO SIMULATION OF HETEROEPITAXY WITHOUT SADDLE POINTS
QUASI-EQUILIBRIUM MONTE-CARLO: OFF-LATTICE KINETIC MONTE CARLO SIMULATION OF HETEROEPITAXY WITHOUT SADDLE POINTS Henry A. Boateng University of Michigan, Ann Arbor Joint work with Tim Schulze and Peter
More informationIntroduction to Accelerated Molecular Dynamics Methods
Introduction to Accelerated Molecular Dynamics Methods Danny Perez and Arthur F. Voter Theoretical Division, T-12 Los Alamos National Laboratory Uncertainty Quantification Workshop April 25-26, 2008 Tucson,
More informationInvestigations of the effects of 7 TeV proton beams on LHC collimator materials and other materials to be used in the LHC
Russian Research Center Kurchatov Institute Investigations of the effects of 7 ev proton beams on LHC collimator materials and other materials to be used in the LHC A.I.Ryazanov Aims of Investigations:
More informationAtomistic Front-End Process Modeling
Atomistic Front-End Process Modeling A Powerful Tool for Deep-Submicron Device Fabrication SISPAD 2001, Athens Martin Jaraiz University of Valladolid, Spain Thanks to: P. Castrillo (U. Valladolid) R. Pinacho
More informationSurface Physics Surface Diffusion. Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science
Surface Physics 008 8. Surface Diffusion Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science Random-Walk Motion Thermal motion of an adatom on an ideal crystal surface: - Thermal excitation the adatom
More informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationKinetic Monte Carlo: Coarsegraining
Kinetic Monte Carlo: Coarsegraining Time and Space Peter Kratzer Faculty of Physics, University Duisburg-Essen, Germany morphology Time and length scales Ga As 2D islands surface reconstruction Methods
More informationNITROGEN CONTAINING ULTRA THIN SiO 2 FILMS ON Si OBTAINED BY ION IMPLANTATION
NITROGEN CONTAINING ULTRA THIN SiO 2 FILMS ON Si OBTAINED BY ION IMPLANTATION Sashka Petrova Alexandrova 1, Evgenia Petrova Valcheva 2, Rumen Georgiev Kobilarov 1 1 Department of Applied Physics, Technical
More informationChristian Ratsch, UCLA
Strain Dependence of Microscopic Parameters and its Effects on Ordering during Epitaxial Growth Christian Ratsch, UCLA Institute for Pure and Applied Mathematics, and Department of Mathematics Collaborators:
More informationStatistical study of defects caused by primary knockon atoms in fcc and bcc metals using molecular. dynamics simulations
Statistical study of defects caused by primary knockon atoms in fcc and bcc metals using molecular dynamics simulations Manoj Warrier, U. Bhardwaj, H. Hemani (CAD, BARC Vizag, India) S. Bukkuru (Nuclear
More informationIntroduction to Semiconductor Physics. Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India
Introduction to Semiconductor Physics 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India http://folk.uio.no/ravi/cmp2013 Review of Semiconductor Physics Semiconductor fundamentals
More informationKinetic Monte Carlo simulation of nucleation on patterned substrates
PHYSICAL REVIEW B, VOLUME 63, 035407 Kinetic Monte Carlo simulation of nucleation on patterned substrates L. Nurminen, A. Kuronen, and K. Kaski Helsinki University of Technology, Laboratory of Computational
More informationMODELING DEFECT MEDIATED DOPANT DIFFUSION IN SILICON. by Brian T. Puchala
MODELING DEFECT MEDIATED DOPANT DIFFUSION IN SILICON by Brian T. Puchala A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Materials Science and
More informationHelium bubbles in bcc Fe and their interactions with irradiation
Loughborough University Institutional Repository Helium bubbles in bcc Fe and their interactions with irradiation This item was submitted to Loughborough University's Institutional Repository by the/an
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 informationAtomistic simulations of point defect diffusion in Si and SiGe
Atomistic simulations of point defect diffusion in Si and SiGe P. Pochet, D. Caliste, K. Rushchanskii, F. Lançon & T. Deutsch CEA UJF INAC Institute for Nanoscience and Cryogenics Partially founded by
More informationFast Monte-Carlo Simulation of Ion Implantation. Binary Collision Approximation Implementation within ATHENA
Fast Monte-Carlo Simulation of Ion Implantation Binary Collision Approximation Implementation within ATHENA Contents Simulation Challenges for Future Technologies Monte-Carlo Concepts and Models Atomic
More informationAndrea Morello. Nuclear spin dynamics in quantum regime of a single-molecule. magnet. UBC Physics & Astronomy
Nuclear spin dynamics in quantum regime of a single-molecule magnet Andrea Morello UBC Physics & Astronomy Kamerlingh Onnes Laboratory Leiden University Nuclear spins in SMMs Intrinsic source of decoherence
More informationHeteroepitaxial growth of Co/ Cu 001 : An accelerated molecular dynamics simulation study
PHYSICAL REVIEW B 72, 035415 2005 Heteroepitaxial growth of Co/ Cu 001 : An accelerated molecular dynamics simulation study Radu A. Miron Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6,
More informationAtomistic simulations on the mobility of di- and tri-interstitials in Si
Atomistic simulations on the mobility of di- and tri-interstitials in Si related publications (since 2001): Posselt, M., Gao, F., Zwicker, D., Atomistic study of the migration of di- and tri-interstitials
More informationWhy Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important
Why Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important The Silicon Age is Ending New materials and growth required Useful in design of Devices
More informationLECTURE 11: Monte Carlo Methods III
1 LECTURE 11: Monte Carlo Methods III December 3, 2012 In this last chapter, we discuss non-equilibrium Monte Carlo methods. We concentrate on lattice systems and discuss ways of simulating phenomena such
More informationMolecular dynamics simulations of the clustering and dislocation loop punching behaviors of noble gas atoms in tungsten
Molecular dynamics simulations of the clustering and dislocation loop punching behaviors of noble gas atoms in tungsten J.Z.Fang, F.Zhou, H.Q.Deng, X.L.Gan, S.F.Xiao, W.Y.Hu Hunan University Contents I,
More informationMultiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models
Multiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models Russel Caflisch 1, Mark Gyure 2, Bo Li 4, Stan Osher 1, Christian Ratsch 1,2, David Shao 1 and Dimitri Vvedensky 3 1 UCLA,
More informationKinetics. Rate of change in response to thermodynamic forces
Kinetics Rate of change in response to thermodynamic forces Deviation from local equilibrium continuous change T heat flow temperature changes µ atom flow composition changes Deviation from global equilibrium
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 informationProgress Report on Chamber Dynamics and Clearing
Progress Report on Chamber Dynamics and Clearing Farrokh Najmabadi, Rene Raffray, Mark S. Tillack, John Pulsifer, Zoran Dragovlovic (UCSD) Ahmed Hassanein (ANL) Laser-IFE Program Workshop May31-June 1,
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary Information Figure S1: (a) Initial configuration of hydroxyl and epoxy groups used in the MD calculations based on the observations of Cai et al. [Ref 27 in the
More informationDopant and Self-Diffusion in Semiconductors: A Tutorial
Dopant and Self-Diffusion in Semiconductors: A Tutorial Eugene Haller and Hughes Silvestri MS&E, UCB and LBNL FLCC Tutorial 1/26/04 1 FLCC Outline Motivation Background Fick s Laws Diffusion Mechanisms
More informationModule 16. Diffusion in solids II. Lecture 16. Diffusion in solids II
Module 16 Diffusion in solids II Lecture 16 Diffusion in solids II 1 NPTEL Phase II : IIT Kharagpur : Prof. R. N. Ghosh, Dept of Metallurgical and Materials Engineering Keywords: Micro mechanisms of diffusion,
More informationFEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH. Yiming Lou and Panagiotis D. Christofides
FEEDBACK CONTROL OF GROWTH RATE AND SURFACE ROUGHNESS IN THIN FILM GROWTH Yiming Lou and Panagiotis D. Christofides Department of Chemical Engineering University of California, Los Angeles IEEE 2003 Conference
More informationDefects and diffusion in metal oxides: Challenges for first-principles modelling
Defects and diffusion in metal oxides: Challenges for first-principles modelling Karsten Albe, FG Materialmodellierung, TU Darmstadt Johan Pohl, Peter Agoston, Paul Erhart, Manuel Diehm FUNDING: ICTP Workshop
More informationPhysics 53. Thermal Physics 1. Statistics are like a bikini. What they reveal is suggestive; what they conceal is vital.
Physics 53 Thermal Physics 1 Statistics are like a bikini. What they reveal is suggestive; what they conceal is vital. Arthur Koestler Overview In the following sections we will treat macroscopic systems
More informationComparison of deuterium retention for ion-irradiated and neutronirradiated
13th International Workshop on Plasma-Facing Materials and Components for Fusion Applications / 1st International Conference on Fusion Energy Materials Science Comparison of deuterium retention for ion-irradiated
More informationDynamical simulations of radiation damage and defect mobility in MgO
PHYSICAL REVIEW B 71, 104102 2005 Dynamical simulations of radiation damage and defect mobility in MgO B. P. Uberuaga, 1 R. Smith, 1, * A. R. Cleave, 2 G. Henkelman, 1 R. W. Grimes, 2 A. F. Voter, 1 and
More informationSimple atomistic modeling of dominant B m I n clusters in boron diffusion
Molecular Simulation, Vol. 31, No. 12, 15 October 2005, 817 824 Simple atomistic modeling of dominant B m I n clusters in boron diffusion J.-H. YOO, C.-O. HWANG, B.-J. KIM and T. WON* Department of Electrical
More informationModelling of the diffusion of self-interstitial atom clusters in Fe-Cr alloys
Modelling of the diffusion of self-interstitial atom clusters in Fe-Cr alloys Dmitry Terentyev, Lorenzo Malerba, Alexander V Barashev To cite this version: Dmitry Terentyev, Lorenzo Malerba, Alexander
More informationOutlook: Application of Positron Annihilation for defects investigations in thin films. Introduction to Positron Annihilation Methods
Application of Positron Annihilation for defects investigations in thin films V. Bondarenko, R. Krause-Rehberg Martin-Luther-University Halle-Wittenberg, Halle, Germany Outlook: Introduction to Positron
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 informationL. Gámez a, B. Gámez a, M.J. Caturla b '*, D. Terentyev c, J.M. Perlado 3 ABSTRACT
Object Kinetic Monte Carlo calculations of irradiated Fe-Cr dilute alloys: The effect of the interaction radius between substitutional Cr and self-interstitial Fe L. Gámez a, B. Gámez a, M.J. Caturla b
More informationThe samples used in these calculations were arranged as perfect diamond crystal of
Chapter 5 Results 5.1 Hydrogen Diffusion 5.1.1 Computational Details The samples used in these calculations were arranged as perfect diamond crystal of a2 2 2 unit cells, i.e. 64 carbon atoms. The effect
More informationFundamental science and synergy of multi-species surface interactions in high-plasma--flux environments
Fundamental science and synergy of multi-species surface interactions in high-plasma--flux environments Our thanks to John Hogan Managed by UT-Battelle for the Department of Energy ReNew PMI, UCLA, March
More informationLecture 1 Modeling and simulation for the growth of thin films
Lecture 1 Modeling and simulation for the growth of thin films Russel Caflisch Mathematics Department Materials Science and Engineering Department UCLA & IPAM www.math.ucla.edu/~material 1 Outline Epitaxial
More informationcollisions of electrons. In semiconductor, in certain temperature ranges the conductivity increases rapidly by increasing temperature
1.9. Temperature Dependence of Semiconductor Conductivity Such dependence is one most important in semiconductor. In metals, Conductivity decreases by increasing temperature due to greater frequency of
More informationStudy of semiconductors with positrons. Outlook:
Study of semiconductors with positrons V. Bondarenko, R. Krause-Rehberg Martin-Luther-University Halle-Wittenberg, Halle, Germany Introduction Positron trapping into defects Methods of positron annihilation
More informationMulti-scale Modeling of Radiation Damage: Large Scale Data Analysis
Journal of Physics: Conference Series PAPER OPEN ACCESS Multi-scale Modeling of Radiation Damage: Large Scale Data Analysis To cite this article: M Warrier et al 2016 J. Phys.: Conf. Ser. 759 012078 View
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 informationIn#uence of microstructure of substrate surface on early stage of thin "lm growth
Vacuum 56 (2000) 185}190 In#uence of microstructure of substrate surface on early stage of thin "lm growth Helin Wei*, Zuli Liu, Kailun Yao Department of Physics, Huazhong University of Science and Technology.
More informationThermally Activated Processes
General Description of Activated Process: 1 Diffusion: 2 Diffusion: Temperature Plays a significant role in diffusion Temperature is not the driving force. Remember: DRIVING FORCE GRADIENT of a FIELD VARIALE
More informationESE 372 / Spring 2013 / Lecture 5 Metal Oxide Semiconductor Field Effect Transistor
Metal Oxide Semiconductor Field Effect Transistor V G V G 1 Metal Oxide Semiconductor Field Effect Transistor We will need to understand how this current flows through Si What is electric current? 2 Back
More informationNucleation theory and the early stages of thin film growth
Nucleation theory and the early stages of thin film growth C. Ratsch a) Department of Mathematics, University of California, Los Angeles, California 90095-1555 J. A. Venables Department of Physics and
More informationPredicting the Dislocation Nucleation Rate as a Function of Temperature and Stress. Abstract
Journal of Materials Research, Invited Feature Article, in press (2011) Predicting the Dislocation Nucleation Rate as a Function of Temperature and Stress Seunghwa Ryu 1, Keonwook Kang 2, and Wei Cai 3
More informationAvailable online at ScienceDirect. Physics Procedia 71 (2015 ) 30 34
Available online at www.sciencedirect.com ScienceDirect Physics Procedia 71 (2015 ) 30 34 18th Conference on Plasma-Surface Interactions, PSI 2015, 5-6 February 2015, Moscow, Russian Federation and the
More informationA dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 15 15 OCTOBER 1999 A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives Graeme Henkelman and Hannes
More informationUpper Critical Dimension for Irreversible Cluster Nucleation and Growth. Abstract
Upper Critical Dimension for Irreversible Cluster Nucleation and Growth Feng Shi, Yunsic Shim, and Jacques G. Amar Department of Physics & Astronomy University of Toledo, Toledo, OH 43606 (Dated: December
More informationMonte Carlo simulation of thin-film growth on a surface with a triangular lattice
Vacuum 52 (1999) 435 440 Monte Carlo simulation of thin-film growth on a surface with a triangular lattice Wei Helin*, Liu Zuli, Yao Kailun Department of Physics, Huazhong University of Science and Technology,
More information31704 Dynamic Monte Carlo modeling of hydrogen isotope. reactive-diffusive transport in porous graphite
31704 Dynamic Monte Carlo modeling of hydrogen isotope reactive-diffusive transport in porous graphite * R. Schneider a, A. Rai a, A. Mutzke a, M. Warrier b,e. Salonen c, K. Nordlund d a Max-Planck-Institut
More informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationFirst passage time Markov chain analysis of rare events for kinetic Monte Carlo: double kink nucleation during dislocation glide
INSTITUTE OF PHYSICS PUBLISHING MODELLING AND SIMULATION IN MATERIALS SCIENCE AND ENGINEERING Modelling Simul. Mater. Sci. Eng. 10 (2002) 581 596 PII: S0965-0393(02)39232-5 First passage time Markov chain
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 informationKinetic Monte Carlo simulation of semiconductor quantum dot growth
Solid State Phenomena Online: 2007-03-15 ISSN: 1662-9779, Vols. 121-123, pp 1073-1076 doi:10.4028/www.scientific.net/ssp.121-123.1073 2007 Trans Tech Publications, Switzerland Kinetic Monte Carlo simulation
More informationComputational Applications in Nuclear Astrophysics using JAVA
Computational Applications in Nuclear Astrophysics using JAVA Lecture: Friday 10:15-11:45 Room NB 6/99 Jim Ritman and Elisabetta Prencipe j.ritman@fz-juelich.de e.prencipe@fz-juelich.de Computer Lab: Friday
More informationPart III. Cellular Automata Simulation of. Monolayer Surfaces
Part III Cellular Automata Simulation of Monolayer Surfaces The concept of progress acts as a protective mechanism to shield us from the terrors of the future. the words of Paul Muad Dib 193 Chapter 6
More informationHigh-Temperature Criticality in Strongly Constrained Quantum Systems
High-Temperature Criticality in Strongly Constrained Quantum Systems Claudio Chamon Collaborators: Claudio Castelnovo - BU Christopher Mudry - PSI, Switzerland Pierre Pujol - ENS Lyon, France PRB 2006
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More informationTheoretical Studies of Self-Diffusion and Dopant Clustering in Semiconductors
phys. stat. sol. (b) 233, No., 24 30 (2002) Theoretical Studies of Self-Diffusion and Dopant Clustering in Semiconductors B. P. Uberuaga )(a), G. Henkelman (b), H. Jónsson (b, c), S. T. Dunham (d), W.
More information