单分子层受限冰-水共存界线毛细波与固-液相变动力学

冰-水界面动力学性质在冰形核、生长、表界面熔化中扮演核心角色,长期以来一直被广泛关注.然而,受限水体系中冰-水界面的动力学性质却鲜有报道.本工作利用平衡态分子动力学模拟方法和受限固-液两相平衡模拟技术,对两种水模型(恒定偶极矩、可极化)描述的单分子层受限冰-水两相平衡体系中的一维固-液界线开展研究.通过对一维受限冰-水...     

纳米受限水的研究进展

水是生命之源,在人类生存和社会生产中扮演了极其重要的角色,然而水的反常性质及在物理、化学、生物过程等领域中的作用和机理却仍存在很多谜团。近年来,水科学研究已逐渐成为科学研究的热点之一。地球上的水大部分是体相水,但在自然界和科学研究中,水同样会以界面/受限水的形式参与到物理、化学过程中。纳米受限水普遍存在于自然及合成的纳米环境中,受限水与体相水的差异主要体现在水的动力学及热力学性质的改变,受限水的存在对材料在生物、环境、地质和传感器等领域的应用也具有深远的影响。本文对纳米受限水的结构进行分析,并归纳了纳米受限水的动力学、热力学以及电学特性,对纳米受限水的研究手段及发展历程进行分类总结,举例介绍了纳米受限水在环境和能源等领域的潜在应用,最后对受限水研究进展及存在问题进行了总结,并对其后续发展进行展望。     

浓差极化的介电模型——复合膜/溶液体系的数值模拟

提出了具有电导率和介电常数线性分布的介质的介电模型, 并导出了其内部电的和结构性质的参数与宏观测量的电容和电导之间定量关系的理论表达式, 以模拟复合膜中的多孔层部分的介电弛豫行为. 大量的模拟计算描述并解释了多孔层介电谱随介电常数分布、厚度等性质而变化的规律. 进一步对具有层状构造的复合膜以及复合膜和溶液相组成的多层体系的弛豫行为进行了数值模拟, 比较了三个体系(多孔层、复合膜、复合膜/液相层状体系)的介电谱, 结果揭示了介电谱对各层性质的依赖关系. 所提出的电导率和介电常数线性分布的多孔层的介电模型, 也可用于具有其他电导率、介电常数分布规律的体系.     

石墨烯狭缝受限孔道中水分子的分子动力学模拟

石墨烯是一种具有广泛应用前景的纳米材料,特别是由石墨烯片层自组装形成的二维纳米通道能够应用于物质的过滤分离.本文采用分子动力学模拟方法研究了原态石墨烯/羟基改性石墨烯狭缝孔道中水分子的微观行为,模拟计算了水的界面结构性质和扩散动力学性质,所研究的石墨烯孔宽为0.6-1.5 nm.模拟结果表明,在石墨烯狭缝孔道中,水分子受限结构呈现层状分布,在超微石墨烯孔道(0.6-0.8 nm)中水分子可形成特殊的环状有序结构,石墨烯表面可诱导产生特殊的水分子界面取向.在石墨烯孔道中,水分子的扩散运动低于主体相水分子的扩散运动,羟基化石墨烯孔道可以促使水分子的扩散能力降低.对于改性石墨烯狭缝孔道,由于羟基的作用,水分子可以自发渗入0.6 nm的石墨烯孔道内.模拟所得到的受限水分子的动力学性质与水分子在石墨烯孔道内的有序结构有关.本文研究结果将有助于分析理解水分子通过石墨烯纳米通道的渗透机理,为设计基于石墨烯的纳米膜提供理论指导.     

双层氧化石墨烯纳米体系中受限水的介电常数

为了深入了解孔径和氧化程度对介电性能的影响, 通过分子动力学模拟的方法, 探究了水在不同孔径和氧化程度的双层石墨烯纳米通道中的介电行为. 实验结果表明, 在较窄的通道环境中, 水分子会表现出更强的有序取向, 受限水的介电常数随纳米通道空间的减小而减小. 随着氧化程度的增加, 较宽的层间距对介电常数的影响大于较窄的纳米通道. 对于最宽的通道(d=1.2 nm), 石墨烯双分子层内水的介电常数随氧化程度的增加而降低, 而对于相对较窄的通道(d=0.6, 0.9 nm), 介电行为呈现出非单调的趋势. 为了理解其背后的物理原因, 分别计算了3个纳米通道的氢键数量. 结果表明, 氢键数量以及动态稳定性（对应最快的衰减率）最低的是1.2 nm宽的纳米通道, 这表明水分子在大的纳米通道中更不稳定, 同时也比在0.6和0.9 nm的纳米通道中更无序.     

