Numerical Considerations

Our numerical task is to solve the IVP given by the eight DEs (4.103) to (4.106) both for j = 1 and j = 2 with the given 26 system parameters and the two rate equations (4.107) and (4.108) for any physically feasible set of initial values b(0) and feed parameters k fee,i numerically. 

A numerical solution of the problem may consist of a plot of all eight profiles h (t), su(t), S2i(t), S3i(f), /12(f), 312(f), 322(f), and 333(f) f°r a certain time interval 0 t Tend. Or it may involve phase plots, such as that of the acetylcholine concentration in compartment (I) versus that in compartment (II). Our aim in the computations that follow is to show the variations in the quality of the solutions for differing values of hf, the hydrogen ion concentration of the feed to compartment (I). Another dependency of the solutions is explored in the exercises. 

Throughout this section we work with the initial values vector 

Here is the MATLAB program neurocycle. m, which upon the user s specification and using MATLAB s (un)commenting feature, which places the symbol % at the start of comment lines, either plots all eight profiles, or one profile and one phase plot, or only one phase plot for relatively high values of the time parameter t, when the system has reached its periodic limit cycle. 

The commenting or uncommenting of MATLAB code line blocks can best be achieved from the MATLAB text editor window for an m file. Simply highlight a block of code lines via a mouse drag in the MATLAB text editor window, then click on the Text entry of the editor s toolbar and click Comment or Uncomment as appropriate. This action makes % commenting marks appear at or disappear from the front of each code line of the highlighted block. 

As mentioned before, the solution of the time-dependent Kohn-Sham equations is an initial value problem. At t = to the system is in some initial state described by the Kohn-Sham orbitals to). In most cases the initial state will be the ground state of the system (i.e., (fi(r,to) will be the solution of the ground-state Kohn-Sham equations). The main task of the computational physicist is then to propagate this initial state until some final time, tf. 

The time-dependent Kohn-Sham equations can be rewritten in the integral form 

Note that Hks is explicitly time-dependent due to the Hartree and xc potentials. It is therefore important to retain the time-ordering propagator, T, in the definition of the operator U. The exponential in expression (4.52) is clearly too complex to be applied directly, and needs to be approximated in some suitable maimer. To reduce the error in the propagation from to to tf, this large interval is usually split into smaller sub-intervals of length A t. The wave-functions are then propagated from to to + At, then from to At to + 2At and so on. 

The simplest approximation to (4.52) is a direct expansion of the exponential in a power series of At 

Unfortunately, the expression (4.53) does not retain one of the most important properties of the Kohn-Sham time-evolution operator unitarity. In other words, if we apply (4.53) to a normalized wave-function the result will no longer be normalized. This leads to an inherently unstable propagation. 

---

Cowen, 1. R., I. P. Castro, and. A. G. Robins. 1997. Numerical considerations for sinuilations of flow and di.spersion around buildings. In 2nd Conference on Computational Wind Engineering, Colorado. 

N. De Leon, M. A. Mehta, andR. Q. Topper, Cylindrical manifolds in phase space as mediators of chemical reaction dynamics and kinetics. II. Numerical considerations and applications to models with two degrees of freedom, J. Chem. Phys. 94, 8329 (1991). 

Certain additional numerical considerations should be satisfied before a spawning attempt is successful. First, in order to avoid unnecessary basis set expansion, we require that the parent of a spawned basis function have a population greater than or equal to Fmln, where the population of the ktU basis function on electronic state / is defined as... 

Over the course of a reaction model, a mineral may dissolve away completely or become supersaturated and precipitate. In either case, the modeling software must alter the basis to match the new mineral assemblage before continuing the calculation. Finally, the basis sometimes must be changed in response to numerical considerations (e.g., Coudrain-Ribstein and Jamet, 1989). Depending on the... 

Steward (S3) proposed an algorithm based on tearing a variable from only one equation at a time and evaluating each tear on the basis of the size of the resulting subsystems of simultaneous equations in the torn system and numerical considerations of the particular equations. Each variable is torn successively from each equation in which it appears and the effectiveness of the tear evaluated. 

However, this limit has some trouble as mentioned in Ref. [30]. For numerical consideration, we do not take this limit. Only the behavior for small is needed. 

For routine drug-excipient interaction studies, DSC and an isothermal stress test are normally employed. According to van Dooren (274), DSC curves are difficult to evaluate and positive conclusions are rarely obtained thus, the latter test is still necessary. Van Dooren claims that numerous considerations must be made before a positive identification of drug-excipient interaction by DSC can be made. 

The development process converts the latent image in the polymer into the final 3-D relief image. This process is perhaps the most complex of resist technology. It can generally be achieved by either liquid development or dry (plasma) development. Numerous considerations are critical to either alternative. We will first focus on the wet development process. Plasma development will be discussed in a later section. 

Numerical considerations dealing with the total number of individual components in a chemical sense, e,g. binary (b), ternary (t), etc. and with the number of particles [host, guest separately, e.g. monomolecular (Im), binuclear (2n), respectively] are also practicable Thus, the clathrate formed between Dianin s compound (5) and chloroform (see Fig. 4) is identified as binary, hexamolecular, and mononuclear. The full description of this inclusion compound applying the complete set of symbols and designations explained above, hence follows as b, 6m, In-cryptato-clathrate . 

---