ISRN Condensed Matter Physics

Volume 2014 (2014), Article ID 719120, 23 pages

http://dx.doi.org/10.1155/2014/719120

Research Article

Projection-Reduction Approach to Optical Conductivities for an Electron-Phonon System and Their Diagram Representation

Nam Lyong Kang¹ and Sang Don Choi²

¹Department of Applied Nanoscience, Pusan National University, Miryang 627-706, Republic of Korea

²Department of Physics, Kyungpook National University, Daegu 702-701, Republic of Korea

Received 7 October 2013; Accepted 16 January 2014; Published 7 April 2014

Academic Editors: D. L. Nika and T. Yanagisawa

Copyright © 2014 Nam Lyong Kang and Sang Don Choi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Utilizing state-dependent projection operators and the Kang-Choi reduction identities, we derive the linear, first, and second-order nonlinear optical conductivities for an electron system interacting with phonons. The lineshape functions included in the conductivity tensors satisfy “the population criterion” saying that the Fermi distribution functions for electrons and Planck distribution functions for phonons should be combined in multiplicative forms. The results also contain energy denominator factors enforcing the energy conservation as well as interaction factors describing electron-phonon interaction properly. Therefore, the phonon absorption and emission processes as well as photon absorption and emission processes in all electron transition processes can be presented in an organized manner and the results can be represented in diagrams that can model the quantum dynamics of electrons in a solid.

1. Introduction

Studies of the optical transitions for electron systems interacting with a background are useful for examining the electronic properties of solids because the absorption lineshapes are sensitive to the type of scattering mechanisms affecting the carrier behavior. One of the approaches to this problem is the projection method [1–17], which is generally based on the response theory. The projection techniques can be divided into two kinds: single-electron projection techniques [1–9] and many-electron projection techniques [10–17].

A proper many-body theory should satisfy the so-called “many-body criteria,” one of the most difficult being “the population criterion,” meaning that the Fermi distribution functions for electrons and the Planck distribution functions for phonons should be combined in multiplicative forms because the electrons and phonons belong to different categories in quantum-statistical physics. The population factors satisfying the population criterion should be combined in a physically acceptable manner with the energy denominator factor to maintain energy conservation and with the interaction factor to interpret the occurrence of local fluctuations due to electron-phonon interactions in such a way that the transition can take place via implicit states. Theories based on the single-electron projection technique starting from the Kubo formalism have failed to meet this criterion.

The many-electron projection techniques are classified into two categories; the first being state-independent or current-dependent [10–13] and the second being state-dependent [14–17]. The conductivity tensors [10, 11] produced using the isolation and projection operators, which are state-independent projection techniques, have limited applications for the following reasons. (a) They do not satisfy the population criterion because the two distribution functions are simply added (or subtracted), such as , where and are the Planck distribution function for a phonon with wave vector and the Fermi distribution function for an electron with energy , respectively. (b) They do not satisfy the “projection criterion,” so they are applicable only to cases where the spacing between the energy levels is constant, such as cyclotron resonance phenomena; that is, they are unsuitable for a system with nonuniform energy separation, such as a quantum square well (see the statements below (51)).

Since Franken et al. [18] pioneered the field of nonlinear optical phenomena by the experimental work on optical wave mixing, it has attracted considerable research attention [19, 20]. Most studies focused on the electronic structures and the electromagnetic properties of media with increasing the available laser power [20–33], but there were some problems. Shen [22] did not have explicit forms of the damping terms (relaxation rates) appearing in the conductivity. Suzuki and Ashikawa [28] reported the explicit forms of the damping term but did not satisfy the population criterion.

The present authors introduced an identity relation (the Kang-Choi reduction identity (KCRI) of the first kind) to derive the lineshape formula for a cyclotron transition [12]. The population factors appearing in this formula satisfy the population criterion (see [12, equation (27)]). This shows that correct use of the KCRI can satisfy the population criterion. Later, the identity was generalized to cover the higher-order nonlinear optical conductivities [13–16]. The second kind (see [13, equation (4.18)]) was used along with current-dependent projectors, which is called the state-independent projectors, to derive the first-order nonlinear optical conductivity [13]. This satisfies the population criterion but has limited applications because it does not satisfy the projection criterion. The corrected extended version of the nonlinear optical conductivity formalism, including the second-order nonlinear optical conductivity, was presented recently by applying the projection-reduction (PR) method by combining the generalized state-dependent projection operator (SDPO) and the KCRI of the 3rd kind (see [14, equation (2.5)]). The axes of the generalized projectors are the creation-annihilation operator pairs (see [14, equation (6.4)]), which makes the formalism applicable to general electron systems; that is, it satisfies the projection criterion.

