CCP5 Summer School 2021

Molecular simulation methods


11th-22nd July 2021, online

registration closed

About


Emails have been sent please check your spam/junk folder if you see nothing in your inbox

Due to the current uncertainty over the COVID situation the school will be taking place online.

Organised by CCP5 and sponsored by CECAM, the School is intended for newcomers to the science of molecular simulation and will provide a comprehensive introduction to the theoretical background as well as practical sessions on computational methods and research seminars to illustrate the versatility of simulation in modern research. There will also be opportunities for participants to present their own research.

The Summer School starts with a two-day programming course, where students can opt to take either Python or modern Fortran. After this preparation, the first five days of the main School will cover the basics of molecular simulation, and the remaining three days will be devoted to more advanced courses with options in mesoscale, ab initio, and biomolecular simulation. Course notes will be provided in electronic format. In addition to the lectures, there will be extensive practical sessions in which students will undertake computational exercises to reinforce and further explore the material.

The school will take place between 11th and 22nd of July 2021 in online.

A fee of £50 to cover part of the expenses will be charged to successful applicants. The school has 50 to 100 places available. Successful applicants from last year’s delayed summer school will have an automatic reserve for this year’s online school and they will be receiving an email about this.

Please note The school can be recognized towards your doctoral training in UK, also upon request we can provide a letter for ECTS credits for your school.

Key dates

  • Application deadline: 15th of April 2021
  • Acceptance decision: 1st of May 2021
  • Fee payments: 1st of June 2021

Organising Committee

  • Dr Colin Freeman, University of Sheffield
  • Prof Neil Allan, University of Bristol
  • Dr Mark Miller, University of Durham
  • Dr Alin Elena, STFC Daresbury Laboratory

Sponsors

Code of Conduct

We value the participation of everyone and want to ensure that everyone has an enjoyable and fulfilling experience, both professionally and personally. Accordingly, all participants of the CCP5 Summer School are expected to always show respect and courtesy to others. The CCP5 and its partners strive to maintain inclusivity in all of our activities. All participants (staff and students) are entitled to a harassment-free experience, regardless of gender identity and expression, sexual orientation, disability, physical appearance, body size, race, age, and/or religion. Harassment in any form is not acceptable for any of us.

We respectfully ask all attendees of the CCP5 Summer School to kindly conform to the following Code of Conduct:

  • Treat all individuals with courtesy and respect.
  • Be kind to others and do not insult or put down other members.
  • Behave professionally. Remember that harassment and sexist, racist, or exclusionary jokes are not appropriate.
  • Harassment includes, but is not limited to, offensive verbal comments related to gender, sexual orientation, disability, physical appearance, body size, race, religion, sexual images in public spaces, deliberate intimidation, stalking, following, harassing photography or recording, sustained disruption of discussions, and unwelcome sexual attention.
  • Participants asked to stop any harassing behaviour are expected to comply immediately.
  • Contribute to communications with a constructive, positive approach.
  • Be mindful of talking over others during presentations and discussion and be willing to hear out the ideas of others.
  • All communication should be appropriate for a professional audience, and be considerate of people from different cultural backgrounds. Sexual language and imagery are not appropriate at any time.
  • Challenge behaviour, action and words that do not support the promotion of equality and diversity.
  • Arrive at the events punctually where possible.
  • Show consideration for the welfare of your friends and peers and, if appropriate, provide advice on seeking help.
  • Seek help for yourself when you need it.

please report any issues to alin-marin.elena@stfc.ac.uk

Registration


Please note the registration is closed!

Lectures


Programming Courses

  • Introduction to Modern Fortran (6 lectures, 4 practical sessions)
  • Introduction to Python

Basic Courses

  • An Overview of Molecular Simulation
  • Statistical Mechanics (2 lectures)
  • Molecular Dynamics (3 lectures)
  • Monte Carlo Methods (3 lectures)
  • Free Energy Methods (3 lectures)
  • Optimisation Methods
  • Introduction to Force Fields
  • Long timescale methods
  • Advanced Free Energy methods
  • Practicals (10 sessions over 5 afternoons)

Advanced Courses

Lecturers

First principles simulations

Mesoscale methods

Simulation of organic and biomolecules

