def METHOD_NAME(self, temps=(100, 200, 300, 400), **kwargs) -> Figure:
    """
    Plot phonon band structure with markers proportional to the occupation
    of each phonon mode for different temperatures.
    Args:
        temps: List of temperatures in Kelvin.
    Return: |matplotlib-Figure|
    """
    temps = np.array(temps)
    ntemp = len(temps)
    # Build plot grid.
    num_plots, ncols, nrows = ntemp, 1, 1
    if num_plots > 1:
        ncols = 2
        nrows = (num_plots // ncols) + (num_plots % ncols)
    ax_list, fig, plt = get_axarray_fig_plt(None, nrows=nrows, ncols=ncols,
                                            sharex=True, sharey=True, squeeze=False)
    ax_list = ax_list.ravel()
    for ax, temp in zip(ax_list, temps.ravel()):
        self.phb_qpath.plot(ax=ax, units="eV", temp=temp, fontsize=8, show=False)
    return fig