Motion Maxwell-Boltzmann distribution

The basic principles are described in many textbooks [24, 26]. They are thus only sketchily presented here. In a conventional classical molecular dynamics calculation, a system of particles is placed within a cell of fixed volume, most frequently cubic in size. A set of velocities is also assigned, usually drawn from a Maxwell-Boltzmann distribution appropriate to the temperature of interest and selected in a way so as to make the net linear momentum zero. The subsequent trajectories of the particles are then calculated using the Newton equations of motion. Employing the finite difference method, this set of differential equations is transformed into a set of algebraic equations, which are solved by computer. The particles are assumed to interact through some prescribed force law. The dispersion, dipole-dipole, and polarization forces are typically included whenever possible, they are taken from the literature. 

The frequency with which the transition state is transformed into products, iT, can be thought of as a typical unimolecular rate constant no barrier is associated with this step. Various points of view have been used to calculate this frequency, and all rely on the assumption that the internal motions of the transition state are governed by thermally equilibrated motions. Thus, the motion along the reaction coordinate is treated as thermal translational motion between the product fragments (or as a vibrational motion along an unstable potential). Statistical theories (such as those used to derive the Maxwell-Boltzmann distribution of velocities) lead to the expression ... 

Transition State Theory [1,4] is the most frequently used theory to calculate rate constants for reactions in the gas phase. The two most basic assumptions of this theory are the separation of the electronic and nuclear motions (stemming from the Bom-Oppenheimer approximation [5]), and that the reactant internal states are in thermal equilibrium with each other (that is, the reactant molecules are distributed among their states in accordance with the Maxwell-Boltzmann distribution). In addition, the fundamental hypothesis [6] of the Transition State Theory is that the net rate of forward reaction at equilibrium is given by the flux of trajectories across a suitable phase space surface (rather a hypersurface) in the product direction. This surface divides reactants from products and it is called the dividing surface. Wigner [6] showed long time ago that for reactants in thermal equilibrium, the Transition State expression gives the exact... 

Begin by considering the translational motion of molecules in a container that has a total concentration of molecules per unit volume [c]. The distribution of velocities in the x, y, and z directions is given by the one-dimensional Maxwell-Boltzmann distribution of... 

In the kinetic theory of gases, the molecules are assumed to be smooth, rigid, and elastic spheres. The only kinetic energy considered is that from the translational motion of the molecules. In addition, the gas is assumed to be in an equilibrium state in a container where the gas molecules are uniformly distributed and all directions of the molecular motion are equally probable. Furthermore, velocities of the molecules are assumed to obey the Maxwell-Boltzmann distribution, which is described in the following section. 

We now proceed to develop a specific expression for the rate constant for reactants where the velocity distributions /a( )(va) and /B(J)(vB) for the translational motion are independent of the internal quantum state (i and j) and correspond to thermal equilibrium.4 Then, according to the kinetic theory of gases or statistical mechanics, see Appendix A.2.1, Eq. (A.65), the velocity distributions associated with the center-of-mass motion of molecules are the Maxwell-Boltzmann distribution, a special case of the general Boltzmann distribution law ... 

That is, the Maxwell-Boltzmann distribution for the two molecules can be written as a product of two terms, where the terms are related to the relative motion and the center-of-mass motion, respectively. After substitution into Eq. (2.18) we can perform the integration over the center-of-mass velocity Vx. This gives the factor y/2iVksTjM (IZo eXP( —ax2)dx = sjnja) and, from the equation above, we obtain the probability distribution for the relative velocity, irrespective of the center-of-mass motion. 

The average velocity for the motion from the left to the right over the barrier is then evaluated. From the one-dimensional Maxwell-Boltzmann distribution of velocities, Eq. (2.26),... 

The well-known Maxwell-Boltzmann distribution for the velocity or momentum associated with the translational motion of a molecule is valid not only for free molecules but also for interacting molecules in a liquid phase (see Appendix A.2.1). The average kinetic energy of a molecule at temperature T is, accordingly, (3/2)ksT. For the molecules to react in a bimolecular reaction they should be brought into contact with each other. This happens by diffusion when the reactants are dispersed in a solution, which is a quite different process from the one in the gas phase. For fast reactions, the diffusion rate of reactant molecules may even be the limiting factor in the rate of reaction. 

This deeper formulation of the question was shelved in Boltzmann s investigations by the soon emerging formulation of the 17-theorem (1872). Boltzmann took up this problem again only to counter Loschmidt s Umkehreinwand (Section 7b) and to obtain a modified formulation of the 17-theorem. He tried to show that, if we consider a motion of unlimited duration, then the Maxwell-Boltzmann distribution very strongly dominates in time over all other distributions, and hence the tendency to approach this particular distribution is quite understandable. 

Ill) If we consider a motion of the gas model which is of unlimited duration, the Maxwell-Boltzmann distribution will predominate overwhelmingly in time over all other appreciably different state distributions. 