Machine Learning and Interatomic Potentials

  • Prof Gábor Csányi (University of Cambridge)
  • David Kovacs (University of Cambridge)
  • Tamas Stenczel (University of Cambridge)

Programming

Timetable


Date Activity Location
July 12 Day 1
9:00 - 9:10 Welcome
9:10 - 10:00 Fortran I/Python I AE/MB
10:00 - 10:10 break
10:10 - 10:30 Icebreaker I
10:30 - 11:30 Fortran II/Python II AE/MB
11:30 - 11:45 break
11:45 - 12:45 Practical
12:45 - 13:45 Lunch
13:45 - 14:45 Fortran III/Python III AE/MB
14:45 - 15:00 break
15:00 - 16:00 Practical
16:00 - 16:15 break
16:15 - 17:15 Student talks
July 13 Day 2
9:00 - 10:00 Fortran IV/Python IV AE/MB
10:00 - 10:10 break
10:10 - 10:30 Icebreaker II
10:30 - 11:30 Fortran V/Python V AE/MB
11:30 - 11:45 Break
11:45 - 12:45 Practical
12:45 - 13:45 Lunch
13:45 - 14:45 Fortran VI/Python VI AE/MB
14:45 - 15:00 break
15:00 - 16:00 Practical
16:00 - 16:15 break
16:15 - 17:15 Dr Chantal Valeriani, Complutense University of Madrid
Trapping active particles up to the limiting case: bacteria enclosed in a biofilm
July 14 Day 3
9:00 - 10:00 Overview of molecular simulations - PC
10:00 - 10:15 break
10:15 - 11:15 Statistical Mechanics 1 - MB
11:15 - 11:30 break
11:30 - 12:30 Statistical Mechanics 2 - MB
12:30 - 13:30 Lunch
13:30 - 14:30 Practical - Stat Mech Problems
14:30 - 14:40 Break
14:40 - 15:10 Flash posters
15:10 - 15:20 break
15:20 - 16:20 Poster Session I
16:20 - 16:30 break
16:30 - 17:30 Professor David Wales, University of Cambridge
Energy landscapes: from molecules and nanodevices to machine learning
July 15 Day 4
9:00 - 10:00 Introduction to force fields - PC
10:00 - 10:15 break
10:15 - 11:15 Monte Carlo 1 - NA
11:15 - 11:30 Virtual Refreshments
11:30 - 12:30 Practical - MC integration
12:30 - 13:30 Lunch
13:30 - 14:30 Monte Carlo 2 - NA
14:30 - 14:40 Virtual Refreshments
14:40 - 15:40 Practicals - Intro to MC
15:40 - 16:20 flash session
16:20 - 16:30 break
16:30 - 17:30 Poster session 2
July 16 Day 5
9:00 - 10:00 Molecular Dynamics 1 - CF
10:00 - 10:15 break
10:15 - 11:15 Molecular Dynamics 2 - CF
11:15 - 11:30 Virtual Refreshments
11:30 - 12:30 Practical - Intro to MD
12:30 - 13:30 Lunch
13:30 - 14:30 Monte Carlo 3 - MA
14:30 - 14:50 Virtual Refreshments
14:50 - 16:00 Practicals - Phase Equilibria
16:00 - 16:15 break
16:15 - 17:15 Prof David Quigley, University of Warwick
Ice interfaces, metastable phases and nucleation – testing hypotheses with free energy calculations
July 17 Day 6
9:00 - 10:00 Molecular Dynamics 3 - MA
10:00 - 10:15 break
10:15 - 11:15 Long timescale methods - JH
11:15 - 11:30 Virtual Refreshments
11:30 - 12:30 Practicals - Stability + accur MD
12:30 - 13:30 Lunch
13:30 - 14:30 Optimisation methods - JH
14:30 - 14:50 Virtual Refreshments
14:50 - 16:00 Practicals - Forcefields optimisation
16:00 - 16:15 break
16:15 - 17:15 free social
** July 18 Free day ** Day 7
July 19 Day 8
9:00 - 10:00 Free energy methods 1 - JA
10;00 - 10:15 break
10:15 - 11:15 Free energy methods 2 - JA
11:15 - 11:30 Virtual Refreshments
11:30 - 12:30 practical Shake
12:30 - 13:30 Lunch
13:30 - 14:30 Advanced Free Energy - LBP
14:30 - 14:50 Virtual Refreshments
14:50 - 16:00 practical Thermostats
16:00 - 16:15 break
16:15 - 17:15 Dr Marina Filip, University of Oxford
Optical excitations in metal-halide perovskite semiconductors from first principles computational modeling
July 20 Day 9
9:00 - 12:30 Advanced Seminar 1-3
12:30 - 13:30 Lunch
13:30 - 16:00 Practicals
16:00 - 16:15 break
16:15 - 17:15 Elin Barrett, Unilever UK
Computational Methods for Next-Generation Risk Assessment of Consumer Goods
July 21 Day 10
9:00 - 12:30 Advanced Seminar 4-6
12:30 - 13:30 Lunch
13:30 - 16:00 Practicals
16:00 - 16:15 break
16:15 - 17:15 Professor Peter Coveney, UCL
Reliability and reproducibility in computational science
July 22 Day 11
9:00 - 12:30 Advanced Seminar 7
12:30 - 13:30 Lunch
13:30 - 15:30 Practicals

