Merge branch 'master' of https://github.com/wehs7661/ensemble_md into…

… add_tests

ensemble_md/data/CB7-10/complex/complex.top

@@ -0,0 +1,15 @@
+#include "amber03.ff/forcefield.itp"
+#include "/ocean/projects/bio230014p/wehs7661/EEXE_experiments/CB7-10/complex/Prep/Topology/itp_files/10atomtypes.itp"
+#include "/ocean/projects/bio230014p/wehs7661/EEXE_experiments/CB7-10/complex/Prep/Topology/itp_files/cucurbit_7_uril_GMX.itp"
+#include "/ocean/projects/bio230014p/wehs7661/EEXE_experiments/CB7-10/complex/Prep/Topology/itp_files/10_GMX.itp"
+#include "amber03.ff/tip3p.itp"
+#include "amber03.ff/ions.itp"
+
+[ system ]
+Complex of CB7 and 10.acpype with solvent in water
+
+[ molecules ]
+HOS     1
+MOL     1
+SOL         2711
+CL               3

ensemble_md/data/CB7-10/complex/fixed_weight.mdp

@@ -0,0 +1,122 @@
+; Run control
+integrator = md-vv
+tinit = 0
+dt = 0.002
+nsteps = 100000000
+nstcomm = 10
+
+; Output control
+nstlog = 1000
+nstcalcenergy = 10
+nstenergy = 1000
+nstxout_compressed = 1000
+
+; Neighborsearching and short-range nonbonded interactions
+nstlist = 10
+ns_type = grid
+pbc = xyz
+rlist = 1.3
+
+; Electrostatics
+coulombtype = PME-switch
+rcoulomb_switch = 0.88
+rcoulomb = 0.9
+
+; van der Waals
+vdw_type = switch
+rvdw_switch = 0.85
+rvdw = 0.9
+
+; Apply long range dispersion corrections for Energy and Pressure
+DispCorr = AllEnerPres
+
+; Spacing for the PME/PPPM FFT grid
+fourierspacing = 0.12
+
+; EWALD/PME/PPPM parameters
+pme_order = 4
+ewald_rtol = 1e-05
+ewald_geometry = 3d
+epsilon_surface = 0
+optimize_fft = yes
+
+; Temperature coupling
+tcoupl = nose-hoover
+nsttcouple = 10
+tc_grps = System
+tau_t = 1.0
+ref_t = 300
+
+; Pressure coupling is on for NPT
+pcoupl = no
+
+gen_vel = yes
+gen_temp = 300
+gen_seed = -1
+
+; options for bonds
+constraints = h-bonds
+
+; Type of constraint algorithm
+constraint_algorithm = shake
+shake_tol = 1e-05
+
+; Free energy calculation
+free_energy = expanded
+calc_lambda_neighbors = -1
+sc_alpha = 0.5
+sc_power = 1
+sc_sigma = 0.5
+couple_moltype = MOL
+couple_lambda0 = vdw-q
+couple_lambda1 = none
+couple_intramol = no
+init_lambda_state = 0
+nstdhdl = 100
+dhdl_print_energy = total
+
+; Seed for Monte Carlo in lambda space
+symmetrized_transition_matrix = no
+nst_transition_matrix = 100000
+; wl-scale                 = 0.8
+; wl-ratio                 = 0.7
+; init-wl-delta            = 10
+
+; expanded ensemble variables
+nstexpanded = 100
+lmc_stats = no
+lmc_move = metropolized-gibbs
+; lmc-weights-equil       = wl-delta
+; weight-equil-wl-delta   = 0.001
+; wl-oneovert             = yes
+
+; lambda-states          = 1      2      3      4      5      6      7      8      9      10     11     12     13     14     15     16     17     18     19     20     21     22     23     24     25     26     27     28     29     30     31     32     33     34     35     36     37     38     39     40
+
+coul_lambdas = 0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
+vdw_lambdas = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.1 0.2 0.3 0.4 0.45 0.5 0.55 0.6 0.63 0.66 0.69 0.72 0.75 0.78 0.81 0.84 0.88 0.92 1.0
+
+; PULL CODE
+pull = yes
+pull_ngroups = 2
+pull_ncoords = 1
+pull_group1_name = HOS
+pull_group2_name = MOL
+pull_pbc_ref_prev_step_com = yes
+
+pull_coord1_groups = 1 2
+pull_coord1_type = umbrella
+pull_coord1_geometry = distance
+pull_coord1_dim = Y Y Y
+pull_coord1_origin = 0.0 0.0 0.0
+pull_coord1_vec = 0.0 0.0 0.0
+pull_coord1_start = yes
+pull_coord1_init = 0
+pull_coord1_rate = 0
+pull_coord1_k = 0
+pull_coord1_kB = 1000
+pull_nstfout = 400000
+pull_nstxout = 1000
+pull_print_ref_value = yes
+
+restraint_lambdas = 0.0 0.0 0.0 0.0 0.0 0.0 0.01 0.05 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.95 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
+init_lambda_weights = 0.0 57.88597 112.71883 163.84425 210.48097 253.80261 294.79849 333.90408 370.82669 406.02515 438.53116 468.53751 496.24649 521.58417 544.57404 565.26697 583.7337 599.60651 613.43958 624.70471 633.95947 638.29785 642.44977 646.33551 649.91626 651.54779 652.93359 654.13263 654.94073 655.13086 655.07239 654.66443 653.68683 652.32123 650.72308 649.2381 647.94586 646.599 645.52063 643.99133

ensemble_md/data/CB7-10/complex/itp_files/10atomtypes.itp

@@ -0,0 +1,9 @@
+; 10_GMX.itp created by acpype (Rev: 7268) on Fri May 17 15:26:08 2013
+
+[ atomtypes ]
+;name   bond_type     mass     charge   ptype   sigma         epsilon       Amb
+ c3       c3          0.00000  0.00000   A     3.39967e-01   4.57730e-01 ; 1.91  0.1094
+ n4       n4          0.00000  0.00000   A     3.25000e-01   7.11280e-01 ; 1.82  0.1700
+ hx       hx          0.00000  0.00000   A     1.95998e-01   6.56888e-02 ; 1.10  0.0157
+ hn       hn          0.00000  0.00000   A     1.06908e-01   6.56888e-02 ; 0.60  0.0157
+