Subsystems

The HILL-SCAN boards contain the following subsystems ... 

For a system composed of two subsystems a and p separated from each other by a diathemiic wall and from the surroundings by adiabatic walls, the equation corresponding to equation (A2.1.12) is... 

It is still necessary to consider the role of entropy m irreversible changes. To do this we return to the system considered earlier in section A2.1.4.2. the one composed of two subsystems in themial contact, each coupled with the outside tliroiigh movable adiabatic walls. Earlier this system was described as a function of tliree independent variables, F , and 0 (or 7). Now, instead of the temperature, the entropy S = +. S P will be... 

Consider two ideal-gas subsystems a and (3 coupled by a movable diatliemiic wall (piston) as shown in figure A2.1.5. The wall is held in place at a fixed position / by a stop (pin) that can be removed then the wall is free to move to a new position / . The total system (a -t P) is adiabatically enclosed, indeed isolated q = w = 0), so the total energy, volume and number of moles are fixed. 

Consider the situation illustrated in figure A2.1.5. with the modifieation that the piston is now an adiabatie wall, so the two temperatures need not be equal. Energy is transmitted from subsystem a to subsystem (3 only in the fomi of work obviously dF = -dF so, in applying equation (A2.1.20), is dlf- P equal to dF = dF or equal todk , or is it something else entirely One ean measure the ehanges in temperature,... 

The paradox involved here ean be made more understandable by introdueing the eoneept of entropy ereation. Unlike the energy, the volume or the number of moles, the entropy is not eonserved. The entropy of a system (in the example, subsystems a or P) may ehange in two ways first, by the transport of entropy aeross the boundary (in this ease, from a to P or vice versa) when energy is transferred in the fomi of heat, and seeond. 

The total change d.S can be detennined, as has been seen, by driving the subsystem a back to its initial state, but the separation into dj.S and dj S is sometimes ambiguous. Any statistical mechanical interpretation of the second law requires that, at least for any volume element of macroscopic size, dj.S > 0. However, the total... 

The assumption (frequently unstated) underlying equations (A2.1.19) and equation (A2.1.20) for the measurement of irreversible work and heat is this in the surroundings, which will be called subsystem p, internal equilibrium (unifomi T, p and //f diroughout the subsystem i.e. no temperature, pressure or concentration gradients) is maintained tliroughout the period of time in which the irreversible changes are... 

Subsystem p may now be called the surroundings or as Callen (see further reading at the end of this article) does, in an excellent discussion of this problem, a source . To fomudate this mathematically one notes that, if dj.S P = 0, one can then write... 

Two subsystems a. and p, in each of which the potentials T,p, and all the p-s are unifonn, are pennitted to interact and come to equilibrium. At equilibrium all infinitesimal processes are reversible, so for the overall system (a + P), which may be regarded as isolated, the quantities conserved include not only energy, volume and numbers of moles, but also entropy, i.e. there is no entropy creation in a system at equilibrium. One now... 

If there are more than two subsystems in equilibrium in the large isolated system, the transfers of S, V and n. between any pair can be chosen arbitrarily so it follows that at equilibrium all the subsystems must have the same temperature, pressure and chemical potentials. The subsystems can be chosen as very small volume elements, so it is evident that the criterion of internal equilibrium within a system (asserted earlier, but without proof) is unifonnity of temperature, pressure and chemical potentials tlu-oughout. It has now been... 

Thus, the neglect of the off-diagonal matrix elements allows the change from mixed states of the nuclear subsystem to pure ones. The motion of the nuclei leads only to the deformation of the electronic distribution and not to transitions between different electronic states. In other words, a stationary distribution of electrons is obtained for each instantaneous position of the nuclei, that is, the elechons follow the motion of the nuclei adiabatically. The distribution of the nuclei is described by the wave function x (R i) in the potential V + Cn , known as the proper adiabatic approximation [41]. The off-diagonal operators C n in the matrix C, which lead to transitions between the states v / and t / are called operators of nonadiabaticity and the potential V = (R) due to the mean field of all the electrons of the system is called the adiabatic potential. 

Let us examine a special but more practical case where the total molecular Hamiltonian, H, can be separated to an electronic part, W,.(r,s Ro), as is the case in the usual BO approximation. Consequendy, the total molecular wave function fl(R, i,r,s) is given by the product of a nuclear wave function X uc(R, i) and an electronic wave function v / (r, s Ro). We may then talk separately about the permutational properties of the subsystem consisting of electrons, and the subsystemfs) formed of identical nuclei. Thus, the following commutative laws Pe,Hg =0 and =0 must be satisfied X =... 

As pointed out in the previous paragraph, the total wave function of a molecule consists of an electronic and a nuclear parts. The electrons have a different intrinsic nature from nuclei, and hence can be treated separately when one considers the issue of permutational symmetry. First, let us consider the case of electrons. These are fermions with spin and hence the subsystem of electrons obeys the Fermi-Dirac statistics the total electronic wave function... 

Let us discuss further the pemrutational symmetry properties of the nuclei subsystem. Since the elechonic spatial wave function t / (r,s Ro) depends parameti ically on the nuclear coordinates, and the electronic spacial and spin coordinates are defined in the BF, it follows that one must take into account the effects of the nuclei under the permutations of the identical nuclei. Of course. 

According to Section VI, the size M of the sub-Hilbert space is determined whether the respective M states form an isolated set of states fulfilling Eqs. (91). In this case, diabatization is always valid for this subsystem. However, it can happen that under certain geometrical situations the size of the sub-Hilbert space for which diabatization is valid is even smaller than this particular M... 

As ab initio MD for all valence electrons [27] is not feasible for very large systems, QM calculations of an embedded quantum subsystem axe required. Since reviews of the various approaches that rely on the Born-Oppenheimer approximation and that are now in use or in development, are available (see Field [87], Merz ]88], Aqvist and Warshel [89], and Bakowies and Thiel [90] and references therein), only some summarizing opinions will be given here. 