The aim of this study is to generalize the KCRI to cover higher-order nonlinear optical conductivities and show how physically acceptable forms of the optical conductivities can be obtained by PR method. In addition, a method is introduced to represent the optical conductivity formulae using diagrams, through which a physical intuition to the quantum dynamics of the electrons in a solid can be achieved.

2. Hamiltonian and Density Operator

To derive the optical conductivities, the Hamiltonian is defined and the density operator is introduced.

2.1. Hamiltonian

We consider the total Hamiltonian, , for a system of electrons interacting with phonons as follows:

where the time-independent part of the system in thermodynamic equilibrium, , is given by

Here, , , and are the electron, phonon, and electron-phonon interaction Hamiltonians, respectively, is the creation (annihilation) operator for an electron in the state with energy , is that for a phonon in the state with energy , , is the phonon wave vector, is the polarization index, and is the electron-phonon interaction matrix element defined as

where is the coupling factor that depends on the mode of phonons and is the electron position vector.

In (1), when a linearly polarized time-dependent electric field with angular frequency , , is applied, the following time-dependent interaction Hamiltonian is considered:

where , , and denote the , , and components, respectively, of the electron position vector, for an arbitrary operator , and . means the complex conjugate.

2.2. Density Operator

To determine the general-order nonlinear optical conductivity, we assume that when the time-dependent electric field is applied to the system, the total density operator can be split into two parts as follows:

where is the density operator for a system in thermal equilibrium and is the perturbed term by the time-dependent external field. Using the Liouville equation

we can obtain

Here, we have used the fact that and is the Liouville operator corresponding to the total Hamiltonian defined as for an arbitrary operator, . can be split into two parts for calculation convenience: where and , respectively, correspond to and .

To obtain , we define the density operator in the Dirac picture as

By differentiating (8) and considering (7),

where we have used . Integrating (9) from to subject to the initial condition , a solution to this equation can be obtained:

By inserting (8) into (10) and replacing by ,

Finally, the following can be obtained by iteration:

where involves ’s times.

Using (12), the ensemble average of the th component of the many-electron current density operator can be determined as follows:

where , , and means the many-body trace. Equation (13) will be used in the following sections to calculate the linear and nonlinear optical conductivities.

3. Optical Conductivity

Using (12) and (13), the linear, first-, and second-order nonlinear optical conductivities are introduced.

3.1. Linear Optical Conductivity

The following can be obtained by inserting in (12) into and considering (4):

where we have used the cyclic property of trace, , and  (). Using , (14) can be written as follows:

Here, and the linear optical conductivity is defined as follows:

where

3.2. First-Order Nonlinear Optical Conductivity

Considering in (12) and , becomes

where and . Using , (18) becomes

Inserting (4) into (19),

where and the first-order nonlinear optical conductivity is defined as follows:

3.3. Second-Order Nonlinear Optical Conductivity

Inserting in (12) into ,

where , , and . Considering and inserting into (22),

where . Therefore, inserting , , and into (23) results in

Finally, considering (4), we obtain

where and the second-order nonlinear optical conductivity is defined as

4. Kang-Choi Reduction Identity

The essential part of this paper is expansion of the optical conductivities (see (16), (21) and (26)) to easily solvable forms. For this purpose, this section introduces various forms of a useful identity, known as the Kang-Choi reduction identity (KCRI).

4.1. KCRI of First Kind (KCRI-I)

For arbitrary operators, and , we consider

Considering and the trace property, , for the first and fourth terms, the KCRI-I of rank 1 can be obtained as follows:

which is useful because and are separated. Similarly, for arbitrary operators, , , and ,

Considering and for the first and eighth terms and adding

which is null, to the right-hand side of (29), the KCRI-I of rank 2 can be obtained as follows:

Therefore, the KCRI-I of rank [] can be inferred as follows:

