Molecular modelling density matrices

The concept of the molecular orbital is, however, not restricted to the Hartree-Fock model. Sets of orbitals can also be constructed for more complex wave functions, which include correlation effects. They can be used to obtain insight into the detailed features of the electron structure. One choice of orbitals are the natural orbitals, which are obtained by diagonalizing the spinless first-order reduced density matrix. The occupation numbers (T ) of the natural orbitals are not restricted to 2, 1, or 0. Instead they fulfill the condition ... 

The next two chapters are devoted to ultrafast radiationless transitions. In Chapter 5, the generalized linear response theory is used to treat the non-equilibrium dynamics of molecular systems. This method, based on the density matrix method, can also be used to calculate the transient spectroscopic signals that are often monitored experimentally. As an application of the method, the authors present the study of the interfadal photo-induced electron transfer in dye-sensitized solar cell as observed by transient absorption spectroscopy. Chapter 6 uses the density matrix method to discuss important processes that occur in the bacterial photosynthetic reaction center, which has congested electronic structure within 200-1500cm 1 and weak interactions between these electronic states. Therefore, this biological system is an ideal system to examine theoretical models (memory effect, coherence effect, vibrational relaxation, etc.) and techniques (generalized linear response theory, Forster-Dexter theory, Marcus theory, internal conversion theory, etc.) for treating ultrafast radiationless transition phenomena. 

In the present article, we report a study concerning the reaction mechanism of a prototype reaction using both static and dynamic approaches to explore a DFT potential surface. The static approach is the standard IRC model, while the dynamic one is based on a Carr-Parrinello method performed with localized (Gaussian) orbitals, the so-called atom-centered density matrix propagation (ADMP) model.25 Our aim is to elucidate the differences, and the common aspects, between the two approaches in the analysis of bond breaking/formation. To this end, we have chosen topological quantities as probe molecular descriptors. 

According to the assumption we have made the change in the density matrix, ARX, due to the coulombic interaction between fragments will be more or less localized. It is tempting to set ARX = XL. By doing that, however, one is forced [8] to split off the local space from the remainder of the system to satisfy the idempotency condition. This results in an ordinary cluster model which does not allow electron transfer to or from the surroundings and, as we will see in Sect. 5, is unsuitable for our purposes. In order to properly embed the cluster we take advantage of the fact that the sum of the occupied and unoccupied molecular orbital (MO) spaces is identical to the total AO space. So, instead of ARX = XL, we write... 

It is neither feasible nor illustrative to detail a general solution of the self-consistency equations. The five input parameters in the crystal field model, Dq, Cso, F(h F2, and F4, determine the five independent density matrix elements when the choice of chemical potential /( and the relative population of the molecular... 

The concept of the molecular orbital and their occupation is, however, not restricted to the HF model. It has much wider relevance and is applicable also for more accurate wave functions. For each wave function we can form the first-order reduced density matrix. This matrix is Hermitian and can be diagonalized. The basis for this diagonal form of the density matrix are the Natural Orbitals first introduced in quantum chemistry by Per-Olof Lowdin [4]. 

The general theory outlined in Chapter 14 has been used in the ROHF model. It can also provide the basis for a whole class of multi-determinant models of molecular electronic structure which are constructed from sets o/doubly occupied orbitals. Two such models are outlined in this chapter Paired-Electron MCSCF (PEMCSCF) and the General Valence Bond (GVB). These pair expansion theories are based on the idea of natural orbitals which diagonalise the one-electron density matrix. 

We have seen earlier that the Hartree-Fock model of molecular electronic structure is basically a density-matrix theory. The physical interpretation and even the formalism of the theory may be expressed in terms of the basic invariant of the theory either expressed as a function of six spatial and two spin variables or as an orbital expansion in terms of some set of MOs Xi ... 

For convenience, we shall classify the molecular models according to their topological dimensionality, p. A molecular conformation defined by the set of nuclear position vectors is a zero-dimensional (OD) model. A one-dimensional (ID) model corresponds to a molecular skeleton, defined by the set of nuclear positions and their connectivity (bond) matrix. Contour surfaces of one-particle molecular properties such as electron density or electrostatic potential are topologically two-dimensional (2D) models embedded in three dimensions. Finally, we find a true three-dimensional (3D) model whenever an entire one-electron property over all space is involved. This model can be regarded as the continuum of all 2D isoproperty surfaces. The difference among the models is summarized in Figure 1. We shall deal with pD models in this work (p = 0,1, 2, 3). Each of them requires a different type of shape descriptor. 

A theory for the ultrafast pump-probe spectroscopy of large polyatomic molecules in condensed phases was developed in the work [15]. A multimode Brownian oscillator model was used to account for high-frequency molecular vibrations and local intermolecular modes as well as collective solvent motions. A semiclassical picture was provided using the density matrix in Liouville space. Conditions for the observation of quantum beats, spectral diffusion, and solvation dynamics (dynamic Stokes shift) are specified. 

---