Advanced Seminars may be structured different depending on the lecturers.

Advanced: Machine Learning

Date Activity Location
July 20 Day 9
9:30 - 10:00 setup, verify that Deepnote works (Tamas Stenczel)
10:00 - 11:00 Lecture: ML potentials and GAP for materials (Gábor Csányi)
11:00 - 11:15 Introduction to practical on materials (Tamas Stenczel)
11:15 - 12:30 practical session 1 (online help provided)
12:30 - 13:30 Lunch
13:30 - 15:00 practical session 2 (online help provided)
15:00 - 15:30 Demonstration of further tools and examples (Tamas Stenczel)
15:30 - 16:00 Miguel Caro (Research Seminar) - Large-scale simulation of carbon materials with GAP
16:00 - 16:15 break
16:15 - 17:15 Elin Barrett, Unilever UK
Computational Methods for Next-Generation Risk Assessment of Consumer Goods
July 21 Day 10
9:00 - 10:00 Lecture: ML potentials for molecules (Gábor Csányi)
10:00 - 10:30 Introduction to practical on molecules (Dávid Kovács)
10:30 - 12:30 practical session 1 (online help provided)
12:30 - 13:30 Lunch
13:30 - 15:00 practical session 2 (online help provided)
15:00 - 15:30 Seminar on molecular fitting with other tools (Dǽvid Kovács)
15:30 - 16:00 Noam Bernstein (Research Seminar) - Machine Learning Interatomic Potentials with and without (much) human labor
16:00 - 16:15 break
16:15 - 17:15 Professor Peter Coveney, UCL
Reliability and reproducibility in computational science
July 22 Day 11
9:00 - 09:45 Volker Deringer - GAP modelling for inorganic materials chemistry
9:45 - 10:30 Johannes Margraf and Simon Wengert - Gap for molecular properties
10:30 - 10:45 break
10:45 - 11:30 Ioan-Bogdan Magdau - Machine Learning Molecular Liquids: Challenges and Solutions
11:30 - 12:15 Albert Bartók-Pártay - Gaussian Approximation Potentials: implementation notes

Advanced: Mesoscale

Date Activity Location
July 20 Day 9
9:00 - 9:30 Introduction to Mesoscale Simulation (Michael Seaton)
9:30 - 10:00 Dissipative Particle Dynamics (Michael Seaton)
10:00 - 10:15 break
10:15 - 11:15 Dissipative Particle Dynamics
11:15 - 11:30 break
11:30 - 12:00 Dissipative Particle Dynamics
12:00 - 12:30 Mesoscale Practical Exercises: Guide to DPD exercises (Michael Seaton)
12:30 - 13:30 Lunch
13:30 - 16:00 Practicals: DPD tutorial exercises (1-4)
16:00 - 16:15 break
16:15 - 17:15 Elin Barrett, Unilever UK
Computational Methods for Next-Generation Risk Assessment of Consumer Goods
July 21 Day 10
9:00 - 10:00 Lattice Boltzmann (Ian Halliday)
10:00 - 10:15 break
10:15 - 11:15 Lattice Boltzmann
11:15 - 11:30 break
11:30 - 12:00 Lattice Boltzmann
12:00 - 12:30 Notes on LBE tutorial exercises (Ian Halliday)
12:30 - 13:30 Lunch
13:30 - 16:00 Practicals: LBE tutorial exercise 1
16:00 - 16:15 break
16:15 - 17:15 Professor Peter Coveney, UCL
Reliability and reproducibility in computational science
July 22 Day 11
9:00 - 10:00 (Multicomponent) Lattice Boltzmann (Ian Halliday)
10:00 - 10:15 break
10:15 - 12:30 Practicals: LBE tutorial exercise 2
12:30 - 13:30 Lunch
13:30 - 15:30 Practicals

