Blaak Ronald, Hansen Jean-Pierre
Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom.
J Chem Phys. 2006 Apr 14;124(14):144714. doi: 10.1063/1.2185092.
We present extensive molecular dynamics simulation results for the structure and the static and dynamical responses of a droplet of 1000 soft spheres carrying extended dipoles and confined to spherical cavities of radii R=2.5, 3, and 4 nm embedded in a dielectric continuum of permittivity epsilon(')>or=1. The polarization of the external medium by the charge distribution inside the cavity is accounted for by appropriate image charges. We focus on the influence of the external permittivity epsilon(') on the static and dynamic properties of the confined fluid. The density profile and local orientational order parameter of the dipoles turn out to be remarkably insensitive to epsilon('). Permittivity profiles epsilon(r) inside the spherical cavity are calculated from a generalized Kirkwood formula. These profiles oscillate in phase with the density profiles and go to a "bulk" value epsilon(b) away from the confining surface; epsilon(b) is only weakly dependent on epsilon('), except for epsilon(')=1 (vacuum), and is strongly reduced compared to the permittivity of a uniform (bulk) fluid under comparable thermodynamic conditions. The dynamic relaxation of the total dipole moment of the sample is found to be strongly dependent on epsilon(') and to exhibit oscillatory behavior when epsilon(')=1; the relaxation is an order of magnitude faster than in the bulk. The complex frequency-dependent permittivity epsilon(omega) is sensitive to epsilon(') at low frequencies, and the zero-frequency limit epsilon(omega=0) is systematically lower than the bulk value epsilon(b) of the static permittivity.
我们展示了关于1000个携带扩展偶极子的软球液滴的结构以及静态和动态响应的广泛分子动力学模拟结果,该液滴被限制在半径R = 2.5、3和4纳米的球形腔内,嵌入介电常数ε(′)≥1的介电连续介质中。腔内电荷分布对外部介质的极化由适当的镜像电荷来考虑。我们关注外部介电常数ε(′)对受限流体静态和动态性质的影响。结果表明,偶极子的密度分布和局部取向序参量对ε(′)非常不敏感。根据广义柯克伍德公式计算球形腔内的介电常数分布ε(r)。这些分布与密度分布同相位振荡,并在远离限制表面处趋于“体相”值ε(b);ε(b)仅微弱地依赖于ε(′),除了ε(′)=1(真空)的情况,并且与在可比热力学条件下均匀(体相)流体的介电常数相比大幅降低。发现样品总偶极矩的动态弛豫强烈依赖于ε(′),并且当ε(′)=1时表现出振荡行为;弛豫比在体相中快一个数量级。复频率依赖介电常数ε(ω)在低频时对ε(′)敏感,并且零频率极限ε(ω = 0)系统地低于静态介电常数的体相值ε(b)。