Dynamic equations, eigenvectors

The L are the eigenvectors of the dynamical equation 5 ( or 7) and can be calculated in the usual way once the force field (F matrix, eq.2) and the geometry (O matrix, eq.3) of the systems are known. This has been done for Polyacetylene (paper I) and for the systems discussed in this paper. 

The quantities U Ix)j and UFx)j in (13) are projections of the eigenvector j along lx- From the above equations, we can interpret these as follows. The term UFx)j is the amount that the transition j received from the total X magnetization, created from the equilibrium state, and (U Ix)j is how much that transition contributes to the observed signal. These two terms may not be equal, as we see in exchanging systems. This general approach forms the basis of the description of dynamic NMR lineshapes. 

The fast stage of relaxation of a complex reaction network could be described as mass transfer from nodes to correspondent attractors of auxiliary dynamical system and mass distribution in the attractors. After that, a slower process of mass redistribution between attractors should play a more important role. To study the next stage of relaxation, we should glue cycles of the first auxiliary system (each cycle transforms into a point), define constants of the first derivative network on this new set of nodes, construct for this new network an (first) auxiliary discrete dynamical system, etc. The process terminates when we get a discrete dynamical system with one attractor. Then the inverse process of cycle restoration and cutting starts. As a result, we create an explicit description of the relaxation process in the reaction network, find estimates of eigenvalues and eigenvectors for the kinetic equation, and provide full analysis of steady states for systems with well-separated constants. 

In the simplest case, the auxiliary discrete dynamical system for the reaction network W is acyclic and has only one attractor, a fixed point. Let this point be A (n is the number of vertices). The correspondent eigenvectors for zero eigenvalue are r = S j and Z = 1. For such a system, it is easy to find explicit analytic solution of kinetic equation (32). 

As it is demonstrated, dynamics of this system approximates relaxation of the whole network in subspace = 0. Eigenvalues for Equation (45) are —k, (i < n), the corresponded eigenvectors are represented by Equations (34), (36) and zero-one multiscale asymptotic representation is based on Equations (37) and (35). 

