_{1}

^{*}

The dielectric constant in Coulomb’s Law, D, can quantify an empirical reduction of force. It can also quantify a reduction of electrostatic field as seen in classical electrostatic theory where the induced charge layer is assumed to be infinitely thin. The two approaches exemplify two traditions that have been used in parallel for decades. They produce Potential Energy Functions (PEFs) that differ by a factor of the permittivity, εr. The classical electrostatic theory result can be incorporated into force field models with an effective dielectric function, Deff, which spans the induced charge layer and accommodates both traditions. The Deff function increases the magnitude of local terms as compared with cumulative long distance terms. It is shown that the Deff function reduces distance dependence of the radial PEF within the induced charge layer and improves computational stability for some systems including substrate in dilute salt solution. End use applications include pharmaceutical development (e.g. protein calculations with docking), materials development, solvation energy calculations and QM/MM calculations.

The dielectric effect [

The dielectric effect is the macroscopic manifestation of induced charge. Electrostatic interactions (and the associated dielectric effect) are foundational in all areas of chemistry including organic, physical, inorganic, biochemistry and electrochemistry. They are broadly used in molecular modeling including models for force fields, reactions, molecular structure, and ions in solution. Accompanying software is used in a wide variety of applications including the development of materials (material science), biotech, nanotech, and pharmaceuticals. Specific applications range from circuit board design to gas chromatography.

Many textbooks [

where F is the force, the dielectric constant, D, is experimentally determined, r is a unit vector pointing from q_{1} to q_{2}, ε_{0} is the vacuum permittivity, and the electrostatic field, E, was defined as the force acting on a probe charge.

Two interpretations of Equation (1) are common in the literature. The “force between two charges” interpretation, depicted in

These two interpretations represent two traditions that have been used in parallel for decades. They roughly correlate with scale where researchers working with interactions within induced charge layers (computational chemistry) tend to use the interpretation of

This paper proposes a function to transition between these two traditions. The transition function improves computational stability for many systems. Section 2 contains background material for a wide range of researchers. It also shows how the parallel traditions are consistent with distinct interpretations of the term “probe charge” used in mid-20th century textbooks. Section 3 derives the potential energy equation (PEF) for the two traditions and shows how they differ by a factor of the dielectric constant. Section 4 examines the distinguishability of the traditions. Section 5 introduces a function to transition between the traditions. Section 6 describes an analysis tool, the radial PEF, which quantifies the effect of transitioning between the traditions. Section 7 shows applica- tions in chemistry and Section 8 is the discussion.

Coulomb’s Law [

When charges are not in vacuum, the surrounding medium polarizes, forming induced charge. The macroscopic manifestation of this induced charge is the dielectric effect. It has been shown that this induced charge can be handled like any other charge ( [

The dielectric effect is due to induced polarization including electronic polarization, field-induced rotation of molecules having permanent dipole moments, and internal rotations and vibrations that produce net dipoles, all of which oppose and reduce the imposed field. The sources of permittivity can be divided into two groups: 1) polarization stemming from rotations of dipoles, ε_{ν} (i.e. due to movement of atoms), and 2) polarization stemming from electronic polarization and intra-molecular vibrations, ε_{∞} (primarily due to movement of electrons). Experimentally, ε_{∞} and ε_{ν} can be distinguished, as reviewed in Gilson and Honig [

In a molecular simulation, the induced charge can be determined by summing all charges within a sphere centered about the ion and averaging over time. As the radius of the sphere increases, the magnitude of the induced charge reaches a constant value, q_{induced}. The total charge is q_{total} = q + q_{induced}. Electrostatic analysis of simulation results is exemplified in Warshel and Levitt [

Quantitative models relating macroscopic permittivity to microscopic polarization were developed by Debye, Onsager, Kirkwood, and Fröhlich among others. Debye developed models relating macroscopic permittivity with the polarizability of a molecule, its permanent dipole moment, and the energy of thermal agitation [

Onsager introduced the concept of a reaction field [_{r}. The internal field in the cavity consists of the field that would be produced by the applied external field and the reaction field set up in the cavity by the polarization which the dipole induces in its surroundings [

Debye and Onsager used macroscopic arguments for calculating the local field which limited method applicability to fluids. Kirkwood and Fröhlich used statistical methods for calculating the local field that extended method applicability. A review is provided by Hill, Vaughan, Price, and Davies [

Relative permittivity, ε_{r}, quantizes the linear reduction in electrostatic field when charge is immersed in media. It can be experimentally measured from the change in capacitance when dielectric material is introduced [

More formally, permittivity, ε, is the ratio of the electric flux density, D, to the electrostatic field, E. It is separated into two factors, the vacuum permittivity, ε_{0}, which incorporates units and conventions and the relative permittivity, ε_{r}, which is unit-less.

Intuitively, D is the electrostatic field due to initial charge; D has the relation to initial charge that E has to total charge. That is, just as the divergence of E is the total charge, the divergence of D is the initial charge ( [

In many cases the relation between D and E may be considered to be linear, which is the first term when D is expanded in powers of the electrostatic field. This linearity results from the normally small size of external electrostatic fields compared to internal molecular fields ( [

It is often possible to determine the magnitude of induced charge. The electrostatic field stems from both initial and induced charge. When the initial charge is known, the induced charge can be back calculated. For example, for a point charge in dielectric media the electrostatic field is linearly reduced by ε_{r}. This must be equal to the field stemming from all charges, initial and induced.

This equality shows the relation between the total charge and the initial charge ( [

In the case of electronic polarization with permittivity of_{r} = 80, the induced charge is 98.75% as large.

The PEF is the potential energy as a function of the position of the charges. The potential energy is the integral of the negative force between charges as the charges are brought in from infinity to the position, r_{12}. For two point charges in infinite dielectric the potential energy is given by V_{PEF}.

In the literature, D is frequently used to quantize the reduction of potential energy. The dielectric media is assumed to be infinite, the charges are point charges that are fully formed at infinity, these charges are surrounded by induced charge at all positions (including infinity), and the zero point energy is defined as the state where the point charges are an infinite distance apart from each other.

In practice, D is set to a constant or a function (e.g. Gelin and Karplus [

On the face of it, the work equation has the same permittivity dependence as the PEF. However, examination of the derivation of the work equation shows that a bare charge is brought in from infinity. In the PEF, all charges are surrounded by induced charge at all positions. The permittivity dependence of the work equation and the PEF are not the same.

Many textbooks use the dielectric version of Coulomb’s Law, Equation (1), where the electrostatic field is defined as the force acting on a probe charge, see Smythe ( [

Two interpretations of the dielectric Coulomb’s Law are depicted in _{r}.

Since Equation (1) holds for a probe charge, the definition of “probe charge” reflects on the intended interpretation of Equation (1). The probe charge can be interpreted either as a small charge with surrounding induced charge, or as an individual charge exclusive of induced charge. That is, while one of the charges in Equation (1) has surrounding induced charge, the other charge does or does not, depending on the interpretation of “probe charge”. Two definitions for “probe charge” are depicted in

Smythe defines a probe charge to be a charge so minute that no redistribution of charge takes place due to its

presence ( [

Classically, a charge always induces charge about itself, regardless of size: The magnitude of induced charge is proportional to the magnitude of the initial charge. As the initial charge becomes smaller, the induced charge becomes proportionally smaller. This interpretation of “small charge” maps onto

Alternately, one can imagine an abstract, hypothetical charge that does not induce charge about itself. This hypothetical charge interpretation is consistent with

These different interpretations of the probe charge are consistent with the different usages of Equation (1) that are seen in the literature. The “small charge with surrounding induced charge” definition is consistent with a “total force between two charges” interpretation. The “small charge without surrounding induced charge” is consistent with the “total force acting on an individual charge” interpretation.

Probe charges are no longer used in textbooks on electromagnetism [

At the molecular scale, some of the approximations and assumptions used in macroscopic analysis do not hold. In particular, at the macroscopic scale the induced charge layer is assumed to be infinitely thin, whereas in chemistry much of the data is based on interactions within the induced charge layer of multiple charges.

In the derivation of the PEF, the force is the total force between charges. The two traditions use different equations for this force as depicted in

In _{1} and q_{2} is given by Equation (1). In molecular modeling, D is generally set to ε_{r}, D = ε_{r}.

The potential energy, V_{Fig}_{ 3}_{(a}_{)}, is found by integration.

In _{1} and q_{2} is determined by substituting the total charge (initial and induced) into Coulomb’s Law, i.e. substitute total charge from Equation (5) in for both charges in Equation (2).

Again, the potential energy is found by integration.

The two traditions in _{r}.

When permittivity is large, the two traditions are empirically similar. For example, in water the permittivity is about 80,

The two traditions can be distinguished in cases where ε_{r} is small resulting in terms that are sufficiently large to be measured. That is, for electronic polarization where ε_{r} = ~2, the PEFs are reduced by a factor of 2 and 4 (Equations (8) and (10) respectively) and this difference can be experimentally measured.

The traditions can also be distinguished in cases where, even though ε_{r} is large, there are a large number of interactions and the sum of the interactions is measureable as discussed below in the section on the radial PEF.

When the dielectric constant is measured in circumstances similar to those where it is used, then the two traditions are similar. For example,

When polarization is approximated to be pairwise separable, a pairwise function can be used to describe the dielectric effect. Then force field models can incorporate classical electrostatic theory with an effective dielectric, D_{eff}.

D_{eff} is the effective dielectric constant at distance, r (See

increases and the effective dielectric increases as the induced charge layer is traversed. At distances outside the induced charge layer, the macroscopic value is used,

Many molecular simulations have terms spanning the range between molecular and classical distances. The effective dielectric handles interactions at all distances and enables the utilization of both traditions. For example a simulation of substrate in condensed phase typically has media that extends to infinity. In this simulation some terms are for interactions within induced charge layers and others are for interactions beyond the induced charge layer.

The question remains, under what circumstances is it important to transition between the traditions?

The magnitude of the potential energy between two charges decreases with distance and long distance terms are very small. However, the number of terms at this distance is large. Is the number of terms large enough for long distance terms to be important? The radial PEF quantifies the net effect of these terms. The radial PEF is related to the radial distribution function. It is independent of approximation and can be applied to an exact wavefunc- tion as well as to a trajectory from a force field model.

For a charge, q, and a distance, r, the radial PEF at r is the sum of potential energy terms for all charges at distance r from q. The integrated radial PEF is the sum over r of the radial PEFs. Examination of the PEF shows that it is equal to the sum of all integrated radial PEF’s divided by 2 since every interaction is counted twice.

In the radial PEF, the magnitude of each charge-charge interaction decreases linearly with inverse distance, but the number of interactions increases by r^{2}, i.e. for a constant density of charges, the number of charges at distance, r, increases in proportion to the surface area of a sphere of radius r centered at q, 4πr^{2}. As a result, the radial PEF increases linearly with distance. The sum over all distances is the integral over r resulting in an inte- grated radial PEF that depends on r^{2}. Ostensibly this means that the PEF also depends on r^{2} and extending to in- finity produces an infinite PEF.

For a neutral system, there are equal numbers of positive and negative charges. In most chemical systems, counter-ions are associated with each other and the terms in a region approximately cancel. That is, for large dis- tances the PEF is the sum of a very large number of positive and negative numbers which can affect computa- tional stability. Algorithmic methods can be used to improve computational stability.

Any charge distribution can be represented as an infinite series consisting of the net charge, a dipole, a qua- drapole, and higher order multipoles ( [

The net charge in a region depends on the characteristics of the system; while it tends toward zero, there may be a net charge in a region (e.g. in a dilute salt solution). The dipole moment is significant within the induced charge layer; indeed the induced charge layer frequently consists of induced dipoles. Outside the induced charge layer, dipole moments return to bulk structure.

The potential energy between two charges varies as

D_{eff} can be modeled as linear with distance. When this D_{eff} is incorporated into the radial PEF, it reduces distance dependence by a factor of r for interactions within the induced charge layer. Then the charge term is independent of distance and the dipole term varies as_{eff} levels out to the macroscopic value further mod- erating the PEF as compared with the PEF produced using the tradition of

When the induced charge layer is thin, then the transfer between traditions occurs within a short distance. An example is electronic polarization. Outside the induced charge layer, macroscopic equations hold and

When the induced charge layer is thick, the number of electrostatic interactions at a given distance becomes important as seen in

The charge term is significant when the sum of charges in a region is not zero (for example, when charges are not closely associated with their counter charges as seen in dilute salt solutions).

The dipole term is a significant portion of the radial PEF for many calculations. As distance increases, dipole moments continue to contribute to the induced charge until the edge of the induced charge layer is reached and the dipoles return to bulk structure. The distance dependence of the dipole radial PEF is_{eff} reaches the macroscopic value.

The quadrapole radial PEF will tend to drop out for large distances due to the

Examination of Coulomb’s Law as it is used in different fields elucidates the dielectric effect in different cir- cumstances. This material can be used to inform the parameterization of force field models, solvent models, and matrix models where distances range from molecular scales to macroscopic scales.

Term | Distance Dependence of Radial PEF inside the Induced Charge Layer | ||
---|---|---|---|

Radial PEF | Radial PEF with D_{eff} | Integrated Radial PEF with D_{eff} | |

Net Charge | r^{ } | Constant | r |

Dipole Moments | Constant | ln(r) | |

Quadrapole Moments |

The dielectric effect is a linear response and different fields of research focus on the linear reduction of different quantities. The dielectric constant can represent a linear reduction of force between charges, consistent with Figures 1(a)-3(a). The dielectric constant can also represent a linear reduction of electrostatic field consistent with Figures 1(b)-3(b).

For the purposes of this Discussion, D quantifies the linear reduction of the potential energy as seen in Equation (6). Setting D = ε_{r} results in Equation (8) and setting _{r}. The traditions have been used in parallel for decades and are self-consistent.

Equation (8) (D = ε_{r}) has been in ubiquitous use in computational chemistry literature for over 20 years. It is used in force field model software for pharmaceuticals, biotech, nanotech, and materials (material science). For example, a web search finds over 12 million cached pages that cite software (AMBER or CHARMM) where users set D = ε_{r}.

Equation (10) (

Smythe’s classic textbook on electromagnetism [

Examination of multiple textbooks from several fields reveals that both traditions are currently in use, possibly stemming from different interpretations of Smythe and his contemporaries. The two traditions roughly correlate with scale where researchers working with interactions within induced charge layers tend to use the interpretation of Figures 1(a)-3(a). Textbooks on macroscopic electromagnetic theory assume the induced charge layer is infinitely thin and are consistent with the interpretation of Figures 1(b)-3(b).

The two traditions may have developed due to different experiences of the magnitude of electrostatic interactions. That is, experience with interactions within the induced charge layer verses experience with interactions outside the induced charge layer. Complexity is compounded if the induced charge layer is thicker than thought: in water a linear D_{eff} and the macroscopic result that D = 6400 would result in an induced charge layer that is thousands of angstroms wide.

The two traditions can be rationalized with a distance dependent function as depicted in _{eff}, is measured at that distance. D_{eff} increases as distance increases and the induced charge layer is traversed. As distance reaches the edge of the induced charge layer, D_{eff} reaches the classical theoretical result,_{eff} function will vary with substrate and matrix.

In PEF calculations, D_{eff} increases the magnitude of local terms as compared to cumulative long distance terms. For example, in water the cumulative long distance PEF is decreased by a factor of 80, reducing the impact of mid and long distance terms on convergence and local effects.

Use of the interpretation in _{eff}. Use of D_{eff} is important wherever cumulative long distance terms are significant.

Special appreciation is extended to Dr. Dana Backman for editorial comment.