COM motion removal during NVT simulation with electric field

User discussions
,
1

Hi,

I am studying the motion of a positively charged protein in the presence of electric field at 300 K. I have observed a conflict between maintaining the physical drift of the protein under the influence of electric field and preventing the COM motion of the system.

I have performed simulations with and without the COM removal: (a) comm-mode = Linear and (b) comm-mode = None.

The following are the observations:

  • In case (a), the displacement of the protein is reduced, ~20 nm, implying that the COM removal is fighting the field-induced motion.

  • In case (b), the protein covers ~ 300 nm in the same time interval. But the COM of the system experiences a drift, (i.e.) COM velocity of the system along the direction of the electric field is ~ 0.009 nm/ps.

  • However the temperature of the protein and the system are at 300 K in both the cases.

Question:
What is the recommended approach to perform NVT simulation, of such a system, in the presence of electric field while preventing the artifacts due to the COM motion?

I would greatly appreciate any insights or references to similar studies.

Thank you.

Best regards,
Afi

2

I don’t know what scientific question you are trying to answer. If you don’t, you need to figure that out.

Putting an infinite periodic array of charged proteins with very short separation in between in an electric field is a very artificial situation.

3

Dear Prof. @hess,

Thank you for your swift response.

To clarify, the ultimate objective is to study protein translocation through a solid-state nano channel under the influence of electric field.

The system mentioned, a protein in water and ions, is a simplified model to understand the effect of COM motion removal on the movement of protein.

Best regards,
Afi