受限于纳米狭缝中的四缔合Lennard-Lones流体的相平衡和界面张力研究

分别用改进的基础测量理论和平均球近似理论表达短程作用和长程作用对四缔合Lennard-Jones流体的过剩自由能的贡献. 在密度函泛理论的框架下, 研究了平均密度等温线, 密度分布, 未缔合分子在平衡汽相和液相中的分布, 相平衡以及平衡时的界面张力等热力学性质. 分析了缔合能量, 流体-固体作用和孔宽对受限于纳米狭缝中的四缔合Lennard-Jones流体相行为的影响.     

一种对大体系局部进行高精度量子化学计算的新方法

提出一种可以对大体系的局部区域进行高精度量子化学计算的新方法.将体系划分为活性区和环境区,采用比较粗略的算法获得大体系的密度矩阵,通过基组的变换将密度矩阵分解为对应于活性区与环境区的子矩阵的直和,活性区和环境区分别拥有整数个电子,成为相对独立的化学键饱和的子体系.然后冻结环境区密度矩阵,同时考虑环境区在活性区所产生势场的作用,对活性区进行高精度的量子化学计算.对一些分子进行计算,得到的电荷分布、键长、键断裂能、电离势与电负性等与对大体系整体进行高精度计算得到的结果相当一致,而计算量减少很多,表明提出的方法可行,并有应用价值.     

用密度泛函理论研究Lennard-Jones 流体在狭缝中的相平衡

用改进的基础度量理论(modified fundamental measure theory, MFMT)和密度Taylor展开分别表达过剩自由能中的短程作用和色散作用. 流体分子与狭缝壁之间的相互作用以10-4-3势能函数表达. 由巨势最小原理确定Lennard-Jones (LJ)流体在狭缝中的密度分布和过剩吸附量, 所得结果与分子模拟数据吻合良好. 根据平衡时两相温度, 化学势及巨势相等, 计算了LJ流体在狭缝中的相平衡.     

用密度泛函理论研究Lennard-Jones 流体在狭缝中的相平衡

用改进的基础度量理论(modified fundamental measure theory, MFMT)和密度Taylor展开分别表达过剩自由能中的短程作用和色散作用. 流体分子与狭缝壁之间的相互作用以10-4-3势能函数表达. 由巨势最小原理确定Lennard-Jones (LJ)流体在狭缝中的密度分布和过剩吸附量, 所得结果与分子模拟数据吻合良好. 根据平衡时两相温度, 化学势及巨势相等, 计算了LJ流体在狭缝中的相平衡.     

含时密度矩阵重正化群的理论与应用

密度矩阵重正化群(DMRG)作为计算低维强关联体系强有力的方法为人熟知,在量子化学电子结构计算中得到广泛应用.最近几年,含时密度矩阵重正化群(TD-DMRG)的理论取得较快发展, TD-DMRG逐渐成为复杂体系量子动力学理论模拟的重要新兴方法之一.本文综述了基于矩阵乘积态(MPS)和矩阵乘积算符(MPO)的DMRG基本理论,并重点介绍了若干最常见的TD-DMRG时间演化算法,包括基于演化再压缩(P&C)的算法、基于含时变分原理(TDVP)的算法和时间步瞄准(TST)算法;还对利用TD-DMRG模拟有限温体系的纯化(Purification)算法和最小纠缠典型量子热态(METTS)算法进行了介绍.最后,对近年TD-DMRG在复杂体系量子动力学中的应用进行了总结.     

过渡金属络合物荷移光谱的研究

在配位场近似下,用不可约张量方法对具有O_4对称性的FeCl_6~(3-)及与其结构相似的过渡金属络合物的Hamilton矩阵进行了矩阵分割,得到了一种计算荷移光谱的简单方法.在计算中采用方案,考虑了纯组态内部的Coulomb作用和旋轨作用,忽略了不同组态间的Coulomb作用.并采取冻结轨道近似.计算结果与实验光谱数据符合较好.     

N2,O2水溶液光谱的分子动力学模拟

用分子动力学模拟方法研究了N2和O2水溶液的光谱性质．给出了能描述分子内部运动的溶质－溶剂相互作用势．对溶质和溶剂原子的速度自相关函数（VACF）作了计算．讨论了所得VACF的性质并计算了其谱密度．溶质分子振动谱出现的红移，与液态N2，O2的Raman实验结果相吻合．模拟得出的转动谱表明了溶剂分子对溶质转动运动的阻滞，模拟结果也表明VACF计算对溶液和液体光谱的研究十分有效．     

