一个计算效率高的子空间波放松算法,用于粗粒度的分子系统,几乎完全的热力学一致性
João V M Pimentel1, Vladimir A Mandelshtam1
1Department of Chemistry, University of California, Irvine, California 92697, USA.
The Journal of chemical physics
|October 1, 2025
概括
这项研究增强了对刚性集群的粗粒度 (CG) 分子建模方法. 算法改进使得CG潜能在计算上与全原子 (AA) 模型可比,使得像蛋白质这样的大型系统的高效模拟成为可能.
科学领域:
- 计算化学计算化学
- 分子动力学分子动力学
- 生物物理学的生物物理.
背景情况:
- 一项先前的研究引入了一种直接粗粒度 (CG) 方法,用于从全原子 (AA) 潜力的刚性分子.
- 最初的CG方法虽然对水和氨集群准确,但由于子空间最小化和正常模式计算,由于计算成本昂贵.
研究的目的:
- 扩大CG硬化方法的适用性,包括对像蛋白质这样的大分子.
- 引入算法改进,以降低CG潜力评估的计算成本.
- 为了实现与分子模拟的AA模型可比的计算成本.
主要方法:
- 开发了对刚性分子集群的直接CG潜力构造的增强配方.
- 实现了算法优化,以减少子空间最小化和正常模式分析的计算开销.
- 在水集群上进行模拟,以证明改进的CG模型的效率和准确性.
主要成果:
- 改进的CG方法显著降低了子空间最小化和正常模式计算的成本.
- 增强的CG模拟需要较少的蒙特卡洛步骤来实现与AA模型相比的可比统计准确性.
- 改进的CG方法的整体计算成本与传统的AA方法具有竞争力.
结论:
- 精细的CG刚性化技术为模拟分子和大型生物分子提供了更具计算可行的途径.
- 算法进步使直接的CG潜力构建成为复杂系统的AA模拟的可行和高效替代方案.
- 这项工作为将高级CG建模应用于更大,更复杂的分子系统铺平了道路.
更多相关视频
相关概念视频
Equilibrium Conditions for a Particle
2.2K
When an object is in equilibrium, it is either at rest or moving with a constant velocity. There are two types of equilibrium: static and dynamic. Static equilibrium occurs when an object is at rest, while dynamic equilibrium occurs when an object is moving with a constant velocity. In both cases, there must be a balance of forces acting on the object.
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
2.2K
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
290
Mechanistic models play a crucial role in algorithms for numerical problem-solving, particularly in nonlinear mixed effects modeling (NMEM). These models aim to minimize specific objective functions by evaluating various parameter estimates, leading to the development of systematic algorithms. In some cases, linearization techniques approximate the model using linear equations.
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
290
Thermal Sigmatropic Reactions: Overview
2.4K
Sigmatropic rearrangements are a class of pericyclic reactions in which a σ bond migrates from one part of a π system to another. These are intramolecular rearrangements where the total number of σ and π bonds remain unchanged.
Sigmatropic shifts are classified based on an order term [i, j ], where i and j indicate the number of atoms across which each end of the σ bond migrates. Below are examples of a [3,3] sigmatropic shift in 1,5-hexadiene, referred...
Sigmatropic shifts are classified based on an order term [i, j ], where i and j indicate the number of atoms across which each end of the σ bond migrates. Below are examples of a [3,3] sigmatropic shift in 1,5-hexadiene, referred...
2.4K
IR Spectroscopy: Hooke's Law Approximation of Molecular Vibration
2.8K
A covalently bonded heteronuclear diatomic molecule can be modeled as two vibrating masses connected by a spring. The vibrational frequency of the bond can be expressed using an equation derived from Hooke's law, which describes how the force applied to stretch or compress a spring is proportional to the displacement of the spring. In this case, the atoms behave like masses, and the bond acts like a spring.
According to Hooke's law, the vibrational frequency is directly proportional to...
According to Hooke's law, the vibrational frequency is directly proportional to...
2.8K
Hybridization of Atomic Orbitals II
48.0K
sp3d and sp3d 2 Hybridization
48.0K
¹H NMR of Conformationally Flexible Molecules: Temporal Resolution
1.3K
At room temperature, the chair conformer of cyclohexane undergoes rapid ring flipping between two equivalent chair conformers at a rate of approximately 105 times per second. These two chair conformers are in equilibrium. The rapid ring flipping results in the interconversion of the axial proton to an equatorial proton and an equatorial to the axial proton. Such interconversions are too rapid and cannot be detected on the NMR timescale. Hence, the NMR spectrometer cannot distinguish between the...
1.3K


