Studies of the optical transitions in electron systems are powerful for examining the electronic structure of solids, because the absorption lineshapes are quite sensitive to the type of scattering mechanism affecting the transport of electrons and to the interaction of electrons with intense laser light. A general method for gaining knowledge on the dynamics of a system is a perturbation-based study. This consists of dividing the Hamiltonian into an exactly soluble part and a nontrivial perturbativepart, the effect of which is studied in order. The most popular method for representing the terms in perturbative expressions is Feynman diagram. This diagrammatic method can be used directly for reasoning and problem solving as well as for representing the perturbative expressions by drawings. The easily recognizable topology of the diagrams makes the diagrammatic method a powerful tool for constructing approximation schemes. Furthermore, the diagrammatic representation can be a suggestive tool providing physical intuition to quantum dynamics by increasing the diagrams to a representation method for possible alternative physical processes.

On the other hand, the method, although invented originally for particle physics, has been adopted in solid-state physics, where the behavior of phonons may be expressed in analogy to that of photons, e.g. in the theory of superconductivity. An electron traveling in a solid distorts the lattice due to the Coulomb interactions with the ions. The lattice distortion in turn has a feedback on the electron dynamics, resulting in an increase in the electron mass and a shortening of the electron lifetime in a particular quasi-particle state. This effect is described in terms of the self-energy that the electron acquires due to the electron-phonon interaction. The real part of selfenergy describes the change in electron energy, and the imaginary part provides information on the electron lifetime. The electron self-energy can be calculated using the standard Feynman diagram.

The standard diagram method can represent the trajectories of particles well in the intermediate states of the scattering processes. However, in the line shape (or selfenergy) function for electron-phonon system, the Fermi distribution functions for electrons and the Planck distribution functions for phonons are simply added [1-7], which violates the “population criterion” in that the population of electrons (fermions) and phonons (bosons) appear independently. In other words, a theory can be said to be proper if the Fermi distributions are multiplied by the Planck distributions in the formalism.

The present authors have developed some projection methods for the optical transitionsin electron-phonon systems and used them to calculate the linewidths in semiconductors [8,9]. Normally, the resolvent factor contained in the conductivity tensor is expanded using projectors, and various formulae can be obtained. Recently, the formalism was improved with the inclusion of nonlinear terms near the resonance points and suggested a meaningful result including the Fermi and Planck distribution functions properly with the proper use of the SDPO and KCRI [10-13]. This paper introduces a method for the nonlinear optical conductivity and line shape functions to represent them in loop diagrams. It can be categorized as another Feynman-like scheme because all the contributions to the line shape functions (or self-energy terms) can be grouped into topologicallydistinct diagrams.

The diagram approach to the nonlinear phenomena is based on the following methods.

For diagrammatic representation of the line shape functions, , we consider the form

where.

1) Let an implicit state exist between state andstate.

2) There is a diagram (Figure 1) connecting the initial state to the implicit state with a dotted line and connecting the implicit state to a state with aclockwise loop. The dotted line and the clockwise loop correspond to and, respectively.

Here called the C-factor and called the P-factor are defined as follows:

where is the coupling factor that depends on the mode of the phonons and is the Planck distribution function for phonons. means that the state is coupled with the state by a phonon with wave vector q and means an electron transition from a state (distribution function:) to a state (distribution function:) with a phonon emission (distribution function:) minus a transition from a state to a state with a phonon absorption. There is an another diagram connecting with a counterclockwise loop, , which is defined

In a loop, the upper and lower half circles correspond to phonon emission () and phonon absorption ().

3) There are diagrams connecting the final state to the implicit state with dotted line and connecting the implicit state to a state with a clockwise (or counterclockwise) loop. Figure 2 corresponds to clockwise loop.

4) Assign for the clockwise loop , and assign for the counterclockwise loop. Here called the G-factors are defined as follows:

By, the energy conservation is satisfied, i.e.,.

5) Let the states connected to the dotted line involve no loops.

6) a) Multiply, , and for the clockwise loop. b) Multiply, , and for the counterclockwise loop.

7) Finally, sum all the diagrams after summing eachdiagram over all q and values for the line shape function.

When an electromagnetic wave of frequency is applied to a system along the l (x, y, z) direction, the linear optical conductivity is given by the following: [10-12]

where for electron states and, j_k is the k component of the single electron current op-

erator, r_l is the l component of the electron position vector, is the Fermi distribution function for an electron with energy, and. Equation (6) can be represented by the diagram in Figure 3.

The energy term in Equation (6),

