def _calc_diffusivity(self):
        """Calculate the diffusivity."""
        self.state.diffusivity = np.array([[ion.diffusivity(self.state.pH)]
                                          for ion in self.state.ions])