where . The present form given by (32) is the general case. The proof for is as follows. The left-hand side (lhs) of (32) for can be expanded as follows:

which is the same as the right-hand side (rhs) of (32) for .

4.2. KCRI of Second Kind (KCRI-II)

The KCRI of the second kind (KCRI-II) is given by the following:

which will be used for derivation of the first-order nonlinear optical conductivity formula. Equation (34) can be proven by applying the identity

on ; that is, the lhs of (34) can be changed into

Applying with and gives the following:

Therefore, applying results in

which is reduced to the rhs of (34) by applying (35) in a reverse manner.

4.3. A Lemma Useful for Derivation of the KC Expansion Relation

A lemma useful to derive the general form of KCRI is introduced. For an analytic function of , ,

for . It suffices to prove the case of , because can be expanded in a Taylor series. The case for is dealt with for illustration. Using ,

Applying , the rhs of (40) becomes

Therefore, applying in a reverse manner results in

completing the proof for . The cases for and can be proven by an iterative operation of (40). Using this lemma, the calculations needed to prove the KCRI for are reduced.

4.4. KCRI of Third Kind (KCRI-III)

The KCRI of the 3rd kind (KCRI-III) can be expressed as

The proof can be obtained as follows. From (37), the lhs of (43) becomes

where the KC expansion relation (see (39)) was applied to the last step. Therefore, applying (37) in a reverse manner results in the right-hand side of (43). KCRI-III (see (43)) is needed in deriving the second-order nonlinear optical conductivity formula.

5. State-Dependent Projection Operators

In this section, we calculate the linear, first-, and second-order optical conductivities using the PR method.

5.1. Linear Optical Conductivity

Using the many-electron current density operator , which is written in terms of the single-electron current operator, (), as

the linear optical conductivity (16) becomes

where

To calculate , the first state-dependent projection operators for an arbitrary operator are introduced as follows:

where

Applying the identity, , to the right-hand side of the Liouville operator in (47) as and using (35) give

which can be rearranged as

Here, we have used and where , , and and are the Liouville operators corresponding to and , respectively.

Note that the projection operators given in (48) are state-dependent; that is, projects an arbitrary operator onto the operator , which depends on the states and . On the other hand, state-independent projection operators project into the state-independent current density operator, . The state-independent projection operator method is applicable only to the case in which or . Cyclotron phenomenon belongs to this category because for the current operator, , , where is the cyclotron frequency. This is called the “projection criterion.” Equation (48) satisfies this criterion because .

Taking the ensemble average of (51) gives

Here

which can be rearranged as

where we have used (35) and in the second term. To calculate (54) further, we make use of the following relation in the second term of (54):

where the first and second terms are canceled out. In this stage, is assumed and the terms up to the second order of are considered, assuming that the electron-phonon interaction is quite weak. Then, the first term in (54) and the fourth term in (55) are neglected because the ensemble averages of and are zeros. Therefore,

where

Equation (56) is an easily calculable form because it contains neither nor the projection operator in the denominator.

5.2. First-Order Nonlinear Optical Conductivity

By inserting (45) into (21), the first-order nonlinear optical conductivity becomes

where

which can be calculated by defining the second state-dependent projection operators for an arbitrary operator, , as follows:

where

Applying the identity, , to the right-hand side of the Liouville operator in (59) as and using the identity (35) give the following after a similar procedure as that used in (52):

Here

where

5.2.1. Calculation of

Applying (35) to (64) with and results in

in the second-order approximation in , where we have used the . Considering ,

where

5.2.2. Calculation of

Applying KCRI-II to (65) results in

Considering that

we obtain, in the second-order approximation in ,

where we have used (37) in the first term of (69),

Therefore, inserting (67) and (71) into (63) results in

in the second-order approximation in , where

Note that values in (74) are already approximated up to the second order of , so is approximated to .

5.3. Second-Order Nonlinear Optical Conductivity

By applying (45) to (26), the second-order nonlinear optical conductivity becomes

where

To calculate (77), the third state-dependent projection operators for an arbitrary operator, , are defined as follows:

where

Applying the identity, , to the right-hand side of the Liouville operator in (77) as and using (35) result in the following after a similar procedure as that used in (52) and (62):

Here

where

5.3.1. Calculation of

Considering (35) with and , (82) can be expanded as follows:

