Boundary thermostat

[Ashley Dickson]

Hello all,

I want to run some collision cascades with a boundary thermostat (so as to not affect the atom velocities in the cascade). I have looked through the cp2k documentation and can't find anything regarding a boundary thermostat. Am I missing something or is this still not implemented?

Many thanks 

Ash

[Krack Matthias]

Hi Ash

 

There is an example input for cascade simulations in UO2 using the force field of Morelon with ZBL. The system size is with 2x2x2 unit cells of course too small and has to be increased significantly (the RCUT range in the &FORCEFIELD section should then be increased to 10-12 A too). You have to uncomment the &THERMAL_REGION section and add a LIST of all atomic indices referring to the boundary thermostat region. I have been using an onion-like order for the UO2 unit cell sequence which allows for an easy definition of these atomic indices by taking advantage of the range feature of the LIST keyword, e.g. LIST 1..5000.

 

HTH

 

Matthias

    

--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+uns...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/0605699a-96ba-4142-a39f-56db50d4931fn%40googlegroups.com.

[Ashley Dickson]

This is very helpful, thank you Matthias!