code: def run_thermal_displacement_matrices(self,
                                            t_min=0,
                                            t_max=1000,
                                            t_step=10,
                                            temperatures=None,
                                            freq_min=None,
                                            freq_max=None):
          
          if self._dynamical_matrix is None:
              msg = ("Dynamical matrix has not yet built.")
              raise RuntimeError(msg)
          if self._mesh is None:
              msg = ("run_mesh has to be done.")
              raise RuntimeError(msg)
          mesh_nums = self._mesh.mesh_numbers
          ir_grid_points = self._mesh.ir_grid_points
          if not self._mesh.with_eigenvectors:
              msg = ("run_mesh has to be done with with_eigenvectors=True.")
              raise RuntimeError(msg)
          if np.prod(mesh_nums) != len(ir_grid_points):
              msg = ("run_mesh has to be done with is_mesh_symmetry=False.")
              raise RuntimeError(msg)
  
          tdm = ThermalDisplacementMatrices(
              self._mesh,
              freq_min=freq_min,
              freq_max=freq_max,
              lattice=self._primitive.get_cell().T)
  
          if temperatures is None:
              tdm.set_temperature_range(t_min, t_max, t_step)
          else:
              tdm.set_temperatures(temperatures)
          tdm.run()
  
          self._thermal_displacement_matrices = tdm
docs: Prepare thermal displacement matrices
  
          Parameters
          ----------
          t_min, t_max, t_step : float, optional
              Minimum and maximum temperatures and the interval in this
              temperature range. Default valuues are 0, 1000, and 10.
          freq_min, freq_max : float, optional
              Phonon frequencies larger than freq_min and smaller than
              freq_max are included. Default is None, i.e., all phonons.
          temperatures : array_like, optional
              Temperature points where thermal properties are calculated.
              When this is set, t_min, t_max, and t_step are ignored.
              Default is None.