\section{Introduction}

 The presence of giant magneto-resistance effects in trilayer systems
consisting of two ferromagnetic metallic layers -- source and drain --
separated by a non-magnetic or anti-ferromagnetic crystalline metallic
spacer is meanwhile well-known \cite{Fert,Gruenberg}.
Recently also non-crystalline and non-metallic spacers
have been considered
 \cite{La1,La2,Buergler}, i.e.\ there is renewed  interest
in spin-polarized tunneling through semi-conducting or
non-conducting spacers \cite{Ted,Jul,Ste},
 based on  the perspective of new applications  in
magneto-electronic technology, e.g.\ for
magnetic-field sensors, spin-valve transistors \cite{Mon} or spin-polarized
field-effect transistors \cite{Da1,Pri}.

 In \cite{Slo,Mat2} spin-valve properties of
ferromagnet/insulator/ferromagnet junctions were studied, and
\cite{But} presents a first $ab$ $initio$
 calculation of the {\it  one-particle} electronic and magnetic properties of FM/NM/FM
tunneling structures, with
FM=Fe, and NM=Ge and GaAs.

However concerning the GMR effect, which is based on {\it two-particle}
  properties, equally accurate calculations for FM/NM/FM systems,
particularly in the CPP geometry, apparently could only be performed without
impurities, \cite{Butler2}. For systems with {\it metallic}
 spacers one should mention at
the same time the calculations of Tsymbal and Pettifor, \cite{Tsymbal}, or
Mathon, \cite{Mathon97}, where also fully realistic calculations have been performed,
but again with phenomenological assumptions for the impurity scattering,
\cite{Tsymbal}, or with no impurities at all, \cite{Mathon97}. Furthermore,
in a two-band tight-binding approximation, two of the present authors,
\cite{KroKrey98}, have already treated systems, where one or two of the
ferromagnet-spacer interfaces were decorated with ultrathin non-metallic layers,
but again with translational invariance within the planes.
 In this approximation, the two bands mentioned above describe the s- and
d-electrons and their hybridization, or the conduction and valence band in
case of a non-metallic spacer.

 In the ferromagnetic metals, the parameters
of the d-band are of course spin-dependent.

Since in \cite{KroKrey98} the method for the calculation of the
resistivities turned out to be extremely flexible and accurate, we extend it
in the  present
paper to a systematic survey of the magnetoresistivity of
trilayers of two ferromagnetic sandwiches separated by different
non-magnetic spacers, metallic and non-metallic ones, with impurities of
various kind, and
without impurities, in CIP and and CPP geometries, but within a
{{simplified  semi-empirical hybridized two-band
model}}.

We stress that our approach is rigorous in principle,   numerically very
 accurate, and avoids the Coherent Potential Approximation (CPA),
 which is often applied to disordered systems, but difficult to control.

As already mentioned, there exists already a large body of theoretical
calculations on the GMR of magnetic multilayer systems.
For systems with metallic spacers and without impurities  our simplified two-band model
calculations, to be presented
below, can of course not { compete}  -- and this is  also not
at all intended -- with the much more ambitious and cumbersome calculations
mentioned above   e.g.\ those of \cite{Mathon97}, who performed fully
realistic calculations of the CPP- and CIP-GMR for Co/Cu/Co/Cu and
Fe/Cr/Fe/Cr multilayers with ideal leads made of Cu rsp.\ Cr, in the
ballistic, Ohmic, and Anderson-localized regimes corresponding to the cases of
 no-disorder rsp.\ certain cases of one-dimensional disorder, i.e.\ with
respect to the stacking of the planes, but without impurities. There exist
more such ambitious calculations, e.g.\ the seminal paper of Schep {\it et
al.}, \cite{Schep}, and other publications, which also in principle would
have earned citing but which are omitted on lack of space. In this respect
the purpose of our paper is, instead, to study in some detail the question,
{\it to which extent} simplified models as the present one are able (or not)
to reproduce those results.

 On the other hand, with impurities, on which we concentrate,
similar extensive calculations
{{of corresponding rigour}
do not yet exist, to our knowledge. In fact,
concerning the influence of impurities, the
{{\it methods}}
 of our calculation have
already been presented in the paper of Asano {\it et al},
\cite{Asano}. However, rigorous calculations with impurities are so
demanding that those authors, although they  wrote already down the formalism for a
{{\it two-band model}}, published numerical results only for
{{\it one-band}} cases.

Here we present results obtained with the two-band model within the
formalism of \cite{Asano}. For a {\it nonmetallic} spacer it is actually necessary
to consider at least two orbitals per atom to simulate the valence and
conduction bands, respectively, but also for ferromagnetic or non-magnetic
3d-{\it metals} it is a natural requirement to include a second band with the
intention to mimic as far as possible the d-electrons in addition to the
s-states. Namely, the d-electrons are not only responsible for the magnetic
properties, but also for a large part of the resistivity, due to
{{\it  Mott's
$s$-$d$-scattering mechanism}}. Therefore it is also important to include not
only two bands, but also realistic values for the $s$-$d$ {\it
hybridization}.

This should give enough motivation to our two-band calculation; but  we
stress again that these calculations should be considered as model
calculations and are not intended to replace fully realistic calculations,
as far as they are possible already now, e.g.\ \cite{Mathon97}, or in future.
In fact, it turns out that the present two-band model calculations will be
able to give estimates of the GMR, which yield the correct order of
magnitude, but not more. Further, concerning the impurity effects, it turns
out that certain contradictory results (see below) concerning essential
trends depend strongly on the model-parameters, so that again, now with
impurities, completely realistic calculations cannot be avoided: Unfortunately, such
 calculations, where the impurities are treated with the same kind
of rigour as the electronic structure of the pure system itself, are
apparently not yet possible,
{{\cite{Mertig1}--\cite{Mertig4}}}.

 On the other hand, one of our most remarkable predictions, which should not
be overlooked and can be stated already from the present model
calculations, is that a large enhancement of the CPP-GMR may happen, if
{\it nonmetallic spacers} of increasing thickness are used. Here it is
{{of course}}
required that the thickness of the nonmetallic spacer remains smaller than
the dephasing length for inelastic scattering; i.e.\ one should remain in
the ballistic limit.

In the following sections we present the basic theory and then the results
of our\, survey, which is followed at the end
 by a section presenting our conclusions.
\section{Basic definitions}
\subsection{Systems}
Our systems are defined in Fig.~1 and Fig.~2 for the CPP and CIP geometries,
respectively.

We start with metallic 'ideal leads' on the left-hand side and
right-hand-side, respectively,  followed  -- in the CPP geometry -- by three monolayers
of ferromagnetic metal $F_1$ rsp.\ $F_2$, with $n_s$ spacer monolayers
inbetween.
In the CPP geometry all layers have a quadratic cross-section of $M\times M$
square-lattice unit cells.

For the CIP geometry, the natural definition of our systems is similar and
explained by Fig.\ 2.

The electronic structure of our systems is described by the following
tight-binding Hamiltonian\ :
\begin{eqnarray*}
 {\cal H}& = &\sum_{l,i}\,E_{l.i}^s\,\sum_{\sigma}
            \hat c^+_{l,i,\sigma}\hat c_{l,i,\sigma}
          +\sum_{l,i}\sum_{l',i'}t^s_{l,i,l',i'}
           \sum_\sigma\hat c^+_{l,i,\sigma}\hat c_{l',i',\sigma} \\
         & + &\sum_{l,i}\sum_{\sigma}\,E^d_{l,i,\sigma}\,
            \hat d^+_{l,i,\sigma}\hat d_{l,i,\sigma}
          +\sum_{l,i}\sum_{l',i'}t^d_{l,i,l',i'}
           \sum_\sigma\hat d^+_{l,i,\sigma}\hat d_{l',i',\sigma} \\
         & + &\sum_{li}\,V^{s,d}_{l,i}\,\sum_{\sigma}
            (\hat d^+_{l,i,\sigma}\hat c_{l,i,\sigma}+\hat c^+_{l,i,\sigma}
\hat d_{l,i,\sigma})\,\,.
\end{eqnarray*}

Here $\hat c^+_{l,i,\sigma}$ and $\hat c_{l,i,\sigma}$  denote creation and
destruction operators of an s-electron occupying the site $l$ in plane $i$
with spin $\sigma$. The corresponding operators for d-electrons are
$\hat d^+_{l,i,\sigma}$ and $\hat d_{l,i,\sigma}$. The
band-energy parameters are $E_{l,i}^s$ for the s-states, which do not depend on
the spin $\sigma$, and $E^d_{l,i,\sigma}$,
which are spin-dependent, for the d-states. The hopping-matrix elements for s-
and d-states  are
$t^s_{l,i,l',i'}$ and $t^d_{l,i,l',i'}$, and the local s-d-hybridization is
given by $ V^{s,d}_{l,i}$.

  Of course, taking more bands into account,
 i.e.\ {{all}} five d-bands, {{ the 4s-band and the
three 4p-}}{{bands, and adding the spin dependence}},
we could have tried  to fit the band
structure of realistic systems as far as possible, see e.g. \cite{Sanvito};
but this is
{{\it not}} our purpose, since we try to remain  semi-quantitative and study
only the essential trends.
 Particularly, in view of the still much simplified band structure, we also
use typical simplifications in the description of the impurities, see below.
But in principle, within our formalism, more rigour in both respects
would be possible, but only at the cost of excessive computing.

For the following, $\cal H$ is written in matrix representation with respect
 to an orthonormal basis by replacing
the creation and annihilation operators by the corresponding
ket- and bra-vectors, respectively, such that ${\cal H}\to \hat
H$, e.g.\ by $\hat d^+_{l,i,\sigma}\hat c_{
l,i\sigma}\to \ket{d,l,i,\sigma}\bra{s,l,i,\sigma}$.
\subsection{Resolvent operators and Kubo formula}
For the calculation of the conductances we use the formalism of
Fisher and Lee, \cite{Fis}. Therefore, we need the following
matrix elements of the advanced and retarded resolvent operators $\hat
 G^{\pm} :=\{ E \pm {\rm i}\,0^+ -{\hat H}\}^{-1}$~:
\begin{equation}
G^\pm(i,i')_{{\alpha},l,\sigma\,;\,{\alpha}',l',\sigma'}=
\bra{{\alpha},i,l,\sigma}\{E\pm {\rm i}\,0^+ -{\hat H} \}^{-1}
\ket{{\alpha}',i',l',\sigma'}\,\,.
\end{equation}
Here the index ${\alpha}$ counts the orbitals, e.g.~${\alpha} = s$ or $d$.
As usual, we also call these matrices 'Green's operators', and ${\rm
i}\,0^+$
is a positive  imaginary infinitesimal,
such that the matrix-inverse in this
equation always  exists in our approach.

 For auxiliary purposes we need the 'left-sided' and
 'right-sided' Green's operators $G^L$ and $G^R$, e.g.\
  \begin{equation}
