Polarisable force fields

This situation, despite the fact that reliability is increasing, is very undesirable. A considerable effort will be needed to revise the shape of the potential functions such that transferability is greatly enhanced and the number of atom types can be reduced. After all, there is only one type of carbon it has mass 12 and charge 6 and that is all that matters. What is obviously most needed is to incorporate essential many-body interactions in a proper way. In all present non-polarisable force fields many-body interactions are incorporated in an average way into pair-additive terms. In general, errors in one term are compensated by parameter adjustments in other terms, and the resulting force field is only valid for a limited range of environments. 

In order to fully explain the way in which molecular level structures and interactions bridge to macroscopic liquid-state properties, Voth stressed that a wide range of length and timescale calculations are required to predict different properties, like the interfacial tension, selfdiffusion and viscosity [30]. Both Voth and Borodin demonstrated the necessity of using polarisable force fields to accurately predict several properties of ionic liquids and their mixtures, namely, the transport of ionic species in the bulk or at metallic interfaces [96]. 

In summary, various approaches were developed and validated for ionic liquids over the last years. It was shown that dispersion corrected KS-DFT approaches allow reliable results for ionic liquids. " Recently, a comparison of trajectories obtained from ab initio molecular dynamics simulations with and without a dispersion correction revealed that the dynamics of the system is more accurately described for the dispersion corrected one, ° which highlights the necessity of dispersion corrected approaches in ab initio molecular dynamics simulations. Polarisable force fields were developed which allow the investigation of various ionic liquids.Cheaper nonpolarisable force fields are also available for a broad range of ionic liquids. However, the latter force fields tend to show too... 

The answers which are discussed in this book are based on the following three concepts. The first one introduces ion specificity through collective dispersion type interactions an ion specificity is thereby obtained by the explicit consideration of the size and the polarisability of the ions. Based on molecular dynamics (MD) simulation with polarisable force fields, Jungwirth and Tobias state that induction interactions close to the free surface may be responsible for the preference of heavier ions at interfacial solvation sites. The asymmetric, incomplete solvation shell induces a sizable dipole on the anion at the interface, which is assumed to be the driving force for the interfacial propensity of the ions. MD simulation provides a very detailed picture of the interfacial architecture however, the results depend strongly on the interaction potentials which are not exactly known. Hence, experiments are needed to verify the predictions. Indeed, this task is challenging and many sophisticated surface analytical techniques, even when pushed to the limits, may still yield only inconclusive results. 

These observations and interpretations defined the textbook picture. This point of view has been challenged by the progress in understanding atmospheric reactions which in turn motivated molecular dynamics (MD) simulations. MD simulations using polarisable force fields predict that soft ions such as halides are enriched at the interface with non-monotonic ion profiles. The book chapter of Pavel Jungwirth covers this in greater detail. 

Let us consider the simplest surface that shows ion-specific adsorption, namely the water-air interface. In a by now classical series of papers, Jungwirth and co-workers have shown that iodide ions do adsorb at the air-water interface, in strong contrast with the traditional view. Those simulations were performed with polarisable force fields, while the non-polarisable force fields employed at that time did not show adsorption of iodide. It was concluded that the polarisability plays a dominant role in the adsorption mechanism. Let us reconsider that problem using our novel thermodynamically optimised force fields discussed in the earlier section. We show results for the potential of mean force of a single ion at an air-water interface, calculated using umbrella sampling and the WHAM method. ... 

But the methods have not really changed. The Verlet algorithm to solve Newton s equations, introduced by Verlet in 1967 [7], and it s variants are still the most popular algorithms today, possibly because they are time-reversible and symplectic, but surely because they are simple. The force field description was then, and still is, a combination of Lennard-Jones and Coulombic terms, with (mostly) harmonic bonds and periodic dihedrals. Modern extensions have added many more parameters but only modestly more reliability. The now almost universal use of constraints for bonds (and sometimes bond angles) was already introduced in 1977 [8]. That polarisability would be necessary was realized then [9], but it is still not routinely implemented today. Long-range interactions are still troublesome, but the methods that now become popular date back to Ewald in 1921 [10] and Hockney and Eastwood in 1981 [11]. 

Finally, there are groups of liquid crystals where, at the current time, force fields are not particularly useful. These include most metal-containing liquid crystals. Some attempts have been made to generalise traditional force fields to allow them to cover more of the periodic table [40, 43]. However, many of these attempts are simple extensions of the force fields used for simple organic systems, and do not attempt to take into account the additional strong polarisation effects that occur in many metal-containing liquid crystals, and which strongly influence both molecular structure and intermolecular interactions. 

We will not discuss the more obvious but very tedious and cumbersome task to improve force fields. Electrostatic interactions are the weakest parts of the force fields polarisability has to be introduced and dependable ways have to be found to include long range effects. 

Thus, in contrast to preceding MM approaches explicit treatment of electronic polarisability is integral to a semi-empirical QM approach and promises excellent prospects for quantitative theoretical modelling of carbohydrates across a range of condensed phase environments. The results of the PM3CARB-1 model do however indicate in line with classical force field approaches [65, 73] that perhaps greater... 

Intermolecular Energy decomposition analyses (EDA) are very useful approaches to calibrate force fields. Indeed, an evaluation of the different physical components of the interaction energy, especially of the many-body induction, is a key issue for the development of polarisable models. 

Trichlorobenzene (TCB), a well known termite exterminator, is prepared selectively from 1,2-dichlorobenzene (DCB) using zeolite K-L as a catalyst and monochloroacetic acid as a promoter. An attempt has been made to apply the combination of molecular graphics, force field calculations and quantum chemical calculations to understand the mechanism of selective chlorination of 1,2-DCB to 1,2,4-TCB over K-L promoted by monochloroacetic acid. It was found that the zeolite lattice plays an important role in polarising the molecules. The peculiar "barrel shaped pore architecture allows zeolite L to act as a reactor vessel where monochloroacetic acid, chlorine and 1,2-DCB can be accommodated on a molecular level. 

---