gdirect ! integral-direct mode basis={ default,aug-cc-pwCVTZ ! orbital basis set,ri;default,aug-cc-pwCVQZ/mp2fit ! RI BASIS set,oep;default,aug-cc-pVDZ/mp2fit ! OEP basis set,dfit;default,aug-cc-pwCV5Z/mp2fit ! density fitting basis } symmetry,nosym ! USCRPA does not use symmetry angstrom geometry={ 1 Li 0.0 0.0 0.0 } spin=1 df-uhf,maxit=0,df_basis=dfit ! HF calculation with 0 iteration {cfit,basis_coul=dfit,basis_exch=dfit} acfd;uscrpa,thr_fai_oep=5d-2,plot_z=1 ! USCRPA calculation