Particle finite, systems

In finite boundary conditions the solute molecule is surrounded by a finite layer of explicit solvent. The missing bulk solvent is modeled by some form of boundary potential at the vacuum/solvent interface. A host of such potentials have been proposed, from the simple spherical half-harmonic potential, which models a hydrophobic container [22], to stochastic boundary conditions [23], which surround the finite system with shells of particles obeying simplified dynamics, and finally to the Beglov and Roux spherical solvent boundary potential [24], which approximates the exact potential of mean force due to the bulk solvent by a superposition of physically motivated tenns. 

In his early survey of computer experiments in materials science , Beeler (1970), in the book chapter already cited, divides such experiments into four categories. One is the Monte Carlo approach. The second is the dynamic approach (today usually named molecular dynamics), in which a finite system of N particles (usually atoms) is treated by setting up 3A equations of motion which are coupled through an assumed two-body potential, and the set of 3A differential equations is then solved numerically on a computer to give the space trajectories and velocities of all particles as function of successive time steps. The third is what Beeler called the variational approach, used to establish equilibrium configurations of atoms in (for instance) a crystal dislocation and also to establish what happens to the atoms when the defect moves each atom is moved in turn, one at a time, in a self-consistent iterative process, until the total energy of the system is minimised. The fourth category of computer experiment is what Beeler called a pattern development... 

The motion of particles of the film and substrate were calculated by standard molecular dynamics techniques. In the simulations discussed here, our purpose is to calculate equilibrium or metastable configurations of the system at zero Kelvin. For this purpose, we have applied random and dissipative forces to the particles. Finite random forces provide the thermal motion which allows the system to explore different configurations, and the dissipation serves to stabilize the system at a fixed temperature. The potential energy minima are populated by reducing the random forces to zero, thus permitting the dissipation to absorb the kinetic energy. 

Fig. 2.2. Average electrostatic potential mc at the position of the methane-like Lennard-Jones particle Me as a function of its charge q. mc contains corrections for the finite system size. Results are shown from Monte Carlo simulations using Ewald summation with N = 256 (plus) and N = 128 (cross) as well as GRF calculations with N = 256 water molecules (square). Statistical errors are smaller than the size of the symbols. Also included are linear tits to the data with q < 0 and q > 0 (solid lines). The fit to the tanh-weighted model of two Gaussian distributions is shown with a dashed line. Reproduced with permission of the American Chemical Society...

The above observation suggests an intriguing relationship between a bulk property of infinite nuclear matter and a surface property of finite systems. Here we want to point out that this correlation can be understood naturally in terms of the Landau-Migdal approach. To this end we consider a simple mean-field model (see, e.g., ref.[16]) with the Hamiltonian consisting of the single-particle mean field part Hq and the residual particle-hole interaction Hph-... 

The investigation of the behaviour of Z(r, t) in a finite system is a difficult problem due to the fact that the boundary conditions for (7.3.3) are unknown for fractal systems although one can use some scaling arguments based on the knowledge of the properties of the infinite system. In the case of the independent production of different particles we obtain... 

Let us consider two parallel plates immersed in a solution of a 1 1 electrolyte and of small finite size particles which can dissociate monovalent counterions, thus acquiring a charge of ze per particle. The system contains two ions of the electrolyte which are assumed to be of negligible size, one kind of charged particles which have a finite size, and their counterions which are of negligible size. Assuming Boltzmann distributions for the ions and charged particles, eq 2 becomes... 

Classical thermodynamics is based on a description of matter through such macroscopic properties as temperature and pressure. However, these properties are manifestations of the behavior of the countless microscopic particles, such as molecules, that make up a finite system. Evidently, one must seek an understanding of the fundamental nature of entropy in a microscopic description of matter. Because of the enormous number of particles contained in any system of interest, such a description must necessarily be statistical in nature. We present here a very brief indication of the statistical interpretation of entropy, t... 

The advantage of the ROPM lies in the fact that due to the Fock form of [n] the self-interaction of the KS-particles is cancelled exactly. This also manifests itself in the asymptotic form of for finite systems,... 

