charge-charge, charge-dipole, dipole-charge, dipole-dipole interaction Masahiro Yamamoto Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University, Kyoto 606-8501, Japan (Dated: December 1, 2008) In the bulk and at the interface of polar solvents and/or ionic liquids, molecules feel the strong electric field and this strong field distort the electron cloud of molecule and induces dipole around atoms. To consider this polarization effect we should consider charge-charge, charge-dipole, dipole-charge, dipole-dipole interaction. I. CHARGE-CHARGE (C-C) INTERACTION The coulomb potential φ from the point charge z i e at r i is given by FIG. 1: Coulomb potential of negative charge φ C (|r − r i |)= z i e 4πϵ 0 1 |r − r i | (1) The electric field is given by the gradient of the potential E i C (r)= −∇φ C (|r − r i |)= z i e 4πϵ 0 r − r i |r − r i | 3 (2) If we put the other point charge z j e at at r j the electrostatic energy V cc is given by V cc = z j eφ C (|r j − r i |)= z i z j e 2 4πϵ 0 1 |r i − r j | (3) II. CHARGE-DIPOLE (C-D) INTERACTION In this interaction we consider two cases, i.e. (I) the coulomb potential interact with dipole µ j = ez j p j and (II) dipole field interact with point charge. In the case of (I), the C-D interaction becomes V cd = z i (−z j )e 2 4πϵ 0 1 |r j − r i | + z i z j e 2 4πϵ 0 1 |r j + p j − r i | (4) r j − r i ≡ r ji , |r ji | = r ji , (1 + x) −1/2 ≅ 1 − x/2+3x 2 /8... (5) V cd = − z i z j e 2 4πϵ 0 1 r ji − 1 q r 2 ji +2r ji · p j + p 2 j = − z i z j e 2 4πϵ 0 r ji " 1 − 1 (1 + 2r ji · p j /r 2 ji + p 2 j /r 2 ji ) 1/2 #
8
Embed
charge-charge, charge-dipole, dipole-charge, dipole-dipole ... · PDF filecharge-charge, charge-dipole, dipole-charge, dipole-dipole interaction Masahiro Yamamoto Department of Energy
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Masahiro YamamotoDepartment of Energy and Hydrocarbon Chemistry,
Graduate School of Engineering, Kyoto University, Kyoto 606-8501, Japan(Dated: December 1, 2008)
In the bulk and at the interface of polar solvents and/or ionic liquids, molecules feel the strong electric field andthis strong field distort the electron cloud of molecule and induces dipole around atoms. To consider this polarizationeffect we should consider charge-charge, charge-dipole, dipole-charge, dipole-dipole interaction.
I. CHARGE-CHARGE (C-C) INTERACTION
The coulomb potential φ from the point charge zie at ri is given by
The electric field is given by the gradient of the potential
EiC(r) = −∇φC(|r − ri|) =
zie
4πϵ0
r − ri
|r − ri|3(2)
If we put the other point charge zje at at rj the electrostatic energy Vcc is given by
Vcc = zjeφC(|rj − ri|) =zizje
2
4πϵ0
1|ri − rj |
(3)
II. CHARGE-DIPOLE (C-D) INTERACTION
In this interaction we consider two cases, i.e. (I) the coulomb potential interact with dipole µj = ezjpj and (II)dipole field interact with point charge.
Here the interaction tensors are given by [check this→OK]
Tij =1rij
(19)
Tαij = ∇αTij = −rij,αr−3
ij (20)
Tαβij = ∇α∇βTij = (3rij,αrij,β − r2
ijδαβ)r−5ij (21)
Tαβγij = ∇α∇β∇γTij = −[15rij,αrij,βrij,γ − 3r2
ij(rij,αδβγ + rij,βδγα + rij,γδαβ)]r−7ij (22)
Here ∇α means ∂∂rij,α
for α = x, y, z. Tαβγij will be used in the force formulation.
V. POLARIZATION
When a strong electric field is applied to an atom i, the electrons around atom i starts to deform and a dipolemoment may be induced.
The dipole moment of atom i in the α direction may be written as the superposition of the electric field EjC generated
by the charge of atom j and that EjD by the dipole of atom j
µi,α = µstati,α + µind
i,α ≅ µindi,α (23)
Usually we take µstatici,α = 0. The electric field at atom i can be written as
E(ri) =∑j(=i)
[EjC(ri) + Ej
D(ri)] (24)
5
������ ++++++d i p o l ei n d u c e d
+ +� �FIG. 5: What’s polarization?
=1
4πϵ0
∑j(=i)
{zje
ri − rj
|ri − rj |3− 1
|ri − rj |3
[µind
j − 3µind
j · (ri − rj)(ri − rj)|ri − rj |2
]}(25)
Eα(ri) =1
4πϵ0
∑j(=i)
−Tαijqj +
∑β
Tαβij µind
j,β
(26)
The induced dipole moment µindi,α is given by
µindi,α = αi
α,βEβ(ri) (27)
Here αiα,β polarizability tensor of atom i. If we assume
[αi] =
αi 0 00 αi 00 0 αi
(28)
then, we have
µindi = αiE(ri), µind
i,α = αiEα(ri) (29)
µindi,α (out) =
αi
4πϵ0
∑j( =i)
−Tαijqj +
∑β
Tαβij µind
j,β (input)
(30)
µindi,x (out) =
αi
4πϵ0
∑j( =i)
−T xijqj +
∑β
T xβij µind
j,β (input)
(31)
=αi
4πϵ0
∑j( =i)
xij
r3ij
qj +∑
β
[3xijrij,β
r5ij
− δxβ
r3ij
]µindj,β (input)
(32)
=αi
4πϵ0
∑j( =i)
xij
r3ij
qj +3xij
r5ij
∑β
rij,βµindj,β (input) −
µindj,x (input)
r3ij
(33)
The last equations should be solved self-consistently.In the fixed charge model the interaction of C-C is calculated by Ewald method. Then when we consider the induced
dipoles at the atom sites we should C-D, D-C, and D-D interaction. The total C-C, C-D, D-C, D-D interaction isgiven by [check this]
4πϵ0Vtot =∑i>j
qiTijqj − qi
∑α
Tαijµ
indj,α +
∑α
µindi,αTα
ijqj −∑α,β
µindi,αTαβ
ij µindj,β
(34)
6
From Bottcher’s book (Theory off Electric Polarization vol 1 Dielectrics in static field 2nd ed. p110), the energy Uof the induced dipole system
U = −µindi · E(ri) + Upol (35)
where Upol is the work of polarization. At equilibrium, the energy will be minimal with an infinitesimal change ofinduced moment
dU = 0, for all dµindi (36)
Then we have
dUpol = −d[−µindi · E(ri)] = E(ri) · dµind
i =µind
i
αi· dµind
i =1
2αid[(µind
i )2] (37)
The induced dipole is formed in a reversible process
Upol =∫
dUpol =1
2αi
∫ µiind
0
d[(µindi )2] =
12αi
(µindi )2 (38)
If we count all the cotribution
Upol =∑
i
12αi
µindi · µind
i (39)
The MD code with this polarization scheme is available, e.g. Amber or Lucretius.