Skip to Main Content
Book cover for CP Violation CP Violation

Contents

Book cover for CP Violation CP Violation

There are many extensions of the standard model in which fermions with non-standard SU(2)U(1) quantum numbers—‘exotic fermions’—naturally occur.62 For example, in the grand unified theory based on the Lie algebra E6, each family of left-handed fermions is in the 27 representation of that algebra. If we consider the SU(3)cSU(2)U(1) subgroup of E6, that representation has the following branching rule:

(24.1)

The 27 contains, apart from the fifteen usual chiral fields in the first line of eqn (24.1)—the doublet of quarks, the up antiquark, the down antiquark, the doublet of leptons, and the positron—, a vector-like isosinglet quark of charge 1/3, a vector-like isodoublet of leptons, and two isosinglet neutrinos. Here, we denote by ‘vector-like fermions’ those fermions whose left-handed and right-handed components transform in the same way under SU(3)cSU(2)U(1).

In this chapter we concentrate our attention on vector-like isosinglet quarks. Our interest in extensions of the SM with those exotic fermions is justified for the following reasons:

1.

They provide a framework for having a naturally small violation of the 3×3 unitarity of the CKM matrix. This leads to non-vanishing, but naturally suppressed, flavour-changing neutral currents (FCNC). The presence of FCNC opens up many interesting possibilities for rare K and B decays, as well as for CP asymmetries in neutral-B decays.

2.

Adding isosinglet quarks to the SM leads to new sources of CP violation. In particular, isosinglet quarks enable one to achieve spontaneous CP violation (SCPV) using a very simple set of scalar fields, with only one non-Hermitian SU(2)U(1)-singlet scalar added to the usual doublet—see § 24.7.

3.

Extensions of the SM with isosinglet quarks may solve the strong CP problem through an implementation of the mechanism of Barr (1984) and Nelson (1984)—see § 27.7.1.

Models with vector-like quarks have been considered extensively in the literature. Early references include the papers by Aguila and Cortés (1985), Fishbane et al. (1985, 1986), Barger et al. (1986), Branco and Lavoura (1986), and Enqvist et al. (1986). A revived interest in these models was generated by the study of the impact of heavy isosinglet quarks on the CP-violating asymmetries in neutral-B decays, both in the limit where Z-exchange gives the dominant contribution to Bd0Bd0¯ mixing (Nir and Silverman 1990) and in the general case where both Z-exchange and the standard box-diagram contributions are taken into account (Branco et al. 1993). Recent analyses of the experimental constraints on these models may be found in the papers by Silverman (1996) and by Barenboim and Botella (1998).

Our extension of the SM has ng doublets of left-handed quarks,

(24.2)

together with the following singlet quark fields: np charge 2/3 left-handed PL,nn charge –1/3 left-handed NL,ng+np charge 2/3 right-handed PR, and ng+nn charge –1/3 right-handed NR:

(24.3)

Therefore, the electromagnetic current is

(24.4)

and the electroweak gauge interactions are

(24.5)

The quark mass terms may in general be written

(24.6)

We shall use the letter x to denote either p or n in equations which hold for both x = p and x = n. The mx are ng×(ng+nx) matrices. They contain the usual ΔT=1/2 mass terms (T is the weak isospin), arising from the VEV of one or more Higgs doublets. Their mass scale should be mv. The Mx are nx×(ng+nx) matrices. They contain ΔT=0 mass terms. Since these are SU(2)U(1)-invariant, they may be present in the Lagrangian prior to the spontaneous breaking of the gauge symmetry. Not being protected by that symmetry, their mass scale M may be significantly larger than v.

We may denote p and n the full mass matrices for the up-type and down-type quarks, respectively:

(24.7)

These are (ng+nx)×(ng+nx) matrices.

We denote the quark mass eigenstates by uα and dk, where α and k run over the ranges in eqns (24.3). The unitary matrices Wp and Wn relate the left-handed-quark weak and mass eigenstates:

(24.8)

The Wx are (ng+nx)×(ng+nx) matrices. It is useful to write them as

(24.9)

where the Ax are rectangular matrices consisting of the first ng rows of Wx, while the Bx consist of the last nx rows of Wx. The unitarity of Wx implies, on the one hand,

(24.10)

and, on the other hand,

(24.11)

The Wx must be chosen such as to diagonalize the Hermitian mass matrices xx:

