Asymmetric nuclear matter at zero temperature is studied using a variational method which is an extension of the methods used by the present authors previously for simpler systems. An approximate expression for the energy per nucleon in asymmetric nuclear matter is derived through a combination of two procedures, one used for symmetric nuclear matter and the other for spin-polarized liquid 3 He with spin polarization replaced by isospin polarization. The approximate expression for the energy is obtained as a functional of various spin-isospin-dependent radial distribution functions, tensor distribution functions, and spin-orbit distribution functions. The Euler-Lagrange equations are derived to minimize this approximate expression for the energy; they consist of 16 coupled integrodifferential equations for various distribution functions. These equations were solved numerically for several values of the nucleon number density ρ and for many degrees of asymmetry ζ [ ζ = (ρ n − ρ p )/ρ, where ρ n (ρ p ) is the neutron (proton) number density]. Unexpectedly, we find that the energies at a fixed density cannot be represented by a power series in ζ 2 . A new energy term, ε 1 (ζ 2 + ζ 2 0 ) 1/2 , where ζ 0 is a small number and ε 1 is a positive coefficient, is proposed. It is shown that if the power series is supplemented with this new term, it reproduces the energies obtained by variational calculations very accurately. This new term is studied in relation to cluster formation in nuclear matter, and some mention is made of a possible similar term in the mass formula for finite nuclei. §1. IntroductionThe variational study 1) of infinite nuclear matter has a long history, and the results have been compared with other many-body calculations, such as the BruecknerHartree-Fock calculations. For asymmetric nuclear matter, which is the main object of our present study, there are fairly many studies employing nonrelativistic and relativistic Brueckner theories. 2) However, most variational calculations are for symmetric nuclear matter and neutron matter. 3) The variational calculations for asymmetric nuclear matter of which the author is aware are very few: the Fermi hypernetted chain (FHNC) calculations carried out in 1981, 4) and the lowest-order constrained variational (LOCV) calculations. 5) Since 1994, the present authors have been developing a different type of variational theory, in which approximate expressions for the energies per nucleon are constructed as functionals of various two-body distribution functions and structure functions. We now have approximate energy expressions for liquid 3 He, 6) neutron matter, and symmetric nuclear matter. 7) For liquid 3 He, we have an expression valid for arbitrary spin-polarized states. 8) In this paper, we extend the theory to treat asymmetric nuclear matter, utilizing the analogy between spin polarization and isospin polarization. Here, we should mention the guiding principle of these varia-