where

Applying (35) to again with and results in

in the second-order approximation in , where we have used the . Equation (87) can be calculated as follows:

where we have used and

KCRI-II is applied to to calculate , resulting in

in the second-order approximation in , where

Therefore, inserting (87) and (90) into (84) gives

5.3.2. Calculation of

Applying KCRI-III to (83) and using

result in

in the second-order approximation in , where

Finally, inserting (93) and (95) into (81) gives

in the second-order approximation in , where

6. Lineshape Functions and Diagram Representations

In this section, the lineshape functions for the electron-phonon system are derived.

6.1. Linear Term

In (52) and (56),

where is the Fermi distribution function for an electron in a state, , with energy . Inserting (101) into (52) gives

Therefore, (46) becomes

where the first lineshape function that contains the effect of the electron-phonon interaction is given as

which can be calculated as follows:

Here,

were used, where is the Planck distribution function for phonons with energy . In (105), the “energy denominator factor” is defined as

and the “population factors” are defined as

The denominator in (103) represents the transition from the initial state to the final state with a photon absorption of frequency and can be represented as Figure 1. If there is no phonon scattering, the lineshape of the conductivity would be like a delta function. On the other hand, the lineshape is broadened as electrons are scattered by phonons, so the lineshape function, , is included. Implicit states exist between the initial and final states because they are connected by , including the implicit states induced by electron-phonon interactions; that is, contains the implicit transitions between the initial (or final) and implicit states. An implicit state should be coupled with the initial (or final) state by a proper interaction factor in the case of an implicit transition between the final (or initial) and implicit states. must also contain an energy denominator factor enforcing energy conservation in each implicit transition and a population factor justifying the implicit transition.

Figure 1: Schematic diagram of the energy denominator, , in .

Although is called the lineshape function, the lineshape must be calculated from (103), not directly from (105). is represented by four diagrams in Figure 2. The two terms correspond to and the other two terms correspond to . is called the “interaction factor,” through which momentum conservation is satisfied.

(a)

(b)

(c)

(d)

Figure 2: Schematic diagrams of (a) , (b) , (c) , and (d) in .

The physical meaning of the first term in (105) is as follows. means that the implicit state is coupled with the initial state by a phonon with a wave vector . is represented by a vertical dotted line (Figure 2(a)). The states denoted by are called implicit states because they are included only in lineshape function (see (105)), not in the conductivity (see (103)). means that the energy of the implicit state is determined by the energy of the final state , photon energy , and phonon energy , so energy conservation can be satisfied; that is, . is represented by a horizontal solid line. The half circle with a rightward (leftward) arrow corresponds to photon absorption (emission) and the upper (lower) half circle corresponds to phonon emission (absorption). means that the reverse implicit transition from the final state to the implicit state with photon emission and phonon absorption () should be subtracted from the forward implicit transition from the implicit state to the final state with photon absorption and phonon emission (). This means that when an electron-phonon interaction is involved in an electron transition, the electron undergoes an implicit transition from an implicit state, which is coupled with the initial state, to the final state (vice versa) with photon absorption (emission) and phonon emission (absorption). The second term in (105) corresponds to phonon absorption (emission) for a forward (backward) process. The third and fourth terms in (105) correspond to the cases that the implicit state is coupled with the final state. Therefore, the transition from the initial state to the final state occurs via two implicit transitions and the implicit state is coupled with the initial or final states. We can represent the result in the diagram because the lineshape function contains an energy denominator, population, and interaction factors properly.

6.2. First-Order Nonlinear Term

In (62), can be calculated as

Considering (102), (109) becomes

On the other hand, (75) can be calculated as follows:

where

was used. Therefore, inserting (110) into (62) and using (74) and (111),

where

Here, we have neglected the values because we are interested in the results near the resonance points, and . Inserting (113) into (58) gives (see appendices (A.1) and (A.2) for details)

Considering (A.1) and (A.2) in the appendix, the lineshape functions are calculated as follows:

and in (115) can be represented by Figure 3. An intermediate state, , exists because of the two photons processes. Therefore, the first-order nonlinear conductivity consists of two parts. The first one, called a () process, is , which contains a photon absorption of frequency in an intermediate transition from the initial state to the intermediate state (Figure 3(a)). The other one, called a () process, is , which is an exchange process with and being exchanged; that is, it contains a photon absorption of frequency in an intermediate transition from the intermediate state to the final state, (Figure 3(b)).

