Third-order integrators

The choices p = q = 1 and p = q = 2 are especially useful. The former yields a third-order integration formula, the latter a fifth-order formula. For example, for p = q = 1, we obtain from Equation 25 ... 

Find the integrated rate equation for a third-order reaction having the rate equation —dc/ ldt = kCf,. ... 

If the hydroxyl and carboxyl group concentrations are equal, both being given by c, Eq. (6) may be replaced by the standard integrated expression for a third-order reaction ... 

If the process is uniformly of the third order without change in velocity constant throughout, the integration constant will be 1/co where co is the initial concentration of the functional groups. It is convenient here, and for many other purposes as well, to introduce a parameter called the extent of reaction and designated by p, which represents the fraction of the functional groups initially present that have undergone reaction at time t. Then c= (1 —p)co, and Eq. (7) may be replaced by... 

The function 7(f) can be chosen for the whole reaction time interval, or two or three subsequent temperature-time data points 7(fi-i), 7(fi), and 7(fi+i) can be approximated by polynomials of second or third order 7,(f), respectively. These polynomials will then be used in a procedure for numerical integration in each integration step i. This method has been successfully applied in a kinetic study of the partial oxidation of hydrocarbons (Skrzypek et al., 1975, Krajewski etai, 1975, 1976, 1977). 

Figure 13 Timing diagram for the clean HMBC experiment with an initial second-order and terminal adiabatic low-pass 7-filter.42,43 The recommended delays for the filters are the same than for a third-order low-pass J filter. <5 and 8 are gradient delays, where 8 — <5 + accounts for the delay of the first point in the 13C dimension. The integral over each gradient pulse G, is H/2yc times the integral over gradient G2 in order to achieve coherence selection. The recommended phase cycle is c/)n = x, x, x, x 3 — 4(x), 4(y), 4( x), 4(—y) with the receiver phase c/)REC = x, x.

The integrals are found numerically. They are tabulated and support the third order mechanism. 

A third order reaction can be the result of the reaction of a single reactant, two reactants or three reactants. If the two or three reactants are involved in the reaction they may have same or different initial concentrations. Depending upon the conditions the differential rate equation may be formulated and integrated to give the rate equation. In some cases, the rate expressions have been given as follows. 

The reaction is first-order with respect to X, second-order with respect to Y and third-order overall k is the rate constant. Note that the order of a reaction does not necessarily have to have integral values. 

Third-order results for closed-shell molecules have average absolute errors of 0.6 - 0.7 eV [31]. Transformed integrals with four virtual indices and OV4 contractions for each value of E are required for the U intermediate, which is needed for ionization energy as well as electron affinity calculations. 

We need to study the numerical integration of only first-order ODEs. Any higher-order equations, say with Mth-order derivatives, can be reduced to N first-order ODEs. For example, suppose we have a third-order ODE ... 

Integration Method or Hit and Trial Method Here known quantities of standard solutions of reactants are mixed in a reaction vessel and the progress of the reaction is determined by determining the amount of reactant consumed after different intervals of time. These values are then substituted in the equations of first, second, third order and so on. The order of the reaction is the order corresponding to that equation which gives a constant value of K. 

Figure 2 Kinetics of gas-phase propylene homometathesis at 0°C, catalyzed by (a) perrhenate/silica-alumina activated by SnMe4 (10 mg, 0.83 wt % Re) and (b) MeReOs on HMDS-capped silica-alumina (10 mg, 1.4 wt % Re). Solid lines are curve-fits to the first-order integrated rate equation. Solid squares first addition solid circles second addition open circles third addition of propylene (30 Torr) to the catalyst.

Note that in die leapfrog method, position depends on the velocities as computed one-half time step out of phase, dins, scaling of the velocities can be accomplished to control temperature. Note also that no force-deld calculations actually take place for the fractional time steps. Eorces (and thus accelerations) in Eq. (3.24) are computed at integral time steps, halftime-step-forward velocities are computed therefrom, and these are then used in Eq. (3.23) to update the particle positions. The drawbacks of the leapfrog algorithm include ignoring third-order terms in the Taylor expansions and the half-time-step displacements of the position and velocity vectors - both of these features can contribute to decreased stability in numerical integration of the trajectoiy. 

Let the value of 8 at node i be represented by bi. Then, we can integrate equation (1) analytically if we assume an appropriate interpolation formula within each element. Linear interpolation was adopted in Nariai and Shigeyama(1984) while third order polynomials were used in Nariai and Ito(1985). Results can be expressed in a form... 

This expression is actually sometimes used to define the first and second order density matrix. It is anyway useful to know that the first order density matrix elements are equal to the coefficients for the corresponding one-electron integrals in the energy expression and similarly for the second order density matrix elements. The definition of the third order density matrix is,... 

We have presented a review of the salient features of nonlinear integrated optics. It appears that nonlinear organic materials can play an important role in second- and third-order guided-wave devices. This field requires a great deal of material characterization and processing, however, before significant advances are realized. 

---