System trajectory

Having specified the interactions (i.e., the model of the system), the actual simulation then constructs a sequence of states (or the system trajectory) in some statistical mechanical ensemble. Simulations can be stochastic (Monte Carlo (MC)) or deterministic (MD), or they can combine elements of both, such as force-biased MC, Brownian dynamics, or generalized Lan-gevin dynamics. It is usually assumed that the laws of classical mechanics (i.e., Newton s second law) may adequately describe the atoms and molecules in the physical system. 

Figures 4.44 and 4.45, best viewed in color, show a benign complication of the problem caused by the Lewis numbers. If, however, we reduce the Lewis number LeA further to 0.07, the system trajectories indicate periodic explosions of the underlying system throughout all time, and the trajectories do not converge to the steady state at all, even with what we thought to be proper feedback. The trajectory that these curves settle at is called a periodic attractor of the system in contradistinction to the earlier encountered point attractor of Figures 4.43 or 4.44, for example. A point attractor, or more accurately a fixed-point attractor, is a more commonly encountered steady state in chemical and biological engineering systems. It could be called a stationary nonequilibrium state to distinguish it from the stationary equilibrium states associated with closed or isolated batch processes.

Instead of enhancing the performance of computers, many theoreticians have tried to enhance the efficiency of computation by improving computational codes. One of the approaches is to reduce the dependence of the computation for the present step on the computational results for previous steps along the system trajectory or increase the parallelism of the computation. The other effective approach is the use of concept of hierarchical coupling of paradigms... 

Let us now examine the behaviour of the solutions for the dynamic system (20) in time and analyze the system trajectories in the phase pattern. This analysis permits us to characterize peculiarities of the unsteady-state behaviour (in particular to establish whether the steady state is stable or unstable), to determine its type (focus, node, saddle, etc.) and to find attraction regions for stable steady states, singular lines, etc. 

Fig. 20. Variation of steady states with temperature on the phase pattern Pq2 = 2.3 x 10 7, Fco = 2.2 x 10 7 Torr. Points stable steady states at given temperatures (K). Lines system trajectories at stepwise temperature variations (533 to 534K and 485 to 484K).

As an introduction to relativistic dynamics, it is of interest to treat time as a dynamical variable rather than as a special system parameter distinct from particle coordinates. Introducing a generic global parameter r that increases along any generalized system trajectory, the function t(r) becomes a dynamical variable. In special relativity, this immediately generalizes to A (r) for each independent particle, associated with spatial coordinates x (r). Hamilton s action integral becomes... 

SA = 0 subject to the energy constraint restates the principle of least action. When the external potential function is constant, the definition of ds as a path element implies that the system trajectory is a geodesic in the Riemann space defined by the mass tensor m . This anticipates the profound geometrization of dynamics introduced by Einstein in the general theory of relativity. 

Predictions can be made about the suitability of different system trajectories on the basis of orbital symmetry conservation rules (207). The most suitable trajectory is an approximation to the reaction path of the reaction under study. The rules can also yield information about the possible structure of the activated complex. The correlation diagram technique has been improved in a series of books by Epiotis et al. (214-216). The method is based on self-consistent field-configuration interaction or valence bond (SCF-CI or VB) (including ionic structures) wave functions. Applications on reactions in the ground states as well as in the excited electronic states are impressive however, the price to be paid for the predictions seems to be rather high. 

Under the copolymerization of more than two monomers Eqs. (5.3) cannot be integrated explicitly, and in order to determine the system trajectories one should need the numerical calculations. Examples of such calculations of the conversional change of composition and structure characteristics of the terpolymers have been reported in Refs. [195-200]. One should pay special attention to Ref. [200] where the programs for the computer realization of such calculations are presented. Under the copolymerization of four or more monomers, the composition drift with the conversion was calculated [7,8] only within the framework of the simplified terminal model described above in Sect. 4.6. 