(24.12)

The real diagonal matrices Dp and Dn contain the masses of the up-type and down-type quarks, respectively.

The electromagnetic current in eqn (24.4) may be written in terms of the mass eigenstates as

(24.13)

The gauge interactions in eqn (24.5) are

(24.14)

where the mixing matrices are

(24.15)

The (generalized) CKM matrix V is an (ng+np)×(ng+nn) rectangular matrix. The CKM matrix in this model is not unitary as in the SM. The mixing matrices for the neutral currents Zx are Hermitian (ng+nx)×(ng+nx) matrices. In general, they are not diagonal, i.e., flavour-changing neutral currents are present.

The generalized CKM matrix is not, in general, unitary. Indeed, if npnn it is not even a square matrix. However, from the unitarity of Wp and Wn it follows (Branco et al. 1993) that there is an (ng+np+nn)×(ng+np+nn) matrix

(24.16)

which is unitary. This may easily be checked using eqns (24.10) and (24.11).

The deviations of V from unitarity are closely related to the presence of FCNC. Indeed, because of the first eqn (24.11),

(24.17)

Thus, if V was unitary then Zp and Zn would be the unit matrix, and FCNC would be absent. From the first eqn (24.11) it also follows that

(24.18)

Equations (24.17) may be elegantly derived without ever referring to the diagonalization of the quark mass matrices (Bamert et al. 1996). From eqns (11.12) and (11.3), we know that the gauge interactions in the SU(2)U(1) gauge theory are given by

(24.19)

In the vector space spanned by the u and the dLk, the matrices T+ and T are given by

(24.20)

cf. eqn (24.14). Indeed, eqns (24.20) may be looked upon as the definition of the CKM matrix. From eqn (24.19), the matrix which determines that part of the neutral current which is not proportional to Jemμ is

(24.21)

This result is very general: it holds in the SU(2)U(1) gauge theory for any standard or exotic Q=2/3 and Q=1/3 quarks. It is valid if one adds to the SM not only vector-like isosinglet quarks, but also vector-like isodoublet quarks and mirror quarks (Lavoura and Silva 1993b).

It is known experimentally that the FCNC are very much suppressed, although the detailed bounds depend on the flavours involved—see for example Lavoura and Silva (1993a) and Silverman (1996). Therefore, in order for the class of models considered here to be plausible, they must have a mechanism for natural suppression of the FCNC. In order to show that such a mechanism exists, we have to explicitly diagonalize the quark mass matrices.

In the discussion that follows, it will be convenient to have p and n in a special form, which we can obtain by using the freedom that one has to make weak-basis transformations (WBT). These are transformations of the quark fields which leave gauge in eqn (24.5) invariant. It is easy to convince oneself that, by making a WBT, one can bring both p and n to the form

(24.22)

where the M̂x are nx×nx matrices, diagonal and real by definition. The matrices Gx and Jx are complex and have dimensions ng×ng and ng×nx, respectively. There is still some freedom left to make a WBT such that either Gp or Gn is made diagonal and real; we shall take Gp to be diagonal and real.

Let us write

(24.23)

The matrices Kx and Tx have dimensions ng×ng and nx×nx, respectively. The matrix Rx is ng×nx, while Sx is nx×ng. Equation (24.10) now reads

(24.24)

We shall also write the matrices Dp and Dn as

(24.25)

Equation (24.12) reads

(24.26)

All the above equations are exact, no approximations have been done yet. We now assume that Gx and Jx are m, while M̂xM, with Mm. One obtains the following solution of eqns (24.24) and (24.26) to leading order in m/M:

(24.27)

together with M¯x2=M̂x2; while Kx is the unitary matrix which diagonalizes GxGx:

(24.28)

In particular, since we have chosen to work in the weak basis where Gp is diagonal, Kp=1ng to leading order.

The quark mixing matrix entering the charged current may now be computed and one obtains

(24.29)

It is clear from eqn (24.29) that the mixing of standard quarks with isosinglet quarks is suppressed by m/M. Also, the matrix Kn, which gives the interactions among the usual quarks, is unitary up to terms m2/M2.

The mixing matrices for the weak neutral current are

(24.30)

The ng×ng upper-left submatrix, which governs the neutral-current interaction among the usual quarks, is the unit matrix up to a correction of the form KxJxM̂x2JxKx. Thus, flavour-changing neutral currents do arise, but they are naturally suppressed by a factor m2/M2, because Jxm while M̂xM.