, represents the transition from the initial state to the final state with a photon absorption of frequency. If there was no phonon scattering, the line shape would be like a delta function. However, as the electrons are scattered by phonons, the shape is broadened, so the line shape function is involved.

In Equation (6), there is no intermediate state between the initial state and final state, so the line shape function is represented by four loop diagrams in Figure 4.

From Figure 4, is obtained as

which is identical to the line shape function given by Equation (4.8) in [11].

The physical meaning of the first term in Equation (7) is as follows (see the first term in Figure 5). Since the absorption power delivered to the system is proportional to the real part of the conductivity [10,11], becomes Dirac’s delta, i.e.,

, since can be replaced by () in the response scheme.

Note that where P means principal value. Therefore, means that the implicit state is determined by the energy of the final state, photon energy, and phonon energy, so that energy conservation can be satisfied, i.e.. means that the implicit state is coupled with the initial state by a phonon with wave vector q. means that the reverse implicit transition from the final state to the intermediate state with phonon absorption should be subtracted from the forward implicit transition from the intermediate state to the final state with a phonon emission. This means that when an electron-phonon interaction is involved in the electron transition, there are local fluctuations, i.e. a transition occurs via implicit states. In other words, an electron undergoes a transition from the implicit state, which is coupled with the initial state, to the final state (or vice versa) with phonon emission (or absorption). The transition forms a loop because the absorption process maintains a balance with the emission process. Although is called the line shape function, the line shape must be calculated from Equation (6) not directly from Equation (7). Therefore, all the states given by are called the implicit states because they are included only in. Therefore, the transition from the initial state to the final state occurs via two implicit transitions, and the implicit state is connected to the initial state by and to the final state by.

Although the implicit transitions are not measured directly, they should be considered in the calculations. Note that the other theories [1-7] cannot provide any diagram representation because they contain the sums of two distribution functions, such as

The other three terms in can be represented in a similar manner.

In Figure 4 or Equation (7), the first two terms are

topologically equivalent because the absorption of photons appears while phonons can be both absorbed and emitted in the transitions. The absorption (emission) of a photon or a phonon in the forward transitions is identical to the emission (absorption) in the reverse transitions, with negative contributions (minus signs) assigned to the reverse transitions. The final two terms are also topologically equivalent for the same reason. A net transition is possible because the implicit states are determined variously according to the energies of the final (or initial) state, photon energies, and phonon energies.

The first-order nonlinear optical conductivity for two incident electromagnetic waves of frequencies and is given by the following [10-12]

where. For the two photon process, there is an intermediate state between the initial state and final state. The first term in Equation (8) means thatthere is an intermediate transition from the intermediatestate to the final state with photon () absorption and line shape function and the direct transition from the initial state to the final state with two photons (,) absorption and line shape function (Figure 5). Note that

.

is called the intermediate state to be discerned from the implicit state because it is explicitly included in the conductivity but not in the line shape function.

Four loop diagrams exist because there is an implicit state between state and state through in the first term in Equation (8) (Figure 5). The first two terms involve in the transition of Figure 6, considering as a pseudo- state through and other two terms involve in the transition of Figure 5.

Therefore can be expressed as follows:

which is identical to the line shape function given by Equation (5.7) in [11]. The physical meaning of the third term in Equation (9) is as follows: means that the implicit state is determined by the energy of the final state, two photon energies, and phonon energy; and means that the implicit state is coupledwith the initial state by a phonon with a wave vector q. There are two forward implicit transitions and two reverse implicit transitions with a center on the implicit state. All the implicit states and types of phonons are included in the processes by the sums in Equation (9). The other three terms in can be interpreted in a similar manner. Similarly, can be obtained. In Figure 6 or Equation (9), the first two terms and last two terms are topologically equivalent pairs for the same reason as in Figure 5 or Equation (7).

In conclusion, this paper introduced a loop diagram approach to the nonlinear optical conductivity formula for an electron-phonon system. It was possible to explain the phonon emissions and absorptions in all electron transition processes in an organized manner because the line shape functions include the electron distribution functions properly as well as the phonon distribution functions. Since the present diagram method is not the one representing the trajectories of particles in the intermediate stages of scattering processes, these diagrams should not be confused with the time-ordered diagrams in the Feynman scheme [14] or with the temperature diagrams

in the Feynman-like scheme [15,16]. However, this method can be classified as another Feynman-like scheme, and shall be called the “KC loop diagrams” because all the contributions to the self-energy terms can be grouped into topologically-distinct loop diagrams based on the electron-phonon population topology originated by proper applications of KCRI’s and SDPO’s. This method can be applied further to other electron transition phenomena.