%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
metalls = pd.read_csv('data/metalls.csv')
metalls
Primer carregarem aquestes dades, tenint en compte que el radi atòmic el tenim en nm i volem la densitat en g/cm³
També definiren el nomnbre d'Avogadre Na
r=metalls['radi']*10**(-7)
A=metalls['A']
rho=metalls['rho']
Na=6.023*10**23
En aquest cas el costat de la cel·la unitat és $a=\frac{4}{\sqrt{2}}\cdot r$
En cada cel·la tenim 4 atoms complets
La densitat es pot calcular com
$\rho = \frac{4}{a^{3}}\cdot \frac{A}{N_{a}}$
n=13
a=4/np.sqrt(2)*r[n]
densitat=4/a**3*A[n]/Na
densitat
Com veiem es tracta d'una bona estimació del valor real
rho[n]
En aquest cas el costat de la cel·la unitat és $a=\frac{4}{\sqrt{3}}\cdot r$
En cada cel·la tenim 2 atoms complets
La densitat es pot calcular com
$\rho = \frac{2}{a^{3}}\cdot \frac{A}{N_{a}}$
n=5
a=4/np.sqrt(3)*r[n]
densitat=2/a**3*A[n]/Na
densitat
Que s'acosta també molt al valor real
rho[n]