When the SM is extended through the addition of vector-like quarks, the CKM matrix is no longer unitary. Counting the independent CP-violating phases in V becomes complicated, especially when there are both Q=1/3 and Q=2/3 vector-like quarks. We shall return to this question later; at this stage, we only want to count the number of phases in V when that matrix is evaluated in the approximation leading to eqn (24.29). Since in this approximation Kn is unitary, the number of physical phases in it is just the same as in the SM, i.e., (ng1)(ng2)/2. In the upper-right block of V, we have the matrix Jn which, being an ng×nn complex matrix, would in general contain ngnn phases. However, nn phases may be eliminated by rephasing the last nn fields dL. Analogously, Jp, in the lower-left block of V, has (ng1)np phases. Therefore, one has altogether

(24.31)

physical phases in V. We shall later show that this is indeed the correct number of CP-violating phases.

The most general CP transformation of the quark fields leaving eqn (24.5) invariant is63

(24.32)

We have suppressed the space-time variables for the sake of clarity. The matrix UL is ng×ng unitary, and the matrices WRx and WLx are (ng+nx)×(ng+nx) and nx×nx unitary, respectively.

Requiring CP invariance of the mass terms in eqn (24.6) leads to the following conditions:

(24.33)

The existence of unitary matrices UL,WRp,WRn,WLp, and WLn which satisfy eqns (24.33) is necessary and sufficient for CP invariance of gauge+M. The fulfilment of this condition in any particular weak basis is equivalent to its fulfilment in any other weak basis, as one easily checks.

In order to find the number of independent CP restrictions, we shall again use the weak basis in which the mass matrices have the form in eqn (24.22). We remind the reader that M̂p,M̂n, and Gp are diagonal and real by definition. In order to eliminate exceptional cases, of null measure in parameter space, we assume that the diagonal matrix elements of each of these three matrices are all different—the matrices are non-degenerate. Under these conditions, we easily find that eqns (24.33) imply

(24.34)

while the matrices WRp and WRn must take the block form

(24.35)

with XL an arbitrary ng×ng unitary matrix. The CP-invariance conditions of eqns (24.33) get reduced to

(24.36)

As XL may be chosen at will, the first eqn (24.36) is in fact equivalent to the simpler condition (see the argument in § 14.2)

(24.37)

Thus, CP invariance imposes restrictions on the complex phases in the matrices Jp,Jn, and GnGn. Remember that the Jx are complex matrices of dimension ng×nx, while GnGn is an ng×ng Hermitian matrix. The number of independent phases in these three matrices is

(24.38)

CP invariance constrains these phases to be equal to differences of ng phases αo,np phases βr, and nn phases γs. The number of independent CP restrictions therefore is

(24.39)

This coincides with the number of CP-violating phases in the CKM matrix, given by eqn (24.31).

From the first two eqns (24.33) one readily derives the following necessary conditions for CP invariance:

(24.40)

where a and b are integers, f is an odd integer, and hxmxmx. The conditions in eqn (24.40) are entirely analogous to those obtained in the SM, cf. eqn (14.21). In models with vector-like quarks, there are necessary conditions for CP invariance of a different type, for example

(24.41)

where HxmxMxMxmx.

The question of finding a set of necessary and sufficient conditions for CP invariance, expressed in terms of weak-basis invariants, is in general very complicated. For the case of a minimal extension in which only one vector-like isosinglet down-type quark is added to the three-generation SM—i.e., ng=3,nn=1, and np=0—it has been shown (Aguila et al. 1998) that

(24.42)

are necessary and sufficient conditions for CP conservation.

These invariant conditions are especially useful for analysing the limit where some of the quark masses are effectively degenerate. This is the case in very-high-energy collisions, where the natural asymptotic states no longer are hadronic states but rather quark jets. At high energies, it should be very difficult, if not impossible, to identify the flavour of the quark jets. In the extreme chiral limit mu=mc=md=ms=0 there is no CP violation in the SM, because there are identical-charge quarks which are degenerate. However, in the model with one down-type vector-like quark, there is CP violation even in that limit. Then, the strength of CP violation is controlled by the second invariant in eqns (24.42), all other invariants being proportional to that one (Aguila et al. 1998).

