code: def get_conversion_factor(structure, species, temperature):
      
      df_sp = get_el_sp(species)
      if hasattr(df_sp, "oxi_state"):
          z = df_sp.oxi_state
      else:
          z = df_sp.full_electronic_structure[-1][2]
  
      n = structure.composition[species]
  
      vol = structure.volume * 1e-24  # units cm^3
      return 1000 * n / (vol * const.N_A) * z ** 2 * (const.N_A * const.e) ** 2 \
          / (const.R * temperature)
docs: Conversion factor to convert between cm^2/s diffusivity measurements and
      mS/cm conductivity measurements based on number of atoms of diffusing
      species. Note that the charge is based on the oxidation state of the
      species (where available), or else the number of valence electrons
      (usually a good guess, esp for main group ions).
  
      Args:
          structure (Structure): Input structure.
          species (Element/Specie): Diffusing species.
          temperature (float): Temperature of the diffusion run in Kelvin.
  
      Returns:
          Conversion factor.
          Conductivity (in mS/cm) = Conversion Factor * Diffusivity (in cm^2/s)