Researcher profile

Gerhard Hummer

Gerhard Hummer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2022arXiv

Rebinding kinetics from single-molecule force spectroscopy experiments close to equilibrium

Analysis of bond rupture data from single-molecule force spectroscopy experiments commonly relies on the strong assumption that the bond dissociation process is irreversible. However, with increased spatiotemporal resolution of instruments it is now possible to observe multiple unbinding-rebinding events in a single pulling experiment. Here, we augment the theory of force-induced unbinding by explicitly taking into account rebinding kinetics, and provide approximate analytic solutions of the resulting rate equations. Furthermore, we use a short-time expansion of the exact kinetics to construct numerically efficient maximum likelihood estimators for the parameters of the force-dependent unbinding and rebinding rates, which pair well with and complement established methods, such as the analysis of rate maps. We provide an open-source implementation of the theory, evaluated for Bell-like rates, which we apply to synthetic data generated by a Gillespie stochastic simulation algorithm for time-dependent rates.

preprint2022arXiv

Small ionic radii limit time step in Martini 3 molecular dynamics simulations

Among other improvements, the Martini 3 coarse-grained force field provides a more accurate description of the solvation of protein pockets and channels through the consistent use of various bead types and sizes. Here, we show that the representation of Na$^+$ and Cl$^-$ ions as "tiny" (TQ5) beads limits the accessible time step to 25 fs. By contrast, with Martini 2, time steps of 30-40 fs were possible for lipid bilayer systems without proteins. This limitation is relevant for, e.g., phase separating lipid mixtures that require long equilibration times. We derive a quantitative kinetic model of time-integration instabilities in molecular dynamics (MD) as a function of time step, ion concentration and mass, system size, and simulation time. With this model, we demonstrate that ion-water interactions are the main source of instability at physiological conditions, followed closely by ion-ion interactions. We show that increasing the ionic masses makes it possible to use time steps up to 40 fs with minimal impact on static equilibrium properties and on dynamical quantities such as lipid and ion diffusion coefficients. Increasing the size of the bead representing the ions (and thus changing their hydration) also permits longer time steps. The use of larger time steps in Martini 3 simulations results in a more efficient exploration of configuration space. The kinetic model of MD simulation crashes can be used to determine the maximum allowed time step whenever sampling efficiency is critical.

preprint2020arXiv

Systematic errors in diffusion coefficients from long-time molecular dynamics simulations at constant pressure

In molecular dynamics simulations under periodic boundary conditions, particle positions are typically wrapped into a reference box. For diffusion coefficient calculations using the Einstein relation, the particle positions need to be unwrapped. Here, we show that a widely used heuristic unwrapping scheme is not suitable for long simulations at constant pressure. Improper accounting for box-volume fluctuations creates, at long times, unphysical trajectories and, in turn, grossly exaggerated diffusion coefficients. We propose an alternative unwrapping scheme that resolves this issue. At each time step, we add the minimal displacement vector according to periodic boundary conditions for the instantaneous box geometry. Here and in a companion paper [J. Chem. Phys. XXX, YYYYY (2020)], we apply the new unwrapping scheme to extensive molecular dynamics and Brownian dynamics simulation data. We provide practitioners with a formula to assess if and by how much earlier results might have been affected by the widely used heuristic unwrapping scheme.

preprint2012arXiv

Pair diffusion, hydrodynamic interactions, and available volume in dense fluids

We calculate the pair diffusion coefficient D(r) as a function of the distance r between two hard-sphere particles in a dense monodisperse suspension. The distance-dependent pair diffusion coefficient describes the hydrodynamic interactions between particles in a fluid that are central to theories of polymer and colloid dynamics. We determine D(r) from the propagators (Green's functions) of particle pairs obtained from discontinuous molecular dynamics simulations. At distances exceeding 3 molecular diameters, the calculated pair diffusion coefficients are in excellent agreement with predictions from exact macroscopic hydrodynamic theory for large Brownian particles suspended in a solvent bath, as well as the Oseen approximation. However, the asymptotic 1/r distance dependence of D(r) associated with hydrodynamic effects emerges only after the pair distance dynamics has been followed for relatively long times, indicating non-negligible memory effects in the pair diffusion at short times. Deviations of the calculated D(r) from the hydrodynamic models at short distances r reflect the underlying many-body fluid structure, and are found to be correlated to differences in the local available volume. The procedure used here to determine the pair diffusion coefficients can also be used for single-particle diffusion in confinement with spherical symmetry.

preprint2009arXiv

A one-dimensional dipole lattice model for water in narrow nanopores

We present a recently developed one-dimensional dipole lattice model that accurately captures the key properties of water in narrow nanopores. For this model, we derive three equivalent representations of the Hamiltonian that together yield a transparent physical picture of the energetics of the water chain and permit efficient computer simulations. In the charge representation, the Hamiltonian consists of nearest-neighbor interactions and Coulomb-like interactions of effective charges at the ends of dipole ordered segments. Approximations based on the charge picture shed light on the influence of the Coulomb-like interactions on the structure of nanopore water. We use Monte Carlo simulations to study the system behavior of the full Hamiltonian and its approximations as a function of chemical potential and system size and investigate the bimodal character of the density distribution occurring at small system sizes.

preprint2008arXiv

Exploration of Effective Potential Landscapes using Coarse Reverse Integration

We describe a reverse integration approach for the exploration of low-dimensional effective potential landscapes. Coarse reverse integration initialized on a ring of coarse states enables efficient "navigation" on the landscape terrain: escape from local effective potential wells, detection of saddle points, and identification of significant transition paths between wells. We consider several distinct ring evolution modes: backward stepping in time, solution arc--length, and effective potential. The performance of these approaches is illustrated for a deterministic problem where the energy landscape is known explicitly. Reverse ring integration is then applied to "noisy" problems where the ring integration routine serves as an outer "wrapper" around a forward-in-time inner simulator. Three versions of such inner simulators are considered: a system of stochastic differential equations, a Gillespie--type stochastic simulator, and a molecular dynamics simulator. In these "equation-free" computational illustrations, estimation techniques are applied to the results of short bursts of "inner" simulation to obtain the unavailable (in closed form) quantities (local drift and diffusion coefficient estimates) required for reverse ring integration; this naturally leads to approximations of the effective landscape.