# Semiconductor Bloch equations

The semiconductor Bloch equations (abbreviated as SBEs) describe the optical response of semiconductors excited by coherent classical light sources, such as lasers. They are based on a full quantum theory, and form a closed set of integro-differential equations for the quantum dynamics of microscopic polarization and charge carrier distribution. The SBEs are named after the structural analogy to the optical Bloch equations that describe the excitation dynamics in a two-level atom interacting with a classical electromagnetic field. As the major complication beyond the atomic approach, the SBEs must address the many-body interactions resulting from Coulomb force among charges and the coupling among lattice vibrations and electrons.

## Background

The optical response of a semiconductor follows if one can determine its macroscopic polarization $\mathbf {P}$ as a function of the electric field $\mathbf {E}$ that excites it. The connection between $\mathbf {P}$ and the microscopic polarization $P_{\mathbf {k} }$ is given by
$\mathbf {P} =\mathbf {d} \,\sum _{\mathbf {k} }P_{\mathbf {k} }+\operatorname {c.c.} \;,$
where the sum involves crystal-momenta $\hbar {\mathbf {k} }$ of all relevant electronic states. In semiconductor optics, one typically excites transitions between a valence and a conduction band. In this connection, $\mathbf {d}$ is the dipole matrix element between the conduction and valence band and $P_{\mathbf {k} }$ defines the corresponding transition amplitude.
The derivation of the SBEs starts from a system Hamiltonian that fully includes the free-particles, Coulomb interaction, dipole interaction between classical light and electronic states, as well as the phonon contributions. Like almost always in many-body physics, it is most convenient to apply the second-quantization formalism after the appropriate system Hamiltonian ${\hat {H}}_{\mathrm {System} }$ is identified. One can then derive the quantum dynamics of relevant observables ${\hat {\mathcal {O}}}$ by using the Heisenberg equation of motion
$\mathrm {i} \hbar {\frac {\mathrm {d} }{\mathrm {d} t}}\langle {\hat {\mathcal {O}}}\rangle =\langle [{\hat {\mathcal {O}}},{\hat {H}}_{\mathrm {System} }]_{-}\rangle \;.$
Due to the many-body interactions within ${\hat {H}}_{\mathrm {System} }$, the dynamics of the observable ${\hat {\mathcal {O}}}$ couples to new observables and the equation structure cannot be closed. This is the well-known BBGKY hierarchy problem that can be systematically truncated with different methods such as the cluster-expansion approach.
At operator level, the microscopic polarization is defined by an expectation value for a single electronic transition between a valence and a conduction band. In second quantization, conduction-band electrons are defined by fermionic creation and annihilation operators ${\hat {a}}_{c,{\mathbf {k} }}^{\dagger }$ and ${\hat {a}}_{c,{\mathbf {k} }}$, respectively. An analogous identification, i.e., ${\hat {a}}_{v,{\mathbf {k} }}^{\dagger }$ and ${\hat {a}}_{v,{\mathbf {k} }}$, is made for the valence band electrons. The corresponding electronic interband transition then becomes
$P_{\mathbf {k} }^{\star }=\langle {\hat {a}}_{c,\mathbf {k} }^{\dagger }{\hat {a}}_{v,\mathbf {k} }\rangle \,,\qquad P_{\mathbf {k} }=\langle {\hat {a}}_{v,\mathbf {k} }^{\dagger }{\hat {a}}_{c,\mathbf {k} }\rangle \,,$
that describe transition amplitudes for moving an electron from conduction to valence band ($P_{\mathbf {k} }^{\star }$ term) or vice versa ($P_{\mathbf {k} }$ term). At the same time, an electron distribution follows from
$f_{\mathbf {k} }^{e}=\langle {\hat {a}}_{c,\mathbf {k} }^{\dagger }{\hat {a}}_{c,\mathbf {k} }\rangle \;.$
It is also convenient to follow the distribution of electronic vacancies, i.e., the holes,
$f_{\mathbf {k} }^{h}=1-\langle {\hat {a}}_{v,\mathbf {k} }^{\dagger }{\hat {a}}_{v,\mathbf {k} }\rangle =\langle {\hat {a}}_{v,\mathbf {k} }{\hat {a}}_{v,\mathbf {k} }^{\dagger }\rangle$
that are left to the valence band due to optical excitation processes.

## Principal structure of SBEs

