QuantumATK Forum

QuantumATK => General Questions and Answers => Topic started by: Shan on December 4, 2017, 15:20

Title: pH value - molecular dynamics
Post by: Shan on December 4, 2017, 15:20
Dear Experts,

I am eager to know if there is any way to adjust the pH value of solution inside the supercell, while performing molecular dynamics.

Thanks.
Title: Re: pH value - molecular dynamics
Post by: Ulrik G. Vej-Hansen on December 5, 2017, 16:27
That depends on how you plan to model your solution, but as far as i am aware, it is quite tricky. It will be easier to guide you how to do it in VNL/ATK, if you can provide a reference describing the method you are using.
Title: Re: pH value - molecular dynamics
Post by: Shan on December 5, 2017, 19:30
Dear Dr. Ulrik,

I want to check the dispersion behaviour of graphene flakes in an aqueous solution at different pH levels.
I want to use "Molecular Dynamics" incorporated in ATK.
All I could not get is, How to vary the pH levels of aqueous solution in ATK-VNL ???

Awaiting your reply,
thank you so much.
Title: Re: pH value - molecular dynamics
Post by: Julian Schneider on December 15, 2017, 09:24
Modelling a given pH value explicitly, i.e. via introducing a certain number of H3O+ ions is normally difficult because, unless you go to very extreme pH values, the H3O+ concentration is so low that you'd have to use a huge number of water molecules. That also requires that you'd need to use DFT or at least a very good reactive force field, to correctly model proton transfer.

Another more implicit approach is, if you know what effect the pH has on your system (e.g. if you have acidic groups) , to manually fix the protonation state of these groups according to the given pH, e.g. if you have an acidic group and you want to simulate neutral pH, the acidic group would most likely be present in the deprotonated state.

A good thing to start with  could be to try and find out what the effect of pH is on the components of your system, for example if water molecules might dissociate at the edges of your graphene flake, leading to adsorbed OH groups, whose behavior might then be pH dependent.