chenzhe 发表于 2023-6-19 16:58
用脚本处理。help里面有关于粒子初速度的设置。
比如下面这两行的意思就是将选定原子初速度设置为-10
fo ...
chenzhe 发表于 2023-6-19 16:58
用脚本处理。help里面有关于粒子初速度的设置。
比如下面这两行的意思就是将选定原子初速度设置为-10
fo ...您好,我刚接触到perl,看了一下HELP,改了之后发现不太对
my $doc = $Documents{"daosan.xsd"};
my $H2 = $doc-> Unitcell-> Sets("H2")->Atoms;
foreach my $atom (@$H2){$atom->Velocity = Point ( X => -10, Y => 0, Z => 0 );}
my $results = Modules->Forcite->Dynamics->Run($doc, Settings(
"3DPeriodicvdWAtomCubicSplineWidth" => 3,
"3DPeriodicNonBond2BodyListBufferWidth" => 1,
CurrentForcefield => "pcff",
Ensemble3D => "NVT",
Pressure => 0.0001,
NumberOfSteps => 1000000,
TrajectoryFrequency => 1000,
InitialVelocities => "Random",
StressXX => -0.0001,
StressYY => -0.0001,
StressZZ => -0.0001));
my $outTrajectory = $results->Trajectory;
我的晶胞是无定形晶胞,不能使用Unitcell,应该怎么改呢,而且这样设置速度的话,是初始速度还是整个过程的速度,谢谢。
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) | Powered by Discuz! X3.3 |