计算液态混合物内能和内压的统计热力学理论

假定二元液态混合物分子间的相互作用势能可以表示成多体相互作用势能的和,分子间的力为短程力,相互作用势能只与分子间的相对距离有关.利用分布函数理论导出了二元液态混合物的过剩内能和内压的公式.二元液态混合物的过剩内能和内压可以表示成体积的幂级数形式,其中的系数可以用多体相互作用势和多体径向分布函数表出.讨论了单元液体的内压和过剩内能的表达式,在两种特殊情形下,过剩内能和内压的表达式分别与Egelstaff的微扰论结果及Frank的实验结果具有相同的形式.讨论了二元混合物内压和内能的两个特例,其一,在特殊情形下,给出了混合液体过剩内能的混合规则的一个证明.其二,给出的二元混合物的过剩内能和内压的表达式与Frank的实验结果具有相同的形式.     

Molecular dynamics simulation of potassium along the liquid-vapor coexistence curve

The applicability of pair potential functions to liquid alkali metals is questionable. On the one hand, some recent reports in the literature suggest the validity of two-parameter pair-wise additive Lennard-Jones (LJ) potentials for liquid alkali metals. On the other hand, there are some reports suggesting the inaccuracy of pair potential functions for liquid metals. In this work, we have performed extensive molecular dynamics simulations of vapor-liquid phase equilibria in potassium to check the validity of the proposed LJ potentials and to improve their accuracy by changing the LJ exponents and taking into account the temperaturedependencies of the potential parameters. We have calculated the orthobaric liquid and vapor densities of potassium using LJ (12–6), LJ (8.5–4) and LJ (5–4), effective pair potential energy functions. The results show that using an LJ (8.5–4) potential energy function with temperature-independent parameters, ε and σ, is inadequate to account for the vapor-liquid coexistence properties of potassium. Taking into account the temperature-dependencies of the LJ parameters, ε(T) and σ(T), we obtained the densities of coexisting liquid and vapor potassium in a much better agreement with experimental data. Changing the magnitude of repulsive and attractive contributions to the potential energy function shows that a two-parameter LJ (5–4) potential can well reproduce the densities of liquid and vapor potassium. The results show that LJ (5–4) potential with temperature-dependent parameters produces the densities of liquid and vapor potassium more accurately, compared to the results obtained using LJ (12–6) and LJ (8.5–4) potential energy functions.     

A new anisotropic soft-core model for the simulation of liquid crystal mesophases

A new anisotropic soft-core model is presented, which is suitable for the rapid simulation of liquid crystal mesophases. The potential is based on a soft spherocylinder, which can be easily tuned to favor different liquid crystal mesophases. The soft-core nature of the potential makes it suitable for long-time step molecular dynamics or dissipative particle dynamics simulations, particularly as a reference model for mesogens or as an anisotropic solvent for use in combination with atomistic models. Results are presented for two variants of the new potential, which show different mesophase behaviors. Variants of the potential can also be linked together to produce more complicated molecular structures. Here, as an example, results are provided for a model multipedal liquid crystal, which has eight liquid crystalline groups linked to a central core via semiflexible chains. Here, despite the complexity of molecular structure, the model succeeds in showing the spontaneous formation of a liquid crystal phase. The results also demonstrate that there is a very strong coupling between the internal structure of the multipedal mesogen and the molecular order of the phase, with the mesogen spontaneously undergoing major structural rearrangement at the transition to the liquid crystal phase.     

A microfluidic electrochemical detection technique for assessing stability of thin films and emulsions

An experimental technique is developed for assessing stability of thin liquid films by application of electric potential to compress the liquid film and to simultaneously measure the electrical properties of the system. The concept involves creating a thin film at the intersection of two microchannels etched onto a glass substrate. A ramped DC potential difference is applied across the film, which develops an electrical stress across the film. Increasing the potential to a critical value leads to the rupture of the film. The critical potential is used to assess the stability of the liquid film. Small channel dimensions in this microfluidic platform allow characterization of thin films formed between micron-sized droplets representing systems with high capillary pressures, analysis of which are typically beyond the scope of conventional thin film characterization techniques. The results of DC potential breakdown of films show that critical potential can be considered as a measure of thin film stability.     

