Fourth-order, generally

Here c[-], which will be called the elastic modulus tensor, is a fourth-order linear mapping of its second-order tensor argument, while b[-], which will be called the inelastic modulus tensor, is a linear mapping of k whose form will depend on the specific properties assigned to k. They depend, in general, on and k. For example, if k consists of a single second-order tensor, then in component form... 

The work of the present section shows that shock-compression experiments provide an effective method for determination of higher-order elastic properties and that, by the same token, the effects of nonlinear elastic response should generally be taken into account in investigations of shock compression (see, e.g., Asay et al. [72A02]). Fourth-order contributions are readily apparent, but few coefficients have been accurately measured. 

This is the general linear equation of motion for an almost planar and rough one-dimensional phase boundary. The fourth-order term in the spatial derivative acts as a stabilizer just like the second-order term, and is not really crucial here. 

As was proven later by Bishop [19], the coefficient A in the expansion (73) is the same for all optical processes. If the expansion (73) is extended to fourth-order [4,19] by adding the term the coefficient B is the same for the dc-Kerr effect and for electric field induced second-harmonic generation, but other fourth powers of the frequencies than are in general needed to represent the frequency-dependence of 7 with process-independent dispersion coefficients [19]. Bishop and De Kee [20] proposed recently for the all-diagonal components yaaaa the expansion... 

In this chapter we described Euler s method for solving sets of ordinary differential equations. The method is extremely simple from a conceptual and programming viewpoint. It is computationally inefficient in the sense that a great many arithmetic operations are necessary to produce accurate solutions. More efficient techniques should be used when the same set of equations is to be solved many times, as in optimization studies. One such technique, fourth-order Runge-Kutta, has proved very popular and can be generally recommended for all but very stiff sets of first-order ordinary differential equations. The set of equations to be solved is... 

If the characteristic polynomial passes the coefficient test, we then construct the Routh array to find the necessary and sufficient conditions for stability. This is one of the few classical techniques that we do not emphasize and the general formula is omitted. The array construction up to a fourth order polynomial is used to illustrate the concept. 

Ctjki is a fourth order tensor that linearly relates a and e. It is sometimes called the elastic rigidity tensor and contains 81 elements that completely describe the elastic characteristics of the medium. Because of the symmetry of a and e, only 36 elements of Cyu are independent in general cases. Moreover only 2 independent rigidity constants are present in Cyti for linear homogeneous isotropic purely elastic medium Lame coefficient A and /r have a stress dimension, A is related to longitudinal strain and n to shear strain. For the purpose of clarity, a condensed notation is often used... 

Here, Flffl are matrix elements of a zeroth-order Hamiltonian, which is chosen as a one-electron operator in the spirit of MP2. is an overlap matrix The excited CFs are not in general orthogonal to each other. Finally, Vf)(i represents the interaction between the excited function and the CAS reference function. The difference between Eq. [2] and ordinary MP2 is the more complicated structure of the matrix elements of the zeroth-order Hamiltonian in MP2 it is a simple sum of orbital energies. Here H is a complex expression involving matrix elements of a generalized Fock operator F combined with up to fourth-order density matrices of the CAS wave function. Additional details are given in the original papers by Andersson and coworkers.17 18 We here mention only the basic principles. The zeroth-order Hamiltonian is written as a sum of projections of F onto the reference function 0)... 

In general, the system (26)-(30) is nonintegrable and its dynamics has to be studied numerically. We examined it with the help of a fourth-order Runge-Kutta... 

Computationally the super-CI method is more complicated to work with than the Newton-Raphson approach. The major reason is that the matrix d is more complicated than the Hessian matrix c. Some of the matrix elements of d will contain up to fourth order density matrix elements for a general MCSCF wave function. In the CASSCF case only third order term remain, since rotations between the active orbitals can be excluded. Besides, if an unfolded procedure is used, where the Cl problem is solved to convergence in each iteration, the highest order terms cancel out. In this case up to third order density matrix elements will be present in the matrix elements of d in the general case. Thus super-CI does not represent any simplification compared to the Newton-Raphson method. 

This is a tensor of fourth order, and in the general case it should be described by a matrix of 81 members (9x9). Since the stress and strain tensors are symmetrical and each has six independent components, the tensor of fourth order derived from them has 6x6 components. 

The equation f(x) = g(x) is a fourth-order algebraic equation, hence to write down the conditions (23) in the explicit form for the general case is difficult. An explicit form of the multiplicity criterion for eqns. (23) solutions can be obtained, e.g. from the simple demand for eqns. (23) to account for the inflexion point x for the f(x) function. Then from f (x ) = 0 we obtain... 

A method similar to the iterative, is the partial closure method [37], It was formulated originally as an approximated extrapolation of the iterative method at infinite number of iterations. A subsequent more general formulation has shown that it is equivalent to use a truncated Taylor expansion with respect to the nondiagonal part of T instead of T-1 in the inversion method. An interpolation of two sets of charges obtained at two consecutive levels of truncations (e.g. to the third and fourth order) accelerates the convergence rate of the power series [38], This method is no longer in use, because it has shown serious numerical problems with CPCM and IEFPCM. 

---