Drift-oscillatory steering with the forward-reverse method for calculating the potential of mean force
Mostafa NategholEslam1, Bryan W Holland, C G Gray
1Department of Physics and Biophysics Interdepartmental Group, University of Guelph, Guelph, Ontario N1G 2W1, Canada. snategho@uoguelph.ca
Abstract:
We present a method that enables the use of the forward-reverse (FR) method of Kosztin et al. on a broader range of problems in soft matter physics. Our method, which we call the oscillating forward-reverse (OFR) method, adds an oscillatory steering potential to the constant velocity steering potential of the FR method. This enables the calculation of the potential of mean force (PMF) in a single unidirectional oscillatory drift, rather than multiple drifts in both directions as required by the FR method. By following small forward perturbations with small reverse perturbations, the OFR method is able to generate a piecewise reverse path that follows the piecewise forward path much more closely than any practical set of paths used in the FR method. We calculate the PMF for four different systems: a dragged Brownian oscillator, a pair of atoms in a Lennard-Jones liquid, a Na(+)-Cl⁻ ion pair in an aqueous solution, and a deca-alanine molecule being stretched in an implicit solvent. In all cases, the PMF results are in good agreement with those published previously using various other methods, and, to our knowledge, we give for the first time PMFs calculated by nonequilibrium methods for the Lennard-Jones and Na(+)-Cl⁻ systems.
Related Concept Videos
Two-Dimensional Force System: Problem Solving
The first step to solving a two-dimensional force system problem is to draw a free-body diagram of the object under consideration. This diagram helps identify all the external forces acting on the object, including their...
Root-Locus Method
This system can be represented by a block diagram,...
Central-Force Motion
Two-Dimensional Force System
Three-Dimensional Force System:Problem Solving
To solve a three-dimensional force system, first resolve each force into its respective scalar components. Do this using...
Time-Domain Interpretation of PD Control
Consider the example of control of motor torque. Initially, a positive...
