We propose an effective method to optimize the working parameters (WPs) of microwave-driven quantum gates implemented with multilevel qubits. We show that by treating transitions between each pair of levels independently, intrinsic gate errors due primarily to population leakage to undesired states can be determined by spectroscopic properties of the qubits and minimized by choosing proper WPs. The validity and efficiency of the approach are demonstrated by applying it to optimize the WPs of two coupled rf SQUID flux qubits for controlled-NOT operation. The result of this independent transition approximation (ITA) is in good agreement with that of dynamic method (DM). The ratio of the speed of ITA to that of DM scales exponentially as 2 n when the number of qubits n increases. DOI: 10.1103/PhysRevLett.95.120501 PACS numbers: 03.67.Lx, 85.25.Dq, 89.70.+c A practical quantum computer would be comprised of a large number of coupled qubits which must be kept in highdegree quantum coherence states for a sufficiently long time. During the past decade, significant progress has been made on quantum computation. High-degree quantum coherence has been demonstrated experimentally in systems such as trapped ions [1,2], nuclear spins [3,4], atoms in optical resonators [5], and photons in microwave cavities [6]. However, it seems quite difficult to realize a large number of coupled qubits using these systems. Meanwhile, solid-state qubits are of particular interest because of their advantages of large-scale integration, flexibility in design, and easy connection to conventional electronic circuits [7]. Of those, qubits based on superconducting devices have recently attracted much attention [8] as manipulation of quantum coherent states has been successfully demonstrated in a variety of single qubits [9-13] and coupled two-qubit systems [14 -18].However, the solid-state qubits demonstrated in experiments so far all have relatively short coherence time and high probability of gate errors [9][10][11][12][13][14][15]17]. One of the causes of these problems is extrinsic gate error arising from interaction between the environment and qubits resulting in decoherence, such as dephasing and relaxation [8]. Another cause is intrinsic gate error resulting from population leakage to undesired states due to the typical multilevel structures of solid-state qubits [19,20]. The intrinsic gate error is crucial since it determines the ultimate performance of the quantum gates and cannot be eliminated by reducing the environmentally caused decoherence.The leakage of a gate can be characterized quantitatively by summing up the maximum transition probabilities to all undesired states of a multilevel qubit [19,20]. It depends strongly on energy level structure and transition matrix elements, i.e., spectroscopic properties of the multilevel qubit determined completely by device parameters (DPs) and external control parameters of the qubit, which we call working parameters (WPs) for simplicity. For instance, inductance (capacitance) of the super...