All chemical process involves the motion of atoms within a molecule. Molecular dynamics (MD), in the broadest sense, is concerned with molecules in motion. It combines the energy calculations from molecular mechanics with equations of motion. Generally, an appropriate starting structure is selected (normally an energy minimized structure). Each atom in the system is then assigned a random velocity that is consistent with the Maxwell-Boltzmann distribution for the temperature of interest. The MM formalism is used to calculate the forces on all the atoms. Once the atom positions are known, the forces, velocities at time t, and the position of the atoms at some new time t + 5t can be predicted. More details about the method can be found in Ref.. ... 

An easy way to find this correction factor is to look at the history of an exit trajectory. This history is followed by starting atx = xb trajectories with velocity sampled from a Maxwell-Boltzmann distribution in the outward direction—these represent the outgoing equilibrium flux, then inverting the velocity (y -> —v) so that the particle is heading into the well, and integrating the equations of motion... 

Thus, although the individual molecular motions are chaotically unpredictable, their average behavior is entirely predictable and satisfies a particular probability distribution (the Maxwell-Boltzmann distribution). Quantities, such as the temperature that appear in Boyle s Law are measures of the average speed of the molecules. If we reran the tape of the history of our gas, we would find essentially the same average behavior, in accord with Boyle s Law, even though the individual trajectories of the molecules would be quite different. 

To obtain hyperpolarizabilities of calibrational quality, a number of standards must be met. The wavefunctions used must be of the highest quality and include electronic correlation. The frequency dependence of the property must be taken into account from the start and not be simply treated as an ad hoc add-on quantity. Zero-point vibrational averaging coupled with consideration of the Maxwell-Boltzmann distribution of populations amongst the rotational states must also be included. The effects of the electric fields (static and dynamic) on nuclear motion must likewise be brought into play (the results given in this section include these effects, but exactly how will be left until Section 3.2.). All this is obviously a tall order and can (and has) only been achieved for the simplest of species He, H2, and D2. Comparison with dilute gas-phase dc-SHG experiments on H2 and D2 (with the helium theoretical values as the standard) shows the challenge to have been met. 

As the next step we further simplify the problem by treating the translational motion classically, i.e., we will replace Po/M by the molecular center of mass velocity Fq and we will assume that between collisions Fq is a constant of the motion. In standard microwave spectroscopy, a Maxwell-Boltzmann distribution may be assumed for the molecular velocities Fq (compare the detailed discussion in Section III). Under these assumptions, Eq. (IV.55f) may be rewritten as... 

When H has reached its minimum value this is the well known Maxwell-Boltzmann distribution for a gas in thermal equilibrium with a uniform motion u. So, argues Boltzmann, solutions of his equation for an isolated system approach an equilibrium state, just as real gases seem to do. Up to a negative factor -k, in fact), differences in H are the same as differences in the thermodynamic entropy between initial and final equilibrium states. Boltzmann thought that his //-theorem gave a foundation of the increase in entropy as a result of the collision integral, whose derivation was based on the Stosszahlansatz. 

No special equilibrium between activated complexes and reactants is assumed. It is supposed however, that within the space between qi and qi+Aqi, q2 and q2+Aq2, configurations (activated complexes) have impulses (motions) between pi and pi+Api,p2 and P2+AP2 respectively. These configurations are computed in accordance with Maxwell-Boltzmann distribution. Recall from the gas laws that an energy profile for molecules can be describe by the Maxwell-Boltzman distribution diagram (Figure 3.3). As the temperature goes up, the population of molecules with more energy also increases. 

Beyond its influence on IMS separation parameters for a fixed geometry, the gas temperature affects ion geometries. As ion-molecule collisions must be sufficientiy frequent for a steady drift, ions are thermalized their internal, rotational, and translational modes are equilibrated at a single temperature. At low E/N where the ion drift is much slower than the Brownian motion of gas molecules, relative ion-molecule velocities conform to the Maxwell-Boltzmann distribution and ion temperature equals T of the gas. 

In a MD simulation, the system is placed within a cell of fixed volume, usually of cubic shape. A set of velocities is assigned, usually drawn from a Maxwell-Boltzmann distribution suitable for the temperature of interest and selected to make the linear momentum equal to zero. Then the trajectories of the particles are calculated by integration of the classical equation of motion. It is also assumed that the particles interact through some forces. 

Once the molecular assembly is constructed, energy minimization is necessary to reduce steric contacts and the system is heated to a desired temperature (corresponding to velocities in the Maxwell-Boltzmann distribution) to generate initial positions and velocities of the atoms. During dynamic runs, changes in position and velocity in response to the force field according to Newton s laws of motion are determined with the calculations repeated every 1 or 2fs. 

In its simplest form MD, considers a box of N particles and monitors their relative positions, velocities and accelerations by solving Newton s laws of motion at regular finite time intervals. Initially the particles are assigned pseudo-random velocities. These are often determined from a Maxwell-Boltzmann distribution and are required to meet certain conditions. These are that the kinetic energy of the system is such that the simulation temperature is fixed and that there is no net translational momentum. The forces acting on each particle, together with their velocities and positions are calculated for all subsequent time steps by considering Newton s Laws of Motion. If the time step is infinitely small then the acceleration, a, of an atom can be calculated from the force. 

---