We present a new scalar relativistic formulation for the full-potential linear-combinationof-atomic-orbitals method based on the density-functional theory. Three approximations are introduced to overcome computational difficulty. The first is to consider only the large component of the four-component spinor, neglecting the small component. The second is to neglect the energy dependence in the Hamiltonian reduced for the large component. The third is to replace the material-dependent potential with the atomic potential in relativistic corrections. After the three approximations, we identify the mass-velocity and Darwin terms and also the spin-orbit coupling, where the latter is to be omitted according to the definition of the scalar relativistic formulation. The computational effort of the present method is reduced considerably in comparison with that of the fully relativistic method, being almost the same as that of the nonrelativistic method. We apply the present method within the local-density approximation to several diatomic molecules with heavy elements, crystalline Au, and crystalline InSb. The results are improved considerably in comparison with the nonrelativistic results. The calculated structural properties are in good agreement with the fully relativistic results and also with the experimental results. The calculated electronic properties are also improved considerably in comparison with the nonrelativistic results and are also in good agreement with the fully relativistic results except for the effect due to the spin-orbit coupling.KEYWORDS: scalar relativistic calculations, full-potential calculations, LCAO method, density-functional theory, band calculations, structure optimizations 1 §1. IntroductionThe structural and electronic properties of materials with heavy elements are strongly affected by relativistic effects. Core electrons experience strong nuclear field and, consequently, the spatial extent of the core-electron orbitals is contracted substantially. As a result, the efficient screening of nuclear charge is caused by the core electrons, inducing the expansion of valence-electron orbitals, especially of d and f orbitals. 1) This effect, which is known as the indirect relativistic effect, strongly affects the structural properties as well as the electronic properties.2) The mass-velocity and Darwin terms are responsible for the indirect relativistic effect. Furthermore, of another importance is the spin-orbit coupling. For example, the splitting due to the spin-orbit coupling in the conduction and valence bands in crystals with heavy elements sometimes has the same order of magnitude as the crystal-field splitting.The strategy which treats the mass-velocity and Darwin terms but neglects the spin-orbit coupling is known as the scalar relativistic method. The method has a practical advantage that the computational effort is reduced considerably in comparison with the fully relativistic method, which considers not only the mass-velocity and Darwin terms but also the spin-orbit coupling. So far...