Growing Hamiltonian

Consider a growing Hamiltonian as in Figure 4 and the following equation, where the superscript k indicates the stage of growth ... 

For simple calculations such as the one given above, a linear parametrization of the Hamiltonian is sufficient. However in most cases, one is interested in growing new atoms or removing atoms. This is the case if for example one is transforming a glycine residue into an alanine. See Fig. 4.14. 

The largest arrays which occur in calculations are of two types. One arises from the electron repulsion integrals and grows in size like the fourth power of the number of basis functions. The other is the configuration interaction hamiltonian matrix which grows like the square of the number of configurations. Many other smaller arrays, whose size is proportional to the square of the number of basis functions, occur throughout the calculation. 

In his interesting paper Professor Nicolis raises the question whether models can be envisioned which lead to a spontaneous spatial symmetry breaking in a chemical system, leading, for example, to the production of a polymer of definite chirality. It would be even more interesting if such a model would arise as a result of a measure preserving process that could mimic a Hamiltonian flow. Although we do not have such an example of a chiral process, which imbeds an axial vector into the polymer chain, several years ago we came across a stochastic process that appears to imbed a polar vector into a growing infinite chain. 

This is also the Hamiltonian of the activated complex. We will encounter it in Eq. (23) with the customary symbol H. Regardless of its stability properties or the size of the nonlinearity, Eq. (12) is always an invariant manifold. However, we are interested in the case when it is of the saddle type with stable and unstable manifolds. If the physical Hamiltonian is of the form of Eq. (1), then a preliminary, local transformation is not required. The manifold (12) is invariant regardless of the size of the nonlinearity. Moreover, it is also of saddle type with respect to stability in the transverse directions. This can be seen by examining Eq. (1). On qn = Pn = 0 the transverse directions, (i.e., q and p ), are still of saddle type (more precisely, they grow and decay exponentially). 

In order to develop the criterion more quantitatively, consider the sequence of phase-space portraits shown in Figs. 5.4(a) - (d). This sequence suggests that, as the control parameter K increases, the diameter of the resonance islands at Z = 0 mod 27t grows in action. In order to predict the touching point of the resonances, we need the widths of the resonances as a function of K. The width of the resonances is derived on the basis of the Hamiltonian (5.2.1). Since the dynamics induced by H is equivalent to the chaotic mapping (5.1.6), the Hamiltonian H itself cannot be treated analytically and has to be simplified. One way is to consider only the average effect of the periodic 6 kicks in (5.2.1). The average perturbation... 

Application of quantum chemical methods to polymeric systems is complicated. The computer hardware capabilities to deal with real polymers are, at this point, inadequate. In general, the CPU requirements grow with the fourth to sixth power of the number of functions used and the system Hamiltonian does not normally contain information about the solvent. The major limitation is the length of time that can be simulated. Important dynamic polymer properties are associated with relaxation times that are many orders of magni-... 

In addition to the analysis of the topology of a conical intersection, the quadratic expansion of the Hamiltonian matrix can be used as a new practical method to generate a subspace of active coordinates for quantum dynamics calculations. The cost of quantum dynamics simulations grows quickly with the number of nuclear degrees of freedom, and quantum dynamics simulations are often performed within a subspace of active coordinates (see, e.g., [46-50]). In this section we describe a method which enables the a priori selection of these important coordinates for a photochemical reaction. Directions that reduce the adiabatic energy difference are expected to lead faster to the conical intersection seam and will be called photoactive modes . The efficiency of quantum dynamics run in the subspace of these reduced coordinates will be illustrated with the photochemistry of benzene [31,51-53]. 

In order to obtain virtual and resonances states, we have to find eigenfunctions of Hamiltonian (101) which grows up exponentially when x —> oo. Using the fact that the potential goes to zero very fast, we can obtain... 

---