Molecular dynamics simulation of the water|nitrobenzene interface

The structure and dynamics of the neat water|nitrobenzene liquid|liquid interface are studied at 300 K using molecular dynamics computer simulations. The water is modeled using the flexible SPC potential, and the nitrobenzene is modeled using an empirically determined nitrobenzene potential energy function. Although nitrobenzene is a polar liquid with a large dielectric constant, the structure of the interface is similar to other water|non-polar organic liquid interfaces. Among the main structural features we describe are an enhancement of interfacial water hydrogen bonds, the specific orientation of water dipoles and nitrobenzene molecules, and a rough surface that is locally sharp. Surface roughness is also characterized dynamically. The dynamics of molecular reorientation are shown to be only mildly modified at the interface. The effect due to the polarizable many-body potential energy functions of both liquids is investigated and is found to affect only mildly the above results.     

Molecular simulation of the phase behaviour of ternary fluid mixtures: the effect of a third component on vapour–liquid and liquid–liquid coexistence

The Gibbs ensemble algorithm is implemented to determine the vapour–liquid and liquid–liquid phase coexistence of dilute ternary fluid mixtures interacting via a Lennard–Jones potential. Calculations are reported for mixtures with a third component characterised by different intermolecular potential energy parameters. Comparison with binary mixture data indicates that the choice of energy parameter for the third component affects the composition range of vapour–liquid substantially. The addition of a third component lowers the energy of liquid phase while slightly increasing the energy of the vapour phase.     

Atom-atomic potentials and the correlation distribution functions for modeling liquid benzene by the molecular dynamics methods

A brief review of the potential functions used in the molecular dynamics modeling of liquid benzene is presented. The structural characteristics of liquid benzene obtained from the correlation distribution functions are discussed. It is demonstrated that, within the framework of this approach, the predicted structure of liquid benzene is virtually independent of the form of the potential used.     

Sedimentation velocity and potential in a concentrated suspension of charged liquid drops

The sedimentation behavior of a concentrated suspension of charged liquid drops is analyzed theoretically at arbitrary surface potential and arbitrary double-layer thickness; that is, the effects of double-layer polarization and double-layer overlapping are taken into account. Kuwabara's unit cell model is employed to model the suspension system, and a pseudospectral method based on the Chebyshev polynomial is adopted to solve the governing electrokinetic equations numerically. Several interesting phenomena, which are of significant influence if the internal flow inside a liquid drop is taken into account, are observed. Key factors are examined such as the thickness of the electric double layer, the magnitude of the surface potential, the volume fraction of liquid drops, and the viscosity of the internal fluid. The results presented here add another dimension to the previous studies, which include concentrated suspensions of rigid particles and mercury drops under low zeta potential, with the consideration of the internal flow of liquid drops and double-layer polarization, characterized by its viscosity and the zeta potential respectively. It is found, among other things, that the smaller the viscosity of the internal fluid is, the higher the sedimentation velocity of liquid drops. The higher the zeta potential is, the larger the decrease in sedimentation velocity. In particular, the sedimentation velocity of an inviscid drop (gas bubble) is about three times higher than that of a rigid one. The decrease in sedimentation velocity resulting from the effect of double-layer polarization achieves about 50% if the zeta potential is sufficiently high.     

Liquid Marble Photosensor

A liquid marble is a liquid droplet coated by a hydrophobic powder. The liquid marble does not wet adjacent surfaces and therefore can be manipulated as a dry soft body. A Belousov-Zhabotinsky (BZ) reaction is an oscillatory chemical reaction exhibiting waves of oxidation. We demonstrate how to make a photo-sensor from BZ medium liquid marbles. We insert electrodes into a liquid marble, prepared from BZ solution and coated with polyethylene powder. The electrodes record a potential difference which oscillates due to oxidation wave-fronts crossing the electrodes. When the BZ marble is illuminated by a light source, the oxidation wave-fronts are hindered and, thus, the electrical potential recorded ceases to oscillate. We characterise several types of responses of BZ marble photosensors to various stimuli, and provide explanations of the recorded activity. BZ liquid marble photosensors may find applications in the fields of liquid electronics, soft robotics and unconventional computing.     