He-LiH体系转动非弹性碰撞的理论研究

在单双迭代耦合簇CCSD(T)势能面的基础上,运用密耦方法讨论了He-LiH体系的转动非弹性碰撞.计算结果表明,对LiH分子j=0→j′跃迁,跃迁截面主要由各向异性的短程相互作用和长程的“软”排斥共同作用的结果,未见明显的长程吸引势贡献.态－态跃迁总截面表现出振荡结构,长程“软”排斥分波只对j=0向j′=1、2、3的跃迁总截面有较大贡献,而j′≥ 4跃迁的积分截面则几乎由各向异性的短程部分贡献.     

熔融LiF、KCl和互易系LiF-KCl的分子动力学模拟计算

我们用分子动力学方法模拟了1200K的LiF、KCl和LiF-KCl熔盐系的结构和性质,包括热力学性质、结构动力学性质和输运性质,并对其瞬态图象进行了分析。计算中取Fumi-Tosi势为离子间有效对势,所用程序MDIONS由中科院理论物理研究所的CPC库提供。模拟系统采用立方元胞,内含216个离子,正负离子各半。互易系     

频域空间密度矩阵重正化群的研究进展

密度矩阵重正化群(DMRG)作为低维强关联体系中电子结构计算的强有力方法被广泛熟知,并被迅速地应用于量子化学,不仅在电子结构计算中发挥重要作用,同时也在近几年迅速地成为复杂体系量子动力学计算的重要方法.在DMRG框架中,衍生出了一系列计算动态响应性质的有效方法,并得到了广泛应用.本文简述了DMRG的基本理论,其矩阵乘积态(MPS)表示有效地扩展了该方法的应用范围.重点介绍了基于线性响应理论的动态DMRG,在频率空间求解系统在零温以及有限温度下响应性质的算法,并介绍其在电子关联问题和电子-声子关联问题中的应用,最后展望了该领域的未来发展方向.     

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.     

Description of self-diffusion coefficients of gases,liquids and fluids at high pressure based on molecular simulation data

The purpose of this work is to evaluate the potential of modeling the self-diffusion coefficient (SDC) of real fluids in all fluid states based on Lennard–Jones analytical relationships involving the SDC, the temperature, the density and the pressure. For that, we generated an equation of state (EOS) that interrelates the self-diffusion coefficient, the temperature and the density of the Lennard–Jones (LJ) fluid. We fit the parameters of such LJ–SDC–EOS using recent wide ranging molecular simulation data for the LJ fluid. We also used in this work a LJ pressure–density–temperature EOS that we combined with the LJ–SDC–EOS to make possible the calculation of LJ–SDC values from given temperature and pressure. Both EOSs are written in terms of LJ dimensionless variables, which are defined in terms of the LJ parameters ɛ and σ. These parameters are meaningful at molecular level. By combining both EOSs, we generated LJ corresponding states charts which make possible to conclude that the LJ fluid captures the observed behavioral patterns of the self-diffusion coefficient of real fluids over a wide range of conditions. In this work, we also performed predictions of the SDC of real fluids in all fluid states. For that, we assumed that a given real fluid behaves as a Lennard–Jones fluid which exactly matches the experimental critical temperature T_c and the experimental critical pressure P_c of the real fluid. Such an assumption implies average true prediction errors of the order of 10% for vapors, light supercritical fluids, some dense supercritical fluids and some liquids. These results make possible to conclude that it is worthwhile to use the LJ fluid reference as a basis to model the self-diffusion coefficient of real fluids, over a wide range of conditions, without resorting to non-LJ correlations for the density–temperature–pressure relationship. The database considered here contains more than 1000 experimental data points. The database practical reduced temperature range is from 0.53 to 2.4, and the practical reduced pressure range is from 0 to 68.4.     

Computer simulation of condensed gas mixtures

Molecular dynamics simulations have been performed for an equimolar mixture of oxygen and nitrogen using a diatomic Lennard—Jones model potential to obtain excess enthalpy and volume and to study the effect of attractive forces on the microscopic structure in the liquid. Simulations have also been performed to study the effect of mixing rules on excess properties in equimolar binary mixtures of atomic and diatomic liquids and to test the existing perturbation theories.     

