Momentum space scheme

In this section, we briefly summarize the major theoretical techniques employed in contemporary pseudopotential calculations. The approach involves the use of initio pseudopotentials3 7 to compute the electron-core interactions and local density function-als° for evaluating electron-electron interactions. Several different basis sets can be used to solve the electron wave equation. The particular choice is determined by the system of interest. Finally, the total energies and forces are calculated using a momentum space scheme. 

Molecular dynamics has also been used to replace the MC moves for conformational advancement [43]. In the molecular dynamics version of parallel tempering, often referred to as replica exchange molecular dynamics, momenta are used in the propagation scheme such that a constant temperature is maintained between the swaps. After the swap in conformational space (with the same acceptance criterion as in the MC implementation), a readjustment in momentum space is also needed. This is done by renewing the momenta for replica i by the transformation... 

In the case of semiconductors, the idea of electron tunneling has been used by Zener [42] to describe the so-called interband tunneling. Such tunneling represents one of the possible mechanisms of semiconductor breakdown. To understand the nature of interband tunneling, we shall first follow Ziman [43] and consider the one-dimensional motion of an electron in a separate band under the influence of an electric field. If we use the scheme of repeated bands, then the electron motion in momentum space is an up and down motion along the OABC periodic curve (Fig. 16). In the coordinate space, the electron, starting from point O, accelerates then slows down as it approaches point A here, the direction of the motion is changed to the... 

According to (5.35), the most fortunate circumstance for the present scheme is a system with heavy mass and parallel potential energy surfaces (A(x) const.). The steepness of the potential difference A(x) is the most crucial parameter it not only affects the validity of this level approximation (5.8) but it also changes the efficiency of excitation according to (5.32). It is obvious that a narrow wavepacket can be relatively easily excited by a quadratically chirped pulse (cf. 5.32). However, a narrow one can easily break the level approximation (5.35) because of the broad distribution in momentum space. The optimal width of a wavepacket can be roughly estimated as... 

The decomposition of the irreducible part of the self-energy wave-function correction term is depicted in Fig. 2. The divergent terms are these with zero and one interaction in the binding potential present, below referred to as zero-potential term and one-potential term , respectively. The charge divergences cancel between both terms. In addition, a mass connter term Sm has to be subtracted to obtain proper mass renormalization similar to the case of the free self energy [47] (for onr schemes see also [44]). The zero- and one-potential terms are then semianalytically evaluated in momentum space (for details cf. 

One simple scheme to effect this interpolation is as follows. Consider a wave function which is sampled by N points. It is first transformed to momentum space. Then the wave function is cast onto a larger grid of M points by adding M — N zeros to the momentum values for fc > ttN/L. A back transform will increase the density of points without adding any new information to the wave function. 

Other algorithms have been developed for speeding up the uniform sampling of phase space points. For example, the efficient microcanonical sampling series of schemes exploits the possibility of sampling independently the spatial coordinates and momenta, simply by weighting the sampled geometries by their associated momentum space density. 

Most schemes that have been proposed to propel starships involve plasmas. Schemes differ both in the selection of matter for propulsion and the way it is energi2ed for ejection. Some proposals involve onboard storage of mass to be ejected, as in modem rockets, and others consider acquisition of matter from space or the picking up of pellets, and their momentum, which are accelerated from within the solar system (184,185). Energy acquisition from earth-based lasers also has been considered, but most interstellar propulsion ideas involve nuclear fusion energy both magnetic, ie, mirror and toroidal, and inertial, ie, laser and ion-beam, fusion schemes have been considered (186—190). 

It is easy to see that this expression has two minima within the Brillouin zone. One minimum is at fc = 0 and gives the correct continuum limit. The other, however, is at k = 7t/a and carries an infinite momentum as the lattice spacing a 0. In other words, discretizing the fermion field leads to the unphysical problem of species doubling. (In fact, since there is a doubling for each space-time dimension, this scheme actually results in 2 = 16 times the expected number of fermions.)... 

Multiparticle collision dynamics describes the interactions in a many-body system in terms of effective collisions that occur at discrete time intervals. Although the dynamics is a simplified representation of real dynamics, it conserves mass, momentum, and energy and preserves phase space volumes. Consequently, it retains many of the basic characteristics of classical Newtonian dynamics. The statistical mechanical basis of multiparticle collision dynamics is well established. Starting with the specification of the dynamics and the collision model, one may verify its dynamical properties, derive macroscopic laws, and, perhaps most importantly, obtain expressions for the transport coefficients. These features distinguish MPC dynamics from a number of other mesoscopic schemes. In order to describe solute motion in solution, MPC dynamics may be combined with molecular dynamics to construct hybrid schemes that can be used to explore a variety of phenomena. The fact that hydrodynamic interactions are properly accounted for in hybrid MPC-MD dynamics makes it a useful tool for the investigation of polymer and colloid dynamics. Since it is a particle-based scheme it incorporates fluctuations so that the reactive and nonreactive dynamics in small systems where such effects are important can be studied. 

As a consequence of the transformation, the equation of motion depends on three extra coordinates which describe the orientation in space of the rotating local system. Furthermore, there are additional terms in the Hamiltonian which represent uncoupled momenta of the nuclear and electronic motion and moment of inertia of the molecule. In general, the Hamiltonian has a structure which allows for separation of electronic and vibrational motions. The separation of rotations however is not obvious. Following the standard scheme of the various contributions to the energy, one may assume that the momentum and angular momentum of internal motions vanish. Thus, the Hamiltonian is simplified to the following form. 

The Verlet scheme propagates the position vector with no reference to the particle velocities. Thus, it is particularly advantageous when the position coordinates of phase space are of more interest than the momentum coordinates, e.g., when one is interested in some property that is independent of momentum. However, often one wants to control the simulation temperature. This can be accomplished by scaling the particle velocities so that the temperature, as defined by Eq. (3.18), remains constant (or changes in some defined manner), as described in more detail in Section 3.6.3. To propagate the position and velocity vectors in a coupled fashion, a modification of Verlet s approach called the leapfrog algorithm has been proposed. In this case, Taylor expansions of the position vector truncated at second order... 

With the introduction of electronic angular momentum, we have to consider how the spin might be coupled to the rotational motion of the molecule. This question becomes even more important when electronic orbital angular momentum is involved. The various coupling schemes give rise to what are known as Hund s coupling cases they are discussed in detail in chapter 6, and many practical examples will be encountered elsewhere in this book. If only electron spin is involved, the important question is whether it is quantised in a space-fixed axis system, or molecule-fixed. In this section we confine ourselves to space quantisation, which corresponds to Hund s case (b). 

---