Centering coordinates for energy calculation

Gary gary... at mail.mcgill.ca
Wed Jul 20 20:59:24 UTC 2016


Ahh I see. Thank you! I assumed it would be in angstrom because those are 
the units for coordinates.

Thanks,

Gary 

On Wednesday, July 20, 2016 at 6:24:20 AM UTC-7, Matt W wrote:
>
> Hi,
>
> this is a bit of a daft one (not you, cp2k). The center_point command 
> takes values in atomic units (I know that you gave the cell in angstrom 
> etc)!! So your molecule will appear off centre. If you convert to Bohr, 
> then you should get what you want (Unfortunately you can't specify units 
> for center_point at the moment).
>
> Matt
>
> On Tuesday, July 19, 2016 at 6:10:34 PM UTC+1, Gary wrote:
>>
>> Hello,
>>
>> I am trying to calculate the electric field above a single non-periodic 
>> molecule system. I am using a cell much larger than the molecule to avoid 
>> interactions between the periodic image (the cell is [26.0, 21.0, 15.0]).
>>
>> As I am only interested in looking at the field above the molecule, I 
>> have performed the geometry optimisation and the electric field calculation 
>> with the keyword
>>
>> &CENTER_COORDINATES
>>   CENTER_POINT 13.0 10.5 0.0
>> &END
>>
>> which should place the molecule in the center bottom of the cell. Looking 
>> at the coordinates of the optimised configuration, this is not the case. 
>> The molecule is in the bottom left bottom of the cell instead, and the 
>> resulting EFIELD_CUBE is cut off the edge. How do I place the molecule 
>> where I want it? Should I use the CENTER_COORDINATES at all during 
>> optimisation? 
>>
>> Attached is the input file.
>>
>> Thanks,
>>
>> Gary
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20160720/2a594488/attachment.htm>


More information about the CP2K-user mailing list