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.
 
From the BCS theory of superconductivity is well known that the superfluid density smoothly decreases with increasing temperature. Annihilated superfluid carriers become normal and lose their momenta on lattice atoms. So if we induce a persistent supercurrent in a ring below Tc and after that slowly increase the temperature, we must observe a decrease in the actual supercurrent, because the density of electron pairs and total supercurrent momentum decrease. However, this supercurrent...
Back
Top