Hence, as the second class of techniques, we discuss adaptive methods for accurate short-term integration (Sec. 4). For this class, it is the major requirement that the discretization allows for the stepsize to adapt to the classical motion and the coupling between the classical and the quantum mechanical subsystem. This means, that we are interested in discretization schemes which avoid stepsize restrictions due to the fast oscillations in the quantum part. We can meet this requirement by applying techniques recently developed for evaluating matrix exponentials iteratively [12]. This approach yields an adaptive Verlet-based exponential integrator for QCMD. 

To begin with, we compare the stepsizes used in the simulations (Fig. 3). As pointed out before, it seems to be unreasonable to equip the Pickaback scheme with a stepsize control, because, as we indeed observe in Fig. 3, the stepsize never increases above a given level. This level depends solely on the eigenvalues of the quantum Hamiltonian. When analyzing the other integrators, we observe that the stepsize control just adapts to the dynamical behavior of the classical subsystem. The internal (quantal) dynamics of the Hydrogen-Chlorine subsystem does not lead to stepsize reductions. 

In the sequel, we assume that the quantum subsystem has been truncated to a finite-dimensional system by an appropriate spatial discretization and a corresponding representation of the wave function by a complex-valued vector Ip C. The discretized quantum operators T, V and H are denoted by T e V(q) E and H q) e respectively. In the following... 

Here we suggest a different approach that propagates the system using multiple step-sizes, i.e., few steps with step-size At are taken in the slow classical part whereas many smaller steps with step-size 5t are taken in the highly oscillatory quantum subsystem (see, for example, [19, 4] for symplectic multiple-time-stepping methods in the context of classical molecular dynamics). Therefore, we consider a splitting of the Hamiltonian H = Hi +H2 in the following way ... 

NAMD was implemented in an object-oriented fashion (Fig. 3). Patches, the encapsulated communication subsystem, the molecular structure, and various output methods were objects. Every patch owned specialized objects... 

Fig. 3. NAMD 1 employs a modular, object-oriented design in which patches communicate via an encapsulated communication subsystem. Every patch owns an integrator and a complete set of force objects for bonded (BondForce), nonbonded (ElectForce), and full electrostatic (DPMTA) calculations.

The classical architecture of an expert system comprises a knowledge base, an inference engine, and some kind of user interface. Most expert systems also include an explanation subsystem and a knowledge acquisition subsystem. This architecture is given in Figure 9-34 and described in more detail below. 

Inference engine The inference engine represents the central problem-solving subsystem. It contains strategies for using the information contained in the... 

Besides these three major components, many expert systems also comprise an explanation subsystem and a knowledge acquisition subsystem,... 

Explanation mhsystenr. The explanation subsystem supplies information about the reasoning process, for example how the conclusion has been drawn and which facts were used to come to the conclusion,... 

Knowledge acquisition subsystem The task of the knowledge acquisition subsystems is to assemble and upgrade the knowledge base. A major task is to verify the data and check for consistency. 

Their symmetry labels can be obtained by vector coupling (see Appendix G) the spin and orbital angular momenta of the two subsystems. The orbital angular momentum coupling... 

---