Thermal diffusion in Lennard–Jones fluids in the frame of the law of the corresponding states

This work is related to the definition of a reduced thermal diffusion coefficient thanks to numerical microscale molecular dynamics simulations. This cross transport process, also called Soret effect, couples mass flux and thermal gradient and is still largely misunderstood. For this study, we have applied a boundary driven non-equilibrium molecular dynamics algorithm on Lennard–Jones spheres mixtures. Simulations have been performed at a constant reduced supercritical state, using a van der Waals’ one fluid approximation in order to fulfil the law of the corresponding states. In binary mixtures, we have studied the molecular parameters and the molar fraction influences on thermal diffusion separately and then combined. It is shown that, on pressure and on thermal conductivity, the corresponding states law is fulfilled for a wide range of molecular parameters ratios. In this frame, we have then constructed simple correlations which relate thermal diffusion factor to the mixture parameters. Combining the relations obtained, a reduced thermal diffusion factor taking into account all the various contributions has been defined. Finally, it is shown that this relation enables us to estimate thermal diffusion in various binary and ternary mixtures of Lennard–Jones spheres representing alkanes with a maximum deviation of 15%.     

Theoretical study of tetrahydrofuran: Comparative investigation of spectroscopic and structural properties between gas and liquid phases

Rovibrational spectroscopic constant of tetrahydrofuran (THF) dimer have been calculated starting from three potential energy curves, each one obtained in a different way: (i) by ab initio calculations at MP2/aug‐cc‐pVDZ level; (ii) using Lennard‐Jones liquid parameters available in the literature, and (iii) from the pair obtained through Monte Carlo Simulation of liquid THF. The comparison among these results allowed the characterization of many solvent effect contributions. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Stability of Metastable Phases and Kinetics of Nucleation in a Simple Single-Component System (Molecular Dynamics Simulation) (A Review)

Russian Journal of General Chemistry - The review presents the results of molecular dynamics simulation of metastable states in the systems of the Lennard–Jones particles. The boundaries of...     

Exact solution of the Schrödinger equation with a Lennard–Jones potential

The Schrödinger equation with a Lennard–Jones potential is solved by using a procedure that treats in a rigorous way the irregular singularities at the origin and at infinity. Global solutions are obtained thanks to the computation of the connection factors between Floquet and Thomé solutions. The energies of the bound states result as zeros of a function defined by a convergent series whose successive terms are calculated by means of recurrence relations. The procedure gives also the wave functions expressed either as a linear combination of two Laurent expansions, at moderate distances, or as an asymptotic expansion, near the singular points. A table of the critical intensities of the potential, for which a new bound state (of zero energy) appears, is also given.     

Viscosity Prediction for Oxygen,Nitrogen and Their Mixtures at Zero and Moderately Dense Regimes via Semi‐Empirically Based Assessment

The viscosity coefficients for the gaseous states of N₂ and O₂ and their mixtures are determined at zero and moderately density regimes. The Lennard‐Jones 12–6 (LJ 12–6) potential energy function is used as the initial model potential required y the technique. The interaction potential energies from the inversion procedure reproduce the viscosity commensurate to the best measurements. The initial density dependence of gaseous viscosity coefficient according to the Rainwater‐Friend theory, which was given by Najafi et al., has been considered for pure N₂ and pure O₂.     

Bohr Sommerfeld quantisation and molecular potentials

We combine, within the Bohr Sommerfeld quantization rule, a systematic perturbation with asymptotic analysis of the action integral for potentials which support a finite number of bound states with E < 0 to obtain an interpolation formula for the energy eigenvalues. We find interpolation formulae for the Morse potential as well as potentials of the form \({V=V_0 \left[ {\left( {\frac{a}{x}} \right)^{2k}-\left( {\frac{a}{x}}\right)^{k}} \right]}\). For k = 6 i.e. the well known Lennard Jones potential this yields results within 1 per cent of the highly accurate numerical values. For the Morse potential this procedure yields the exact answer. We find that the result for the Morse potential which approaches zero exponentially is the \({k\rightarrow\infty}\) limit of the Lennard Jones class of potentials.