The linear approximation

In the simplest case, one can assume that there is neither global nor local anisotropy, which means that 

The non-diagonal terms of the matrixes Hai and Gai are connected with mutual influence of the particles of the chain. One can admit that, in accordance with the works by Edwards and Freed (1974) and Freed and Edwards (1974, 1975), the hydrodynamic interaction in the system between the particles of the chain becomes negligible, and one can introduce a diagonal matrix of external resistance, but, in virtue of relation (3.5), one cannot introduce non-zero diagonal matrix of internal resistance, so that the simplest forms of the matrixes are 

The symmetrical numerical matrix Ga r represents the influence of motion of the particle 7 on the motion of the particle a. The only general requirement one ought to put on matrix G 7 in the last relation is the following in normal co-ordinates it has a zero eigenvalue. The simplest forms, satisfying the requirement (3.5), can be written as 

equation (3.4) in the simplest case can be specified in the form 

The equation (3.11) is the equation for the dynamics of a single macromolecule in the case of linear dependence on the co-ordinates and velocities.4 Let us note that, if memory functions (3(s) and p(s) turn into -functions, equation (3.11) becomes identical to the equation of motion of the macromolecule in a viscous liquid, which was used in Chapter 2 to describe the dynamics of a macromolecule in this case. 

In the first approximation, neglecting the term proportional to the square of the flattening, we obtain from Equations (2.184 and 2.185) 

Formulas (2.187 and 2.188) constitute the Clairaut s theorem. The last one defines the flattening of the spheroid in terms of the parameter m, known with a sufficient accuracy, and the coefficient/ . The first one. Equation (2.187), gives the law of a distribution of the normal field on the surface of the spheroid. 

Now we evaluate a contribution of terms in Equation (2.188). The gravitational field varies from the equator to pole by approximately 5 Gal. Then, 

It is instructive to derive this equation in a different way. By definition, the equation of the earth s surface is 

Neglecting the square of the parameter /, this equation becomes 

---

A linear stability analysis of (A3.3.57) can provide some insight into the structure of solutions to model B. The linear approximation to (A3.3.57) can be easily solved by taking a spatial Fourier transfomi. The result for the Ml Fourier mode is... 

The limitations and range of validity of the linear theory have been discussed in [17, 23, 24]- The linear approximation to equation (A3.3.54) and equation (A3.3.57) assumes that the nonlinear temis are small compared to the linear temis. As t[increases with time, at some crossover time i the linear... 

In the linear approximation there is a direct Fourier relationship between the FID and the spectrum and, in the great majority of experunents, the spectrum is produced by Fourier transfonnation of the FID. It is a tacit assumption that everything behaves in a linear fashion with, for example, imifonn excitation (or effective RF field) across the spectrum. For many cases this situation is closely approximated but distortions may occur for some of the broad lines that may be encountered in solids. The power spectrum P(v) of a pulse applied at Vq is given by a smc fiinction 18]... 

This treatment illustrates several important aspects of relaxation kinetics. One of these is that the method is applicable to equilibrium systems. Another is that we can always generate a first-order relaxation process by adopting the linearization approximation. This condition usually requires that the perturbation be small (in the sense that higher-order terms be negligible relative to the first-order term). The relaxation time is a function of rate constants and, often, concentrations. 

The quasi-one-dimensional model allows analyzing the behavior of the vapor-liquid system, which undergoes small perturbations. In the frame of the linear approximation the effect of physical properties of both phases, the wall heat flux and the capillary sizes, on the flow instability is studied, and a scenario of the development of a possible processes at small and moderate Peclet number is considered. 

A first approach to the definition of the confidence regions in parameter space follows the linear approximation to the parameter joint distribution that we have already used If the estimates are approximately normally distributed around 9 with dispersion [U. U.] then an approximate 100(1 - a)%... 

Without loss of generality y = y can be assumed. If the dipole moment can be assumed to be a linear function of coordinate within the spread of the frozen Gaussian wave packet, the matrix element (gy,q,p, Pjt(r) Y,q, p ) can be evaluated analytically. Since the integrand in Eq. (201) has distinct maxima usually, we can introduce the linearization approximation around these maxima. Namely, the Taylor expansion with respect to bqp = Qq — Qo and 8po = Po — Po is made, where qj, and pj, represent the maximum positions. The classical action >5qj, p , ( is expanded up to the second order, the final phase-space point (q, p,) to the first order, and the Herman-Kluk preexponential factor Cy pj to the zeroth order. This approximation is the same as the ceUularization procedure used in Ref. [18]. Under the above assumptions, various integrations in U/i(y, q, p ) can be carried out analytically and we have... 

Here I — ma is moment of inertia, a angular acceleration, and z the resultant moment. Note that we have neglected attenuation but in reality, of course, it is always present. This equation characterizes a motion for any angle a, but we consider only the vicinity of points of equilibrium. For this reason, the resultant moment in the linear approximation can be represented as... 

Next let us turn our attention to models described by a set of ordinary differential equations. We are interested in establishing confidence intervals for each of the response variables y, j=l,...,/w at any time t=to. The linear approximation of the output vector at time to,... 

Because the degrees of freedom decouple in the linear approximation, it is easy to describe the dynamics in detail. There is the motion across a harmonic barrier in one degree of freedom and N — 1 harmonic oscillators. Phase-space plots of the dynamics are shown in Fig. 1. The transition from the reactant region at q <0 to the product region at q >0 is determined solely by the dynamics in (pi,qi), which in the traditional language of reaction dynamics is called the reactive mode. 

So far, the discussion of the dynamics and the associated phase-space geometry has been restricted to the linearized Hamiltonian in eq. (5). However, in practice the linearization will rarely be sufficiently accurate to describe the reaction dynamics. We must then generalize the discussion to arbitrary nonlinear Hamiltonians in the vicinity of the saddle point. Fortunately, general theorems of invariant manifold theory [88] ensure that the qualitative features of the dynamics are the same as in the linear approximation for every energy not too high above the energy of the saddle point, there will be a NHIM with its associated stable and unstable manifolds that act as separatrices between reactive and nonreactive trajectories in precisely the manner that was described for the harmonic approximation. 

With the identification of the TS trajectory, we have taken the crucial step that enables us to carry over the constructions of the geometric TST into time-dependent settings. We now have at our disposal an invariant object that is analogous to the fixed point in an autonomous system in that it never leaves the barrier region. However, although this dynamical boundedness is characteristic of the saddle point and the NHIMs, what makes them important for TST are the invariant manifolds that are attached to them. It remains to be shown that the TS trajectory can take over their role in this respect. In doing so, we follow the two main steps of time-independent TST first describe the dynamics in the linear approximation, then verify that important features remain qualitatively intact in the full nonlinear system. 

Because in an autonomous system many of the invariant manifolds that are found in the linear approximation do not remain intact in the presence of nonlinearities, one should expect the same in the time-dependent case. In particular, the separation of the bath modes will not persist but will give way to irregular dynamics within the center manifold. At the same time, one can hope to separate the reactive mode from the bath modes and in this way to find the recrossing-free dividing surfaces and the separatrices that are of importance to TST. As was shown in Ref. 40, this separation can indeed be achieved through a generalization of the normal form procedure that was used earlier to treat autonomous systems [34]. 

The scaling (58) assigns order s° to the leading term of the autonomous Hamiltonian Hsys. The leading term of the time-dependent part is assigned order s-1. This term is treated exactly in the linear approximation (33). It should therefore have order s°. To achieve this, we also have to scale... 

For the second point of the linear approximation the maximum A T has to be found. Assuming that the highest temperature Tout max can be reached after a complete desorption and that all of the water vapor within the air stream will be adsorbed ( Ax = x n), Tmax can be written as... 

This point can be easily calculated from the adsorption equilibrium of each adsorbent. As shown in Figures 2 and 3 the linear approximation gives a sufficient accurate estimates ( 5% of the experimental values) within a range of realistic conditions to predict Tout and xout under given de- and adsorption conditions for each adsorbent [4],... 

During an early time period in which approximately less than 5% of the mass is transported, the kinetics are effectively linear. Subsequently, the linear approximations to Eqs. (19) and (22) are... 

I presented a group of subroutines—CORE, CHECKSTEP, STEPPER, SLOPER, GAUSS, and SWAPPER—that can be used to solve diverse theoretical problems in Earth system science. Together these subroutines can solve systems of coupled ordinary differential equations, systems that arise in the mathematical description of the history of environmental properties. The systems to be solved are described by subroutines EQUATIONS and SPECS. The systems need not be linear, as linearization is handled automatically by subroutine SLOPER. Subroutine CHECKSTEP ensures that the time steps are small enough to permit the linear approximation. Subroutine PRINTER simply preserves during the calculation whatever values will be needed for subsequent study. 

For small overpotentials the linear approximations of Eqs. (5.4) and (5.5) should be sufficient, but at high overpotentials higher-order terms are expected to contribute. 

The error in this linear approximation approaches zero proportionally to (Ax)2 as Ax approaches zero. Given initial values for the variables, all nonlinear functions in the problem are linearized and replaced by their linear Taylor series approximations at this initial point. The variables in the resulting LP are the Ax/s, representing changes from the base vaiues. In addition, upper and lower bounds (called step bounds) are imposed on these change variables because the linear approximation is reasonably accurate only in some neighborhood of the initial point. 

Note that as the line search process continues and the total step from the initial point gets larger, the number of Newton iterations generally increases. This increase occurs because the linear approximation to the active constraints, at the initial point (0.697,1.517), becomes less and less accurate as we move further from that point. 

Combining equations (5.4.86) to (5.4.88) gives the linear approximation of the matrix ffSx pftT as... 

---