A polarizable ion model for the structure of molten AgI

The results are reported of the molecular dynamics simulations of the coherent static structure factor of molten AgI at 923 K using a polarizable ion model. This model is based on a rigid ion potential, to which the many body interactions due to the anions induced polarization are added. The calculated structure factor is in better agreement with recent neutron diffraction data than that obtained by using simple rigid ion pair potentials. The Voronoi-Delaunay method has been applied to study the relationship between voids in the spatial distribution of cations and the prepeak of the structure factor.     

Molecular dynamics study of polarization effects on AgI

Three different models of AgI are studied by molecular dynamics simulations. The first one is the rigid ion model (RIM) with the effective pair potential of the Vashishta and Rahman form and the parametrization proposed by Shimojo and Kobayashi. The other two are polarizable ion models in which the induced polarization effects have been added to the RIM effective pair potential. In one of them (PIM1), only the anions are assumed to be polarizable by the local electric field. In the other one (PIM2s), the silver polarization is also included, and a short-range overlap-induced polarization opposes the electrically induced dipole moments. This short-range polarization is proved to be necessary to avoid overpolarization when both species are assumed to be polarizable. The three models reproduce the superionic character of alpha-AgI at 573 K and the liquid behavior of molten AgI at 923 K. The averaged spatial distribution of the cations in the alpha-phase obtained for PIM1 appears to be in better agreement with experimental data analysis. The PIM1 also reproduces the structure factor prepeak at about 1 A(-1) observed from neutron diffraction data of molten AgI. The three models retain in the liquid phase the superionic character of alpha-AgI, as the mobility of the cations is significantly larger than that for the anions. The ionic conductivity for the polarizable ion models is in better agreement with experimental data for alpha-AgI and molten AgI.     

The structure of molten AgCl, AgI and their eutectic mixture as studied by molecular dynamics simulations of polarizable ion model potentials

The structure of molten AgCl, AgI, and their eutectic mixture Ag(Cl(0.43)I(0.57)) is studied by means of molecular dynamics simulations of polarizable ion model potentials. The corresponding static coherent structure factors reproduce quite well the available neutron scattering data. The qualitative behavior of the simulated partial structure factors and radial distribution functions for molten AgCl and AgI is that predicted by the reverse Monte Carlo modeling of the experimental data. The AgI results are also in qualitative agreement with those calculated from ab initio molecular dynamics.     

A polarizable ion model for the structure of molten CuI

The results are reported of the molecular dynamics simulations of the coherent static structure factor of molten CuI at 938 K using a polarizable ion model. This model is based on a rigid ion potential to which the many body interactions due to the anions induced polarization are added. The calculated structure factor reproduces the clear sharp prepeak observed in neutron diffraction data. The corresponding partial structure factors and the related radial distribution functions calculated by molecular dynamics are compared with those found in the literature derived from a combination of neutron and x-ray diffraction data with the aid of the reverse Monte Carlo simulation technique, as well as those calculated by ab initio MD simulations.     

Molecular dynamics study of polarizable ion models for molten AgBr

Three different polarizable ion models for molten AgBr have been studied by molecular dynamics simulations. The three models are based on a rigid ion model (RIM) with a pair potential of the type proposed by Vashishta and Rahman for alpha-AgI, to which the induced dipole polarization of the ions is added. In the first (PIM1) the dipole moments are only induced by the local electric field, while in the other two (PIM1s and PIM2s) a short-range overlap induced polarization opposes the electrically induced dipole moments. In the PIM1 and the PIM1s only the anions are assumed polarizable, while in the PIM2s both species are polarizable. Long molecular dynamics simulations show that the PIM2s is an unphysical model since, for some improbable but possible critical configurations, the ions become infinitely polarized. The results of using the PIM1, the PIM1s, as well as those of the simple RIM, have been compared for the static structure and ionic transport properties. The PIM1 reproduces the broad main peak of the total structure factor present in the neutron diffraction data, although the smoothed three-peak feature of this broad peak is slightly overestimated. The structural results for the PIM1s are intermediate between those for the RIM and the PIM1, but fail to reproduce the experimental features within the broad principal peak. Concerning the ionic transport properties, the value of the conductivity obtained using PIM1 is in good agreement with experimental values, while the self-diffusion coefficients and the conductivity for the PIM1s are lower than the corresponding values using the PIM1 and the RIM.     

