The Cholesky decomposition (CD) technique, used to approximate the two-electron repulsion integrals (ERIs), is applied to the valence bond self-consistent field (VBSCF) method. Test calculations on ethylene, C2 n H2 n +2 , and C2 n H4 n -2 molecules (n = 1-7) show that the performance of the VBSCF method is much improved using the CD technique, and thus, the integral transformation from basis functions to VB orbitals is no longer the bottleneck in VBSCF calculations. The errors of the CD-based ERIs and of the total energy are controlled by the CD threshold, for which a value of 10(-6) ensures to control the total energy error within 10(-6) Hartree. © 2016 Wiley Periodicals, Inc.