(a)

(b)

Figure 3: Schematic diagrams of the energy denominators in (a) and (b) .

The first term, , contains two transition processes. The first is an intermediate transition from the initial state, , to the intermediate state, , with a photon absorption of frequency and lineshape function, , and the other is a direct transition from the initial state to the final state, , with two photon absorption of frequencies and and lineshape function . Note that . The dotted half circle in Figure 3(a) means that the state is not connected to the state by any lineshape function. Therefore, no implicit transition exists between the and states. The second term, , can be interpreted in a similar manner, where the state is not connected to the state by any lineshape function. Figures 4 and 5 show schematic diagrams of the lineshape functions and , respectively. The linear part included in has the same form as (105).

(a)

(b)

(c)

(d)

Figure 4: Schematic diagrams of (a) , (b) , (c) , and (d) in .

(a)

(b)

(c)

(d)

Figure 5: Schematic diagrams of (a) , (b) , (c) , and (d) in .

Equation (116) can be interpreted as follows (Figure 4). Note that an intermediate state, , exists between the initial state, , and final state, , in the first term of (115) (Figure 3(a)). With reasoning similar to that for the linear conductivity, there are four diagrams; two terms corresponding to and the other two terms corresponding to . The first term in (116) (Figure 4(a)) corresponds to a transition from the implicit state to the intermediate state with a photon absorption of frequency . (horizontal solid line) represents the intermediate transition and (vertical dotted line) means that the implicit state is coupled with the initial state. The third term in (116) (Figure 4(c)) corresponds to a transition from the initial state to the implicit state with two photon absorption of frequencies and . represents the transition and means that the implicit state is coupled with the final state. Note that (dotted cicle) should be subtracted from because the intermediate state is not connected to the final state by any lineshape function (see Figure 3(a)). The second and fourth terms (Figures 4(b) and 4(d)) in (116) and the four terms in (117) (Figure 5) can be interpreted in a similar manner. Note that the implicit state in (117) is not connected to the initial state by any lineshape function (see Figure 3(b)).

6.3. Second-Order Nonlinear Term

In (80), can be calculated as follows:

Considering (113), (118) becomes

On the other hand, (100) can be calculated as

where we have used (112). Applying (112) to (120) gives

Therefore, inserting (119) into (80) and using (99) and (121) yield

where the four lineshape functions are given as follows:

Here, , , and values were neglected because the resonance phenomena are considered.

Finally, inserting (122) into (76) gives (see appendices (A.3), (A.4), and (A.5) for details)

where, considering (A.3)–(A.5) in the appendix, the lineshape functions are calculated as follows:

The physical meaning of (124) is clear. Two intermediate states exist because there are three photon processes, and thus, four terms in the second-order nonlinear conductivity exist (Figure 6). With reasoning similar to that for the first-order nonlinear conductivity, there are and processes, which consist of two processes each. The first term, , contains a photon absorption process of frequency and the process (Figure 6(a)). The third term, , contains the process exchanged with a photon absorption process of frequency (Figure 6(c)). The second and fourth terms, and , contain the process (Figures 6(b) and 6(d)).

(a)

(b)

(c)

(d)

Figure 6: Schematic diagrams of the energy denominators in (a) , (b) , (c) , and (d) .

The first term, , contains three transition processes (Figure 6(a)). The first transition is an intermediate transition from the initial state to an intermediate state , with a photon absorption of frequency and lineshape function . The second transition is an intermediate transition from the initial state to an intermediate state , with two photons absorption of frequencies and and lineshape function . The third transition is a direct transition from the initial state to the final state , with three photon absorption of frequencies , , and and lineshape function . Note that . The dotted half circle in Figure 6(a) means that the state is not connected to the and states by any lineshape function. On the other hand, the and states in the second term, , are not connected to the and states, respectively, by any lineshape function. The other terms in (124) can be interpreted in a similar manner. Figures 7, 8, 9, and 10 show schematic diagrams of the four lineshape functions (see (125)–(128)) in (124). The linear part, , included in has the same form as (105) and the first-order nonlinear parts, , has the same forms as (116) or (117).

