A Loop Diagram Approach to the Nonlinear Optical Conductivity for an Electron-Phonon System ()
1. Introduction
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.
2. Methods
For diagrammatic representation of the line shape functions,
, we consider the form
(1)
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:
(2)
(3)
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
(4)
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:
(5)
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.
3. Linear Optical Conductivity
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]
(6)
where
for electron states
and
, jk is the k component of the single electron current op-
Figure 1. The dotted line and the clockwise loop correspond to
and
, respectively.
Figure 2. The dotted line and the clockwise loop correspond to
and
, respectively.
erator, rl 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
(7)
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 ![](https://www.scirp.org/html/21-7500528\c0f86344-6edd-4001-8029-7bd43690e2d9.jpg)
The other three terms in
can be represented in a similar manner.
In Figure 4 or Equation (7), the first two terms are
Figure 4. The dotted straight line means that the implicit state λ is coupled with the initial state α by a phonon in the first diagram. Two implicit transitions occur between the implicit state λ and final state β, which forms a loop. The other three diagrams can be interpreted in a similar manner.
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.
4. First-Order Nonlinear Optical Conductivity
The first-order nonlinear optical conductivity for two incident electromagnetic waves of frequencies
and
is given by the following [10-12]
(8)
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:
(9)
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).
5. Concluding Remarks
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
Figure 6. The dotted straight line means that the implicit state λ is coupled with the initial state α by a phonon in the third diagram. There are two forward implicit transitions and two reverse implicit transitions centered on the implicit state λ, which form loops. The other three diagrams can be interpreted in a similar manner.
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.