Molecules , ensembles

Figure 8.1 Simplified representation of slice of free-energy hypersurface for C10H22 is shown. Free energy vs. configuration (conformation) for two molecules, ensembles of which represent two isomeric compounds n-decane (A) and isodecane (B) are given.

In the constitutional model of Ugi, rather than molecules, "ensemble of molecules (EM) are used in which the molecules can be either chemically different or identical. Like molecules, an EM has an empirical formula, which is the sum of the empirical formulae of the constiment molecules and describes the collection A of atoms within the EM under consideration. All the EM s which can be formed from A have the same empirical formula . Therefore, an EM(A) consists of one or more molecules which can be obtained from A using each atom which belongs to A only once. Moreover, a FIEM(A) or a family of isomeric EM, is the collection of all EM(A) and it is determined by the empirical formula . On the other hand, a chemical reaction, or a sequence of chemical reactions, is the conversion of an EM into an isomeric EM, and therefore a FIEM contains all EMs which are chemically interconvertible, as far as stoichiometry is concerned. In summary, a FIEM(A) contains, at least in principle, the whole chemistry of the collection A of atoms and since any collection of atoms may be chosen here, Ugi concludes that a theory of FIEM is, in fact, a theory of all chemistry. 

PML can also represent polymers at different levels of certainty where the structure is known, it can be encoded and where this information is not available, a polymer can be codified in terms of other concepts, such as the monomers it was prepared from. Furthermore, PML provides essentially a coarse-grained representation larger structural fragments can be mapped back to fully atomistic fragments, if desired. PML, therefore, allows data to be associated with a polymer representation at the atom, molecular fragment, molecule and molecule ensemble level, and. 

For fuzzy virtual screening purposes, highly active molecules with different scaffolds are combined into an MTree model. By combining the information of remotely related actives into a single model, efficient database searches with molecule ensembles are possible. 

The continuity and momentum equations can be solved by tracking the movements of molecule ensembles through the evolution of the distribution function [20] using the popular LBM. The lattice Boltzmann equation (LBE) can be derived from the Boltzmann equation. For the flows with external forces, the continuous Boltzmann-BGK equation with an external force term, F, is... 

The microcanonical ensemble is a set of systems each having the same number of molecules N, the same volume V and the same energy U. In such an ensemble of isolated systems, any allowed quantum state is equally probable. In classical thennodynamics at equilibrium at constant n (or equivalently, N), V, and U, it is the entropy S that is a maximum. For the microcanonical ensemble, the entropy is directly related to the number of allowed quantum states C1(N,V,U) ... 

In this chapter, the foundations of equilibrium statistical mechanics are introduced and applied to ideal and weakly interacting systems. The coimection between statistical mechanics and thennodynamics is made by introducing ensemble methods. The role of mechanics, both quantum and classical, is described. In particular, the concept and use of the density of states is utilized. Applications are made to ideal quantum and classical gases, ideal gas of diatomic molecules, photons and the black body radiation, phonons in a hannonic solid, conduction electrons in metals and the Bose—Einstein condensation. Introductory aspects of the density... 

No system is exactly unifomi even a crystal lattice will have fluctuations in density, and even the Ising model must pemiit fluctuations in the configuration of spins around a given spin. Moreover, even the classical treatment allows for fluctuations the statistical mechanics of the grand canonical ensemble yields an exact relation between the isothemial compressibility K j,and the number of molecules Ain volume V ... 

RRKM theory assumes a microcanonical ensemble of A vibrational/rotational states within the energy interval E E + dE, so that each of these states is populated statistically with an equal probability [4]. This assumption of a microcanonical distribution means that the unimolecular rate constant for A only depends on energy, and not on the maimer in which A is energized. If N(0) is the number of A molecules excited at / =... 

In the above discussion it was assumed that the barriers are low for transitions between the different confonnations of the fluxional molecule, as depicted in figure A3.12.5 and therefore the transitions occur on a timescale much shorter than the RRKM lifetime. This is the rapid IVR assumption of RRKM theory discussed in section A3.12.2. Accordingly, an initial microcanonical ensemble over all the confonnations decays exponentially. However, for some fluxional molecules, transitions between the different confonnations may be slower than the RRKM rate, giving rise to bottlenecks in the unimolecular dissociation [4, ]. The ensuing lifetime distribution, equation (A3.12.7), will be non-exponential, as is the case for intrinsic non-RRKM dynamics, for an mitial microcanonical ensemble of molecular states. 

