Matrix differential propagation

Once obtained, the differential propagation Jones matrix, N, can be used to con-... 

Equation (2.52) satisfies the requirement that the ordering of the lamellae will not affect the final result. Furthermore, the differential propagation Jones matrix for a composite material is simply the sum of the matrices, N., for each of the separate optical effects. The N matrices for different optical properties can be derived using equation (2.44). As an example, consider a birefringent material with retardation S = (InAn z) /X. One form of the Jones matrix for this material is... 

Formally, if one has the experimental values of the dielectric tensor e, the magnetic permeability tensor /jl, and the optical rotation tensors p and p for the substrate, one can construct first the optical matrix M, then the differential propagation matrix A, and C, which, to repeat, is the x component of the wavevector of the incident wave. Once A is known, the law of propagation (wave equation) for the generalized field vector ift (the components of E and H parallel to the x and y axes) is specified by Eq. (2.15.18). Experimentally, one travels this path backwards. 

Consider the relationship between A and the dielectric tensor e. In ellipsometry, there is reflection and transmission by the surface (z = 0) of a semi-infinite anisotropic substrate (biaxial crystal) into an isotropic ambient (air, for z<0). Suppose that this semi-infinite anisotropic medium (the crystal) is homogeneous and that its optical matrix M is independent of z (if A does depend on z—that is, on how far into the crystal one goes—then the problem becomes much more difficult). If the optical matrix M of the substrate is independent of z, then so is the differential propagation matrix A if A is independent of z and has a value (, to be found below, the solution of Eq. (2.15.25) is given by... 

So far, however, one still needs an expression for the reflection matrix that shows how to extract from it the tensor elements for the refractive index tensor of the biaxial medium. We seek the reflection matrix R for the semi-infinite anisotropic biaxial medium. Using Eq. (2.15.8) and Eq. (2.15.21), we can relate the 4x4 differential propagation matrix A to the dielectric tensor e from Eqs. (2.15.21) and (2.15.24). Then it can be shown that... 

We shall call Q z) the differential propagation matrix. When (z) does not vary appreciably over an interval h, an integral of this equation is... 

The 2-independent term Po(A) is the propagation matrix corresponding to the invariant part,, of the differential propagation matrix. It may be written in the following exact, closed form. 

The correction to the relaxing density matrix can be obtained without coupling it to the differential equations for the Hamiltonian equations, and therefore does not require solving coupled equations for slow and fast functions. This procedure has been successfully applied to several collisional phenomena involving both one and several active electrons, where a single TDHF state was suitable, and was observed to show excellent numerical behavior. A simple and yet useful procedure employs the first order correction F (f) = A (f) and an adaptive step size for the quadrature and propagation. The density matrix is then approximated in each interval by... 

The Eik/TDDM approximation can be computationally implemented with a procedure based on a local interaction picture for the density matrix, and on its propagation in a relax-and-drive perturbation treatment with a relaxing density matrix as the zeroth-order contribution and a correction due to the driving effect of nuclear motions. This allows for an efficient computational procedure for differential equations coupling functions with short and long time scales, and is of general applicability. 

This is the simplest of the models where violation of the Flory principle is permitted. The assumption behind this model stipulates that the reactivity of a polymer radical is predetermined by the type of bothjts ultimate and penultimate units [23]. Here, the pairs of terminal units MaM act, along with monomers M, as kinetically independent elements, so that there are m3 constants of the rate of elementary reactions of chain propagation ka ]r The stochastic process of conventional movement along macromolecules formed at fixed x will be Markovian, provided that monomeric units are differentiated by the type of preceding unit. In this case the number of transient states Sa of the extended Markov chain is m2 in accordance with the number of pairs of monomeric units. No special problems presents writing down the elements of the matrix of the transitions Q of such a chain [ 1,10,34,39] and deriving by means of the mathematical apparatus of the Markov chains the expressions for the instantaneous statistical characteristics of copolymers. By way of illustration this matrix will be presented for the case of binary copolymerization ... 

It is required to calculate the Mueller elements by first calculating the differential elements for the slab using equation (2.3) and integrating along the path of light propagation as indicated in equation (5.51). The Mueller matrix elements are... 

Baluja, K., Burke, P.G. and Morgan, L.A. (1982). R-matrix propagation program for solving coupled second-order differential equations, Comput. Phys. Commun. 27, 299-307. 

As to future directions, the problem of the canonical density matrix, or equivalently the Feynman propagator, for hydrogen-like atoms in intense external fields remain an unsolved problem of major interest. Not unrelated, differential equations for the diagonal element of the canonical density matrix, the important Slater sum, are going to be worthy of further research, some progress having already been made in (a) intense electric fields and (b) in central field problems. Finally, further analytical work on semiclassical time-dependent theory seems of considerable interest for the future. 

---