                 GGA = PE                    # GGA of Perdew, Burke and Ernzerhof 1997.
               ISPIN = 1                     # Controls whether the calculation is spin-polarized (2) or not (1)
               LASPH = .True.                # Include non-spherical contributions from the gradient corrections inside the PAW spheres
               ENCUT = 500.00                # The energy cutoff of the plane-wave basis (eV)
              ISTART = 0                     # Whether to read WAVECAR: 0: from scratch; 1: constant E basis; 2: constant basis.
                ALGO = Fast                  # The electronic minimization algorithm and/or type of GW calculation
             ISEARCH = 1                     # 0: use legacy linesearch; 1 for new.
                NELM = 60                    # The maximum number of electronic self-consistency steps
              NELMIN = 2                     # The minimum number of electronic self-consistency steps. 2 is the default, 4-8 recommended for minimization and MD
               EDIFF = 1.00E-05              # Convergence criterion for electronic self-consistency
                PREC = normal                # The precision of the FFT grids the real-space projections.
              ISMEAR = 2                     # the methfessel-paxton method
               SIGMA = 0.05                  # The width of the electron smearing (eV)
              IBRION = -1                    # How to change the structure during the calculation
                ISIF = 1                     # Determines whether the stress is calculated, and degrees of freedom.
                 NSW = 0                     # The maximum number of ionic steps
              EFERMI = MIDGAP                # Defines how the Fermi energy is calculated, or gives the value
                 LH5 = .False.               # Use HDF5 files instead of text files
              LORBIT = 11                    # Project the Kohn-Sham orbitals onto local quantum number (lm) to get on-site charge and magnetic moments
               NCORE = 2                     # Number of compute cores that work on each orbital
                KPAR = 2                     # Number of compute cores that work on each k-point
              LPLANE = .True.                # Use plane-wise data distribution in real space
               LREAL = .False.               # Evaluate the projection operators in real space (Auto) or reciprocal (.False.)
                NSIM = 4                     # Number of bands to optimize simultaneously when using RMM-DIIS orthe blocked Davidson algorithm (NSIM)
          LSCALAPACK = .True.                # Use ScaLAPACK for orthonormalization and subspace rotations.
              LSCALU = .False.               # Use ScaLAPACK for the subspace rotations.