在分子动力学模拟中计算库伦相互作用:Evjen方法的修订
1Department of Physics, Royal Holloway, University of London, Egham, Surrey TW20 0EX, United Kingdom.
The Journal of chemical physics
|January 7, 2026
概括
一种新的直接总和 (DS) 方法准确地计算了库伦相互作用在分子动力学模拟的散装和薄膜系统. 这种实体空间方法为复杂系统提供了对Ewald方法的直接替代方案.
科学领域:
- 计算化学的计算化学
- 材料科学 材料科学 材料科学
- 物理 物理学 物理
背景情况:
- 准确计算库伦相互作用对于分子动力学模拟至关重要.
- 像Ewald方法这样的现有方法在实施和适用性方面存在局限性.
- 埃维恩的方法提供了一个早期的真实空间方法,但需要扩展.
研究的目的:
- 引入一种新的直接总和 (DS) 方法来计算库伦相互作用.
- 将DS方法扩展到处理非立方单位细胞和板状系统.
- 评估DS方法的准确性和效率与Ewald方法相比.
主要方法:
- 直接总和 (DS) 方法是Evjen方法的延伸,具有单元细胞双极校正.
- 在单位电池基础上对截断的点电荷相互作用的实时空间总和.
- 将DS理论扩展到具有3D电荷和2D周期性的非立方单位细胞和板状系统.
主要成果:
- 对于中等规模的系统,DS方法的准确性和计算效率与Ewald方法相比相当.
- 对于散装阶段和板式系统,DS方法是有效的.
- 该方法在形式上是准确的,没有可调节的参数.
结论:
- DS方法为分子动力学中的库伦相互作用提供了强大而准确的方法.
- 它的真实空间性质和简单的实施使其成为Ewald方法的有价值的替代方案.
- 对于具有非立方单位细胞和板状结构的系统,DS方法特别有用.
更多相关视频
相关概念视频
Van der Waals Interactions
70.1K
Atoms and molecules interact with each other through intermolecular forces. These electrostatic forces arise from attractive or repulsive interactions between particles with permanent, partial, or temporary charges. The intermolecular forces between neutral atoms and molecules are ion–dipole, dipole–dipole, and dispersion forces, collectively known as van der Waals forces.
70.1K
Coulomb's Law
11.3K
Experiments with electric charges have shown that if two objects each have an electric charge, they exert an electric force on each other. The magnitude of the force is linearly proportional to the net charge on each object and inversely proportional to the square of the distance between them. The direction of the force vector is along the imaginary line joining the two objects and is dictated by the signs of the charges involved.
Newton's third law applies to the Coulomb force — the...
Newton's third law applies to the Coulomb force — the...
11.3K
Van der Waals Equation
6.2K
The ideal gas law is an approximation that works well at high temperatures and low pressures. The van der Waals equation of state (named after the Dutch physicist Johannes van der Waals, 1837−1923) improves it by considering two factors.
First, the attractive forces between molecules, which are stronger at higher densities and reduce the pressure, are considered by adding to the pressure a term equal to the square of the molar density multiplied by a positive coefficient a. Second, the volume...
First, the attractive forces between molecules, which are stronger at higher densities and reduce the pressure, are considered by adding to the pressure a term equal to the square of the molar density multiplied by a positive coefficient a. Second, the volume...
6.2K
Calculations of Electric Potential II
2.2K
An electric dipole is a system of two equal but opposite charges, separated by a fixed distance. This system is used to model many real-world systems, including atomic and molecular interactions. One of these systems is the water molecule, but only under certain circumstances. These circumstances are met inside a microwave oven, where electric fields with alternating directions make the water molecules change orientation. This vibration is equivalent to heat at the molecular level.
Consider a...
Consider a...
2.2K
Calculations of Electric Potential I
2.5K
Consider a ring of radius R with a uniform charge density λ. What will the electric potential be at point M, which is located on the axis of the ring at a distance x from the center of the ring?
The ring is divided into infinitesimal small arcs such that point M is equidistant from all the arcs. Here, the cylindrical coordinate system is used to calculate the electric potential at point M. A general element of the arc between angles θ and θ + dθ is of the length Rdθ and has a charge of...
The ring is divided into infinitesimal small arcs such that point M is equidistant from all the arcs. Here, the cylindrical coordinate system is used to calculate the electric potential at point M. A general element of the arc between angles θ and θ + dθ is of the length Rdθ and has a charge of...
2.5K
Calculating the Equilibrium Constant
37.4K
The equilibrium constant for a reaction is calculated from the equilibrium concentrations (or pressures) of its reactants and products. If these concentrations are known, the calculation simply involves their substitution into the Kc expression.
For example, gaseous nitrogen dioxide forms dinitrogen tetroxide according to this equation:
For example, gaseous nitrogen dioxide forms dinitrogen tetroxide according to this equation:
37.4K


