Coulson (ZINDO), Mulliken (MP2/6-31G*) and Natural (MP2/6-31G*) population analyses of several large molecules were performed by the Fragment Reassociation (FR) method. The agreement between the conventional ZINDO (or conventional MP2) and FR-ZINDO (or FR-MP2) charges of these molecules was excellent. The standard deviations of the FR-ZINDO net atomic charges from the conventional ZINDO net atomic charges were 0.0008 for C 10 H 22 (32 atoms), 0.0012 for NH 2-C 16 O 2 H 28-COOH (53 atoms), 0.0014 for NH 3 +-C 16 O 2 H 28-COOH (54 atoms), 0.0017 for NH 2-C 16 O 2 H 28-COO − (52 atoms), 0.0019 for NH 3 +-C 16 O 2 H 28-COO − (53 atoms), 0.0024 for a conjugated model (O=CH-(CH=CH) 15-C=O-(CH=CH) 12-CH=CH 2), 118 atoms), 0.0038 for aglycoristocetin (C 60 N 7 O 19 H 52 + , 138 atoms), 0.0023 for a polypropylene model complexed with a zirconocene catalyst (C 68 H 121 Zr + , 190 atoms) and 0.0013 for magainin (C 112 N 29 O 28 SH 177 , 347 atoms), respectively. The standard deviations of the FR-MP2 Mulliken (or Natural) partial atomic charges from the conventional ones were 0.0016 (or 0.0016) for C 10 H 22 , 0.0019 (or 0.0018) for NH 2-C 16 O 2 H 28-COOH and 0.0033 (or 0.0023) for NH 3 +-C 16 O 2 H 28-COO − , respectively. These errors were attributed to the shape of molecules, the choice of fragments and the degree of ionic characters of molecules as well as the choice of methods. The CPU time of aglycoristocetin, conjugated model, polypropylene model complexed with zirconocene and magainin computed by the FR-ZINDO method was respectively 2, 4, 6 and 21 times faster than that by the normal ZINDO method. The CPU time of NH 2-C 16 O 2 H 28-COOH and NH 3 +-C 16 O 2 H 28-COO − computed by the FR-MP2 method was, respectively, 6 and 20 times faster than that by the normal MP2 method. The largest molecule calculated by the FR-ZINDO method was B-DNA (766 atoms). These results will enable us to compute atomic charges of huge molecules near future.