Advanced: First Principles

Date Activity Location
July 20 Day 9
9:00 - 10:00 Intro to DFT (SC)
10:00 - 10:15 break
10:15 - 11:15 Intro to Crystal (BS)
11:15 - 11:30 break
11:30 - 12:30 Intro to Castep (KR)
12:30 - 13:30 Lunch
13:30 - 16:00 Practicals
16:00 - 16:15 break
16:15 - 17:15 Elin Barrett, Unilever UK
Computational Methods for Next-Generation Risk Assessment of Consumer Goods
July 21 Day 10
9:00 - 10:00 Quantum Chemistry (BS)
10:00 - 10:15 break
10:15 - 11:15 Reciprocal space (SC)
11:15 - 11:30 break
11:30 - 12:30 Plane wave machinery (SC)
12:30 - 13:30 Lunch
13:30 - 16:00 Practicals
16:00 - 16:15 break
16:15 - 17:15 Professor Peter Coveney, UCL
Reliability and reproducibility in computational science
July 22 Day 11
9:00 - 10:00 Convergence (KR)
10:00 - 10:15 break
10:15 - 12:30 Practicals
12:30 - 13:30 Lunch
13:30 - 15:30 Practicals

Research Seminar Speakers


Campus information


it is online

Student Events


Student Seminar (13+2 minutes)

  • Alexander von Wedelstedt, Entropic Separation Effects in Adsorption Processes, HTW University of Applied Sciences, Dresden, Germany
  • Rebecca Walters, Interactive, Flexible Docking to the SARS-CoV-2 Main Protease in Virtual Reality, University of Bristol, UK
  • Mohammed El ismaili,Activity-modulated phase transition in a two-dimensional mixture of active and passive colloids. Univerisite Lorraine, France
  • Samira Anker, Novel approach to control pharmaceutical nucleation using heterogeneous nucleation, University of Strathclyde, Glasgow, UK

CCP5 Student Poster Prize

will be virtual.

  1. Charlotte E. Breakwell - Atomistic Simulations of Sulfonated Polymer of Intrinsic Microporosity (PIM) Membranes - Imperial College London, UK

  2. Adele Hardie - Markov State Modelling for Allosteric Drug Discovery  - University of Edinburgh, UK

  3. Anna Katarina Antonovic - Rational Design of Modulators of Skeletal Myosin - Queen Mary University London, UK

  4. Federica Lodesani - Driving Crystallization in Highly Viscous Binary Silicate Melts Through Metadynamics Simulations, Univerisy of Modena, Italy

  5. Wojciech Lipinski - Make Less Noise Controlling protein aggregation by compartmentalisation - Radboud University, Netherlands

Contact


For more information do not hesitate to contact Alin M Elena alin-marin.elena@stfc.ac.uk

Overview basic lectures


An Overview of Molecular Simulation

An overview of the current state of molecular simulation with examples of special interest taken from the literature.

Introduction to force fields

Statistical Mechanics 1

In this lecture we will begin with an important question: why bother with statistical thermodynamics? We will progress to basic statistical quantities and concepts such as averages, fluctuations and correlations and how to use them in practice to calculate the physical properties of systems. This will lead us to the determination of the true statistical error for system properties obtained by simulation. We will apply these ideas to commonly calculated properties such as diffusion, radial distribution functions and velocity autocorrelation, while also examining the physical meaning of these properties. We will conclude with a look at distribution functions: how they arise and what they mean.

Statistical Mechanics 2

In the second lecture we shall look at more theoretical aspects of statistical mechanics. Beginning with the Lagrange and Hamiltonian description of classical mechanics we shall progress to the idea of phase space and the concept of a probability distribution function. This will be followed by basic applications (and associated mathematical manipulations) of the distribution function to obtain various physical properties of a system. We will examine the common ensembles (NVE, NVT and NPT) and discuss their application and interrelation. Finally we shall look at time dependence, beginning with the Liouville Equation and its connection with other time dependent equations. We shall conclude with the fluctuation-dissipation theorem.

