We propose a torque method for the theoretical determination of the magnetocrystalline anisotropy ͑MCA͒ energy for systems with uniaxial symmetry. While the dependence of the total energy on the angle between the magnetization and the normal axis ( ) can be expressed as E( )ϭE 0 ϩK 2 sin 2 ( )ϩK 4 sin 4 ( ), we show that the MCA energy ͓defined as E MCA ϭE( ϭ90°)ϪE( ϭ0°)ϭK 2 ϩK 4 ͔ can be easily evaluated through the expectation value of the angular derivative of the spin-orbit coupling Hamiltonian ͑torque͒ at an angle of ϭ45 o . Unlike other procedures, the proposed method is independent of the validity of the MCA force theorem, or of the absolute accuracy of two total energy calculations. Calculated MCA energies for the free Fe monolayer with different lattice constants are analyzed and compared with results of other ab initio calculations, especially those obtained with our previously reported state tracking method. ͓S0163-1829͑96͒02926-8͔Recent developments in magnetic thin films and overlayers that show perpendicular magnetic anisotropy 1 with their tremendous practical implications for high-density magnetooptical storage media, have stimulated ab initio theoretical determinations of the magnetocrystalline anisotropy energy ͑MAE͒. 2-12 Unfortunately, despite the great advances in local-spin-density electronic structure theory and computational power in the past decades, the accurate determination of MAE still remains difficult and computationally demanding. In the traditional approach, the value of MAE is calculated through comparing the total energies of a given system for two different magnetization orientations ͑in-plane and perpendicular directions for surface-interface systems͒. Since the spin-orbit coupling ͑SOC͒ is very weak in 3d transition metals, the so-called magnetocrystalline anisotropy ͑MCA͒ force theorem 6 is usually adopted and the MAE is calculated by merely comparing the band energies between the two magnetic orientations. ͑We call them direct approaches.͒ The main difficulty associated with the direct approaches concerns the numerical stability of calculating a very small difference of two large numbers. In fact, in order to eliminate numerical fluctuations, [3][4][5][6][7]9 extremely fine sampling meshes are required for the k-space integrations.In this paper, we propose a torque method for the determination of MCA. The main advantage that this method offers is that one only needs to calculate MAE at one particular magnetic orientation and to do the k-space integration with the single Fermi surface at this orientation. In the following, we first present the method, then give results for free standing Fe monolayers as a test case, and compare them with the results of other methods. [3][4][5]9 To demonstrate the idea, let us recall that the total energy of an uniaxial system can be well approximated in the formwhere is the angle between the magnetization and the normal axis. If we define the torque, T( ), as the angular derivative of the total energy, i.e.,it is easy to show that M AEϵE͑ ϭ...