YouTube video by Mike Ladd - Topic

Cu3BiS3 (CBS) has emerged as a promising earth-abundant absorber for thin-film photovoltaics, offering a sustainable alternative to conventional technologies. However, ab initio studies on its optoelectronic properties remain scarce and often yield contradictory results. This study systematically examines the influence of two density functional theory (DFT) methodologies, linear combination of atomic orbitals (LCAO) and projector augmented wave (PAW), on the structural and electronic properties of CBS, aiming to establish a reliable computational framework for future research. With this in mind, we also assessed the impact of a wide range of exchange-correlation (XC) functionals within both methods, including 6 from the local density approximation (LDA) family (HL, PW, PZ, RPA, Wigner, XA), 10 from the generalized gradient approximation (GGA) family (BLYP, BP86, BPW91, GAM, KT2, PBE, PBEsol, PW91, RPBE, XLYP), 2 meta-GGA functionals (SCAN, R2SCAN), and the hybrid HSE06 functional. Both LCAO and PAW consistently predict an indirect bandgap for CBS across all XC functionals, aligning with most previous DFT studies but contradicting experimental reports of a direct transition. The LDA and meta-GGA functionals systematically underestimated the CBS bandgap (<1 eV), with further reductions upon structural relaxation. GGA functionals performed better, with BLYP and XLYP yielding the most experimentally consistent results. The hybrid HSE06 functional substantially overestimated the bandgap (1.9 eV), with minimal changes after relaxation. The calculated hole and electron effective masses reveal strong anisotropy along the X, Y, and Z crystallographic directions. Additionally, CBS exhibits an intrinsic p-type nature, as the Fermi level consistently lies closer to the valence band maximum across all methods and functionals. However, the PAW method generally predicted more accurate lattice parameters than LCAO; the best agreement with experimental values was achieved using the PW91 (1.2% deviation) and HSE06 (0.9% deviation) functionals within LCAO. Based on these findings, we recommend the PW91 functional with LCAO for structural optimizations in large supercell studies of CBS dopants and/or defects and BLYP/XLYP for electronic properties.

arXiv abstract link

arXiv abstract link

We present a domain-based local pair natural orbital Møller--Plesset second order perturbation theory (DLPNO-MP2) for periodic systems, working within an LCAO formalism within the Tubromole program package. This approach, Megacell-DLPNO-MP2, embeds a supercell correlation treatment within a megacell and does not involve periodic image summation for the Coulomb integrals. Working in a basis of well-localised Wannier functions, periodicity is instead imposed through rigorous translational symmetry of Hamiltonian integrals and wavefunction parameters. The accuracy of the method is validated through comparison with a complementary periodic DLPNO-MP2 method that employs Born--von K{á}rm{á}n boundary conditions, described in paper I of this series. The PNO approximations are shown to be equivalent in the two approaches and entirely consistent with molecular DLPNO-MP2 calculations. The Megacell-DLPNO-MP2 method displays sub-linear scaling with respect to supercell size at the asymptotic limit and example calculations are presented with up to 15000 basis functions in the correlation treatment.