异质结型AgI/AgBr/TiO2催化剂的制备及其光催化性能

首先以沉积-沉淀法制备AgBr/TiO₂复合催化剂,然后采用离子交换法制备出新型的异质结型AgI/AgBr/TiO₂光催化剂.利用XRD和UV-Vis对AgI/AgBr/TiO₂光催化剂进行了表征.以甲基橙为染料模型,在可见光条件下(500 W、λ>420 nm)研究了AgI的含量对AgI/AgBr/TiO₂催化活性的影响.结果表明,AgI拓展了催化剂的吸收光谱范围;AgI生成量为AgBr的5%时,AgI/AgBr/TiO₂的催化活性最高.AgI/AgBr异质结的形成有利于光生电子和空穴的分离,提高AgI/AgBr/TiO₂的催化活性.     

Internal mobilities and diffusion in an ionic liquid mixture

The internal mobility gives the rate at which one ionic species moves relative to the other species present in an ionic mixture, it mirrors the differential strength of the interactions between different ionic species. In this work we examine the dependence of the internal mobilities of the Li(+) and K(+) ions on the composition in molten mixtures of LiF and KF. We compare them to the behaviour of the individual diffusion coefficients and the self-exchange velocities, which measure the rate at which an ion separates from its nearest-neighbour coordination shell. The examination is made using molecular dynamics simulations with polarizable, first-principles parameterised interaction potentials which are shown to reproduce the limited available experimental data on the transport properties of these mixtures extremely well. The results confirm that the composition-dependence of the internal mobilities in LiF/KF follows the unusual type-II behaviour, which is not reflected in that of the diffusion coefficients or the self-exchange velocities.     

Molecular dynamics study of polarizable point dipole models for molten sodium iodide

The structure, the ionic transport properties, and the dynamics of long-wavelength charge-density fluctuations, for two polarizable point dipole models of molten NaI, have been studied by molecular dynamics simulations. These models are based on a rigid ion potential to which the induced dipole polarization of the anions is added. The polarization is added in such a way that point dipoles are induced on the anions by both local electric field and short-range damping interactions that oppose the electrically induced dipole moments. The two polarizable ion models differ only in the range of the damping polarization interactions. The influence of the induced anion polarization on the different properties of simulated molten NaI is discussed.     

Equilibrium and nonequilibrium molecular dynamics simulations of the thermal conductivity of molten alkali halides

The thermal conductivity of molten NaCl and KCl was calculated through the Evans-Gillan nonequilibrium molecular dynamics (NEMD) algorithm and Green-Kubo equilibrium molecular dynamics (EMD) simulations. The EMD simulations were performed for a "binary" ionic mixture and the NEMD simulations assumed a pure system for reasons discussed in this work. The cross thermoelectric coefficient obtained from Green-Kubo EMD simulations is discussed in terms of the homogeneous thermoelectric power or Seebeck coefficient of these materials. The thermal conductivity obtained from NEMD simulations is found to be in very good agreement with that obtained through Green-Kubo EMD simulations for a binary ionic mixture. This result points to a possible cancellation between the neglected "partial enthalpy" contribution to the heat flux associated with the interdiffusion of one species through the other and that part of the thermal conductivity related to the coupled fluxes of charge and heat in "binary" ionic mixtures.     

熔融ZnCl2结构的分子动力学模拟研究

