B3LYP is the most famous hybrid density functional theory model, which includes Hartree-Fock exchange, local exchange, gradient exchange correction, local correlation, and gradient correlation correction. Historically, the relative weight of each component in B3LYP, which is controlled by three empirical parameters (a 0 , a x , a c ), has not been optimized. In this work, we perform global optimization against accurate experimental reference, optimal empirical parameters, and the better version of B3LYP are obtained and denoted as OpB3LYP. The performance of OpB3LYP is widely tested over many species and chemical properties, the results show that the computational accuracy is significantly improved as compared to original B3LYP and the serious size dependence of B3LYP is remarkably overcome by the employment of OpB3LYP. The comparative assessment of OpB3LYP and other prevalent functionals indicates that OpB3LYP is a promising functional for large molecules.