(a)

(b)

(c)

(d)

Figure 7: Schematic diagrams of (a) , (b) , (c) , and (d) in .

(a)

(b)

(c)

(d)

Figure 8: Schematic diagrams of (a) , (b) , (c)  , and (d)   in .

(a)

(b)

(c)

(d)

Figure 9: Schematic diagrams of (a) , (b) , (c)  , and (d)   in .

(a)

(b)

(c)

(d)

Figure 10: Schematic diagrams of (a) , (b) , (c) , and (d) in .

Equation (125) (Figure 7) can be interpreted by (116) and (126) can be interpreted as follows (Figure 8). Like the linear and first-order nonlinear conductivities, there are four diagrams. The first two terms correspond to and the other two terms correspond to . The first term (Figure 8(a)) represents a transition from the implicit state to an intermediate state , with two photon absorption of frequencies and . (horizontal solid line) represents the transition and (vertical dotted line) means that the implicit state is coupled with the initial state, . Note that (dotted circle) should be subtracted from because the intermediate state, , is not connected to the initial state by any lineshape function (see Figure 6(b)) in process. The third term (Figure 8(c)) corresponds to a transition from the initial state to the final state, , with three photon absorption of frequencies , , and . Note that the intermediate state, , is not connected to the final state by any lineshape function (see Figure 6(b)). Therefore, should be subtracted from , not from , because the state is also not connected to the initial state, , by any lineshape function (see Figure 6(b)). Equations (127) and (128) can be interpreted in a similar manner (Figures 9 and 10).

7. Concluding Remarks

So far, we presented a method for deriving the nonlinear optical conductivity tensors for a system of electrons interacting with phonons using the PR method. The standard inductive inference will be required for a more rigorous proof of the general form of KCRI. The use of this method in the present problem, however, was beyond the scope of this study because the general-order term contains a large number of elementary factors. Consequently, this study used a simplified version of the induction instead, which states that the proposition KCRI of the th kind is valid for all positive integers, , if it holds for = 1, 2, and 3. This method is justifiable because the structure of KCRI and the population factor are iterative.

Nonlinear optics is concerned with the phenomena resulting from light-induced changes in the optical properties of a system, the first-order nonlinear part being the simplest example. This study considers the second harmonic generation (SHG) as a proper example, in which the second parts of and play an important role with . Note that a small fraction is converted to light at if an incident beam at frequency is focused on a proper sample. The two fractions can be filtered technically using a prism. Although there have been many studies on SHG in systems of electrons interacting with phonons [34–36], less attention has been paid to how the absorption and emission of photons and phonons occur in each process. The present work aims to clear up the physical cause and effect of the processes. These results will be improved if radiation-induced inharmonic phonon modes are included in the consideration. This part with a derivation of the high-order nonlinear term will be examined in future.

The projection operators used in this study are generally applicable because they satisfy the projection criterion (see statements below (51)). The lineshape functions in the linear, first-, and second-order nonlinear optical conductivities include the electron distribution functions to satisfy the population criterion as well as energy denominator factors enforcing energy conservation and the interaction factors coupling the implicit state with the initial, final, or intermediate states. Therefore, the phonon emissions and absorptions as well as photon absorptions and emissions in all electron transition processes can be presented in an organized manner and the results can be represented in diagrams, through which a physical intuition to quantum dynamics of electron transition in a solid can be obtained. These diagrams are quite different from the temperature diagrams with imaginary times in the Feynman-like scheme [37–42]. In the present scheme, the contribution of each term to the lineshape function (or the damping term) is the sum of all weighted topologically distinct diagrams. Therefore, the present scheme can be categorized as another Feynman-like scheme. In conclusion, we claim that a theory, which is generally applicable and satisfies the population criterion, can be obtained using the PR method. Formulations of other electron transition phenomena might be obtained using this method.

Appendix

Calculation of Values

Using (106) properly, (68), (72), (89), (91), and (96) can be calculated. Here, we neglect the terms including and because we consider only the electron-phonon interaction given by (2) and (3):

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

R. Kubo, “Statistical-mechanical theory of irreversible processes. I: general theory and simple applications to magnetic and conduction problems,” Journal of the Physical Society of Japan, vol. 12, no. 6, pp. 570–586, 1957.