有机液体的慢极化介电谱

高绝缘液体存在一种特殊类型的极化响应; 其性质很不同于普通的介电极化响应。后者可用频域方法描述, 称为快响应;前者不能用频域而只能用时域方法描述, 故区别地称为慢响应。在微分时域介电谱中, 每种快或慢极化机构贡献一个峰。苯的时域谱有一个快峰和两个慢峰。室温下石蜡只有一个快峰和一个慢峰; 在加热熔解过程中慢峰分裂为两个峰。慢极化响应和液体中分子的局部束缚空间电荷有关; 慢极化介电谱可以给出分子动力的许多信息。     

压力驱动膜过程中浓差极化的介电模型

提出了电导率呈指数型分布的浓差极化层的介电模型. 利用传递方程研究了压力驱动膜过程的介电谱, 并导出层内电导率分布与宏观测量的电容和电导之间定量关系的理论表达式, 得到介电谱高低频特征值与浓差极化模数之间的关系. 计算结果表明, 较之线性分布电导率, 指数型分布电导率能更好地描述浓差极化层的介电谱. 研究表明, 介电谱及其特征值和弛豫强度随浓差极化模数和层厚度显著变化. 对由分离膜、主体溶液和浓差极化层组成的复合体系的介电谱进行了数值模拟, 发现电导较之电容其频率依存关系对浓差极化模数的变化更为敏感, 上述研究为利用介电谱在线监测膜分离过程提供了有价值的理论参考.     

模型分子筛对水分子的吸附

用PM3半经验量子化学方法研究了模型分子筛对水分子吸附的性质.通过模拟计算确定了两个低能吸附位点及相应的吸附热,结果与实验值非常接近.同时,从两个不同低能吸附位点可以看出,大笼内的吸附流动性比β笼内的流动性好,即对水分子的吸附有较宽的区域.     

模型分子筛对水分子的吸附

The adsorptive process of H2O molecule in model zeolite have been studied by PM3 semi empirical quantum chemistry method. The variation of heat of formation and O-H bond lengths with the adsorption distance have been calculated .The results show that there are two active adsorption site on the model zeolite. The average heats of adsorption is about -44 kJ•mol-1,in accordance with experiment results (-51 kJ•mol-1).     

膜/液界面浓度极化现象的介电解析

对由强、弱荷电膜和溶液构成的膜／液非均匀体系，在10－107Hz频率范围进行了介电测量，在直流偏压下，该体系显示了二个显著的介电弛豫．利用具有电导率分布相，即浓度极化相的介电理论对强行电的离子交换膜的结果进行了介电解析，从实验上测得的介电参数求出了反映膜／液界面浓度极化层构造的参数，讨论了该体系产生介电弛像的原因、说明了介电理论的合理性，并提出了膜／液界面体系产生介电弛豫现象的可能机理．     

限域单分子层冰-水两相平衡体系的分子动力学模拟

限域水因有极其丰富的结构物相变化而成为近年来水科学研究的一个热点.然而,不同相限域水之间的相平衡结构与性质却鲜有报道.论文提出一套分子动力学模拟技术,可实现纳米尺度限域条件下冰和水的不同结构相间形成的低维固-液界面（线）的平衡态模拟.应用此模拟技术,我们探索了0.65 nm限域尺寸、5000 bar限域压强条件下,单分子层厚度的冰-水（固-液）两相平衡,计算了该平衡体系一系列热力学量在界线附近的分布.平衡态的分子模拟结果直观地展示了粗糙型固-液界线的热毛细涨落、界线固-液结构转变的微观机制、以及缺陷在固-液相变区附近的形成与输运.各种热力学量分布函数呈现了二维限域冰-水共存界面（线）的特殊性质,如：相平衡区域的尺寸异于块体材料固-液界面,固-液界线处于切向压缩状态等.     

Hydration Dynamics of Water near an Amphiphilic Model Peptide at Low Hydration Levels: A Dielectric Relaxation Study

A dielectric relaxation study of aqueous solutions of the amphiphilic model peptide N‐acetyl‐leucine amide (NALA) at 298 K over a wide range of hydration levels is presented. The experiments range from states where water builds up several hydration layers to states where single water molecules or small water clusters are shared by several NALA molecules. The dielectric spectra reveal two modes on the 10 and 100 ps timescales. These are largely broadened with regard to the Lorentzian shape caused by simple Debye‐type relaxation, and are well described by the Kohlrausch–Williams–Watts stretched exponential function. The fast mode is assigned to water reorientation comprising bulk water as well as hydration water. Even when all water molecules are in contact with the solute, this fast component is dominant, and its mean relaxation time is retarded by less than a factor of two relative to neat water. The amplitude of the slow process is far higher than expected for the dipolar reorientation of the solute. The observations are consistent with results from molecular dynamics simulations for a similar model peptide reported in the literature. They suggest that the slow relaxation mode is mainly founded in peptide–water dipolar couplings, with some additional contribution from slowly reorienting hydration water molecules. The results are discussed with regard to the hydration dynamics of proteins and the interpretation of dielectric spectra of protein solutions.     

