CHM 532 Notes on Angular Momentum Eigenvalues and Eigenfunctions

Total Page:16

File Type:pdf, Size:1020Kb

CHM 532 Notes on Angular Momentum Eigenvalues and Eigenfunctions CHM 532 Notes on Angular Momentum Eigenvalues and Eigenfunctions In your textbooks, the eigenfunctions and eigenvalues of the angular momentum operators are determined using differential equations methods. A more powerful approach is to solve the angular momentum eigenfunction problem using operator methods analogous to the creation and annihilation operators we have used in the harmonic oscillator problem. These notes provide details about the operator approach. 1 Preliminaries We begin by reviewing the angular momentum operators and their commutation relations. The detailed derivation of these preliminary results can be found in your textbooks. Recall that the angular momentum is a vector with three components, and there is a Hermitian operator associated with each component h¯ ∂ ∂ ! Lˆ = y − z , (1) x i ∂z ∂y h¯ ∂ ∂ ! Lˆ = z − x , (2) y i ∂x ∂z and h¯ ∂ ∂ ! Lˆ = x − y . (3) z i ∂y ∂x Additionally, it is convenient to define the operator corresponding to the square of the length of the angular momentum vector ˆ2 ˆ2 ˆ2 ˆ2 L = Lx + Ly + Lz. (4) Because the angular momentum is a conserved quantity for systems having spherically sym- metric potentials, it is important to express these operators as well in spherical polar coor- dinates h¯ ∂ ∂ ! Lˆ = − sin φ − cot θ cos φ , (5) x i ∂θ ∂φ 1 h¯ ∂ ∂ ! Lˆ = cos φ − cot θ sin φ , (6) y i ∂θ ∂φ h¯ ∂ Lˆ = , (7) z i ∂φ and 1 ∂ ∂ 1 ∂2 ! Lˆ2 = −h¯2 sin θ + . (8) sin θ ∂θ ∂θ sin2 θ ∂φ2 It is important to recognize that in spherical polar coordinates, all the angular momentum operators are independent of the coordinate r, and the eigenfunctions must be functions only of the angular coordinates; i.e. θ and φ. Using the operator expressions in either coordinate system (it is most easy to use Carte- sian coordinates), it is possible to demonstrate the following commutators (see the text for details) ˆ ˆ ˆ [Lx, Ly] = ih¯Lz, (9) ˆ ˆ ˆ [Ly, Lz] = ih¯Lx, (10) ˆ ˆ ˆ [Lz, Lx] = ih¯Ly, (11) and ˆ ˆ2 [Li, L ] = 0 (12) for i = x, y or z. From the commutation relations, it is clear that we can know only one component of the angular momentum and the square of the angular momentum vector at ˆ the same time. Because Lz has a particularly simple form in spherical polar coordinates, we 2 ˆ choose to know Lz and L simultaneously; i.e. we seek simultaneous eigenfunctions of Lz and Lˆ2. 2 The Eigenvalues In the case of the harmonic oscillator, we discovered the eigenvalues of the Hamiltonian by introducing creation and annihilation operators. In analogy, we make the following defini- tions: Definition 1 The step-up operator is defined by ˆ ˆ ˆ L+ = Lx + iLy. (13) Definition 2 The step-down operator is defined as the adjoint of the step-up operator; i.e. ˆ ˆ† ˆ ˆ L− = L+ = Lx − iLy. (14) 2 We next introduce and prove a series of lemmas from which we can extract the eigenvalues of ˆ ˆ2 Lz and L . We let {Yl,m} represent the common complete orthonormal set of eigenfunctions ˆ ˆ2 of Lz and L with m and l respectively the quantum numbers associated with each operator. We write ˆ LzYl,m = mhY¯ l,m (15) and ˆ2 2 L Yl,m = f(l)¯h Yl,m (16) where f(l) is some function of the l quantum number the explicit expression for which is derived below. It is easily verified that the units ofh ¯ are units of angular momentum, and we have included factors ofh ¯ explicitly in Eqs.(15) and (16). The orthonormality condition is written hYl,m|Yl0,m0 i = δl,l0 δm,m0 (17) We now introduce and provide proofs for the necessary lemmas. Lemma 1 ˆ ˆ ˆ ˆ LzL+ = L+(Lz +h ¯). (18) Proof: ˆ ˆ ˆ ˆ ˆ LzL+ = Lz(Lx + iLy) (19) ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ = LxLz + [Lz, Lx] + i(LyLz + [Lz, Ly]) (20) ˆ ˆ ˆ ˆ ˆ ˆ = LxLz + ih¯Ly + i(LyLz − ih¯Lx) (21) ˆ ˆ ˆ ˆ ˆ ˆ ˆ = (Lx + iLy)Lz +h ¯(Lx + iLy) = L+(Lz +h ¯). (22) 2 Lemma 2 ˆ ˆ ˆ ˆ LzL− = L+(Lz − h¯). (23) Proof: ˆ ˆ ˆ ˆ ˆ LzL− = Lz(Lx − iLy) (24) ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ = LxLz + [Lz, Lx] − i(LyLz + [Lz, Ly]) (25) ˆ ˆ ˆ ˆ ˆ ˆ = LxLz + ih¯Ly − i(LyLz − ih¯Lx) (26) ˆ ˆ ˆ ˆ ˆ ˆ ˆ = (Lx − iLy)Lz − h¯(Lx − iLy) = L−(Lz − h¯). (27) 2 Lemma 3 ˆ L+Yl,m = c+(l, m)Yl,m+1 (28) where c+(l, m) is a function only of the quantum numbers l and m. 3 ˆ ˆ2 ˆ Proof: We must show that L+Yl,m is an eigenfunction of L and Lz with the value of the z ˆ ˆ2 quantum number increased by one unit ofh ¯. Notice that L+ commutes with L so that the ˆ l quantum number must be unaffected by the action of L+. Next ˆ ˆ ˆ ˆ Lz(L+Yl,m) = L+(Lz +h ¯)Yl,m (29) ˆ ˆ = L+(m + 1)¯hYl,m = (m + 1)¯h(L+Yl,m). (30) 2 Lemma 4 ˆ L−Yl,m = c−(l, m)Yl,m−1 (31) where c−(l, m) is a function only of the quantum numbers l and m. ˆ ˆ2 ˆ Proof: We must show that L−Yl,m is an eigenfunction of L and Lz with the value of the z ˆ ˆ2 quantum number decreased by one unit ofh ¯. Notice that L− commutes with L so that the ˆ l quantum number must be unaffected by the action of L−. Next ˆ ˆ ˆ ˆ Lz(L−Yl,m) = L−(Lz − h¯)Yl,m (32) ˆ ˆ = L−(m − 1)¯hYl,m = (m − 1)¯h(L−Yl,m). (33) 2 Lemma 5 f(l) ≥ m2. (34) Proof: We begin with ˆ ˆ ˆ ˆ hYl,m|L−L+|Yl,mi = hL+Yl,m|L+Yl,mi ≥ 0. (35) Next, ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 ˆ2 ˆ ˆ ˆ ˆ L−L+ = (Lx − iLy)(Lx + iLy) = Lx + Ly + i(LxLy − LyLx) (36) ˆ2 ˆ2 ˆ ˆ2 ˆ2 ˆ = Lx + Ly − h¯Lz = L − Lz − h¯Lz. (37) Then ˆ ˆ 2 2 2 2 2 hYl,m|L−L+|Yl,mi = f(l)¯h − m h¯ − h¯ m = [f(l) − m(m + 1)]¯h ≥ 0 (38) or f(l) − m(m + 1) ≥ 0. (39) Next, ˆ ˆ ˆ ˆ hYl,m|L+L−|Yl,mi = hL−Yl,m|L−Yl,mi ≥ 0. (40) Also, ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 ˆ2 ˆ ˆ ˆ ˆ L+L− = (Lx + iLy)(Lx − iLy) = Lx + Ly − i(LxLy − LyLx) (41) ˆ2 ˆ2 ˆ ˆ2 ˆ2 ˆ = Lx + Ly +h ¯Lz = L − Lz +h ¯Lz. (42) 4 Then ˆ ˆ 2 2 2 2 2 hYl,m|L+L−|Yl,mi = f(l)¯h − m h¯ +h ¯ m = [f(l) − m(m − 1)]¯h ≥ 0 (43) or f(l) − m(m − 1) ≥ 0. (44) Adding Eqs.(39) and (44), we obtain f(l) − m2 ≥ 0 (45) or f(l) ≥ m2. (46) 2 An implication of Eq.(46) is that for a given l, there exits a maximum value of m which we denote m and a minimum value of m which we denote m. The next lemma relates m to m. Lemma 6 m = −m (47) Proof: ˆ ˆ L+Yl,m = L−Yl,m = 0. (48) ˆ ˆ ˆ2 ˆ2 ˆ 2 hYl,m|L−L+|Yl,mi = hYl,m|L − Lz − h¯Lz|Yl,mi = [f(l) − m(m + 1)]¯h = 0. (49) ˆ ˆ ˆ2 ˆ2 ˆ 2 hYl,m|L+L−|Yl,mi = hYl,m|L − Lz +h ¯Lz|Yl,mi = [f(l) − m(m − 1)]¯h = 0. (50) Then f(l) = m(m + 1) = m(m − 1). (51) Solving for m in terms of m, we obtain m = −m. (52) 2 From Eq.(49), we have f(l) = m(m + 1). (53) We define l = m; i.e. the maximum value of m for a given l, and we can write f(l) = l(l + 1). (54) In the next lemma, we identify what values l can have. Lemma 7 The quantum number l can be an integer or a 1/2 integer. 5 Proof: We have the following series of equations ˆ L−Yl,l = c−(l, l)Yl,l−1 (55) ˆ L−Yl,l−1 = c−(l, l − 1)Yl,l−2 (56) . ˆ L−Y (l, −l + 1) = c−(l, −l + 1)Yl,−l. (57) We connect m = l to m = −l by an integer number k of operations. Then l − (−l) = 2l = k (58) or k l = (59) 2 an integer or a half integer. 2 We summarize the results of these lemmas with the following set of equations and statements: ˆ LzYl,m = mhY¯ l,m (60) ˆ2 2 L Yl,m = l(l + 1)¯h Yl,m (61) −l ≤ m ≤ l (62) where m and l can either be integers or half integers. Notice for a given l, there are gl = 2l+1 different m-states. It is interesting that the algebraic commutation relations for the angular momentum operators allows both integer and half integer quantum numbers for the angular momentum. When we solve for the eigenfunctions in a subsequent section of these notes, we find that the half integer quantum numbers are excluded.
Recommended publications
  • Quaternion Zernike Spherical Polynomials
    MATHEMATICS OF COMPUTATION Volume 84, Number 293, May 2015, Pages 1317–1337 S 0025-5718(2014)02888-3 Article electronically published on August 29, 2014 QUATERNION ZERNIKE SPHERICAL POLYNOMIALS J. MORAIS AND I. CAC¸ AO˜ Abstract. Over the past few years considerable attention has been given to the role played by the Zernike polynomials (ZPs) in many different fields of geometrical optics, optical engineering, and astronomy. The ZPs and their applications to corneal surface modeling played a key role in this develop- ment. These polynomials are a complete set of orthogonal functions over the unit circle and are commonly used to describe balanced aberrations. In the present paper we introduce the Zernike spherical polynomials within quater- nionic analysis ((R)QZSPs), which refine and extend the Zernike moments (defined through their polynomial counterparts). In particular, the underlying functions are of three real variables and take on either values in the reduced and full quaternions (identified, respectively, with R3 and R4). (R)QZSPs are orthonormal in the unit ball. The representation of these functions in terms of spherical monogenics over the unit sphere are explicitly given, from which several recurrence formulae for fast computer implementations can be derived. A summary of their fundamental properties and a further second or- der homogeneous differential equation are also discussed. As an application, we provide the reader with plot simulations that demonstrate the effectiveness of our approach. (R)QZSPs are new in literature and have some consequences that are now under investigation. 1. Introduction 1.1. The Zernike spherical polynomials. The complex Zernike polynomials (ZPs) have long been successfully used in many different fields of optics.
    [Show full text]
  • Vectors, Matrices and Coordinate Transformations
    S. Widnall 16.07 Dynamics Fall 2009 Lecture notes based on J. Peraire Version 2.0 Lecture L3 - Vectors, Matrices and Coordinate Transformations By using vectors and defining appropriate operations between them, physical laws can often be written in a simple form. Since we will making extensive use of vectors in Dynamics, we will summarize some of their important properties. Vectors For our purposes we will think of a vector as a mathematical representation of a physical entity which has both magnitude and direction in a 3D space. Examples of physical vectors are forces, moments, and velocities. Geometrically, a vector can be represented as arrows. The length of the arrow represents its magnitude. Unless indicated otherwise, we shall assume that parallel translation does not change a vector, and we shall call the vectors satisfying this property, free vectors. Thus, two vectors are equal if and only if they are parallel, point in the same direction, and have equal length. Vectors are usually typed in boldface and scalar quantities appear in lightface italic type, e.g. the vector quantity A has magnitude, or modulus, A = |A|. In handwritten text, vectors are often expressed using the −→ arrow, or underbar notation, e.g. A , A. Vector Algebra Here, we introduce a few useful operations which are defined for free vectors. Multiplication by a scalar If we multiply a vector A by a scalar α, the result is a vector B = αA, which has magnitude B = |α|A. The vector B, is parallel to A and points in the same direction if α > 0.
    [Show full text]
  • Laplacian Eigenmodes for the Three Sphere
    Laplacian eigenmodes for the three sphere M. Lachi`eze-Rey Service d’Astrophysique, C. E. Saclay 91191 Gif sur Yvette cedex, France October 10, 2018 Abstract The vector space k of the eigenfunctions of the Laplacian on the 3 V three sphere S , corresponding to the same eigenvalue λk = k (k + 2 − k 2), has dimension (k + 1) . After recalling the standard bases for , V we introduce a new basis B3, constructed from the reductions to S3 of peculiar homogeneous harmonic polynomia involving null vectors. We give the transformation laws between this basis and the usual hyper-spherical harmonics. Thanks to the quaternionic representations of S3 and SO(4), we are able to write explicitely the transformation properties of B3, and thus of any eigenmode, under an arbitrary rotation of SO(4). This offers the possibility to select those functions of k which remain invariant under V a chosen rotation of SO(4). When the rotation is an holonomy transfor- mation of a spherical space S3/Γ, this gives a method to calculates the eigenmodes of S3/Γ, which remains an open probleme in general. We illustrate our method by (re-)deriving the eigenmodes of lens and prism space. In a forthcoming paper, we present the derivation for dodecahedral space. 1 Introduction 3 The eigenvalues of the Laplacian ∆ of S are of the form λk = k (k + 2), + − k where k IN . For a given value of k, they span the eigenspace of ∈ 2 2 V dimension (k +1) . This vector space constitutes the (k +1) dimensional irreductible representation of SO(4), the isometry group of S3.
    [Show full text]
  • Chapter 5 ANGULAR MOMENTUM and ROTATIONS
    Chapter 5 ANGULAR MOMENTUM AND ROTATIONS In classical mechanics the total angular momentum L~ of an isolated system about any …xed point is conserved. The existence of a conserved vector L~ associated with such a system is itself a consequence of the fact that the associated Hamiltonian (or Lagrangian) is invariant under rotations, i.e., if the coordinates and momenta of the entire system are rotated “rigidly” about some point, the energy of the system is unchanged and, more importantly, is the same function of the dynamical variables as it was before the rotation. Such a circumstance would not apply, e.g., to a system lying in an externally imposed gravitational …eld pointing in some speci…c direction. Thus, the invariance of an isolated system under rotations ultimately arises from the fact that, in the absence of external …elds of this sort, space is isotropic; it behaves the same way in all directions. Not surprisingly, therefore, in quantum mechanics the individual Cartesian com- ponents Li of the total angular momentum operator L~ of an isolated system are also constants of the motion. The di¤erent components of L~ are not, however, compatible quantum observables. Indeed, as we will see the operators representing the components of angular momentum along di¤erent directions do not generally commute with one an- other. Thus, the vector operator L~ is not, strictly speaking, an observable, since it does not have a complete basis of eigenstates (which would have to be simultaneous eigenstates of all of its non-commuting components). This lack of commutivity often seems, at …rst encounter, as somewhat of a nuisance but, in fact, it intimately re‡ects the underlying structure of the three dimensional space in which we are immersed, and has its source in the fact that rotations in three dimensions about di¤erent axes do not commute with one another.
    [Show full text]
  • Solving the Geodesic Equation
    Solving the Geodesic Equation Jeremy Atkins December 12, 2018 Abstract We find the general form of the geodesic equation and discuss the closed form relation to find Christoffel symbols. We then show how to use metric independence to find Killing vector fields, which allow us to solve the geodesic equation when there are helpful symmetries. We also discuss a more general way to find Killing vector fields, and some of their properties as a Lie algebra. 1 The Variational Method We will exploit the following variational principle to characterize motion in general relativity: The world line of a free test particle between two timelike separated points extremizes the proper time between them. where a test particle is one that is not a significant source of spacetime cur- vature, and a free particles is one that is only under the influence of curved spacetime. Similarly to classical Lagrangian mechanics, we can use this to de- duce the equations of motion for a metric. The proper time along a timeline worldline between point A and point B for the metric gµν is given by Z B Z B µ ν 1=2 τAB = dτ = (−gµν (x)dx dx ) (1) A A using the Einstein summation notation, and µ, ν = 0; 1; 2; 3. We can parame- terize the four coordinates with the parameter σ where σ = 0 at A and σ = 1 at B. This gives us the following equation for the proper time: Z 1 dxµ dxν 1=2 τAB = dσ −gµν (x) (2) 0 dσ dσ We can treat the integrand as a Lagrangian, dxµ dxν 1=2 L = −gµν (x) (3) dσ dσ and it's clear that the world lines extremizing proper time are those that satisfy the Euler-Lagrange equation: @L d @L − = 0 (4) @xµ dσ @(dxµ/dσ) 1 These four equations together give the equation for the worldline extremizing the proper time.
    [Show full text]
  • A User's Guide to Spherical Harmonics
    A User's Guide to Spherical Harmonics Martin J. Mohlenkamp∗ Version: October 18, 2016 This pamphlet is intended for the scientist who is considering using Spherical Harmonics for some appli- cation. It is designed to introduce the Spherical Harmonics from a theoretical perspective and then discuss those practical issues necessary for their use in applications. I expect great variability in the backgrounds of the readers. In order to be informative to those who know less, without insulting those who know more, I have adopted the strategy \State even basic facts, but briefly." My dissertation work was a fast transform for Spherical Harmonics [6]. After its publication, I started receiving questions about Spherical Harmonics in general, rather than about my work specifically. This pamphlet aims to answer such general questions. If you find mistakes, or feel that important material is unclear or missing, please inform me. 1 A Theory of Spherical Harmonics In this section we give a development of Spherical Harmonics. There are other developments from other perspectives. The one chosen here has the benefit of being very concrete. 1.1 Mathematical Preliminaries We define the L2 inner product of two functions to be Z hf; gi = f(s)¯g(s)ds (1) R 2π whereg ¯ denotes complex conjugation and the integral is over the space of interest, for example 0 dθ on the R 2π R π 2 p 2 circle or 0 0 sin φdφdθ on the sphere. We define the L norm by jjfjj = hf; fi. The space L consists of all functions such that jjfjj < 1.
    [Show full text]
  • Geodetic Position Computations
    GEODETIC POSITION COMPUTATIONS E. J. KRAKIWSKY D. B. THOMSON February 1974 TECHNICALLECTURE NOTES REPORT NO.NO. 21739 PREFACE In order to make our extensive series of lecture notes more readily available, we have scanned the old master copies and produced electronic versions in Portable Document Format. The quality of the images varies depending on the quality of the originals. The images have not been converted to searchable text. GEODETIC POSITION COMPUTATIONS E.J. Krakiwsky D.B. Thomson Department of Geodesy and Geomatics Engineering University of New Brunswick P.O. Box 4400 Fredericton. N .B. Canada E3B5A3 February 197 4 Latest Reprinting December 1995 PREFACE The purpose of these notes is to give the theory and use of some methods of computing the geodetic positions of points on a reference ellipsoid and on the terrain. Justification for the first three sections o{ these lecture notes, which are concerned with the classical problem of "cCDputation of geodetic positions on the surface of an ellipsoid" is not easy to come by. It can onl.y be stated that the attempt has been to produce a self contained package , cont8.i.ning the complete development of same representative methods that exist in the literature. The last section is an introduction to three dimensional computation methods , and is offered as an alternative to the classical approach. Several problems, and their respective solutions, are presented. The approach t~en herein is to perform complete derivations, thus stqing awrq f'rcm the practice of giving a list of for11111lae to use in the solution of' a problem.
    [Show full text]
  • Rotation Matrix - Wikipedia, the Free Encyclopedia Page 1 of 22
    Rotation matrix - Wikipedia, the free encyclopedia Page 1 of 22 Rotation matrix From Wikipedia, the free encyclopedia In linear algebra, a rotation matrix is a matrix that is used to perform a rotation in Euclidean space. For example the matrix rotates points in the xy -Cartesian plane counterclockwise through an angle θ about the origin of the Cartesian coordinate system. To perform the rotation, the position of each point must be represented by a column vector v, containing the coordinates of the point. A rotated vector is obtained by using the matrix multiplication Rv (see below for details). In two and three dimensions, rotation matrices are among the simplest algebraic descriptions of rotations, and are used extensively for computations in geometry, physics, and computer graphics. Though most applications involve rotations in two or three dimensions, rotation matrices can be defined for n-dimensional space. Rotation matrices are always square, with real entries. Algebraically, a rotation matrix in n-dimensions is a n × n special orthogonal matrix, i.e. an orthogonal matrix whose determinant is 1: . The set of all rotation matrices forms a group, known as the rotation group or the special orthogonal group. It is a subset of the orthogonal group, which includes reflections and consists of all orthogonal matrices with determinant 1 or -1, and of the special linear group, which includes all volume-preserving transformations and consists of matrices with determinant 1. Contents 1 Rotations in two dimensions 1.1 Non-standard orientation
    [Show full text]
  • Coordinate Transformation
    Coordinate Transformation Coordinate Transformations In this chapter, we explore mappings – where a mapping is a function that "maps" one set to another, usually in a way that preserves at least some of the underlyign geometry of the sets. For example, a 2-dimensional coordinate transformation is a mapping of the form T (u; v) = x (u; v) ; y (u; v) h i The functions x (u; v) and y (u; v) are called the components of the transforma- tion. Moreover, the transformation T maps a set S in the uv-plane to a set T (S) in the xy-plane: If S is a region, then we use the components x = f (u; v) and y = g (u; v) to …nd the image of S under T (u; v) : EXAMPLE 1 Find T (S) when T (u; v) = uv; u2 v2 and S is the unit square in the uv-plane (i.e., S = [0; 1] [0; 1]). Solution: To do so, let’s determine the boundary of T (S) in the xy-plane. We use x = uv and y = u2 v2 to …nd the image of the lines bounding the unit square: Side of Square Result of T (u; v) Image in xy-plane v = 0; u in [0; 1] x = 0; y = u2; u in [0; 1] y-axis for 0 y 1 u = 1; v in [0; 1] x = v; y = 1 v2; v in [0; 1] y = 1 x2; x in[0; 1] v = 1; u in [0; 1] x = u; y = u2 1; u in [0; 1] y = x2 1; x in [0; 1] u = 0; u in [0; 1] x = 0; y = v2; v in [0; 1] y-axis for 1 y 0 1 As a result, T (S) is the region in the xy-plane bounded by x = 0; y = x2 1; and y = 1 x2: Linear transformations are coordinate transformations of the form T (u; v) = au + bv; cu + dv h i where a; b; c; and d are constants.
    [Show full text]
  • Unit 5: Change of Coordinates
    LINEAR ALGEBRA AND VECTOR ANALYSIS MATH 22B Unit 5: Change of Coordinates Lecture 5.1. Given a basis B in a linear space X, we can write an element v in X in a unique 3 way as a sum of basis elements. For example, if v = is a vector in X = 2 and 4 R 1 1 2 B = fv = ; v = g, then v = 2v + v . We say that are the B 1 −1 2 6 1 2 1 B 3 coordinates of v. The standard coordinates are v = are assumed if no other 4 basis is specified. This means v = 3e1 + 4e2. n 5.2. If B = fv1; v2; ··· ; vng is a basis of R , then the matrix S which contains the vectors vk as column vectors is called the coordinate change matrix. Theorem: If S is the matrix of B, then S−1v are the B coordinates of v. 1 1 6 −1 5.3. In the above example, S = has the inverse S−1 = =7. We −1 6 1 1 compute S−1[3; 4]T = [2; 1]T . Proof. If [v]B = [a1; : : : ; an] are the new coordinates of v, this means v = a1v1 + ··· + −1 anvn. But that means v = S[v]B. Since B is a basis, S is invertible and [v]B = S v. Theorem: If T (x) = Ax is a linear map and S is the matrix from a basis change, then B = S−1AS is the matrix of T in the new basis B. Proof. Let y = Ax. The statement [y]B = B[x]B can be written using the last theorem as S−1y = BS−1x so that y = SBS−1x.
    [Show full text]
  • COORDINATE TRANSFORMATIONS Members of a Structural System Are Typically Oriented in Differing Directions, E.G., Fig
    COORDINATE TRANSFORMATIONS Members of a structural system are typically oriented in differing directions, e.g., Fig. 17.1. In order to perform an analysis, the ele- ment stiffness equations need to be expressed in a common coor- dinate system – typically the global coordinate system. Once the element equations are expressed in a common coordinate system, the equations for each element comprising the structure can be Figure 17.1 – Frame Structure 1 (Kassimali, 1999) 2 assembled. Coordinate Transformations: Frame Elements Consider frame element m of Fig. 17.7. Y e x y m b X (a) Frame Q , u 6 6 Q , u m 4 4 Q3, u3 Q , u Q5, u5 1 1 Figure 17.7: Local – Global Q2, u2 Coordinate Relationships (b) Local Coordinate End Forces and Displacements 3 4 1 Figure 17.7 shows that the local – QcosFsinFxXY global coordinate transformations QsinFcosFyXY (17.9) can be expressed as QFzZ x = cos X + sin Y y = -sin X + cos Y where x, X = 1 or 4; y, Y = 2 or 5; and z, Z = 3 or 6. and since z and Z are parallel, this coordinate transformation is Utilizing (17.9) for all six member expressed as force components and expressing z = Z the resulting transformations in Using the above coordinate matrix form gives transformations, the end force and QFb [t] [0] b displacement transformations can QFee[0] [t] be expressed as or 5 {Q} = [T] {F} (17.11) T {Q} Q Q Q T where b 123 = member; {Q} = <<Q>b <Q>e> = beginning node local coordinate element local coordinate force T T force vector; {Q}e = <Q4 Q5 Q6> = vector; {F} = <<F>b <F>e> = end node local coordinate force element global coordinate force T vector; {F}b = <F1 F2 F3> = [t] [0] beginning node global coordinate vector; and [T] = = T [0] [t] force vector; {F}e = <F4 F5 F6> = end node global coordinate force element local to global coordinate cos sin 0 transformation matrix.
    [Show full text]
  • 1 Vectors & Tensors
    1 Vectors & Tensors The mathematical modeling of the physical world requires knowledge of quite a few different mathematics subjects, such as Calculus, Differential Equations and Linear Algebra. These topics are usually encountered in fundamental mathematics courses. However, in a more thorough and in-depth treatment of mechanics, it is essential to describe the physical world using the concept of the tensor, and so we begin this book with a comprehensive chapter on the tensor. The chapter is divided into three parts. The first part covers vectors (§1.1-1.7). The second part is concerned with second, and higher-order, tensors (§1.8-1.15). The second part covers much of the same ground as done in the first part, mainly generalizing the vector concepts and expressions to tensors. The final part (§1.16-1.19) (not required in the vast majority of applications) is concerned with generalizing the earlier work to curvilinear coordinate systems. The first part comprises basic vector algebra, such as the dot product and the cross product; the mathematics of how the components of a vector transform between different coordinate systems; the symbolic, index and matrix notations for vectors; the differentiation of vectors, including the gradient, the divergence and the curl; the integration of vectors, including line, double, surface and volume integrals, and the integral theorems. The second part comprises the definition of the tensor (and a re-definition of the vector); dyads and dyadics; the manipulation of tensors; properties of tensors, such as the trace, transpose, norm, determinant and principal values; special tensors, such as the spherical, identity and orthogonal tensors; the transformation of tensor components between different coordinate systems; the calculus of tensors, including the gradient of vectors and higher order tensors and the divergence of higher order tensors and special fourth order tensors.
    [Show full text]