A relativistic low-energy electron-diffraction theory has been modified to allow computations for energies up to 200 eV. Its application to W(001) yields intensity profiles in reasonable correspondence with experiment, and pronounced polarization profiles which show encouraging agreement with novel and as yet limited experimental data. A high sensitivity of the polarization profiles to displacement of the topmost layer suggests that spin polarization measurements could be valuable for structure determination of surfaces involving heavy atoms.Theoretical studies of spin-polarization effects in low-energy electron diffraction (LEED) from the (001) surface of tungsten 1 * 2 were recently followed by an experiment, 3 in which significant degrees of spin polarization were found for electron energies from 45 to 190 eV at various small angles of incidence. Comparison of these novel data with the theoretical predictions 1 ' 2 is not possible, however, since the latter are confined to energies below about 40 eV. It is the aim of this Letter to report the first theoretical results at higher energies, to compare them with the recent experiment, 3 and to draw conclusions regarding the value of a polarization analysis in surface-structure determination.The theory underlying the present calculations follows to a large extent a relativistic LEED theory described earlier. 4 In particular, the crystal is assumed to consist of a finite number of monatomic layers, for which the Dirac-equation boundary problem is solved to yield the fourspinor amplitudes and thence the intensities and the spin polarization vectors of the reflected beams as functions of the energy, the polar and azimuthal angles of incidence, and the polarization vector of the (normalized) primary beam. A modification was, however, made in the treatment of the monolayer part of the problem. The Korringa-Kohn-Rostoker-Ziman (KKRZ) 5 type relativistic pseudopotential used in Ref. 4 depends on radii R t -associated with the spin-up and spin-down atomic phase shifts 5^-which have to be optimized individually in a numerical application of the method. Since calculations up to 200 eV require phase shifts up to about l-l (cf. also Van Hove and Tong), 6 I decided to avoid the resulting radii parameter problem by developing and applying the Dirac-equation-based analog 7 of a KKR-type Schrodinger-equation-based method, 8,9 which has been used with great success for LEED intensity calculations. 10 The ef-fect of thermal lattice vibrations is taken into account by replacing the actual (real) spin-up and spin-down atomic phase shifts 6j* by effective (complex) phase shifts fy*, which are obtained from the fy* by averaging the corresponding scattering amplitudes over a Debye spectrum. 11 In the computational application to W(001) the following specific model assumptions have been made. Effective phase shifts oV, obtained from Mattheiss's muffin-tin potential 12 and with use of the bulk Debye temperature of 380°K, are included up to 1 = 1. The number of surface reciprocal lattice ...