Monte Carlo 1

Basics: The system. Random sampling. Importance sampling. Detailed balance. Metropolis algorithm in the canonical ensemble. Isothermal-isobaric ensemble. Grand-canonical ensemble. Which ensemble?

Monte Carlo 2

Practicalities: Finite-size effects. Random number generators. Tuning the acceptance rate. Equilibration. Configurational temperature. Ergodicity and free-energy barriers. Measuring ensemble averages. Examples (showing ensemble independence for the Lennard-Jones fluid)

Monte Carlo 3

(Free) Energy Barriers: Quasi non-ergodicity. Vapour-liquid phase transition as an example. Removing the interface by Gibbs ensemble MC. Free-energy barrier in the grand-canonical ensemble. Multicanonical preweighting. Histogram reweighting. Parallel tempering

Molecular Dynamics 1

Molecular dynamics: the basic methodology. Integration algorithms and their derivation. Static properties: thermodynamics and structure. Dynamic properties: correlation functions and collective properties

Molecular dynamics 2

Practical aspects of molecular dynamics - Verlet neighbour list, link cell algorithm. Calculating pressure: the virial theorem and the thermodynamic method. Estimating statistical errors: the blocking method. Symplectic algorithms and the Tuckerman-Berne-Martyna approach.

Molecular dynamics 3

Extended systems: canonical (NVT) and isothermal-isobaric (NPT) ensembles. Rigid Bodies, SHAKE, RATTLE.

Free energy methods 1

Free energy, chemical potential & thermodynamics. Applications. Essential statistical mechanics. Ensemble averages, probability distributions & simulations. Free energy, the challenge. Particle insertion & removal. Energy density distributions. The perturbation method.

Free energy methods 2

Review essential statistical mechanics. Thermodynamic integration. Potential of mean force calculations. Umbrella sampling. Absolute free energies. Free energy of liquids.Free energy of solids.

Optimization Methods

The energy landscape, geometrical optimisation and saddle points. Minimisation methods (steepest descent, conjugate gradient, genetic algorithm). Saddle-points (transition state theory, harmonic theory, nudged elastic band, dimer method).

Long timescale methods

Long timescales simulations - the problems. Transition state theory and kinetic Monte Carlo. Temperature accelerated hyperdynamics. Metadynamics.

Advanced Free Energy Methods

TBD

First-principles simulation


First-principles simulation has grown to become one of the most influential and important techniques for modelling at the atomic level. With nuclei and electrons as the basic ingredients the system is modelled at a deeper level of physics than with atoms and interatomic potentials. By explicitly including the electrons in the model and treating their interactions using quantum-mechanical laws, chemical bonding arises as an emergent phenomenon of the model. All kinds of bonding - ionic, covalent, metallic, hydrogen can be treated using the same method. The price of this accurate Hamiltonian is a computational cost orders of magnitude higher than atomic potential models. Nevertheless it is possible and convenient with modern parallel computers to simulate systems of hundreds of atoms, and perform optimization and molecular dynamics in a variety of ensembles.

In this advanced course I will provide a rapid introduction to the “nuts and bolts” of first-principles simulation. In accordance with the philosophy of the CCP5 Summer School, the aim is to attempt to open up the “black box” and explain the concepts and algorithms used. The presentation will assume a familiarity with wave mechanics at the undergraduate level and Dirac notation.

In the practicals you will be able to try for yourself using an advanced density functional code. You should be capable of running realistic calculations by the end of the course, and aware of the major aspects of setup and testing that are vital ingredients for success. The practicals will consist of a series of guided exercises with the CASTEP and CRYSTAL codes.

Synopsis

An Introduction to First-Principles Simulation

  • Motivation
  • Quantum-Mechanical approaches
  • Density-Functional Theory
  • Excited states: TD-DFT
  • Electronic Structure of Condensed Phases
  • Total-energy calculations
  • Basis sets
  • Plane waves and pseudopotentials
  • How to solve the equations
  • Ab-initio simulations

Practical calculations using first-principles QM: Convergence, convergence, convergence

  • Convergence
  • Structural Calculations
  • Lattice Dynamics
  • Exchange and Correlation Functionals
  • Summary