With the addition of vector-like quarks to the SM, parametrizing the CKM matrix V becomes rather involved, essentially due to the fact that it no longer is a unitary matrix. Still, in some cases one may use the fact that V is a submatrix of the larger unitary matrix Y in eqn (24.16). Various approaches to this problem have been proposed, including:

1.

parametrization through Euler angles and phases (Branco and Lavoura 1986);

2.

parametrization through the moduli of matrix elements and the arguments of quartets (Branco et al. 1993);

3.

Wolfenstein-type parametrization (Lavoura and Silva 1993a).

Here we shall only describe the first of these approaches, which is interesting because it yields a different way of counting the number of CP-violating phases in V.

This type of parametrization is possible when there are vector-like quarks of only one electric charge, either Q=1/3 or Q=2/3. In this case V consists of the first ng rows of a unitary matrix which, in principle, does not have any zero matrix elements. This is no longer true when both np and nn are non-zero; then, the unitary matrix Y defined in eqn (24.16) has a zero submatrix, and its parametrization through Euler angles and phases is awkward.

For definiteness, let us assume that there are nn down-type vector-like quarks, but no up-type vector-like quarks, i.e., that np=0. In this case, one may choose, without loss of generality, a weak basis in which the up-type-quark mass matrix is diagonal and real non-negative. Then, V consists of the first ng rows of the (ng+nn)×(ng+nn) unitary matrix Wn which diagonalizes nn. The task then consists of finding a parametrization of Wn with only (ng1)(ng/21+nn) phases in its first ng rows, according to eqns (24.31) and (24.39).

In order to construct such a parametrization, let us introduce the orthogonal matrices

(24.43)

and the rephasing matrices

(24.44)

We may write Wn as a product of ‘complex rotations’ Ωij:

(24.45)

where N=ng+nn, and the unitary matrices Ωij mix the ith and jth rows and columns, and contain one rotation angle and three phases:

(24.46)

Most of the rephasings Ij(δk) in the product in eqn (24.45) may be exchanged in position in such a way that they are brought out of the product of rotations Oij, and then they become equivalent to rephasings of the quark fields. If we omit writing down those trivial rephasings, it can be shown (Anselm et al. 1985) that the matrix Wn may be written as

(24.47)

Explicit computation of the N×N unitary matrix in eqn (24.47) yields that it has (ng1)(Nng/21) phases in its upper ng rows. This coincides with eqn (24.31), as nn=Nng.

Let us consider the example of E6 with three families in the 27 representation. In each 27 there is one standard quark doublet together with one vector-like isosinglet down-type quark; therefore, ng=nn=3 and np=0. The quark mixing matrix V consists of the first three rows of the 6×6 unitary matrix Wn, in the weak basis in which the up-type-quark mass matrix is diagonal. We parametrize

(24.48)

The distribution of phases among the rows of this matrix turns out to be

(24.49)

i.e., the first row of V is real, the second row depends on four phases, the third row depends on seven phases, and so on. This coincides with the formula (ng1)(Nng/21) for N=6.

Models with vector-like quarks provide one of the simplest scenarios for spontaneous CP violation. We shall illustrate this feature by considering a minimal model proposed by Bento et al. (1991). This consists of the SM supplemented with only one charge 1/3 vector-like quark and one non-Hermitian scalar singlet S. For later use, we shall also introduce a discrete symmetry under which S,NL, and NR4 change sign. This discrete symmmetry will not, however, play any role in the discussion of SCPV; it will only be important in the discussion of the strong CP problem in § 27.7.1

The scalar potential and the pattern of vacuum expectation values (VEVs) for this model have been discussed in § 23.6. One obtains the VEVs for ϕ and S in eqns (23.38), with a non-trivial phase α. However, we have emphasized in § 23.6 that, in the context of the SM with standard quarks only, the non-trivial vacuum phase does not lead to spontaneous CP breaking. In the extension of the SM with vector-like quarks the situation is different, because the presence of extra interactions means that there is less freedom in the choice of a CP transformation and, in particular, eqns (23.41) are no longer a valid CP transformation. This can readily be seen by examining the most general quark Yukawa couplings and the mass term consistent with the SU(2)U(1) gauge symmetry and with the discrete symmetry:

(24.50)

