MSE 485/CSE 485/PHYS 466: Atomic-Scale Simulations

Fall 2019


TR 9:30-10:50am in 3217 Everitt Laboratory

Course content



Students will be able to

  1. explain fundamentals of molecular dynamics: integration algorithms, static and dynamic correlations functions, and their connection to order and transport
  2. explain fundamentals of Monte Carlo and random walks: variance reduction, Metropolis algorithms, kinetic Monte Carlo, heat diffusion, Brownian motion
  3. apply molecular dynamics and Monte Carlo to relevant material problems, such as
    1. phase transitions (melting/freezing, calculating free energies)
    2. polymers (growth and equilibrium structure)
  4. determine proper parameters for a simulation;
  5. determine if a simulation result is reliable based on the underlying physical principles in the simulation;
  6. systematically analyze data to extract meaningful materials quantities.


If you have not passed a prerequisite course or are concerned about your background (especially if you have no prior programming experience), please see the instructor before continuing.


Dallas R. Trinkle (; @DallasTrinkle; 308 MSEB in the west stairwell).

Teaching assistant: Han-Yi Chou (grading;

Teaching approach: Lecture + Hands-on project-based learning

An “active learning” approach where we focus on the course objectives: understanding and applying computational modeling tools to simulate material behavior, determine if a simulation is reliable, and use a systematic approach to computation. You will:

Over the course of the semester, you will implement statistical analysis, molecular dynamics, and Monte Carlo techniques to atomic-scale simulations.



Homework will be assigned through the MSE 485 website and assignments are due 11:59 pm the day posted on the MSE 485 website. Late submissions will be penalized by 50% for each day late, unless excused in advance. Your reports must be submitted electronically via GradeScope. The only format for the report that will be accepted for submission is a single, properly-ordered PDF, in portrait format; your name must be printed legibly on the top of the first page.

These homework assignments require you to write code that is bug free and actually works. This code will be automatically graded and checked using the PrairieLearn web site, which will provide you with immediate feedback, allowing you to fix buggy code before using it to prepare the data needed to write your reports. The TA grades only your report.

The written reports are assigned to practice the communication of scientific concepts in writing. They will be graded based on presentation, neatness, correct use of symbols, quality of drawings and diagrams, and clarity of explanation. Reports should be neat and organized! More details and an example will be provided. Tables and graphical representations of results need to be generated using some software program such as Python, Excel, TecPlot, MatLab, etc., rather than being hand-drawn. Correct interpretation and implementation of the problem and correct final answers are important.


In the next few weeks, we will form teams that balance interests, programming ability, and experience. The team will be given (i) a collective grade for a status report and its presentation in class, (ii) a collective grade for the final report, and (iii) a collective grade for the presentation of the final results. For your status report and final report we will use peer review, which will also be part of your grade.

We expect the project itself to take into account:

Prior to the status report, two slides used for this presentation need to be submitted. It needs to outline a problem (Scientific research) and explain what the team will do to solve it (Algorithm development), according to the criteria given above. Also the final reports and the final presentations need to be submitted electronically. Late submissions will be penalized by 50 % for each day late. If you have any questions about the suitability of your project please get in touch with the instructor.



Letter grades:


To succeed in this class, you will need to

Obtaining help: The main two ways to obtain help are online at Canvas or in person at office hours. Please do not send email directly to TAs or professors for routine help or absences! In cases of emergencies related to assignments (e.g., illness) you should contact your professor at the earliest possible opportunity.

Formal and Informal Accommodations

I am committed to assisting students requiring special accommodations for circumstances that are registered with the DRES Student Services Department. These formal accommodations should be discussed with me as early as possible in the semester or as soon after DRES approval as possible.

If you are not formally registered with DRES and have anxiety, depression, learning disabilities, or other issues that affect your ability to fully participate and learn in this class, you are encouraged to check-in with me so we can determine together the kind of support you need to thrive in this class. Please set up a meeting with me through Canvas.

Inclusion and Diversity

I value all students regardless of their background, race, religion (creed), ethnicity, gender, gender expression, age, country of origin, disability status, marital status, sexual orientation, or military status, etc., and am committed to providing a climate of excellence and inclusiveness within all aspects of the course. If there are aspects of your culture or identity that you would like to share with me as they relate to your success in this class, I am happy to meet to discuss. Likewise, if you have any concerns in this area of facing any special issues or challenges, you are encouraged to discuss the matter with me (set up a meeting via Canvas) with an assurance of full confidentiality (only exception being mandatory reporting of academic integrity / code violation and sexual harassment). Harassment or discrimination of any kind will not be tolerated.

Academic Integrity

You are bound by the University Honor Code in this course. Any violation of the Honor Code will result in disciplinary action.

Students are responsible for producing their own code, reports, and projects. Collaborative interaction is encouraged, but each student must do their own implementation, perform all calculations themselves, and write their own reports. Plagiarism will not be tolerated, and verified incidents will result in all parties receiving a zero and formal academic sanctions. Students are responsible for familiarizing themselves with the definition of and penalties for plagiarism in Section I-401 of the UIUC Student Code. Note that plagiarism includes “copying another student’s paper or working with another person when both submit similar papers without authorization to satisfy an individual assignment.”

Changes to syllabus

May occur as deemed necessary by the professor; they will be announced via Canvas and the updated syllabus posted on the course website.


Suggested readings are from

Lecture day Reading Topic notes
T Aug 27 Orientation, Introduction, Statistics I
R Aug 29 Statistical errors, Frenkel Smit (Chapter 4 and Appendix D), LeSar (Appendix G) Statistics and Errors
T Sep 03 Stat. mech, Frenkel Smit (Chapter 2), LeSar (Appendix D, G) Correlation, Bias; Python and example notebook
R Sep 05 Statistical mechanics
T Sep 10 MD 1, MD 2, Frenkel Smit (Chapter 4), LeSar (Chapter 3) Molecular dynamics, boundary conditions and measurement; forces, time propagation
R Sep 12 Potentials, Frenkel Smit (Appendix F), LeSar (Chapter 5) Molecular dynamics: Code, force fields HW1/pdf due 9/13 upload
T Sep 17 no class
R Sep 19 Molecular dynamics: force fields II
T Sep 24 Molecular dynamics: correlation functions
R Sep 26 Correlations Pair correlation functions HW2/pdf due 9/27 upload
T Oct 01 Dynamics, MD ensembles, Frenkel Smit (Chapter 6, Appendix B, C, E), LeSar (Chapter 6) Thermostats
R Oct 03 Monte Carlo intro Intro to Monte Carlo
T Oct 08 Sampling Random numbers and psuedorandom number generators
R Oct 10 Non-uniform distributions HW3/pdf due 10/11 upload
T Oct 15 Importance Sampling, Random walk, Frenkel Smit (Chapter 3), LeSar (Chapter 7) Variance reduction
R Oct 17 Alternate MC schemes, Frenkel Smit (Chapter 13, 17) Alternate Monte Carlo schemes
T Oct 22 Brownian dynamics
R Oct 24 Kinetic MC, Frenkel Smit (Chapter 8), LeSar (Chapter 9) Kinetic Monte Carlo HW4/pdf due 10/25 upload
T Oct 29 Lattice MC, Advanced MC, LeSar (Chapter 7) MC and Ising model
R Oct 31 Phase transitions and critical slowing down
T Nov 05 Frankel Smit (Chapter 15) Ewald summation
R Nov 07 Free energy, Frenkel Smit (Chapter 7, 11) Free energy techniques HW5/pdf due 11/08 upload
T Nov 12 Constraints in MD + Constraints: polymers
R Nov 14 Quantum mechanics for electrons Project proposal draft due 11/15
T Nov 19 Quantum mechanics for electrons
R Nov 21 Project proposal presentations Peer review due 11/22
T Nov 26 no class (fall break)
R Nov 28 no class (fall break)
T Dec 03 no class (fall MRS)
R Dec 05 no class (fall MRS) HW6 due 12/06 upload
T Dec 10 Review and questions Final report due 12/11
R Dec 12 no class (reading day)
T Dec 17 1:30-4:30pm Final presentations (exam week) slides due 10:00am
R Dec 19 no class Peer review due 12/19

Accessing library resources

The Univ. Illinois library has access to a huge variety of electronic resources; this plus additional online resources will be our references. Many can be accessed from the library’s website, or via the campus VPN. Alternatively, you can take advantage of the library proxy. This is done by appending to the web address; when reloaded, you will be asked for Univ. Illinois authentication, and then will be able to access the resource as if you were on campus. In general, this authentication is required only once per session. So, the website

would become

Alternatively, install the Proxy Bookmarklet which makes it extremely easy to use the proxy. I highly recommend this method.