Dear all,
I am wondering if I can modify the syntax applied used in the custom dump file. The purpose is for some vibrationnal density of states calculations : I need to extract the velocity per atom every N timestep, to compute the autocorrelation function and FT it, using an external C code with fftw.The syntax between the dump file and the readable type of file for the external code are a bit different (I need to write the elapsed time at each timestep, to erase the “ITEM : etc” at every timestep, and some time to convert some component of the velocity from int to float - eg 5 to 5.0).
For the example, this is my current dump_custom line :
variable f file tmp.times
variable s equal next(f)
dump veloc all custom 5 crystal_v.vel vx vy vz
dump_modify veloc sort id
dump_modify veloc every v_s
For the moment I am preparing the simulations, so in the first place I figured I could just use a small matlab routine. But I wondered if it was possible to do it directly from lammps? Without modifying the dump_custom c++ file for controlling directly the printings, because I work on a cluster.
I have another general question, concerning VDAOS simulations : as a first guess, I would use a fix NVT to obtain the veolcity field, and pay attention to the frequency at witch it is dumped. Could somebody share his experience?
regards,
Arthur France-Lanord