<div dir="ltr">Dear molpro users/developers,<div><br></div><div style>I'm trying to perform a AIMS trajectory specifying initial momenta instead of doing standard samplings. The problem is that the initial momenta given is not being read. The corresponding files are as follow:</div>

<div style><br></div><div style><b>Geometry.dat:</b></div><div style><br></div><div style><div>   UNITS=ANG</div><div>        5</div><div>X1  0.226707876  1.54860805  7.12387356</div><div>X2 -0.201399539  1.10690738  7.45922101</div>

<div>X3 -0.135985033  1.89880715  9.26828747</div><div>X4 -0.064964325  1.78959108  8.36970196</div><div>X5  0.175641021  2.65608635  8.80329899</div><div>Momenta:</div><div>X1  5.2742467920  0.4623106270  0.1090012330</div>

<div>X2 -0.3497877210 -0.2012876050  0.4642698440</div><div style>X3  0.0815010823  0.3953913210 -0.4911153070</div><div style>X4  0.0947156466 -0.1953009890 -0.0691696517</div><div style>X5 -0.1006758000 -0.4611133530 -0.0129861184</div>

</div><div style><br></div><div style>--------------------------------------------------------------------------</div><div style><br></div><div style><b>Control.dat:</b></div><div style><br></div><div style><div> &control</div>

<div><br></div><div>!------ Initialization -------!</div><div> InitState=1</div><div> iRndSeed=101</div><div> InitialCond='NoSample'</div><div> Temperature=10.0</div><div><br></div><div>!------ Simulation Properties -------!</div>

<div> IMethod=0              !Single- state dynamics - no couplings.</div><div> NumParticles=5         !Number of atoms</div><div> NumStates=1            !number of electronic states</div><div> IRestart=0</div><div><br></div>

<div>!-------- Integration and numerics --------!</div><div> TimeStep=2.06706          !timestep for integration</div><div> SimulationTime=206.706  !integration time</div><div>!MinTimeStep=1.0        !minimum adaptive timestep</div>

<div>!------- Ouput options ---------!</div><div style>....</div><div style><br></div><div style>-------------------------------------------------------------------</div><div style><br></div><div style><b>Input:</b></div>

<div style><br></div><div style><div>memory,400,M</div><div>SYMMETRY,nosym</div><div>ORIENT,noorient</div><div>GEOMTYP=xyz</div><div>GEOMETRY={</div><div>   5</div><div>QCISD(T)/xxx  ENERGY=?</div><div>X1  0.226707876  1.54860805  7.12387356</div>

<div>X2 -0.201399539  1.10690738  7.45922101</div><div>X3 -0.135985033  1.89880715  9.26828747</div><div style>X4 -0.064964325  1.78959108  8.36970196</div><div style>X5  0.175641021  2.65608635  8.80329899</div><div>}</div>

<div>BASIS=xxx</div><div><br></div><div>MPINIT={</div><div>{hf;wf,X,0,0}</div><div>qcisd(t)</div><div>}</div><div><br></div><div>MPCALC={</div><div> ! Calculate Energy and Forces</div><div> show,itimestep,itrj;</div><div>

 {hf;wf,X,0,0}</div><div> qcisd(t)</div><div> ener(1)=energy(1)</div><div> pop;density</div><div> force</div><div>}</div><div><br></div><div>aims,dir='.'</div><div><br></div><div>-------------------------------------------</div>

<div><br></div><div style>In the output E.dat, the EKineticCl is always 0.0 for the initial step, and the all the momenta outputed TrajDump.1 are 0.0 for time 0, regardless my input in Geometry.dat.</div><div style><br></div>

<div style>Any help will be very appreciated.</div><div style><br></div><div style>Best wishes,</div><div style><br></div><div style>Gabriel</div></div></div></div>