DOCSLIB.ORG

Lecture #8: Monte Carlo Method

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text

Load more

Lecture #8: Monte Carlo method ENGG304: Uncertainty, Reliability and Risk Edoardo Patelli Institute for Risk and Uncertainty E: [email protected] W: www.liv.ac.uk/risk-and-uncertainty T: +44 01517944079 A MEMBER OF THE RUSSELL GROUP Edoardo Patelli University of Liverpool 18 March 2019 1 / 81 Lecture Outline 1 Introduction 2 Monte Carlo method Random Number Generator Sampling Methods Buffon’s experiment Monte Carlo integration Probability of failure 3 Summary 4 Computer based class Some useful slides Assignments Edoardo Patelli University of Liverpool 18 March 2019 2 / 81 Introduction Programme ENGG304 1 28/01/2019 Introduction 2 08/02/2019 Human error 3 11/02/2019 Qualitative risk assessment: Safety analysis 4 18/02/2019 Qualitative risk assessment: Event Tree and Fault Tree 5 25/02/2019 Tutorial I 6 04/03/2019 Model of random phenomena 7 11/03/2019 Structural reliability 8 18/03/2019 Monte Carlo simulation I + Hands-on session 9 25/03/2019 Tutorial II 10 01/04/2019 Monte Carlo simulation II + Tutorial III (Hands-on session) Edoardo Patelli University of Liverpool 18 March 2019 3 / 81 Introduction Summary Lecture #7 Safety Margin Fundamental problem Performance function defined as “Capacity” - “Demand” Safety margin: M = C − D = g(x) 2 For normal and independent random variables: M ∼ N(µM ; σM ) q 2 2 µM = µC − µD σM = σC + σD Reliability index β β = µM /σM represents the number of standard deviations by which the mean value of the safety margin M exceeds zero Edoardo Patelli University of Liverpool 18 March 2019 4 / 81 Introduction Summary Lecture #7 Safety Margin Fundamental problem Performance function defined as “Capacity” - “Demand” Safety margin: M = C − D = g(x) 2 For normal and independent random variables: M ∼ N(µM ; σM ) q 2 2 µM = µC − µD σM = σC + σD Reliability index β β = µM /σM represents the number of standard deviations by which the mean value of the safety margin M exceeds zero Edoardo Patelli University of Liverpool 18 March 2019 4 / 81 Introduction Summary Lecture #7 Reliability index Geometrical Interpretation/1 Edoardo Patelli University of Liverpool 18 March 2019 5 / 81 Monte Carlo method Introduction Lecture Outline 1 Introduction 2 Monte Carlo method Random Number Generator Sampling Methods Buffon’s experiment Monte Carlo integration Probability of failure 3 Summary 4 Computer based class Some useful slides Assignments Edoardo Patelli University of Liverpool 18 March 2019 6 / 81 Monte Carlo method Introduction What is the Monte Carlo method? Edoardo Patelli University of Liverpool 18 March 2019 7 / 81 Monte Carlo method Introduction Monte Carlo method Computation technique based on random numbers Numerical experiment by generating a random sequence of number with prescribed probability distribution. Collecting quantity of interest Edoardo Patelli University of Liverpool 18 March 2019 8 / 81 Monte Carlo method Introduction Monte Carlo method How did Monte Carlo simulation get its name? Name derived from the Principality of Monaco Edoardo Patelli University of Liverpool 18 March 2019 9 / 81 Monte Carlo method Introduction Monte Carlo method How did Monte Carlo simulation get its name? Name derived from the Principality of Monaco Edoardo Patelli University of Liverpool 18 March 2019 9 / 81 Monte Carlo method Introduction Monte Carlo method Origins 1777 Comte de Buffon - earliest documented use of random sampling Edoardo Patelli University of Liverpool 18 March 2019 10 / 81 Monte Carlo method Introduction Monte Carlo method Origins 1777 Comte de Buffon - earliest documented use of random sampling 2 ∗ L P(needle intersects the grid) = πt Edoardo Patelli University of Liverpool 18 March 2019 10 / 81 Monte Carlo method Introduction Monte Carlo method Origins 1777 Comte de Buffon - earliest documented use of random sampling 2 ∗ L P(needle intersects the grid) = πt Edoardo Patelli University of Liverpool 18 March 2019 10 / 81 Monte Carlo method Introduction Monte Carlo method Origins 1777 Comte de Buffon - earliest documented use of random sampling 2 ∗ L P(needle intersects the grid) = πt Edoardo Patelli University of Liverpool 18 March 2019 10 / 81 Monte Carlo method Introduction Monte Carlo method Origins 1777 Comte de Buffon - earliest documented use of random sampling 2 ∗ L P(needle intersects the grid) = πt 1786 Laplace suggested to estimate π by random sampling Edoardo Patelli University of Liverpool 18 March 2019 10 / 81 Monte Carlo method Introduction Monte Carlo method Origins 1777 Comte de Buffon - earliest documented use of random sampling 2 ∗ L P(needle intersects the grid) = πt 1786 Laplace suggested to estimate π by random sampling http://mste.illinois.edu/reese/ buffon/bufjava.html Edoardo Patelli University of Liverpool 18 March 2019 10 / 81 Monte Carlo method Introduction Monte Carlo method Modern time 1930s, Enrico Fermi first experimented with the Monte Carlo Fermiac tools (transport of neutrons) Edoardo Patelli University of Liverpool 18 March 2019 11 / 81 Monte Carlo method Introduction Monte Carlo method Manhattan Project 1942-1946 Calculation of the explosive yield of atomic bomb Stanislaw Ulam suggested to used ENIAC for “method of statistic trials” Metropolis suggested the name “Monte Carlo” Von Neumann developed the first computer code (1947) Edoardo Patelli University of Liverpool 18 March 2019 12 / 81 Monte Carlo method Introduction Is it really used? Edoardo Patelli University of Liverpool 18 March 2019 13 / 81 Monte Carlo method Examples Monte Carlo applications Solution of integrals, differential equation, complex systems etc... Simulation random events Cryptography, Decision-Making Games Edoardo Patelli University of Liverpool 18 March 2019 14 / 81 Monte Carlo method Examples Monte Carlo applications Solution of integrals, differential equation, complex systems etc... Simulation random events Cryptography, Decision-Making Games Edoardo Patelli University of Liverpool 18 March 2019 14 / 81 Monte Carlo method Examples Monte Carlo applications Solution of integrals, differential equation, complex systems etc... Simulation random events Cryptography, Decision-Making Games Edoardo Patelli University of Liverpool 18 March 2019 14 / 81 Monte Carlo method Examples Monte Carlo applications Solution of integrals, differential equation, complex systems etc... Simulation random events Cryptography, Decision-Making Games Edoardo Patelli University of Liverpool 18 March 2019 14 / 81 Monte Carlo method Main components Monte Carlo method Main components Random number generator Probability distribution functions (describing the model) Sampling rules (how to sample from PDFs) Error estimator Variance reduction technique Computational resources Edoardo Patelli University of Liverpool 18 March 2019 15 / 81 Monte Carlo method Random Number Generator Lecture Outline 1 Introduction 2 Monte Carlo method Random Number Generator Sampling Methods Buffon’s experiment Monte Carlo integration Probability of failure 3 Summary 4 Computer based class Some useful slides Assignments Edoardo Patelli University of Liverpool 18 March 2019 16 / 81 Monte Carlo method Random Number Generator Random Number Generator Definition A tool or method able to generate an (automatic) sequence (of number) in which each term is unpredictable Edoardo Patelli University of Liverpool 18 March 2019 17 / 81 Monte Carlo method Random Number Generator Type of Random Number Generators Definitions True random Each element has equal probability of being chosen from the set Pseudo-random A finite set that display qualities of random numbers Edoardo Patelli University of Liverpool 18 March 2019 18 / 81 Monte Carlo method Random Number Generator Type of Random Number Generators True Random Number Generator vs Pseudo Random Number Generator Pseudo-Random True Random Number Number Generators Generators Characteristic Excellent efficiency Poor efficiency Deterministic Non-deterministic Periodic Aperiodic Applications Edoardo Patelli University of Liverpool 18 March 2019 19 / 81 Monte Carlo method Random Number Generator Type of Random Number Generators True Random Number Generator vs Pseudo Random Number Generator Pseudo-Random True Random Number Number Generators Generators Characteristic Excellent efficiency Poor efficiency Deterministic Non-deterministic Periodic Aperiodic Applications Simulation and Mod- elling, Video Games, Security Edoardo Patelli University of Liverpool 18 March 2019 19 / 81 Monte Carlo method Random Number Generator Type of Random Number Generators True Random Number Generator vs Pseudo Random Number Generator Pseudo-Random True Random Number Number Generators Generators Characteristic Excellent efficiency Poor efficiency Deterministic Non-deterministic Periodic Aperiodic Applications Simulation and Mod- Lottery and Draws, elling, Video Games, Random Sampling, Security Games and Gambling Edoardo Patelli University of Liverpool 18 March 2019 19 / 81 Monte Carlo method Random Number Generator Pseudo-Random Number Generator Anyone who considers arithmetical methods of producing random digits is, of course, in a state of sin (John Von Neumann, 1951) Edoardo Patelli University of Liverpool 18 March 2019 20 / 81 Monte Carlo method Random Number Generator Pseudo-Random Number Generator Characteristics Approximate the PDF Statistically independent Should pass a battery of statistical tests Be aperiodic Difficult to prove whether a sequence of numbers is random Edoardo Patelli University of Liverpool 18 March 2019 21 / 81 Monte Carlo method Random Number Generator Pseudo-Random Number Generator Characteristics Approximate
Recommended publications
  • The Convergence of Markov Chain Monte Carlo Methods 3

    The Convergence of Markov Chain Monte Carlo Methods 3

    The Convergence of Markov chain Monte Carlo Methods: From the Metropolis method to Hamiltonian Monte Carlo Michael Betancourt From its inception in the 1950s to the modern frontiers of applied statistics, Markov chain Monte Carlo has been one of the most ubiquitous and successful methods in statistical computing. The development of the method in that time has been fueled by not only increasingly difficult problems but also novel techniques adopted from physics. In this article I will review the history of Markov chain Monte Carlo from its inception with the Metropolis method to the contemporary state-of-the-art in Hamiltonian Monte Carlo. Along the way I will focus on the evolving interplay between the statistical and physical perspectives of the method. This particular conceptual emphasis, not to mention the brevity of the article, requires a necessarily incomplete treatment. A complementary, and entertaining, discussion of the method from the statistical perspective is given in Robert and Casella (2011). Similarly, a more thorough but still very readable review of the mathematics behind Markov chain Monte Carlo and its implementations is given in the excellent survey by Neal (1993). I will begin with a discussion of the mathematical relationship between physical and statistical computation before reviewing the historical introduction of Markov chain Monte Carlo and its first implementations. Then I will continue to the subsequent evolution of the method with increasing more sophisticated implementations, ultimately leading to the advent of Hamiltonian Monte Carlo. 1. FROM PHYSICS TO STATISTICS AND BACK AGAIN At the dawn of the twentieth-century, physics became increasingly focused on under- arXiv:1706.01520v2 [stat.ME] 10 Jan 2018 standing the equilibrium behavior of thermodynamic systems, especially ensembles of par- ticles.
  • Computer History a Look Back Contents

    Computer History a Look Back Contents

    Computer History A look back Contents 1 Computer 1 1.1 Etymology ................................................. 1 1.2 History ................................................... 1 1.2.1 Pre-twentieth century ....................................... 1 1.2.2 First general-purpose computing device ............................. 3 1.2.3 Later analog computers ...................................... 3 1.2.4 Digital computer development .................................. 4 1.2.5 Mobile computers become dominant ............................... 7 1.3 Programs ................................................. 7 1.3.1 Stored program architecture ................................... 8 1.3.2 Machine code ........................................... 8 1.3.3 Programming language ...................................... 9 1.3.4 Fourth Generation Languages ................................... 9 1.3.5 Program design .......................................... 9 1.3.6 Bugs ................................................ 9 1.4 Components ................................................ 10 1.4.1 Control unit ............................................ 10 1.4.2 Central processing unit (CPU) .................................. 11 1.4.3 Arithmetic logic unit (ALU) ................................... 11 1.4.4 Memory .............................................. 11 1.4.5 Input/output (I/O) ......................................... 12 1.4.6 Multitasking ............................................ 12 1.4.7 Multiprocessing .........................................
  • RN Random Numbers

    RN Random Numbers

    Random Numb ers RN Copyright C June Computational Science Education Pro ject Remarks Keywords Random pseudorandom linear congruential lagged Fib onacci List of prerequisites some exp osure to sequences and series some abilitytowork in base arithmetic helpful go o d background in FORTRAN or some other pro cedural language List of computational metho ds pseudorandom numb ers random numb ers linear congruential generators LCGs lagged Fib onacci generators LFGs List of architecturescomputers any computer with a FORTRAN compiler and bitwise logical functions List of co des supplied ranlcfor linear congruential generator ranlffor lagged Fib onacci generator randomfor p ortable linear congruential generator after Park and Miller getseedfor to generate initial seeds from the time and date implemented for Unix and IBM PCs Scop e lectures dep ending up on depth of coverage Notation Key a integer multiplier a m c integer constant c m D length of needle in Buons needle exp eriment F fraction of trials needle falls within ruled grid in Buons needle exp eriment k lag in lagged Fib onacci generator k lag in lagged Fib onacci generator k M numb er of binary bits in m m integer mo dulus m N numb er of trials N numb er of parallel pro cessors p P p erio d of generator p a prime integer R random real numb er R n n S spacing b etween grid in Buons needle exp eriment th X n random integer n Subscripts and Sup erscripts denotes initial value or seed th j denotes j value th denotes value th k denotes k value th n denotes n value
  • Markov Chain Monte Carlo Methods

    Markov Chain Monte Carlo Methods

    Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos Contents Markov Chain Monte Carlo Methods Sampling • Rejection • Importance • Hastings-Metropolis • Gibbs Markov Chains • Properties Particle Filtering • Condensation 2 Monte Carlo Methods 3 Markov Chain Monte Carlo History A recent survey places the Metropolis algorithm among the 10 algorithms that have had the greatest influence on the development and practice of science and engineering in the 20th century (Beichl&Sullivan, 2000). The Metropolis algorithm is an instance of a large class of sampling algorithms, known as Markov chain Monte Carlo (MCMC) MCMC plays significant role in statistics, econometrics, physics and computing science. There are several high-dimensional problems, such as computing the volume of a convex body in d dimensions and other high-dim integrals, for which MCMC simulation is the only known general approach for providing a solution within a reasonable time. 4 Markov Chain Monte Carlo History While convalescing from an illness in 1946, Stanislaw Ulam was playing solitaire. It, then, occurred to him to try to compute the chances that a particular solitaire laid out with 52 cards would come out successfully (Eckhard, 1987). Exhaustive combinatorial calculations are very difficult. New, more practical approach: laying out several solitaires at random and then observing and counting the number of successful plays. This idea of selecting a statistical sample to approximate a hard combinatorial problem by a much simpler problem is at the heart of modern Monte Carlo simulation. 5 FERMIAC FERMIAC: The beginning of Monte Carlo Methods Developed in the Early 1930’s The Monte Carlo trolley, or FERMIAC, was an analog computer invented by physicist Enrico Fermi to aid in his studies of neutron transport.
  • THE BEGINNING of the MONTE CARLO METHOD by N

    THE BEGINNING of the MONTE CARLO METHOD by N

    THE BEGINNING of the MONTE CARLO METHOD by N. Metropolis he year was 1945. Two earth- Some Background first electronic computer—the ENIAC— shaking events took place: the at the University of Pennsylvania in Phil- successful test at Alamogordo Most of us have grown so blase about adelphia. Their mentors were Physicist and the building of the first elec- computer developments and capabilities First Class John Mauchly and Brilliant tronicT computer. Their combined impact -even some that are spectacular—that Engineer Presper Eckert. Mauchly, fa- was to modify qualitatively the nature of it is difficult to believe or imagine there miliar with Geiger counters in physics global interactions between Russia and was a time when we suffered the noisy, laboratories, had realized that if electronic the West. No less perturbative were the painstakingly slow, electromechanical de- circuits could count, then they could do changes wrought in all of academic re- vices that chomped away on punched arithmetic and hence solve, inter alia, dif- search and in applied science. On a less cards. Their saving grace was that they ference equations—at almost incredible grand scale these events brought about a continued working around the clock, ex- speeds! When he’d seen a seemingly renascence of a mathematical technique cept for maintenance and occasional re- limitless array of women cranking out known to the old guard as statistical sam- pair (such as removing a dust particle firing tables with desk calculators, he’d pling; in its new surroundings and owing from a relay gap). But these machines been inspired to propose to the Ballistics to its nature, there was no denying its new helped enormously with the routine, rela- Research Laboratory at Aberdeen that an name of the Monte Carlo method.
  • Benchmarking OLTARIS for Deep Space Dose Analyses Using MCNP6

    Benchmarking OLTARIS for Deep Space Dose Analyses Using MCNP6

    Benchmarking OLTARIS For Deep Space Dose Analyses Using MCNP6 By John Daniel Baunach Thesis Submitted to the Faculty of the Graduate School of Vanderbilt University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE in Physics May, 2015 Nashville, Tennessee Approved: Michael G. Stabin, Ph.D. Todd E. Peterson, Ph.D. Copyright © 2015 by John Daniel Baunach All Rights Reserved ii ACKNOWLEDGEMENTS I much first acknowledge the Nuclear Regulatory Commission, for without their Nuclear Education Grant Program I would not have had the chance to follow this research wherever it led. This opportunity for graduate-level research is something I never expected to have, so on behalf of all students who are supported by these funds -- thank you. An immeasurable debt of gratitude is owed to my advisor, Dr. Michael Stabin. His endless patience in teaching the most stubborn and frustrating of students is an inspiration as a person, as a teacher, and as an aspiring health physicist. The lessons I have learned under his tutelage far exceed any mere memorization of facts, and have taught me to expect a high standard in my own research and professional relationships. For this and much more -- thank you. Much debt is also owed to Dr. Robert Singleterry of NASA's Langley Research Center. I could not ask for a more knowledgeable guide for the content of this thesis, and for impressing on me the true challenges of any potential interplanetary spaceflight. For allowing this student to even be in the room with professionals I have looked up to my entire life -- thank you.
  • Stanislaw M. Ulam Papers 1916-1984 Mss.Ms.Coll.54

    Stanislaw M. Ulam Papers 1916-1984 Mss.Ms.Coll.54

    Stanislaw M. Ulam Papers 1916-1984 Mss.Ms.Coll.54 American Philosophical Society 2002 105 South Fifth Street Philadelphia, PA, 19106 215-440-3400 [email protected] Stanislaw M. Ulam Papers 1916-1984 Mss.Ms.Coll.54 Table of Contents Summary Information ................................................................................................................................. 3 Background note ......................................................................................................................................... 4 Scope & content ..........................................................................................................................................7 Administrative Information .........................................................................................................................9 Related Materials ........................................................................................................................................ 9 Indexing Terms ........................................................................................................................................... 9 Bibliography ..............................................................................................................................................11 Collection Inventory ..................................................................................................................................12 Series I. Professional correspondence, 1932-1986................................................................................12
  • History of Monte Carlo

    History of Monte Carlo

    The History of the Monte Carlo Methods Hamed Jafari Research Division of Computer Graphics Institute of Visual Computing & Human-Centered Technology TU Wien, Austria The Monte Carlo Integration The History of Monte Carlo Methods 2 1930s: Enrico Fermi Italian Physicist (1901-1954) „Architect of nuclear age“ Creator of world‘s 1st reactor Developed the idea of statistical sampling Left Italy and became a scientist on the Manhattan Project at Los Alamos in New Mexico The History of the Monte Carlo Methods 3 Early 1940s: The idea of the first digital Computer by John Mauchly and Presper Eckert At University of Pennsylvania They Proposed the idea of electronic computer to Ballistic Research Lab in Aberdeen Presper Eckert (center) Lead to “Project PX” The History of the Monte Carlo Methods 4 The ENIAC The first programmable, electronic, general-purpose digital computer Electronic Numerical Integrator and Computer 27 tons, 167 square meters, 18000 vacuum tubes Completed in 1946 – operational from 1947 until 1955 First review tests were firing table problems for Aberdeen vacuum tubes function tables The History of the Monte Carlo Methods 5 John von Neumann Hungarian-American mathematician, computer scientist and physicist. 1903-1957 1940s: Consultant to both Los Alamos and Aberdeen That connection lead to his involvement with ENIAC The History of the Monte Carlo Methods 6 Nicholas Metropolis Greek-American physicist (1915-1999) Was also working at Los Alamos Known for Metropolis-Hasting algorithm and as one of the creators of MANIAC The History
  • ENIAC– the First Electronic Computer, University of Pennsylvania

    ENIAC– the First Electronic Computer, University of Pennsylvania

    Operated by Los Alamos National Security, LLC for the U.S. Department of Energy's NNSA The Monte Carlo Method and MCNP – A Brief Review of Our 40 Year History Presentation to the International Topical Meeting on Industrial Radiation and Radioisotope Measurement Applications Conference Avneet Sood, PhD XCP-3 Group Leader 10 July 2017 Operated by Los Alamos National Security, LLC for the U.S. Department of Energy's NNSA LA-UR-17-25633 Outline • Origins of the Monte Carlo method – Development of electronic computers and Monte Carlo method occur simultaneously – Ulam, Von Neumann, Richtmeyer, Metropolis, Fermi • Growth and usage of Monte Carlo codes – 1950’s, 1960’s, and 1970’s • Monte Carlo becomes mainstream; nuclear criticality and reactor • Emergence of MCNP • MCNP’s history • MCNP’s upcoming future *nearly all references can be found at: https://laws.lanl.gov/vhosts/mcnp.lanl.gov/references.shtml Los Alamos National Laboratory 7/14/2017 | 3 Abstract The Monte Carlo method for radiation particle transport has its origins at LANL dating back to the 1940’s. The creators of these methods were Drs. Stanislaw Ulam, John von Neumann, Robert Richtmyer, and Nicholas Metropolis. Monte Carlo methods for particle transport have been driving computational developments since the beginning of modern computers; this continues today. In the 1950’s and 1960’s, these new methods were organized into a series of special-purpose Monte Carlo codes, including MCS, MCN, MCP, and MCG. These codes were able to transport neutrons and photons for specialized LANL applications. In 1977, these separate codes were combined to create the first generalized Monte Carlo radiation particle transport code, MCNP.
  • Le Prime Macchine Da Calcolo Italiane E I Loro Inventori1

    Le Prime Macchine Da Calcolo Italiane E I Loro Inventori1

    Le prime macchine da calcolo italiane e i loro inventori1 Silvio Hénin Introduzione Se si tien conto di tutti i tipi di strumenti matematici inventati, fin dalle tavole aritmetiche sumere (ca. 2500 a.C.) e dagli abachi cinesi e greci (V – IV secolo a.C.) si può dedurre che dispositivi per alleviare la noiosa fatica del calcolo numerico apparvero molto presto nella storia delle civiltà [Campbell-Kelly 2003, Needham 1959, Schärlig 2001]. Le tavole delle moltiplicazioni e gli abachi erano più che sufficienti per il limitato fabbisogno di calcolo di mercanti, contabili e amministratori, almeno fino al XVII secolo, quando la rivoluzione scientifica e la crescita delle burocrazie governative richiesero calcoli più complessi e una maggior precisione, tanto da indurre il filosofo e matematico Gottfried Leibniz a esclamare “è indegno che uomini eccellenti debbano sprecare il loro tempo nel calcolo, quando un qualunque contadino potrebbe farlo altrettanto bene con l’aiuto di una macchina” e a inventare un calcolatore meccanico nel 1672, il terzo della storia dopo quelli di Schickard (1623) e di Pascal (1645). Nei successivi due secoli, numerosi scienziati e semplici artigiani si dedicarono all’invenzione e alla costruzione di macchine per il calcolo aritmetico, trasformando lentamente questi congegni da ‘inutili giocattoli […] destinati al fallimento’ [Williams 1990] ‒ seppure ingegnosi ed esteticamente attraenti ‒ in oggetti utili, producibili in massa e venduti in gran numero. Fu l’avvento del capitalismo industriale alla fine del XIX secolo, che causò una ‘crisi del controllo’ [Beninger 1986] che a sua volta generò un mercato ricettivo per macchine calcolatrici veloci e precise, necessarie alla elaborazione di dati numerici e statistici.
  • Monte Carlo Methods: Early History and the Basics

    Monte Carlo Methods: Early History and the Basics

    MCMs: Early History and The Basics Monte Carlo Methods: Early History and The Basics Prof. Michael Mascagni Department of Computer Science Department of Mathematics Department of Scientific Computing Florida State University, Tallahassee, FL 32306 USA E-mail: [email protected] or [email protected] URL: http://www.cs.fsu.edu/∼mascagni MCMs: Early History and The Basics Outline of the Talk Early History of Probability Theory and Monte Carlo Methods Early History of Probability Theory The Stars Align at Los Alamos The Problems The People The Technology Monte Carlo Methods The Birth General Concepts of the Monte Carlo Method References MCMs: Early History and The Basics Early History of Probability Theory and Monte Carlo Methods Early History of Probability Theory Early History of Probability Theory I Probability was first used to understand games of chance 1. Antoine Gombaud, chevalier de Méré, a French nobleman called on Blaise Pascal and Pierre de Fermat were called on to resolve a dispute 2. Correspondence between Pascal and Fermat led to Huygens writing a text on “Probability" 3. Jacob Bernoulli, Abraham de Moivre, and Pierre-Simon, marquis de Laplace, led development of modern “Probability" 4. 1812: Laplace, Théorie Analytique des Probabilités MCMs: Early History and The Basics Early History of Probability Theory and Monte Carlo Methods Early History of Probability Theory Early History of Monte Carlo: Before Los Alamos I Buffon Needle Problem: Early Monte Carlo (experimental mathematics) 1. Problem was first stated in 1777 by Georges-Louis Leclerc, comte de Buffon 2. Involves dropping a needle on a lined surface and can be used to estimate 3.
  • Computer Meets Theoretical Physics the New Frontier of Molecular Simulation

    Computer Meets Theoretical Physics the New Frontier of Molecular Simulation

    THE FRONTIERS COLLECTION Giovanni Battimelli Giovanni Ciccotti Pietro Greco COM PU T ER MEETS THEORETICAL PHYSICS The New Frontier of Molecular Simulation THE FRONTIERS COLLECTION Series Editors Avshalom C. Elitzur, Iyar, Israel Institute of Advanced Research, Rehovot, Israel Zeeya Merali, Foundational Questions Institute, Decatur, GA, USA Thanu Padmanabhan, Inter-University Centre for Astronomy and Astrophysics (IUCAA), Pune, India Maximilian Schlosshauer, Department of Physics, University of Portland, Portland, OR, USA Mark P. Silverman, Department of Physics, Trinity College, Hartford, CT, USA Jack A. Tuszynski, Department of Physics, University of Alberta, Edmonton, AB, Canada Rüdiger Vaas, Redaktion Astronomie, Physik, bild der wissenschaft, Leinfelden-Echterdingen, Germany The books in this collection are devoted to challenging and open problems at the forefront of modern science and scholarship, including related philosophical debates. In contrast to typical research monographs, however, they strive to present their topics in a manner accessible also to scientifically literate non-specialists wishing to gain insight into the deeper implications and fascinating questions involved. Taken as a whole, the series reflects the need for a fundamental and interdisciplinary approach to modern science and research. Furthermore, it is intended to encourage active academics in all fields to ponder over important and perhaps controversial issues beyond their own speciality. Extending from quantum physics and relativity to entropy, consciousness,