From first principles calculations to statistical physics and simulation of magnetic materials
|
|
- Rosamond Knight
- 5 years ago
- Views:
Transcription
1 From first principles calculations to statistical physics and simulation of magnetic materials Markus Eisenbach Oak Ridge National Laboratory This work was sponsored in parts by the Center for Defect Physics, an Energy Frontier Research Center funded by the US DOE Office of Basic Energy Sciences and by the US DOE Offices of Advanced Scientific Computing Research and Basic Energy Sciences. This research used resources of the Oak Ridge Leadership Computing Facility at ORNL, which is supported by the US DOE, Office of Science.
2 Technological Opportunities in Nano-Magnetism Motivation Ø Magnetic recording/memory: Blocking temperature Super-paramagnetism Ø Magnetic refrigeration: Curie temperature High entropy change Brass bound lodestone Ferrite block NdFeB magnet Ø Permanent magnets: Curie temperature Magnetic moment Large Energy product Issue Ø Understanding and controlling magnetism at the nanoscale: R Harris, A.J. Williams Materials World, Vol. 7 No. 8 pp , August Managed by UT-Battelle for the U.S. Department of Energy
3 Oak Ridge National Laboratory Tennessee Spallation Neutron Source Oak Ridge Leadership Computing Facility OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY 3
4 Opportunities for Modeling Magnetic Nanostructures Spallation Neutron Source Center for Nanophase Materials Sciences Synthesis Characterization Fe 5 nm ~12,000 atoms ~4,000 surface atoms Managed by UT-Battelle for the U.S. Department of Energy
5 Opportunities for Modeling Magnetic Nanostructures Spallation Neutron Source Center for Nanophase Materials Sciences Opportunity Synthesis Fe Characterization Simulation size Real device size Need First Principles approaches capable 5 nm of simulating ~10 4 atoms Methods for dealing with complex magnetic ground states Approaches for dealing ~12,000 with finite atoms temperature statistical physics ~4,000 surface atoms Managed by UT-Battelle for the U.S. Department of Energy
6 Length & Time Scales in Magnetism Time Scale (s) Loss of detailed physics Increasing computational complexity Micromagnetics Spin Models Quantum Length Scale (m) Managed by UT-Battelle for the U.S. Department of Energy
7 Surface Magnetic Nanostructures Co-chains at a Pt(111) Step Edge Ø STM image & magnetic measurements Magnetic moment is canted towards step θ θ Expt =43 o Expt: P. Gambardella et al., Nature 416, 301 (2002) t=1 Managed by UT-Battelle for the U.S. Department of Energy
8 Surface Magnetic Nanostructures Co-chains at a Pt(111) Step Edge Ø STM image & magnetic measurements Magnetic moment is canted towards step θ θ Expt =43 o Expt: P. Gambardella et al., Nature 416, 301 (2002) Energy Scales Ø Exchange and Magneto-crystalline anisotropy treated on same footing E Tot ~10 3 Ry E Coh ~10-1 Ry t=1 E Exch ~10-3 Ry E MAE ~10-6 Ry Managed by UT-Battelle for the U.S. Department of Energy
9 Modeling surface nanostructures Extended impurity Green s function embedding Ø Fully relativistic theory (Dirac equation) Ø Screened KKR layer approach for pristine Pt (111)-surface c - Sites in region to be replaced Ø Green s function embedding method for extended impurity Vacuum Semi-Infinite Solid Managed by UT-Battelle for the U.S. Department of Energy
10 Modeling surface nanostructures Extended impurity Green s function embedding Ø Fully relativistic theory (Dirac equation) Ø Screened KKR layer approach for pristine Pt (111)-surface c - Sites in region to be replaced Ø Green s function embedding method for extended impurity Vacuum Extended Impurity Semi-Infinite Solid Managed by UT-Battelle for the U.S. Department of Energy
11 Modeling surface nanostructures Extended impurity Green s function embedding Ø Fully relativistic theory (Dirac equation) Ø Screened KKR layer approach for pristine Pt (111)-surface c - Sites in region to be replaced Ø Green s function embedding method for extended impurity Vacuum Extended Impurity Semi-Infinite Solid Model geometry Ø 7 Co atoms embedded into a cavity of the Pt(111) surface layer Managed by UT-Battelle for the U.S. Department of Energy
12 Modeling surface nanostructures Extended impurity Green s function embedding Ø Fully relativistic theory (Dirac equation) Ø Screened KKR layer approach for pristine Pt (111)-surface c - Sites in region to be replaced Ø Green s function embedding method for extended impurity Vacuum Extended Impurity Semi-Infinite Solid Model geometry Ø 7 Co atoms embedded into a cavity of the Pt(111) surface layer Managed by UT-Battelle for the U.S. Department of Energy
13 Modeling surface nanostructures Extended impurity Green s function embedding Ø Fully relativistic theory (Dirac equation) Ø Screened KKR layer approach for pristine Pt (111)-surface c - Sites in region to be replaced Ø Green s function embedding method for extended impurity Vacuum Extended Impurity Semi-Infinite Solid Model geometry Ø 7 Co atoms embedded into a cavity of the Pt(111) surface layer Managed by UT-Battelle for the U.S. Department of Energy
14 Modeling surface nanostructures Extended impurity Green s function embedding Ø Fully relativistic theory (Dirac equation) Ø Screened KKR layer approach for pristine Pt (111)-surface c - Sites in region to be replaced Ø Green s function embedding method for extended impurity Vacuum Extended Impurity Semi-Infinite Solid Model geometry Ø 7 Co atoms embedded into a cavity of the Pt(111) surface layer z y θ y φ y x Managed by UT-Battelle for the U.S. Department of Energy
15 Managed by UT-Battelle for the U.S. Department of Energy
16 Managed by UT-Battelle for the U.S. Department of Energy
17 Co-Chain at a Pt(111) Step: SD t=1 Managed by UT-Battelle for the U.S. Department of Energy
18 Co-Chain at a Pt(111) Step: SD t=1 t=10 t=25 t=100 t=800 Managed by UT-Battelle for the U.S. Department of Energy
19 Co-Chain at a Pt(111) Step: SD Full convergence in about 800 time steps t=1 t=10 t=25 t=100 t=800 Managed by UT-Battelle for the U.S. Department of Energy
20 Co-Chain at a Pt(111) Step: SD Full convergence in about 800 time steps t=1 t=10 t=25 t=100 t=800 θ θ Calc =42 o θ Expt =43 o Expt: P. Gambardella et al., Nature 416, 301 (2002) Managed by UT-Battelle for the U.S. Department of Energy
21 LSDA and Finite Temperature Magnetism LSDA is theory of ground states Correct magnetic ground state FM AFM NM Moment (µ B ) Fe Ni Calc. LSDA Expt ~0.5eV
22 LSDA and Finite Temperature Magnetism LSDA is theory of ground states Correct magnetic ground state FM AFM NM Moment (µ B ) Fe Ni Calc. LSDA Expt LSDA magnetism at 7inite temperature ~0.5eV T=0K; FM
23 LSDA and Finite Temperature Magnetism LSDA is theory of ground states Correct magnetic ground state FM AFM NM Moment (µ B ) Fe Ni Calc. LSDA Expt LSDA magnetism at 7inite temperature ~0.5eV T=0K; FM
24 LSDA and Finite Temperature Magnetism LSDA is theory of ground states Correct magnetic ground state FM AFM NM Moment (µ B ) Fe Ni Calc. LSDA Expt LSDA magnetism at 7inite temperature T C (K) Fe Ni Stoner 20,000 10,000 Expt ~0.5eV T=0K; FM
25 LSDA and Finite Temperature Magnetism LSDA is theory of ground states Correct magnetic ground state FM AFM NM Moment (µ B ) Fe Ni Calc. LSDA Expt LSDA magnetism at 7inite temperature T C (K) Fe Ni Stoner 20,000 10,000 Expt PM Disordered Local Moment (DLM) state ~0.5eV T=0K; FM
26 LSDA and Finite Temperature Magnetism LSDA is theory of ground states Correct magnetic ground state FM AFM NM Moment (µ B ) Fe Ni Calc. LSDA Expt LSDA magnetism at 7inite temperature T C (K) Fe Ni Stoner 20,000 10,000 PM Disordered Local Moment (DLM) state Expt T C (K) Fe Ni DLM Expt ~0.5eV T=0K; FM
27 LSDA and Finite Temperature Magnetism LSDA is theory of ground states Correct magnetic ground state FM AFM NM Moment (µ B ) Fe Ni Calc. LSDA Expt ~0.5eV LSDA magnetism at 7inite temperature T C (K) Fe Ni Stoner 20,000 10,000 PM Disordered Local Moment (DLM) state Expt T C (K) Fe Ni DLM Expt Missing Longitudinal Fluctuations T=0K; FM
28 Constrained Local Moments LSDA is theory of ground states Ø General non-collinear state is not a ground state Ø Local moment orientation New degrees of freedom e i i dr m M ( r) Ω = i dr m ( r) Ω i Ø Define general orientational state Set of orientations Selfconsistentcy of the constraint is enforced by ρ δ E i M ( r) = Tr ˆ ρ( r) ( ) = Trσρ ˆ ˆ ( ) M r e e r con " # i $ % & i OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY ρ i e $ i % & ' ( r) ; M( r) ; B ( r) $ con % & ' δ e = 0, i Orientational State e " # i $ % & 11
29 A parallel implementation and scaling of the LSMS method: perfectly scalable at high performance j j m i k l N-2 N-1 N Need only block i of A C B D 1 = (A BD 1 C) 1 m k Calculation dominated by ZGEMM Sustained performance similar to Linpack l
30 13 Wang-Landau Method F. Wang and Landau, PRL 86, 2050 (2001) Conventional Monte Carlo methods calculate expectation values by sampling with a weight given by the Bolzmann distribution In the Wang-Landau Method we rewrite the partition function in terms of the density of states which is calculated by this algorithm Z( ) = { i } To derive an algorithm to estimate g(e) we note that states are randomly generated with a probability proportional to 1/g(E) each energy interval is visited with the same frequency (flat histogram) OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY e H({ i }) = g 0 g(e)e E de
31 Metropolis Method Wang-Landau Method Metropolis et al, JCP 21, 1087 (1953) F. Wang and Landau, PRL 86, 2050 (2001) Z = e E[x]/kBT dx Z Z = g(e)e E/kBT de Compute partition function and other averages with configurations that are weighted with a Boltzmann factor Sample configuration where Boltzmann factor is large. 1. Select configuration E i = E[x i ] 2. Modify configuration (move) E f = E[x f ] 3. Accept move with probability A i f = min{1,e β(e i E f ) } If configurations are accepted with probability 1/g all energies are visited equally (flat histogram) 1. Begin with prior estimate, eg 2. Propose move, accepted with probability 3. Increase DOS at new energy 4. Iterate 2 & 3 until histogram is flat f f = f g 0 (E) =1 A if =min{1,g 0 (E i )/g 0 (E f )} g 0 (E f ) g 0 (E f ) f f>1 5. Reduce and go back to 2
32 Metropolis Method Wang-Landau Method Z Z = e E[x]/k BT dx Z = g(e)e E/k BT de Sample configuration space with probability e E[x]/k BT 1/g(E[x]) Samples mainly regions around energy minima Samples all energies equally -
33 16 Wang-Landau vs Metropolis [same CPU time] Specific heat, Cv / N Metropolis (10 8 trials each run) Wang-Landau Effective Temperature, k B T / ε HH OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY HP Protein model [Ying Wai Li]
34 Not quite embarrassingly parallel Metropolis MC acceptance: A i f =min{1,e β(e i E f ) } random walker 1 random walker 2
35 Not quite embarrassingly parallel Metropolis MC acceptance: A i f =min{1,e β(e i E f ) } Wang-Landau acceptance:a if =min{1,e [g(x f ) g(x i )] } random walker 1... global update of joint DOS at every MC step random walker 2
36 Not quite embarrassingly parallel Metropolis MC acceptance: A i f =min{1,e β(e i E f ) } Wang-Landau acceptance:a if =min{1,e [g(x f ) g(x i )] } random walker 1... global update of joint DOS at every MC step random walker 2 limited by latency ~ microseconds local calculation of energy and observable ~ millisecond to minutes
37 Organization of the WL-LSMS code using a master-slave approach Communicate moment directions and energy Master/driver node controlling WL acceptance, DOS, and histogram LSMS running on N processors to compute energy of particular spinconfigurations
38 Sustained performance of WL-LSMS on Cray XT Petaflop/s on 223,232 cores Petaflop/s on 223,752 cores 829 Teraflop/s on 100,008 cores
39 What it takes to compute a converged DOS on a Cray XT atoms (2x2x2 unit cells) atoms (5x5x5 unit cells) ln g(e) ln g(e) E [Ry] E [Ry] 16 atoms 250 atoms WL walkers total cores 3, ,008 WL samples 23, ,000 CPU-core hours 12,300 4,885,720 Just 2 days
40 Wang-Landau Magnetic Phase Transition: Fe 5x5x5 cell: Curie temperature of Fe from first principles. T C =980K (Exp: 1043K) The internal energy jumps at the Curie temp. leading to a spike in the specific heat. Managed by UT-Battelle for the U.S. Department of Energy
41 22 Cementite (Fe 3 C) Important constituent of carbon steel and cast iron Orthorombic (Pnma) Unit cell: 16 atoms: 4 C Fe Picture: Wikipedia Commons OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY
42 23 Cementite (Fe 3 C) 600 Wang-Landau DOS Specific Heat [Ry/supercell/K] ln g(e) Speci c Heat E [Ry] Temperature [K] Cementite 128 atom unit cell leads to a Tc of 425K OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY
43 24 Comparison with Cementite (Fe3C) exp Curie Temperature calc Curie Temperature Fe 1050K 980K Fe3C 480K 425K Tc(Fe)/ Tc(Fe3C) J. Appl. Phys. 109, 07E138 (2011) OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY
44 Conclusions (I) It is now possible to compute free energies in nanoscale systems - using ab initio methods based on Density Functional Theory - fully taking into account entropy Ab initio calculation of ferromagnetic transition temperature in Fe that does not rely on mean-field approximation - LDA answer based on WL-LSMS underestimate Tc by (only) 7% Future: Improve WL sampling methods - bin-less Wang-Landau - perturbative improvements ( warping ) of model based density of states using first principles calculations. Perturbation calculation of thermodynamic density of states, G. Brown, T. C. Schulthess, D. M. Nicholson, M. Eisenbach, G. M. Stocks, Phys Rev. E 84, (2011). - Convergence improvements of the WL sampling Convergence for the Wang-Landau Density of States, G. Brown, Kh. Odbadrakh, D. M. Nicholson, M. Eisenbach, Phys. Rev. E 84, (R) (2011).
45 26 Effect of lattice vibration on magnetic phase transition in BCC Iron Although BCC Fe is not an ideal Heisenberg magnet [1], Critical exponents Classical Heisenberg Model [2] Experiment [3] α (11) -0.12(1) β (32) 0.363(4) γ 1.385(10) 1.33(2) [1] V. Heine etc, Europhys. Lett. 12 (1990) 545; A.V. Ruban etc, Phys. Rev. B 75.(2007 ) [2] K. Chen etc, Phys. Rev. B 48 (1993) 3249 [3] L.W. Shacklette, Phys. Rev. B 9 (1974) 3789; N. Stusser etc, 31 (1985) 5905 OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY
46 Model MD generated configuration with displacement at T = 900K Sample magnetic configurations using WL-LSMS method [1] Fit J values from J ij = 1 4 (Eij + Eī ȷ Eīj E i ȷ ) [1] M. Eisenbach etc, J. Appl. Phys. 109 (2011) 07E138
47 Model H = U({r n }) X J ij ({r n })ŝ i ŝ j i6=j J(r ) is a cubic function of separation distance. J(r, v, d) is a function of both separation distance and atomic volume v and volume distortion d. d i = X k2nn(i) r k r i r 0 k + r 0 i d is the displacement of nearest neighbors from the ideal crystal position
48 29 Specific Heat OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY
49 30 Magnetization OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY
50 Magnetism at dislocations Dislocations and other defects are important for the understanding of the behavior of structural materials (e.g. Steel, Ni alloys) Magnetism might play a role in the dynamics of dislocations and its effect has not been investigated in detail yet. (Most force fields used in molecular dynamics don t include explicit magnetism) Screw dislocation OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY 31
51 increase in the moment is a coordination effect, not a volumetric effect. We see in Fig. 4 that the atomic volume of atoms near the core is very near the average atomic volume, a.u. In addition to the coordination effect at the core, there is a volume or strain effect between the dislocations. The effect is either compressive or dilational, depending on the handedness of the neighboring dislocations. The moments under compressive strain are reduced and those in dilation are enhanced, as is expected from the typical rela- Screw Dislocation in bcc Fe with LSMS magnetic moments and local volumes 07E159-3 Odbadrakh et al. J. Appl. Phys. 109, 07E159 (2011) bulk tiona from cates lated caus valu core trost addr in re Fully axis held ques from ACK tural the Offic Ridg Labo Depa 1 M. S S. L 3 C. W 4 B. U 5 Y. W 6 S. Is 7 S. Pl 8 M. I 2 3. (Color) local moments are shown for the quadrapole of screw OAKFIG. RIDGE NThe ATIONAL LABORATORY dislocations. U. S. DEPARTMENT OF ENERGY FIG. 5. (Color online) Distribution of local moments from the center of the simulation supercell. FIG. 4. (Color) The Voronii volumes of each site are shown for the quadrapole of screw dislocations. 32
52 Local Magnetization N = 9240 T = 965 K M bulk =
53 Combined MD- SD Simula1ons of BCC Iron
54 Combined MD- SD Simula1ons of BCC Iron
55 Equilibra1on via la:ce thermostat Combined MD- SD Simula1ons of BCC Iron
56 Conclusion (II) l Distribution of local magnetic moment is not uniform but affected by atomic environment parameters. l In BCC Fe, the effect of atomic distance and local volume cancel each other, and an uniform interaction coupling already works well. l Possible applications to dislocation systems.
57 Team Oak Ridge National Laboratory Markus Eisenbach Ying Wai Li University of North Carolina at Asheville Donald M. Nicholson University of Tennessee Khorgolkhuu Odbadrakh Junqi Yin Florida State University Gregory Brown University of Georgia David Landau Dilina Perera
A scalable method for ab initio computation of free energies in nanoscale systems
A scalable method for ab initio computation of free energies in nanoscale systems M. Eisenbach, C.-G. Zhou, D. M. Nicholson, G. Brown, J. Larkin, and T. C. Schulthess Oak Ridge National Laboratory, Oak
More informationParallel Tempering Algorithm in Monte Carlo Simulation
Parallel Tempering Algorithm in Monte Carlo Simulation Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Metropolis Algorithm on Ising Model
More informationEFFECTIVE MAGNETIC HAMILTONIANS: ab initio determination
ICSM212, Istanbul, May 3, 212, Theoretical Magnetism I, 17:2 p. 1 EFFECTIVE MAGNETIC HAMILTONIANS: ab initio determination Václav Drchal Institute of Physics ASCR, Praha, Czech Republic in collaboration
More informationIgor A. Abrikosov Department of Physics, Chemistry, and Biology (IFM), Linköping University, Sweden
Correlation between electronic structure, magnetism and physical properties of Fe-Cr alloys: ab initio modeling Igor A. Abrikosov Department of Physics, Chemistry, and Biology (IFM), Linköping University,
More informationSurface effects in frustrated magnetic materials: phase transition and spin resistivity
Surface effects in frustrated magnetic materials: phase transition and spin resistivity H T Diep (lptm, ucp) in collaboration with Yann Magnin, V. T. Ngo, K. Akabli Plan: I. Introduction II. Surface spin-waves,
More informationOptimized statistical ensembles for slowly equilibrating classical and quantum systems
Optimized statistical ensembles for slowly equilibrating classical and quantum systems IPAM, January 2009 Simon Trebst Microsoft Station Q University of California, Santa Barbara Collaborators: David Huse,
More informationImpact of magnetism upon chemical interactions in Fe alloys
Impact of magnetism upon chemical interactions in Fe alloys A.V. Ruban KTH, Stockholm 2009 M. Ohr, 1985 (Pettifor) C. Wolverton, 2005 Multi-scale engineering design FEM CAD-CAM Engineering design Density
More informationComputational Physics (6810): Session 13
Computational Physics (6810): Session 13 Dick Furnstahl Nuclear Theory Group OSU Physics Department April 14, 2017 6810 Endgame Various recaps and followups Random stuff (like RNGs :) Session 13 stuff
More informationTemperature-dependence of magnetism of free Fe clusters
Temperature-dependence of magnetism of free Fe clusters O. Šipr 1, S. Bornemann 2, J. Minár 2, S. Polesya 2, H. Ebert 2 1 Institute of Physics, Academy of Sciences CR, Prague, Czech Republic 2 Universität
More informationJ ij S i S j B i S i (1)
LECTURE 18 The Ising Model (References: Kerson Huang, Statistical Mechanics, Wiley and Sons (1963) and Colin Thompson, Mathematical Statistical Mechanics, Princeton Univ. Press (1972)). One of the simplest
More informationStructure and Curie temperature of Y 2 Fe 17 x Cr x
Vol. 46 No. 4 SCIENCE IN CHINA (Series G) August 2003 Structure and Curie temperature of Y 2 Fe 17 x Cr x HAO Shiqiang ( ) 1 & CHEN Nanxian ( ) 1,2 1. Department of Physics, Tsinghua University, Beijing
More informationWORLD SCIENTIFIC (2014)
WORLD SCIENTIFIC (2014) LIST OF PROBLEMS Chapter 1: Magnetism of Free Electrons and Atoms 1. Orbital and spin moments of an electron: Using the theory of angular momentum, calculate the orbital
More informationSkyrmions à la carte
This project has received funding from the European Union's Horizon 2020 research and innovation programme FET under grant agreement No 665095 Bertrand Dupé Institute of Theoretical Physics and Astrophysics,
More informationFinite-temperature magnetism of ultrathin lms and nanoclusters PhD Thesis Booklet. Levente Rózsa Supervisor: László Udvardi
Finite-temperature magnetism of ultrathin lms and nanoclusters PhD Thesis Booklet Levente Rózsa Supervisor: László Udvardi BME 2016 Background of the research Magnetic materials continue to play an ever
More informationHanoi 7/11/2018. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.
Hanoi 7/11/2018 Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam. Finite size effects and Reweighting methods 1. Finite size effects 2. Single histogram method 3. Multiple histogram method 4. Wang-Landau
More informationMonte Carlo study of the Baxter-Wu model
Monte Carlo study of the Baxter-Wu model Nir Schreiber and Dr. Joan Adler Monte Carlo study of the Baxter-Wu model p.1/40 Outline Theory of phase transitions, Monte Carlo simulations and finite size scaling
More informationPresentation Groupmeeting June 3 rd, sorry 10 th, 2009 by Jacques Klaasse
Presentation Groupmeeting June 3 rd, sorry 10 th, 2009 by Jacques Klaasse Spin Density Waves This talk is based on a book-chapter on antiferromagnetism, written by Anthony Arrott in Rado-Suhl, Volume IIB,
More informationTheory of magnetic interactions in real materials. Mikhail Katsnelson
Theory of magnetic interactions in real materials Mikhail Katsnelson Outline 1. Introduction 2. Exchange interactions from first principles 3. Beyond DFT: correlated systems and LDA+DMFT 4. Applications:
More informationarxiv: v1 [physics.comp-ph] 6 Mar 2018
Combined Molecular and Spin Dynamics Simulation of BCC Iron with Vacancy Defects Mark Mudrick, 1, Markus Eisenbach, Dilina Perera, 3 and David P. ndau 1 1 Center for Simulational Physics, he University
More informationExtreme scale simulations of high-temperature superconductivity. Thomas C. Schulthess
Extreme scale simulations of high-temperature superconductivity Thomas C. Schulthess T [K] Superconductivity: a state of matter with zero electrical resistivity Heike Kamerlingh Onnes (1853-1926) Discovery
More informationJ 12 J 23 J 34. Driving forces in the nano-magnetism world. Intra-atomic exchange, electron correlation effects: Inter-atomic exchange: MAGNETIC ORDER
Driving forces in the nano-magnetism world Intra-atomic exchange, electron correlation effects: LOCAL (ATOMIC) MAGNETIC MOMENTS m d or f electrons Inter-atomic exchange: MAGNETIC ORDER H exc J S S i j
More informationOutline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary
Outline Introduction: graphene Adsorption on graphene: - Chemisorption - Physisorption Summary 1 Electronic band structure: Electronic properties K Γ M v F = 10 6 ms -1 = c/300 massless Dirac particles!
More informationHidden Interfaces and High-Temperature Magnetism in Intrinsic Topological Insulator - Ferromagnetic Insulator Heterostructures
Hidden Interfaces and High-Temperature Magnetism in Intrinsic Topological Insulator - Ferromagnetic Insulator Heterostructures Valeria Lauter Quantum Condensed Matter Division, Oak Ridge National Laboratory,
More informationVSOP19, Quy Nhon 3-18/08/2013. Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam.
VSOP19, Quy Nhon 3-18/08/2013 Ngo Van Thanh, Institute of Physics, Hanoi, Vietnam. Part III. Finite size effects and Reweighting methods III.1. Finite size effects III.2. Single histogram method III.3.
More informationPhysics 115/242 Monte Carlo simulations in Statistical Physics
Physics 115/242 Monte Carlo simulations in Statistical Physics Peter Young (Dated: May 12, 2007) For additional information on the statistical Physics part of this handout, the first two sections, I strongly
More informationNumerical Analysis of 2-D Ising Model. Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011
Numerical Analysis of 2-D Ising Model By Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011 Contents Abstract Acknowledgment Introduction Computational techniques Numerical Analysis
More informationMagnetism at finite temperature: molecular field, phase transitions
Magnetism at finite temperature: molecular field, phase transitions -The Heisenberg model in molecular field approximation: ferro, antiferromagnetism. Ordering temperature; thermodynamics - Mean field
More informationScaling Theory. Roger Herrigel Advisor: Helmut Katzgraber
Scaling Theory Roger Herrigel Advisor: Helmut Katzgraber 7.4.2007 Outline The scaling hypothesis Critical exponents The scaling hypothesis Derivation of the scaling relations Heuristic explanation Kadanoff
More informationParallel Tempering Algorithm in Monte Carlo Simula5on
Parallel Tempering Algorithm in Monte Carlo Simula5on Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Monte Carlo Algorithms Mo5va5on: Idea:
More informationAn introduction to magnetism in three parts
An introduction to magnetism in three parts Wulf Wulfhekel Physikalisches Institut, Karlsruhe Institute of Technology (KIT) Wolfgang Gaede Str. 1, D-76131 Karlsruhe 0. Overview Chapters of the three lectures
More informationParallel Tempering And Adaptive Spacing. For Monte Carlo Simulation
Parallel Tempering And Adaptive pacing For Monte Carlo imulation Authors: iu Wun Cheung Department of Mathematics, CUHK Zhao, Yiwei Department of Mathematics, CUHK upervisors: Dr. Markus isenbach ORNL
More informationMONTE CARLO METHODS IN SEQUENTIAL AND PARALLEL COMPUTING OF 2D AND 3D ISING MODEL
Journal of Optoelectronics and Advanced Materials Vol. 5, No. 4, December 003, p. 971-976 MONTE CARLO METHODS IN SEQUENTIAL AND PARALLEL COMPUTING OF D AND 3D ISING MODEL M. Diaconu *, R. Puscasu, A. Stancu
More informationshows the difference between observed (black) and calculated patterns (red). Vertical ticks indicate
Intensity (arb. unit) a 5 K No disorder Mn-Pt disorder 5 K Mn-Ga disorder 5 K b 5 K Observed Calculated Difference Bragg positions 24 28 32 2 4 6 8 2 4 2θ (degree) 2θ (degree) Supplementary Figure. Powder
More informationFrom electronic structure to magnetism
From electronic structure to magnetism Olle Eriksson, Department of Physics and Astronomy, Uppsala University And School of Science and Technology, Örebro University Collaborators Abrikosov, Bergman, Bergqvist,
More informationUnderstanding the Magnetic Ground States for Improper Multiferroic Materials
Wright State University CORE Scholar Physics Seminars Physics 4-11-2013 Understanding the Magnetic Ground States for Improper Multiferroic Materials Jason T. Haraldsen Follow this and additional works
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 informationMetropolis Monte Carlo simulation of the Ising Model
Metropolis Monte Carlo simulation of the Ising Model Krishna Shrinivas (CH10B026) Swaroop Ramaswamy (CH10B068) May 10, 2013 Modelling and Simulation of Particulate Processes (CH5012) Introduction The Ising
More informationChapter 3. Magnetic Model. 3.1 Magnetic interactions
Chapter 3 Magnetic Model In this chapter, the micromagnetic model for the description of the magnetic properties of a laterally nanostructured film during growth is presented. The main physical idea of
More informationTransition Elements. pranjoto utomo
Transition Elements pranjoto utomo Definition What is transition metal? One of which forms one or more stable ions which have incompletely filled d orbitals. 30Zn? Definition Zink is not transition elements
More informationAb Initio Calculation of Exchange Interactions, Adiabatic Spin-Waves, and Curie Temperature of Itinerant Ferromagnets
Ab Initio Calculation of Exchange Interactions, Adiabatic Spin-Waves, and Curie Temperature of Itinerant Ferromagnets Patrick Bruno Max-Planck-Institut für Mikrostrukturphysik, Halle, Germany People involved:
More informationarxiv: v1 [cond-mat.mtrl-sci] 15 Jan 2008
New High-T c Half-Heusler Ferromagnets NiMnZ (Z = Si, P, Ge, As) Van An DINH and Kazunori Sato and Hiroshi Katayama-Yoshida The Institute of Scientific and Industrial Research, Osaka University, Mihogaoka
More informationInterface tension of the 3d 4-state Potts model using the Wang-Landau algorithm
Interface tension of the 3d 4-state Potts model using the Wang-Landau algorithm CP3-Origins and the Danish Institute for Advanced Study DIAS, University of Southern Denmark, Campusvej 55, DK-5230 Odense
More informationElectronic, magnetic and spectroscopic properties of free Fe clusters
Electronic, magnetic and spectroscopic properties of free Fe clusters O. Šipr 1, M. Košuth 2, J. Minár 2, S. Polesya 2 and H. Ebert 2 1 Institute of Physics, Academy of Sciences of the Czech Republic,
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationThe Mott Metal-Insulator Transition
Florian Gebhard The Mott Metal-Insulator Transition Models and Methods With 38 Figures Springer 1. Metal Insulator Transitions 1 1.1 Classification of Metals and Insulators 2 1.1.1 Definition of Metal
More informationEvaluation of Wang-Landau Monte Carlo Simulations
2012 4th International Conference on Computer Modeling and Simulation (ICCMS 2012) IPCSIT vol.22 (2012) (2012) IACSIT Press, Singapore Evaluation of Wang-Landau Monte Carlo Simulations Seung-Yeon Kim School
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 25 Mar 2004
Ab-initio spin dynamics applied to nanoparticles: canted magnetism of a finite Co chain along a Pt(111) surface step edge arxiv:cond-mat/43637v1 [cond-mat.mtrl-sci] 25 Mar 24 B. Újfalussy1, B. Lazarovits
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationSpin Superfluidity and Graphene in a Strong Magnetic Field
Spin Superfluidity and Graphene in a Strong Magnetic Field by B. I. Halperin Nano-QT 2016 Kyiv October 11, 2016 Based on work with So Takei (CUNY), Yaroslav Tserkovnyak (UCLA), and Amir Yacoby (Harvard)
More informationMagnetic recording technology
Magnetic recording technology The grain (particle) can be described as a single macrospin μ = Σ i μ i 1 0 1 0 1 W~500nm 1 bit = 300 grains All spins in the grain are ferromagnetically aligned B~50nm Exchange
More informationHigh Temperature High Pressure Properties of Silica From Quantum Monte Carlo
High Temperature High Pressure Properties of Silica From Quantum Monte Carlo K.P. Driver, R.E. Cohen, Z. Wu, B. Militzer, P. Lopez Rios, M. Towler, R. Needs, and J.W. Wilkins Funding: NSF, DOE; Computation:
More informationFirst steps toward the construction of a hyperphase diagram that covers different classes of short polymer chains
Cent. Eur. J. Phys. 12(6) 2014 421-426 DOI: 10.2478/s11534-014-0461-z Central European Journal of Physics First steps toward the construction of a hyperphase diagram that covers different classes of short
More informationKosterlitz-Thouless Transition
Heidelberg University Seminar-Lecture 5 SS 16 Menelaos Zikidis Kosterlitz-Thouless Transition 1 Motivation Back in the 70 s, the concept of a phase transition in condensed matter physics was associated
More informationThe XY-Model. David-Alexander Robinson Sch th January 2012
The XY-Model David-Alexander Robinson Sch. 08332461 17th January 2012 Contents 1 Introduction & Theory 2 1.1 The XY-Model............................... 2 1.2 Markov Chains...............................
More informationLuigi Paolasini
Luigi Paolasini paolasini@esrf.fr LECTURE 5: MAGNETIC STRUCTURES - Mean field theory and magnetic order - Classification of magnetic structures - Collinear and non-collinear magnetic structures. - Magnetic
More information7. Basics of Magnetization Switching
Beyond CMOS computing 7. Basics of Magnetization Switching Dmitri Nikonov Dmitri.e.nikonov@intel.com 1 Outline Energies in a nanomagnet Precession in a magnetic field Anisotropies in a nanomagnet Hysteresis
More informationAvailable online at ScienceDirect. Physics Procedia 57 (2014 ) 82 86
Available online at www.sciencedirect.com ScienceDirect Physics Procedia 57 (2014 ) 82 86 27th Annual CSP Workshops on Recent Developments in Computer Simulation Studies in Condensed Matter Physics, CSP
More informationarxiv: v1 [physics.comp-ph] 21 Jul 2017
A histogram-free multicanonical Monte Carlo algorithm for the basis expansion of density of states Ying Wai Li National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee
More informationSkyrmion à la carte. Bertrand Dupé. Skyrmion à la carte Bertrand Dupé. Institute of Physics, Johannes Gutenberg University of Mainz, Germany
Skyrmion à la carte Bertrand Dupé Institute of Physics, Johannes Gutenberg University of Mainz, Germany 1 Acknowledgement Charles Paillard Markus Hoffmann Stephan von Malottki Stefan Heinze Sebastian Meyer
More informationMonte Carlo Simulation of the Ising Model. Abstract
Monte Carlo Simulation of the Ising Model Saryu Jindal 1 1 Department of Chemical Engineering and Material Sciences, University of California, Davis, CA 95616 (Dated: June 9, 2007) Abstract This paper
More informationl μ M Right hand Screw rule
Magnetic materials Magnetic property The response of the materials to external magnetic field All the materials are magnetic, only the degree of response varies, which is measured in terms of their magnetization
More informationChapter 2. Theoretical background. 2.1 Itinerant ferromagnets and antiferromagnets
Chapter 2 Theoretical background The first part of this chapter gives an overview of the main static magnetic behavior of itinerant ferromagnetic and antiferromagnetic materials. The formation of the magnetic
More informationClassical Monte Carlo Simulations
Classical Monte Carlo Simulations Hyejin Ju April 17, 2012 1 Introduction Why do we need numerics? One of the main goals of condensed matter is to compute expectation values O = 1 Z Tr{O e βĥ} (1) and
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 informationMagnetic ordering, magnetic anisotropy and the mean-field theory
Magnetic ordering, magnetic anisotropy and the mean-field theory Alexandra Kalashnikova kalashnikova@mail.ioffe.ru Ferromagnets Mean-field approximation Curie temperature and critical exponents Magnetic
More informationSupporting Information
Electronic Supplementary Material (ESI) for Journal of Materials Chemistry C. This journal is The Royal Society of Chemistry 2016 Supporting Information Metal-Free Half-Metallicity in a High Energy Phase
More informationThe Magnetic Properties of Superparamagnetic Particles by a Monte Carlo Method
The Magnetic Properties of Superparamagnetic Particles by a Monte Carlo Method D. A. Dimitrov and G. M. Wysin Department of Physics Kansas State University Manhattan, KS 6656-261 (June 19, 1996) We develop
More informationNi 8 Cu n Ni 9. Lectue 4 Trilayers a prototype of multilayers. for FM1 and FM2 interlayer exchange coupling IEC, J inter
Lectue 4 Trilayers a prototype of multilayers Ni 8 Cu n Ni 9 Important parameters: K anisotropy, E band for FM1 and FM2 interlayer exchange coupling IEC, J inter 1 4a Optical and acoustic modes in the
More informationSubir Sachdev. Yale University. C. Buragohain K. Damle M. Vojta
C. Buragohain K. Damle M. Vojta Subir Sachdev Phys. Rev. Lett. 78, 943 (1997). Phys. Rev. B 57, 8307 (1998). Science 286, 2479 (1999). cond-mat/9912020 Quantum Phase Transitions, Cambridge University Press
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
More informationMagnetic Materials. The inductor Φ B = LI (Q = CV) = L I = N Φ. Power = VI = LI. Energy = Power dt = LIdI = 1 LI 2 = 1 NΦ B capacitor CV 2
Magnetic Materials The inductor Φ B = LI (Q = CV) Φ B 1 B = L I E = (CGS) t t c t EdS = 1 ( BdS )= 1 Φ V EMF = N Φ B = L I t t c t B c t I V Φ B magnetic flux density V = L (recall I = C for the capacitor)
More informationNeutron Diffraction Study of Antiferromagnetic Phase Transitions in an Ordered Pt 3 Fe(111) Film
Neutron Diffraction Study of Antiferromagnetic Phase Transitions in an Ordered Pt 3 Fe(111) Film G. J. Mankey, V. V. Krishnamurthy, and I. Zoto MINT Center, The University of Alabama, Tuscaloosa, AL 35487-0209
More informationLectures 16: Phase Transitions
Lectures 16: Phase Transitions Continuous Phase transitions Aims: Mean-field theory: Order parameter. Order-disorder transitions. Examples: β-brass (CuZn), Ferromagnetic transition in zero field. Universality.
More informationAnisotropic Magnetic Structures in Iron-Based Superconductors
Anisotropic Magnetic Structures in Iron-Based Superconductors Chi-Cheng Lee, Weiguo Yin & Wei Ku CM-Theory, CMPMSD, Brookhaven National Lab Department of Physics, SUNY Stony Brook Another example of SC
More informationPhase transitions and finite-size scaling
Phase transitions and finite-size scaling Critical slowing down and cluster methods. Theory of phase transitions/ RNG Finite-size scaling Detailed treatment: Lectures on Phase Transitions and the Renormalization
More informationSupporting information. Realizing Two-Dimensional Magnetic Semiconductors with. Enhanced Curie Temperature by Antiaromatic Ring Based
Supporting information Realizing Two-Dimensional Magnetic Semiconductors with Enhanced Curie Temperature by Antiaromatic Ring Based Organometallic Frameworks Xingxing Li and Jinlong Yang* Department of
More informationFerromagnetism and Metal-Insulator Transition in Hubbard Model with Alloy Disorder
Ferromagnetism and Metal-Insulator Transition in Hubbard Model with Alloy Disorder Krzysztof Byczuk Institute of Physics, Augsburg University Institute of Theoretical Physics, Warsaw University October
More informationOptimization Methods via Simulation
Optimization Methods via Simulation Optimization problems are very important in science, engineering, industry,. Examples: Traveling salesman problem Circuit-board design Car-Parrinello ab initio MD Protein
More informationSupplementary Figures
Supplementary Figures Supplementary Figure 1: Region mapping. a Pristine and b Mn-doped Bi 2 Te 3. Arrows point at characteristic defects present on the pristine surface which have been used as markers
More informationFinite-temperature magnetism in bcc Fe under compression. Xianwei Sha* and R. E. Cohen
Finite-temperature magnetism in bcc Fe under compression Xianwei Sha* and R. E. Cohen Carnegie Institution of Washington, 5251 Broad Branch Road, NW, Washington, D. C. 20015, U. S. A. We investigate the
More informationUniversity of Bristol. 1 Naval Research Laboratory 2 II. Physikalisches Institut, Universität zu Köln
Charge ordering as alternative to Jahn-Teller distortion In collaboration with Michelle Johannes 1, Daniel Khomskii 2 (theory) and Mohsen Abd-Elmeguid et al 2, Radu Coldea et al 3 (experiment) 1 Naval
More informationExchange bias in core/shell magnetic nanoparticles: experimental results and numerical simulations
Exchange bias in core/shell magnetic nanoparticles: experimental results and numerical simulations Xavier Batlle, A. Labarta, Ò. Iglesias, M. García del Muro and M. Kovylina Goup of Magnetic Nanomaterials
More informationBand calculations: Theory and Applications
Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation
More informationQuantum critical itinerant ferromagnetism
Quantum critical itinerant ferromagnetism Belitz et al., PRL 2005 Cavendish Laboratory University of Cambridge Two types of ferromagnetism Localised ferromagnetism: moments localised in real space Ferromagnet
More informationDepartment of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions
Department of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions Friday May 15th, 2014, 14-17h Examiners: Prof. J. Cline, Prof. H. Guo, Prof. G. Gervais (Chair), and Prof. D. Hanna INSTRUCTIONS
More informationApplications to simulations: Monte-Carlo
Applications to simulations: Monte-Carlo A.C. Maggs ESPCI, Paris June 2016 Summary Algorithms Faster/simpler codes Thermodynamics of Electric fields Partition function of electric field Thermal Casimir/Lifshitz
More informationFrustrated diamond lattice antiferromagnets
Frustrated diamond lattice antiferromagnets ason Alicea (Caltech) Doron Bergman (Yale) Leon Balents (UCSB) Emanuel Gull (ETH Zurich) Simon Trebst (Station Q) Bergman et al., Nature Physics 3, 487 (007).
More informationIntroduction to the Renormalization Group
Introduction to the Renormalization Group Gregory Petropoulos University of Colorado Boulder March 4, 2015 1 / 17 Summary Flavor of Statistical Physics Universality / Critical Exponents Ising Model Renormalization
More informationThe Physics of Ferromagnetism
Terunobu Miyazaki Hanmin Jin The Physics of Ferromagnetism Springer Contents Part I Foundation of Magnetism 1 Basis of Magnetism 3 1.1 Basic Magnetic Laws and Magnetic Quantities 3 1.1.1 Basic Laws of
More informationTuning magnetic anisotropy, Kondo screening and Dzyaloshinskii-Moriya interaction in pairs of Fe adatoms
Tuning magnetic anisotropy, Kondo screening and Dzyaloshinskii-Moriya interaction in pairs of Fe adatoms Department of Physics, Hamburg University, Hamburg, Germany SPICE Workshop, Mainz Outline Tune magnetic
More informationMonte Caro simulations
Monte Caro simulations Monte Carlo methods - based on random numbers Stanislav Ulam s terminology - his uncle frequented the Casino in Monte Carlo Random (pseudo random) number generator on the computer
More informationMarkov Chain Monte Carlo Method
Markov Chain Monte Carlo Method Macoto Kikuchi Cybermedia Center, Osaka University 6th July 2017 Thermal Simulations 1 Why temperature 2 Statistical mechanics in a nutshell 3 Temperature in computers 4
More informationPhase Transitions in Condensed Matter Spontaneous Symmetry Breaking and Universality. Hans-Henning Klauss. Institut für Festkörperphysik TU Dresden
Phase Transitions in Condensed Matter Spontaneous Symmetry Breaking and Universality Hans-Henning Klauss Institut für Festkörperphysik TU Dresden 1 References [1] Stephen Blundell, Magnetism in Condensed
More informationCompeting Ferroic Orders The magnetoelectric effect
Competing Ferroic Orders The magnetoelectric effect Cornell University I would found an institution where any person can find instruction in any study. Ezra Cornell, 1868 Craig J. Fennie School of Applied
More informationMagnetism (FM, AFM, FSM)
Magnetism (FM, AFM, FSM) Karlheinz Schwarz Institute of Materials Chemistry TU Wien Localized vs. itinerant systems In localized systems (e.g. some rare earth) the magnetism is mainly governed by the atom
More informationTrivially parallel computing
Parallel Computing After briefly discussing the often neglected, but in praxis frequently encountered, issue of trivially parallel computing, we turn to parallel computing with information exchange. Our
More informationSpins and spin-orbit coupling in semiconductors, metals, and nanostructures
B. Halperin Spin lecture 1 Spins and spin-orbit coupling in semiconductors, metals, and nanostructures Behavior of non-equilibrium spin populations. Spin relaxation and spin transport. How does one produce
More informationCollective behavior, from particles to fields
978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
More information