熔融ZnCl2作为一种离子性共价性参半的典型熔盐, 其近邻结构在实验测量和分子动力学模拟方面均作过一些研究。本文依据新近EXAFS实验结果, 比较了不同的有效势下模拟得到的径向分布函数,表明KDR势可作为一种实用势。并进一步在KDR势模拟产生的瞬态构型基础上, 使用键序参数方法研究了晶态和熔融态ZnCl2中的近邻结构。结果表明, 和晶态ZnCl2一样, 在熔融ZnCl2中存在稳定的Zn/Cl正四面体结构, 但熔态和晶态Zn/Cl近邻结构热波方差σ不同。计算表明300K晶态σ=5.0℃, 613K熔融态σ=12.2℃。也对熔融ZnCl2的网络状结构和宏观输运性质进行了讨论。     

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.     

First principles molecular dynamics of molten NaI: structure, self-diffusion, polarization effects, and charge transfer

The structure and self-diffusion of NaI and NaCl at temperatures close to their melting points are studied by first principles Hellmann-Feynman molecular dynamics (HFMD). The results are compared with classical MD using rigid-ion (RI) and shell-model (ShM) interionic potentials. HFMD for NaCl was reported before at a higher temperature [N. Galamba and B. J. Costa Cabral, J. Chem. Phys. 126, 124502 (2007)]. The main differences between the structures predicted by HFMD and RI MD for NaI concern the cation-cation and the anion-cation pair correlation functions. A ShM which allows only for the polarization of I- reproduces the main features of the HFMD structure of NaI. The inclusion of polarization effects for both ionic species leads to a more structured ionic liquid, although a good agreement with HFMD is also observed. HFMD Green-Kubo self-diffusion coefficients are larger than those obtained from RI and ShM simulations. A qualitative study of charge transfer in molten NaI and NaCl was also carried out with the Hirshfeld charge partitioning method. Charge transfer in molten NaI is comparable to that in NaCl, and results for NaCl at two temperatures support the view that the magnitude of charge transfer is weakly state dependent for ionic systems. Finally, Hirshfeld charge distributions indicate that differences between RI and HFMD results are mainly related to polarization effects, while the influence of charge transfer fluctuations is minimal for these systems.     

硝酸钾电解质溶液水化结构的分子动力学模拟

Molecular dynamics simulations were carried out to study the structure of ion clusters and hydration properties of KNO3 solution. The water molecule was treated as a simple-point-charge (SPC) model, and a four-site model for the nitrate ion was adopted. Both the Coulomb and Lennard-Jones interactions between all the charged sites were considered, and the long-range Coulomb electrostatic interaction was treated using Ewald summation techniques. The configuration of ionic pairs, the radial distribution function of the solution, and the effect of solution concentration on ionic hydration were studied in detail. It was found that there are ionic association phenomena in KNO3 solution and that the dimeric, triplet, solvent-separated ion pairs, and other complex clusters can be observed at high ionic concentration condition. As the concentration of solution decreases, the ionic hydration number increases, 5-7 for cation K+ and 3.5-4.7 for anion NO3-, which is in good agreement with former Monte Carlo and time-of-flight neutron diffraction results.     

The hall effect in molten salts

Abstract

The Hall Effect is measured in the molten ionic salts CuCl, AgCl and AgBr. Values of 0.023, 0.031 and 0.027 cm²/volt sec are obtained for the effective mobility. The values obtained are one or two orders of magnitude larger than drift mobilities and so indicate there can be no correspondence between Hall and drift mobilities in molten salts. These measurements appear to be the first of the Hall Effect in molten ionic salts     

LiCl熔盐急冷形成非晶固体的分子动力学计算机模拟研究

卤化物玻璃目前已成为引人注目的光纤新材料,用分子动力学方法研究液态急冷形成非晶态的过程,对于卤化物玻璃的形成过程研究也应是有用的。鉴于碱金属卤化物是最简单的熔盐,其动态结构亦很清楚。用分子动力学方法研究其急冷以形成玻     

Raman spectra of ionic liquids: a simulation study of LaCl3 and its mixtures with alkali chlorides

