Molecular heat baths

In both cases, because of restrictions imposed on the excitation process (e.g. optical selection rules), the initially excited state is not an exact eigenstate of the molecular Hamiltonian (see below). At the same time, if the molecule is large enough, this initially prepared zero-order excited state is embedded in a bath of a very large number of other states. Interaction between these zero-order states results from residual molecular interactions such as corrections to the Bom Oppenheimer approximation in the first example and anharmonic corrections to nuclear potential surfaces in the second. These exist even in the absence of interactions with other molecules, giving rise to relaxation even in isolated (large) molecules. The quasi-continuous manifolds of states are sometimes referred to as molecular heat baths. The fact that these states are initially not populated implies that these baths are at zero temperature. 

In a molecular dynamics calculation, you can add a term to adjust the velocities, keeping the molecular system near a desired temperature. During a constant temperature simulation, velocities are scaled at each time step. This couples the system to a simulated heat bath at Tq, with a temperature relaxation time of "r. The velocities arc scaled bv a factor X. where... 

For a constant temperature simulation, a molecular system is coupled to a heat bath via a Bath relaxation constant (see Temperature Control on page 72). When setting this constant, remember that a small number results in tight coupling and holds the temperature closer to the chosen temperature. A larger number corresponds to weaker coupling, allowing more fluctuation in temper-... 

The simplest method that keeps the temperature of a system constant during an MD simulation is to rescale the velocities at each time step by a factor of (To/T) -, where T is the current instantaneous temperature [defined in Eq. (24)] and Tq is the desired temperamre. This method is commonly used in the equilibration phase of many MD simulations and has also been suggested as a means of performing constant temperature molecular dynamics [22]. A further refinement of the velocity-rescaling approach was proposed by Berendsen et al. [24], who used velocity rescaling to couple the system to a heat bath at a temperature Tq. Since heat coupling has a characteristic relaxation time, each velocity V is scaled by a factor X, defined as... 

Concluding this section, one should mention also the method of molecular dynamics (MD) in which one employs again a bead-spring model [33,70,71] of a polymer chain where each monomer is coupled to a heat bath. Monomers which are connected along the backbone of a chain interact via Eq. (8) whereas non-bonded monomers are assumed usually to exert Lennard-Jones forces on each other. Then the time evolution of the system is obtained by integrating numerically the equation of motion for each monomer i... 

G. Grest, K. Kremer. Molecular dynamics simulation for polymers in the presence of heat bath. Phys Rev A 55 3628-3631, 1986. 

A higher level of understanding would require a knowledge of molecular dynamics and presently represents a rather distant goal. In addition to reliable knowledge of the shapes of potential energy hypersurfaces, it would also require information such as vibronic coupling elements, densities of vibrational states, detailed mechanism of the action of the heat bath, etc. 

The book covers a variety of questions related to orientational mobility of polar and nonpolar molecules in condensed phases, including orientational states and phase transitions in low-dimensional lattice systems and the theory of molecular vibrations interacting both with each other and with a solid-state heat bath. Special attention is given to simple models which permit analytical solutions and provide a qualitative insight into physical phenomena. 

We now present results from molecular dynamics simulations in which all the chain monomers are coupled to a heat bath. The chains interact via the repiflsive portion of a shifted Lennard-Jones potential with a Lennard-Jones diameter a, which corresponds to a good solvent situation. For the bond potential between adjacent polymer segments we take a FENE (nonhnear bond) potential which gives an average nearest-neighbor monomer-monomer separation of typically a 0.97cr. In the simulation box with a volume LxL kLz there are 50 (if not stated otherwise) chains each of which consists of N -i-1... 

Once the boundary conditions have been implemented, the calculation of solution molecular dynamics proceeds in essentially the same manner as do vacuum calculations. While the total energy and volume in a microcanonical ensemble calculation remain constant, the temperature and pressure need not remain fixed. A variant of the periodic boundary condition calculation method keeps the system pressure constant by adjusting the box length of the primary box at each step by the amount necessary to keep the pressure calculated from the system second virial at a fixed value (46). Such a procedure may be necessary in simulations of processes which involve large volume changes or fluctuations. Techniques are also available, by coupling the system to a Brownian heat bath, for performing simulations directly in the canonical, or constant T,N, and V, ensemble (2,46). 

Open quantum systems have attracted much attention over the last decades. While most of the studies dealt with systems coupled to bosonic heat baths, recently systems coupled to fermionic reservoirs describing for example molecular wires have been in the focus of many investigations. This chapter will not try to give a concise overview of the available literature but will focus on a particular approach time-local (TL) quantum master equations (QMEs) and in particular their combination with specific forms of the spectral density. 

Thus, the key point is [13-15] to circumvent the Boltzmann distribution (i.e., the temperature ) by replacing the heat bath by a controllable electromagnetic field. This generates distributions over the molecular states that are inaccessible when energy is supplied in the form of heat. In addition, when the electromagnetic field is generated by a laser there is a more fundamental difference laser excitation is a coherent excitation of states this is explained in more detail below. 

Although a theoretical approach has been desecrated as to how one can apply the generalized coupled master equations to deal with ultrafast radiationless transitions taking place in molecular systems, there are several problems and limitations to the approach. For example, the number of the vibrational modes is limited to less than six for numerical calculations. This is simply just because of the limitation of the computational resources. If the efficient parallelization can be realized to the generalized coupled master equations, the limitation of the number of the modes can be relaxed. In the present approach, the Markov approximation to the interaction between the molecule and the heat bath mode has been employed. If the time scale of the ultrashort measurements becomes close to the characteristic time of the correlation time of the heat bath mode, the Markov approximation cannot be applicable. In this case, the so-called non-Markov treatment should be used. This, in turn, leads to a more computationally demanding task. Thus, it is desirable to develop a new theoretical approach that allows a more efficient algorithm for the computation of the non-Markov kernels. Another problem is related to the modeling of the interaction between the molecule and the heat bath mode. In our model, the heat bath mode is treated as... 

Apart from the heat bath mode, the harmonic potential surface model has been used for the molecular vibrations. It is possible to include the generalized harmonic potential surfaces, i.e., displaced-distorted-rotated surfaces. In this case, the mode coupling can be treated within this model. Beyond the generalized harmonic potential surface model, there is no systematic approach in constructing the generalized (multi-mode coupled) master equation that can be numerically solved. The first step to attack this problem would start with anharmonicity corrections to the harmonic potential surface model. Since anharmonicity has been recognized as an important mechanism in the vibrational dynamics in the electronically excited states, urgent realization of this work is needed. 

We have derived a formula for the molecular partition function by considering a system containing many molecules at equilibrium with a heat bath. We can generalize our statistical mechanics by a gedanken experiment of considering a large number of identical systems, each with volume V and number of particles N at equilibrium with the heat bath at temperature T. Such a supersystem is called a canonical ensemble. Our derivation is the same the fraction of systems that are in a state with energy Et is... 

---