# diffusion coefficient in cm^2/s from the slope of the MSD in Bohr^2
# and time step size in au_t
diffusion <- function(slope,dt) { return(0.192946*slope/dt) }