凤凰时尚网移动版

$atom->Z);
        }   
    }
    #close trajectory file
    $xmolFile->Close;
}
else {
    print "The " . $doc->Name . " is not a multiframe trajectory file \n";
}

使用时先将这些内容复制到一个文本文件里，后缀名改为.pl。然后在MS里将这个.pl加入到项目中。要转换哪个目录下的xtd文件就把这pl文件挪到哪个目录中，并且把my $doc = $Documents{"./benzene.xtd"}; 当中的文件名改成要转换的文件名。之后，保持此脚本文件窗口处于激活状态，选tools-scripting-debug（或者直接按F5，或者按工具栏的蓝色三角按钮）就开始对xtd文件进行转换，转换结束后在当前目录下会输出trj.txt文件。将其后缀改为.xyz之后就能被VMD等程序直接读取了。

debug模式对于大体系、帧数较多的轨迹转换起来颇慢，可以用tools-scripting-Run on server模式来运行，这样转换速度明显快得多，trj.txt将会生成到新的目录，当前目录下的其它文件也会被强行复制过去一份。

此脚本默认转换所有帧。如果想转换指定帧数范围，就把my $framebegin=1;和my $frameend=10;改成自定的起止帧号就行了，需要先将my $frameend=10;前面的注释去掉。

如果是周期性体系，那么在MS当中看起来轨迹是什么样转换过去就是什么样。比如，如果在display style-lattice中在某个方向上多显示一个周期，那么转换出的轨迹在相应方向上也会多出一倍原子。Default、In-Cell、Original的显示模式下转换出的原子坐标也会相应地可能有所不同。

上面的这个名为xtd2xyz脚本只能转换全原子模拟的轨迹，如Forcite的xtd轨迹，却不能转换粗粒化模拟的轨迹，如Mesocite的以bead描述粒子的xtd轨迹。我将之修改成下面的xtdbead2xyz脚本，专门用来转换粗粒化模拟的轨迹（但不能转换全原子的），用法同前。

#!perl
# XTDbead2XYZ - Convert the XTD files containing beads into XYZ format
# Creator: Sobereva (sobereva@sina.com)
# Date:    2012-May-23

use strict;
use MaterialsScript qw(:all);

#open the multiframe trajectory structure file or die
my $doc = $Documents{"./bilayer.xtd"};

if (!$doc) {die "no document";}

my $trajectory = $doc->Trajectory;

if ($trajectory->NumFrames>1) {

    print "Found ".$trajectory->NumFrames." frames in the trajectory\n";
    # Open new xmol trajectory file
    my $xmolFile=Documents->New("trj.txt");
   
    #get atoms in the structure
    my $Beads = $doc->DisplayRange->Beads;
    my $NBeads=@$Beads;

    # loops over the frames
    my $framebegin=1;
    my $frameend=$trajectory->NumFrames;
#    my $frameend=10;
    for (my $frame=$framebegin; $frame<=$frameend; ++$frame){
        $trajectory->CurrentFrame = $frame;
        #write header xyz
        $xmolFile->Append(sprintf "%i \n", $NBeads);
        $xmolFile->Append(sprintf "%s %i \n", "Frame",$frame);
        foreach my $Bead (@$Beads) {
            # write atom symbol and x-y-z- coordinates
            $xmolFile->Append(sprintf "%s %f  %f  %f \n",$Bead->Name, $Bead->X, $Bead->Y, $Bead-