anotações.txt

Equações adimensionais para MD
m' = adimensional
m* = valor dimensional de referência
[] vetor

m' = m/m* (massa)
sigma' = sigma/sigma* (distância)
epsilon' = epsion/epsilon* (energia = força*[distância]) No caso do LJ, epsilon*[distância]/distância^2 = força
[xi'] = [xi]/sigma* (distância)
Dimensões:

comprimento    sigma (LJ)
energia        epsilon (LJ)
massa          m
Velocidade     sqrt(epsilon/m)
força          epsilon/sigma 
pressão        epsilon/sigma^3
temperatura    epsilon*k_b

Adimensionais

[rij] = [rij]/sigma*
T' = T*k_b/epsilon
E' = E/epsilon*
V' = V/epsilon* (energia potencial) 
t' = t/alpha com alpha = sqrt( (m^2*sigma^2/epsilon*) )

O que dá pra fazer aqui para admensionalizar. 
Escolher uma das particulas como referência, esta terá valores de
sigma, m, epsilon = 1 
A velocidade das partículas determinada por T ou E será o que fará as particulas
se aglutinarem mais ou menos. 



Ek =  T'* epsilon 
T' = Ek/epsilon (numero grande) = Ek'


Raios atômicos:
https://en.wikipedia.org/wiki/Atomic_radii_of_the_elements_(data_page)
para o carbono usar uma média entre sigma e raio = 0.335, 0.170 (van der Waals)
https://www.sciencedirect.com/topics/chemistry/rhombohedral-graphite

Usar o fator de empacotamento (em geral de 0.64 a 0.74 BCC, FCC, HCF)

Ideias: 

Variar densidade (número)
Variar Temperatura

m_C = 12
m_Ar = 40
fe = (.74*(.335+.127)) # ou 0.64*
rs = 5
rC = 0.07
sigma_Ar = 0.341
nC = (sigma_Ar*rs)**2*fe/(rC**2) # numero de atomos na região
m_Cparticle = nC*m_C/m_Ar # massa adimensional do carbono
print(nC)
print(m_Cparticle)