Normal Mode Analysis ©David Ronis Mcgill University

Total Page:16

File Type:pdf, Size:1020Kb

Normal Mode Analysis ©David Ronis Mcgill University Chemistry 365: Normal Mode Analysis ©David Ronis McGill University 1. Quantum Mechanical Treatment Our starting point is the Schrodinger wav e equation: N −2 ∂2 − h + → → Ψ → → = Ψ → → Σ → U(r1,...,r N ) (r1,...,r N ) E (r1,...,r N ), (1.1) = ∂ 2 i 1 2mi ri where N is the number of atoms in the molecule, mi is the mass of the i’th atom, and → → U(r1,...,r N )isthe effective potential for the nuclear motion, e.g., as is obtained in the Born- Oppenheimer approximation. If the amplitude of the vibrational motion is small, then the vibrational part of the Hamil- tonian associated with Eq. (1.1) can be written as: N −2 ∂2 N ↔ → → ≈− h + + 1 ∆ ∆ Hvib Σ →2 U0 Σ Ki, j: i j,(1.2) i=1 2mi ∂∆ 2 i, j=1 i → ∆ ≡ → − → → where U0 is the minimum value of the potential energy, i ri Ri, Ri is the equilibrium posi- tion of the i’th atom, and 2 ↔ ∂ ≡ U Ki, j → → (1.3) ∂ ∂ → ri r j → r k = Rk is the matrix of (harmonic) force constants. Henceforth, we will shift the zero of energy so as to = make U0 0. Note that in obtaining Eq. (1.2), we have neglected anharmonic (i.e., cubic and higher order) corrections to the vibrational motion. The next and most confusing step is to change to matrix notation. We introduce a column vector containing the displacements as: ∆≡ ∆x ∆y ∆z ∆x ∆y ∆z T [ 1 , 1, 1,..., N , N , N ] ,(1.4) where "T"denotes a matrix transpose. Similarly,wecan encode the force constants or masses into 3N × 3N matrices, and thereby rewrite Eq. (1.2) as −2 † h ∂ − ∂ 1 H =− M 1 + ∆†K∆,(1.5) vib 2 ∂∆ ∂∆ 2 where all matrix quantities are emboldened, the superscript "†" denotes the Hermitian conjugate (matrix transpose for real matrices), and where the mass matrix, M,isadiagonal matrix with the masses of the givenatoms each appearing three times on the diagonal. The Hamiltonian givenbyEq. (1.5) is the generalization of the usual harmonic oscillator Hamiltonian to include more particles and to allowfor "springs" between arbitrary particles. Unless K is diagonal (and it usually isn’t) Eq. (1.5) would seem to suggest that the vibrational problem for a polyatomic is non-separable (why?); nonetheless, as we nowshow, itcan be sepa- rated. This is first shown using using the quantum mechanical framework we’ve just set up. Winter Term, 2019 Chemistry 365 -2- Normal Mode Analysis Later,wewill discuss the separation using classical mechanics. This is valid since, for harmonic oscillators, you get the same result for the vibration frequencies. → → ∆ → −1/2∆ We now makethe transformation i mi i,which allows Eq. (1.5) to be reexpressed in the newcoordinates as −h 2 ∂ † ∂ 1 H =− + ∆†K˜ ∆,(1.6) vib 2 ∂∆ ∂∆ 2 where − − K˜ ≡ M 1/2KM 1/2.(1.7) Since the matrix K˜ is Hermitian (or symmetric for real matrices), it is possible to find a uni- tary (also referred to as an orthogonal matrix for real matrices) matrix which diagonalizes it; i.e., you can find a matrix P which satisfies P†KP˜ = λ or K˜ = PλP†,(1.8) where λ is diagonal (with real eigenvalues) and where PP† = 1(1.9) where 1 is the identity matrix. Note that P is the unitary matrix who’scolumns are the normal- ized eigenvectors. (See agood quantum mechanics book or anylinear algebra text for proofs of these results). We now makethe transformation ∆→P∆ (1.10) in Eq. (1.5). This gives −h 2 ∂ † ∂ 1 H =− P†P + ∆†λ˜ ∆,(1.11) vib 2 ∂∆ ∂∆ 2 Finally,Eq. (1.8) allows us to cancel the factors of P in this last equation and write 3N 1 ∂2 1 H = Σ − −h 2 + λ ∆2 .(1.12) vib ∂∆2 i i i=1 2 i 2 Thus the twotransformations described above hav e separated the Hamiltonian into 3N uncoupled harmonic oscillator Hamiltonians, and are usually referred to as a normal mode transformation. Remember that in general the normal mode coordinates are not the original displacements from the equilibrium positions, but correspond to collective vibrations of the molecule. 2. Normal Modes in Classical Mechanics The starting point for the normal mode analysis in classical mechanics are Newton’sequa- tions for a system of coupled harmonic oscillators. Still, using the notation that led to Eq. (1.12), we can write the classical equations of motion as: M∆¨ (t) =−K∆(t), (2.1) where components of the left-hand-side of the equation are the rates of change of momentum of the nuclei, while the right-hand-side contains the harmonic forces. Winter Term, 2019 Chemistry 365 -3- Normal Mode Analysis Since we expect harmonic motion, we’ll look for a solution of the form: ∆(t) = cos(ω t + φ)Y,(2.2) where φ is an arbitrary phase shift. If Eq. (2.2) is used in Eq. (2.1) it follows that: (Mω 2 − K)Y = 0, (2.3) where remember that Y is a column vector and M,and K are matrices. We want to solvethis homogeneous system of linear equations for Y.Ingeneral, the only way to get a nonzero solu- tion is to makethe matrix [Mω 2 − K]singular; i.e., we must set − − det(Mω 2 − K) = det(ω 2 − M 1/2KM 1/2)det(M) = 0. (2.4) = 3 ≠ Since, det(M) (M1 M2...M N ) 0wesee that the second equality is simply the characteristic equation associated with the eigenvalue problem: Ku˜ = λu,(2.5) with λ ≡ ω 2 and with K˜ givenbyEq. (2.7) above.The remaining steps are equivalent to what wasdone quantum mechanically. 3. Force Constant Calculations Here is an example of a force constant matrix calculation. We will consider a diatomic molecule, where the twoatoms interact with a potential of the form: 1 2 U(r , r ) ≡ K |r − r | − R ;(3.1) 1 2 2 1 2 0 i.e., a simple Hookian spring. It is easy to takethe various derivativesindicated in the preceding sections; here, however, wewill explicitly expand the potential in terms of the atomic displace- ments, ∆i.Bywriting, ri = Ri +∆i (where Ri is the equilibrium position of the i’th nucleus), Eq., (3.1) can be rewritten as: 2 1 1/2 = 2 + ⋅ ∆ −∆ + ∆ −∆ 2 − U(r1, r2) K R12 2R12 ( 1 2) | 1 2| R0 ,(3.2) 2 ≡ − = where R12 R1 R2.Clearly,the equilibrium will have |R12| R0.Moreover, weexpect that the vibrational amplitude will be small, and thus, the terms in the ∆’s inthe square root in Eq. (3.2) will be small compared with the first term. The Taylor expansion of the square root implies that B √ √ A+ B ≈ √ A + +..., (3.3) 2 √ A we can write ⋅ ∆ −∆ + ∆ −∆ 2 2 = 1 + 2R12 ( 1 2) | 1 2| − U(r1, r2) KR12 R0 ,(3.4) 2 2R12 which can be rewritten as 1 U(r , r ) = K[Rˆ ⋅ (∆ −∆ )]2,(3.5) 1 2 2 12 1 2 Winter Term, 2019 Chemistry 365 -4- Normal Mode Analysis where the ˆ denotes a unit vector and where all terms smaller than quadratic in the nuclear dis- placements have been dropped. If the squareisexpanded, notice the appearance of cross terms in the displacements of 1 and 2. It is actually quite simple to finish the normal mode calculation in this case. To doso, define the equilibrium bond to point in along the x axis. Equation (3.5) shows that only x-com- ponents of the displacements cost energy,and hence, there will no force in thy y or z directions (thereby resulting in 4 zero eigen-frequencies). Forthe x components, Newton’sequations become: ∆¨x − ∆x m1 0 ⋅ 1 =− K K ⋅ 1 ∆x − ∆x,(3.6) 0 m2 2 K K 2 where mi is the mass of the i’th nucleus. This in turn leads to the following characteristic equa- tion for the remaining frequencies: m 0 K −K = 1 ω 2 − 0 det − 0 m2 K K = ω 2 − ω 2 − − 2 (m1 K)(m2 K) K and = ω 2(µω 2 − K), ≡ + where µ m1m2/(m1 m2)isthe reduced mass. Thus we pick up another zero frequencyand a nonzero one with ω = √ K/µ,which is the usual result. (Note that we don’tcount ± roots twice--WHY?). 4. Normal Modes in Crystals The main result of the preceding sections is that the characteristic vibrational frequencies are obtained by solving for eigenvalues, cf., Eq. (2.5). For small to mid-size molecules this can be done numerically,onthe other hand, this is not practical for crystalline solids where the matri- ces are huge (3N × 3N with N∼O(N A))! Crystalline materials differ from small molecules in one important aspect; theyare periodic structures made up of identical unit cells; specifically,each unit cell is placed at → = → + → + → R n1a1 n2a2 n3a3,(4.1) → where ni is an integer and ai is a primitive lattice vector.The atoms are positioned within each unit cell at positions labeled by an index, α .Thus, the position of anyatom in entire crystal is determined by specifying R and α .Hence, Eq. (2.5) can be rewritten as λ = ˜ uR,α ΣΣ KR,α ;R′,α ′uR′,α ′.(4.2) R′ α ′ The periodicity of the lattice implies that only the distance between unit cells can matter; i.e., ˜ = ˜ KR,α ;R′,α ′ KR−R′;α ,α ′,and this turns Eq. (4.2) into a discrete convolution that can be simplified by introducing a discrete Fourier transform, i.e., we let ≡ ik⋅R u˜ k,α Σ e uR,α .(4.3) R Winter Term, 2019 Chemistry 365 -5- Normal Mode Analysis When applied to both sides of Eq.
Recommended publications
  • Resonance Beyond Frequency-Matching
    Resonance Beyond Frequency-Matching Zhenyu Wang (王振宇)1, Mingzhe Li (李明哲)1,2, & Ruifang Wang (王瑞方)1,2* 1 Department of Physics, Xiamen University, Xiamen 361005, China. 2 Institute of Theoretical Physics and Astrophysics, Xiamen University, Xiamen 361005, China. *Corresponding author. [email protected] Resonance, defined as the oscillation of a system when the temporal frequency of an external stimulus matches a natural frequency of the system, is important in both fundamental physics and applied disciplines. However, the spatial character of oscillation is not considered in the definition of resonance. In this work, we reveal the creation of spatial resonance when the stimulus matches the space pattern of a normal mode in an oscillating system. The complete resonance, which we call multidimensional resonance, is a combination of both the spatial and the conventionally defined (temporal) resonance and can be several orders of magnitude stronger than the temporal resonance alone. We further elucidate that the spin wave produced by multidimensional resonance drives considerably faster reversal of the vortex core in a magnetic nanodisk. Our findings provide insight into the nature of wave dynamics and open the door to novel applications. I. INTRODUCTION Resonance is a universal property of oscillation in both classical and quantum physics[1,2]. Resonance occurs at a wide range of scales, from subatomic particles[2,3] to astronomical objects[4]. A thorough understanding of resonance is therefore crucial for both fundamental research[4-8] and numerous related applications[9-12]. The simplest resonance system is composed of one oscillating element, for instance, a pendulum. Such a simple system features a single inherent resonance frequency.
    [Show full text]
  • Normal Modes of the Earth
    Proceedings of the Second HELAS International Conference IOP Publishing Journal of Physics: Conference Series 118 (2008) 012004 doi:10.1088/1742-6596/118/1/012004 Normal modes of the Earth Jean-Paul Montagner and Genevi`eve Roult Institut de Physique du Globe, UMR/CNRS 7154, 4 Place Jussieu, 75252 Paris, France E-mail: [email protected] Abstract. The free oscillations of the Earth were observed for the first time in the 1960s. They can be divided into spheroidal modes and toroidal modes, which are characterized by three quantum numbers n, l, and m. In a spherically symmetric Earth, the modes are degenerate in m, but the influence of rotation and lateral heterogeneities within the Earth splits the modes and lifts this degeneracy. The occurrence of the Great Sumatra-Andaman earthquake on 24 December 2004 provided unprecedented high-quality seismic data recorded by the broadband stations of the FDSN (Federation of Digital Seismograph Networks). For the first time, it has been possible to observe a very large collection of split modes, not only spheroidal modes but also toroidal modes. 1. Introduction Seismic waves can be generated by different kinds of sources (tectonic, volcanic, oceanic, atmospheric, cryospheric, or human activity). They are recorded by seismometers in a very broad frequency band. Modern broadband seismometers which equip global seismic networks (such as GEOSCOPE or IRIS/GSN) record seismic waves between 0.1 mHz and 10 Hz. Most seismologists use seismic records at frequencies larger than 10 mHz (e.g. [1]). However, the very low frequency range (below 10 mHz) has also been used extensively over the last 40 years and provides unvaluable information on the whole Earth.
    [Show full text]
  • 22.51 Course Notes, Chapter 9: Harmonic Oscillator
    9. Harmonic Oscillator 9.1 Harmonic Oscillator 9.1.1 Classical harmonic oscillator and h.o. model 9.1.2 Oscillator Hamiltonian: Position and momentum operators 9.1.3 Position representation 9.1.4 Heisenberg picture 9.1.5 Schr¨odinger picture 9.2 Uncertainty relationships 9.3 Coherent States 9.3.1 Expansion in terms of number states 9.3.2 Non-Orthogonality 9.3.3 Uncertainty relationships 9.3.4 X-representation 9.4 Phonons 9.4.1 Harmonic oscillator model for a crystal 9.4.2 Phonons as normal modes of the lattice vibration 9.4.3 Thermal energy density and Specific Heat 9.1 Harmonic Oscillator We have considered up to this moment only systems with a finite number of energy levels; we are now going to consider a system with an infinite number of energy levels: the quantum harmonic oscillator (h.o.). The quantum h.o. is a model that describes systems with a characteristic energy spectrum, given by a ladder of evenly spaced energy levels. The energy difference between two consecutive levels is ∆E. The number of levels is infinite, but there must exist a minimum energy, since the energy must always be positive. Given this spectrum, we expect the Hamiltonian will have the form 1 n = n + ~ω n , H | i 2 | i where each level in the ladder is identified by a number n. The name of the model is due to the analogy with characteristics of classical h.o., which we will review first. 9.1.1 Classical harmonic oscillator and h.o.
    [Show full text]
  • Normal Modes (Free Oscillations)
    Introduction to Seismology: Lecture Notes 22 April 2005 SEISMOLOGY: NORMAL MODES (FREE OSCILLATIONS) DECOMPOSING SEISMOLOGY INTO SUBDISCIPLINES Seismology can be decomposed into three representative subdisciplines: body waves, surface waves, and normal modes of free oscillation. Technically, these domains form a continuum, each pertaining to particular frequency bands, spatial scales, etc. In all cases, these representations satisfy the wave equation, but each is subject to different boundary conditions and simplifying assumptions. Each is therefore relevant to particular types of subsurface investigation. Below is a table summarizing the salient characteristics of the three. Boundary Seismic Domains Type Application Data Conditions Body Waves P-SV SH High frequency travel times; waveforms unbounded dispersion; group c(w) & Surface Waves Rayleigh Love Lithosphere phase u(w) velocities interfaces spherical Normal Modes Spheroidal Modes Toroidal Modes Global power spectra earth As the table suggests, the normal modes provide a framework for representing global seismic waves. Typically, these modes of free oscillation are of extremely low frequency and are therefore difficult to observe in seismograms. Only the most energetic earthquakes are capable of generating free oscillations that are readily apparent on most seismograms, and then only if the seismograms extend over several days. NORMAL MODES To understand normal modes, which describe the modes of free oscillation of a sphere, it’s instructive to consider the 1D analog of a vibrating string fixed at both ends as shown in panel Figure 1b. This is useful because the 3D case (Figure 1c), similar to the 1D case, requires that Figure 1 Figure by MIT OCW. 1 Introduction to Seismology: Lecture Notes 22 April 2005 standing waves ‘wrap around’ and meet at a null point.
    [Show full text]
  • Waves & Normal Modes
    Waves & Normal Modes Matt Jarvis February 2, 2016 Contents 1 Oscillations 2 1.0.1 Simple Harmonic Motion - revision . 2 2 Normal Modes 5 2.1 Thecoupledpendulum.............................. 6 2.1.1 TheDecouplingMethod......................... 7 2.1.2 The Matrix Method . 10 2.1.3 Initial conditions and examples . 13 2.1.4 Energy of a coupled pendulum . 15 2.2 Unequal Coupled Pendula . 18 2.3 The Horizontal Spring-Mass system . 22 2.3.1 Decouplingmethod............................ 22 2.3.2 The Matrix Method . 23 2.3.3 Energy of the horizontal spring-mass system . 25 2.3.4 Initial Condition . 26 2.4 Vertical spring-mass system . 26 2.4.1 The matrix method . 27 2.5 Interlude: Solving inhomogeneous 2nd order di↵erential equations . 28 2.6 Horizontal spring-mass system with a driving term . 31 2.7 The Forced Coupled Pendulum with a Damping Factor . 33 3 Normal modes II - towards the continuous limit 39 3.1 N-coupled oscillators . 39 3.1.1 Special cases . 40 3.1.2 General case . 42 3.1.3 N verylarge ............................... 44 3.1.4 Longitudinal Oscillations . 47 4WavesI 48 4.1 Thewaveequation ................................ 48 4.1.1 TheStretchedString........................... 48 4.2 d’Alambert’s solution to the wave equation . 50 4.2.1 Interpretation of d’Alambert’s solution . 51 4.2.2 d’Alambert’s solution with boundary conditions . 52 4.3 Solving the wave equation by separation of variables . 54 4.3.1 Negative C . 55 i 1 4.3.2 Positive C . 56 4.3.3 C=0.................................... 56 4.4 Sinusoidalwaves ................................
    [Show full text]
  • The 1D Schrödinger Equation for a Free Particle
    Lecture 3 Phys 3750 Two Coupled Oscillators / Normal Modes Overview and Motivation: Today we take a small, but significant, step towards wave motion. We will not yet observe waves, but this step is important in its own right. The step is the coupling together of two oscillators via a spring that is attached to both oscillating objects. Key Mathematics: We gain some experience with coupled, linear ordinary differential equations. In particular we find special solutions to these equations, known as normal modes, by solving an eigenvalue problem. I. Two Coupled Oscillators Let's consider the diagram shown below, which is nothing more than 2 copies of an harmonic oscillator, the system that we discussed last time. We assume that both oscillators have the same mass m and spring constant k s . Notice, however, that because there are two oscillators each has it own displacement, either q1 or q2 . ks m m ks q1 q2 q1 = 0 q2 = 0 Based on the discussion last time you should be able to immediately write down the equations of motion (one for each oscillating object) as ~2 q&&1 +ω q1 = 0, and (1a) ~ 2 q&&2 + ω q2 = 0 , (1b) ~ 2 where ω = k s m . As we saw last time, the solution to each of theses equations is harmonic motion at the (angular) frequency ω~ . As should be obvious from the D M Riffe -1- 1/4/2013 Lecture 3 Phys 3750 picture, the motion of each oscillator is independent of the other oscillator. This is also reflected in the equation of motion for each oscillator, which has nothing to do with the other oscillator.
    [Show full text]
  • What Is a Photon? Foundations of Quantum Field Theory
    What is a Photon? Foundations of Quantum Field Theory C. G. Torre June 16, 2018 2 What is a Photon? Foundations of Quantum Field Theory Version 1.0 Copyright c 2018. Charles Torre, Utah State University. PDF created June 16, 2018 Contents 1 Introduction 5 1.1 Why do we need this course? . 5 1.2 Why do we need quantum fields? . 5 1.3 Problems . 6 2 The Harmonic Oscillator 7 2.1 Classical mechanics: Lagrangian, Hamiltonian, and equations of motion . 7 2.2 Classical mechanics: coupled oscillations . 8 2.3 The postulates of quantum mechanics . 10 2.4 The quantum oscillator . 11 2.5 Energy spectrum . 12 2.6 Position, momentum, and their continuous spectra . 15 2.6.1 Position . 15 2.6.2 Momentum . 18 2.6.3 General formalism . 19 2.7 Time evolution . 20 2.8 Coherent States . 23 2.9 Problems . 24 3 Tensor Products and Identical Particles 27 3.1 Definition of the tensor product . 27 3.2 Observables and the tensor product . 31 3.3 Symmetric and antisymmetric tensors. Identical particles. 32 3.4 Symmetrization and anti-symmetrization for any number of particles . 35 3.5 Problems . 36 4 Fock Space 38 4.1 Definitions . 38 4.2 Occupation numbers. Creation and annihilation operators. 40 4.3 Observables. Field operators. 43 4.3.1 1-particle observables . 43 4.3.2 2-particle observables . 46 4.3.3 Field operators and wave functions . 47 4.4 Time evolution of the field operators . 49 4.5 General formalism . 51 3 4 CONTENTS 4.6 Relation to the Hilbert space of quantum normal modes .
    [Show full text]
  • Normal Modes and Quality Factors of Spherical Dielectric Resonators: I – Shielded Dielectric Sphere
    PRAMANA °c Indian Academy of Sciences Vol. 62, No. 6 | journal of June 2004 physics pp. 1255{1271 Normal modes and quality factors of spherical dielectric resonators: I { Shielded dielectric sphere R A YADAV and I D SINGH Spectroscopy Laboratory, Department of Physics, Banaras Hindu University, Varanasi 221 005, India E-mail: [email protected] MS received 16 July 2002; revised 27 October 2003; accepted 1 April 2004 Abstract. Electromagnetic theoretic analysis of shielded homogeneous and isotropic di- electric spheres has been made. Characteristic equations for the TE and TM modes have been derived. Dielectric spheres of radii of the order of ¹m size are found suitable for the optical frequency region whereas for the microwave region radii of the order of mm size are found suitable. Parameters suitable for their application in the optical and microwave frequency ranges have been used to compute the frequencies corresponding to the normal modes for the TE and TM modes. Expressions for the quality factors for realistic res- onators, i.e., for a dielectric sphere with a non-zero conductivity and a metal shield with a ¯nite conductivity have also been derived for the TE and TM modes. Computations of the quality factors have been made for resonators with parameters suitable for the optical and the microwave regions. Keywords. Eigenmodes; spherical resonators; spherical dielectric resonators; quality factors. PACS No. 42.50.Dv 1. Introduction The earliest reported study on spherical resonators seems to be that of Debye [1] wherein he has studied normal modes of a conducting sphere embedded in a perfect dielectric medium.
    [Show full text]
  • Calculation of Molecular Vibrational Normal Modes
    Calculation of Molecular Vibrational Normal Modes Benjamin Rosman 0407237H September 4, 2008 Supervisor Dr Alex Welte Abstract Normal mode analysis provides a vital key to understanding the dynamics of a complicated system. In this case, this is the motion and vibrations of the atoms in a molecule. It is shown in several test cases that the algorithm successfully detects every normal mode of the molecule, as well as all rigid body rotations and translations. The frequencies of these vibrations are also computed to within 6% of the values reported in the literature. 1 Introduction The energy and wave function describing the states of a molecule are computed by solving the time- independent Schrodinger¨ Equation. However, for large molecules this can involve finding the solution to partial differential equations of several hundred variables, which can become completely infeasible. This calculation can be simplified by use of the Born-Oppenheimer approximation, which allows for the electronic and nuclear problems to be solved independently [Griffiths 2005]. Research is currently being conducted into different ways of solving the electronic problem. The work presented in this document serves to complement that, by discussing the development of tools which would allow for solving the nuclear component of the molecular energy problem, assuming that the electronic problem had already been resolved. The output from the electronic problem is a grid of potentials, which is not something that is directly measurable in a molecule. Instead observables such as frequencies are calculated by this software. Section 2 provides a discussion of normal modes and their particular relevance to quantum chemistry, including the method used in their computation.
    [Show full text]
  • 8 Surface Waves and Normal Modes
    8 Surface waves and normal modes Our treatment to this point has been limited to body waves, solutions to the seismic wave equation that exist in whole spaces. However, when free surfaces exist in a medium, other solutions are possible and are given the name surface waves. There are two types of surface waves that propagate along Earth’s surface: Rayleigh waves and Love waves. For laterally homogeneous models, Rayleigh waves are radially polarized (P/SV) and exist at any free surface, whereas Love waves are transversely polarized and require some velocity increase with depth (or a spherical geometry). Surface waves are generally the strongest arrivals recorded at teleseismic distances and they provide some of the best constraints on Earth’s shallow structure and low-frequency source properties. They differ from body waves in many respects – they travel more slowly, their amplitude decay with range is generally much less, and their velocities are strongly frequency dependent. Surface waves from large earthquakes are observable for many hours, during which time they circle the Earth multiple times. Constructive interference among these orbiting surface waves, to- gether with analogous reverberations of body waves, form the normal modes,or free oscillations of the Earth. Surface waves and normal modes are generally ob- served at periods longer than about 10 s, in contrast to the much shorter periods seen in many body wave observations. 8.1 Love waves Love waves are formed through the constructive interference of high-order SH surface multiples (i.e., SSS, SSSS, SSSSS, etc.). Thus, it is possible to model Love waves as a sum of body waves.
    [Show full text]
  • Normal Mode and Surface Wave Observations
    1.04 Theory and Observations: Normal Mode and Surface Wave Observations G Laske, Scripps Institution of Oceanography, University of California San Diego, La Jolla, CA, USA R Widmer-Schnidrig, Black Forest Observatory, Wolfach, Germany ã 2015 Elsevier B.V. All rights reserved. 1.04.1 Introduction 117 1.04.1.1 Free Oscillations 117 1.04.1.2 Surface Waves 118 1.04.1.3 What We See in Seismograms: The Basics 119 1.04.2 Free Oscillations 120 1.04.2.1 Modes of a Spherically Symmetric Earth 120 1.04.2.2 Modes of a Heterogeneous Earth 122 1.04.2.2.1 Mode splitting 122 1.04.2.2.2 Mode coupling 126 1.04.2.3 Measuring Mode Observables 128 1.04.2.3.1 Multiplet stripping and degenerate mode frequencies 130 1.04.2.3.2 Singlet stripping and receiver stripping 132 1.04.2.3.3 Retrieving the splitting matrix with the matrix autoregressive technique 133 1.04.2.3.4 Retrieving the splitting matrix with iterative spectral fitting 133 1.04.2.3.5 Observed mode coupling 135 1.04.2.4 Example of a Mode Application: Inner-Core Rotation 136 1.04.2.5 Example of a Mode Application: Earth’s Hum 139 1.04.3 Surface Waves 141 1.04.3.1 Standing Waves and Traveling Waves 141 1.04.3.2 The Measurement of Fundamental Mode Dispersion 144 1.04.3.2.1 Group velocity 148 1.04.3.2.2 The study of ambient noise 149 1.04.3.2.3 Phase velocity 149 1.04.3.2.4 Time-variable filtering 151 1.04.3.3 Other Surface Wave Observables 152 1.04.3.4 Higher Modes 154 1.04.3.4.1 Higher-mode dispersion and waveform modeling 154 1.04.3.4.2 Love waves and overtones 157 1.04.3.5 Surface Waves and Structure at Depth 157 1.04.3.6 The Validity of the Great-Circle Approximation and Data Accuracy 160 1.04.4 Concluding Remarks 160 Acknowledgments 161 References 161 1.04.1 Introduction Valdivia, Chile, earthquake.
    [Show full text]
  • Demonstration of Ultra Low Dissipation Optomechanical Resonators on a Chip
    Demonstration of Ultra Low Dissipation Optomechanical Resonators on a Chip G. Anetsberger, R. Rivi`ere,A. Schliesser, O. Arcizet & T. J. Kippenberg∗ Max-Planck-Institut f¨urQuantenoptik, Hans-Kopfermann-Str.1, 85748 Garching Cavity-enhanced radiation-pressure coupling of optical and mechanical degrees of freedom gives rise to a range of optomechanical phenomena, in particular providing a route to the quantum regime of mesoscopic mechanical oscillators. A prime challenge in cavity optomechanics has however been to realize systems which simultaneously maximize optical finesse and mechanical quality. Here we demonstrate for the first time independent control over both mechanical and optical degree of freedom within one and the same on-chip resonator. The first direct observation of mechanical normal mode coupling in a micromechanical system allows for a quantitative understanding of mechanical dissipation. Subsequent optimization of the resonator geometry enables intrinsic material loss limited mechanical Q-factors, rivalling the best values reported in the high MHz frequency range, while simultaneously preserving the resonators' ultra-high optical finesse. Besides manifesting a complete understanding of mechanical dissipation in microresonator based optomechanical systems, our results provide an ideal setting for cavity optomechanics. Over the past years it has become experimentally pos- of both is important for applications such as low loss, sible to study the coupling of optical and mechanical narrow-band \photonic clocks"[2, 27] and indispensable modes via cavity enhanced radiation pressure, which for fundamental studies aiming at approaching and de- gives rise to a diverse set of long anticipated optome- tecting quantized motion in mesoscopic optomechanical chanical phenomena[1] such as radiation pressure driven systems[11, 12, 28].
    [Show full text]