离子液体在多孔碳电极上的电化学性能

制备了数种离子液体及离子液体有机溶液电解质,用线性电位扫描法测试了它们的电化学窗口;并通过循环伏安、交流阻抗、电势阶跃等电化学方法,对它们在多孔碳电极中的电化学性能进行了研究.实验结果表明:溶剂对1-乙基-3-甲基咪唑类离子液体有机溶液电解质的电化学窗口有较大的影响.离子液体及离子液体有机溶液电解质在多孔碳电极上的电化学性能与其电导率密切相关,电导率越大,充电时间常数越小,比容量越大;但比容量降低的倍率远小于电导率降低的倍率.     

三维纳米液膜沉积结构的模拟研究

基于三维动力学蒙特卡罗模型,再现了纳米液膜干燥后形成的双尺度沉积结构,探究了液体化学势、纳米粒子移动速率、液体临界蒸发率以及化学势锐度等参数对纳米液膜沉积结构的影响。结果表明,纳米液膜干燥时,液膜中的纳米粒子随着三相线移动自组装形成各种各样的沉积结构沉积在基板上。纳米液膜蒸发过程中,当液体蒸发达到临界蒸发率时液体化学势将发生突变,使其形成两种沉积结构。随液体初始化学势的增大,沉积结构逐渐变为均匀分布的密集网状结构。液体临界蒸发率越小,液体化学势突变后形成的沉积结构越明显。纳米粒子的移动速率越快,沉积结构中的枝状结构越少。化学势锐度对双尺度沉积结构的差异有很大的影响,锐度越大,两种沉积结构的差异越大。     

Investigation of dynamic streaming potential by dimensional analysis

The theory of streaming potential at sinusoidal flow of liquid in a porous medium is a convenient and fruitful tool for determination of the interface properties of materials and also for construction of apparatus for zeta potential measurements and electrokinetic transducers. An investigation of the dynamic streaming potential by the method of dimensional analysis is presented. This method provides a wider approach to the problem under consideration. As a result, relationships between streaming potential in a porous medium and mechanical quantities are established. These quantities include pressure gradient in a liquid inside pores and capillaries, acceleration of capillaries, and the solid part of a porous medium, and the viscous friction force the liquid exerts on the solid part. The corresponding formulas for streaming potential are presented. The relationship between the streaming potential and viscous friction force does not depend on the frequency of oscillation and pore size. All these formulas in particular cases are transformed to known formulas for the streaming potential.     

Interfacial tension measurements at the interface between two highly immiscible electrolyte solutions: the trapped air bubble method

A dynamic method for the measurement of interfacial tension at a liquid–liquid interface under controlled interfacial potential difference is described. The interface was formed as a drop suspended at the tip of a liquid filled glass micro-syringe into which a trapped air bubble of known volume had been introduced. Changes in volume of the air bubble allow changes in the surface tension to be calculated. Application to measurement of the variation in interfacial tension with cyclic variation in interfacial potential difference (‘cyclic volttensiometry’) is demonstrated.     

Second order of perturbation theory correction to the radial distribution function of a liquid with the interaction potential of hard spheres plus a square-well

A liquid with the interaction potential of hard spheres plus a square-well is analyzed using the Monte-Carlo technique. Numerical results for the perturbation theory series over a square-well potential are obtained in the form of the Barker and Henderson discrete representation. Approximating expressions for the correction to a liquid radial distribution function in the second order of perturbation theory are presented. The obtained results allow us to define this correction with a root-mean-square deviation of about 0.007. It is shown that the given approach provides a complete calculation in the second order of perturbation theory, and also the determination of the third order correction to the free energy for a liquid interacting with the potential of the Lennard-Jones type.     

液体的内压和内能的统计热力学理论