In this section we are going to develop a different approach to the calculation of excitation energies which is based on TDDFT [69, 84, 152]. Similar ideas were recently proposed by Casida [223] on the basis of the one-particle density matrix. To extract excitation energies from TDDFT we exploit the fact that the frequency-dependent linear density response of a finite system has discrete poles at the excitation energies of the unperturbed system. The idea is to use the formally exact representation (156) of the linear density response n j (r, cu), to calculate the shift of the Kohn-Sham orbital energy differences coj (which are the poles of the Kohn-Sham response function) towards the true excitation energies Sl in a systematic fashion. 

Phase transitions in statistical mechanical calculations arise only in the thermodynamic limit, in which the volume of the system and the number of particles go to infinity with fixed density. Only in this limit the free energy, or any thermodynamic quantity, is a singular function of the temperature or external fields. However, real experimental systems are finite and certainly exhibit phase transitions marked by apparently singular thermodynamic quantities. Finite-size scaling (FSS), which was formulated by Fisher [22] in 1971 and further developed by a number of authors (see Refs. 23-25 and references therein), has been used in order to extrapolate the information available from a finite system to the thermodynamic limit. Finite-size scaling in classical statistical mechanics has been reviewed in a number of excellent review chapters [22-24] and is not the subject of this review chapter. 

In statistical mechanics, the existence of phase transitions is associated with singularities of the free energy per particle in some region of the thermodynamic space. These singularities occur only in the thermodynamic limit [68,69] in this limit the volume (V) and particle number (TV) go to infinity in such a way that the density (p = N/V) stays constant. This fact could be understood by examining the partition function. For a finite system, the partition function is a finite sum of analytical terms, and therefore it is itself an analytical function. It is necessary to take an infinite number of terms in order to obtain a singularity in the thermodynamic limit [68,69]. 

If fl = 1, every atom in the slider has the same velocity at every instant of time, once steady state (not necessarily smooth sliding) has been reached. Hence the problem is reduced to the motion of a single particle, for which one obtains Fj = 1. This provides an upper bound of Fj for arbitrary a. If the walls are incommensurate or disordered, one can again make use of the argument that the motion of all atoms relative to their preferred positions is the same up to temporal shifts once steady state has been reached. Owing to the incommensurability, the distribution of these temporal shifts with respect to a reference trajectory cannot change with time in the thermodynamic limit, and the instantaneous value of Fk is identical to Fk at all times. This gives a lower bound for Fj for arbitrary a. The static friction for arbitrary commensurability and/or finite systems lies in between the upper and the lower bound. 

An important difference between quantum and classical mechanics is that in classical mechanics stationary states exist at all energies while in quantum mechanics of finite systems the spectrum is discrete as shown in the examples above. This difference disappears when the system becomes large. Even for a single particle system, the spacings between allowed energy levels become increasingly smaller as the size of accessible spatial extent of the system increases, as seen, for example. 

Still, as discussed in Section 2,8.1, normalization is in some sense still a useful concept even for such processes. As we saw in Section 2.8.1, we may think of an infinite system as a Q oo limit of a finite system of volume Q. Intuition suggests that a scattering process characterized by a short range potential should not depend on system size. On the other hand the normalization condition dx fl/(x)p = 1 implies that scattering wavefunctions will vanish everywhere like as Q CX3. We have noted (Section 2.8) that physically meaningful results are associated either with products such as A V<(x) 2 or yo V (x)p, where jV, the total number of particles, and p, the density of states, are both proportional to Q. Thus, for physical observables the volume factor cancels. 

Thus Fxc(rs,. s) is a measure of the enhancement in the energy per particle over local exchange. In the rest of this subsection, we plot curves of Fxc(r,s) for different approximations for several values of The energies of real systems can contain significant contributions from s up to about 3, and r, up to about 18. Valence electrons in solid metals have s < 2 and 1 rs < 6. In the core of an atom, s < 1 and rs 1. In the limit r —> oo for a finite system, rs and s grow exponentially. 

An outstanding problem concerns itself with the structure of a hard sphere phase. This is a special instance of the more. general difficulty of the specification of the structure of infinitely extended random media. These questions will perhaps be the subject of a future mathematical discipline-stochastic geometry. The pair correlation function g(r), even if it is known, hardly suffices to specify uniqudy the stochastic metric properties of a random structure. For a finite N and V finite) system in equilibrium in thermal contact with a heat reservoir at temperature T, the density in the configuration space of the N particles [Eq. (2)]... 

---