Here, we have reserved the notation NR only for NR4, the singlet quark which changes sign under the discrete symmetry, while for the other PR and NR fields we have returned to the SM notation and denoted them pR and nR, respectively. We assume the CP transformation in eqn (23.36), together with

(24.51)

Hence, the 1×3 matrices of Yukawa couplings F and F are real, just as the 3×3 matrices Γ and Δ, and the ΔT=0 mass μ. On the other hand, the transformation in eqn (23.41) would be a symmetry of the Lagrangian only if F were equal to F. We thus see that the presence of extra Yukawa interactions ensures that exactly the same vacuum structure now leads to CP violation, while in § 23.6 it did not.

Since the model has no charge 2/3 vector-like quarks, we may choose a weak basis in which the up-type-quark mass matrix is diagonal and real. The down-type-quark mass matrix is

(24.52)

where

(24.53)

(Notice that would be real if F were equal to F.) The matrix in eqn (24.52) is not in the standard form of eqn (24.22), but we can bring it to that form by means of a unitary transformation of the right-handed quark fields. Let us define

(24.54)

We may then construct the 4×4 unitary matrix

(24.55)

where X is a 3×3 matrix and Y is a 1×3 matrix. The unitarity of U implies, in particular,

(24.56)

and

(24.57)

We then effect the weak-basis transformation

(24.58)

where we have used eqns (24.54) and (24.57). We thus obtain a matrix of the form in eqn (24.22). One may now use the results of § 24.3, provided that one assumes that Mv and M2vV64 The CKM matrix is a 3×4 matrix; its block connecting the standard quarks is a 3×3 matrix K which, in a first approximation, is the unitary matrix which diagonalizes

(24.59)

where we have used eqn (24.56). Because of the presence of the complex matrix in eqn (24.59), it is clear that, if V and M are of the same order of magnitude—or, equivalently, if V and μ are of the same order of magnitude—then the matrix K will in general contain a CP-violating phase and, moreover, this phase will not be suppressed by small mass ratios (Bento et al. 1991).

It is worth recalling the main features of the model that we have described. The scalar sector consists of the standard doublet and a non-Hermitian singlet S. In the fermion sector, only one charge 1/3 isosinglet quark N is added. One imposes CP invariance of the Lagrangian. CP is spontaneously broken through the phase of 0S0. The crucial role played by the Yukawa couplings connecting N with the standard quarks should be emphasized. On the one hand, those couplings render the phase of 0S0 genuinely CP-violating, by restricting the allowed CP transformations of the Lagrangian; on the other hand, it is through these couplings that that phase leads to the complexity of the 3×3 block of the CKM matrix connecting the standard quarks.

It is important to stress that this model may also be important because, if μ and V are made higher than the cosmological-inflation scale, then domain walls, which are a problem which typically plagues models of SCPV, will be absent. Spontaneous CP breaking can occur at such a high energy scale that domain walls disappear during inflation, while CP violation in the CKM matrix remains unsuppressed.

As we have emphasized, one of the salient features of models with vector-like quarks is the fact that the 3×3 CKM matrix connecting the standard quarks is no longer unitary and, as a result, there are FCNC coupling to the Z.

We shall first present some experimental constraints on FCNC in the down-type-quark sector. The relevant Lagrangian is given by

(24.60)

The flavour-changing parameters Zkj are closely related to the deviation of the CKM matrix V from unitarity, through the relations

(24.61)

The strongest constraint on Zds stems from the experimental upper bound

(24.62)

As in Appendix D, we compare the process K+π+νν¯ with the charged-current process K+π0e+νe, which has

(24.63)

and using flavour-SU(3) symmetry, one obtains (Lavoura and Silva 1993b)

(24.64)

where the factor 3 corresponds to the sum over the three neutrino flavours. From eqns (24.62)–(24.64) one obtains

(24.65)

There are other limits, on Re(Zds2) and on Im(Zds2), arising from ΔmK and from ϵK, but eqn (24.65) is the stringest bound on Zds.

The best limit on Zdb and on Zsb is derived from the experimental bound (Particle Data Group 1996)

(24.66)

The FCNC contributes to the decay Bμ+μX, where X is an arbitrary set of particles. A straightforward calculation gives (Parada 1996)

(24.67)

where the function f has been given in eqn (15.18). From eqns (24.66) and (24.67) one obtains

(24.68)

