3D分子动画革命：用Manim构建化学键与反应可视化

2026-02-04 04:05:52作者：俞予舒Fleming

你是否曾为化学课件中的静态分子结构图感到枯燥？是否想让学生直观看到化学键的形成与断裂过程？Manim（数学动画引擎）提供了强大的3D建模能力，能将抽象的分子结构转化为生动的动态演示。本文将带你从零开始创建第一个分子动画，掌握原子建模、化学键连接和反应过程动画的核心技术。

分子建模基础：从柏拉图立体到原子模型

Manim的三维建模模块为分子可视化提供了坚实基础。通过manim/mobject/three_d/polyhedra.py中定义的多面体类，我们可以轻松创建各种原子模型。以甲烷分子(CH₄)为例，碳原子通常用正四面体表示，这可以通过Tetrahedron类实现：

from manim import *
from manim.mobject.three_d.polyhedra import Tetrahedron

class MethaneMolecule(ThreeDScene):
    def construct(self):
        self.set_camera_orientation(phi=75*DEGREES, theta=30*DEGREES)
        
        # 创建碳原子（红色正四面体）
        carbon = Tetrahedron(edge_length=1.5)
        carbon.set_color(RED)
        carbon.set_fill(opacity=0.8)
        
        # 创建氢原子（白色球体）
        hydrogens = [
            Sphere(radius=0.3).set_color(WHITE).shift(UP*1.5),
            Sphere(radius=0.3).set_color(WHITE).shift(DOWN*1.5),
            Sphere(radius=0.3).set_color(WHITE).shift(LEFT*1.5),
            Sphere(radius=0.3).set_color(WHITE).shift(RIGHT*1.5)
        ]
        
        # 创建化学键（灰色圆柱体）
        bonds = [
            Cylinder(radius=0.1, height=1.5, color=GRAY).shift(UP*0.75),
            Cylinder(radius=0.1, height=1.5, color=GRAY).shift(DOWN*0.75),
            Cylinder(radius=0.1, height=1.5, color=GRAY).shift(LEFT*0.75),
            Cylinder(radius=0.1, height=1.5, color=GRAY).shift(RIGHT*0.75)
        ]
        
        self.add(carbon, *hydrogens, *bonds)
        self.begin_3dillusion_camera_rotation()  # 启动3D旋转
        self.wait(5)

这段代码利用了Manim的ThreeDScene场景和Polyhedron类，创建了具有空间立体感的分子模型。其中begin_3dillusion_camera_rotation()方法来自manim/scene/three_d_scene.py，实现了自动旋转视角，让观众能从不同角度观察分子结构。

进阶技巧：动态化学键与反应动画

静态分子模型只是起点，Manim真正的强大之处在于创建动态反应过程。以下是水分解反应（2H₂O → 2H₂ + O₂）的动画实现，展示了如何模拟化学键的断裂与形成：

class WaterDecomposition(ThreeDScene):
    def construct(self):
        self.set_camera_orientation(phi=60*DEGREES, theta=45*DEGREES)
        
        # 创建两个水分子
        def create_water_molecule(position):
            # 氧原子（红色球体）
            oxygen = Sphere(radius=0.4, color=RED).shift(position)
            # 氢原子（白色球体）
            hydrogen1 = Sphere(radius=0.25, color=WHITE).shift(position + LEFT*0.8 + UP*0.5)
            hydrogen2 = Sphere(radius=0.25, color=WHITE).shift(position + RIGHT*0.8 + UP*0.5)
            # 化学键（灰色圆柱体）
            bond1 = Cylinder(radius=0.08, height=1.0, color=GRAY)
            bond1.rotate(PI/4, axis=OUT)
            bond1.shift(position + LEFT*0.4 + UP*0.25)
            bond2 = Cylinder(radius=0.08, height=1.0, color=GRAY)
            bond2.rotate(-PI/4, axis=OUT)
            bond2.shift(position + RIGHT*0.4 + UP*0.25)
            return VGroup(oxygen, hydrogen1, hydrogen2, bond1, bond2)
        
        # 创建两个水分子并放置在不同位置
        water1 = create_water_molecule(LEFT*3)
        water2 = create_water_molecule(RIGHT*3)
        self.add(water1, water2)
        self.wait(1)
        
        # 水分子分解动画
        self.play(
            # 化学键断裂
            water1[3].animate.scale(0),  # 第一个键消失
            water1[4].animate.scale(0),  # 第二个键消失
            water2[3].animate.scale(0),
            water2[4].animate.scale(0),
            # 氢原子分离
            water1[1].animate.shift(LEFT*1.5 + DOWN*1),
            water1[2].animate.shift(RIGHT*1.5 + DOWN*1),
            water2[1].animate.shift(LEFT*1.5 + UP*1),
            water2[2].animate.shift(RIGHT*1.5 + UP*1),
            run_time=3
        )
        
        # 形成氢气分子
        def create_hydrogen_molecule(position):
            h1 = Sphere(radius=0.25, color=WHITE).shift(position + LEFT*0.5)
            h2 = Sphere(radius=0.25, color=WHITE).shift(position + RIGHT*0.5)
            bond = Cylinder(radius=0.08, height=1.0, color=GRAY).shift(position)
            return VGroup(h1, h2, bond)
        
        h2_1 = create_hydrogen_molecule(DOWN*2)
        h2_2 = create_hydrogen_molecule(UP*2)
        o2 = VGroup(
            water1[0], water2[0],  # 两个氧原子
            Cylinder(radius=0.08, height=2.0, color=GRAY).shift(ORIGIN)  # 氧氧键
        )
        
        self.play(
            Create(h2_1[2]), Create(h2_2[2]), Create(o2[2]),
            run_time=2
        )
        self.add(h2_1, h2_2, o2)
        self.wait(2)