R. Zwanzig, “Ensemble method in the theory of irreversibility,” The Journal of Chemical Physics, vol. 33, no. 5, pp. 1338–1341, 1960.

H. Mori, “A continued-fraction representation of the time-vorrelation functions,” Progress of Theoretical Physics, vol. 34, no. 3, pp. 399–416, 1965.

A. Kawabata, “Theory of cyclotron resonance line width,” Journal of the Physical Society of Japan, vol. 23, no. 5, pp. 999–1006, 1967.

P. N. Argyres and J. L. Sigel, “Projection techniques in the theory of electrical resistivity,” Physical Review Letters, vol. 31, no. 23, pp. 1397–1400, 1973.

S. D. Choi and O. H. Chung, “Theory of cyclotron resonance lineshape for electron-phonon systems in two coupling schemes,” Solid State Communications, vol. 46, no. 9, pp. 717–720, 1983.

S. N. Yi, O. H. Chung, and S. D. Choi, “Theory of cyclotron resonance lineshape revisited,” Progress of Theoretical Physics, vol. 77, no. 4, pp. 891–898, 1987.

J. Y. Ryu and S. D. Choi, “Quantum-statistical theory of high-field transport phenomena,” Physical Review B, vol. 44, no. 20, pp. 11328–11338, 1991.

J. Y. Sug and S. D. Choi, “Quantum transport theory based on the equilibrium density projection technique,” Physical Review E, vol. 55, no. 1, pp. 314–321, 1997.

S. Badjou and P. N. Argyres, “Theory of cyclotron resonance in an electron-phonon system,” Physical Review B, vol. 35, no. 12, pp. 5964–5968, 1987.

Y. J. Cho and S. D. Choi, “Theory of cyclotron-resonance line shapes based on the isolation-projection technique,” Physical Review B, vol. 47, no. 15, pp. 9273–9278, 1993.

N. L. Kang, Y. J. Cho, and S. D. Choi, “A many-body theory of quantum-limit cyclotron transition line-shapes in electron-phonon systems based on projection technique,” Progress of Theoretical Physics, vol. 96, no. 2, pp. 307–316, 1996.

H. J. Lee, N. L. Kang, J. Y. Sug, and S. D. Choi, “Calculation of the nonlinear optical conductivity by a quantum-statistical method,” Physical Review B, vol. 65, no. 19, Article ID 195113, 2002.

N. L. Kang and S. D. Choi, “Derivation of nonlinear optical conductivity by using a reduction identity and a state-dependent projection method,” Journal of Physics A, vol. 43, no. 16, Article ID 165203, 2010.

N. L. Kang, Y. J. Lee, and S. D. Choi, “Generalization of the nonlinear optical conductivity formalism by using a projection-reduction method,” Journal of the Korean Physical Society, vol. 58, no. 3, pp. 538–544, 2011.

N. L. Kang and S. D. Choi, “Validity of nth-order conductivity by the general projection-reduction method,” Progress of Theoretical Physics, vol. 125, no. 5, pp. 1011–1019, 2011.

N. L. Kang and S. D. Choi, “Derivation of second-order nonlinear optical conductivity by the projection diagram method,” AIP Advances, vol. 2, Article ID 012161, 2012.

P. A. Franken, A. E. Hill, C. W. Peters, and G. Weinreich, “Generation of optical harmonics,” Physical Review Letters, vol. 7, no. 4, pp. 118–119, 1961.

J. A. Armstrong, N. Bloembergen, J. Ducuing, and P. S. Pershan, “Interactions between light waves in a nonlinear dielectric,” Physical Review, vol. 127, no. 6, pp. 1918–1939, 1962.

N. Bloembergen and Y. R. Shen, “Quantum-theoretical comparison of nonlinear susceptibilities in parametric media, lasers, and Raman lasers,” Physical Review, vol. 133, no. 1, pp. A37–A49, 1964.

L. V. Keldysh, “Ionization in the Field of a strong electromagnetic wave,” Soviet Physics JETP, vol. 20, no. 5, pp. 1307–1314, 1965.

Y. R. Shen, The Principles of Nonlinear Optics, chapter 2, Willey, New York, USA, 1984.

B. S. Wherrett, “Scaling rules for multiphoton interband absorption in semiconductors,” Journal of the Optical Society of America B, vol. 1, no. 1, pp. 67–72, 1984.