G^L(i_0)_{{\alpha},l,\sigma\,;\,{\alpha}',l',\sigma'}\, :=\,
\bra{{\alpha},i_0,l,\sigma}\{E+{\rm i}\,0^+-\hat H_L^{i_0}\}^{-1}
\ket{{\alpha}',i_0,l',\sigma'}
\,\,.\end{equation}
Here $\hat H_L^{i_0}$ denotes the Hamiltonian for a system, in which all
planes $i > i_0$ are deleted; $G^R(i_0)$ is defined in a similar way with
$E+{\rm i}\,0^+-\hat H_R^{i_0}$, where $ H_R^{i_0}$
 is the Hamiltonian of a system with all
planes $i < i_0$  deleted.
 $G^L$ and $G^R$ are obtained by {\it recursion}, e.g.\ $G^L$ by recursion
 from the left, i.e.\
\begin{equation}
G^L(i_0)=[\hat g(i_0)^{-1}-\hat TG^L(i_0-1)\hat T]^{-1}\,\,.
\end{equation}
$\hat g(i_0)$ is  the resolvent of the isolated plane $i_0$; $\hat T$ is
the matrix describing the 'hopping' from plane to plane.

 From this one gets finally the two desired resolvent operators
\begin{equation}
G(i,i)=[\hat g(i)^{-1}-\hat TG^L(i-1)\hat T-\hat TG^R(i+1)
\hat T]^{-1}
\end{equation}
and
\begin{equation}
G(i,i+1)=G(i,i)\hat TG^R(i+1)\,,
\end{equation}
which are needed for the conductance, see below.
Here the plane $i$ ($=i_0$) is arbitrary, since under stationary
 conditions through any of the
planes the same amount of current is flowing.

Finally, with the definition
\begin{equation}
\tilde G(i,i')\,:=\frac{1}{2{\rm i}}[G^-(i,i')-G^+(i,i')]
\end{equation}

one gets the conductance from the
{ Kubo-formula}, \cite{Fis}:

\begin{equation}\label{eqKubo}
\Gamma_\sigma=\frac{4e^2}{h}{\rm Tr}_{l\,,\,\alpha }
\,[\tilde G(i,i)\hat T\tilde G(i+1,i+1)
{\hat T}
-{\hat T}\tilde G(i,i-1)\hat T\tilde G(i,i-1)]\,,
\end{equation}
where the trace ${\rm Tr}_{l\,,\,\alpha }$
is performed with respect to the orbitals ${\alpha} =s,d$,
and with respect to the sites $l$ of the atoms of a given plane, whereas
the spin-index $\sigma$  ($=\sigma'$) is kept fixed. Here $e$ is the
electronic charge, and $h$ is Planck's constant.

 The GMR effect is  defined by the equation
\begin{equation} {\rm GMR}\,=\,
\frac{\Gamma^{++}_{\sigma=\uparrow}+\Gamma^{++}_{\sigma=\downarrow}}
{\Gamma^{+-}_{\sigma=\uparrow}+\Gamma^{+-}_{\sigma=\downarrow}} -1
\,=\,\frac{\Gamma^{++}}
{\Gamma^{+-}} -1\,\,.
\end{equation}
Here the superscripts ($++$) resp.\ ($+-$) denote
{{\it
mutually parallel}} rsp.\
{{\it
antiparallel}} magnetizations of the two ferromagnetic sandwiches. Since
 we neglect spin-flip scattering altogether, the conductances are of course
additively composed out of the separate contributions from majority and minority
spins, respectively, as denoted in the formula.

We have performed real-space (x-space) and k-space calculations. In the
x-space calculations each layer consists of $M\times M=$10$\times$10 atoms
with free boundary conditions, {{\cite{REMrealspace}}}.
 Here the matrices are of size 200$\times$200, again for given
$\sigma$ ($=\sigma'$). Impurities can be added at will, which is important
below.

In our $k$-space calculations, which are only performed in the CPP case
and without impurities, the planes are infinitely extended.
 Here, in the trace of equation (\ref{eqKubo}), in the $k$-space
 representation, the sum over $l$ is replaced by the corresponding
  sum over $\vec k$-vectors in the two-dimensional Brillouin zone,
since the two-dimensional vector $\vec {k}_{||}$ is now 'a good quantum number', i.e.\
there is translational invariance in the planes.   This summation, or
the corresponding integral, is approximated extremely
accurately by a Cunningham formula, \cite{Cun}, with $\sim 10^6$
$\vec{k}_{||}$-points, for which the matrix elements of the Green's operators
can be calculated separately without difficulty and extremely accurately. (All
numerical evaluations have been performed with  Cunningham's accuracy parameter
$m=10$ or $11$).
\subsection{Model parameters}
For the models considered, we have chosen four parameter sets: Two parameter
sets, (i) and (ii), correspond to a {\it metallic} spacer and mimic the
main features of trilayers composed of Co as ferromagnetic metals and Cu as
nonmagnetic metallic spacer material (see below), whereas the third and
fourth parameter sets, (iii) and (iv), correspond again to Co ferromagnets, but
with a {\it non-metallic} spacer.

Concerning the  sets (i) and (ii), we note that, as already mentioned, they are
constructed to reflect the two main features of Co and Cu as far as
possible in our simplified two-band model: The two features are
\begin{itemize} \item the similarity of the s-bands of Co and Cu, and
\item the similarity of the d-bands of Cu and the 'majority-spin' d-bands of
Co.
\end{itemize}
This  qualitative
similarity is obvious from  Fig.~3, which presents
results obtained by Mathon {\it et al.}, \cite{Mat3}.
 The parameter sets (i) and (ii) are in fact derived from rough
 fits to {\it the overall band-structure of Co} --
parameter set (i), see Fig.\ 4a --  respectively
 to {\it the overall structure of the DOS for Co} -- parameter set (ii), see
Fig.\ 4b.
But these  fits are rather crude anyway, and so the
derivation  is not emphasized. The values are given in
Table 1 and Table 2, respectively, \cite{REMgroupvelocity}.

 In Tables 3 and 4, the parameter sets (iii) and (iv)  are
given, which correspond
to {\it nonmetallic} spacers ('Isolator
 1 and 2', respectively). The parameters of the ferromagnet correspond to
case (i), with $E_F=-2.8$.

\noindent
The notions 'indirect energy gap'  and 'direct
energy gap', which are used in the table captions of Table 3 and 4,
respectively, refer to Fig.~5:

\noindent For parameter set (iii) there is an '{\it indirect energy gap}', since the
 valence-band maximum is at $k_{[111]}=(\pi/a)\sqrt{3}$ with an energy -2.81,
which is slightly below the Fermi energy $E_F= -2.80$, whereas the minimum
of the conduction-band energy is at $k_{[111]}=0$ with E=-2.79. These
numbers refer to the first line of the three parameter sets of Table 3.

 On the other hand, with parameter set (iv), also with $E_F=-2.8$, there is
 a {\it direct} energy gap at $k_{[111]}= 0.814$,
 i.e.\ at the {\it arrow} in the r.h.s.\ of Fig.\ 5. In both figures, and
in the results presented below, the uppermost line, i.e.~with the smallest
gap, of Table 3 and Table 4 has been used.

\subsection{Impurity models}
We consider a) bulk impurities and b) interface impurities.

In case a), we have considered {\it bulk impurities} only in the spacer,
{{and exclusively nonmagnetic ones, which}}
 are modelled by adding to the parameter $E_d$ of the single-site
d-band energy a {\it Gaussian spatial noise} proportional to a 'disorder
strength'  $\sigma_r$, namely
\begin{equation} \label{eqBulk_impurities} E_d\to E_d + \sigma_r\cdot n_l^{\rm Gauss}\,.\end{equation}
Here the $n_l^{\rm Gauss}$, where $l$ enumerates the different spacer atoms,
 are independent random numbers  distributed according to a 'Gaussian' with
average 0 and variance 1;
the parameters of the s-bands remain unchanged.

These assumptions for our bulk impurities are of course rather schematic,
i.e.\ these are the same simplified assumptions e.g.\ of Asano {\it et al.},
\cite{Asano}. But principally one might have been more ambitious and might
have considered e.g.\ specific magnetic impurities. Such calculations are of
course always possible within our formalism,
{{\cite{REMgeom,Garcia}}},
 and have partially been performed in
course of the diploma work of the first one of the present authors;
 however, as a first step, and since the results with the magnetic
impurities were not particularly exciting and have not been documented,
 also because of lack of space, we keep at present to our simplified
description of the impurities. This seems also justified in view of the
simplifications made for the band structure and for the geometry of our
systems.

Concerning case b), the {\it interface impurities} are defined as follows~:
For every ferromagnet-spacer interface, in each of the two boundary
layers of the interface there are $n$ atoms (out of $M\times M$) selected at
random and replaced by atoms of the other kind.

For impure systems we have of course {\it averaged}
 over a large number of different samples, see the figures below; this leads
to the 'averages with error bars' in the results presented in the following
section.

\section{Results}
\subsection{CPP-GMR}
 \subsubsection{CPP-GMR; no impurities}
 At first we describe the results obtained without impurities obtained
 in the CPP case ('current perpendicular to the planes'). In Fig.\ 6
and Fig.\ 7
the results of k-space calculations with model parameters
 (i) and (ii) are presented,
 for variable thickness $n_s$ of the 'metallic spacer'. Interestingly,
in both cases there is a slight but significant spatial oscillation
of the GMR as a function of $n_s$;  in case (i) the GMR oscillates between
$\sim 70$ and $\sim 80$ \% with a quasiperiod of $\lambda \approx 5$
 monolayers (ML); in case (ii)  the oscillation of the GMR is between $\sim 47$
and $\sim 51$ \% with $\lambda \approx 3$. Here one can explicitly see that
 the conductance is {\it ballistic} from the fact that on average the
conductances hardly depend on $n_s$ in the limit $n_s\to\infty$; in
contrast, in case of {\it Ohmic} behaviour, one would expect that the
conductances  converge to 0, namely $ \propto n_s^{-1}$, as soon as
$n_s$ becomes larger than the relevant ''scattering length'' for diffusive
elastic or inelastic scattering, \cite{Mathon97}.

The reason for the spatially oscillating behaviour of the GMR are tiny
oscillations of the conductances $\Gamma^{++}_\uparrow$ and
$\Gamma^{++}_\downarrow$, which are visible in the plots on the r.h.s.\
 of the preceding figures. These oscillations are caused by our treatment of
the ideal leads: In fact, for the ideal leads we have always used the same
parameters as within the 'Cu' spacer, but usually -- with exceptions mentioned
below -- for the pure system the s-d-hybridization has been neglected in the
 ideal leads (but not in the spacer). The reason for this neglection  is
mainly technical, since in case of neglected hybridization we can
directly express the matrix elements of the Green's operators $G^L({\vec k}, i)$
 and $G^R({\vec k},i)$ at the boundary layers from the ideal lead to the
ferromagnetic by simple {\it analytical} expressions, whereas the
 numerical calculation
 is much harder if the s-d-hybrization is taken into account already
in the ideal leads, by which the above-mentioned tiny oscillations vanish.
 What is, however, more important:
 In this way we have made explicit that the properties
of the ideal lead really matter, which is not astonishing with
 'ballistic electrons', where the 'contact resistance' from the ideal
leads to the system under study plays an important role.

\noindent{\it
As a first consequence one can learn from this study that it is a mistake
to speak of 'the ideal lead' as a uniquely defined entity:
In ballistic experiments there are {\rm
different} ideal leads, and it is necessary to characterize them, too,
as already stressed in a recent book on electronic transport in mesoscopic
systems}, \cite{Da2}.

It is also interesting at this place to contrast the different, i.e.\
parameter-dependent, results of our Fig.~6 and Fig.~7 with the corresponding
results found by the much more ambitious {\it ab initio} calculations of
\cite{Mathon97} and \cite{Schep}. These authors find values of the CPP-GMR
for Co/Cu-multilayers amounting to $\sim 150$ \%.
 Moreover, for the conductances $\Gamma/M^2[e^2/h]$, instead of our
values 0.8; 0.3; and 0.2 (Fig.\ 6) rsp.\ 0.45; 0.2; and 0.15 (Fig.\ 7) for
the cases (++,$\uparrow$), (++,$\downarrow$), and (+-), respectively, they
find values of
$\Gamma$ corresponding to 0.43$\cdot$10$^{15}$, 0.25$\cdot$10$^{15}$, and
0.165$\cdot$10$^{15}$
$\Omega^{-1}$ $m^{-2}$, i.e.\ $\Gamma/M^2[e^2/h] = 0.67$, 0.39,
and 0.26.

{As a conclusion, our values are correct concerning the order of the three
conductances, and their order of magnitude, but not more. This is exactly what
can be expected.}
\subsubsection{CPP-GMR; bulk impurities}
We would  like to stress at this place
 that with impurities, all calculations have been
performed in the {\it real} lattice space, with cross sections of
$M\times M=$10$\times$10 sites, {\it free} boundary conditions, in the
perpendicular directions,
{{instead
of {\it periodic} ones}}, and with complete
s-d-hybridization everywhere, i.e.\ now always {\it including} the ideal leads.

In  Fig.~8 and Fig.~9, for $n_s=3$ and the parameters (i) and (ii),
respectively, our results for the GMR and the conductances are presented over the
strength of the disorder, i.e.~over the standard deviation
 $\sigma_r$ in equation (\ref{eqBulk_impurities}). For the parameter set
 (i), Fig.\ 8, one obtains first a rapid {\it decrease} of the GMR from GMR $\approx
 1.15$   to ${\rm GMR} \approx 0.45$, when  $\sigma_r$
 increases from 0 to 1, followed
  by  a plateau behaviour with further increase of $\sigma_r$. Concerning
  the  behaviour of $\Gamma_\sigma^{++}$,  $\Gamma_\sigma^{--}$, and
  $\Gamma_\sigma^{+-}$, one can see from the r.h.s.~of Fig.~8 that the
  GMR reflects essentially the behaviour of   $\Gamma_\uparrow^{++}$.

For the parameter set (ii), Fig.\ 9, there is also at first a rapid decrease
 from ${\rm GMR} \approx 0.37$ at $\sigma_r =0$ down to ${\rm GMR}\approx 0.15$
at $\sigma_r \approx 1$, but with a further increase of $\sigma_r$ there is
now a significant {\it re-increase} of ${\rm GMR}$
 back to $\sigma_r \approx 0.22 $ for $\sigma_r
\approx 6$. In fact, in this case  the behaviour of the conductances,
which are plotted on the r.h.s.\ of Fig.\ 9, is significantly
different and more subtle than that of Fig.\ 8.  In particular, the results
of Fig.\ 9 differ from what in {\cite{Asano}} the authors have obtained from a
simple one-band model.

 {\it So from these two figures, Fig.\ 8 and Fig.\ 9, concerning the
dependence on  the strength of disorder $\sigma_r$, one may only conclude that at first a
rapid decrease of ${\rm GMR}(\sigma_r)$ from its value at $\sigma_r=0$ to
 ${\rm GMR}(1)\,\approx 0.4\cdot {\rm
GMR}(0)$
 would be typical; this qualitative behaviour
seems to be essentially
parameter-independent; but for $\sigma_r >1$ a quantitative
calculation is hardly avoidable}, \cite{REMgmr,Moser}.

The reason for the common '{\it decreasing behaviour}'
 for $\sigma_r\widetilde < 1$ is the
following: In this case, an increase of $\sigma_r$ leads to an effective
reduction of the number of  open channels for coherent motion of
an electron from one plane to the next. This reduction concerns mainly
$\Gamma^{++}_\uparrow$ and therefore leads to a {\it decrease} of the GMR, since
the majority-spin d-electrons, {\it without} impurities, would not feel any
potential-variation at all.

 In this connection, a very recent paper of Bruno {\it
et
al.}, \cite{Bruno98}, should be mentioned: There it is shown that
 the conductance of
systems with impurities can be decomposed into two parts, a '{\it
diffusive}' part, arising from scattering processes, where the in-plane $\vec
k_{||}$-vector is not conserved, as is generically always the case for the present
systems with impurities and free-boundary conditions,
 and a '{\it ballistic}' part, to which only those processes
contribute, where the $\vec k_{||}$-vector is the same for
the incoming and outgoing waves. The diffusive part ( rsp.\ ballistic part)
always dominates the conductance, as long as the spacer length is much
shorter (rsp.\ much longer) than the {\it scattering length} $l_s$ of the
system. According to \cite{Bruno98}, $l_s$ should be of the order
of several 100 lattice spacings under similar conditions as in the present
work. So in our systems
with impurities, where $n_s$ is very small ($n_s=3$),
 always the {\it diffusive} part of the conductance
dominates.

\subsubsection{CPP-GMR; interface impurities}
In the following Fig.~10 and Fig.~11 the number $n$ of impurities
is varied between $n=0$  and $n=80$ for every interface layer
(remember that per interface, two layers are involved, with 100 sites per
layer; so we have $2\times n$ impurity sites out of 200 per interface, occupied by a
'wrong atom').

For the parameter set (i), Fig.\ 10, the ${\rm GMR}$ {\it decreases} with
increasing number of impurities, whereas for parameter set (ii),
Fig.~11, there is again the opposite behaviour: Here the ${\rm GMR}$
{\it increases} with increasing number of impurities.

\noindent This difference contains again an important statement in itself, which
should be particularly relevant to the experimentalist\  :

{\it Without a
precise calculation for a specific model one cannot predict whether an
{\rm increase} of the number of interface impurities will lead to an {\rm
increase} or a {\rm decrease} of the ${\rm GMR}$.}

\subsection{CIP-GMR}
\subsubsection{CIP-GMR; no impurities} The formalisms of the preceding
sections can also be applied on equal footings to the situation, where the
current is in the direction of the planes (CIP geometry, see Fig.~2). The
ideal leads have cross-sections of $\Delta z\times \Delta y=N\times M$
with
$N=6+n_s$, see below, while the CIP-conductor consists again of $\Delta z=N$
 monolayers of length $L=\Delta x=10$ and width $\Delta y=M=10$, namely two
ferromagnetic sandwiches of heigth $\Delta z=3$ and inbetween a spacer
sandwich of $\Delta z=n_s$. Again we consider at first the situation without
impurities; however in this geometry we have applied the x-space calculation
right from the beginning, i.e.~with cross sections of size $10\times 10$ and
with {{\it free}} boundaries, although with respect to the coordinate $y$
one could have still performed a Fourier transform. Moreover, as for the CPP
geometry without impurities, this time the s-d-hybridization has been
switched off for the lead wires.

 In Fig.~12 and Fig.~13, the ballistic CIP conductances are plotted as a
 function of $n_s$, for the parameter sets (i) and (ii), respectively.
In the first case, with increasing $n_s$ the conductances $\Gamma/M^2$
 rise, because the spacer cross section increases,
 the limiting value given by the ''shunting'' through the Cu spacer. The GMR shows
irregular oscillations as a function of $n_s$; these oscillations change
with the boundary conditions, the geometry, and the material,
  and the results are more than one order of magnitude smaller than for the
CPP geometry, \cite{REMcip,Gijs,Theeuwen}.
Note that with parameter set (i), the GMR is even  negative on average over
$n_s$, which is not observed experimentally. By discrete Fourier
transformation of GMR  with respect to $n_s$, \cite{Gebele}, we would get
maxima in the Fourier spectrum at periods of 2.5, 4.6 and 8.33 monolayers
(ML) from parameter set (i), Fig.~12, and 2.27 and 5 ML from (ii), Fig.~13.

\subsubsection{CIP-GMR; bulk impurities}
In case of bulk impurities, we have considered a system with $n_s=4$, i.e.\
$N=L=10$. The results are presented as a function of $\sigma_r$, the
 parameter characterizing the 'strength of the disorder' according to
 Eq.~($\ref{eqBulk_impurities}$). Here, in case of the parameter sets (i)
and (ii), we find again quite different behaviour in Fig.\ 14 and 15,
respectively: In case of Fig.\ 14, i.e.\ for parameter set (i), the GMR
{\it increases} at first from 0 to $0.015$, and then for $\sigma_r > 1$ it
roughly  remains constant with increasing $\sigma_r$, whereas in case of
the parameter set (ii), i.e.\ in Fig.\ 15, there is at first a drastic
{\it decrease} from GMR $\sim 0.035$ to GMR $\sim 0$, and then again a roughly
constant value of GMR $\sim 0$, {{\cite{REMcip}}}.
 In both cases the {\it conductances} decrease at first rather fast, and
then there is a slow re-increase; but this behaviour is slightly different
for the up- and down-spin channels in the ($++$)-configuration, and for the
($+-$)-configuration.

 {\it So again the behaviour of the GMR as a function of $\sigma_r$ is
 hardly predictable without an extensive and completely realistic
 calculation for a specific model.}
\subsubsection{CIP-GMR; interface impurities}
A somewhat different conclusion refers to the influence of {\it interface
impurities} for the CIP geometry.  Here in Fig.\ 16 and 17 we use the same
geometry as in the two preceding figures. The GMR rises significantly up to
$50$ \% impurities - i.e.\ this time we have the same trend for both models
(i) and (ii) --, and then it falls again, which is natural since a system of
two {\it three-ML} ferromagnets separated by a 4-ML-spacer and $x>50$ \%
interface impurities is equivalent with two {\it four-ML} ferromagnets
separated by a 2-ML spacer with $x'\,=\,(1-x) < 50$ \% interface impurities.
Here it should be noted that the majority-spin carriers are hardly
influenced by our interface disorder, which is clear since we have assumed
that the s-bands and the majority d-bands of our metallic ferromagnets and
of the metallic spacers are identical. But it should also be noted that
 at the above-mentioned maximum  the
CIP-GMR with {\it interface} scattering reaches much higher values (e.g.\ 9\%
and 8\% respectively, for the parameter sets of (i) and (ii)) as with
our nonmagnetic impurities in the {\it bulk}.
\section{Non-metallic spacers, CPP geometry}
Finally, we come to our most startling result, for the non-metallic spacers,
where we have considered metallic ferromagnets as in model (i), but
{\it {{semi-conducting}} spacer layers} as in the first line of Tables 3 and 4. The
geometrical situation corresponds to the CPP geometry of Fig.\ 1, with
metallic ferromagnets (two-times 3 monolayers) and a variable number $n_s$
of semi-conducting
 spacer layers, i.e.\ with a conduction band and a valence band, inbetween. Of course one now expects a strong exponential
decrease of all conductances with increasing $n_s$, but GMR
=($\Gamma^{++}/\Gamma^{+-})\, -\,1$ may be well defined and even {\it
increase} with increasing $n_s$, as long as the  spacer thickness does not
become larger then the 'dephasing scattering length' or 'spin-flip
scattering length'.

The calculations have been performed by the highly accurate $\vec k$-space
method in reliable numerical accuracy, and we have plotted the conductances
and the GMR against the number
$n_s$ of spacer layers, calculated in the ballistic limit of perfect phase
coherence. In Fig.\ 18 and Fig.~19 there is now the expected drastic exponential
decrease of the conductance,
which for our model
is {\it weaker} for $\Gamma^{++}$ than for $\Gamma^{+-}$. (The decay is
weakest for $\Gamma^{++}_\downarrow$, which is reasonable, since the
 $\downarrow$ d-electrons have a much lower energy gap to the lower edge of
the semi-conductor conduction band than the $\uparrow$ d-electrons,
 and with the $\Gamma^{+-}$ case, in contrast to the $\Gamma^{++}$ case, there
 are the additional 'matching' problems at the interface to the second
ferromagnet. Note
that due to the s-d-hybridization, and with the finite value of $t_d$, the
d-electrons do in any case contribute to the tunneling current.)
 So this means that in our
case, and in the recent calculation of MacLaren {\it et al}, \cite{Butler2},
 to be discussed below, and
 in contrast to Julliere's assumptions, \cite{Jul},  the
tunneling matrix through the semi-conducting spacer is naturally strongly
spin-dependent. (The fact that here for $\Gamma^{+-}$ the exponential decay
with the spacer thickness is stronger than for $\Gamma^{++}$, depends of
course on our model:
for other assumptions on the semi-conductor it
might be just the other-way round.)

  As a consequence, one gets a
positive GMR, which {\it strongly increases with} $n_s$, e.g.~up to GMR $\sim 300$
\% for $n_s\widetilde > \,5$ in case of model (iv), Fig.\ 18, i.e.~for the
model with the 'direct energy gap'. For the case of model (iii) with the
'indirect gap', Fig.\ 19, the increase is at first similar, i.e.~for $n_s
\widetilde < 4$, but then the GMR increases even further to {\it collosal}
values, e.g.~GMR $\sim 2500$ \% for $n_s=7$. However because of the strong
exponential decrease of the conductance, and in view of the fact that we
have neglected indirect transitions involving phonons, this region is
probably beyond experimental realizability, and the results of Fig.\ 18 and
Fig.\ 19 should not be taken literally, although the numerical accuracy of
our calculations seems to be sufficient  (In Fig.\ 18 and Fig.\ 19, the Cunningham
accuracy parameter, see \cite{Cun}, was $m=10$).

In fact, however, in the already cited recent paper of MacLaren {\it et al.},
\cite{Butler2}, which contains an ambitious realistic LKKR calculation of
the ballistic CPP conductances and the corresponding GMR for bcc (001)
Fe/ZnSe/Fe trilayers, the conductances of Fig.\ 3 of that paper resemble
very much -- almost quantitatively -- to the results of our simple model
calculation in Fig.19. One only must take into account that in the paper of
MacLaren {\it et al.}, one Zn-Se bcc double-layer is replaced by two monolayers in
our simple {\it sc} two-band model calculation with the parameter set (iii).
So here we have an example where a fully realistic ab-initio calculation may be
sometimes not necessary after all.

 Finally we present at this place additional comments  on the seminal
paper of Julliere, \cite{Jul}, written already in 1975. Julliere describes the
tunneling conductance from the  ferromagnetic metallic {\it  source} to a metallic
ferromagnetic {\it drain} through an isolating Al$_2$O$_3$ spacer as a
product of  the densities of states (DOS) of the ferromagnetic
source and drain, multiplied by an
effective tunneling matrix element. This should be essentially equivalent to
our approach, which involves two imaginary parts of Green's function, which
 represent in principle densities of states. But from the relevant densities
of states of the magnetic metals, the contributions of {\it confined} electrons, which
contribute to the DOS, but not to the conductance, should be excluded; such
confined states exist e.g.\
 with the
Co down-spin electrons, which are strongly reflected at the Cu interfaces;
so instead of Juliere's DOS, an 'effective DOS' should be used, and since the
reduction  factor of the DOS to this 'effective DOS' is spin-dependent, the
remaining 'effective tunneling matrix elements', if one uses the original
DOS, i.e.\ Julliere's formula,
would in any case depend on the spin, contrary to Julliere's assumption,
 if applied to systems as Co/Al$_2$O$_3$/Co, \cite{MathonPriv}.

In principle, however, our approach is even more explicit, and more demanding,
 since it also involves the
transfers to (and  from) the ideal leads to the ferromagnetic metal, and
from the reservoirs to the ideal leads.

In Julliere's  TMR experiments, and similar experiments performed at
present, the high values of the TMR
predicted in the present paper and also in the paper of MacLaren {\it et al.},
\cite{Butler2}, have {\it not} been observed hitherto. However one should stress
that our insulating layers are extremely thin (from 1 to 15 monolayers only)
with ideal interfaces, whereas the experimental thicknesses of the
Al$_2$O$_3$ spacers are larger and with rougher interfaces, so that in that case
spin-channel mixing might play a role, particularly if on the metallic side
of the interface the transverse spin components are almost as probable as the
$z$-components. The band-structure, e.g.\ that of the semi-conductors
 involved, is also crucial, since the
exponential decay with the spacer-thickness should be weaker for $\Gamma^{++}$
than for $\Gamma^{+-}$, i.e.\ the semi-conducting compounds  involved should have
 properties similar to those of our present models, or to the Fe/ZnSe/Fe-system
 studied by Mac\-Laren {\it et al.}, \cite{Butler2}. In any case, one can only
 speculate that perhaps ballistic point contacts as 'nanocontacts' might
ultimately work, with particular
 few-$\AA$-thick
semi-conducting nanospacers, and with
very clean interfaces.

The suggested experiment concerning the TMR
 would therefore be analogous to the very recent set-up of Garcia {\it et al.},
\cite{Garcia}, but with a small semi-conducting tip of the above-mentioned
kind.

  \section{Conclusions}
We have performed a systematic study of  the CPP-GMR and CIP-GMR (i.e.\ with
{\it Current Perpendicular to the Planes} or {\it Current in the Planes},
respectively), for systems consisting of two 'ideal leads' attached to two
ferromagnetic metallic slabs, which are three monolayers thick, separated by
$n_s$ monolayers of non-magnetic (metallic or non-metallic) spacer material.
A simple two-band model is used throughout, where in the metallic case
 one band mimics the s-bands, and the other one   the
d-bands
 which are spin-dependent in the
ferromagnetic metals, whereas for nonmetallic spacers the
two bands describe
 the valence band and the conduction band, respectively. The
s-d-hybridization is taken into account, and is important for our results.
For the metals, and also for the semi-conducting spacer, we use two
{{specific}} parameter sets,
 corresponding in the metallic case very roughly to the situation
of a Co/Cu/Co-trilayer system with Cu leads, whereas in the non-metallic
case the two parameter sets of the spacer distinguish between situations
corresponding to an 'indirect' and a 'direct' energy gap, respectively.

The calculations have been performed by very accurate
 Green's function methods and
applied to systems without impurities, with {{{nonmagnetic}}
bulk impurities in the spacer,
and with {{magnetic and non-}} {{magnetic}} interface impurities
 produced by mutual interdiffusion of the atoms near the
interface. Within our formalism other situations could also have been
treated, e.g.\ more complicated geometries as (i) the 'constrictions'
studied  experimentally by Garcia {\it el al.}, \cite{Garcia}, or (ii)
 quasi-one-dimensional problems, where the potentials are
constant within a layer, but random from one layer to the next. The
last-mentioned problem has already been extensively studied by J.\ Mathon,
\cite{Mathon97}, and in two-dimensional systems by one of the present
authors, \cite{Boehm99}. Also the recent work of Sanvito {\it et al.} should
be mentioned in this broader context.

 In all cases the influence of the strength of
disorder leads to strong effects, sometimes to an enhanced GMR and sometimes
to a reduction. However for different model parameters the outcome can be
quite different, which is not astonishing after all, since three
conductances are involved in the changes, namely $\Gamma^{(++)}_\uparrow$,
 $\Gamma^{(++)}_\downarrow$, and $\Gamma^{(+-)}_\sigma$, and the GMR effect is --
after all -- a difference effect out of these quantities.

So our
 calculation shows among other results that for the
GMR one can hardly avoid extensive calculations for very realistic models of
the specific systems considered.
Only for the CPP case with impurities in the bulk of a metallic spacer, the GMR
seems to be strongly reduced with increasing strength $\sigma_r$ of the
disorder, but only as long as $\sigma_r$ remains $\widetilde < \, 1$ (Fig.\ 8
and Fig.\ 9). In contrast, interface impurities may lead in the same case to
a reduction of the GMR for one parameter set (Fig.\ 10), but to an
enhancement for another set (Fig.\ 11).

 One of our main results is perhaps the {{suggested}} existence of
a drastic {\it increase} of the GMR with increasing spacer thickness in case
of ballistic conductance through a {\it nonmetallic} spacer, see section 4
above,
where we also  suggest an experimental realization
along the lines of the recent experiments of Garcia {\it et al.},
\cite{Garcia}.

\subsection*{Acknowledgements}
We would like to thank Prof.\ Dr. P.\ Bruno,
Prof. Dr. P. Weinberger, and Prof.\ Dr.\ G.\ Bayreuther, {{and the
experimentalists around him}}, for stimulating discussions.
S.K. thanks the Humboldt foundation for financial support and the University
of Regensburg for its hospitality, and U.K.  does the same  w.r.\ to  the
Polish Academy of Sciences and the Institute of Molecular Physics  in
Pozna\'n.


\section{Introduction\label{Intro}}

At the time of the previous (1994) International Congress of
Mathematicians, steady, but slow, progress was being made on the
classification of simply connected closed smooth 4-manifolds. In
particular, the Donaldson invariants had begun to take a particularly nice
form \cite{KM} (also \cite{FSstructure}), their computations were
becoming more routine \cite{rat}, and their behavior under blowing up
(i.e. taking connected sum with $\overline{\mathbf{CP}}{}^{\hspace{.005in} 2}$)
 was well understood \cite{blowup}. Due to the complexity of the Donaldson
invariants, great hope was held out that an even better understanding of
these invariants would close the books on the classification of
simply connected 4-manifolds.

A few short months after the 1994 ICM, the 4-manifold community was
blind-sided by the introduction of the now famous Seiberg-Witten equations
\cite{W}. Most of the results obtained by using  Donaldson theory were
found to
have quicker, and sometimes more general, counterparts using the Seiberg-Witten
technology. The potential applications of the difficult Donaldson
technology became much more transparent using these new equations. As of
July 1998, there is good news as well as bad news. The good news is that
many of the earlier focus problems have been solved. In particular, the
Thom conjecture \cite{KM2} and  its natural generalizations have been
verified \cite{MST,OS}; also the study of symplectic 4-manifolds
has taken a more central role \cite{T1,T2,TSWG,TCPH}.
The bad news is that recent constructions and computations indicate that
the Seiberg-Witten and Donaldson theories are too weak to distinguish
simply connected smooth 4-manifolds \cite{KL4M}.  It is these
latter constructions and computations that we will discuss at this 1998
International Congress of Mathematicians. It is becoming more
apparent that we are seeing only a small constellation of 4-dimensional
manifolds. More seriously, we are lacking a reasonable conjectural
classification of simply connected closed smooth 4-manifolds.

Current technology has given us many more
4-manifolds than had been expected in 1994. The authors hope that during
the 2002 ICM the
construction of large classes of new 4-manifolds will be discussed; in
particular, they hope that a sufficiently large collection of 4-manifolds
will have been
discovered so as to allow for
some general patterns to emerge and, at least, a conjectural classification
to again be on the books.

\section{The knot surgery construction \label{knots}}

Let $X$ be a simply connected oriented smooth closed 4-manifold. Its most basic
invariant
is its intersection form
\[ Q_X:H_2(X;\mathbf{Z})\otimes H_2(X;\mathbf{Z})\to\mathbf{Z}\]
defined by counting signed transverse intersections of embedded oriented
surfaces
representing given homology classes. It is a famous theorem of M. Freedman
\cite{Fr} that $Q_X$ determines the homeomorphism type of $X$, and an equally
renowned theorem of S.K. Donaldson \cite{Donpoly} that $Q_X$ is not
sufficient to
determine the diffeomorphism type of $X$. In this section we shall discuss
geometric operations on a given smooth 4-manifold which preserve the underlying
topological structure and alter its smooth structure. In particular, we shall
consider the following construction: Let $X$ be a simply connected smooth
4-manifold which contains a smoothly embedded torus
$T$ of self-intersection 0. Given a knot $K$ in $S^3$, we replace a
tubular neighborhood of $T$ with $S^1\times (S^3\setminus K)$ to
obtain the {\it knot surgery manifold} $X_K$.

More formally, this procedure is accomplished by performing 0-framed
surgery on $K$ to obtain the 3-manifold $M_K$. The meridian $m$ of $K$ can be
viewed as a circle in $M_K$; so in $S^1\times M_K $ we have the smooth torus
$T_m=S^1\times m$ of self-intersection~0. Since a neighborhood of
$m$ has a canonical framing in $M_K$, a neighborhood of the torus
$T_m$ in $S^1\times M_K$ has a canonical identification with
$T_m\times D^2$. The knot surgery manifold $X_K$ is given by the
fiber sum \[ X_K=X\#_{T=T_m}S^1\times M_K=(X\setminus T\times D^2)
\cup(S^1\times M_K\setminus T_m\times D^2)\] where the two pieces are
glued together so as to preserve the homology class
$[{\rm{pt}}\times \partial D^2]$. This latter condition does not, in
general, completely determine the isotopy type of the gluing, and
$X_K$ is taken to be any manifold constructed in this fashion.

Because $S^1\times (S^3\setminus K)$ has the same homology as a tubular
neighborhood of $T$ in $X$ (and because the gluing preserves
$[{\rm{pt}}\times \partial D^2]$)
the homology and intersection form of
$X_K$ will agree
with that of $X$. If it is also assumed that $X\setminus T$ is simply
connected,
then $\pi_1(X_K)=1$; so $X_K$ will be homeomorphic to $X$.

In order to distinguish the diffeomorphism types of the $X_K$, we rely on
Seiberg-Witten invariants. We view the Seiberg-Witten invariant of a smooth
4-manifold as a multivariable (Laurent) polynomial. To do this, recall
that the Seiberg-Witten
invariant  of a smooth closed oriented $4$-manifold
$X$ with $b_2 ^+(X)>1$ is an integer-valued function which is defined on
the set of $spin ^{\, c}$ structures over $X$ (cf. \cite{W}). In case
$H_1(X,\mathbf{Z})$
has no 2-torsion (for example, as here where $X$ is simply connected) there
is a
natural identification of the $spin ^{\, c}$ structures of
$X$ with the characteristic elements of $H_2(X,\mathbf{Z})$ (i.e. those elements
$k$ whose
Poincar\'e duals $\hat{k}$ reduce mod~2 to $w_2(X)$).
In this case we view the
Seiberg-Witten invariant as
\[ \text{SW}_X: \lbrace k\in H_2(X,\mathbf{Z})|\hat{k}\equiv w_2(TX)\pmod2)\rbrace
\rightarrow \mathbf{Z}. \]
The sign of $\mathcal{SW}_X$
depends on an orientation of
$H^0(X,\mathbf{R})\otimes\det H_+^2(X,\mathbf{R})\otimes \det H^1(X,\mathbf{R})$. If $\text{SW}_X(\beta)\neq
0$, then $\beta$ is called a {\it{basic class}} of $X$. It is a fundamental
fact that the set of
basic classes is
finite. Furthermore, if $\beta$ is a basic class, then so is $-\beta$ with
$\text{SW}_X(-\beta)=(-1)^{(\text{e}+\text{sign})(X)/4}\,\text{SW}_X(\beta)$ where
$\text{e}(X)$ is
the Euler number and $\text{sign}(X)$ is the signature of $X$.

Now let
$\{\pm \beta_1,\dots,\pm \beta_n\}$ be the set of nonzero basic classes for $X$.
Consider variables $t_{\beta}=\exp(\beta)$ for each $\beta\in H^2(X;\mathbf{Z})$ which
satisfy the
relations $t_{\alpha+\beta}=t_{\alpha}t_{\beta}$. We may then view the Seiberg-Witten
invariant
of $X$ as the Laurent polynomial
\[\mathcal{SW}_X = \text{SW}_X(0)+\sum_{j=1}^n \text{SW}_X(\beta_j)\cdot
(t_{\beta_j}+(-1)^{(\text{e}+\text{sign})(X)/4}\, t_{\beta_j}^{-1}).\]
As an example of this notational device, consider the simply connected
minimally elliptic
surface
$E(n)$ with holomorphic Euler characteristic $n$ and no multiple fibers. Its
Seiberg-Witten invariant is
$\mathcal{SW}_{E(n)}=(t-t^{-1})^{n-2}$ where $t=t_F$ for $F$ the fiber class. Thus,
$\text{SW}_{E(n)}((n-2m)F)=(-1)^{m-1}\binom{n-2}{m-1}$ for
$m=1,\dots,n-1$ and $\text{SW}_{E(n)}(\alpha)=0$ for any other $\alpha$.
When $b^+(X)>1$, the Laurent polynomial $\mathcal{SW}_X$ is a diffeomorphism
invariant of $X$.

For our theorem, we need to place a mild hypothesis on the embedded
torus $T$. We say that a smoothly embedded torus representing a
nontrivial homology class $[T]$ is {\it{c-embedded}} if there is a
neighborhood $N$ of $T$ in $X$ and a diffeomorphism $\varphi:N\to U$ where
$U$ is a neighborhood of a cusp fiber in an elliptic surface and
$\varphi(T)$ is a
smooth elliptic fiber in $U$. Equivalently, $T$ is c-embedded if it
contains two
simple closed curves which generate $\pi_1(T)$ and which bound  vanishing
cycles
in $X$. Note that a c-embedded torus has self-intersection~0.

\begin{thm}[\cite{KL4M}]\label{knotthm} Let $X$ be a simply connected
oriented smooth
4-manifold with $b^+>1$. Suppose that $X$ contains a c-embedded
torus $T$ with $\pi_1(X\setminus T)=1$, and let $K$ be any knot in $S^3$.
Then the knot
surgery manifold $X_K$ is homeomorphic to $X$ and has Seiberg-Witten invariant
\[ \mathcal{SW}_{X_K}=\mathcal{SW}_{X}\cdot\Delta_K(t) \]
where $\Delta_K(t)$ is the symmetrized Alexander polynomial of $K$ and
$t=\exp(2[T])$.
\end{thm}

For example, the theorem applies to the K3-surface $E(2)$ where $T$ is a
smooth elliptic
fiber, and since $\mathcal{SW}_{E(2)}=1$, we have $\mathcal{SW}_{E(2)_K}=\Delta_K(t)$. It is a
theorem of
Seifert that any Laurent polynomial of the
form $P(t)=a_0+\sum\limits_{j=1}^na_j(t^j+t^{-j})$ with coefficient sum
$P(1)=\pm1$ is the Alexander polynomial of some
knot in $S^3$.
Call such a Laurent polynomial an {\it {$A$-polynomial}}. It follows that
if $(X,T)$
satisfies the hypothesis of Theorem~\ref{knotthm}, then for any
$A$-polynomial $P(t)$, there is a smooth simply connected 4-manifold
$X_{P}$ which is
homeomorphic to $X$
and has Seiberg-Witten invariant $\mathcal{SW}_{X_{P}}=\mathcal{SW}_X\cdot P(t)$
where $t=\exp(2[T])$. In particular, for each $A$-polynomial $P(t)$, there
is a manifold
homeomorphic to the K3-surface with $\mathcal{SW}=P(t)$.

The relationship between Seiberg-Witten type invariants and the Alexander
polynomial was
first discovered by Meng and Taubes. In \cite{MT} they showed that the
3-manifold
Seiberg-Witten invariant is related to Milnor torsion.

If one starts with a fibered knot $K$, then $S^1\times M_K$ is a surface bundle
over a torus and thus carries a symplectic structure \cite{Th} for which
$T_m$ is a
symplectic submanifold. Thus if
$X$ is a symplectic 4-manifold
containing a c-embedded symplectic torus $T$, then
$X_K=X\#_{T=T_m}S^1\times M_K$ is also
symplectic \cite{Gompf, MW}. In a fashion similar
to the treatment of the Seiberg-Witten invariant as a Laurent polynomial,
one can
view the Gromov invariant of a symplectic 4-manifold $X$ as a polynomial
$\mathcal{G}r_X = \sum\text{Gr}_X(\beta)\,t_{\beta}$
where $\text{Gr}_X(\beta)$ is the usual Gromov invariant of $\beta$. Let
$\text{A}_K(t)=t^{d}\Delta_K(t)$
denote the normalized Alexander polynomial, where $d$ is the degree of
$\Delta_K(t)$.
As a corollary to Theorem~\ref{knotthm} and the theorems of Taubes relating the
Seiberg-Witten and Gromov invariants of a symplectic
$4-$manifold \cite{TSWG,TCPH} we have:

\begin{cor}[\cite{KL4M}] \label{symcor} Let $X$ be a symplectic
$4$-manifold with
$b^+>1$ containing a symplectic c-embedded torus $T$. If
$K$ is a fibered knot, then $X_K$ is a symplectic $4$-manifold whose Gromov
invariant is \
$\mathcal{G}r_{X_K}=\mathcal{G}r_X\cdot\text{A}_K(\tau)$ \
where $\tau=\exp([T])$.
\end{cor}

This last calculation can also be made purely within the realm of
symplectic topology
\cite{IP,L}. Our interest is directed more to the opposite situation. The
Alexander
polynomial of a fibered knot is monic; i.e. its top coefficient is $\pm1$.
On the other
hand:

\begin{cor}[\cite{KL4M}] \label{nonsymp} If $\Delta_K(t)$ is not monic, then
$X_K$ does
not admit a symplectic structure. Furthermore, if
$X$  contains a homologically nontrivial surface $\Sigma_g$ of genus $g$
disjoint from
$T$  with $[\Sigma_g]^2 < 2-2g$ if $g>0$ or $[\Sigma_g]^2 < 0$ if
$g=0$, then $X_K$ with the opposite orientation does not admit a symplectic
structure.
\end{cor}

Until the summer of 1996, it was still a plausible conjecture (sometimes
called the
`minimal conjecture') that each irreducible simply connected 4-manifold
should admit a
symplectic structure with one of its orientations. The first
counterexamples to this
conjecture were constructed by Z. Szabo \cite{S}.
The knot surgery construction gives a multitude of examples of simply
connected irreducible
`nonsymplectic' 4-manifolds. In fact, if $X$ is simply connected with
$\mathcal{SW}_X\ne0$ and if
$X$ contains a c-embedded torus $T$ with $\pi_1(X\setminus T)=1$, then
Theorem~\ref{knotthm} and  Corollary~\ref{nonsymp} imply that there are
infinitely many
distinct nonsymplectic smooth 4-manifolds $X_K$ homeomorphic to $X$.

If $K_1$ and $K_2$ have the same Alexander polynomial, Seiberg-Witten
invariants are
not able to distinguish $X_{K_1}$ from $X_{K_2}$. For example, take
$X=E(2)$. Then $X_K$ has
a self-intersection 0 homology class $\sigma$ satisfying $\sigma \cdot [T]=1$ which is
represented by an embedded surface of genus $g(K)+1$ where $g(K)$ is the
genus of $K$. One
might hope that these classes could be used to distinguish $X_{K_1}$ from
$X_{K_2}$ when
$g(K_1)\ne g(K_2)$.

\begin{conj} For $X=E(2)$, the manifolds  $X_{K_1}$ and $X_{K_2}$ are
diffeomorphic if
and only if $K_1$ and $K_2$ are equivalent knots.
\end{conj}

The proof of Theorem~\ref{knotthm} proceeds by successively simplifying the
manifold $X_K$
in a fashion which mimics the calculation of the Alexander polynomial of $K$
via skein
relations. Recall that $\Delta_K(t)$ can be calculated via the relation
\begin{equation} \Delta_{K_+}(t)=\Delta_{K_-}(t)+(t^{1/2}-t^{-1/2})\cdot\Delta_{K_0}(t)
\label{boogie}\end{equation}
where $K_+$ is an oriented knot or link, $K_-$ is the result of
changing a single oriented positive (right-handed) crossing in $K_+$ to a
negative (left-handed) crossing, and $K_0$ is the result of resolving the
crossing
as shown in Figure~1.

\centerline{\unitlength .65cm
\begin{picture}(9,4)
\put (1,1){\vector(1,2){1.12}}
\put (2,1){\line(-1,2){0.45}}
\put (1.45,2.1){\vector(-1,2){.575}}
\put (1.25,0.5){$K_+$}
\put (4,1){\line(1,2){0.45}}
\put (4.55,2.1){\vector(1,2){.575}}
\put (5,1){\vector(-1,2){1.12}}
\put (4.25,0.5){$K_-$}
\put (7,1){\line(1,4){.255}}
\put (8,1){\line(-1,4){.255}}
\put (7.255,2){\vector(-1,4){.3}}
\put (7.745,2){\vector(1,4){.3}}
\put (7.25,0.5){$K_0$}
\put (3.3,-.15){Figure 1}
\end{picture}}
\medskip

The point of using \eqref{boogie} to calculate $\Delta_K$ is that $K$ can be
simplified to an unknot via a sequence of crossing changes. One builds a
`resolution tree' starting from $K$ and at each stage adding the
bifurcation of Figure~2,
where each $K_+$, $K_-$, $K_0$ is a knot or 2-component link, and
so that at the bottom of the tree, there are only unknots, and split links.
Then,
because the Alexander polynomial of an unknot is $1$, and is $0$ for a
split link
(of more than one component) one can can work backwards using
\eqref{boogie} to calculate
$\Delta_K(t)$.

\centerline{\unitlength .5cm
\begin{picture}(5.5,5.5)
\put (2.25,4){\line(-3,-4){1.5}}
\put (2.75,4){\line(3,-4){1.5}}
\put (2.25,4.3){$K_+$}
\put (0.25,1.2){$K_-$}
\put (4,1.2){$K_0$}
\put (1.4,.3){Figure 2}
\end{picture}}

The manifold $X_{K_+}$ can be obtained from $X_{K_-}$ by means of a
$(+1)$-log transform on
a nullhomologous torus in $X_{K_-}$, and then the gluing theorems of
\cite{MMS} show that
$\mathcal{SW}_{X_{K_+}}$ can be computed in terms of the Seiberg-Witten invariants
of $X_{K_-}$  and
a manifold $X_{K_-,0}$ obtained by a  0-log transform on $X_{K_-}$. With
some work, this
leads to a related resolution diagram of 4-manifolds where each knot $K'$
corresponds
to $X_{K'}$, and this diagram can be used to prove Theorem~\ref{knotthm}.

We conclude this section by pointing out that the knot surgery construction
can be
generalized to manifolds with $b^+=1$ and to links in $S^3$ of more than one
component in a
more-or-less obvious way. One glues the complements of c-embedded tori in
4-manifolds to
the product of $S^1$ with the link complement. See \cite{KL4M} for details.
For example, if
to each boundary component of $S^1\times (S^3\setminus N(L))$ we glue $E(1)$
minus the
neighborhood of a smooth elliptic fiber, we obtain a manifold with
$\mathcal{SW}=\Delta_L(t_1,\dots,t_n)$, the multivariable Alexander polynomial of the
link. Szabo's
examples in \cite{S} can be obtained from this construction.

\section{Embeddings of surfaces in 4-manifolds \label{surf}}

Knot surgery can also be used to change the embedding of a surface in a
fixed 4-manifold.
To motivate the construction, note that one can tie a knot in the core $\{
0\}\times I$ of a
cylinder $D^2\times I$ by removing a tubular neighborhood of a meridian circle
and replacing it
with a knot complement $S^3\setminus K$. We shall perform a parametrized
version of this
construction in the 4-manifold setting.  Consider an oriented surface $\Sigma$ of
genus $g>0$ which is smoothly embedded in a simply connected 4-manifold
$X$. Let $\alpha$ be a
simple closed curve on $\Sigma$ which is part of a symplectic basis, and let
$\alpha\times I$ be an
annular neighborhood of $\alpha$ in $\Sigma$. In $X$ we see the neighborhood
$D^2\times\alpha\times I$. For a
fixed knot $K$ in $S^3$, we parametrize the above construction so as to
perform it
on each of the cylinders
$D^2\times\{ y\}\times I$, $y\in\alpha$, to obtain an embedded surface $\Sigma_K$. This
is equivalent to
performing knot surgery on the (nullhomologous) rim torus $R=\partial D^2\times\alpha$.
We call this
operation {\it{rim surgery}}.

\begin{thm}[\cite{surfaces}]\label{homeo} Let $X$ be a simply connected
smooth 4-manifold
with an embedded surface $\Sigma$ of positive genus. Suppose that
$\pi_1(X\setminus\Sigma)=1$. Then for each knot
$K$ in
$S^3$, rim surgery produces a surface $\Sigma_K$, and there is a
homeomorphism $(X, \Sigma)\cong (X,\Sigma_K)$.
\end{thm}

The Seiberg-Witten invariant can be used to study these embeddings, but
first, an
auxilliary construction is needed. For each positive integer $g$, let $Y_g$ be
the union of the Milnor fiber of the $(2,2g+1,4g+1)$ Brieskorn singularity
and a
generalized nucleus consisting of the 4-manifold obtained as the trace of
the 0-framed
surgery on $(2,2g+1)$ torus knot in $\partial B^4$ and a $-1$ surgery on a
meridian. Then $Y_g$
is a Kahler surface and admits a holomorphic fibration over
${\mathbf{CP}}^{1}$ with
generic fiber a surface $S_g$ of genus $g$.

Let $(X,\Sigma)$ be as in Theorem~\ref{homeo}, and suppose that the
self-intersection $\Sigma^2=0$.  We call $(X,\Sigma)$ an {\it SW-pair}\/ if
satisfies the property that $\mathcal{SW}_{X\#_{\Sigma=S_g} Y_g}\ne 0$. (In general, if
$\Sigma^2=n>0$, one makes this definition by first blowing up $n$ times.)
For example, if $X$ is symplectic and
$\Sigma$ is a
symplectic submanifold (of square 0), then $X\#_{\Sigma=S_g} Y_g$ is
symplectic, and it
follows that
$(X,\Sigma)$ is an SW-pair. In $X\#_{\Sigma=S_g} Y_g$, the rim torus $R$
becomes homologically
essential and is c-embedded. We can use Theorem~\ref{knotthm} to calculate
Seiberg-Witten
invariants:
\[\mathcal{SW}_{X\#_{\Sigma_K=S_g}Y_g}=\mathcal{SW}_{(X\#_{\Sigma=S_g}Y_g)_K}=
\mathcal{SW}_{X\#_{\Sigma=S_g}Y_g}\cdot\Delta_K(r)\]
where $r=\exp(2[R])$, viewing $[R]$ as a class in the fiber sum. We have:
\begin{thm}[\cite{surfaces}]\label{ss} Consider any SW-pair $(X,\Sigma)$ with
$\Sigma^2\ge0$. If
$K_1$ and $K_2$ are two knots in $S^3$ and if there is a diffeomorphism of
pairs
$(X,\Sigma_{K_1})\cong (X,\Sigma_{K_2})$, then  $\Delta_{K_1}(t)=\Delta_{K_2}(t)$.
\end{thm}

As a special case:
\begin{thm}[\cite{surfaces}] Let $X$ be a simply connected symplectic
4-manifold and $\Sigma$ a symplectically embedded surface of positive genus
and nonnegative self-intersection.
Assume also that $\pi_1(X\setminus\Sigma)=1$. If
$K_1$ and $K_2$ are knots in $S^3$ and if
$(X, \Sigma_{K_1})\cong (X,\Sigma_{K_2})$, then  $\Delta_{K_1}(t)=\Delta_{K_2}(t)$.
Furthermore, if $\Delta_K(t)\ne 1$, then
$\Sigma_K$ is not smoothly ambient isotopic to a symplectic submanifold of $X$.
\end{thm}

The second part of the theorem holds because if $\Sigma_K$ were symplectic,
$X\#_{\Sigma_K=S_g}Y_g$ would be a symplectic manifold. The symplectic form
$\omega$ on this
manifold is inherited from the forms on $X$ and $Y_g$; so $\langle\omega,R\rangle=0$.
But $\mathcal{SW}_{X\#_{\Sigma_K=S_g}Y_g}= \mathcal{SW}_{X\#_{\Sigma=S_g}Y_g}\cdot\Delta_K(r)$, and
it follows that the among the basic classes $k$ of $X\#_{\Sigma_K=S_g}Y_g$,
more than one has $\langle\omega,k\rangle$ maximal. This contradicts the fact that, for a
symplectic manifold, the maximality of $\langle\omega,K\rangle$
characterizes the
canonical class among all basic classes \cite{T2}.

\section{Fiber sums of holomorphic Lefschetz fibrations \label{Lefschetz}}
In this section we shall construct for every integer $g \ge 3$ a pair
$(X_g, X'_g)$ of simply connected complex surfaces carrying holomorphic
genus $g$
Lefschetz fibrations with the property that their fiber sum (along a
regular fiber)
is a symplectic 4-manifold $Z_g$  which supports no complex structure;
in fact $Z_g$ is not even homeomorphic to a complex manifold.

Let $T(p,q)$ denote the $(p,q)$ torus knot in $S^3$ and let $N(p,q)$
denote the 4-manifold obtained by attaching a 2-handle to the 4-ball along
$T(p,q)$
with $0$-framing. It is well known that $N(p,q)$ is a Lefschetz fibration over
$D^2$ with generic fiber a Riemann surface of genus
$g(p,q)=(p-1)(q-1)/2$. Let
$W(p,q)$ denote the canonical resolution of the Brieskorn singularity
$\Sigma(p,q,pq)$, the Seifert-fibered 3-manifold with three exceptional
fibers of order $p$, $q$, and $pq$, and with $H_1=\mathbf{Z}$. It is known that
$W(p,q)$ also supports the structure of a genus $g(p,q)$ Lefschetz
fibration
over $D^2$ with a singular fiber over $0$ which is a sequence of $2$-spheres
plumbed    according to the resolution diagram of $\Sigma(p,q,pq)$. Finally, let
\[Z(p,q)= W(p,q) \cup N(p,q).\]
The manifold $Z(p,q)$ is a rational
surface  which is diffeomorphic to the connected sum of ${\mathbf{CP}}^{\hspace{.005in} 2}$ and $r(p,q)$
copies
of $\overline{\mathbf{CP}}{}^{\hspace{.005in} 2}$ for some computable integer $r(p,q)$. Furthermore, $Z(p,q)$
supports the
structure of a holomorphic Lefschetz fibration whose fiber has genus $g(p,q)$.

Now consider nontrivial torus knots $T(p,q)$ and $T(p',q')$ with the
property that
$g(p,q)=g(p',q')$. (This is  possible for every $g(p,q) \ge 3$.) Let
$F(p,q;p',q')$ denote the fiber sum along a regular fiber of
$Z(p,q)$ with $Z(p',q')$. Then $F(p,q;p',q')$ is a simply connected symplectic
$4$-manifold with
\[c_1^2=10+8g(p,q)-r(p,q)-r(p',q'), \;\;\chi=(b^++1)/2=1+g(p,q).\]

Furthermore, $F(p,q;p',q')$ supports the structure of a Lefschetz fibration
with fiber
of genus
$g(p,q)$. A computation of the Seiberg-Witten invariants of
$F(p,q;p',q')$  shows that, up to sign, there is a unique Seiberg-Witten basic
class. It follows that $F(p,q;p',q')$ is minimal.

\begin{conj} $F(p,q;p',q')$ supports the structure of a complex
4-manifold if and only if $\{p,q\}=\{p',q'\}$.
\end{conj}

\noindent As evidence, consider the pairs $(2, 2n+1)$ and $(3,n+1)$, $n
\not\equiv 2 \mod 3$. For
$F(2,2n+1;3,n)$ one can show that $r(2,2n+1)=4n+4$ and  $r(3,n+1)=3n+7$
so that
\[c_1^2=n-2, \;\;\chi=n+1. \] Thus, $c_1^2=\chi-3$, which violates the
Noether inequality $c_1^2 \ge 2\chi-6$. This means that
$F(2,2n+1;3,n)$ is a minimal symplectic 4-manifold that is not even
homotopy equivalent to a complex manifold.  In fact, it can be shown that
the fiber sum of $Z(2,2n+1)$ with itself is the elliptic surface
$E(n+1)$ and that the fiber sum of $Z(3,n+1)$ with itself is a Horikawa
surface with
$\chi = n+1$. Furthermore $F(2,2n+1;3,n)$ can be obtained from
$E(n+1)$ by removing from $Z(2,2n+1)\setminus F\subset E(n+1)$, $F$ a
regular fiber, the regular neighborhood of the configuration of  ($n-2$)
2-spheres:

\centerline{\unitlength 1cm
\begin{picture}(5,1.15)
\put(.9,.2){$\bullet$}
\put(1,.3){\line(1,0){1.3}}
\put(2.2,.2){$\bullet$}
\put(2.3,.3){\line(1,0){.75}}
\put(3.3,.3){.}
\put(3.5,.3){.}
\put(3.7,.3){.}
\put(4,.38){\line(1,0){.75}}
\put(4.65,.3){$\bullet$}
\put(.35,.6){$-(n+1)$}
\put(2.1,.6){$-2$}
\put(4.55,.6){$-2$}
\end{picture}}
\noindent whose boundary is the lens space $L((n-1)^2,-n)$ and replacing
it with the rational ball that this lens space  bounds. (See \cite{rat}
for all the details concerning this rational blowdown procedure.) Thus
$F(2,2n+1;3,n)$ is the manifold $Y(n)$ constructed in Lemma 7.5 of
\cite{rat}.

\section{Homeomorphic but non-diffeomorphic $4$-manifolds with the same
Seiberg-Witten invariants
\label{sameSW}}

In this section we construct examples of a pair $(X_1,X_2)$ of symplectic
4-manifolds with $X_1$ homeomorphic to $X_2$, $\mathcal{SW}_{X_1}=\mathcal{SW}_{X_2}$, but
$X_1$ is not
diffeomorphic to $X_2$. To do this choose a pair of fibered 2-bridge knots
$K(\alpha,\beta_1)$ and
$K(\alpha,\beta_2)$ with the same Alexander polynomials; for example
$K_1=K(105,64)$ and
$K_2=K(105,76)$ with Alexander polynomial
\[\Delta_K(t)=t^{-4}-5t^{-3}+13t^{-2}-21t^{-1}+25-21t+13t^2-5t^3+t^4.\]
Although these knots have the same Alexander polynomial, they can be
distinguished
by the fact that their branch covers are the lens spaces $L(\alpha,\beta_1)$
and
$L(\alpha,\beta_2)$ which are distinct; in our specific case $L(105,64)$ is
not diffeomorphic to $L(105,76)$. These knots are also
distinguished by their dihedral linking numbers; let
$S_{K_1}$ and
$S_{K_2}$ denote the 2-fold covers of
$S^3$ branched over
$K_1$ and $K_2$, with lifted branched loci $\tilde{K_1}$ and
$\tilde{K_2}$, respectively. Thus we have knots
$\tilde{K_i}$ in $S_{K_i}=L(\alpha,\beta_i)$. Take the $\alpha$-fold covers
of these
lens spaces to obtain links
$L_i=\{K_1^{(i)},\dots,K_\alpha^{(i)}\}$ which are
the lifts of the branch loci $\tilde{K_i}$. The linking numbers
of the links $L_1$ and $L_2$ are known as the `dihedral linking numbers' of
the 2-bridge knot $K(\alpha,\beta)$.

Now perform the knot surgery construction of \S\ref{knots} on the $K3$
surface, replacing $T^2\times D^2$ with
$S^1\times(S_{K_j}\setminus\tilde{K_j})$. The resulting 4-manifolds are
the manifolds $X_i$. Either by adapting the arguments of
\cite{KL4M} or by using \cite{IP} or \cite{L}, it can be shown that
$\mathcal{SW}_{X}=\mathcal{SW}_{Y}=\Delta_K(t)\cdot \Delta_K(-t)$.
Unfortunately, the $X_i$ are not simply connected (but are
homeomorphic). In particular,
$\pi_1(X_1)=\pi_1(X_2)=\mathbf{Z}_{\alpha}$, and the
$\alpha$-fold covers $\tilde{X_1}$ and $\tilde{X_2}$ of $X_1$ and $X_2$ are not
diffeomorphic. To see this, observe that $\tilde{X_i}$ is
obtained as our link construction in \cite{KL4M} (cf. \S~\ref{knots}) by
gluing one
copy of
$E(2)$ minus a neighborhood of a smooth elliptic fiber to every boundary
component of
$S^1\times (S^3\setminus L_i)$. It follows from \cite{KL4M} that
\[ \text{SW}_{\tilde{X_i}}=\Delta_{L_i}(t_1,\dots,t_\alpha)\cdot
    \prod_{j=1}^{\alpha} (t_j^{1/2}-t_j^{-1/2})\]
\noindent Since the linking numbers of the links $L_1$ and $L_2$ are
different, it can be shown the Hosokawa polynomials of the links
$L_1$ and $L_2$, when evaluated at $1$ are distinct \cite{FSHos}. Thus their
Alexander polynomials are different and $\tilde{X_1}$ is not diffeomorphic to
$\tilde{X_2}$.

There is a lesson to be learned from these examples. One must consider the
Seiberg-Witten invariants of a 4-manifold $X$ together with those of all of
its covers as the appropriate invariant for $X$.

\section{Nonsymplectic 4-manifolds with one basic class}

Recall from \S~\ref{knots}, that if $k$ is a basic class of $X$, so is $-k$.
Because of this, we say that $X$ {\it{has $n$ basic classes}} if the set $\{k\,
|\,\text{SW}_X(k)\ne0\}/ \{\pm1\}$ consists of $n$ elements.
There are abundant examples of $4$-manifolds with one basic class. Minimal
nonsingular
algebraic surfaces of general type have one basic class (the canonical class)
\cite{W}. The authors and others have constructed many examples of minimal
symplectic
manifolds with one basic class and $\chi -3\le {c_1}^2 < 2\chi -6$. (These
manifolds
cannot admit complex structures due to the geography of complex surfaces.)
However, the
examples described here are the first nonsymplectic manifolds with one basic
class.

Let $X=E(2)$ and $T$ a smooth elliptic fiber. For a knot $K$ of genus $g$
form the knot
surgery construction to obtain $X_K$. In $X_K$ there is a surface
$\Sigma$ of genus $g+1$ with $[\Sigma]^2=0$ and $[\Sigma]\cdot [T]=1$. Let $M$ be the
3-manifold obtained from 0-surgery on the trefoil knot. Then $S^1\times M$ is a
$T^2$-fiber
bundle over $T^2$. The fiber sum of $g+1$ copies of the fiber bundle gives
a 4-manifold
$Y$ which is an $F=T^2$-bundle over a surface of genus $g+1$, and it is
easily seen that
there is a section $C$. Furthermore, $Y$ is a symplectic 4-manifold
with $c_1(Y)=-2g[F]$. Our example, corresponding to the genus $g$ knot $K$ is
$ Z_K=X_K\#_{\Sigma=C}Y$.
We perform this fiber sum so that $Z_K$ is a spin 4-manifold \cite{Gompf}.
It can be
seen to be simply connected.

Write the symmetrized Alexander polynomial of $K$ as
$\Delta_K(t)=a_0+\sum\limits_{n=1}^da_n(t^n+t^{-n})$, and call $d$ the
{\it{degree}} of
$\Delta_K(t)$. Since the genus of $K$ is $g$, we have
$d\le g$. If $K$ is an alternating knot, for example, then $d=g$. Say that the
Alexander polynomial of $K$ has {\it{maximal degree}} if $d=g$. Using
techniques of
\cite{MST} we calculate:

\begin{thm}[\cite{1basic}] Let $K$ be a knot in $S^3$ whose Alexander
polynomial has
maximal degree. Then $Z_K$ has one basic class, $k$, with
$|\text{SW}_{Z_K}(k)|=a_d$, the top coefficient of
$\Delta_K(t)$. When $|a_d|>1$, $Z_K$ is nonsymplectic.
\end{thm}