Mean-field solutions computation

In this Section we introduce a stochastic alternative model for surface reactions. As an application we will focus on the formation of NH3 which is described below, equations (9.1.72) to (9.1.76). It is expected that these stochastic systems are well-suited for the description via master equations using the Markovian behaviour of the systems under study. In such a representation an infinite set of master equations for the distribution functions describing the state of the surface and of pairs of surface sites (and so on) arises. As it was told earlier, this set cannot be solved analytically and must be truncated at a certain level. The resulting equations can be solved exactly in a small region and can be connected to a mean-field solution for large distances from a reference point. This procedure is well-suited for the description of surface reaction systems which includes such elementary steps as adsorption, diffusion, reaction and desorption.The numerical part needs only a very small amount of computer time compared to MC or CA simulations. 

This chapter is organized as follows. The thermodynamics of the critical micelle concentration are considered in Section 3.2. Section 3.3 is concerned with a summary of experiments characterizing micellization in block copolymers, and tables are used to provide a summary of some of the studies from the vast literature. Theories for dilute block copolymer solutions are described in Section 3.4, including both scaling models and mean field theories. Computer simulations of block copolymer micelles are discussed in Section 3.5. Micellization of ionic block copolymers is described in Section 3.6. Several methods for the study of dynamics in block copolymer solutions are sketched in Section 3.7. Finally, Section 3.8 is concerned with adsorption of block copolymers at the liquid interface. 

From the rigorous treatment of the double-layer problem on the molecular level, it becomes clear that the Gouy-Chapman theory of the interface is equivalent to a mean field solution of a simple primitive model (PM) of electrolytes at the interface (6). To consider the correlation between ions, integral equations that describe the PM are devised and solved in different approximations. An exact solution of the PM of the electrolyte can be obtained from the computer simulations. This solution can be compared with the solutions obtained from different integral equations. For detailed discussion of this topic, refer to the review by Camie and Torrie (6). In many cases, the molecular description of the solvent must be introduced into the theory to explain the complexity of the observed phenomena. The analytical treatment in such cases is very involved, but initial success has already been achieved. Some of the theoretical developments along these lines were reviewed by Blum (7). 

The minimum value in EKs p(r)] corresponds to the exact ground-state dectron density. To determine the actual energy, variations in E p(r)] must be optimized with respect to variations in p, subject to the orthonormality constraints. In KS-DFT, an arriiidal reference system of noninteracting dectrons is constracted, which has exactly the same electron density as the real molecular system. Therefore, from a computational viewpoint, the KS version of DFT leads to a mean-field solution, which is only an approximation. 

The above theoretical framework was developed in real space. However, computing the correlation functions in real space can be carried out only for simple geometries, such as the lamellar phase [31]. In order to apply the theory to more complex structures, efficient methods other than the direct real-space computation have to be developed. One particularly useful method is the reciprocal-space technique, which utilizes the symmetries of the ordered phases. The key observation is that the mean-field solution w (r) = is a periodic... 

A central issue in statistical thermodynamic modelling is to solve the best model possible for a system with many interacting molecules. If it is essential to include all excluded-volume correlations, i.e. to account for all the possible ways that the molecules in the system instantaneously interact with each other, it is necessary to do computer simulations as discussed above, because there are no exact (analytical) solutions to the many-body problems. The only analytical models that can be solved are of the mean-field type. 

The determination of the ground state energy and the ground state electron density distribution of a many-electron system in a fixed external potential is a problem of major importance in chemistry and physics. For a given Hamiltonian and for specified boundary conditions, it is possible in principle to obtain directly numerical solutions of the Schrodinger equation. Even with current generations of computers, this is not feasible in practice for systems of large total number of electrons. Of course, a variety of alternative methods, such as self-consistent mean field theories, also exist. However, these are approximate. 

For this extended IBM, exact solutions are not possible, and we are forced to use a mean field approximation. The system is cranked to generate the high spin states. This HB approximation is simpler than the corresponding Hartree Fock lowest state all the bosons (of each type) are in a single boson condensate rather than in a set of occupied states as in the HF case. Finding the self consistent solution is a simple numerical problem and requires little computer time. The ground state is of the form ... 