The new contribution to B0B0¯ mixing arising from the Z-mediated ΔB=2 tree-level diagram can have a significant impact on CP asymmetries in B0 decays. Let us write the complete matrix element M12 of Bd0Bd0¯ mixing as

(24.69)

where M12SM is the standard-model box-diagram contribution, while (Barenboim et al. 1998)

(24.70)

where xt=mt2/mW2, the Inami–Lim (1981) function S0 has been defined in eqn (B.16), and

(24.71)

is another Inami–Lim function. The first term in eqn (24.70) corresponds to the box diagram; the second term results from the Z-exchange tree-level diagram; the third term arises from a one-loop Z-vertex correction, which is relevant for small values of Zbd (Barenboim and Botella 1998). Analogous expressions hold in the case of Bs0Bs0¯ mixing. From eqn (24.68) it can be readily verified (Branco et al. 1993) that in the case of Bd0Bd0¯ mixing the Z contribution can be the dominant one, while in the case of Bs0Bs0¯ mixing it can at most compete with the usual box diagram.

The study of CP asymmetries should provide a very sensitive probe of FCNC. The measurement of those asymmetries with the expected experimental uncertainties may detect FCNC effects (Branco et al. 1993; Barenboim et al. 1998) even for rather small values of Zbd, at the level

(24.72)

It has been shown (Branco et al. 1998) that, in a model with an extra singlet complex scalar, the first-order electroweak phase transition can be strong enough to avoid the baryon-asymmetry washout by sphalerons. The crucial point is the fact that in this model there is CP violation even in the extreme chiral limit where md=ms=mu=mc=0. As a result, and contrary to what happens in the SM, there is no strong suppression of the baryon asymmetry by the ratio of light-quark masses over the critical temperature Tc. The dominant contribution to the baryon asymmetry was estimated to be (Branco et al. 1998)

(24.73)

where vw is the bubble-wall velocity. The appearance of the weak-basis (WB) invariant in the numerator of eqn (24.73) was to be expected, since it is the invariant which controls the strength of CP violation in the extreme chiral limit. That WB invariant can be expressed in terms of quark masses and mixing angles, and one obtains in leading order

(24.74)

when there is only one vector-like Q=1/3 quark D, with mass mD. If mD=200GeV and Im(VtbV4DVtDV4b)=104 one obtains the right order of magnitude for nB/s.

Models with vector-like quarks have new sources of CP violation and provide a well-defined framework to study deviations from unitarity of the 3×3 CKM matrix, as well as flavour-changing neutral currents mediated by the Z. Those models have a rich phenomenology, which may be tested through the search for rare K and B decays, and through the study of CP asymmetries at B factories.

With the addition of at least one vector-like quark to the SM, spontaneous CP violation can be achieved with a very simple Higgs system, consisting of only one non-Hermitian scalar, singlet under SU(2)U(1). CP violation originates in the phase of the VEV of the singlet scalar. Through the mixing of the vector-like quarks with standard quarks, this phase generates an unsuppressed CP-violating phase in the 3×3 block of the CKM matrix connecting the usual quarks.

We are now in a position to add another item to the summary of models in § 23.7:

ONE HIGGS DOUBLET AND ONE SCALAR SINGLET

(CP SYMMETRY IMPOSED ON THE LAGRANGIAN)

No extra fermions: The vacuum contains one irremovable phase, but that phase does not lead to spontaneous CP breaking. The scalar sector does not lead to any new sources of CP violation.

Isosinglet quarks: The physical phase in the VEV of the singlet non-Hermitian scalar leads to spontaneous CP violation—it generates a complex CKM matrix. This is a minimal realization of the Nelson–Barr mechanism.

Notes
62

Generic properties of models with exotic quarks have been studied, for example, by Langacker and London (1988), Joshipura (1992), Lavoura and Silva (1993a,  b), and Nardi et al. (1995).

63

We have eliminated the inconsequential phase ξW from the CP transformation of W± in eqn (13.9).

64

The assumption is indeed that the vacuum expectation value V and the ΔT=0 mass μ are of the same order of magnitude, which is much larger than v. Then, both X and Y have matrix elements 1.

Close
This Feature Is Available To Subscribers Only

Sign In or Create an Account

Close

This PDF is available to Subscribers Only

View Article Abstract & Purchase Options

For full access to this pdf, sign in to an existing account, or purchase an annual subscription.

Close