Further Study

The lecture notes from the CASTEP workshop held in 2007 are available from http://www.castep.org. Links to a number of ab-initio methods and resources are available at http://electronicstructure.org/.

Mesoscale Methods


Mesoscale methods of modelling are capable of tackling larger length and time scales than those available using atomistic methods. By using particles considerably larger than atoms and appropriate choices of interactions between them, these techniques can readily model bulk materials and large structures at the cost of omitting some fine atomic detail. Hydrodynamics start to become more important at these scales: these modelling techniques are thus designed to ensure correct (emergent) fluid behaviour. A mesoscale model can be set up either using a ‘bottom-up’ approach from atomistic models, a ‘top-down’ approach from continuum fluid models, or both.

In this advanced course we will provide an introduction to two mesoscale methods: Dissipative Particle Dynamics (DPD) and the Lattice Boltzmann Equation (LBE) method. We will explain the origins, concepts and algorithms of both methods, as well as their applications, continuing developments and how they can be related to material models at smaller and larger scales (including those covered by the basic lectures).

In the practicals, you will be able to try out DPD and LBE using both simple ‘hackable’ codes and the general-purpose mesoscale modelling package DL_MESO. By the end of the course, you will gain insight into the capabilities of both mesoscale modelling methods. The practicals will consist of a series of guided exercises using the provided codes.

Synopsis

Introduction to the Mesoscale

  • Techniques
  • Physical scales
  • Mesoscale simulation strategies

Dissipative Particle Dynamics (DPD)

  • DPD algorithm
  • Fokker-Planck formulation
  • Application to simple/complex fluids
  • Boundary conditions
  • Thermodynamics and DPD
  • Molecular dynamics and DPD

Lattice Boltzmann Equation (LBE)

  • Classical Boltzmann/Boltzmann Bhatnagar-Gross-Krook (BGK) Equations
  • Lattice Gas Cellular Automata (LGCA)
  • Multiple component or “diphasic” LGCA
  • Lattice Boltzmann Equation method
  • Lattice Boltzmann BGK Equation and kinetic theory
  • LBE for multi-component flow

Simulation of Organic and Bio Molecules


Biomolecular systems can include proteins, DNA, lipids and the small molecules that interact with them. Individual residues, such as amino acids or nucleic acids, combine to form large complex macromolecules. Here we will focus on how simulation tools like those you have learned about in the summer school can be used to study the structure and function of biomolecules.

This advanced course will cover everything from setting up a biomolecular system for simulation to analysing the results. In addition to the standard molecular dynamics, enhanced sampling methods (including metadynamics), free energy methods and multiscale methods will be explained. One lecture will be devoted to nucleic acids. There will be hands-on practical sessions to accompany each of the lectures.

Synopsis

Introduction

  • Biomolecules
  • Molecular Dynamics Software
  • Force Fields for Biomolecules

Set-up / Analysis

  • Errors/problems in PDB files
  • Checking/choosing protonation states
  • Solvation
  • Analysis: assessing convergence and sampling

Enhanced Sampling

  • Replica exchange methods
  • Biased sampling: methods based on modified Hamiltonians
  • Biased sampling: methods based on unmodified Hamiltonians

Nucleic Acids

  • DNA

Multiscale Modelling

  • QM/MM
  • Coarse-graining
  • Code Coupling

Free Energy Methods

  • Ligand Binding
  • FEP
  • Alchemical Perturbations

Machine Learning for Interatomic Potentials


Synopsis

The two-and-a-half-day Machine Learning for Molecular Simulation activities consist of introductory lectures, practicals, and advanced seminars. The first two days are structure similarly, starting with a lecture with a high level overview and some theoretical background, followed by an introduction to the practicals. Students will use Jupyter notebooks on the Deepnote.com cloud provider, and the Gaussian Approximation Potential framework for fitting machine learned interatomic potentials. On the first day, the practical is about hard materials, the example is silicon carbide. The second day again starts with a higher level lecture, this time focussing on fitting molecular potential energy surfaces, and the following practical sessions, again using the GAP framework take the students through fitting a potential for toluene. On the last half-day, a series of advanced research seminars will be given by practitioners at the forefront of the field.

Campus Maps


to be added soon