def METHOD_NAME(
    energies,
    volumes,
    force_constants,
    structure,
    t_min,
    t_step,
    t_max,
    mesh,
    eos,
    pressure=0,
):
    """
    Compute QHA thermal expansion coefficient using the phonopy interface.
    Args:
        energies (list):
        volumes (list):
        force_constants (list):
        structure (Structure):
        t_min (float): min temperature
        t_step (float): temperature step
        t_max (float): max temperature
        mesh (list/tuple): reciprocal space density
        eos (str): equation of state used for fitting the energies and the volumes.
            options supported by phonopy: vinet, murnaghan, birch_murnaghan
        pressure (float): in GPa, optional.
    Returns:
        (numpy.ndarray, numpy.ndarray): thermal expansion coefficient, Temperature
    """
    # quasi-harmonic approx
    phonopy_qha = get_phonopy_qha(
        energies,
        volumes,
        force_constants,
        structure,
        t_min,
        t_step,
        t_max,
        mesh,
        eos,
        pressure=pressure,
    )
    # thermal expansion coefficient and temperature
    max_t_index = phonopy_qha._qha._len
    alpha = phonopy_qha.get_thermal_expansion()[:max_t_index]
    T = phonopy_qha._qha._temperatures[:max_t_index]
    return alpha, T