We classify the allowed order parameter symmetries in multilayer cupra
tes and discuss their physical consequences, with emphasis on Josephso
n tunneling and impurity scattering. Our solutions to the gap equation
are based on highly nonspecific forms for the inter- and intraplane p
airing interactions in order to arrive at the most general conclusions
. Within this framework, the bilayer (N = 2) case is discussed in deta
il with reference to Y-Ba-Cu-O (YBCO) and Bi-Sr-Ca-Cu-O (BSCCO) and th
e related Landau-Ginzburg free energy functional. Particular attention
is paid to the role of small orthorhombic distortions as would derive
from the chains in YBCO and from superlattice effects in BSCCO, which
give rise to a rich and complex behavior of the multilayer order para
meter. This order parameter has N components associated with each of t
he N bands or layers. Moreover, these components have specific phase r
elationships. In the orthorhombic bilayer case the (s, - s) State is o
f special interest, since for a wide range of phase space, this state
exhibits rr phase shifts in corner Josephson junction experiments. In
addition, its transition temperature is found to be insensitive to non
magnetic interplane disorder, as would be present at the rare earth si
te in YBCO, for example. Of particular interest, also, are the role of
van Hove singularities which are seen to stabilize states with d(x)2(
-y)2-like Symmetry (as well as nodeless s states) and to elongate the
gap functions along the four van Hove points, thereby leading to a sub
stantial region of gaplessness. We find that for these rather general
models of the pairing interaction the d(x)2(-y)2-like states are the m
ost stable solutions in a large region of parameter space. In this way
, they should not be specifically associated with a spin fluctuation d
riven pairing mechanism.