Dear Matthew:
good afternoon.
Sorry to bother once more.
I took Gly and Ala oriented the way that their COO-Ca-N plane is parallel to XZ plane
and modeled charge dynamics propagating under such 5 field settings
{Ex = 0, Ey = 0} no field
{Ex=0, Ey} polarized along Y
{Ex, Ey = 0} polarized along X
{Ex, Ey + P/2 shifted} left circular polarized
{Ex, Ey - P/2 shifted} right circular polarized.
The initial geometries are proximal to the global minimum optimized under DFT.
I attach an image where I sum-up the results.
When there are no fields, apparently, there is just the equillibration: energy flows into kinetic populating vibrations,
and there is the intra-molecular H dynamics between the terminals.
Field polarized along X (parallel to the
COO-Ca-N plane
) makes thermalization slowe.
Field polarized along Y (perpendicular to
the COO-Ca-N plane) introduce strong anticorrelation of side groups charges for Gly,
but it does not slow down the thermalization for Gly.
Instead, in the case of Ala, the temperature and the kinetic energy explode at a certain point - the shake up start to break the molecule.
At the same time, the shake up for the potential energy is about the same as when field is parallel to X.
This is peculiar since one could think that shaking between the terminals should have larger effect.
Using CD polarized light, phase of charge modulations of the side groups, only, reflect "chirality" of the drive.
In the case of Ala, cp2k simulates both,
the temperature and the kinetic energy explode at a certain point.
I simulate using energy constraining NVE and instructing T=300K.
Potential energy does change smooth, always.
I wish to ask your opinion: the increase of temperature as cp2k simulate is it a signature
of investing energy into molecule that it start to break it apart "properly",
or it is a "proper" mis-performance of the NVE part of the equations,
or it is a signature of the equations collapse?
I doubt about the third because I tested increasing the amplitude of the drive: at certain point
the potential energy start to show a ripple.
The data I attach here does show potential energy to vary smoothly.
Therefore, I think it is likely that the molecule is getting very hot.
Thank you.
With best,
Victor