Effective force field

Effective force fields that include the new potential functions based on pair interaction energies are used to calculate the dehydration of hexahydrate complexes and they are compared to DFT results for the same reaction. It is concluded that the new functions may appropriately describe high temperature and low density states, where repulsive forces start to be significant in comparison to attractive forces. 

MD simulations are used to further test the new effective force fields. The main features observed are given by changes in the dynamic behavior. The new functions yield shorter characteristic lifetimes for water molecules in the first ionic shell than those found when the 12-6 function is used to represent the short-range interactions. Although more studies are required, these simulations indicate that combinations of short-range interaction exponents such as those proposed here may be more suitable for the representation of aqueous electrolyte solutions at high temperatures. 

Twenty years ago Car and Parrinello introduced an efficient method to perform Molecular Dynamics simulation for classical nuclei with forces computed on the fly by a Density Functional Theory (DFT) based electronic calculation [1], Because the method allowed study of the statistical mechanics of classical nuclei with many-body electronic interactions, it opened the way for the use of simulation methods for realistic systems with an accuracy well beyond the limits of available effective force fields. In the last twenty years, the number of applications of the Car-Parrinello ab-initio molecular d3mam-ics has ranged from simple covalent bonded solids, to high pressure physics, material science and biological systems. There have also been extensions of the original algorithm to simulate systems at constant temperature and constant pressure [2], finite temperature effects for the electrons [3], and quantum nuclei [4]. 

Medium-ring monocyclic diamines (lOa-e) have pA H.i values 12 in aqueous solution, so they cannot be measured by potentiometric methods <88CC1528>. Studies in DMSO with the proton sponge, 2,7-dimethoxy-l,8-bis(dimethylamino)naphthalene, indicate that (lOe) has the highest pATn value in this series. This results from internal hydrogen bonding as indicated by NMR measurements and effective force-field calculations. 

Force match (FM) method for multiscale simulation was initially introduced by Voth et al. Flowever the FM method was proposed as an extension of least square force match approach originally developed by Ercolessi and Adams. But the method proposed by Voth et al. can determine pair wise effective force field from a given trajectory and force data. The origin of this force data can be obtained from any type of calculations such as ab-initio MD simulation, path integral MD or from the coarse graining of atomistic data. ... 

The results for the unsaturated systems PBD and PIP are less easily assessed than for aPP, PE, and PIB. The simulation results are partly in good agreement with experiment, and partly defective. Especially for the unsaturated carbons, large deviations from experimental data are found. To identify the origin of the difficulties (correlation effects , force field ), and for a final decision on the validity of the simulation method for this class of systems, additional investigations will be necessary. 

The choice of the adjustable parameters used in conjunction with classical potentials can result to either effective potentials that implicitly include the nuclear quantization and can therefore be used in conjunction with classical simulations (albeit only for the conditions they were parameterized for) or transferable ones that attempt to best approximate the Born-Oppenheimer PES and should be used in nuclear quantum statistical simulations. Representative examples of effective force fields for water consist of TIP4P (Jorgensen et al. 1983), SPC/E (Berendsen et al. 1987) (pair-wise additive), and Dang-Chang (DC) (Dang and Chang 1997) (polarizable, many-body). The polarizable potentials contain - in addition to the pairwise additive term - a classical induction (polarization) term that explicitly (albeit approximately) accounts for many-body effects to infinite order. These effective potentials are fitted to reproduce bulk-phase experimental data (i.e., the enthalpy at T = 298 K and the radial distribution functions at ambient conditions) in classical molecular dynamics simulations of liquid water. Despite their simplicity, these models describe some experimental properties of liquid... 

The present high cost of full CASSCF direct dynamics means that it is not possible to use such calculations to run large numbers of trajectories. As a result it cannot be used to build up experience of the types of effects to be found in dynamical studies of organic photochemistry, and in their interpretation. This problem can be remedied by performing calculations using the MMVB force field [63,64]. 

These numbers are not negligible. At present such effects are on the average compensated by other force field terms through empirical parametrization. 

Conformational Adjustments The conformations of protein and ligand in the free state may differ from those in the complex. The conformation in the complex may be different from the most stable conformation in solution, and/or a broader range of conformations may be sampled in solution than in the complex. In the former case, the required adjustment raises the energy, in the latter it lowers the entropy in either case this effect favors the dissociated state (although exceptional instances in which the flexibility increases as a result of complex formation seem possible). With current models based on two-body potentials (but not with force fields based on polarizable atoms, currently under development), separate intra-molecular energies of protein and ligand in the complex are, in fact, definable. However, it is impossible to assign separate entropies to the two parts of the complex. 

Large stepsizes result in a strong reduction of the number of force field evaluations per unit time (see left hand side of Fig. 4). This represents the major advantage of the adaptive schemes in comparison to structure conserving methods. On the right hand side of Fig. 4 we see the number of FFTs (i.e., matrix-vector multiplication) per unit time. As expected, we observe that the Chebyshev iteration requires about double as much FFTs than the Krylov techniques. This is due to the fact that only about half of the eigenstates of the Hamiltonian are essentially occupied during the process. This effect occurs even more drastically in cases with less states occupied. 

In this model of electrostatic in teraction s, two atoms (i and j) have poin t charges tq and qj. The magnitude of the electrostatic energy (V[. , [ ) varies inversely with the distance between the atoms, Rjj. fh e effective dielectric constant is . For in vacuo simulations or simulation s with explicit water rn olecules, the den om in a tor equals uRjj, In some force fields, a distance-dependent dielectric, where the denominator is uRjj Rjj, represen is solvent implicitly. 

United Atom force fieldsare used often for biological polymers. In th esc m oleciiles, a reduced ii nm ber of explicit h ydrogen s can have a notable effect on the speed of the calculation. Both the BlOn and OPLS force fields are United Atom force fields. AMBER con tain s both aU nited and an All Atom force field. 

An atom lhai hassp h yhridi/alion lends lo be coplanar wilh its attached aloms. This effect isaccoiinled for by improper torsions in Olher force fields and by oiil-of-plane-bending inlcraclions in... 

---