def METHOD_NAME(input_file, library_flag):
    """
    Estimate thermo for a list of species using RMG and the settings chosen inside a thermo input file.
    """
    rmg = RMG()
    rmg.load_thermo_input(input_file)
    rmg.database = RMGDatabase()
    path = os.path.join(settings['database.directory'])
    # forbidden structure loading
    rmg.database.load_thermo(os.path.join(path, 'thermo'), rmg.thermo_libraries, depository=False)
    if rmg.solvent:
        rmg.database.load_solvation(os.path.join(path, 'solvation'))
        Species.solvent_data = rmg.database.solvation.get_solvent_data(rmg.solvent)
        Species.solvent_name = rmg.solvent
    for species in rmg.initial_species:
        submit(species)
    if library_flag:
        library = ThermoLibrary(name='Thermo Estimation Library')
        for species in rmg.initial_species:
            library.load_entry(
                index=len(library.entries) + 1,
                label=species.label,
                molecule=species.molecule[0].to_adjacency_list(),
                thermo=species.get_thermo_data().to_thermo_data(),
                shortDesc=species.get_thermo_data().comment,
            )
        library.save(os.path.join(rmg.output_directory, 'ThermoLibrary.py'))
    # Save the thermo data to chemkin format output files and dictionary, with no reactions    
    save_chemkin_file(os.path.join(rmg.output_directory, 'chem_annotated.inp'), species=rmg.initial_species, reactions=[])
    save_species_dictionary(os.path.join(rmg.output_directory, 'species_dictionary.txt'), species=rmg.initial_species)