受限于不同螺旋性的纳米碳管中水的分子动力学模拟

近年来将纳米碱米碳管引入到与生命过程息息相关的离子通道膜的研究逐渐成 为热点，而其中的关键就是要了解受限于膜孔道（碳管）中水分子的行为。采用分 子动力学模拟在300 K和1.01 * 10~5 Pa下对受限于（6，6）armchair型和（10， 0）zigzag型纳米碳管中的水进行了研究，得到了水分子在碳管中的局部密度分布 等静态性质以及水分子在碳管中的传递等动态性质，并对不同势能模型的模拟结果 作了比较。结果表明选择不同的势能模型并没有改变此体系的固有性质，即水分子 不仅能够进入到憎水性的（6，6）碳管中而且能形成一条稳定的由氢键相连的纵列 （single file），而且在管中以纵列的形式进行同歇传递。此外，碳管螺旋性对 受限水的静态性质影响不大但对动态性质则有一定程度的影响，水分子在（10，0 ）zigzag型碳管中的传递能力要强于在（6，6）armchair型碳管中的能力。     

Molecular Simulation on Transport Properties of Confined Water

The diffusivity and viscosity of water confined in micropores were studied by molecular dynamics simulations. The effects of pore width and density were analyzed at pore widths from 0.9 to 2.6nm. The diffusivity in micropores is lower than that of the bulk, and it decreases as pore width decreases and as density increases. But the viscosity in micropores is much larger than that of the bulk, and it increases as pore width decreases and as density increases. The diffusivity in channel parallel direction is obviously larger than that in channel perpendicular directions.     

Structure and Thermodynamics of Mg:Phosphate Interactions in Water: A Simulation Study

The association of Mg²⁺ and H₂PO₄^? in water can give insights into Mg:phosphate interactions in general, which are very widespread, but for which experimental data is surprisingly sparse. It is studied through molecular dynamics simulations (>100 ns) by using the polarizable AMOEBA force field, and the association free energy is computed for the first time. Explicit consideration of outer‐sphere and two types of inner‐sphere association provides considerable insight into the dynamics and thermodynamics of ion pairing. After careful assessment of the computational approximations, the agreement with experimental values indicates that the methodology can be extended to other inorganic and biological Mg:phosphate interactions in solution.     

Water Confined in the Local Field of Ions

Interionic distances are shorter in concentrated ionic solutions, thus instigating the interaction and overlap of hydration shells, as ions become separated by only one or two layers of water molecules. The simultaneous interaction of water with two oppositely charged ions has, so far, only been investigated by computer simulation studies, because the isolated vibrational spectroscopic signature of these molecules remains undetected. Our combined near‐infrared spectroscopic and molecular dynamics simulation studies of alkali halide solutions present a distinct spectral feature, which is highly responsive to depletion of bulk water and merging of hydration shells. The analysis of this spectral feature demonstrates that absorption trends are in good agreement with the law of matching affinities, thus providing the first successful vibrational spectroscopic treatment of this topic. Combined with commonly observed near‐infrared bands, this feature provides a spectral pattern that describes some relevant aspects of ionic hydration.     

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl₂) solutions at different concentrations (0.6–2.8 mol kg⁻¹) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl₂ has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.     

Anisotropic Dynamics of Binary Particles in Confined Geometries

The diffusion dynamics of colloidal particles in a good solvent confined between two parallel quartz walls have been studied within the framework of dynamical density functional theory. The highly ordered density layers induced by interfacial effects give rise to the oscillating dynamics, resulting in position-dependent structural relaxations and diffusivities. Further investigation reveals that particle size, particle-wall interaction, and slitpore width play different roles in affecting the oscillating behaviors along different directions. As a result, the theory yields the local mean square displacements in perpendicular and parallel directions, which agree remarkably well with prior experimental measurements. The results indicate that the mean square displacements can be quantitatively predicted based on the knowledge of inhomogeneous thermodynamics and dynamics. The local and averaged free energy evolution of the binary particles has been described and presented to understand their dynamic mechanism in confined geometry.