def infos_on_system(self):
        """
        Returns infos on initial parameters as in the log file of Fiesta
        """

        o = []
        o.append("=========================================")
        o.append("Reading infos on system:")
        o.append("")
        o.append(" Number of atoms = {} ; number of species = {}".format(
            int(self._mol.composition.num_atoms), len(self._mol.symbol_set)))
        o.append(" Number of valence bands = {}".format(
            int(self._mol.nelectrons / 2)))
        o.append(" Sigma grid specs: n_grid = {} ;  dE_grid = {} (eV)".format(
            self.correlation_grid['n_grid'], self.correlation_grid['dE_grid']))
        if int(self.Exc_DFT_option['rdVxcpsi']) == 1:
            o.append(" Exchange and correlation energy read from Vxcpsi.mat")
        elif int(self.Exc_DFT_option['rdVxcpsi']) == 0:
            o.append(" Exchange and correlation energy re-computed")

        if self.COHSEX_options['eigMethod'] == "C":
            o.append(
                " Correcting  {} valence bands and   {} conduction bands at COHSEX level".format(
                    self.COHSEX_options['nv_cohsex'],
                    self.COHSEX_options['nc_cohsex']))
            o.append(" Performing   {} diagonal COHSEX iterations".format(
                self.COHSEX_options['nit_cohsex']))
        elif self.COHSEX_options['eigMethod'] == "HF":
            o.append(
                " Correcting  {} valence bands and   {} conduction bands at HF level".format(
                    self.COHSEX_options['nv_cohsex'],
                    self.COHSEX_options['nc_cohsex']))
            o.append(" Performing   {} diagonal HF iterations".format(
                self.COHSEX_options['nit_cohsex']))

        o.append(" Using resolution of identity : {}".format(
            self.COHSEX_options['resMethod']))
        o.append(
            " Correcting  {} valence bands and  {} conduction bands at GW level".format(
                self.GW_options['nv_corr'], self.GW_options['nc_corr']))
        o.append(
            " Performing   {} GW iterations".format(self.GW_options['nit_gw']))

        if int(self.BSE_TDDFT_options['do_bse']) == 1:
            o.append(" Dumping data for BSE treatment")

        if int(self.BSE_TDDFT_options['do_tddft']) == 1:
            o.append(" Dumping data for TD-DFT treatment")
        o.append("")
        o.append(" Atoms in cell cartesian A:")
        symbols = []
        for syb in self._mol.symbol_set:
            symbols.append(syb)

        for site in self._mol:
            o.append(" {} {} {} {}".format(site.x, site.y,
                                           site.z, int(
                    symbols.index(site.specie.symbol)) + 1))

        o.append("=========================================")

        return str(o)