Polarisation forces

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. 

The following discussion is limited to a couple of other dipoles, in order to illustrate some of the less straightforward cases. The balance of factors leading to a particular chemo- and regioselectivity is often close—the choice of which pair of frontier orbitals to take is sometimes difficult, the fact that some frontier orbitals are not strongly polarised forces us to judge each case carefully on its merits, and the outcome is not quite always in agreement with the experimental results. There is more discussion of these and all the other cases in Houk s... 

London(11), it appears that at large internuclear distances, where the repulsive forces of the first approximation have aheady practically gone to zero, weak attractive forces occur even for these pairs of atoms, giving rise to a flat minimum of the energy curve there. These forces may be termed polarisation forces or van der Wools farces since they are the cause of the van der Waals correction in the equation of state of the vapour which makes possible the liquefaction at sufficiently... 

It is also noteworthy that on the whole the activation energy for diffusion and for conductivity (Table 64) in ionic lattices is smaller than the activation energy for diffusion in metallic lattices. It is probable that this difference may be traced to the greater influence of polarisation forces in the case of salts, and to the greater density and so closer atomic packing in the case of many metals. 

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 existence of such forces also has to be assumed in quantum mechanics. For hydrogen atoms Wang has already calculated them, neglecting the chemical interactions. Since the polarisation forces l/R ) vanish more slowly at infinity... 

Especially for He we have to point out a special circumstance. He is special in that its polarisation forces and its mass are very small, so that for the motion of one atom in the field of another already the zero-point vibration fills a very considerable portion of the potential well. To give a rough indication, our potential for He has been presented in Fig. 1 as it looks in the boundaries (13) given above. Also a level Eo has been drawn, which has been determined based on the quantum condition... 

If, next to the dipole effects, we take into account the polarisation effects of the substituents, then nothing substantial changes in the energetics. We can estimate from the polarisabilities, which are known from the refraction equivalents, the dipole moments of the groups and their distances that the polarisation forces will be smaller than the pure dipole interactions. For instance, in the case of dichloroethylene, taking the polarisability into account, the maximum for the cis-position is somewhat lower, and in the case of compound (4) the two minima are deepened (order of magnitude after the change about 500 and 300 cal respectively). 

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. ... 

Lamoureux G, Roux B. (2006) Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarisable force BAA. J Phys Chem B WO 3308-3322. 

---