Theoretical Raman spectra of the elpasolite-structured crystal Cs2NaLaCl6 and of molten mixtures of LaCl6 with NaCl and CsCl have been obtained from computer simulations in order to examine how the Raman spectra reflect the coordination structure around the La3+ ions. This system is a model for many other trivalent metal halides and for examining how the network structure of the pure melts is broken down by the addition of alkali halides with different structure-breaking powers. The results suggest a way of reconciling the conclusions of Raman studies about the structures of the melts with those of neutron and x-ray-diffraction studies, which have already been examined with the same simulation methods. The Raman spectra, both polarized and depolarized, are calculated from a model for the dependence of the polarizability of the system on the ionic coordinates which was inspired by electronic structure calculations of the polarizabilities of ions in a condensed phase environment. Some results on the lifetimes of the coordination complexes responsible for the appearance of the discrete Raman bands are discussed.     

The sulfur dioxide-1-butyl-3-methylimidazolium bromide interaction: drastic changes in structural and physical properties

The contact between minor amounts of SO2 and crystalline 1-butyl-3-methyl-imidazolium bromide (BMIBr) causes the almost immediate melting of the ionic liquid (mp 45 degrees C) as well as a dramatic decrease in its viscosity in comparison to the pure molten phase. The same behavior was observed for other halide ionic liquids of higher melting points (70, 135, and 220 degrees C). The Raman spectrum of BMIBr-SO2 clearly indicates a specific charge transfer interaction involving SO2 and the halide. The measurements of ionic conductivity and diffusion coefficients obtained for the neat BMIBr (molten phase) and for the BMIBr-SO2 strongly suggest a higher degree of ionic association in the presence of SO2. Molecular dynamic simulations indicate that although the cation-anion distance is preserved in the short range, there is a variation in the interionic distances in the second shell, leading to a less organized structure in the long range. The modulation of the structural and physical properties of ionic liquids by SO2 and the convenient choice of the ions for gas absorption are suggested.     

掺杂正八面体AgBr乳剂介电频谱的研究——I-离子的影响

本文研究了用不同量的I~-(1×10~(-3)-4×10~(-2)mol I~-/mol AgBr)进行表面掺杂的正八面体AgBr乳剂的介电吸收频谱, 并用强X射线光源相应作了多晶X射线物相分析。当I~-的加入量小于1×10~(-2)mol/mol AgBr时, 介电吸收峰随加入I~-量的增加而逐渐向高频方向位移。吸收峰弥散, 分布很宽。当I~-的加入量大于1×10~(-2)mol/mol AgBr时, 样品的介电吸收峰不再明显向高频方向移动, 峰形亦相对比较尖锐。已有证据表明, 在掺杂I~-以后, 在AgBr微晶的表明层内除生成Ag(Br, I)混晶外, 还附生有一薄层β-AgI。分散相颗粒表面层大约2.0 nm范围内的组成与结构对非均匀电介质中的界面极化效应有重要影响, 此外, 对界面极化效应应用的可能性作了初步探讨。     

Effective interaction potentials for alkali and alkaline earth metal ions in SPC/E water and prediction of mean ion activity coefficients

The potential of mean force (PMF) acting between two simple ions surrounded by SPC/E water have been determined by molecular dynamics (MD) simulations using a spherical cavity approach. Such effective ion-ion potentials were obtained for Me-Me, Me-Cl-, and Cl(-)-Cl- pairs, where Me is a Li+, Na+, K+, Mg2+, Ca2+, Sr2+, and Ba2+ cation. The ionic sizes estimated from the effective potentials are not pairwise additive, a feature in the frequently used primitive model for electrolytes. The effective potentials were used in Monte Carlo (MC) simulations with implicit water to calculate mean ion activity coefficients of LiCl, NaCl, KCl, MgCl2, CaCl2, SrCl2, and BaCl2. Predicted activities were compared with experimental ones in the electrolyte concentration range 0.1-1 M. A qualitative agreement for LiCl and a satisfactory agreement for NaCl were found, whereas the predictions for KCl by two K+ models were less coherent. In the case of alkaline earth metal ions, all experimental activities were successfully reproduced at c = 0.1 M. However, at higher concentrations, similar deviations occurred for all divalent cations, suggesting that the dependence of the permittivity on the salt concentration and the polarization deficiency arising from the ordering of water molecules in the ion hydration shells are important in such systems.