利用分布函数理论导出了液体的内能和内压公式.液体的内压和过剩内能可以表示成体积的幂级数形式,其中的系数可以用多体相互作用势和多体径向分布函数表出,它们仅仅与温度有关.讨论了液体仅存在第n次多体相互作用势情形的内压和过剩内能的表达式,结果与Egelstaff的微扰理论结果具有相同的形式,不仅给出了相应参数的表达式而且适用于多体相互作用较强的情形.定义了物性参数α(T)和m,得到的液体过剩内能和内压的表达式与Frank实验结果具有相同的形式,其结果不仅给出了参数α(T)和m的表达式,而且指出了Frank的过剩内能和内压公式只适用于参数α(T)和m与体积无关的液体.     

液体接界电势公式分析

以不同浓度的硝酸银接触为例,讨论了液体接界电势产生的原因,阐明了正负离子迁移的吉布斯总能量变化与过剩电荷的吉布斯能量差异,推导了液体接界电势的计算公式,并说明了液体接界电势与电势正负侧的活度比关系。     

Validation of intermolecular pair potential model of SiH4: molecular-dynamics simulation for saturated liquid density and thermal transport properties

We demonstrate a validation of the intermolecular pair potential model of SiH(4), which is constructed from ab initio molecular-orbital calculations and expressed as the sum of the exponential and the London dispersion terms. The saturated liquid densities of SiH(4) are calculated for temperatures from 100 to 225 K by molecular-dynamics (MD) simulation. The average deviation between the experiment and the MD simulation using the present potential model is 3.9%, while the deviations exceed 10% for other well-known potential models such as the five-center Lennard-Jones (LJ) model. Subsequently, the shear viscosity, the thermal conductivity, and the self-diffusion coefficient of liquid SiH(4) are calculated by an equilibrium MD simulation with the Green-Kubo formula from 100 to 225 K. The average deviations from experiment are 11.8% and 13.7% for the shear viscosity and the thermal conductivity, respectively. Comparing the present model with an empirical one-center LJ model, it turns out that the rotational energy transfer through the intermolecular potential energy, which comes from the anisotropic potential energy, plays an important role in the thermal conductivity of liquid SiH(4). These results indicate that the present intermolecular potential model has an ability to give realistic pictures for liquid SiH(4) through molecular simulations.     

Temperature and density dependence of the structure of liquid lithium from an improved pseudopotential theory

Leribaux, H.R., and Lemarchand, J.L., 1978. Temperature and density dependence of the structure of liquid lithium from an improved pseudopotential theory. Fluid Phase Equilibria, 2: 79–90.The structure of liquid lithium and alloys has recently been measured more accurately near their melting point. We present for liquid lithium a modification of a previous pseudopotential for liquid metals, obtained by including explicitly the Xα exchange energy in the electron pseudopotential. The structure factor S(q), obtained from this model potential via the effective pair potential, and calculated by our solution of the Percus—Yevick equation, is found to have the right oscillations compared with the experimental structure factor, the first peak also having the right magnitude. A Percus—Yevick hard-sphere structure factor is also obtained via thermodynamic perturbation and is found to be surprisingly good even for lithium, except for the high-q oscillations. We present computed structures and radial distributions for several higher temperatures involved typically in the liquid lithium alloys, up to the highest temperatures at which saturated density is measured. Our model potential predicts also the correct cohesive energy for liquid lithium and liquid lithium—magnesium.     

DNA translocation through polyelectrolyte‐modified nanopores: An analytical approximation

An analytical model for the electrophoretic speed of DNA translocating through nanopore functionalized with polyelectrolyte (PE) brush is developed for the first time. The electrophoretic speed depends on DNA surface potential, applied electric field, viscosity, and permittivity of solution along with velocity and electrostatic potential at liquid–polyelectrolyte layer (PEL) interface where the interface seemed to behave similar to that of a solid‐state nanopore wall. Under the limit of Debye–Hückel linearization, the electrostatic potential at liquid–PEL interface and at DNA surface have been calculated. Velocity at liquid–PEL interface has been estimated by assuming a linear variation of hydrodynamic frictional force within the PEL. It is observed that velocity and electrostatic potential at liquid–PEL interface strongly depend on PE charge density and softness parameter. Present analytical results show excellent agreement with exact numerical results (i.e., without any approximation) at a higher salt concentration where Debye–Hückel linearization is applicable.