The second case shows very different behavior The relative concentrations of the degenerate master sequences are subject to random drift, and the dominant eigenvector of W represents at best a time average of the mutant distribution. Then the dynamics can be modelled only by a stochastic process requiring careful choice of the appropriate mathematical technique and approximations in a hierarchy of equations (see refs. 48 and 51 and Section V.2). One difficulty here is that even very distant mutants contribute if sufficiently neutral. The results of Section III.2 indicate that there is (almost... 

We now introduce the ideas of Weyl to distinguish between pure states and mixtures. Pure states were mathematically represented by eigenvectors of observables, which described the properties of a particle or a dynamic state. On the other hand, mixtures were composed of pure states of a certain mixing relationship. These aspects are clearly important to chemists and obviously to the electrochemists too. The canonical variables, G and H [19], have to satisfy the canonical or Heisenberg commutation relation, derived from Equations 3.12 and 3.13 ... 

An efficient way to treat such a system is to assemble all coefficients of the different terms of the mass-balance equations in a matrix and to apply methods of matrix algebra to solve the system for steady-state concentrations (level III) or for the concentrations as functions of time (level IV) [19]. We denote the matrix of coefficients (the fate matrix ) by S, the vector of concentrations in all boxes of the model by c, and the vector of all source terms by q. The set of mass-balance equations describing the temporal changes of the concentrations in all boxes then reads c = -S c + q. The steady-state solution is obtained by setting c equal to zero and solving for c. This leads to ss -1. j obtain the steady-state concentrations the emission vector has to be multiplied by the inverse of the matrix S. For the dynamic solutions of the system, the eigenvalues and eigenvectors of S have to be determined. 

Other variants are due to Fano [76], Anderson [77], Lee [78], and Friedrichs [79] and have been successfully applied to study, for example, autoionization, photon emission, or cavities coupled to waveguides. The dynamics can be solved in several ways, using coupled differential equations for the time-dependent amplitudes and Laplace transforms or finding the eigenstates with Feshbach s (P,Q) projector formalism [80], which allows separation of the inner (discrete) and outer (continuum) spaces and provides explicit expressions ready for exact calculation or phenomenological approaches. For modern treatments with emphasis on decay, see Refs. [31, 81]. Writing the eigenvector as [31, 76]... 

This decomposition can be interpreted as a PCA on the bandpass-filtered process. An inverse Fourier transform using the first few eigenvectors provides a reduced representation of the original multivariate time series [25]. When anarrowband frequency-domain structure is present in the dynamics, it can be shown that Equation 5.20 provides stronger optimal decomposition of the multivariate time series than does conventional PCA [25]. 

A well-known technique for dynamic analysis of structures is the Modal Analysis. It consists of obtaining the natural frequencies and their associated modeshapes. Mathematically, modal analysis is performed by obtaining the eigenvalues and eigenvectors of the homogeneous equation ... 

It is widely accepted that vibrational dynamics of atoms and molecules are reasonably well represented by harmonic force fields. The resolution of the secular equation transforms a set of (say N) coupled oscillators into (N) independent oscillators along orthogonal (normal) coordinates. Eigenvalues of the dynamical matrix are normal frequencies and eigenvectors give atomic displacements for each normal mode [3,4,13]. If band intensities cannot be frilly exploited, as it is normally the case for infrared and Raman spectra, these vectors are unknown and force fields refined with respect to observed frequencies only are largely underdetermined. For complex systems, symmetry consideration or/and isotopic substitutions may remove only partially this under determination. 

Excitation energies are thus computed as poles of the dynamic polarizability, that is, as the values of co leading to zero eigenvalues on the left-hand side of the matrix of Eq. 1.12. In the framework of the above equations, an efficient fast iterative solution for the lowest eigenvalue/excitation energies can be attained [62]. Oscillator strengths can also be obtained by the eigenvectors of Eq. 1.12, as explained by Casida [17]. 

The problem of the calculation of eigenvalues and eigenvectors becomes even more complex when the translational symmetry is destroyed by the introduction of defects. The concept of phonon waves is lost and the k dependency in equation 1 and 2 is lost. One has then to calculate eigenvalues and eigenvectors of very large dynamical matrices corresponding to the size of the piece of crystal or segment of polymer chain one has chosen. Even if drastic structural simplifications are introduced, the problem still remains untractable by the standard numerical procedures. 

The solution of the dynamic problem, which is modeled by Eq, (2.18), would require the evaluation of the characteristic values (eigenvalues) X, and characteristic vectors (eigenvectors) jC( of the matrix K. It is shown in Chap. 5 that the solution of a set of linear ordinary differential equations can be obtained by using Eq. (5.53) ... 

This text of the two-volume treatment contains most of the theoretical background necessary to understand experiments in the field of phonons. This background is presented in four basic chapters. Chapter 2 starts with the diatomic linear chain. In the classical theory we discuss the periodic boundary conditions, equation of motion, dynamical matrix, eigenvalues and eigenvectors, acoustic and optic branches and normal coordinates. The transition to quantum mechanics is achieved by introducing the Sohpddingev equation of the vibrating chain. This is followed by the occupation number representation and a detailed discussion of the concept of phonons. The chapter ends with a discussion of the specific heat and the density of states. 

We start the discussion by formulating the Hamiltonian of the system and the equations of motion. The concept of force constants needs further examination before it can be applied in three dimensions. We shall discuss the restrictions on the atomic force constants which follow from infinitesimal translations of the whole crystal as well as from the translational symmetry of the crystal lattice. Next we introduce the dynamical matrix and the eigenvectors this will be a generalization of Sect.2.1.2. In Sect.3.3, we introduce the periodic boundary conditions and give examples of Brillouin zones for some important structures. In strict analogy to Sect.2.1.4, we then introduce normal coordinates which allow the transition to quantum mechanics. All the quantum mechanical results which have been discussed in Sect.2.2 also apply for the three-dimensional case and only a summary of the main results is therefore given. We then discuss the den-... 

---