[论文解读] Dynamics of charged molecular strands

David Ellis, François Gay–Balmaz+3 more|ArXiv.org|Jan 19, 2009

Porphyrin and Phthalocyanine Chemistry参考文献 23被引用 7

一句话总结

本论文通过将屏蔽库仑或伦纳德-琼斯相互作用纳入几何杆理论的对流表示，推导出带电分子链（如DNA）动力学的非局部欧拉-庞加莱方程。主要贡献是一个统一的变分框架——采用改进的哈密顿-庞特里亚金和欧拉-庞加莱公式——得到仿射欧拉-庞加莱方程，并揭示其与复杂流体动力学的联系，其中四元数形式使非局部力矩和力的直接计算成为可能。

ABSTRACT

Euler-Poincare equations are derived for the dynamical folding of charged molecular strands (such as DNA) modeled as flexible continuous filamentary distributions of interacting rigid charge conformations. The new feature is that the equations of motion for the dynamics of such molecular strands are nonlocal when the screened Coulomb interactions, or Lennard-Jones potentials between pairs of charges are included. These nonlocal dynamical equations are derived in the convective representation of continuum motion by using modified Euler-Poincare and Hamilton-Pontryagin variational formulations that illuminate the various approaches within the framework of symmetry reduction of Hamilton's principle for exact geometric rods. In the absence of nonlocal interactions, the equations recover the classical Kirchhoff theory of elastic rods in the spatial representation. The motion equations in the convective representation are shown to be affine Euler-Poincare equations relative to a certain cocycle. This property relates the geometry of the molecular strands to that of complex fluids. An elegant change of variables allows a direct passage from the infinite dimensional point of view to the covariant formulation in terms of Lagrange-Poincare equations. In another revealing perspective, the convective representation of the nonlocal equations of molecular strand motion is transformed into quaternionic form.

研究动机与目标

• 开发一个一致的连续体模型，用于柔性带电分子链，同时考虑局部弹性形变和非局部静电相互作用。
• 解决长期存在的挑战：将非局部相互作用产生的长程力矩贡献整合进经典基尔霍夫杆理论。
• 在几何和对称框架下，统一多种变分公式——欧拉-庞加莱、仿射和协变拉格朗日-庞加莱——用于非局部杆动力学。
• 将几何杆理论的适用性扩展至包含屏蔽库仑和伦纳德-琼斯等非局部势能。
• 通过哈密顿结构和李-泊松括号，实现对带电链动态折叠和静止状态的解析洞察。

提出的方法

• 通过引入基于空间分离的卷积积分来包含非局部势能，利用改进的哈密顿-庞特里亚金变分原理推导运动方程。
• 应用改进的欧拉-庞加莱方法推导方程的对流表示，表明其构成一个具有特定上循环结构的仿射欧拉-庞加莱系统。
• 引入坐标变换，将动力学分解为水平和垂直分量，通过主丛上的拉格朗日-庞加莱约化实现协变公式。
• 采用非局部相互作用核的四元数表示，用四元数乘积替代向量叉积，简化非局部力矩的计算。
• 利用仿射欧拉-庞加莱方程的半直积公式推导的李-泊松括号，建立系统的哈密顿结构。
• 通过将非局部势能定义为沿链上电荷构型相对位置和取向的泛函，将理论应用于带电链的情形。

实验结果

研究问题

• RQ1 如何在几何杆框架内一致地将长程静电相互作用纳入柔性分子链的动力学？
• RQ2 当同时包含非局部相互作用产生的力和力矩时，所得非局部运动方程的数学结构是什么？
• RQ3 在非局部几何杆的背景下，欧拉-庞加莱、仿射和协变拉格朗日-庞加莱公式之间有何关系？
• RQ4 非局部动力学能否以更易计算的形式重新表述，例如使用四元数？
• RQ5 在固定纤维极限下，非局部系统中哪些守恒律和对称性得以保持？

主要发现

• 带电分子链的非局部运动方程被推导为相对于特定上循环的仿射欧拉-庞加莱方程，揭示了其与复杂流体动力学之间深刻的几何联系。
• 在无非局部相互作用时，方程退化为经典基尔霍夫弹性杆理论，验证了该框架作为自然扩展的合理性。
• 变量变换使得可直接从无穷维观点过渡到协变拉格朗日-庞加莱公式，明确揭示了动力学的几何结构。
• 非局部相互作用核以四元数形式表达，力矩贡献由势能对相对四元数位置导数的虚部给出。
• 通过推广经典形式的李-泊松括号建立哈密顿结构，包含非局部项，同时保持环流和动量守恒。
• 该框架可推广至高维情形，方程在多维参数化下保持其几何形式，但需满足交叉导数的可积性条件。