ESPResSo 4.0 released

The ESPResSo team is proud to announce that we have just released version 4.0!

Download ESPResSo 4.0

ESPResSo 4.0 is the first release of ESPResSo with a Python interface. I.e., the language in which simulation scripts are written is now Python, and support for the TCL-interface has been dropped. We recommend all users switch to ESPResSo 4.0, there will be no more fixes to TCL-based versions.

Transitioning from TCL-based versions

Please note that, starting from ESPResSo 4.0, enabling additional features in myconfig.hpp no longer changes the behavior of a simulation scripts. Features such as particle rotation and virtual sites have to be activated explicitly in the simulation scripts.

New features

  • Reaction ensemble methods for the modelling of chemical reactions, including constant pH and Wang-Landau schemes.
  • Polarisable particles via the Drude-oscillator scheme.
  • Steepest-descent energy minimization scheme.
  • Methods for active particles (swimmer_reaction, dipolar swimmers).
  • GPU-accelerated P3M method for electrostatic interactions.
  • GPU-accelerated direct summation and Barnes-Hutt schemes for dipolar interactions under open boundary conditions.
  • Support for the electrostatic and magnetostatics method provided by the ScaFaCoS library.
    Immersed boundary method for soft immersed objects in an LB fluid.
  • Virtual sites acting as inertialess tracers in an LB fluid usable as part of the immersed boundary method or independently.
  • New bonded quartic and Coulomb interactions.
  • Possibility to tune the skin parameter.
  • Support for saving to the h5md file format for molecular dynamics data.
  • Connection to the MD Analysis Python package.
  • A rotate_system command was added, which rotates the particle positions around the system’s center of mass.

User-visible changes

  • When the ROTATION feature is compiled in, the rotation of a particle around each of its Cartesian axes in the body-fixed frame can be enabled.
  • The rotational degrees of freedom of virtual sites is now integrated. Virtual sites no longer use the quaternion attributes to store their relative orientation.
  • Alternatively, the orientation of the virtual site can be specified relative to that of the particle it is derived from.
  • The random number generator was changed to Mersenne Twister as implemented by C++11.
  • Espresso now depends on the Boost library including Boost-Mpi and a C++11 capable compiler
  • The build system is now based on CMake.

Removed features

The following features will not be part of the ESPResSo 4.0 release, but support will be re-added in a subsequent release:
  • Generalized hybrid Monte Carlo thermostat.
  • Virtual sites located at the center of mass of a group of particles.
  • Metadynamics, umbrella sampling and parallel tempering.
  • Non-equilibrium molecular dynamics (shear boundary conditions implemented by moving slabs of particles at the boundaries).
  • The memd/maggs electrostatic solver.
  • The com_force feature.
  • The Shan Chen lattice Boltzmann extension.
  • Lees-Edwards (shear) boundary conditions.
  • The chain rdf analysis.
The following functionality is removed permanently:
  • Blockfile support, which is replaced by support for the h5md format as well as checkpointing based on the pickle Python module.
  • Some of the convenience functions for particle creation such as salt and counterions, as these can be replaced by a very few lines of Python in the simulation script.
  • The lj_angle and angledist (not-quite)-pair potentials.
  • The following analysis methods:
    diffusion_profile, cwvac, p_inst, …_mol, cel_gpb, dipmom_normal, MSD, angularmomentum, cluster_size_dist, mol, lipid_orient_order, get_lipid_orients, current, wall_stuff, necklace, bilayer, modes2d.
  • The uwerr error estimator

ESPResSo 4.0 in action

The following python script demonstrates how a simple particle simulation with Lennard-Jones interaction and Langevin thermostat is set up and visualized in ESPResSo .


import espressomd
from espressomd import visualization
import numpy as np

# Check that the Lennard-jones feature is compiled in
required_features = ['LENNARD_JONES']
espressomd.assert_features(required_features)

# System setup
system = espressomd.System(box_l = [10, 10, 10], time_step = 0.0001)

# Use a random seed
system.set_random_state_PRNG()

# Add particles at random positions
for i in range(100):
   system.part.add(pos = np.random.random(3) * system.box_l, type = 0)

# Setup the Lennard-Jones interaction between particles of type 0
system.non_bonded_inter[0, 0].lennard_jones.set_params(
   epsilon=1.0,
   sigma=1.0,
   cutoff=2**(1./6.),
   shift='auto')

# Steepest descent energy minimization for warmup
Etot = system.analysis.energy()['total']
print('Energy before minimization: {} '.format(Etot))
system.minimize_energy.init(
   f_max=1.0,
   gamma=10.0,
   max_steps=100,
   max_displacement=0.01)
system.minimize_energy.minimize()
Etot = system.analysis.energy()['total']
print('Energy after minimization: {} '.format(Etot))

# Activate the Langevin thermostat
system.thermostat.set_langevin(kT=1.0, gamma=1.0)

# Initialize the OpenGL visualizer and start the simulation loop
print('Starting...')
visualizer = visualization.openGLLive(system, background_color = [1, 1, 1])
visualizer.run(1)

Invitation to the ESPResSo Summer School 2018

Our cecam school Particle-based Simulations for Hard and Soft Matter will take place on October 8-12, 2018 at the Institute for Computational Physics, University of Stuttgart.

The school focuses on the introduction of particle based coarse-grained Molecular Dynamics simulation techniques for Hard and Soft Matter systems with the freely available software packages ESPResSo and IMD.  It is intended to introduce the participants into specific methods like Lattice-Boltzmann, Smooth Particle Hydrodynamics, advanced electrostatic algorithms, force field development, finding optimal transition states, and coarse-grained model building, as well as visualization methods with the program MegaMol. This visualization middleware focuses on GPU-accelerated visualization of particle-based data with the scope of interactive exploratory analysis of molecular dynamics simulations.

In addition to lectures, hands-on-tutorials will help to illustrate the usage of all packages as well as to understand the discussed topics. The school will cover examples from soft matter as well as from solid state physics and engineering applications. The MegaMol framework includes also scenarios from thermodynamics and structural biology. It is planned that the participants will present their research by posters. Moreover, invited speakers will give talks on topics like scientific programming with Python, the methodology of multiscale methods, and advanced scientific visualization. Therefore the main aim of the school is given by a broad introduction into modern simulation techniques for Soft and Hard Matter and related topics in terms of practical tutorials, discussion rounds as well as lectures.

For details and for registration, please see: 

https://www.cecam.org/workshop-1605.html

More information about the software packages can be found at:

www.espressomd.org (ESPResSo)

http://imd.itap.physik.uni-stuttgart.de/ (IMD)

https://megamol.org/ (MegaMol)

 

ESPResSo Summer School 2017 Slides

Here are the slides of the ESPResSo summer school 2017:

Online Documentation

We are switching from the document based User Guide  to a Online Documentation with examples, development instructions and cross references to the python interface code.

Soft Matter, Issue 18, 2017

The Cover of Soft Matter, Issue 18, 2017 shows a

polyelectrolyte gel simulated with ESPResSo .

Find details in the article:

“A self-consistent mean-field model for polyelectrolyte gels”
by Oleg Rud, Tobias Richter, Oleg Borisov, Christian Holm
and Peter Košovan

Recent publications citing ESPResSo

Check the Publications section for highlight publications and the complete Scopus RSS feed.