Chris Lorenz

Bio: Chris Lorenz is an academic researcher from King's College London. The author has contributed to research in topics: Dissipative system & Pentamidine. The author has an hindex of 9, co-authored 13 publications receiving 1675 citations.

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

Abstract: The authors present a very general integration scheme for simulating realistic dissipative systems with non-Markovian dynamics and colored noise. The key is a simple, efficient algorithm for the memory kernel.

Modulation of dipalmitoylphosphatidylcholine monolayers by dimethyl sulfoxide.

Abstract: The action of the penetration-enhancing agent, dimethyl sulfoxide (DMSO), on phospholipid monolayers was investigated at the air-water interface using a combination of experimental techniques and molecular dynamics simulations. Brewster angle microscopy revealed that DPPC monolayers remained laterally homogeneous at subphase concentrations up to a mole fraction of 0.1 DMSO. Neutron reflectometry of the monolayers in combination with isotopic substitution enabled the determination of solvent profiles as a function of distance perpendicular to the interface for the different DMSO subphase concentrations. These experimental results were compared to those obtained from molecular dynamic (MD) simulations of the corresponding monolayer systems. There was excellent agreement found between the MD-derived reflectivity curves and the measured data for all of the H/D contrast variations investigated. The MD provide a detailed description of the distribution of water and DMSO molecules around the phosphatidylcholine headgroup, and how this distribution changes with increasing DMSO concentrations. Significantly, the measurements and simulations that are reported here support the hypothesis that DMSO acts by dehydrating the phosphatidylcholine headgroup, and as such provide the first direct evidence that it does so primarily by displacing water molecules bound to the choline group.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

Abstract: The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Fast parallel algorithms for short-range molecular dynamics

Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Biomolecular Simulation: A Computational Microscope for Molecular Biology

Abstract: Molecular dynamics simulations capture the behavior of biological macromolecules in full atomic detail, but their computational demands, combined with the challenge of appropriately modeling the relevant physics, have historically restricted their length and accuracy. Dramatic recent improvements in achievable simulation speed and the underlying physical models have enabled atomic-level simulations on timescales as long as milliseconds that capture key biochemical processes such as protein folding, drug binding, membrane transport, and the conformational changes critical to protein function. Such simulation may serve as a computational microscope, revealing biomolecular mechanisms at spatial and temporal scales that are difficult to observe experimentally. We describe the rapidly evolving state of the art for atomic-level biomolecular simulation, illustrate the types of biological discoveries that can now be made through simulation, and discuss challenges motivating continued innovation in this field.

Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

Abstract: GROMACS is a widely used package for biomolecular simulation, and over the last two decades it has evolved from small-scale efficiency to advanced heterogeneous acceleration and multi-level parallelism targeting some of the largest supercomputers in the world. Here, we describe some of the ways we have been able to realize this through the use of parallelization on all levels, combined with a constant focus on absolute performance. Release 4.6 of GROMACS uses SIMD acceleration on a wide range of architectures, GPU offloading acceleration, and both OpenMP and MPI parallelism within and between nodes, respectively. The recent work on acceleration made it necessary to revisit the fundamental algorithms of molecular simulation, including the concept of neighborsearching, and we discuss the present and future challenges we see for exascale simulation - in particular a very fine-grained task parallelism. We also discuss the software management, code peer review and continuous integration testing required for a project of this complexity.