View Javadoc
1   // ******************************************************************************
2   //
3   // Title:       Force Field X.
4   // Description: Force Field X - Software for Molecular Biophysics.
5   // Copyright:   Copyright (c) Michael J. Schnieders 2001-2024.
6   //
7   // This file is part of Force Field X.
8   //
9   // Force Field X is free software; you can redistribute it and/or modify it
10  // under the terms of the GNU General Public License version 3 as published by
11  // the Free Software Foundation.
12  //
13  // Force Field X is distributed in the hope that it will be useful, but WITHOUT
14  // ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
15  // FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
16  // details.
17  //
18  // You should have received a copy of the GNU General Public License along with
19  // Force Field X; if not, write to the Free Software Foundation, Inc., 59 Temple
20  // Place, Suite 330, Boston, MA 02111-1307 USA
21  //
22  // Linking this library statically or dynamically with other modules is making a
23  // combined work based on this library. Thus, the terms and conditions of the
24  // GNU General Public License cover the whole combination.
25  //
26  // As a special exception, the copyright holders of this library give you
27  // permission to link this library with independent modules to produce an
28  // executable, regardless of the license terms of these independent modules, and
29  // to copy and distribute the resulting executable under terms of your choice,
30  // provided that you also meet, for each linked independent module, the terms
31  // and conditions of the license of that module. An independent module is a
32  // module which is not derived from or based on this library. If you modify this
33  // library, you may extend this exception to your version of the library, but
34  // you are not obligated to do so. If you do not wish to do so, delete this
35  // exception statement from your version.
36  //
37  // ******************************************************************************
38  package ffx.openmm.amoeba;
39  
40  import ffx.openmm.Context;
41  import ffx.openmm.Force;
42  
43  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_addParticle;
44  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_create;
45  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_destroy;
46  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setAwater;
47  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setDispoff;
48  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setEpsh;
49  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setEpso;
50  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setParticleParameters;
51  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setRminh;
52  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setRmino;
53  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setShctd;
54  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_setSlevy;
55  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_AmoebaWcaDispersionForce_updateParametersInContext;
56  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_KJPerKcal;
57  import static edu.uiowa.jopenmm.OpenMMAmoebaLibrary.OpenMM_NmPerAngstrom;
58  
59  /**
60   * Weeks-Chandler-Andersen Dispersion Force.
61   */
62  public class WcaDispersionForce extends Force {
63  
64    /**
65     * Create a new Amoeba WCA dispersion force.
66     */
67    public WcaDispersionForce() {
68      pointer = OpenMM_AmoebaWcaDispersionForce_create();
69    }
70  
71    /**
72     * Add a particle to the force field term.
73     *
74     * @param radius  The radius of the particle.
75     * @param epsilon The well depth of the particle.
76     */
77    public void addParticle(double radius, double epsilon) {
78      OpenMM_AmoebaWcaDispersionForce_addParticle(pointer, radius, epsilon);
79    }
80  
81    /**
82     * Set the parameters for a particle.
83     *
84     * @param index   The index of the particle to set.
85     * @param radius  The radius of the particle.
86     * @param epsilon The well depth of the particle.
87     */
88    public void setParticleParameters(int index, double radius, double epsilon) {
89      OpenMM_AmoebaWcaDispersionForce_setParticleParameters(pointer, index, radius, epsilon);
90    }
91  
92  
93    /**
94     * Set the water oxygen epsilon parameter.
95     *
96     * @param epso The water oxygen epsilon parameter.
97     */
98    public void setEpso(double epso) {
99      OpenMM_AmoebaWcaDispersionForce_setEpso(pointer, epso * OpenMM_KJPerKcal);
100   }
101 
102   /**
103    * Set the water hydrogen epsilon parameter.
104    *
105    * @param epsh The water hydrogen epsilon parameter.
106    */
107   public void setEpsh(double epsh) {
108     OpenMM_AmoebaWcaDispersionForce_setEpsh(pointer, epsh * OpenMM_KJPerKcal);
109   }
110 
111   /**
112    * Set the water oxygen radius parameter.
113    *
114    * @param rmino The water oxygen radius parameter.
115    */
116   public void setRmino(double rmino) {
117     OpenMM_AmoebaWcaDispersionForce_setRmino(pointer, rmino * OpenMM_NmPerAngstrom);
118   }
119 
120   /**
121    * Set the water hydrogen radius parameter.
122    *
123    * @param rminh The water hydrogen radius parameter.
124    */
125   public void setRminh(double rminh) {
126     OpenMM_AmoebaWcaDispersionForce_setRminh(pointer, rminh * OpenMM_NmPerAngstrom);
127   }
128 
129   /**
130    * Set the dispersion offset.
131    *
132    * @param dispoff The dispersion offset.
133    */
134   public void setDispoff(double dispoff) {
135     OpenMM_AmoebaWcaDispersionForce_setDispoff(pointer, dispoff * OpenMM_NmPerAngstrom);
136   }
137 
138   /**
139    * Set the water density parameter.
140    *
141    * @param awater The water density parameter.
142    */
143   public void setAwater(double awater) {
144     OpenMM_AmoebaWcaDispersionForce_setAwater(pointer, awater / (OpenMM_NmPerAngstrom * OpenMM_NmPerAngstrom * OpenMM_NmPerAngstrom));
145   }
146 
147   /**
148    * Set the Levy parameter.
149    *
150    * @param slevy The Levy parameter.
151    */
152   public void setSlevy(double slevy) {
153     OpenMM_AmoebaWcaDispersionForce_setSlevy(pointer, slevy);
154   }
155 
156   /**
157    * Set the overlap factor.
158    *
159    * @param shctd The overlap factor.
160    */
161   public void setShctd(double shctd) {
162     OpenMM_AmoebaWcaDispersionForce_setShctd(pointer, shctd);
163   }
164 
165   /**
166    * Update the parameters in the context.
167    *
168    * @param context The OpenMM context.
169    */
170   public void updateParametersInContext(Context context) {
171     if (context.hasContextPointer()) {
172       OpenMM_AmoebaWcaDispersionForce_updateParametersInContext(pointer, context.getPointer());
173     }
174   }
175 
176   /**
177    * Destroy the force.
178    */
179   public void destroy() {
180     if (pointer != null) {
181       OpenMM_AmoebaWcaDispersionForce_destroy(pointer);
182       pointer = null;
183     }
184   }
185 }