A novel procedure to generate initial broken-symmetry solutions is proposed. Conventional methods for the initial broken-symmetry solutions are the MO alter, HOMO-LUMO mixing and fragment methods. These procedures, however, are quite complex. Our new approach is efficient, automatic and highly practical especially for large QM systems. This approach, called the LNO method, is applied to the following four typical open-shell systems: H-2, dicarbene and two iron-sulfur clusters of Rieske-type [2Fe-2S] and [4Fe-4S]. The performance and the efficiency as an electronic structural analysis are discussed. The LNO method will be applicable for general systems in the complicated broken symmetry states. (c) 2014 Elsevier B.V. All rights reserved.