Source author record

Berk Hess

Berk Hess appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.comp-ph physics.flu-dyn Computational Engineering, Finance, and Science physics.chem-ph

Catalog footprint

What is connected

5works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Nanoscale sheared droplet: Volume-of-Fluid, phase-field and no-slip molecular dynamics

The motion of the three-phase contact line between two immiscible fluids and a solid surface arises in a variety of wetting phenomena and technological applications. One challenge in continuum theory is the effective representation of molecular phenomena close to the contact line. Here, we characterize the molecular processes of the moving contact line to assess the accuracy of two different continuum two-phase models. Specifically, molecular dynamics (MD) simulations of a two-dimensional droplet between two moving plates are used to create reference data for different capillary numbers and contact angles. We use a simple-point-charge/extended (SPC/E) water model with particle-mesh Ewald electrostatics treatment. This model provides a very small slip and a more realistic representation of the molecular physics than Lennards-Jones models. The Cahn-Hilliard phase-field model and the Volume-of-Fluid model are calibrated against the drop displacement from MD reference data. It is demonstrated that the calibrated continuum models can accurately capture droplet displacement and droplet breakup for different capillary numbers and contact angles. However, we also observe differences between continuum and atomistic simulations in describing the transient and unsteady droplet behavior, in particular, close to dynamical wetting transitions. The molecular dynamics of the sheared droplet provide insight of the line friction experienced by the advancing and receding contact lines and evidence of large-scale temporal "stick-slip" like oscillations. The presented results will serve as a stepping stone towards developing accurate continuum models for nanoscale hydrodynamics.

preprint2020arXiv

Electrowetting diminishes contact line friction in molecular wetting

We use large-scale molecular dynamics to study dynamics at the three-phase contact line in electrowetting of water and electrolytes on no-slip substrates. Under the applied electrostatic potential the line friction at the contact line is diminished. The effect is consistent for droplets of different sizes as well as for both pure water and electrolyte solution droplets. We analyze the electric field at the contact line to show how it assists ions and dipolar molecules to advance the contact line. Without an electric field, the interaction between a substrate and a liquid has a very short range, mostly affecting the bottom, immobilized layer of liquid molecules which leads to high friction since mobile molecules are not pulled towards the surface. In electrowetting, the electric field attracts charged and polar molecules over a longer range which diminishes the friction.

preprint2015arXiv

Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

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.

preprint2014arXiv

Accelerated weight histogram method for exploring free energy landscapes

Calculating free energies is an important and notoriously difficult task for molecular simulations. The rapid increase in computational power has made it possible to probe increasingly complex systems, yet extracting accurate free energies from these simulations remains a major challenge. Fully exploring the free energy landscape of, say, a biological macromolecule typically requires sampling large conformational changes and slow transitions. Often, the only feasible way to study such a system is to simulate it using an enhanced sampling method. The accelerated weight histogram (AWH) method is a new, efficient extended ensemble sampling technique which adaptively biases the simulation to promote exploration of the free energy landscape. The AWH method uses a probability weight histogram which allows for efficient free energy updates and results in an easy discretization procedure. A major advantage of the method is its general formulation, making it a powerful platform for developing further extensions and analyzing its relation to already existing methods. Here, we demonstrate its efficiency and general applicability by calculating the potential of mean force along a reaction coordinate for both a single dimension and multiple dimensions. We make use of a non-uniform, free energy dependent target distribution in reaction coordinate space so that computational efforts are not wasted on physically irrelevant regions. We present numerical results for molecular dynamics simulations of lithium acetate in solution and chignolin, a 10-residue long peptide that folds into a $β$-hairpin. We further present practical guidelines for setting up and running an AWH simulation.

preprint2013arXiv

A flexible algorithm for calculating pair interactions on SIMD architectures

Calculating interactions or correlations between pairs of particles is typically the most time-consuming task in particle simulation or correlation analysis. Straightforward implementations using a double loop over particle pairs have traditionally worked well, especially since compilers usually do a good job of unrolling the inner loop. In order to reach high performance on modern CPU and accelerator architectures, single-instruction multiple-data (SIMD) parallelization has become essential. Avoiding memory bottlenecks is also increasingly important and requires reducing the ratio of memory to arithmetic operations. Moreover, when pairs only interact within a certain cut-off distance, good SIMD utilization can only be achieved by reordering input and output data, which quickly becomes a limiting factor. Here we present an algorithm for SIMD parallelization based on grouping a fixed number of particles, e.g. 2, 4, or 8, into spatial clusters. Calculating all interactions between particles in a pair of such clusters improves data reuse compared to the traditional scheme and results in a more efficient SIMD parallelization. Adjusting the cluster size allows the algorithm to map to SIMD units of various widths. This flexibility not only enables fast and efficient implementation on current CPUs and accelerator architectures like GPUs or Intel MIC, but it also makes the algorithm future-proof. We present the algorithm with an application to molecular dynamics simulations, where we can also make use of the effective buffering the method introduces.