Conformational probability distribution

It is clear from Figures 4 and 5 that the chain dimensions and probabilities have essentially reached their asymptotic values at modest temperatures. We should expect, therefore, that in melts of semiflexible PLCs, unless there are extreme extenuating circumstances (i. e., perfect core orientational ordering), the dimensions of a spacer chain and the conformer probability distributions within a spacer chain should be similar to those quantities in an alkane liquid. The attachment of a spacer at both ends will of course perturb the dimensions and conformer probabilities of the spacer. Clearly the magnitude of this perturbation of the spacer is the key to understanding the role of the spacer in a PLC. 

The basic idea of these theories is to look at the distribution of conformations of a chain molecule attached to the surface. The conformational probability distribution function is written in terms of the non-local interaction field induced by the other chain molecules. This field is anisotropic, i.e., it depends on the direction perpendicular to the surface, because the presence of the surface and the inhomogeneous variation of the density of polymer segments and solvent molecules as a function of the distance from the surface. The non-local mean-field is determined by packing constraints that take into account the fact that the volume (at all distances from the surface) must be filled by polymer segments or solvent molecules. These self-consistent criteria represent the incompressibility assumption at all distances from the surface. 

Application of Eq. (30) corrects the free energies of the endpoints but not those of the intermediate conformations. Therefore, the above approach yields a free energy profile between qp and q-g, that is altered by the restraint(s). In particular, the barrier height is not that of the namral, unrestrained system. It is possible to correct the probability distributions P,. observed all along the pathway (with restraints) to obtain those of the unrestrained system [8,40]. Erom the relation P(q)Z, = P,(q)Z, cxp(UJkT) and Eqs. (6)-(8), one obtains... 

The portion of the polymer consisting of molecules terminated by transfer will conform to the most probable distribution, its average degree of polymerization being... 

At the gel point, (3 —l) = l/p, which with the foregoing expression gives Eq. (14), thus establishing equivalence of the two procedures. The primary molecules in a condensation polymer must almost invariably conform to a most probable distribution (see Chap. VIII). The random cross-linking of primary molecules otherwise distributed in size has no counterpart in polyfunctional condensation, therefore. 

Here p (I) is the distribution of conformational states that arises from a simulation using the biased potential. The tricky point with this method comes from the fact that we ultimately need to integrate the work function over a series of windows, and the integration constant for each window is undefined. In practice, this problem is addressed using clever approaches that attempt to match up the probability distributions on consecutive intervals. 

The next step is to generate all possible and allowed conformations, which leads to the full probability distribution F). The normalisation of this distribution gives the number of molecules of type i in conformation c, and from this it is trivial to extract the volume fraction profiles for all the molecules in the system. With these density distributions, one can subsequently compute the distribution of charges in the system. The charges should be consistent with the electrostatic potentials, according to the Poisson equation ... 

The probability distribution of isomeric conformations in PDMS is investigated by both conformational energy considerations and by molecular dynamics simulations. A comparatively smooth distribution of isomeric states is obtained from both approaches. A new RIS treatment, compatible with the molecular mechanics and dynamics considerations, is introduced for describing the conformational statistics of PDMS. 

It should be noted that in the limit where Tw is very short in comparison to the time scale of conformational dynamics, R remains fixed during the time interval Tw and E(Tw) therefore reduces to its instantaneous value,E = 1 + [i /i o]6 -1 In such a case, there is a direct and exact relationship between the probability distribution of E and P(R) ... 

Fig. 2.14 (a) Molecular dynamic simulation results of the variation of the free energy change (A G/R) of the chair-to-chair conformational change against of the number of carbon atoms of the spacer group, (b) Probability distribution of the torsional angle obtained at 1000 K for PCHMA... 

In principle, once the probability distribution function is available, bulk solution properties can be evaluated by averaging appropriate functions of conformation space and time. From the Kuhn and Grun analysis leading to equation (7.24) for the refractive index tensor, we are particularly interested in the second moment tensor,... 

Using 112 nodes of the Earth Simulator, we performed a REMD simulation of this system with 224 replicas. The REMD simulation was successful in the sense that we observed a random walk in potential energy space, which suggests that a wide conformational space was sampled. In Fig. 4.10 we show the canonical probability distributions of the total potential energy at the corresponding 224 temperatures ranging from 250 to 700 K. 

First, the current state of affairs is remarkably similar to that of the field of computational molecular dynamics 40 years ago. While the basic equations are known in principle (as we shall see), the large number of unknown parameters makes realistic simulations essentially impossible. The parameters in molecular dynamics represent the force field to which Newton s equation is applied the parameters in the CME are the rate constants. (Accepted sets of parameters for molecular dynamics are based on many years of continuous development and checking predictions with experimental measurements.) In current applications molecular dynamics is used to identify functional conformational states of macromolecules, i.e., free energy minima, from the entire ensemble of possible molecular structures. Similarly, one of the important goals of analyzing the CME is to identify functional states of areaction network from the entire ensemble of potential concentration states. These functional states are associated with the maxima in the steady state probability distribution function p(n i, no, , hn). In both the cases of molecular dynamics and the CME applied to non-trivial systems it is rarely feasible to enumerate all possible states to choose the most probable. Instead, simulations are used to intelligently and realistically sample the state space. 

Knowledge of the conformations of molecules is important for making structure-activity comparison and in developing a drug-receptor site theory. The preferred conformation may indicate a lock and key arrangement or a probable distribution of charge on the surface of the molecule. 

---