The quantum dynamics of optical excitations yields an integro-differential equations that constitute the SBEs
$\mathrm {i} \hbar {\frac {\partial }{\partial t}}P_{\mathbf {k} }={\tilde {\varepsilon }}_{\mathbf {k} }P_{\mathbf {k} }-\left[1-f_{\mathbf {k} }^{e}(t)-f_{\mathbf {k} }^{h}(t)\right]\Omega _{\mathbf {k} }+\mathrm {i} \hbar \left.{\frac {\partial }{\partial t}}P_{\mathbf {k} }\right|_{\mathrm {scatter} }\,,$
$\hbar {\frac {\partial }{\partial t}}f_{\mathbf {k} }^{e}=2\operatorname {Im} \left[\Omega _{\mathbf {k} }^{\star }P_{\mathbf {k} }\right]+\hbar \left.{\frac {\partial }{\partial t}}f_{\mathbf {k} }^{e}\right|_{\mathrm {scatter} }\,,$
$\hbar {\frac {\partial }{\partial t}}f_{\mathbf {k} }^{h}=2\operatorname {Im} \left[\Omega _{\mathbf {k} }^{\star }P_{\mathbf {k} }\right]+\hbar \left.{\frac {\partial }{\partial t}}f_{\mathbf {k} }^{h}\right|_{\mathrm {scatter} }\;.$
These contain the renormalized Rabi energy
$\Omega _{\mathbf {k} }=\mathbf {d} \cdot \mathbf {E} +\sum _{\mathbf {k} '\neq \mathbf {k} }V_{\mathbf {k} -\mathbf {k} '}P_{\mathbf {k} '}$
as well as the renormalized carrier energy
${\tilde {\varepsilon }}_{\mathbf {k} }=\varepsilon _{\mathbf {k} }-\sum _{\mathbf {k} '\neq \mathbf {k} }V_{\mathbf {k} -\mathbf {k} '}\left[f_{\mathbf {k} '}^{e}+f_{\mathbf {k} '}^{h}\right]\,,$
where $\varepsilon _{\mathbf {k} }$ corresponds to the energy of free electron–hole pairs and $V_{\mathbf {k} }$ is the Coulomb matrix element, given here in terms of the carrier wave vector $\mathbf {k}$.
The symbolically denoted $\left.\cdots \right|_{\mathrm {scatter} }$ contributions stem from the hierarchical coupling due to many-body interactions. Conceptually, $P_{\mathbf {k} }$, $f_{\mathbf {k} }^{e}$, and $f_{\mathbf {k} }^{h}$ are single-particle expectation values while the hierarchical coupling originates from two-particle correlations such as polarization-density correlations or polarization-phonon correlations. Physically, these two-particle correlations introduce several nontrivial effects such as screening of Coulomb interaction, Boltzmann-type scattering of $f_{\mathbf {k} }^{e}$ and $f_{\mathbf {k} }^{h}$ toward Fermi–Dirac distribution, excitation-induced dephasing, and further renormalization of energies due to correlations.
All these correlation effects can be systematically included by solving also the dynamics of two-particle correlations. At this level of sophistication, one can use the SBEs to predict optical response of semiconductors without phenomenological parameters, which gives the SBEs a very high degree of predictability. Indeed, one can use the SBEs in order to predict suitable laser designs through the accurate knowledge they produce about the semiconductor's gain spectrum. One can even use the SBEs to deduce existence of correlations, such as bound excitons, from quantitative measurements.
The presented SBEs are formulated in the momentum space since carrier's crystal momentum follows from $\hbar \mathbf {k}$. An equivalent set of equations can also be formulated in position space. However, especially, the correlation computations are much simpler to be performed in the momentum space.

## Interpretation and consequences

The $P_{\mathbf {k} }$ dynamic shows a structure where an individual $P_{\mathbf {k} }$ is coupled to all other microscopic polarizations due to the Coulomb interaction $V_{\mathbf {k} }$. Therefore, the transition amplitude $P_{\mathbf {k} }$ is collectively modified by the presence of other transition amplitudes. Only if one sets $V_{\mathbf {k} }$ to zero, one finds isolated transitions within each ${\mathbf {k} }$ state that follow exactly the same dynamics as the optical Bloch equations predict. Therefore, already the Coulomb interaction among $P_{\mathbf {k} }$ produces a new solid-state effect compared with optical transitions in simple atoms.
Conceptually, $P_{\mathbf {k} }$ is just a transition amplitude for exciting an electron from valence to conduction band. At the same time, the homogeneous part of $P_{\mathbf {k} }$ dynamics yields an eigenvalue problem that can be expressed through the generalized Wannier equation. The eigenstates of the Wannier equation is analogous to bound solutions of the hydrogen problem of quantum mechanics. These are often referred to as exciton solutions and they formally describe Coulombic binding by oppositely charged electrons and holes.
However, a real exciton is a true two-particle correlation because one must then have a correlation between one electron to another hole. Therefore, the appearance of exciton resonances in the polarization does not signify the presence of excitons because $P_{\mathbf {k} }$ is a single-particle transition amplitude. The excitonic resonances are a direct consequence of Coulomb coupling among all transitions possible in the system. In other words, the single-particle transitions themselves are influenced by Coulomb interaction making it possible to detect exciton resonance in optical response even when true excitons are not present.
Therefore, it is often customary to specify optical resonances as excitonic instead of exciton resonances. The actual role of excitons on optical response can only be deduced by quantitative changes to induce to the linewidth and energy shift of excitonic resonances.
The solutions of the Wannier equation produce valuable insight to the basic properties of a semiconductor's optical response. In particular, one can solve the steady-state solutions of the SBEs to predict optical absorption spectrum analytically with the so-called Elliott formula. In this form, one can verify that an unexcited semiconductor shows several excitonic absorption resonances well below the fundamental bandgap energy. Obviously, this situation cannot be probing excitons because the initial many-body system does not contain electrons and holes to begin with. Furthermore, the probing can, in principle, be performed so gently that one essentially does not excite electron–hole pairs. This gedanken experiment illustrates nicely why one can detect excitonic resonances without having excitons in the system, all due to virtue of Coulomb coupling among transition amplitudes.

## Extensions

The SBEs are particularly useful when solving the light propagation through a semiconductor structure. In this case, one needs to solve the SBEs together with the Maxwell's equations driven by the optical polarization. This self-consistent set is called the Maxwell–SBEs and is frequently applied to analyze present-day experiments and to simulate device designs.
At this level, the SBEs provide an extremely versatile method that describes linear as well as nonlinear phenomena such as excitonic effects, propagation effects, semiconductor microcavity effects, four-wave-mixing, polaritons in semiconductor microcavities, gain spectroscopy, and so on. One can also generalize the SBEs by including excitation with terahertz (THz) fields that are typically resonant with intraband transitions. One can also quantize the light field and investigate quantum-optical effects that result. In this situation, the SBEs become coupled to the semiconductor luminescence equations.