A Brief History of the Use of the Metropolis Method at LANL in the 1950s
|
|
- Arlene Boone
- 6 years ago
- Views:
Transcription
1 A Brief History of the Use of the Metropolis Method at LANL in the 1950s William W. Wood Los Alamos National Laboratory Los Alamos, New Mexico Abstract. A personal view of the use of the Metropolis Algorithm in statistical mechanics calculations at Los Alamos during the 1950s will be presented, based on [1] and [2]. INTRODUCTION I came to Los Alamos in the fall of 1950 to work in the explosives division under Duncan P. Macdougall and Eugene H. Eyster. I had nothing to do with the invention of the Metropolis Algorithm, and indeed did nothing with it until several years later. I had done my doctoral work with Professor John G. Kirkwood at Cal Tech, who was then a consultant to the explosives division. He and I worked closely together until his death in 1959, and I came to admire him greatly. The first time I heard about what later would become known as the Metropolis Algorithm was when Marshall Rosenbluth gave me a preprint of [3] in early 1953, which led me to mention it briefly to Kirkwood in a letter dated March 3, But he may well have heard about it the previous summer during his consulting visit. Kirkwood quickly realized the potential of the method and urged that we should use it to calculate the equation of state of Lennard Jones molecules, that is to say, the equation of state of particles interacting with the Lennard Jones 6 12 potential V r ε r r 12 r r We had already been using that potential for some calculations based on the Lennard Jones cell theory, for the equation of state of detonation products. We began our work with the Metropolis Algorithm in late 1953 or early 1954, I do not recall which. (Many of my research notes of that period later fell victim to bureaucratic enthusiasm for reducing the amount of material in the files.) It was, of course, quite natural for us to begin our work in the belief that the Rosenbluths had done as well as could be done on the equation of state of hard spheres [4] with the available calculators of that time (MANIAC, in their case). It was a simple matter to modify a Lennard Jones code to calculate for hard spheres, and we did indeed make such a run as a check on the correctness of our program. Unfortunately, we happened to choose a density somewhat above what would later become known as the transition region for our IBM 704 run, and found reasonable agreement with the Rosenbluths results. 6 (1)
2 In the meantime, Kirkwood s enthusiasm for the Metropolis Algorithm had temporarily waned. This was due to the fact that the Rosenbluths results for the equation of state of hard spheres [4] disagreed badly with the results from the work that he and his students, including Alder, had done using the superposition approximation in solving the Born Green Kirkwood integral equation, which was a difficult calculation in its own right on the then-available calculators. To make matters worse, the simple cell theory gave somewhat better agreement with the Rosenbluths results than did the integral equation results. To see the nature of his objections, we will need to go slightly into the details of the Metropolis Algorithm. THE METROPOLIS ALGORITHM The essence of the Metropolis Algorithm is the generation of a certain kind of random walk in the configuration space of the system. The simplest kind of system to consider is one in which there is a finite number of possible configuration states, say M. Then the method consists of the generation of the type of random walk known as a Markov chain. Suppose we want to estimate the canonical average of some function f i of the state i of the system, f M 1 f i exp Q i 1 M βe i Q exp βe i (3) i 1 with β 1 k B T,andE i being the energy of state i. The realization average over S steps of the Markov chain is S S 1 f S f i n (4) n 1 i n being the state of the system on step n. Then S f (2) f (5) in probability as S tends to, provided p i j exp βe i p j i exp βe j (6) (microscopic reversibility), along with other important conditions regarding the ergodicity of the chain; p i j is the one-step probability that, if the system is in state i on step n, then it will be in state j on step n 1. One way in which the microscopic reversibility condition can be satisfied is to let p i j A i j π i j j i (7) with A i j being a symmetric stochastic matrix (i.e., M j 1 A i j 1, for all i),
3 and 1 if E π i j exp β E 0 if E 0; here E E j E i. Then the necessary normalization condition can be satisfied by letting for each i. (8) M p i i 1 j p 1! i j (9) i KIRKWOOD S OBJECTIONS For a system of hard spheres Eqs.(8 9) become and π i j 1 if E 0 0 if E p i i 1 f i (11) where f i is the fraction of states which can be reached from state i, but which result in an overlap between two (or more) hard spheres. It follows that one must count state i n again in Eq.(4) in order to satisfy the essential normalization condition, but it was this prescription to which Kirkwood first objected, in a letter to me on November 16, 1953, in which he attributed the opposite procedure (keep trying until a successful move is found, without counting state i n again) to Frankel and Lewinson. There then ensued a lively correspondence between Kirkwood and me, and between Kirkwood and Marshall Rosenbulth. Marshall wrote up a more detailed proof of the method, and the two Rosenbluths did a calculation for a one dimensional system of hard rods, finding good agreement with the (exact) Tonks equation of state. I did a calculation, both analytically and numerically, for a system of one hard rod confined between hard walls. These efforts temporarily convinced Kirkwood of the validity of the procedure, according to a letter to me on December 21,1953. But when I visited Kirkwood in New Haven in early January, 1954, he raised a second objection, this time to the bias in Eq.(8) in favor of a state of lower energy. He proposed that one should use instead the more symmetric π i j exp βe j exp βe i exp βe j (12) Both Marshall and I agreed that Eq.(12) was a valid procedure, but we intuitively thought Eq.(8) was a better choice. Finally, on April 24, 1954, Kirkwood wrote to me I finally was able to find the time to devote a weekend to thinking about the Monte Carlo problem. He now agreed that Eq.(8) was one of the correct ones, but not to let us (10)
4 off too lightly, he continued with Marshall s arguments were at no time helpful, and your modes of presentation were opaque. From that time on he was an active proponent of the method and was urgently prompting me to publish our results as soon as possible. By October, 1954, we had results for LJ systems on the k B T ε 2 74 isotherm over the fluid range V V $ 0 95, and by April, 1955, we had seen the first indication of the fluid solid phase transition, in the form of the now familiar two level structure of the pressure versus Markov chain time. Those results were not published until 1957 [5], because we got distracted by the disagreement of the Alder and Wainwright molecular dynamics results and the Rosenbluths results for hard spheres, as will be described in the next section. THE HARD SPHERE TRANSITION In August, 1956, there was a symposium on the theory of transport processes in Brussels, which I did not attend, at which Alder and Wainwright [6] presented a variety of results from their molecular dynamics calculations. As far as I am aware, this was the first public mention of the molecular dynamics method. There they called attention to the significant differences between their results and those of the Rosenbluths [4]. I was unaware of those differences until I received a letter from Alder, written on October 8, 1956, which mentioned it. As soon as I became aware of the discrepancy, I suspected that the disagreement was due to the shortness of the runs made by the Rosenbluths, typically only a few hundred attempted moves per particle, as well as to their choice of system size. We had already noticed that in order to obtain reliable results for Lennard Jones particles, something on the order of several thousand moves per particle were required. In November, 1956, Alder and Wainwright visited Los Alamos, and our conversations resulted in my agreeing that we would repeat the Rosenbluths calculations in the region around V V 0 1 6(V 0 being the close packed volume) where the disagreement was largest. There then ensued a period of intense collaboration between us at Los Alamos and Alder and Wainwright at Livermore. It was of course an easy matter to convert our Lennard Jones code to calculate for hard spheres, although the code was not very efficient for systems of more than N 32 particles, since it examined all N 1 interactions of the displaced particle. So we concentrated our efforts on systems of 32 hard spheres, until a more efficient code could be developed [7]. By mid-december, 1956, we had made Monte Carlo runs at V V 0 = 1.5, 1.6, and 1.7 for systems of N=32 and 108 particles and had observed the transient melting of the initial lattice configuration, interpreting the results as evidence of the existence of a first order phase transition in the hard sphere system. Runs started from the initial lattice configuration tended to remain in the initial ordered arrangement, with all particles remaining relatively near their initial positions, suggestive of the solid state, for some number of attempted moves per particle; but then, over a small number of moves per particle, the initial lattice would suddenly disorder into a state of relatively rapid diffusion, suggestive of the fluid state, with an accompanying increase in the computed pressure. Those results were communicated to Alder and Wainwright in letters written on
5 December 4, 1956, and January 4, 1957, and also, of course, to Kirkwood. Kirkwood was delighted, as he had previously predicted [8] such a transition. Alder and Wainwright were skeptical of this interpretation at first, because they were not able to observe the disordering of the initial lattice configuration for the larger systems they were using, and because they were using a somewhat slower UNIVAC calculator than the IBM 704. Both groups of course realized that such results were very far from definitively proving the existence of such a phase transition. In January, 1957, a symposium on the many body problem was held at the Stevens Institute of Technology. I was not present, but Kirkwood presented our results, including the phase transition interpretation. Alder argued against it, at the time. The proceedings were not published until 1963 [9], after Kirkwood s death in 1959, at which time Alder and Wainwright remarked in their published paper... but it is clear that some transition is occurring in that region, presumably a first order one from a fluid to a solid phase. In March, 1957, the molecular dynamics results and the Monte Carlo results at V V 0 = 1.6 and N = 32 were in serious disagreement, the latter giving higher values of the pressure, after discarding the initial, solid-like portion of the run. It turned out that the molecular dynamics runs were too short to disorder the initial lattice, as became clear after the two groups interchanged their final configurations, which were then used as the initial configurations for additional runs by the other group. In April, 1957, I arranged an invitation for Alder to come to Los Alamos that summer. He accepted, and brought with him a molecular dynamics code for the IBM 704. Between June and August, 1957, he used about 180 hours of time on those machines, during which time we prepared our presentations for the IUPAP meeting in Varenna, in September, which were published in Nuovo Cimento [10],[11] and as adjacent Letters to the Editor in the Journal of Chemical Physics [12],[13]. EPILOGUE Until recently, I had always recalled this period of friendly collaboration with Alder and Wainwright as a very pleasant one. Thus I was quite surprised to learn, after a colleague pointed it out to me, that Alder, at least, recalls it very differently: In the proceedings of a meeting in Alghero (Sardinia) in July 15 17, 1991, he writes [14]: So then we went back to our first love which was the hard sphere transition. I had to know whether Bill Wood had done his job right by Monte-Carlo and it turned out that he had not. By molecular dynamics we found a phase transition and subsequently Wood also found it and that s a long story. REFERENCES 1. Wood, W. W., Early History of Computer Simulations in Statistical Mechanics in Molecular Dynamics Simulation of Statistical Mechanics Systems, edited by G. Ciccotti and W. G. Hoover, North Holland, New York, 1986, pp Wood, W. W., Chapter 36, pp , in Monte Carlo and Molecular Dynamics of Condensed Matter Systems, edited by K. Binder and G. Ciccotti, Italian Physical Society, 1996.
6 3. Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., and Teller, E., J. Chem. Phys. 21, (1953). 4. Rosenbluth, M. N., and Rosenbluth, A. W., J. Chem. Phys. 22, (1954). 5. Wood, W. W., and Parker, F. R., J. Chem. Phys (1957). 6. Alder, B. J., and Wainwright, T. E., Molecular Dynamics by Electronic Computers, in International Symposium on the Statistical Mechanical Theory of Transport Processes, Brussels, 1956, edited by I. Prigogine, Interscience, New York, 1958, pp Wood, W. W., and Jacobson, J. D., Monte Carlo Calculations in Statistical Mechanics, in Proceedings of the Western Joint Computer Conference, San Francisco, 1959, pp Kirkwood, J. G., Crystallization as a Cooperative Phenomenon, in Phase Transformations in Solids, edited by R. Smoluchowski, J. E. Mayer and W. A. Weyl (Wiley, New York, 1951), pp Alder, B. J., and Wainwright, T. E., Investigation of the Many Body Problem by Electronic Computers, Chapter 29 in The Many Body Problem, edited by J. K. Percus, Interscience, New York, 1963, pp Wood, W. W., Parker, F. R., and Jacobson, J. D., Supplement to Nuovo Cimento 9, (Series 10), (1958). 11. Alder, B. J., and Wainwright, T. E., Supplement to Nuovo Cimento 9, (Series 10), (1958). 12. Wood, W. W., and Jacobson, J. D., J. Chem. Phys. 27, (1957). 13. Alder, B. J., and Wainwright, T. E., J. Chem. Phys. 27, (1957). 14. Alder, B. J., in Microscopic Simulations of Complex Hydrodynamic Phenomena, edited by M. Mareschal and B. L. Holian, Plenum Press, New York, 1992, pp
Copyright 2001 University of Cambridge. Not to be quoted or copied without permission.
Course MP3 Lecture 4 13/11/2006 Monte Carlo method I An introduction to the use of the Monte Carlo method in materials modelling Dr James Elliott 4.1 Why Monte Carlo? The name derives from the association
More informationTHE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM
Molecular Simulation, 1987, Vol. 1, pp. 87-93 c Gordon and Breach Science Publishers S.A. THE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM M. MEZEI Department of Chemistry, Hunter
More informationRandom Walks A&T and F&S 3.1.2
Random Walks A&T 110-123 and F&S 3.1.2 As we explained last time, it is very difficult to sample directly a general probability distribution. - If we sample from another distribution, the overlap will
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 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 informationNew Observations Show Carbon Emissions Don t Cause Global Warming
New Observations Show Carbon Emissions Don t Cause Global Warming Background Dr David Evans (david.evans@sciencespeak.com) 11 December 2007 Current Summary of Crucial Evidence A paper I wrote that briefly
More informationAn extended summary of the NCGR/Berkeley Double-Blind Test of Astrology undertaken by Shawn Carlson and published in 1985
From: http://www.astrodivination.com/moa/ncgrberk.htm An extended summary of the NCGR/Berkeley Double-Blind Test of Astrology undertaken by Shawn Carlson and published in 1985 Introduction Under the heading
More informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
More informationMetropolis/Variational Monte Carlo. Microscopic Theories of Nuclear Structure, Dynamics and Electroweak Currents June 12-30, 2017, ECT*, Trento, Italy
Metropolis/Variational Monte Carlo Stefano Gandolfi Los Alamos National Laboratory (LANL) Microscopic Theories of Nuclear Structure, Dynamics and Electroweak Currents June 12-30, 2017, ECT*, Trento, Italy
More informationOPTIMIZATION BY SIMULATED ANNEALING: A NECESSARY AND SUFFICIENT CONDITION FOR CONVERGENCE. Bruce Hajek* University of Illinois at Champaign-Urbana
OPTIMIZATION BY SIMULATED ANNEALING: A NECESSARY AND SUFFICIENT CONDITION FOR CONVERGENCE Bruce Hajek* University of Illinois at Champaign-Urbana A Monte Carlo optimization technique called "simulated
More informationSession 3A: Markov chain Monte Carlo (MCMC)
Session 3A: Markov chain Monte Carlo (MCMC) John Geweke Bayesian Econometrics and its Applications August 15, 2012 ohn Geweke Bayesian Econometrics and its Session Applications 3A: Markov () chain Monte
More informationFibonacci Numbers. November 7, Fibonacci's Task: Figure out how many pairs of rabbits there will be at the end of one year, following rules.
Fibonacci Numbers November 7, 2010 Fibonacci's Task: Figure out how many pairs of rabbits there will be at the end of one year, following rules. Rules: 1. Start with a pair of new rabbits, born in December.
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 informationMITOCW MITRES_18-007_Part3_lec5_300k.mp4
MITOCW MITRES_18-007_Part3_lec5_300k.mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationComputational statistics
Computational statistics Markov Chain Monte Carlo methods Thierry Denœux March 2017 Thierry Denœux Computational statistics March 2017 1 / 71 Contents of this chapter When a target density f can be evaluated
More informationVordiplom (German B.S.), Physics, Summer 2001; GPA: 3.7 Universität des Saarlandes, Saarbrücken, Germany
Home address Steven Diehl 621 Carriage Hill Drive Athens, OH 45701, U.S. Phone: 740-594-1025 diehl@helios.phy.ohiou.edu http://www.phy.ohiou.edu/ diehl Work address Astrophysical Institute Department of
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 informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationA new algorithm for Reverse Monte Carlo simulations
A new algorithm for Reverse Monte Carlo simulations Fernando Lus B. da Silva, Bo Svensson, Torbjörn Åkesson, and Bo Jönsson Citation: The Journal of Chemical Physics 109, 2624 (1998); doi: 10.1063/1.476861
More informationDynamics of Solitary Waves Induced by Shock Impulses in a Linear Atomic Chain*
Dynamics of Solitary Waves Induced by Shock Impulses in a Linear Atomic Chain* PHUOC X. TRAN, DONALD W. BRENNER, and C. T. WHITE Naval Research Laboratory, Washington, DC 20375-5000 Abstract The propagation
More informationPhase transitions of quadrupolar fluids
Phase transitions of quadrupolar fluids Seamus F. O Shea Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada, T1K 3M4 Girija S. Dubey Brookhaven National Laboratory, Upton, New
More informationIsaac Newton Benjamin Franklin Michael Faraday
Isaac Newton (4 January 1643 31 March 1727) was born and raised in England. He was a greater thinker and made many discoveries in physics, mathematics, and astronomy. Newton was the first to describe the
More informationCHAPTER 2--THE DEVELOPMENT OF EVOLUTIONARY THEORY
CHAPTER 2--THE DEVELOPMENT OF EVOLUTIONARY THEORY Student: 1. In Europe during the Middle Ages, it was believed that. A. all species had evolved from a common ancestor B. evolution was the result of natural
More informationEquation of State of Dense Helium
Iowa State University From the SelectedWorks of Richard Alan Lesar October 31, 1988 Equation of State of Dense Helium Richard Alan Lesar, Los Alamos National Laboratory Available at: https://works.bepress.com/richard_lesar/27/
More informationAn Introduction to Two Phase Molecular Dynamics Simulation
An Introduction to Two Phase Molecular Dynamics Simulation David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: April 19, 2016 date last updated: April
More informationLecture 3: From Random Walks to Continuum Diffusion
Lecture 3: From Random Walks to Continuum Diffusion Martin Z. Bazant Department of Mathematics, MIT February 3, 6 Overview In the previous lecture (by Prof. Yip), we discussed how individual particles
More informationParadoxes of special relativity
Paradoxes of special relativity Today we are turning from metaphysics to physics. As we ll see, certain paradoxes about the nature of space and time result not from philosophical speculation, but from
More informationTime-Dependent Statistical Mechanics 1. Introduction
Time-Dependent Statistical Mechanics 1. Introduction c Hans C. Andersen Announcements September 24, 2009 Lecture 1 9/22/09 1 Topics of concern in the course We shall be concerned with the time dependent
More informationAstronomy 10 Introduction to the Solar System
Astronomy 10 Introduction to the Solar System http://lpc1.clpccd.cc.ca.us/lpc/murray/astro10/a10intro.htm This is your source for all things in this class: 1. Grading policies. 2. What topics we ll covering.
More informationVerification Of January HDD Forecasts
Verification Of January HDD Forecasts W2020 / Average HDD stands for Heating Degree Day. A Heating Degree Day is zero if the average temperature is 65 degrees. An HDD of -30 would mean an average temperature
More informationThe Dynamics of Continued Fractions
The Dynamics of Continued Fractions Evan O Dorney May 3, 20 The Story I was first introduced to the Intel Science Talent Search in ninth grade. I knew I would have no trouble entering this contest, as
More informationCARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores
CARBON Providence, Rhode Island Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores D. D. Do* and H. D. Do, University of Queensland, St. Lucia, Qld 7, Australia
More informationAn Interview with ARTHUR W. BURKS OH 136. Conducted by William Aspray. 20 June Ann Arbor, MI
An Interview with ARTHUR W. BURKS OH 136 Conducted by William Aspray on 20 June 1987 Ann Arbor, MI Charles Babbage Institute The Center for the History of Information Processing University of Minnesota,
More informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
More informationCOMA CLUSTER OF GALAXIES
COMA CLUSTER OF GALAXIES Learn the basics of galaxy classification and grouping, using actual Hubble Space Telescope images. Keely Finkelstein, McDonald Observatory Curriculum topic Galaxies Big idea of
More informationIntroducing Proof 1. hsn.uk.net. Contents
Contents 1 1 Introduction 1 What is proof? 1 Statements, Definitions and Euler Diagrams 1 Statements 1 Definitions Our first proof Euler diagrams 4 3 Logical Connectives 5 Negation 6 Conjunction 7 Disjunction
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 informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
More informationChemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. # 26 Problem solving : Heterogeneous reactions Friends, in last few
More informationWhen you throw a ball, increasing the speed of your hand increases the speed of the ball. Does the same apply to a wave pulse?
Tutorial Finding mechanism: Pulses Name I. One traveling pulse: The effect of flicking the spring faster When you throw a ball, increasing the speed of your hand increases the speed of the ball. Does the
More informationDetection ASTR ASTR509 Jasper Wall Fall term. William Sealey Gosset
ASTR509-14 Detection William Sealey Gosset 1876-1937 Best known for his Student s t-test, devised for handling small samples for quality control in brewing. To many in the statistical world "Student" was
More informationQuantum Entanglement. Chapter Introduction. 8.2 Entangled Two-Particle States
Chapter 8 Quantum Entanglement 8.1 Introduction In our final chapter on quantum mechanics we introduce the concept of entanglement. This is a feature of two-particle states (or multi-particle states) in
More informationRemembering Adelchi Fabrocini ( )
Remembering Adelchi Fabrocini (1951 2006) Born on August 27, 1951 in Sepino (Campobasso) June 1977: Laurea in Physics at the University of Pisa with the degree summa cum laude. Thesis title: Variational
More informationScientific Revolution
Age of Revolutions Scientific Revolution Scientific Revolution Period of time in which a new way of thinking came about. The beliefs held by many for so long were now being questioned. Use logic and reason
More informationAlthough we rarely see a shock wave in everyday life,
Shear Viscosity under Shock-Loading Conditions Zhuohan Liang 1 and Steven M. Valone 1 T-4, Los Alamos Summer School MST-8, Los Alamos National Laboratory Although we rarely see a shock wave in everyday
More informationUse of Eigen values and eigen vectors to calculate higher transition probabilities
The Lecture Contains : Markov-Bernoulli Chain Note Assignments Random Walks which are correlated Actual examples of Markov Chains Examples Use of Eigen values and eigen vectors to calculate higher transition
More informationMarkov Chain Monte Carlo The Metropolis-Hastings Algorithm
Markov Chain Monte Carlo The Metropolis-Hastings Algorithm Anthony Trubiano April 11th, 2018 1 Introduction Markov Chain Monte Carlo (MCMC) methods are a class of algorithms for sampling from a probability
More informationPENGUIN READERS. Five Famous Fairy Tales
PENGUIN READERS Five Famous Fairy Tales Introduction Jacob and Wilhelm Grimm the Brothers Grimm were good friends. Jacob was a quiet man and sometimes sad. Wilhelm was often very ill but he was a happier
More informationMinicourse on: Markov Chain Monte Carlo: Simulation Techniques in Statistics
Minicourse on: Markov Chain Monte Carlo: Simulation Techniques in Statistics Eric Slud, Statistics Program Lecture 1: Metropolis-Hastings Algorithm, plus background in Simulation and Markov Chains. Lecture
More informationMath 456: Mathematical Modeling. Tuesday, April 9th, 2018
Math 456: Mathematical Modeling Tuesday, April 9th, 2018 The Ergodic theorem Tuesday, April 9th, 2018 Today 1. Asymptotic frequency (or: How to use the stationary distribution to estimate the average amount
More informationBias-Variance Error Bounds for Temporal Difference Updates
Bias-Variance Bounds for Temporal Difference Updates Michael Kearns AT&T Labs mkearns@research.att.com Satinder Singh AT&T Labs baveja@research.att.com Abstract We give the first rigorous upper bounds
More informationWhen Is A Linear Operator Diagonalizable? Marco Abate
When Is A Linear Operator Diagonalizable? Marco Abate 0 Introduction As it often happens, everything began with a mistake I was teaching for the third year in a row a linear algebra course to engineering
More informationFILE NAME: PYTHAGOREAN_TRIPLES_012-( WED).DOCX AUTHOR: DOUG JONES
FILE NAME: PYTHAGOREAN_TRIPLES_01-(0090107WED.DOCX AUTHOR: DOUG JONES A. BACKGROUND & PURPOSE THE SEARCH FOR PYTHAGOREAN TRIPLES 1. Three positive whole numbers ( a,b,c which are such that a + b = c are
More informationA Mathematical Analysis of the Collapse of the World Trade Center
by Kenneth A. Kuhn Oct. 31, 2008, rev. Nov. 8, 2008 Introduction There has been much debate concerning the speed at which the twin towers collapsed on Sept. 11, 2001. Many have felt that the near free
More informationIntroduction. Model DENSITY PROFILES OF SEMI-DILUTE POLYMER SOLUTIONS NEAR A HARD WALL: MONTE CARLO SIMULATION
169 DENSITY PROFILES OF SEMI-DILUTE POLYMER SOLUTIONS NEAR A HARD WALL: MONTE CARLO SIMULATION WAN Y. SHIH, WEI-HENG SHIH and ILHAN A. AKSAY Dept. of Materials Science and Engineering University of Washington,
More informationarxiv:gr-qc/ v2 30 Oct 2005
International Journal of Modern Physics D c World Scientific Publishing Company arxiv:gr-qc/0505111v2 30 Oct 2005 ENTROPY AND AREA IN LOOP QUANTUM GRAVITY JOHN SWAIN Department of Physics, Northeastern
More informationThere are self-avoiding walks of steps on Z 3
There are 7 10 26 018 276 self-avoiding walks of 38 797 311 steps on Z 3 Nathan Clisby MASCOS, The University of Melbourne Institut für Theoretische Physik Universität Leipzig November 9, 2012 1 / 37 Outline
More informationLab 1: Electric Interactions I
Page 1 of 6 LEARNING GOALS OF THE LAB: 1. Learn how to explain the interactions of charged and uncharged objects. 2. Learn how to test a 3. Find how the strength of the electric interaction depends on
More informationHugh Everett III s Many Worlds
236 My God, He Plays Dice! Hugh Everett III s Many Worlds Many Worlds 237 Hugh Everett III s Many Worlds Hugh Everett III was one of John Wheeler s most famous graduate students. Others included Richard
More informationProtein Folding Challenge and Theoretical Computer Science
Protein Folding Challenge and Theoretical Computer Science Somenath Biswas Department of Computer Science and Engineering, Indian Institute of Technology Kanpur. (Extended Abstract) September 2011 Almost
More informationBay Area Scientists in Schools Presentation Plan
Bay Area Scientists in Schools Presentation Plan Lesson Name Presenter(s) Grade Level 3rd Melting, Freezing, and More!: Phase Transitions Steven Scroggins, Ailey Crow, Tom Holcombe, and Terence Choy California
More informationThe Quantum-Classical Transition and Wave Packet Dispersion. C. L. Herzenberg
The Quantum-Classical Transition and Wave Packet Dispersion C. L. Herzenberg Abstract Two recent studies have presented new information relevant to the transition from quantum behavior to classical behavior,
More informationFeatured Alumna Sarah Caudill ( 06)
Featured Alumna Sarah Caudill ( 06) As a high school student applying to colleges, I had a choice between Stetson and the University of Florida. I reasoned that I would receive a more personalized education
More informationChapter 13 - Inverse Functions
Chapter 13 - Inverse Functions In the second part of this book on Calculus, we shall be devoting our study to another type of function, the exponential function and its close relative the Sine function.
More informationSpecial Theory Of Relativity Prof. Shiva Prasad Department of Physics Indian Institute of Technology, Bombay
Special Theory Of Relativity Prof. Shiva Prasad Department of Physics Indian Institute of Technology, Bombay Lecture - 6 Length Contraction and Time Dilation (Refer Slide Time: 00:29) In our last lecture,
More informationMICROPROCESSOR REPORT. THE INSIDER S GUIDE TO MICROPROCESSOR HARDWARE
MICROPROCESSOR www.mpronline.com REPORT THE INSIDER S GUIDE TO MICROPROCESSOR HARDWARE ENERGY COROLLARIES TO AMDAHL S LAW Analyzing the Interactions Between Parallel Execution and Energy Consumption By
More informationQ7K~ 45-ti FINAL REPORT. ~Cmm~STUDIES. V11 w rte; 1- Contract AF 49(638)-965 IN HYPERSONIC FLOW THEORY
Q7K~ 45-ti o FINAL REPORT Contract AF 49(638)-965 ~Cmm~STUDIES IN HYPERSONIC FLOW THEORY -- T V11 w rte; 1- Department of Aeronautics and Astronautics Stanford University Stanford, California Milton Van
More informationLecture 2+3: Simulations of Soft Matter. 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications
Lecture 2+3: Simulations of Soft Matter 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications D. Frenkel, Boulder, July 6, 2006 What distinguishes Colloids from atoms or
More informationMotivation, Basic Concepts, Basic Methods, Travelling Salesperson Problem (TSP), Algorithms
Motivation, Basic Concepts, Basic Methods, Travelling Salesperson Problem (TSP), Algorithms 1 What is Combinatorial Optimization? Combinatorial Optimization deals with problems where we have to search
More informationMarkov Processes. Stochastic process. Markov process
Markov Processes Stochastic process movement through a series of well-defined states in a way that involves some element of randomness for our purposes, states are microstates in the governing ensemble
More informationPeople. The Shadow Shadow People. The Shadow People A Reading A Z Level O Leveled Book Word Count: 874 LEVELED BOOK O.
The Shadow People A Reading A Z Level O Leveled Book Word Count: 874 LEVELED BOOK O The Shadow Shadow The People People Written by Ann Weil Illustrated by Stephen Marchesi Visit www.readinga-z.com for
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 informationSupplementary Note on Bayesian analysis
Supplementary Note on Bayesian analysis Structured variability of muscle activations supports the minimal intervention principle of motor control Francisco J. Valero-Cuevas 1,2,3, Madhusudhan Venkadesan
More informationPressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster
University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University
More information13.7 ANOTHER TEST FOR TREND: KENDALL S TAU
13.7 ANOTHER TEST FOR TREND: KENDALL S TAU In 1969 the U.S. government instituted a draft lottery for choosing young men to be drafted into the military. Numbers from 1 to 366 were randomly assigned to
More informationCHAPTER 1. Introduction
CHAPTER 1 Introduction A typical Modern Geometry course will focus on some variation of a set of axioms for Euclidean geometry due to Hilbert. At the end of such a course, non-euclidean geometries (always
More informationStatistical Mechanics
42 My God, He Plays Dice! Statistical Mechanics Statistical Mechanics 43 Statistical Mechanics Statistical mechanics and thermodynamics are nineteenthcentury classical physics, but they contain the seeds
More informationComplex Numbers For High School Students
Complex Numbers For High School Students For the Love of Mathematics and Computing Saturday, October 14, 2017 Presented by: Rich Dlin Presented by: Rich Dlin Complex Numbers For High School Students 1
More informationThe PROMYS Math Circle Problem of the Week #3 February 3, 2017
The PROMYS Math Circle Problem of the Week #3 February 3, 2017 You can use rods of positive integer lengths to build trains that all have a common length. For instance, a train of length 12 is a row of
More informationMONTE CARLO METHOD. Reference1: Smit Frenkel, Understanding molecular simulation, second edition, Academic press, 2002.
MONTE CARLO METHOD Reference1: Smit Frenkel, Understanding molecular simulation, second edition, Academic press, 2002. Reference 2: David P. Landau., Kurt Binder., A Guide to Monte Carlo Simulations in
More informationMarkov Chain Monte Carlo
Chapter 5 Markov Chain Monte Carlo MCMC is a kind of improvement of the Monte Carlo method By sampling from a Markov chain whose stationary distribution is the desired sampling distributuion, it is possible
More informationPhysics 562: Statistical Mechanics Spring 2002, James P. Sethna Prelim, due Wednesday, March 13 Latest revision: March 22, 2002, 10:9
Physics 562: Statistical Mechanics Spring 2002, James P. Sethna Prelim, due Wednesday, March 13 Latest revision: March 22, 2002, 10:9 Open Book Exam Work on your own for this exam. You may consult your
More informationUNIT 2 The Particulate Nature of Matter
UNIT 2 The Particulate Nature of Matter Take a moment to think about all of the different properties which can be exhibited by matter. The list seems endless. Matter can be found in a variety of forms,
More informationDepartment of Physics & Astronomy Trent University
Department of Physics & Astronomy Trent University PHYS 2610H : INTRODUCTORY QUANTUM PHYSICS 2017FA PETERBOROUGH Instructor: Balaji Subramanian Campus: Peterborough balajisubramanian@trentu.ca SC 322 15063375166
More information12. LOCAL SEARCH. gradient descent Metropolis algorithm Hopfield neural networks maximum cut Nash equilibria
Coping With NP-hardness Q. Suppose I need to solve an NP-hard problem. What should I do? A. Theory says you re unlikely to find poly-time algorithm. Must sacrifice one of three desired features. Solve
More informationThe problem of transition from school to university mathematics Student Survey E. Krause, F. Wetter, C. Nguyen Phuong (2015)
Demographic Data: Please answer the following questions. 1. Sex: Male Female 2. Age (years): 3. Semester at university: 4. In which country / state, or in which provinces did you go to school? 5. What
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 informationMonte Carlo Simulations for a Soft Sphere Fluid
Monte Carlo Simulations for a Soft Sphere Fluid Patrick Kreitzberg Advisor: Dr. David Roundy May 8 th, 2015 1 Contents 1 2 3 5 7 Abstract Introduction Methods Results and Discussion Conclusion Acknowledgments
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 informationBIBLIOTHECA MATHEMATICA
BIBLIOTHECA MATHEMATICA A series of Monographs on Pure and Applied Mathematics Volume VIII Edited with the cooperation of THE 'MATHEMATISCH CENTRUM' and THE 'WISKUNDIG GENOOTSCHAP' at Amsterdam Editors:
More informationImperfect Gases. NC State University
Chemistry 431 Lecture 3 Imperfect Gases NC State University The Compression Factor One way to represent the relationship between ideal and real gases is to plot the deviation from ideality as the gas is
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationSUBTRACTIVE BLACK HOLES AND BLACK LOOPS
Texas College Mathematics Journal Volume 2, Number 1, Pages 1-9 Article electronically published on August 17, 2005 SUBTRACTIVE BLACK HOLES AND BLACK LOOPS MURRAY H. SIEGEL, WASIN SO AND PETER A. COOPER
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 informationOne day an ant was drinking at a small stream and fell in. She made desperate
(A) One day an ant was drinking at a small stream and fell in. She made desperate efforts to reach the side, but made no progress at all. The poor ant almost exhausted was still bravely doing her best
More informationChapter 16: DC Circuits
Chapter 16: DC Circuits Why might all the electrical devices in your house suddenly turn off if you simultaneously turn too many on? How can you use an electric circuit to model the circulatory system?
More informationThe Next Generation Science Standards (NGSS)
The Next Generation Science Standards (NGSS) CHAPTER 3, LESSON 1: WHAT IS DENSITY? MS-PS1-1. Develop models to describe the atomic composition of simple molecules and extended structures. DISCIPLINARY
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 informationHow the maximum step size in Monte Carlo simulations should be adjusted
Physics Procedia Physics Procedia 00 (2013) 1 6 How the maximum step size in Monte Carlo simulations should be adjusted Robert H. Swendsen Physics Department, Carnegie Mellon University, Pittsburgh, PA
More information( ) ( ) ( ) ( ) Simulated Annealing. Introduction. Pseudotemperature, Free Energy and Entropy. A Short Detour into Statistical Mechanics.
Aims Reference Keywords Plan Simulated Annealing to obtain a mathematical framework for stochastic machines to study simulated annealing Parts of chapter of Haykin, S., Neural Networks: A Comprehensive
More information