Re: [AMBER-Developers] XMIN and writing restart and trajectory files

From: Ben Roberts <roberts.qtp.ufl.edu>
Date: Mon, 10 May 2010 09:40:38 -0400

Hi Istvan (and Andy),

On 8/5/2010, at 4:07 a.m., istvan.kolossvary.hu wrote:

> Hi, Ben, Andy,
>
> Tell me exactly what you want. You would write out coords from the xmin wrapper, but the exact place may depend on what your goal is.
>
>> Hi Ben,
>>
>> Ben Roberts wrote:
>>> Hi all,
>>>
>>> I'm wondering if someone could give me a few pointers here.
>>>
>>> I'm trying to include writing restart and trajectory files in the XMIN algorithm, initially for my use in QM/MM calculations, so I can save files at each step instead of having to go back to the start if the calculation fails (since each step may take an hour or more).
>>>
>>> My initial solution, which works passably well, was to call minrit and corpac (based on ntwr and ntwx) at the end of the main do loop in run_xmin. (At the moment, this is only done if run_xmin is called directly, not from run_lmod.) This is based on n_force_calls. The problem is that some steps do more than one force call, so that a trajectory written on every force call contains more steps than the corresponding mdout file. Ideally, the steps counted to determine ixdump and itdump would be those that are counted for mdout-writing purposes.
>>
>> Each xmin step does a line search which requires several gradient
>> evaluations (force calls). You don't want to write coordinates during
>> the line search but only for each new xmin step.
>>
>>> Inside the do loop, I can see a rather complicated case switch, and it's not obvious to me which branches start a new "step", and which ones do another force evaluation within the same "step". Could anyone clarify which is which? I would guess that "case (CALCENRG_NEWNBL, CALCGRAD_NEWNBL, CALCBOTH_NEWNBL)" does not do a new step and that "case (CALCENRG_OLDNBL, CALCGRAD_OLDNBL, CALCBOTH_OLDNBL)" does, but I'm not certain...
>>
>> As far as I can see, a new xmin step is done for
>>
>> case ( CALCENRG_NEWNBL, CALCGRAD_NEWNBL, CALCBOTH_NEWNBL )
>>
>> which is also when the minimization progress is printed. This is where
>> you want to write the coordinates.
>>
>> All the best,
>> Andy

I'm trying to use XMIN to do QM/MM minimisations with PUPIL. If I'm doing, for example, a large QM region (70+ atoms) at the DFT level, a single force evaluation may take an hour or more. Hence, I would like to be able to save a restart file often, up to every step.

Furthermore, I would like to be able to write trajectories (mdcrd or netcdf), again up to once per mdout step, so I can see what's happening in the system over "time".

At the moment, as outlined, it's doing every force evaluation. But the problem there is that a single step (as counted in the mdout file) may have multiple force evaluations, and it's not obvious when one step begins and ends. Nor is it obvious that a restart file generated in such an "intermediate" force evaluation is safe to pick up and run with if things bomb out.

So I would like to have a way to write restrt and mdcrd frames based on steps as counted in the mdout file.

I hope that helps clarify things.

Cheers,
Ben
_______________________________________________
AMBER-Developers mailing list
AMBER-Developers.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber-developers
Received on Mon May 10 2010 - 07:00:04 PDT
Custom Search