Abstract In this chapter we review recent advances which have been achieved in the theoretical description and understanding of polyelectrolyte solutions. We will discuss an improved density functional approach to go beyond mean-field theory for the cell model and an integral equation approach to describe stiff and flexible polyelectrolytes in good solvents and compare some of the results to computer simulations. Then we review some recent theoretical and numerical advances in the theory of poor solvent polyelectrolytes. At the end we show how to describe annealed polyelectrolytes in the bulk and discuss their adsorption properties. 

It is not necessary to consider spheres only. For any particular system the dielectric-response functions and excess ion densities can be measured or formulated as functions of the size, shape, and charge of the suspended particles as well as from ionic properties of the bathing solution. It is a separate procedure to compute the excess numbers Tv from mean-field theory. They are used here as given quantities. 

The purpose of this chapter is to get a better insight of the first problem listed above, i.e., the polarization of interfaces (colloidal particles) during their interaction. Because of tutorial reasons, the electrolyte solution will be described using a rather simple, mean field approximation, that, however, allows to obtain an analytical solution of the problem. It is clear that this elaboration can easily be followed, and one can extend our model on more advanced situations. This model is identical with so-called weak-coupling theory for point ions treated in a course of the Debye-Hiickel approximation. Before going to make an elaboration for two interacting macrobodies immersed into an electrolyte solution, we would like to introduce a method, which is usually used to model this polarization, and to compute the electrical field next to a polarized medium. Then we will also discuss consequences of the polarization for the ion distribution at the particle-solution interface. 

Despite all this, DFT is credited of many successes and can be very useful, being computationally tractable for complex systems. DFT results with GGA functionals are often of much better quality than Hartree-Fock. Although DFT is often seen as a mean field approach (from its structure and since setting Ex c = 0 amounts to the mean field Hartree solution of the electronic structure problem), approximate functionals like LDA or GGA s incorporate terms that are non mean field contributions. Furthermore they do so while fulfilling some exact relations (sum rules etc.). 

This name covers all polymer chains (diblocks and others) attached by one end (or end-block) at ( external ) solid/liquid, liquid/air or ( internal ) liquid/liq-uid interfaces [226-228]. Usually this is achieved by the modified chain end, which adsorbs to the surface or is chemically bound to it. Double brushes may be also formed, e.g., by the copolymers A-N, when the joints of two blocks are located at a liquid/liquid interface and each of the blocks is immersed in different liquid. A number of theoretical models have dealt specifically with the case of brush layers immersed in polymer melts (and in solutions of homopolymers). These models include scaling approaches [229, 230], simple Flory-type mean field models [230-233], theories solving self-consistent mean field (SCMF) equations analytically [234,235] or numerically [236-238]. Also first computer simulations have recently been reported for brushes immersed in a melt [239]. 

In recent years, there have been many attempts to combine the best of both worlds. Continuum solvent models (reaction field and variations thereof) are very popular now in quantum chemistry but they do not solve all problems, since the environment is treated in a static mean-field approximation. The Car-Parrinello method has found its way into chemistry and it is probably the most rigorous of the methods presently feasible. However, its computational cost allows only the study of systems of a few dozen atoms for periods of a few dozen picoseconds. Semiempirical cluster calculations on chromophores in solvent structures obtained from classical Monte Carlo calculations are discussed in the contribution of Coutinho and Canuto in this volume. In the present article, we describe our attempts with so-called hybrid or quantum-mechanical/molecular-mechanical (QM/MM) methods. These concentrate on the part of the system which is of primary interest (the reactants or the electronically excited solute, say) and treat it by semiempirical quantum chemistry. The rest of the system (solvent, surface, outer part of enzyme) is described by a classical force field. With this, we hope to incorporate the essential influence of the in itself uninteresting environment on the dynamics of the primary system. The approach lacks the rigour of the Car-Parrinello scheme but it allows us to surround a primary system of up to a few dozen atoms by an environment of several ten thousand atoms and run the whole system for several hundred thousand time steps which is equivalent to several hundred picoseconds. 

---