そのメカニズムの解明に有効な手法の1つに分子動力学法（MD法）がある．MD法によってVACNTの圧縮座屈シミュレーションを行い，不安定性解析を行うことでVACNTの圧縮座屈挙動を明らかにできる．不安定性解析にはMD法で用いるポテンシャル関数の2次微分項からなるヘッシアン行列の固有値解析をすれば良い．しかし，ポテンシャル関数が非常に複雑であるため，この2次微分項を求めることは非常に困難である．そこで，本研究では1次微分項を数値的に微分することで対応する．数値微分には精度に問題がある可能性があるため，まずはその有効性を検討する．そして，圧縮変形を受けるVACNTの不安定性解析を行い，座屈点および 座屈モードを明らかにする．