The first classical trajectory study of iinimoleciilar decomposition and intramolecular motion for realistic anhannonic molecular Hamiltonians was perfonned by Bunker [12,13], Both intrinsic RRKM and non-RRKM dynamics was observed in these studies. Since this pioneering work, there have been numerous additional studies [9,k7,30,M,M, ai d from which two distinct types of intramolecular motion, chaotic and quasiperiodic [14], have been identified. Both are depicted in figure A3,12,7. Chaotic vibrational motion is not regular as predicted by tire nonnal-mode model and, instead, there is energy transfer between the modes. If all the modes of the molecule participate in the chaotic motion and energy flow is sufficiently rapid, an initial microcanonical ensemble is maintained as the molecule dissociates and RRKM behaviour is observed [9], For non-random excitation initial apparent non-RRKM behaviour is observed, but at longer times a microcanonical ensemble of states is fonned and the probability of decomposition becomes that of RRKM theory. 

Apparent non-RRKM behaviour occurs when the molecule is excited non-randomly and there is an initial non-RRKM decomposition before IVR fomis a microcanonical ensemble (see section A3.12.2). Reaction patliways, which have non-competitive RRKM rates, may be promoted in this way. Classical trajectory simulations were used in early studies of apparent non-RRKM dynamics [113.114]. 

As discussed in section A3.12.2. intrinsic non-RRKM behaviour occurs when there is at least one bottleneck for transitions between the reactant molecule s vibrational states, so drat IVR is slow and a microcanonical ensemble over the reactant s phase space is not maintained during the unimolecular reaction. The above discussion of mode-specific decomposition illustrates that there are unimolecular reactions which are intrinsically non-RRKM. Many van der Waals molecules behave in this maimer [4,82]. For example, in an initial microcanonical ensemble for the ( 211 )2 van der Waals molecule both the C2H4—C2H4 intennolecular modes and C2H4 intramolecular modes are excited with equal probabilities. However, this microcanonical ensemble is not maintained as the dimer dissociates. States with energy in the intermolecular modes react more rapidly than do those with the C2H4 intramolecular modes excited [85]. 

A diagrannnatic approach that can unify the theory underlymg these many spectroscopies is presented. The most complete theoretical treatment is achieved by applying statistical quantum mechanics in the fonn of the time evolution of the light/matter density operator. (It is recoimnended that anyone interested in advanced study of this topic should familiarize themselves with density operator fonnalism [8, 9, 10, H and f2]. Most books on nonlinear optics [13,14, f5,16 and 17] and nonlinear optical spectroscopy [18,19] treat this in much detail.) Once the density operator is known at any time and position within a material, its matrix in the eigenstate basis set of the constituents (usually molecules) can be detennined. The ensemble averaged electrical polarization, P, is then obtained—tlie centrepiece of all spectroscopies based on the electric component of the EM field. 

Consider an ensemble composed of constituents (such as molecules) per unit volume. The (complex) density operator for this system is developed perturbatively in orders of the applied field, and at. sth order is given by The (complex). sth order contribution to the ensemble averaged polarization is given by the trace over the eigenstate basis of the constituents of the product of the dipole operator, N and = Tr A pp... 

The nonlinear response of an individual molecule depends on die orientation of the molecule with respect to the polarization of the applied and detected electric fields. The same situation prevails for an ensemble of molecules at an interface. It follows that we may gamer infonnation about molecular orientation at surfaces and interfaces by appropriate measurements of the polarization dependence of the nonlinear response, taken together with a model for the nonlinear response of the relevant molecule in a standard orientation. 

Here the ijk coordinate system represents the laboratory reference frame the primed coordinate system i j k corresponds to coordinates in the molecular system. The quantities Tj, are the matrices describing the coordinate transfomiation between the molecular and laboratory systems. In this relationship, we have neglected local-field effects and expressed the in a fomi equivalent to simnning the molecular response over all the molecules in a unit surface area (with surface density N. (For simplicity, we have omitted any contribution to not attributable to the dipolar response of the molecules. In many cases, however, it is important to measure and account for the background nonlinear response not arising from the dipolar contributions from the molecules of interest.) In equation B 1.5.44, we allow for a distribution of molecular orientations and have denoted by () the corresponding ensemble average ... 

An important point for all these studies is the possible variability of the single molecule or single particle studies. It is not possible, a priori, to exclude bad particles from the averaging procedure. It is clear, however, that high structural resolution can only be obtained from a very homogeneous ensemble. Various classification and analysis schemes are used to extract such homogeneous data, even from sets of mixed states [69]. In general, a typical resolution of the order of 1-3 mn is obtained today. 

A consideration of the transition probabilities allows us to prove that microscopic reversibility holds, and that canonical ensemble averages are generated. This approach has greatly extended the range of simulations that can be perfonned. An early example was the preferential sampling of molecules near solutes [77], but more recently, as we shall see, polymer simulations have been greatly accelerated by tiiis method. 

For a multicomponent system, it is possible to simulate at constant pressure rather than constant volume, as separation into phases of different compositions is still allowed. The method allows one to study straightforwardly phase equilibria in confined systems such as pores [166]. Configuration-biased MC methods can be used in combination with the Gibbs ensemble. An impressive demonstration of this has been the detennination by Siepmaim et al [167] and Smit et al [168] of liquid-vapour coexistence curves for n-alkane chain molecules as long as 48 atoms. 

Yoon K, Chae D G, Ree T and Ree F H 1981 Computer simulation of a grand canonical ensemble of rodlike molecules J. Chem. Phys. 74 1412-23... 

The development of tunable, narrow-bandwidtli dye laser sources in tire early 1970s gave spectroscopists a new tool for selectively exciting small subsets of molecules witliin inhomogeneously broadened ensembles in tire solid state. The teclmique of fluorescence line-narrowing [1, 2 and 3] takes advantage of tire fact tliat relatively rigid chromophoric... 

---