I Direct Gibbs energy calculation from molecular dynamics

ab_kein
Messages
3
Reaction score
0
I've seen a lot of articles, where people calculate free Gibbs energy of system using several methods in some sense indirect methods, but I've never seen it is being calculated using its definition:
$$ \Phi=U+pV-TS $$
Imagine a big MD system at its equilibrium in a box with periodic boundary conditions. Now consider a smaller sub-box with walls transparent for particles (atoms), where are k particles right now. Suppose we know each particle's position ##r_i##, velocity ##v_i##, potential energy ##\pi_i##, force acting on it ##f_i## and the sub-system's temperature ##T##. And, mainly, we know per-atom entropy ##s_i##. Can then the sub-system's free Gibbs energy be calculated as
$$\Phi_k=U_k+(pV)_k+T_kS_k;\quad U_k=\frac{m}{2}\sum_{i=1}^{k}v_i^2 + \sum_{i=1}^{k}\pi_i;\quad (pV)_k=kT_k-\frac{1}{3}\sum_{i=1}^k r_i\cdot f_i;\quad S_k=\sum_{i=1}^k s_i$$
 
Physics news on Phys.org
Looks like the last equation (entropy, S_k) was cut off.
 
Hi. I have got question as in title. How can idea of instantaneous dipole moment for atoms like, for example hydrogen be consistent with idea of orbitals? At my level of knowledge London dispersion forces are derived taking into account Bohr model of atom. But we know today that this model is not correct. If it would be correct I understand that at each time electron is at some point at radius at some angle and there is dipole moment at this time from nucleus to electron at orbit. But how...
Back
Top