INTRODUCTION TO VECTORS AND TENSORS
Volume 2
Ray M. Bowen Mechanical Engineering Texas A&M University College Station, Texas
and
C.-C. Wang Mathematical Sciences Rice University Houston, Texas
Copyright Ray M. Bowen and C.-C. Wang (ISBN 0-306-37509-5 (v. 2))
______PREFACE
To Volume 2
This is the second volume of a two-volume work on vectors and tensors. Volume 1 is concerned with the algebra of vectors and tensors, while this volume is concerned with the geometrical aspects of vectors and tensors. This volume begins with a discussion of Euclidean manifolds. The principal mathematical entity considered in this volume is a field, which is defined on a domain in a Euclidean manifold. The values of the field may be vectors or tensors. We investigate results due to the distribution of the vector or tensor values of the field on its domain. While we do not discuss general differentiable manifolds, we do include a chapter on vector and tensor fields defined on hypersurfaces in a Euclidean manifold.
This volume contains frequent references to Volume 1. However, references are limited to basic algebraic concepts, and a student with a modest background in linear algebra should be able to utilize this volume as an independent textbook. As indicated in the preface to Volume 1, this volume is suitable for a one-semester course on vector and tensor analysis. On occasions when we have taught a one –semester course, we covered material from Chapters 9, 10, and 11 of this volume. This course also covered the material in Chapters 0,3,4,5, and 8 from Volume 1.
We wish to thank the U.S. National Science Foundation for its support during the preparation of this work. We also wish to take this opportunity to thank Dr. Kurt Reinicke for critically checking the entire manuscript and offering improvements on many points.
Houston, Texas R.M.B. C.-C.W.
iii
______CONTENTS
Vol. 2 Vector and Tensor Analysis
Contents of Volume 1………………………………………………… vii
PART III. VECTOR AND TENSOR ANALYSIS
Selected Readings for Part III………………………………………… 296
CHAPTER 9. Euclidean Manifolds………………………………………….. 297
Section 43. Euclidean Point Spaces……………………………….. 297 Section 44. Coordinate Systems…………………………………… 306 Section 45. Transformation Rules for Vector and Tensor Fields…. 324 Section 46. Anholonomic and Physical Components of Tensors…. 332 Section 47. Christoffel Symbols and Covariant Differentiation…... 339 Section 48. Covariant Derivatives along Curves………………….. 353
CHAPTER 10. Vector Fields and Differential Forms………………………... 359
Section 49. Lie Derivatives……………………………………….. 359 Section 5O. Frobenius Theorem…………………………………… 368 Section 51. Differential Forms and Exterior Derivative………….. 373 Section 52. The Dual Form of Frobenius Theorem: the Poincaré Lemma……………………………………………….. 381 Section 53. Vector Fields in a Three-Dimensiona1 Euclidean Manifold, I. Invariants and Intrinsic Equations…….. 389 Section 54. Vector Fields in a Three-Dimensiona1 Euclidean Manifold, II. Representations for Special Class of Vector Fields………………………………. 399
CHAPTER 11. Hypersurfaces in a Euclidean Manifold
Section 55. Normal Vector, Tangent Plane, and Surface Metric… 407 Section 56. Surface Covariant Derivatives………………………. 416 Section 57. Surface Geodesics and the Exponential Map……….. 425 Section 58. Surface Curvature, I. The Formulas of Weingarten and Gauss…………………………………………… 433 Section 59. Surface Curvature, II. The Riemann-Christoffel Tensor and the Ricci Identities……………………... 443 Section 60. Surface Curvature, III. The Equations of Gauss and Codazzi 449
v vi CONTENTS OF VOLUME 2
Section 61. Surface Area, Minimal Surface………...... 454 Section 62. Surfaces in a Three-Dimensional Euclidean Manifold. 457
CHAPTER 12. Elements of Classical Continuous Groups
Section 63. The General Linear Group and Its Subgroups……….. 463 Section 64. The Parallelism of Cartan……………………………. 469 Section 65. One-Parameter Groups and the Exponential Map…… 476 Section 66. Subgroups and Subalgebras…………………………. 482 Section 67. Maximal Abelian Subgroups and Subalgebras……… 486
CHAPTER 13. Integration of Fields on Euclidean Manifolds, Hypersurfaces, and Continuous Groups
Section 68. Arc Length, Surface Area, and Volume……………... 491 Section 69. Integration of Vector Fields and Tensor Fields……… 499 Section 70. Integration of Differential Forms……………………. 503 Section 71. Generalized Stokes’ Theorem……………………….. 507 Section 72. Invariant Integrals on Continuous Groups…………... 515
INDEX………………………………………………………………………. x
______
CONTENTS
Vol. 1 Linear and Multilinear Algebra
PART 1 BASIC MATHEMATICS
Selected Readings for Part I………………………………………………………… 2
CHAPTER 0 Elementary Matrix Theory…………………………………………. 3
CHAPTER 1 Sets, Relations, and Functions……………………………………… 13
Section 1. Sets and Set Algebra………………………………………... 13 Section 2. Ordered Pairs" Cartesian Products" and Relations…………. 16 Section 3. Functions……………………………………………………. 18
CHAPTER 2 Groups, Rings and Fields…………………………………………… 23
Section 4. The Axioms for a Group……………………………………. 23 Section 5. Properties of a Group……………………………………….. 26 Section 6. Group Homomorphisms…………………………………….. 29 Section 7. Rings and Fields…………………………………………….. 33
PART I1 VECTOR AND TENSOR ALGEBRA
Selected Readings for Part II………………………………………………………… 40
CHAPTER 3 Vector Spaces……………………………………………………….. 41
Section 8. The Axioms for a Vector Space…………………………….. 41 Section 9. Linear Independence, Dimension and Basis…………….….. 46 Section 10. Intersection, Sum and Direct Sum of Subspaces……………. 55 Section 11. Factor Spaces………………………………………………... 59 Section 12. Inner Product Spaces………………………..………………. 62 Section 13. Orthogonal Bases and Orthogonal Compliments…………… 69 Section 14. Reciprocal Basis and Change of Basis……………………… 75
CHAPTER 4. Linear Transformations……………………………………………… 85
Section 15. Definition of a Linear Transformation………………………. 85 Section 16. Sums and Products of Linear Transformations……………… 93
viii CONTENTS OF VOLUME 2
Section 17. Special Types of Linear Transformations…………………… 97 Section 18. The Adjoint of a Linear Transformation…………………….. 105 Section 19. Component Formulas………………………………………... 118
CHAPTER 5. Determinants and Matrices…………………………………………… 125
Section 20. The Generalized Kronecker Deltas and the Summation Convention……………………………… 125 Section 21. Determinants…………………………………………………. 130 Section 22. The Matrix of a Linear Transformation……………………… 136 Section 23 Solution of Systems of Linear Equations…………………….. 142
CHAPTER 6 Spectral Decompositions……………………………………………... 145
Section 24. Direct Sum of Endomorphisms……………………………… 145 Section 25. Eigenvectors and Eigenvalues……………………………….. 148 Section 26. The Characteristic Polynomial………………………………. 151 Section 27. Spectral Decomposition for Hermitian Endomorphisms…….. 158 Section 28. Illustrative Examples…………………………………………. 171 Section 29. The Minimal Polynomial……………………………..……… 176 Section 30. Spectral Decomposition for Arbitrary Endomorphisms….….. 182
CHAPTER 7. Tensor Algebra………………………………………………………. 203
Section 31. Linear Functions, the Dual Space…………………………… 203 Section 32. The Second Dual Space, Canonical Isomorphisms…………. 213 Section 33. Multilinear Functions, Tensors…………………………..….. 218 Section 34. Contractions…...... 229 Section 35. Tensors on Inner Product Spaces……………………………. 235
CHAPTER 8. Exterior Algebra……………………………………………………... 247
Section 36. Skew-Symmetric Tensors and Symmetric Tensors………….. 247 Section 37. The Skew-Symmetric Operator……………………………… 250 Section 38. The Wedge Product………………………………………….. 256 Section 39. Product Bases and Strict Components……………………….. 263 Section 40. Determinants and Orientations………………………………. 271 Section 41. Duality……………………………………………………….. 280 Section 42. Transformation to Contravariant Representation……………. 287
INDEX…………………………………………………………………………………….x
______PART III
VECTOR AND TENSOR ANALYSIS
Selected Reading for Part III
BISHOP, R. L., and R. J. CRITTENDEN, Geometry of Manifolds, Academic Press, New York, 1964 BISHOP, R. L., and S. I. GOLDBERG, Tensor Analysis on Manifolds, Macmillan, New York, 1968. CHEVALLEY, C., Theory of Lie Groups, Princeton University Press, Princeton, New Jersey, 1946 COHN, P. M., Lie Groups, Cambridge University Press, Cambridge, 1965. EISENHART, L. P., Riemannian Geometry, Princeton University Press, Princeton, New Jersey, 1925. ERICKSEN, J. L., Tensor Fields, an appendix in the Classical Field Theories, Vol. III/1. Encyclopedia of Physics, Springer-Verlag, Berlin-Gottingen-Heidelberg, 1960. FLANDERS, H., Differential Forms with Applications in the Physical Sciences, Academic Press, New York, 1963. KOBAYASHI, S., and K. NOMIZU, Foundations of Differential Geometry, Vols. I and II, Interscience, New York, 1963, 1969. LOOMIS, L. H., and S. STERNBERG, Advanced Calculus, Addison-Wesley, Reading, Massachusetts, 1968. MCCONNEL, A. J., Applications of Tensor Analysis, Dover Publications, New York, 1957. NELSON, E., Tensor Analysis, Princeton University Press, Princeton, New Jersey, 1967. NICKERSON, H. K., D. C. SPENCER, and N. E. STEENROD, Advanced Calculus, D. Van Nostrand, Princeton, New Jersey, 1958. SCHOUTEN, J. A., Ricci Calculus, 2nd ed., Springer-Verlag, Berlin, 1954. STERNBERG, S., Lectures on Differential Geometry, Prentice-Hall, Englewood Cliffs, New Jersey, 1964. WEATHERBURN, C. E., An Introduction to Riemannian Geometry and the Tensor Calculus, Cambridge University Press, Cambridge, 1957.
______Chapter 9
EUCLIDEAN MANIFOLDS
This chapter is the first where the algebraic concepts developed thus far are combined with ideas from analysis. The main concept to be introduced is that of a manifold. We will discuss here only a special case cal1ed a Euclidean manifold. The reader is assumed to be familiar with certain elementary concepts in analysis, but, for the sake of completeness, many of these shall be inserted when needed.
Section 43 Euclidean Point Spaces
Consider an inner produce space V and a set E . The set E is a Euclidean point space if there exists a function f :EE×→ V such that:
(a) fff(,xy )=+ (,) xz (, zy ), xyz ,, ∈E and (b) For every x ∈E and v ∈V there exists a unique element y ∈E such that f (,x y )= v .
The elements of E are called points, and the inner product space V is called the translation space. We say that f (,x y ) is the vector determined by the end point x and the initial point y . Condition b) above is equivalent to requiring the function fx :EV→ defined by ffx (y )= (x ,y ) to be one to one for each x . The dimension of E , written dimE , is defined to be the dimension of V . If V does not have an inner product, the set E defined above is called an affine space.
A Euclidean point space is not a vector space but a vector space with inner product is made a Euclidean point space by defining f (,vv12 )≡ v 1− v 2 for all v ∈V . For an arbitrary point space the function f is called the point difference, and it is customary to use the suggestive notation
f (,x y )= x − y (43.1)
In this notation (a) and (b) above take the forms
297 298 Chap. 9 • EUCLIDEAN MANIFOLDS
x − y =−+−xzzy (43.2) and
x − y = v (43.3)
Theorem 43.1. In a Euclidean point space E
()i xx0−= ()ii xy−=−− ( yx ) (43.4) (iii ) if x −=−y x 'y ', then xx − ' =−yy '
Proof. For (i) take x ==y z in (43.2); then
xxxxxx− =−+− which implies xx0−=. To obtain (ii) take y = x in (43.2) and use (i). For (iii) observe that
x −=−+−=−+−y '''''x yyy xx x y from (43.2). However, we are given x − y =−x 'y ' which implies (iii).
The equation
x − y = v has the property that given any v and y , x is uniquely determined. For this reason it is customary to write
x = y + v (43.5)
Sec. 43 • Euclidean Point Spaces 299 for the point x uniquely determined by y ∈E and v ∈V .
The distance from x to y , written d (x ,y ) , is defined by
1/2 d(,xy )=−= xy{()() xy −⋅− xy} (43.6)
It easily fol1ows from the definition (43.6) and the properties of the inner product that
dd(,x yy )= (,)x (43.7) and
ddd(,x y )≤ (,)xz+ (, zy ) (43.8) for all x , y , and z in E . Equation (43.8) is simply rewritten in terms of the points x , y , and z rather than the vectors x −−y ,xz, and z − y . It is also apparent from (43.6) that
dd(,x y )≥=⇔= 0 and (,x y ) 0 x y (43.9)
The properties (43.7)-(43.9) establish that E is a metric space.
There are several concepts from the theory of metric spaces which we need to summarize. For simplicity the definitions are sated here in terms of Euclidean point spaces only even though they can be defined for metric spaces in general.
In a Euclidean point space E an open ball of radius ε > 0 centered at x0 ∈E is the set
Bd(,)xxxx00ε = { (,)< ε} (43.10) and a closed ball is the set
300 Chap. 9 • EUCLIDEAN MANIFOLDS
Bd(,)xxxx00ε = { (,)≤ ε} (43.11)
A neighborhood of x ∈E is a set which contains an open ball centered at x . A subset U of E is open if it is a neighborhood of each of its points. The empty set ∅ is trivially open because it contains no points. It also follows from the definitions that E is open.
Theorem 43.2. An open ball is an open set.
Proof. Consider the open ball B(,)x0 ε . Let x be an arbitrary point in B(,)x0 ε . Then
ε −>d(,xx0 ) 0 and the open ball Bd(,xxxε − (,0 )) is in B(,)x0 ε , because if y ∈−Bd(,xxxε (,0 )), then dd(,)y xxx<−ε (0 ,) and, by (43.8), ddd(,)x00y ≤ (,)xx+ (,) xy , which yields d(,)xy0 < ε and, thus, yx∈ B(,)0 ε .
A subset U of E is closed if its complement, EU, is open. It can be shown that closed balls are indeed closed sets. The empty set, ∅ , is closed because EE= ∅ is open. By the same logic E is closed since EE=∅ is open. In fact ∅ and E are the only subsets of E which are both open and closed. A subset UE⊂ is bounded if it is contained in some open ball. A subset UE⊂ is compact if it is closed and bounded.
Theorem 43.3. The union of any collection of open sets is open.
Proof. Let {Uα α ∈ I} be a collection of open sets, where I is an index set. Assume that x ∈∪ . Then x must belong to at least one of the sets in the collection, say . Since is α∈IUα Uα0 Uα0 open, there exists an open ball B(,)x ε ⊂⊂∪ . Thus, B (x ,ε ) ⊂ ∪ . Since x is UUα0 αα∈I α∈IUα arbitrary, ∪α∈IUα is open.
Theorem43.4. The intersection of a finite collection of open sets is open.
n Proof. Let {UU1,..., α } be a finite family of open sets. If ∩ii=1U is empty, the assertion is trivial. n n Thus, assume ∩ii=1U is not empty and let x be an arbitrary element of ∩ii=1U . Then x ∈Ui for in= 1,..., and there is an open ball B (x ,εii ) ⊂U for in= 1,..., . Let ε be the smallest of the n n positive numbers ε1,...,εn . Then xx∈⊂B(,)ε ∩ii=1U . Thus ∩ii=1U is open.
Sec. 43 • Euclidean Point Spaces 301
It should be noted that arbitrary intersections of open sets will not always lead to open sets. The standard counter example is given by the family of open sets of R of the form ()−1,1nn, ∞ n = 1,2,3.... The intersection ∩ n=1 ()−1,1nn is the set {0} which is not open.
By a sequence in E we mean a function on the positive integers {1,2,3,...,n ,...} with values in E . The notation {xxx123, , ..., xn ,...} , or simply {xn}, is usually used to denote the values of the sequence. A sequence {xn} is said to converge to a limit x ∈E if for every open ball B (x ,ε ) centered at x , there exists a positive integer n0 ()ε such that xn ∈ B (x ,ε ) whenever nn≥ 0 ()ε .
Equivalently, a sequence {xn} converges to x if for every real number ε > 0 there exists a positive integer n0 ()ε such that d (xxn , ) < ε for all nn> 0 ()ε . If {xn} converges to x , it is conventional to write
xx=→→∞ limnn or xx as n n→∞
Theorem 43.5. If {xn} converges to a limit, then the limit is unique.
Proof. Assume that
xx= limnn ,y = lim x nn→∞ →∞
Then, from (43.8)
dd(,x y )≤ (,xxnn )+ d ( x ,y ) for every n . Let ε be an arbitrary positive real number. Then from the definition of convergence of {xn}, there exists an integer n0 ()ε such that nn≥ 0 ()ε implies d (xx ,n ) < ε and d (xyn , ) < ε . Therefore,
d(,xy )≤ 2ε for arbitrary ε . This result implies d(,xy )= 0 and, thus, x = y . 302 Chap. 9 • EUCLIDEAN MANIFOLDS
A point x ∈E is a limit point of a subset UE⊂ if every neighborhood of x contains a point of U distinct from x . Note that x need not be in U . For example, the sphere
{xxxd(,0 )= ε} are limit points of the open ball B(,)x0 ε . The closure of UE⊂ , written U , is the union of U and its limit points. For example, the closure of the open ball B(,)x0 ε is the closed ball B(,)x0 ε . It is a fact that the closure of U is the smallest closed set containing U . Thus U is closed if and only if UU= .
The reader is cautioned not to confuse the concepts limit of a sequence and limit point of a subset. A sequence is not a subset of E ; it is a function with values in E . A sequence may have a limit when it has no limit point. Likewise the set of pints which represent the values of a sequence may have a limit point when the sequence does not converge to a limit. However, these two concepts are related by the following result from the theory of metric spaces: a point x is a limit point of a set U if and only if there exists a convergent sequence of distinct points of U with x as a limit.
A mapping f :'UE→ , where U is an open set in E and E ' is a Euclidean point space or an inner product space, is continuous at x0 ∈U if for every real number ε > 0 there exists a real number δ (,)0x0 ε > such that d(,xx00 )< δ (,ε x ) implies df'( (xx0 ), f ( )) < ε . Here d ' is the distance function for E ' . When f is continuous at x0 , it is conventional to write
limfff (xxxxx ) or ( )→→ (00 ) as xx→ 0
The mapping f is continuous on U if it is continuous at every point of U . A continuous mapping is called a homomorphism if it is one-to-one and if its inverse is also continuous. What we have just defined is sometimes called a homeomorphism into. If a homomorphism is also onto, then it is called specifically a homeomorphism onto. It is easily verified that a composition of two continuous maps is a continuous map and the composition of two homomorphisms is a homomorphism.
A mapping f :UE→ ' , where U and E ' are defined as before, is differentiable at x ∈U if there exists a linear transformation Ax ∈LV ( ; V ') such that
ff()()xv+= x + Avx + o (,) xv (43.12)
Sec. 43 • Euclidean Point Spaces 303 where
o(,xv ) lim = 0 (43.13) v0→ v
In the above definition V ' denotes the translation space of E ' .
Theorem 43.6. The linear transformation Ax in (43.12) is unique.
Proof. If (43.12) holds for Ax and Ax , then by subtraction we find
()AAvxx−=oo(, xv ) − (, xv )
By (43.13), ()AAexx− must be zero for each unit vector e , so that AAxx= .
If f is differentiable at every point of U , then we can define a mapping gradf :ULVV→ ( ; ') , called the gradient of f , by
gradf (xA )= x , x∈U (43.14)
If grad f is continuous on U , then f is said to be of class C1 . If grad f exists and is itself of class C1 , then f is of class C 2 . More generally, f is of class C r , r > 0 , if it is of class C r−1 and its (1)r − st gradient, written gradr−1 f , is of class C1 . Of course, f is of class C 0 if it is continuous on U . If f is a C r one-to-one map with a C r inverse f −1 defined on f ()U , then f is called a C r diffeomorphism.
If f is differentiable at x , then it follows from (43.12) that
ffd()()xv+−τ x Avx ==+limf ( xτ v ) (43.15) τ →0 ττd τ =0
304 Chap. 9 • EUCLIDEAN MANIFOLDS for all v ∈V . To obtain (43.15) replace v by τ v , τ > 0 in (43.12) and write the result as
ff()()xv+−τ xo(,xvτ ) Av=− (43.16) x ττ
By (43.13) the limit of the last term is zero as τ → 0, and (43.15) is obtained. Equation (43.15) holds for all v ∈V because we can always choose τ in (43.16) small enough to ensure that xv+τ is in U , the domain of f . If f is differentiable at every x ∈U , then (43.15) can be written
d ()gradff (xv )=+ ( xτ v ) (43.17) dτ τ =0
A function f :UR→ , where U is an open subset of E , is called a scalar field. Similarly, f :UV→ is a vector field, and f :UTV→ q ( ) is a tensor field of order q . It should be noted that the term field is defined here is not the same as that in Section 7.
Before closing this section there is an important theorem which needs to be recorded for later use. We shall not prove this theorem here, but we assume that the reader is familiar with the result known as the inverse mapping theorem in multivariable calculus.
r Theorem 43.7. Let f :'UE→ be a C mapping and assume that gradf (x0 ) is a linear isomorphism. Then there exists a neighborhood U1 of x0 such that the restriction of f to U1 is a C r diffeomorphism. In addition
−1 −1 gradff ( (xx00 ))= () grad f ( ) (43.18)
This theorem provides a condition under which one can asert the existence of a local inverse of a smooth mapping.
Exercises
43.1 Let a sequence {xn} converge to x . Show that every subsequence of {xn} also converges to x . Sec. 43 • Euclidean Point Spaces 305
43.2 Show that arbitrary intersections and finite unions of closed sets yields closed sets. 43.3 Let f :UE→ ', where U is open in E , and E ' is either a Euclidean point space or an inner produce space. Show that f is continuous on U if and only if f −1()D is open in E for all D open in f (U ) . 43.4 Let f :UE→ ' be a homeomorphism. Show that f maps any open set in U onto an open set in E ' . 43.5 If f is a differentiable scalar valued function on LV ( ; V ) , show that the gradient of f at A ∈LV( ; V ) , written
∂f ()A ∂A
is a linear transformation in LV(;) V defined by
⎛⎞∂fdfT tr⎜⎟ (AB ) =+() Aτ B ⎝⎠∂A dτ τ =0
for all B ∈LV ( ; V ) 43.6 Show that
∂∂μ μ T 1 ()AI== andN () A() adj A ∂∂AA
306 Chap. 9 • EUCLIDEAN MANIFOLDS
Section 44 Coordinate Systems
Given a Euclidean point space E of dimension N , we define a Cchartr - at x ∈E to be a pair ()U , xˆ , where U is an open set in E containing x and xˆ :UR→ N is a C r diffeomorphism. Given any chart ()U , xˆ , there are N scalar fields xˆi :UR→ such that
xxˆˆ(xx )= ( 1 ( ),..., x ˆN ( x )) (44.1) for all x ∈U . We call these fields the coordinate functions of the chart, and the mapping xˆ is also called a coordinate map or a coordinate system on U . The set U is called the coordinate neighborhood.
N N Two charts xˆ :UR1 → and yˆ :UR2 → , where UU12∩ ≠∅, yield the coordinate −1 −1 transformation yxˆˆ :( x ˆUU12∩→ )y ˆ ( UU 12 ∩ ) and its inverse xyˆˆ :(y ˆUU12∩→ ) x ˆ ( UU 12 ∩ ). Since
yyˆˆ(xx )= ( 1 ( ),...,y ˆN ( x )) (44.2)
The coordinate transformation can be written as the equations
()yyˆˆ111(xx ),...,NN ( )= yxxx ˆˆˆˆ − ( ( xx ),..., ( )) (44.3) and the inverse can be written
()xxˆˆ111(xx ),...,NN ( )= xyyy ˆˆˆˆ − ( ( xx ),..., ( )) (44.4)
The component forms of (44.3) and (44.4) can be written in the simplified notation
y jj=≡yx(1 ,..., x N ) yx jk ( ) (44.5)
Sec. 44 • Coordinate Systems 307 and
x jj=≡xy(1 ,..., y N ) xy jk ( ) (44.6)
The two N-tuples (y1 ,...,y N ) and (x1 ,...,x N ) , where yyjj= ˆ ()x and xxjj= ˆ ()x , are the coordinates of the point x ∈∩UU12. Figure 6 is useful in understanding the coordinate transformations. R N
UU12∩
U1 U2
xˆ yˆ xˆ()∩ UU12 yˆ()UU12∩ yˆˆ x−1
xˆˆ y−1
Figure 6
It is important to note that the quantities x1,..., x N , y1,..., y N are scalar fields, i.e., real- valued functions defined on certain subsets of E . Since U1 and U2 are open sets, UU12∩ is open, r and because xˆ and yˆ are C diffeomorphisms, xˆ (UU12∩ ) and yˆ (UU12∩ ) are open subsets of R N . In addition, the mapping yˆˆ x−1 and xˆˆ y−1 are C r diffeomorphisms. Since xˆ is a diffeomorphism, equation (44.1) written in the form
xˆ(x )= (xx1 ,...,N ) (44.7) can be inverted to yie1d
308 Chap. 9 • EUCLIDEAN MANIFOLDS
xx== (x1 ,...,xxNj ) x ( ) (44.8) where x is a diffeomorphism x :()xˆ UU→ .
r r A C- atlas on E is a family (not necessarily countable) of C- charts {()Uαα, xˆ α ∈ I} , where I is an index set, such that
EU= ∪ α (44.9) α∈I
r Equation (44.9) states that E is covered by the family of open sets {Uα α ∈ I}. A C -Euclidean manifold is a Euclidean point space equipped with a Cr - atlas . A Catlas∞ - and a C ∞ -Euclidean manifold are defined similarly. For simplicity, we shall assume that E is C ∞ .
A C ∞ curve in E is a C ∞ mapping λ : (ab , ) →E , where (ab , ) is an open interval of R . ∞ A C curve λ passes through x0 ∈E if there exists a cab∈ ( , ) such that λ()c = x0 . Given a chart th ()U , xˆ and a point x0 ∈U , the j coordinate curve passing through x0 is the curve λ j defined by
11jj−+ j 1 N λ j (txxxtxx )=+x (0000 ,..., , , ,..., 0 ) (44.10)
11jj−+ j 1 N k for all t such that (xxxtxx0000 ,..., ,+∈ , ,..., 0 ) xˆ (U ) , where ()xx00= ˆ ()x . The subset of U obtained by requiring
xxjj==ˆ ()x const (44.11) is called the jth coordinate surface of the chart
Euclidean manifolds possess certain special coordinate systems of major interest. Let {ii1,..., N } be an arbitrary basis, not necessarily orthonormal, for V . We define N constant vector fields i j :EV→ , jN= 1,..., , by the formulas
Sec. 44 • Coordinate Systems 309
ixjj()= i , x∈U (44.12)
The use of the same symbol for the vector field and its value will cause no confusion and simplifies the notation considerably. If 0E denotes a fixed element of E , then a Cartesian coordinate system on E is defined by the N scalar fields zzˆˆ12, ,..., z ˆN such that
jj j zz==−⋅ˆ ()xx0i (E ) , x ∈E (44.13)
If the basis {ii1,..., N } is orthonormal, the Cartesian system is called a rectangular Cartesian system. The point 0E is called the origin of the Cartesian coordinate system. The vector field defined by
rx()= x− 0E (44.14)
for all x ∈E is the position vector field relative to 0E . The value rx ( ) is the position vector of x . 1 N If {ii1,..., N } is the basis reciprocal to {ii,..., } , then (44.13) implies
1 Njj x0−=EE x (zz ,..., ) −= 0 zˆ ( xi ) jj = z i (44.15)
Defining constant vector fields ii1,..., N as before, we can write (44.15) as
j ri= zˆ j (44.16)
j The product of the scalar field zˆ with the vector field i j in (44.16) is defined pointwise; i.e., if f is a scalar field and v is a vector field, then fv is a vector field defined by
ffvx()= ()() xvx (44.17) for all x in the intersection of the domains of f and v . An equivalent version of (44.13) is
zˆ jj= ri⋅ (44.18) 310 Chap. 9 • EUCLIDEAN MANIFOLDS
where the operator ri⋅ j between vector fields is defined pointwise in a similar fashion as in (44.17). As an illustration, let {ii12, ,..., iN } be another basis for V related to the original basis by
jjk ii= Qk (44.19)
and let 0E be another fixed point of E . Then by (44.13)
z jjj=−()()(x0 ⋅=− i x0 ⋅+ i 0 − 0 ) ⋅ i j EEEE (44.20) jk jjkj =−⋅+−⋅=+QQzckk()(x0EEE i 0 0 ) i
where (44.19) has been used. Also in (44.20) the constant scalars ck , kN= 1,..., , are defined by
kk c = ()00iEE−⋅ (44.21)
12 N 12 N ⎡⎤j If the bases {ii, ,..., i} and {ii, ,..., i} are both orthonormal, then the matrix ⎣⎦Qk is orthogonal. Note that the coordinate neighborhood is the entire space E .
th The j coordinate curve which passes through 0E of the Cartesian coordinate system is the curve
λ()tt= i0j + E (44.22)
12 N Equation (44.22) follows from (44.10), (44.15), and the fact that for x0= E , zz= =⋅⋅⋅= z = 0 . th As (44.22) indicates, the coordinate curves are straight lines passing through 0E . Similarly, the j coordinate surface is the plane defined by Sec. 44 • Coordinate Systems 311
j ()x0−⋅=E i const
z3
i3
z2 0E i2
i1
z1
Figure 7
Geometrically, the Cartesian coordinate system can be represented by Figure 7 for N=3. The coordinate transformation represented by (44.20) yields the result in Figure 8 (again for N=3). Since every inner product space has an orthonormal basis (see Theorem 13.3), there is no loss of generality in assuming that associated with every point of E as origin we can introduce a rectangular Cartesian coordinate system. 312 Chap. 9 • EUCLIDEAN MANIFOLDS
3 z 3 z
r r i 3 i3
2 z 0 0 i 1 E 2 E 2 i i i 1 z 2 z 1 z1
Figure 8
Given any rectangular coordinate system (zzˆˆ1 ,...,N ) , we can characterize a general or a curvilinear coordinate system as follows: Let (U , xˆ) be a chart. Then it can be specified by the coordinate transformation from zˆ to xˆ as described earlier, since in this case the overlap of the coordinate neighborhood is UEU=∩. Thus we have
(zz111 ,...,NN )= zxxxˆ ˆ − ( ,..., ) (44.23) (x111 ,...,xxzzzNN )= ˆ ˆ− ( ,..., ) where zxˆˆ ˆˆ−1 :() xUU→ z () and xzˆˆ ˆˆ−1 :( zUU )→ x ( ) are diffeomorphisms. Equivalent versions of (44.23) are
zzxxjj==(1 ,..., N ) zx jk ( ) (44.24) x jj==xz(1 ,..., z N ) xz jk ( )
Sec. 44 • Coordinate Systems 313
As an example of the above ideas, consider the cylindrical coordinate system1. In this case N=3 and equations (44.23) take the special form
(,zzz123 , )= ( x 1 cos, xx 21 sin, xx 2 3 ) (44.25)
In order for (44.25) to qualify as a coordinate transformation, it is necessary for the transformation functions to be C ∞ . It is apparent from (44.25) that zxˆ ˆ −1 is C ∞ on every open subset of R 3 . Also, by examination of (44.25), zxˆ ˆ −1 is one-to-one if we restrict it to an appropriate domain, say (0,∞× ) (0,2π ) ×−∞∞ ( , ) . The image of this subset of R 3 under zxˆ ˆ −1 is easily seen to be the set R 31231{(,zzz , ) z≥= 0, z 2 0} and the inverse transformation is
2 ⎛⎞1/2 z (,x123xx , )=+⎡⎤ () z 12 () z 22 ,tan− 1 , z 3 (44.26) ⎜⎟⎣⎦1 ⎝⎠z which is also C ∞ . Consequently we can choose the coordinate neighborhood to be any open subset U in E such that
zzzzzzˆ()UR⊂≥=31231{ (, , ) 0, 2 0} (44.27) or, equivalently,
xˆ()U ⊂∞× (0,)(0,2)(π ×−∞∞ ,)
Figure 9 describes the cylindrical system. The coordinate curves are a straight line (for x1 ), a circle lying in a plane parallel to the (,zz12 ) plane (for x2 ), and a straight line coincident with z3 (for x3 ). The coordinate surface x1 = const is a circular cylinder whose generators are the z3 lines. The remaining coordinate surfaces are planes.
1 The computer program Maple has a plot command, coordplot3d, that is useful when trying to visualize coordinate curves and coordinate surfaces. The program MATLAB will also produce useful and instructive plots. 314 Chap. 9 • EUCLIDEAN MANIFOLDS
zx33=
z2
x2 z1 x1 Figure 9
Returning to the general transformations (44.5) and (44.6), we can substitute the second into the first and differentiate the result to find
∂∂yxik (xx11 ,...,NNi ) ( yy ,..., ) = δ (44.28) ∂∂xykj j
By a similar argument with x and y interchanged,
∂∂yxki (xx11 ,...,NNi ) ( yy ,..., ) = δ (44.29) ∂∂xyjk j
Each of these equations ensures that
⎡⎤∂yi 1 det (xx1 ,...,N )= ≠ 0 (44.30) ⎢⎥k j ⎣⎦∂x ⎡⎤∂x 1 N det⎢⎥l (yy ,..., ) ⎣⎦∂y
⎡⎤jk123 1 For example, det⎣⎦()()∂∂zxxxx , , =≠ x 0 for the cylindrical coordinate system. The ⎡⎤jk1 N determinant det⎣⎦()()∂∂xyy ,..., y is the Jacobian of the coordinate transformation (44.6).
Sec. 44 • Coordinate Systems 315
Just as a vector space can be assigned an orientation, a Euclidean manifold can be oriented by assigning the orientation to its translation space V . In this case E is called an oriented Euclidean manifold. In such a manifold we use Cartesian coordinate systems associated with positive basis only and these coordinate systems are called positive. A curvilinear coordinate system is positive if its coordinate transformation relative to a positive Cartesian coordinate system has a positive Jacobian.
Given a chart ()U , xˆ for E , we can compute the gradient of each coordinate function xˆi and obtain a C ∞ vector field on U . We shall denote each of these fields by gi , namely
gii= grad xˆ (44.31) for iN= 1,..., . From this definition, it is clear that gi (x ) is a vector in V normal to the ith coordinate surface. From (44.8) we can define N vector fields gg1,..., N on U by
gxi = [grad ] (0,...,1 ,...,0) (44.32) i
Or, equivalently,
11iN N xxx (xxtx ,...,+− ,..., ) ( xx ,..., ) ∂1 N gi =≡lim (x ,...,x ) (44.33) t→0 tx∂ i
th for all x ∈U . Equations (44.10) and (44.33) show that gi (x ) is tangent to the i coordinate curve. Since
xˆiNi(x (xx1 ,..., )) = x (44.34)
The chain rule along with the definitions (44.31) and (44.32) yield
ii gx()⋅ gjj () x= δ (44.35) as they should. 316 Chap. 9 • EUCLIDEAN MANIFOLDS
The values of the vector fields {gg1,..., N } form a linearly independent set of vectors
{g1(x ),...,g N (x )} at each x ∈U . To see this assertion, assume x ∈U and
12 N λλg12()x + g ()x +⋅⋅⋅+ λgN ()x0 = for λλ1,..., N ∈R . Taking the inner product of this equation with g j (x ) and using equation (44.35) , we see that λ j = 0 , jN= 1,..., which proves the assertion.
Because V has dimension N, {g1(x ),...,g N (x )} forms a basis for V at each x ∈U . 1 N Equation (44.35) shows that {g (x ),...,g (x )} is the basis reciprocal to {g1(x ),...,g N (x )}. Because 1 N of the special geometric interpretation of the vectors {g1(x ),...,g N (x )} and {g (x ),...,g (x )} mentioned above, these bases are called the natural bases of xˆ at x . Any other basis field which cannot be determined by either (44.31) or (44.32) relative to any coordinate system is called an 1 N anholonomic or nonintegrable basis. The constant vector fields {ii1,..., N } and {ii,..., } yield the natural bases for the Cartesian coordinate systems.
If ()U1, xˆ and ()U2 , yˆ are two charts such that UU12∩ ≠∅, we can determine the transformation rules for the changes of natural bases at x ∈UU12∩ in the following way: We shall let the vector fields h j , j = 1,..., N , be defined by
h jj= grad yˆ (44.36)
Then, from (44.5) and (44.31),
∂y j hxjj( )== gradyxxxˆˆ ( x ) (1 ,..., Ni )grad ( x ) ∂xi (44.37) ∂y j = (xx1 ,...,Ni )gx ( ) ∂xi
for all x ∈∩UU12. A similar calculation shows that
Sec. 44 • Coordinate Systems 317
∂xi hx( )= (yy1 ,...,N )g (x ) (44.38) ji∂y j
for all x ∈∩UU12. Equations (44.37) and (44.38) are the desired transformations.
2 ij Given a chart ()U1, xˆ , we can define 2N scalar fields gij :UR→ and g :UR→ by
ggij()x = g i ()x ⋅=g j ()xx ji () (44.39) and
ggij()x =⋅g i ()x g j ()xx = ji () (44.40) for all x ∈U . It immediately follows from (44.35) that
−1 ⎡⎤ij ⎡ ⎤ ⎣⎦gg()xx= ⎣ ij ()⎦ (44.41) since we have
iij gg= g j (44.42) and
i ggjji= g (44.43)
where the produce of the scalar fields with vector fields is defined by (44.17). If θij is the angle between the ith and the jth coordinate curves at x ∈U , then from (44.39)
gij ()x cosθij = 1/2 (no sum) (44.44) ⎡⎤ ⎣⎦ggii()xx jj () 318 Chap. 9 • EUCLIDEAN MANIFOLDS
Based upon (44.44), the curvilinear coordinate system is orthogonal if gij = 0 when ij≠ . The symbol g denotes a scalar field on U defined by
⎡ ⎤ gg()xx= det⎣ ij ()⎦ (44.45) for all x ∈U .
At the point x ∈U , the differential element of arc ds is defined by
ds2 = dxx⋅ d (44.46) and, by (44.8), (44.33), and (44.39),
∂xx ∂ ds2 =⋅(()) xˆˆxx (()) x dxij dx ∂∂xxij (44.47) ij = gdxdxij ()x
If ()U1, xˆ and ()U2 , yˆ are charts where UU12∩ ≠∅, then at x ∈UU12∩
hij()xhxhx=⋅ i () j () ∂∂xxkl (44.48) = (yy11 ,...,NN ) ( yyg ,..., ) (x ) ∂∂yyijkl
Equation (44.48) is helpful for actual calculations of the quantities gij (x ) . For example, for the transformation (44.23), (44.48) can be arranged to yield
∂∂zzkk gxxxx(x )= (11 ,...,NN ) ( ,..., ) (44.49) ij ∂∂xxij
since iikl⋅=δ kl. For the cylindrical coordinate system defined by (44.25) a simple calculation based upon (44.49) yields Sec. 44 • Coordinate Systems 319
⎡100⎤ ⎡⎤gx()x = ⎢ 0 (12 ) 0⎥ (44.50) ⎣⎦ij ⎢ ⎥ ⎣⎢001⎦⎥
Among other things, (44.50) shows that this coordinate system is orthogonal. By (44.50) and (44.41)
⎡100⎤ ⎡⎤gxij ()x = ⎢ 0 1(12 ) 0⎥ (44.51) ⎣⎦⎢ ⎥ ⎣⎢001⎦⎥
And, from (44.50) and (44.45),
gx()x = (12 ) (44.52)
It follows from (44.50) and (44.47) that
ds212122232=+()()()() dx x dx + dx (44.53)
Exercises
44.1 Show that
∂x ii==andkk grad zˆ k ∂zk
for any Cartesian coordinate system zˆ associated with {i j}. 44.2 Show that
Irx= grad ( )
320 Chap. 9 • EUCLIDEAN MANIFOLDS
44.3 Spherical coordinates (,x123xx , ) are defined by the coordinate transformation
zx11= sin x 2 cos x 3 zx21= sin x 2 sin x 3 zx31= cos x 2
relative to a rectangular Cartesian coordinate system zˆ . How must the quantity (,x123xx , ) be restricted so as to make zxˆ ˆ −1 one-to-one? Discuss the coordinate curves and the coordinate surfaces. Show that
⎡10 0 ⎤ ⎡⎤gx()x = ⎢ 0 (12 ) 0 ⎥ ⎣⎦ij ⎢ ⎥ 122 ⎣⎢00(sin)xx⎦⎥
44.4 Paraboloidal coordinates (,x123xx , ) are defined by
zxxx1123= cos zxxx2123= sin 1 zxx31222=−()() () 2
Relative to a rectangular Cartesian coordinate system zˆ . How must the quantity (,x123xx , ) be restricted so as to make zxˆ ˆ −1 one-to-one. Discuss the coordinate curves and the coordinate surfaces. Show that
⎡⎤()xx12+ () 22 0 0 ⎢⎥ ⎡⎤ 12 22 ⎣⎦gxxij ()x =+⎢⎥ 0 ( ) ( ) 0 ⎢⎥122 ⎣⎦00()xx
44.5 A bispherical coordinate system (,x123xx , ) is defined relative to a rectangular Cartesian coordinate system by Sec. 44 • Coordinate Systems 321
axsin23 cos x z1 = coshx12− cos x
axsin23 sin x z2 = coshx12− cos x
and
axsinh 1 z3 = coshx12− cos x
where a > 0 . How must (,x123xx , ) be restricted so as to make zxˆ ˆ −1 one-to-one? Discuss the coordinate curves and the coordinate surfaces. Also show that
⎡⎤2 ⎢⎥a 2 00 ⎢⎥()coshxx12− cos ⎢⎥ ⎢⎥a2 ⎡⎤gij ()x = 0 0 ⎣⎦⎢⎥122 ⎢⎥()coshxx− cos ⎢⎥ ax222(sin ) ⎢⎥00 ⎢⎥122 ⎣⎦()coshxx− cos
44.6 Prolate spheroidal coordinates (,x123xx , ) are defined by
za1123= sinh x sin x cos x za2123= sinh x sin x sin x za312= cosh x cos x
relative to a rectangular Cartesian coordinate system zˆ , where a > 0 . How must (,x123xx , ) be restricted so as to make zxˆ ˆ −1 one-to-one? Also discuss the coordinate curves and the coordinate surfaces and show that
322 Chap. 9 • EUCLIDEAN MANIFOLDS
⎡⎤axx22122()cosh− cos 0 0 ⎢⎥ ⎢⎥22122 ⎡⎤gaxxij (x )=− 0() cosh cos 0 ⎣⎦⎢⎥ ⎢⎥0 0axx22122 sinh sin ⎣⎦
44.7 Elliptical cylindrical coordinates (,x123xx , ) are defined relative to a rectangular Cartesian coordinate system by
za112= cosh x cos x za212= sinh x sin x zx33=
where a > 0 . How must (,x123xx , ) be restricted so as to make zxˆ ˆ −1 one-to-one? Discuss the coordinate curves and coordinate surfaces. Also, show that
⎡⎤axx22122()sinh+ sin 0 0 ⎢⎥ ⎢⎥22122 ⎡⎤gaxxij ()x =+ 0() sinh sin 0 ⎣⎦⎢⎥ ⎢⎥001 ⎣⎦
44.8 For the cylindrical coordinate system show that
22 gii112=+(cosxx ) (sin ) 12 1 2 g212=−xx(sin )ii + x (cos x )
gi33=
44.9 At a point x in E , the components of the position vector rx ( ) = x− 0E with respect to the 1 N basis {ii1,..., N } associated with a rectangular Cartesian coordinate system are zz,..., . This observation follows, of course, from (44.16). Compute the components of rx() with
respect to the basis {g123(),(),()x g x g x } for (a) cylindrical coordinates, (b) spherical coordinates, and (c) parabolic coordinates. You should find that
Sec. 44 • Coordinate Systems 323
13 rx()=+xx g13 () x g () x for (a) 1 rx()= x g1 () x for (b) 11 rx()=+xx12 g () x g () x for (c) 2212
44.10 Toroidal coordinates (,x123xx , ) are defined relative to a rectangular Cartesian coordinate system by
axxsinh13 cos z1 = coshx12− cos x
axxsinh13 sin z2 = coshx12− cos x
and
axsin 2 z3 = coshx12− cos x
where a > 0 . How must (,x123xx , ) be restricted so as to make zxˆ ˆ −1 one to one? Discuss the coordinate surfaces. Show that
⎡ 2 ⎤ ⎢ a ⎥ 2 00 ⎢()coshxx12− cos ⎥ ⎢ ⎥ ⎢ a2 ⎥ ⎡⎤gij ()x = 0 0 ⎣⎦⎢ 122 ⎥ ⎢ ()coshxx− cos ⎥ ⎢ ⎥ ax222sinh ⎢ 00 ⎥ ⎢ 122 ⎥ ⎣ ()cosxx− cos ⎦
324 Chap. 9 • EUCLIDEAN MANIFOLDS
Section 45. Transformation Rules for Vectors and Tensor Fields
In this section, we shall formalize certain ideas regarding fields on E and then investigate the transformation rules for vectors and tensor fields. Let U be an open subset of E ; we shall denote by F ∞ (U ) the set of C ∞ functions f :UR→ . First we shall study the algebraic structure ∞ ∞ ∞ of F (U ) . If f1 and f2 are in F (U ) , then their sum f12+ f is an element of F (U ) defined by
(ff12+ )()xxx=+ f 1 () f 2 () (45.1)
∞ and their produce f12f is also an element of F (U ) defined by
(ff12 )()xxx= f 1 () f 2 () (45.2) for all x ∈U . For any real number λ ∈R the constant function is defined by
λ()x = λ (45.3) for all x ∈U . For simplicity, the function and the value in (45.3) are indicated by the same symbol. Thus, the zero function in F ∞ (U ) is denoted simply by 0 and for every fF∈ ∞ (U )
f + 0= f (45.4)
It is also apparent that
1f = f (45.5)
In addition, we define
− f =−(1)f (45.6)
Sec. 45 • Transformation Rules 325
It is easily shown that the operations of addition and multiplication obey commutative, associative, and distributive laws. These facts show that F ∞ (U ) is a commutative ring (see Section 7).
An important collection of scalar fields can be constructed as follows: Given two charts 2 ij1 N ()U1, xˆ and ()U2 , yˆ , where UU12∩≠∅, we define the N partial derivatives ()()∂∂y xx,..., x 1 N 2 ∞ at every ()xx,..., ∈∩ xˆ ()UU12. Using a suggestive notation, we can define N C functions ij ∂∂yx:UU12 ∩→ R by
∂∂yyii ()xx= xˆ () (45.7) ∂∂xxjj
for all x ∈∩UU12.
As mentioned earlier, a C ∞ vector field on an open set U of E is a C ∞ map v :UV→ , where V is the translation space of E . The fields defined by (44.31) and (44.32) are special cases of vector fields. We can express v in component forms on UU12∩ ,
ij v ==υ ggijυ (45.8)
i υ jji= g υ (45.9)
i As usual, we can computer υ :UU∩→1 R by
ii υ ()xvxgx=⋅ () (), x ∈∩UU12 (45.10)
and υi by (45.9). In particular, if ()U2 , yˆ is another chart such that UU12∩ ≠∅, then the component form of h j relative to xˆ is
∂xi h = g (45.11) ji∂y j
326 Chap. 9 • EUCLIDEAN MANIFOLDS
With respect to the chart ()U2 , yˆ , we have also
k vh=υ k (45.12)
k where υ :UU∩→2 R. From (45.12), (45.11), and (45.8), the transformation rule for the components of v relative to the two charts (U1, xˆ) and (U2 , yˆ ) is
∂xi υ ij= υ (45.13) ∂y j
for all x ∈∩∩UUU12 .
As in (44.18), we can define an inner product operation between vector fields. If v11:UV→ and v22:UV→ are vector fields, then vv12⋅ is a scalar field defined on UU12∩ by
vvx12⋅=⋅() vxvx 1 () 2 (), x ∈∩UU 1 2 (45.14)
Then (45.10) can be written
υ ii= v ⋅g (45.15)
∞ Now let us consider tensor fields in general. Let Tq ()U denote the set of all tensor fields of ∞ ∞ order q defined on an open set U in E . As with the set F (U ) , the set Tq ()U can be assigned ∞ an algebraic structure. The sum of A :()UTV→ q and B :()UTV→ q is a C tensor field
AB+→:()UTVq defined by
()AB+=+()xxx A () B () (45.16)
∞ ∞ ∞ for all x∈U . If fF∈ (U ) and A ∈Tq ()U , then we can define fTA ∈ q ()U by
Sec. 45 • Transformation Rules 327
ffA()xxx= ()A () (45.17)
Clearly this multiplication operation satisfies the usual associative and distributive laws with ∞ ∞ respect to the sum for all x∈U . As with F (U ) , constant tensor fields in Tq ()U are given the same symbol as their value. For example, the zero tensor field is 0 :UTV→ q ( ) and is defined by
0()x0= (45.18) for all x∈U . If 1 is the constant function in F ∞ (U ) , then
−A =−(1)A (45.19)
∞ ∞ The algebraic structure on the set Tq ()V just defined is called a module over the ring F (U ) .
q The components of a tensor field A :UTV→ q ( ) with respect to a chart ()U1, xˆ are the N scalar fields A : ∩→ defined by ii1...q UU 1 R
A (xx )= A ( )(g (x ),...,g (x )) (45.20) ii11... qq i i
for all x∈∩UU1 . Clearly we can regard tensor fields as multilinear mappings on vector fields with values as scalar fields. For example, A(gg ,..., ) is a scalar field defined by ii1 q
AA(gg ,..., )(xx )= ( )(g (x ),...,g (x )) ii11qq i i
for all x∈∩UU1 . In fact we can, and shall, carry over to tensor fields the many algebraic operations previously applied to tensors. In particular a tensor field A :UTV→ q ( ) has the representation
i A =⊗⋅⋅⋅⊗A ggi1 q (45.21) ii1... q
328 Chap. 9 • EUCLIDEAN MANIFOLDS
for all x∈∩UU1 , where U1 is the coordinate neighborhood for a chart (U1, xˆ) . The scalar fields A are the covariant components of A and under a change of coordinates obey the ii1... q transformation rule
i ∂∂xxi1 q AAkk...= ⋅⋅⋅ k ii ... (45.22) 11qq∂yk1 ∂y q
Equation (45.22) is a relationship among the component fields and holds at all points x∈∩∩UUU12 where the charts involved are (U1, xˆ) and (U2 , yˆ ) . We encountered an example of
(45.22) earlier with (44.48). Equation (44.48) shows that the gij are the covariant components of a tensor field I whose value is the identity or metric tensor, namely
ij jj ij I =ggijgggggg ⊗=⊗=⊗= j j gg i ⊗ j (45.23)
for points x∈U1 , where the chart in question is (U1, xˆ) . Equations (45.23) show that the components of a constant tensor field are not necessarily constant scalar fields. It is only in Cartesian coordinates that constant tensor fields have constant components.
Another important tensor field is the one constructed from the positive unit volume tensor
E . With respect to an orthonormal basis {i j }, which has positive orientation, E is given by (41.6), i.e.,
E =∧⋅⋅⋅∧=iiε i ⊗⋅⋅⋅⊗ i (45.24) 1 Niii11⋅⋅⋅ NN i
ˆ Given this tensor, we define as usual a constant tensor field E :()ETV→ N by
EE()x = (45.25)
for all x ∈E . With respect to a chart (U1, xˆ) , it follows from the general formula (42.27) that
E =EEggi1 ⊗⋅⋅⋅⊗iiiNN =1⋅⋅⋅ gg ⊗⋅⋅⋅⊗ (45.26) ii11⋅⋅⋅ NN i i
Sec. 45 • Transformation Rules 329 where E and E ii1⋅⋅⋅ N are scalar fields on defined by ii1⋅⋅⋅ N U1
Eeg= ε (45.27) ii11⋅⋅⋅NN ii ⋅⋅⋅ and
e Eii11⋅⋅⋅NN= ε ii ⋅⋅⋅ (45.28) g
where, as in Section 42, e is + 1 if {gi ()x } is positively oriented and − 1 if {gi ()x } is negatively ⎡ ⎤ oriented, and where g is the determinant of ⎣gij ⎦ as defined by (44.45). By application of (42.28), it follows that
EggEii1⋅⋅⋅ NNN=ij11 ⋅⋅⋅ ij (45.29) jj1⋅⋅⋅ N
An interesting application of the formulas derived thus far is the derivation of an expression for the differential element of volume in curvilinear coordinates. Given the position vector r defined by (44.14) and a chart ()U1, xˆ , the differential of r can be written
i ddrxgx==i () dx (45.30)
where (44.33) has been used. Given N differentials of r , ddrr12,,,⋅⋅⋅ d rN , the differential volume element dυ generated by them is defined by
dddυ =E( rr1,, ⋅⋅⋅ N ) (45.31) or, equivalently,
ddυ =rr1 ∧⋅⋅⋅∧ dN (45.32)
330 Chap. 9 • EUCLIDEAN MANIFOLDS
12 N If we select ddxddxddxrgx11=() , rgx 2 = 2 () , ⋅⋅⋅ , rNN = gx () , we can write (45.31)) as
12 N ddxdxdxυ =E(gx1(),, ⋅⋅⋅ gN () x) ⋅⋅⋅ (45.33)
By use of (45.26) and (45.27), we then get
E()gx112(),,⋅⋅⋅ gNN () x =Eeg⋅⋅⋅ = (45.34)
Therefore,
dυ =⋅⋅⋅ gdx12 dx dx N (45.35)
For example, in the parabolic coordinates mentioned in Exercise 44.4,
22 dxxxxdxdxdxυ =+(()1212123 ()) (45.36)
Exercises
45.1 Let v be a C ∞ vector field and f be a C ∞ function both defined on U an open set in E . We define v f :UR→ by
vxvx ff()= ()grad(),⋅∈ x xU (45.37)
Show that
vvv ()λμf +=gfg λ( ) + μ( )
and
Sec. 45 • Transformation Rules 331
vv ( fgfgfg) =+( ) ( v )
For all constant functions λ , μ and all C ∞ functions f and g . In differential geometry, an operator on F ∞ (U ) with the above properties is called a derivation. Show that, conversely, every derivation on F ∞ (U ) corresponds to a unique vector field on U by (45.37). 45.2 By use of the definition (45.37), the Lie bracket of two vector fields v :UV→ and u :UV→ , written [uv, ] , is a vector field defined by
[uv, ] f =− u( vff) v( u ) (45.38)
For all scalar fields fF∈ ∞ (U ) . Show that [uv, ] is well defined by verifying that (45.38) defines a derivation on [uv, ]. Also, show that
[uv,] =−( grad vu) ( grad u) v
and then establish the following results: (a) [uv,,] =−[ vu]
(b) ⎣⎦⎣⎦⎣⎦⎡⎤⎡⎤⎡⎤vuw,,[ ] ++= uwv ,[ ,] wvu ,,[ ] 0 for all vector fields uv , and w .
ˆ ⎡⎤ (c) Let ()U1, x be a chart with natural basis field {gi} . Show that ⎣⎦ggij , = 0 . The results (a) and (b) are known as Jacobi’s identities. 45.3 In a three-dimensional Euclidean space the differential element of area normal to the plan
formed from dr1 and dr2 is defined by
dddσ = rr12×
Show that
jki dσ = Eijk dx12 dx g ()x 332 Chap. 9 • EUCLIDEAN MANIFOLDS
Section 46. Anholonomic and Physical Components of Tensors
In many applications, the components of interest are not always the components with j respect to the natural basis fields {gi } and {g }. For definiteness let us call the components of a ∞ tensor field A∈Tq ()U is defined by (45.20) the holonomic components of A . In this section, we shall consider briefly the concept of the anholonomic components of A ; i.e., the components of A taken with respect to an anholonomic basis of vector fields. The concept of the physical components of a tensor field is a special case and will also be discussed.
Let U1 be an open set in E and let {ea } denote a set of N vectors fields on U1 , which are ∞ linearly independent, i.e., at each x∈U1 , {ea } is a basis of V . If A is a tensor field in Tq ()U where UU1 ∩≠∅, then by the same type of argument as used in Section 45, we can write
a A =A eea1 ⊗⋅⋅⋅⊗ q (46.1) aa12... aq or, for example,
A =Abb1... q ee ⊗⋅⋅⋅⊗ (46.2) bb1 q
a where {e } is the reciprocal basis field to {ea } defined by
aa exex()⋅ bb ()= δ (46.3)
for all x∈U1 . Equations (46.1) and (46.2) hold on UU∩ 1 , and the component fields as defined by
A = A(ee ,..., ) (46.4) aa12... aqq a 1 a and
bb... b A 1 qq= A(eeb1 ,..., ) (46.5) Sec. 46 • Anholonomic and Physical Components 333
are scalar fields on UU∩ 1 . These fields are the anholonomic components of A when the bases a {e } and {ea } are not the natural bases of any coordinate system.
Given a set of N vector fields {ea } as above, one can show that a necessary and sufficient condition for {ea } to be the natural basis field of some chart is
[eeab, ] = 0 for all ab ,= 1,..., N, where the bracket product is defined in Exercise 45.2. We shall prove this important result in Section 49. Formulas which generalize (45.22) to anholonomic components can easily be derived. If {eˆ a } is an anholonomic basis field defined on an open set U2 such that
UU12∩≠∅, then we can express each vector field eˆ b in anholonomic component form relative to the basis {ea }, namely
a eeˆ bba= T (46.6)
a where the Tb are scalar fields on UU12∩ defined by
aa Tbb= eeˆ ⋅
The inverse of (46.6) can be written
ˆ b eeaab= T ˆ a (46.7) where
ˆ ba b TTac()xx ()= δ c (46.8) and
334 Chap. 9 • EUCLIDEAN MANIFOLDS
acˆ a TTcb()xx ()= δ b (46.9)
for all x∈∩UU12. It follows from (46.4) and (46.7) that
b ATTA=ˆˆb1 ⋅⋅⋅ q ˆ (46.10) aa11...qqq a a bb 1 ... where
Aˆ = A(eeˆˆ ,..., ) (46.11) bb11... qq b b
Equation (46.10) is the transformation rule for the anholonomic components of A . Of course,
(46.10) is a field equation which holds at every point of UUU12∩ ∩ . Similar transformation rules for the other components of A can easily be derived by the same type of argument used above.
We define the physical components of A , denoted by A , to be the anholonomic aa1 ,..., q A components of relative to the field of orthonomal basis {g i } whose basis vectors g i are unit vectors in the direction of the natural basis vectors gi of an orthogonal coordinate system. Let
()U1, xˆ by such a coordinate system with gijij = 0, ≠ . Then we define
gxi ()= gxgxii () () (no sum) (46.12)
at every x∈U1 . By (44.39), an equivalent version of (46.12) is
12 gxi ()= gxiii ()(g ()) x (no sum) (46.13)
Since {gi } is orthogonal, it follows from (46.13) and (44.39) that {g i } is orthonormal:
ggab⋅ = δab (46.14)
Sec. 46 • Anholonomic and Physical Components 335 as it should, and it follows from (44.41) that
⎡⎤1()0g11 x ⋅⋅⋅ 0 ⎢⎥01()g x 0 ⎢⎥22
ij ⎢⎥⋅⋅⋅ ⎡⎤⎣⎦g ()x = ⎢⎥ (46.15) ⎢⎥⋅⋅⋅ ⎢⎥⋅⋅⋅ ⎢⎥ ⎣⎦00⋅⋅⋅ 1()gNN x
This result shows that g i can also be written
gxi () gx()== (g ()) x12 gxi () (no sum) (46.16) i (())gii x 12 ii
Equation (46.13) can be viewed as a special case of (46.7), where
⎡⎤10g ⋅⋅⋅ 0 ⎢⎥11 ⎢⎥01g22 0 ⎢⎥⋅⋅⋅ ⎡⎤Tˆ b = ⎢⎥ (46.17) ⎣⎦a ⎢⎥⋅⋅⋅ ⎢⎥ ⎢⎥⋅⋅⋅ ⎢⎥ ⎣⎦00⋅⋅⋅ 1gNN By the transformation rule (46.10), the physical components of A are related to the covariant components of A by
−12 AggA≡=⋅⋅⋅A(gg , ,..., g ) (no sum) (46.18) aa... aaa12 a ( aa11 aqq a) a 1... a q 12 qq
Equation (46.18) is a field equation which holds for all x∈U . Since the coordinate system is orthogonal, we can replace (46.18)1 with several equivalent formulas as follows:
336 Chap. 9 • EUCLIDEAN MANIFOLDS
12 AggA= ⋅⋅⋅ aa1 ... q aa... a ()aa11 aqq a 12 q
−12 12 =⋅⋅⋅gggAaa2 ... q ()aa11() aa 2 2 aqq a a 1 ⋅ (46.19) ⋅ ⋅ −12 12 =gg ⋅⋅⋅ g Aaq ()aa11 aqq−− 1 a 1() a qq a a 1... a q − 1
In mathematical physics, tensor fields often arise naturally in component forms relative to product bases associated with several bases. For example, if {ea } and {eˆ b} are fields of bases, possibly anholonomic, then it might be convenient to express a second-order tensor field A as a field of linear transformations such that
bˆ Aeeaab==AaNˆ , 1,..., (46.20)
In this case A has naturally the component form
baˆ A = A abeeˆ ⊗ (46.21)
a a Relative to the product basis {eeˆ b ⊗ } formed by {eˆ b} and {e }, the latter being the reciprocal a basis of {ea } as usual. For definiteness, we call {eeˆ b ⊗ } a composite product basis associated a bˆ with the bases {eˆ b} and {e }. Then the scalar fields A a defined by (46.20) or (46.21), may be called the composite components of A , and they are given by
bbˆ A aa=⊗A()eeˆ (46.22)
Similarly we may define other types of composite components, e.g.,
AA==AA(,),eeˆˆbaˆ (,) ee b a (46.23) baˆ ba etc., and these components are related to one another by Sec. 46 • Anholonomic and Physical Components 337
AAgAgAgg===cccdˆˆˆˆ (46.24) baˆˆˆˆacada cbˆˆ b bc etc. Further, the composite components are related to the regular tensor components associated with a single basis field by
AATAT==ccˆ ˆ ,etc. (46.25) baˆˆˆca b bcˆ a where T a and Tˆ aˆ are given by (46.6) and (46.7) as before, In the special case where e and bˆ b { a }
{eˆ b} are orthogonal but not orthonormal, we define the normalized basis vectors e a and eˆ a as before. Then the composite physical components A of A are given by baˆ,
A = A(,)eeˆ (46.26) baˆˆ, b a
and these are related to the composite components by
12 12 AgAg= ˆ baaˆ (no sum) (46.27) baˆ, ()bbˆˆ a ( )
Clearly the concepts of composite components and composite physical components can be defined for higher order tensors also.
Exercises
46.1 In a three-dimensional Euclidean space the covariant components of a tensor field A
relative to the cylindrical coordinate system are Aij . Determine the physical components of A .
46.2 Relative to the cylindrical coordinate system the helical basis {ea } has the component form
338 Chap. 9 • EUCLIDEAN MANIFOLDS
eg1 = 1
eg2 =+()cosαα23() sin g (46.28)
e3 =−()sinααgg23 +() cos
where α is a constant called the pitch, and where {g α } is the orthonormal basis associated
with the natural basis of the cylindrical system. Show that {ea } is anholonomic. Determine the anholonomic components of the tensor field A in the preceding exercise relative to the helical basis. 46.3 Determine the composite physical components of A relative to the composite product basis
{e ab⊗ g }
Sec. 47 • Christoffel Symbols, Covariant Differentiation 339
Section 47. Christoffel Symbols and Covariant Differentiation
In this section we shall investigate the problem of representing the gradient of various tensor fields in components relative to the natural basis of arbitrary coordinate systems. We consider first the simple case of representing the tangent of a smooth curve in E . Let λ :,()ab →E be a smooth curve passing through a point x , say x = λ (c ) . Then the tangent vector of λ at x is defined by
dtλ() λ = (47.1) x dt tc=
Given the chart ()U , xˆ covering x , we can project the vector equation (47.1) into the natural basis
{gi} of xˆ . First, the coordinates of the curve λ are given by
xˆ()λλ(ttt )==()λλ11 ( ),...,NN ( ) ,( t) x ( λλ ( tt ),..., ( )) (47.2)
for all t such that λ ()t ∈U . Differentiating (47.2)2 with respect to t , we get
∂x dλ j λ = (47.3) x j ∂x dt tc=
By (47.3) this equation can be rewritten as
dλ j λ = g ()x (47.4) x j dt tc=