def METHOD_NAME(self) -> List[str]:
    """Formatted result as a list of strings"""
    lines = []
    lines.append("=== GROUND STATE ===")
    lines.append(" ")
    lines.append(
        "* Vibrational ground state energy "
        f"(cm^-1): {np.round(self.computed_vibrational_energies[0], self.formatting_precision)}"
    )
    if len(self.num_occupied_modals_per_mode) > 0:
        lines.append("The number of occupied modals for each mode is: ")
        for i, m in enumerate(self.num_occupied_modals_per_mode[0]):
            lines.append(f"- Mode {i}: {np.round(m, self.formatting_precision)}")
    if (
        self.computed_vibrational_energies is not None
        and len(self.computed_vibrational_energies) > 1
    ):
        lines.append(" ")
        lines.append("=== EXCITED STATES ===")
        lines.append(" ")
        for idx, vib_energy in enumerate(self.computed_vibrational_energies[1:]):
            lines.append(
                f"* {(idx + 1): 3d}: Vibrational excited state energy "
                f"(cm^-1): {np.round(vib_energy, self.formatting_precision)}"
            )
            if idx < len(self.num_occupied_modals_per_mode):
                lines.append("The number of occupied modals for each mode is")
                for i, m in enumerate(self.num_occupied_modals_per_mode[idx]):
                    lines.append(f"- Mode {i}: {np.round(m, self.formatting_precision)}")
            lines.append(" ")
    return lines