Figure 5. The distributions of the recrossing trajectories over configurational surface S qi = 0) at time t = 0 on the phase-space planes (pf (p,q), (p,q)) at E = 0.5e, where most modes are strongly chaotic—except 4i(p,q). (a) First and (b) second orders The circle and triangle symbols denote the system trajectories having negative and positive incident momenta p (t = 0) on the S(qi = 0), and the open and filled symbols denote those whose final states were predicted correctly and falsely by Eq. (11), respectively [45].

Symbolic dynamics is a powerful, but qualitative, tool. We can easily imagine a system trajectory that is close to the trajectory shown in Fig. 2.11(a), but not quite identical. Both trajectories axe then characterized by the same word, ABC AC. Thus, symbol sequences do not specify system trajectories uniquely. In other words, there is no one-to-one correspondence between the impact parameters b and the words constructed from the alphabet A, B, C. ... 

The quantities pa are called generahzed momenta. They can be used together with the coordinates qa to define a system trajectory. The system trajectory evolves in the 2/-dimensional space spanned by the / coordinates q and the / coordinates p. This space plays a central role in analytical mechanics. It is called the phase space of the system. A point in phase space uniquely defines the mechanical state of a system. In connection with Poincare s method of surfaces of section, the phase space is also an important vehicle for the visuahzation of the quahtative behaviour of a given dynamical system. An example is presented in Section 3.2. 

For a given set of initial conditions Zh t = 0), k = 1,. ..,4, the solution of (3.2.6) defines a system trajectory. Since the energy is conserved, this trajectory winds in a three-dimensional sub-space of the four-dimensional phase space of the pendulum. It is not easy to imagine the motion in such a high-dimensional space, and we need to devise a visualization method that gives some clues as to the qualitative nature of the system trajectories. One particularly useful method, the method of the surface of... 

A model of competition between plasmid-bearing and plasmid-free organisms in a chemostat was proposed by Stephanopoulos and Lapidus [SLa], who give a local analysis of various models. A global analysis of the behavior of system trajectories was presented in [HWW], and a portion of that analysis is sketched here. The major remaining problem is discussed after the description of the known results. There are other models of plasmid loss (and conjugation), for example, Stewart and Levin [SL2] and Macken, Levin, and Waldstatter [MLW]. The survey article of Si-monsen [Si] contains a discussion of the experiments and the theory. 

Figure 5. Schematic variation in the electron potential energy as a function of the nuclear coordinates (qa and qp for the reactant and product systems, respectively) and of the electron coordinate, (a) Forbidden electron transfer (b) allowed electron transfer (c) projection of the system trajectory on the electron coordinate-nuclear coordinates plane, q and qp are the values of the nuclear coordinate at the equilibrium for the reactant or product systems, respectively, and q that at the transition state. Solid lines, variations of the nuclear coordinates ( ->) electron tunneling.

This observation has an important practical consequence In numerical simulation we usually follow a single-system trajectory in time, and the system temperature can be obtained from such an equilibrium trajectory using Eq. (5.14)." Note that... 

It is evident that the system (1.1) for a < 0 (or a > 0) is structurally stable with respect to changes in the parameter a a small change of a (not altering its sign) does not change the form of a trajectory on the phase plane [x(t), y(t)]. On the other hand, the system (1.1) for a = 0 is not structurally stable — an arbitrarily small perturbation of the parameter a change the form of the system trajectory. 

The Lie techniques may provide us with the physical footings or analytical means to elucidate dynamical correlations among successive saddle crossings by enabling us to scrutinize connectivity of manifolds from and to the sequential saddle points and extent of volume of the region of a junction of manifolds in terms of the backward system trajectories initiated from S (gi(p, q) = 0) at one saddle point and the forward from the other S (gi(p,q) = 0) at the previous saddle point, through which the system has passed before reaching the first [80]. 

The mechanical system being controlled is referred to as the plant. The configuration of the system at any instant in time comprises the plant states. The devices that power the system are called actuators. The signals driving the actuators are called controls. The controller encompasses processes by which the controls are generated. The time histories of the plant states in response to the control signals are referred to the system trajectory. 

Integrating the eqiaations of motion(14.II)at the above conditions, the classical trajectories x (t) (or the interatomic distances) as functions of time and the system trajectory in configuration space... 

---