Properties of the Density Matrix

The term density matrix arises by analogy to classical statistical mechanics, where the state of a system consisting of N molecules moving in a real three-dimensional space is described by the density of points in a 6N-dimensional phase space, which includes three orthogonal spatial coordinates and three conjugate momenta for each of the N particles, thus giving a complete description of the system at a particular time. In principle, the density matrix for a spin system includes all the spins, as we have seen, and all the spatial coordinates as well. However, as we discuss subsequently we limit our treatment to spins. For simplicity we deal only with application to systems of spin % nuclei, but the formalism also applies to nuclei of higher spin. 

From Eq. 11.5 it is apparent that ( , = p m, which defines a Hermitian matrix, so the density matrix is seen to be Hermitian. For a system of N spins /2, p is 2N X 2n in size. From the fact that the probabilities p f in Eq. 11.5 must sum to unity and the basis functions used are orthonormal, it can be shown that the trace of the density matrix is 1. 

We can simplify the mathematical expressions by recognizing that the sum in Eq. 11.5 is just an ensemble average, so this equation can be rewritten as 

As we observed previously, the phases are random hence there is no correlation between the values of a, and ciy and we may average separately to give 

0 only if both the average value of c c and the average value of exp[i(an — , )] are nonzero. The latter requirement means that there must be a coherence between the phases am and an for the corresponding element of p to be nonzero. 

---

Husimi, K., Proc. Phys.-Math. Soc. Japan 22, 264, "Some formal properties of the density matrix." Introduction of the concept of reduced density matrix. Statistical-mechanical treatment of the Hartree-Fock approximation at an arbitrary temperature and an alternative method of obtaining the reduced density matrices are discussed. 

K. Husimi, Some formal properties of the density matrix. Proc. Phys. Math. Soc. Japan 22, 264 (1940). 

A few remarks are necessary, concerning the general properties of the density matrix. Keeping in mind Eq. (17) and that the density matrix is Hermitian (p = P a) we obtain the following symmetry relations ... 

In this case the transformation properties of the density matrix... 

The convergence properties of the density matrix-based equations, i.e., the number of iterations to converge P , are similar to the ones encountered for a solution in the MO space, so that the advantage of using sparse multiplications within the density-based approach allows us to reduce the scaling property of the computational effort in an efficient manner. In this way, NMR chemical shift calculations with linear-scaling effort become possible and systems with 1000 and more atoms can be treated at the HF or DFT level on today s computers.Extensions to other molecular properties can be formulated in a similar fashion. 

It is also of interest to study the "inverse" problem. If something is known about the symmetry properties of the density or the (first order) density matrix, what can be said about the symmetry properties of the corresponding wave functions In a one electron problem the effective Hamiltonian is constructed either from the density [in density functional theories] or from the full first order density matrix [in Hartree-Fock type theories]. If the density or density matrix is invariant under all the operations of a space CToup, the effective one electron Hamiltonian commutes with all those elements. Consequently the eigenfunctions of the Hamiltonian transform under these operations according to the irreducible representations of the space group. We have a scheme which is selfconsistent with respect to symmetty. 

This equation relates the temporal concentration of a diffusing chemical to its location in space. In real soil and aquifer materials, the diffusion coefficient can be affected by the temperature and properties of the solid matrix, such as mineral composition (which affects sorption, a process that can be difficult to separate from diffusion), bulk density, and critically, water content. 

The problem is that this exchange correlation potential does not in a simple fashion depend on only the density or the diagonal part of the one-matrix. The antisymmetry property of the two-matrix, which yields the last term in Eq. (33) as well as the apperance of the coefficents rkj kj, prevent this. Only some further (ad hoc) assumptions will allow an exchange-correlation potential that is simply a function of the electron density. This can be achieved in a number of ways. For a full discussion on this subject see [11]. 

For multichannel scattering where there are two or more open channels, the S matrix is a true matrix with elements Sy and the cross section for the transition from channel i to channel j is proportional to 5y - Sy 2. The symmetry of collision processes with respect to the time reversal leads to the symmetric property of the S matrix, ST = S, which, in turn, leads to the principle of detailed balance between mutually reverse processes. The conservation of the flux of probability density for a real potential and a real energy requires that SSf = SfS = I, i.e., S is unitary. For a complex energy or for a complex potential, in general, the flux is not conserved and S is non-unitary. 

We have written the operator Fl(x) as a function of the combined space-spin coordinates X, because while the spin summations can be carried out in Jl(x) before calculating matrix elements, Kl(x) may connect spin-orbitals that are off-diagonal in the spin wavefunctions however in the special case of the density matrix p (xi, Xa) arising from a wavefunction that is a spin singlet (5 = 0) one can show that must also be diagonal. This leads to a useful simplification here since we can usually assume this property for Wlo, and it means that Vl(x) reduces to a (non-local) function of the space variable r only we can therefore consistently parameterize the matrix elements for the whole potential, (/bI Vl(x) j) without having to decompose them into different spin combinations for the Coulomb and exchange potentials. 

The presence of filler in the rubber as well as the increase of the surface ability of the Aerosil surface causes an increase in the modulus. The temperature dependence of the modulus is often used to analyze the network density in cured elastomers. According to the simple statistical theory of rubber elasticity, the modulus should increase twice for the double increase of the absolute temperature [35]. This behavior is observed for a cured xmfilled sample as shown in Fig. 15. However, for rubber filled with hydrophilic and hydrophobic Aerosil, the modulus increases by a factor of 1.3 and 1.6, respectively, as a function of temperature in the range of 225-450 K. It appears that less mobile chain units in the adsorption layer do not contribute directly to the rubber modulus, since the fraction of this layer is only a few percent [7, 8, 12, 21]. Since the influence of the secondary structure of fillers and filler-filler interaction is of importance only at moderate strain [43, 47], it is assumed that the change of the modulus with temperature is mainly caused by the properties of the elastomer matrix and the adsorption layer which cause the filler particles to share deformation. Therefore, the moderate decrease of the rubber modulus with increasing temperature, as compared to the value expected from the statistical theory, can be explained by the following reasons a decrease of the density of adsorption junctions as well as their strength, and a decrease of the ability of filler particles to share deformation due to a decrease of elastomer-filler interactions. 

Charge density and related nominal capacity Properties of the support matrix, including composition and pore diameter... 

Because T operates on each element of a matrix it is called a superoperator. In fact, the Hilbert-space formulation of quantum mechanics leading to the von Neumann equation of motion of the density matrix can be simplified considerably by introduction of a superoperator notation in the so-called Liouville space. Furthermore, for the analysis of NMR experiments with complicated pulse sequences it is of great help to expand the density matrix into products of operators, where each product operator exhibits characteristic transformation properties under rotation [Eml]. 

In order to introduce some notation, we first recall a few of the well-known properties of the interaction of light pulses with molecules in the linear approximation. Frequently, the signals in nonlinear optical experiments are expressed in terms of the polarization induced in the medium by the incident pulses. The complex linear polarization P t) vector for a distribution of identical two-level systems is obtained from an elementary calculation of the density matrix using the Liouville equation of a system perturbed by an electric field and proceeding as follows ... 

The relations between expectation values of properties and the density matrix will be the same as in statistical thermodynamics and, therefore, will conform to the laws of thermodynamics. 

The present picture is based on speciahzation of the density matrix amplitude (x, Tj x, such as to become uniformly in the valence shell properties, such as the space-time Bohr-Slater quantification on orbits (in atomic units m = h = e l4ns = 1), see (Bohr, 1921 White, 1934)... 

---