Hamiltonian systems overlap

In the recent past, analytical research in Celestial Mechanics has centred on KAM theory and its applications to the dynamics of low dimensional Hamiltonian systems. Results were used to interpret observed solutions to three body problems. Order was expected and chaos or disorder the exception. Researchers turned to the curious exception, designing analytical models to study the chaotic behaviour at resonances and the effects of resonant overlaps. Numerical simulations were completed with ever longer integration times, in attempts to explore the manifestations of chaos. These methods improved our understanding but left much unexplained phenomena. 

Various numerical methods have been proposed for collisional Hamiltonian systems [136, 176, 184, 263, 348, 352]. Typically, these schemes rely on the Verlet method to propagate the system between collisions, with collisions detected either (i) by checking for overlap at the end of the step, (ii) checking for overlap during the step, or (iii) approximating the time to collision before the step. Collisions lead to momentum exchange between particles according the principle outlined above. 

The system that we are interested in is divided into three partitions a molecular part (M) and left/right electrodes (L/R). The total Hamiltonian (//) and overlap (5) matrix within the localized atomic orbital basis representation are written as follows ... 

With the use of symmetries, the Hamiltonian and overlap matrices can be reduced to a minimum number of independent parameters [20] that are calculated once forever for a prototype system or averaging over a wide class of molecules. 

For bound state systems, eigenfunctions of the nuclear Hamiltonian can be found by diagonalization of the Hamiltonian matiix in Eq. (11). These functions are the possible nuclear states of the system, that is, the vibrational states. If these states are used as a basis set, the wave function after excitation is a superposition of these vibrational states, with expansion coefficients given by the Frank-Condon overlaps. In this picture, the dynamics in Figure 4 can be described by the time evolution of these expansion coefficients, a simple phase factor. The periodic motion in coordinate space is thus related to a discrete spectrum in energy space. 

The systems discussed in this chapter give some examples using different theoretical models for the interpretation of, primarily, UPS valence band data, both for pristine and doped systems as well as for the initial stages of interface formation between metals and conjugated systems. Among the various methods used in the examples are the following semiempirical Hartree-Fock methods such as the Modified Neglect of Diatomic Overlap (MNDO) [31, 32) and Austin Model 1 (AMI) [33] the non-empirical Valence Effective Hamiltonian (VEH) pseudopotential method [3, 34J and ab initio Hartree-Fock techniques. 

The principles involved are illustrated well by a system consisting of two particles. In the first instance it may be assumed that the particles are so weakly interacting that their respective wave functions do not overlap significantly. The Hamiltonian may then be assumed to consist of two parts... 

Fig. 2. The quantum mechanics of the two-state prpblem provide a paradigm for the much more extensive electronic state space of a real molecular or macromolecular system. The eigenvectors c, of the Hamiltonian are symmetric and antisymmetric linear combinations of the localized basis vectors with an eigenvalue splitting of 2A, where s is the overlap integral and A is the direct coupling (the only kind possible in this case)...

So far we have assumed that the electronic structure of the crystal consists of one band derived, in our approximation, from a single atomic state. In general, this will not be a realistic picture. The metals, for example, have a complicated system of overlapping bands derived, in our approximation, from several atomic states. This means that more than one atomic orbital has to be associated with each crystal atom. When this is done, it turns out that even the equations for the one-dimensional crystal cannot be solved directly. However, the mathematical technique developed by Baldock (2) and Koster and Slater (S) can be applied (8) and a formal solution obtained. Even so, the question of the existence of otherwise of surface states in real crystals is diflBcult to answer from theoretical considerations. For the simplest metals, i.e., the alkali metals, for which a one-band model is a fair approximation, the problem is still difficult. The nature of the difficulty can be seen within the framework of our simple model. In the first place, the effective one-electron Hamiltonian operator is really different for each electron. If we overlook this complication and use some sort of mean value for this operator, the operator still contains terms representing the interaction of the considered electron with all other electrons in the crystal. The Coulomb part of this interaction acts in such a way as to reduce the effect of the perturbation introduced by the existence of a free surface. A self-consistent calculation is therefore essential, and the various parameters in our theory would have to be chosen in conformity with the results of such a calculation. 

Equation (4.150) expresses the score P as an overlap between the gradient P and the change of system state Ap. In order to find expressions for P in terms of physically insightful quantities, we decompose the total Hamiltonian into system, bath, and interaction parts. 

Alternatively, Eq. (4.160) can be written as an overlap of system and bath matrices, which allows a more direct physical interpretation. To do so, we assume /-dimensional Hilbert space and expand the interaction Hamiltonian as a sum of products of system and bath operators. 

In the perturbative "transfer Hamiltonian approach developed by Bardeen 58), the tip and sample are treated as two non-interacting subsystems. Instead of trying to solve the problem of the combined system, each separate component is described by its wave function, i tip and i/zj, respectively. The tunneling current is then calculated by considering the overlap of these in the tunnel junction. This approach has the advantage that the solutions can be found, for many practical systems, at least approximately, by solution of the stationary Schrodinger equation. 

Show that if the overlap between torsional-vibration wave functions corresponding to oscillation about different equilibrium configurations is neglected, the perturbation-theory secular equation (1.207) for internal rotation in ethane has the same form as the secular equation for the Hiickel MOs of the cyclopropenyl system, thereby justifying (5.96)-(5.98). Write down an expression (in terms of the Hamiltonian and the wave functions) for the energy splitting between sublevels of each torsional level. 

Even at high n s one needs to follow the system for many orbital periods if one is to mimic the experimental results. The difficulty is compounded if one measures the time in units of periods of the core motion. This suggests that the time evolution be characterized using the stationary states of the Hamiltonian rather than propagating the initial state. We have done so, but our experience is that in the presence of DC fields of experimental magnitude (which means that Stark manifolds of adjacent n values overlap), and certainly so in the presence of other ions that break the cylindrical symmetry and hence mix the m/ values, the size of the basis required for convergence is near the limit of current computers. In our experience, truncating the quan-... 

---