M. Sheik-Bahae, D. C. Hutchings, D. J. Hagan, and E. W. Van Stryland, “Dispersion of bound electronic nonlinear refraction in solids,” IEEE Journal of Quantum Electronics, vol. 27, no. 6, pp. 1296–1309, 1991.

W. S. Fann, R. Storz, H. W. K. Tom, and J. Bokor, “Electron thermalization in gold,” Physical Review B, vol. 46, no. 20, pp. 13592–13595, 1992.

C.-K. Sun, F. Vallée, L. H. Acioli, E. P. Ippen, and J. G. Fujimoto, “Femtosecond-tunable measurement of electron thermalization in gold,” Physical Review B, vol. 50, no. 20, pp. 15337–15348, 1994.

R. H. M. Groeneveld, R. Sprik, and A. Lagendijk, “Femtosecond spectroscopy of electron-electron and electron-phonon energy relaxation in Ag and Au,” Physical Review B, vol. 51, no. 17, pp. 11433–11445, 1995.

A. Suzuki and M. Ashikawa, “Quantum-statistical theory of nonlinear optical conductivity for an electron-phonon system,” Physical Review E, vol. 58, no. 4, pp. 4307–4320, 1998.

N. Del Fatti, C. Voisin, M. Achenmann, S. Tzortzakis, D. Christofilos, and F. Vallée, “Nonequilibrium electron dynamics in noble metals,” Physical Review B, vol. 61, no. 24, pp. 16956–16966, 2000.

Y. Hamanaka, J. Kuwabata, I. Tanahashi, S. Omi, and A. Nakamura, “Ultrafast electron relaxation via breathing vibration of gold nanocrystals embedded in a dielectric medium,” Physical Review B, vol. 63, no. 10, Article ID 104302, 2001.

A. Zarifi, C. Fisker, and T. G. Pedersen, “Theoretical study of quadratic electro-optic effect in semiconducting zigzag carbon nanotubes,” Physical Review B, vol. 76, no. 4, Article ID 045403, 2007.

N. Olivier and E. Beaurepaire, “Third-harmonic generation microscopy with focus-engineered beams: a numerical study,” Optics Express, vol. 16, no. 19, pp. 14703–14708, 2008.

F. Munhoz, H. Rigneault, and S. Brasselet, “High order symmetry structural properties of vibrational resonances using multiple-field polarization coherent anti-stokes Raman spectroscopy microscopy,” Physical Review Letters, vol. 105, no. 12, Article ID 123903, 2010.

C.-C. Wu and J. Tsai, “Coupling mechanisms of acoustic second-harmonic generation in piezoelectric semiconductor,” Journal of Physics C, vol. 20, no. 10, pp. 1527–1533, 1987.

Y.-M. Chang and N.-A. Chang, “Observation of coherent optical phonons in GaAs/AlxGa1-xAs single heterojunctions,” Journal of Applied Physics, vol. 93, no. 4, pp. 2015–2018, 2003.

D. H. Son, J. S. Wittenberg, U. Banin, and A. P. Alivisatos, “Second harmonic generation and confined acoustic phonons in highly excited semiconductor nanocrystals,” Journal of Physical Chemistry B, vol. 110, no. 40, pp. 19884–19890, 2006.

R. P. Feynman, “Space-time approach to quantum electrodynamics,” Physical Review, vol. 76, no. 6, pp. 769–789, 1949.

A. L. Fetter and J. D. Walecka, Quantum Theory of Many Particle Systems, McGraw Hill,, New York, NY, USA, 1971.

R. D. Mattuch, A Guide to Feynman Diagrams in the Many Body Problem, McGraw Hill, New York, NY, USA, 2nd edition, 1976.

G. D. Mahan, Many Particle Physics, chapter 3, Plenum, New York, NY, USA, 2000.

F. Giustino, M. L. Cohen, and S. G. Louie, “Small phonon contribution to the photoemission kink in the copper oxide superconductors,” Nature, vol. 452, no. 7190, pp. 975–978, 2008.

A. R. Wright and C. Zhang, “Dynamic conductivity of graphene with electron-LO-phonon interaction,” Physical Review B, vol. 81, no. 16, Article ID 165413, 2010.