A computer-aided molecular design method (CAMD) based on the higher-order group contribution method (GC+) and the conductor like screening model-segment activity coefficient (COSMO-SAC) model is proposed. First, based on the GC+ method and the COSMO-SAC model, the GC+-COSMO method is constructed, which correlates the molecular group combination and the surface screening charge density distribution [σ-profiles, p(σ)] as well as the volume of molecular cavity (Vc), so as to achieve high-throughput prediction of them. Then, the simplified molecular input line entry system (SMILES) based isomer generation algorithm and the GC+-COSMO method are combined to realize the recognition and property differentiation of solvent molecular isomers by CAMD technology. Finally, the solvent design problem is established by the mixed integer nonlinear programming model (MINLP) composed of objective function and constraint equations, and the decomposition algorithm is further used to optimize the solution to achieve the solvent optimization design goal. Based on the above models and methods, the Diels-Alder (DA) competitive reaction solvent design case is carried out, and the feasibility and effectiveness of the proposed method are verified. © 2021, Chemical Industry Press Co., Ltd.. All rights reserved.