A phonovoltaic (pV) cell converts vibrational (phonons) energy into a direct current much like the photovoltaic effect in a photovoltaic (PV) cell converts light (photon) into power. That is, it uses a p-n junction to separate the electrons and holes generated as valence electrons absorb optical phonons more energetic than the band gap, and then collects them in the metallic contacts for use in a circuit. The pV cell is an application of heat transfer physics and competes with other thermal energy harvesting devices like the thermoelectric generator. While the thermoelectric generator converts heat, a broad spectrum of phonon and electron energy, to electricity, the pV cell converts only a narrow band of phonon energy, i.e., only the most energetic optical phonon modes. A narrow band of excited optical phonons has much less entropy than heat. Thus, the pV cell can exceed the thermoelectric efficiency. However, exciting and harvesting the optical phonon poses a challenge.
Satisfying the laws of thermodynamics By the first law of thermodynamics, the excitation driving electron generation in both photo- and phonovoltaic cells, i.e., the photon or phonon, must have more energy than the semiconductor band gap. For a PV cell, many materials are available with a band gap ( Δ E e , g {\displaystyle \Delta E_{e,g}} ) well matched to the solar photon spectrum, like Silicon or Gallium Arsenide. For a pV cell, however, no current semiconducting materials have a band gap smaller than the energy of their most energetic (optical) phonon modes ( E p , O {\displaystyle E_{p,\mathrm {O} }} ). Thus, novel materials are required with both energetic optical phonon modes ( E p , O ≥ 100 {\displaystyle E_{p,\mathrm {O} }\geq 100} meV, e.g., graphene, diamond, or boron nitride) and a small band gap ( Δ E e , g < E p , O {\displaystyle \Delta E_{e,g}<E_{p,\mathrm {O} }} , e.g., graphene). By the second law of thermodynamics, the excitation must be "hotter" than the cell for power generation to occur. In a PV, the light comes from an outside source, for example, the sun, which is nearly 6000 kelvins, whereas the PV is around 300 kelvins. Thus, the second law is satisfied and energy conversion is possible. However, the crystal vibrations driving power generation in a pV are intrinsic to the material itself. As such, they can not be imported from an outside source like the sun, but must instead be excited by some other process until they are hotter than the cell. The temperature of the optical phonon population is calculated by comparing the number of optical phonon to the number expected at a given temperature, which comes from the Bose–Einstein statistics.
Non-equilibrium optical phonon population and the nanoscale requirement
There are a number of ways to excite a population of vibrations, i.e., create a hot optical phonon population. For example, if the electron population is excited, using a laser or electric field, they will typically relax by emitting optical phonons. Additionally, a hot molecular gas can impart its vibrations to a crystal when chemisorbed. Regardless of method, the conversion efficiency is limited by the optical phonon temperature achieved as compared to the electron temperature within the device due to Carnot's theorem. In a nanoscale device, this temperature is approximately equal to the temperature of the device itself. However, in a macroscale device the generated electrons accumulate faster than they are collected. Thus, the electron population is heated up to the optical phonon temperature and further generation is inhibited. The down-conversion is simultaneously inhibited as the acoustic phonon population is heated to the optical phonon temperature. Thus, the large pV cell develops a near-equilibrium state where it is heated. At best, it will act like a thermoelectric generator and exhibit thermoelectric effects. Such a device is called a thermovoltaic, rather than a phonovoltaic.
Entropy generation and efficiency Entropy generation and inefficiency in a PV cell is the result of photons more energetic than the band gap producing electrons with kinetic energy in addition to the potential energy provided by the band gap. Similarly, optical phonon energy in excess of the band gap generates an entropy flow in the pV cell, rather than electric power. The energy efficiency ( η ϕ {\displaystyle \eta _{\phi }} ) is quantified by the ratio of the band gap and optical phonon energy, that is
In addition this typical inefficiency, hot optical phonon populations tend to downconvert into multiple low-energy, acoustic phonon modes (whereas photons typically do not downconvert into low energy infrared waves). This efficiency ( η Q E {\displaystyle \eta _{QE}} ) is quantified by the tendency of a hot optical phonon to downconvert rather than generate an electron-hole pair, that is
where γ ˙ e − p {\displaystyle {\dot {\gamma }}_{e-p}} is the rate of generation and γ ˙ p − p {\displaystyle {\dot {\gamma }}_{p-p}} is the rate of downconversion, i.e., the rate at which an optical phonon produces multiple low-energy, acoustic phonons. This provides a second entropy flow reducing the efficiency of a pV cell. Finally, entropy is generated in both pV and PV cells due to the inefficient separation of the generated electrons and holes. This efficiency ( η p n {\displaystyle \eta _{pn}} ) is limited by the Carnot efficiency given by
where T p V {\displaystyle T_{\mathrm {pV} }} is the temperature of the pV cell and T p , O {\displaystyle T_{p,\mathrm {O} }} is the temperature of the optical phonon population, as dictated by the Bose–Einstein statistics. This efficiency is reduced the smaller the band gap is in comparison to the thermal energy ( k B T {\displaystyle k_{\mathrm {B} }T} , where k B {\displaystyle k_{\mathrm {B} }} is the Boltzmann constant and T {\displaystyle T} is the temperature). Indeed, the p-n junction efficiency is approximately
Thus, the overall efficiency ( η p V {\displaystyle \eta _{\mathrm {pV} }} ) is
where the temperature independent terms become the material figure of merit ( Z p V {\displaystyle Z_{\mathrm {pV} }} ),
If the band gap and optical phonon mode are resonant, and the optical phonon tends to generate electrons, the phonovoltaic cell can approach the Carnot limit as T p v → 0 {\displaystyle T_{pv}\rightarrow 0} .
The electron-phonon coupling
The electron-phonon coupling is responsible for electron generation in the pV cell. In this phenomenon, the phonon leads to ion motion which perturbs the highest occupied valence state (HOS). This state begins to overlap with the lowest unoccupied conduction state (LUS), and the electron can switch states if energy and momentum are conserved. If it does, an electron-hole pair is generated. Using a taylor expansion of the change in electron potential, φ {\displaystyle \varphi } , due to the ionic displacement of a phonon provides a matrix element for use in Fermi's golden rule, and the derivation of a generation rate. This Taylor expansion gives the following matrix element
⟨ f | H e − p ′ | i ⟩ = M e − p = ( ℏ 2 ⟨ m ⟩ ω κ p , α ) 1 / 2 ⟨ κ e + κ p , j | ∂ φ ∂ d κ p , α | κ e , i ⟩ , {\displaystyle {\begin{aligned}\langle f|\mathrm {H} _{e-p}'|i\rangle =M_{e-p}=({\frac {\hbar }{2\langle m\rangle \omega _{\kappa _{p},\alpha }}})^{1/2}\langle {\boldsymbol {\kappa }}_{e}+{\boldsymbol {\kappa }}_{p},j|{\frac {\partial \varphi }{\partial \mathbf {d} _{{\boldsymbol {\kappa }}_{p},\alpha }}}|{\boldsymbol {\kappa }}_{e},i\rangle ,\end{aligned}}}
where ⟨ m ⟩ {\displaystyle \langle m\rangle } is the average atomic mass, ω κ p , α {\displaystyle \omega _{\kappa _{p},\alpha }} and d κ p , α {\displaystyle \mathbf {d} _{{\boldsymbol {\kappa }}_{p},\alpha }} are the frequency and atomic displacement due to a phonon with polarization α {\displaystyle \alpha } and momentum κ p {\displaystyle {\boldsymbol {\kappa }}_{p}} , and | κ e , i ⟩ {\displaystyle |{\boldsymbol {\kappa }}_{e},i\rangle } is the electron wavefunction for an electron with momentum κ e {\displaystyle {\boldsymbol {\kappa }}_{e}} in band i. From Fermi's golden rule
γ ˙ e − p = 2 π ℏ | M e − p | 2 δ ( E e , i , κ e − E e , j , κ e + κ p ± ℏ ω κ p , α ) f e , i , κ e ( 1 − f e , i , κ e + κ p ) ( 1 2 ∓ 1 2 + f p , α , κ p ) {\displaystyle {\begin{aligned}{\dot {\gamma }}_{e-p}={\frac {2\pi }{\hbar }}|M_{e-p}|^{2}\delta (E_{e,i,\mathbf {\kappa } _{e}}-E_{e,j,\mathbf {\kappa } _{e}+\mathbf {\kappa } _{p}}\pm \hbar \omega _{\mathbf {\kappa } _{p},\alpha })f_{e,i,\mathbf {\kappa } _{e}}(1-f_{e,i,\mathbf {\kappa } _{e}+\mathbf {\kappa } _{p}})({\frac {1}{2}}\mp {\frac {1}{2}}+f_{p,\alpha ,\mathbf {\kappa } _{p}})\end{aligned}}}
where E e , i , κ e {\displaystyle E_{e,i,\mathbf {\kappa } _{e}}} is the energy of an electron in band i and momentum κ e {\displaystyle \mathbf {\kappa } _{e}} , f e , i , κ e {\displaystyle f_{e,i,\mathbf {\kappa } _{e}}} is the corresponding electron occupation, and f p , α , κ p {\displaystyle f_{p,\alpha ,\mathbf {\kappa } _{p}}} is the phonon occupancy.
The phonon-phonon coupling Competing with the generation of electrons is the downconversion of optical phonons into multiple acoustic phonons. The coupling arises from the crystal Hamiltonian (H) expanded in terms of the ionic displacement ( d i α {\displaystyle d_{i\alpha }} ) from the equilibrium position ( r i {\displaystyle r_{i}} ) of atom i in direction α {\displaystyle \alpha } in direction, i.e.,
H = H o + 1 2 ∑ i , j ∑ α , β ∂ 2 ⟨ φ ⟩ ∂ d i α ∂ d j β | o d i α d j β + 1 6 ∑ i , j , k ∑ α , β , γ ∂ 3 ⟨ φ ⟩ ∂ d i α ∂ d j β ∂ d k γ | o d i α d j β d k γ + ⋯ = H o + ∑ i , j ∑ α , β Γ i j , α β d i α d j β + ∑ i , j , k ∑ α , β , γ Ψ i j k , α β γ d i α d j β d k γ + ⋯ , {\displaystyle {\begin{aligned}\mathrm {H} &=\mathrm {H_{o}} +{\frac {1}{2}}\sum _{i,j}\sum _{\alpha ,\beta }{\frac {\partial ^{2}\langle \varphi \rangle }{\partial d_{i\alpha }\partial d_{j\beta }}}|_{\mathrm {o} }d_{i\alpha }d_{j\beta }+{\frac {1}{6}}\sum _{i,j,k}\sum _{\alpha ,\beta ,\gamma }{\frac {\partial ^{3}\langle \varphi \rangle }{\partial d_{i\alpha }\partial d_{j\beta }\partial d_{k\gamma }}}|_{\mathrm {o} }d_{i\alpha }d_{j\beta }d_{k\gamma }+\cdots \\&=\mathrm {H_{o}} +\sum _{i,j}\sum _{\alpha ,\beta }\Gamma _{ij,\alpha \beta }d_{i\alpha }d_{j\beta }+\sum _{i,j,k}\sum _{\alpha ,\beta ,\gamma }\Psi _{ijk,\alpha \beta \gamma }d_{i\alpha }d_{j\beta }d_{k\gamma }+\cdots ,\end{aligned}}}
where H o {\displaystyle \mathrm {H_{o}} } is the ground-state Hamiltonian, the linear term vanishes (as the ground state is found by minimizing the energy in terms of the ionic position), and Γ i j , α β {\displaystyle \Gamma _{ij,\alpha \beta }} and Ψ i j k , α β γ {\displaystyle \Psi _{ijk,\alpha \beta \gamma }} are the second- and third-order force constants between atoms i, j, and k when moved in along coordinate α {\displaystyle \alpha } , β {\displaystyle \beta } , and γ {\displaystyle \gamma } . The second order term is primarily responsible for the phonon dispersion, while the anharmonic (third order and higher) terms are responsible for thermal expansion as well as the phonon up- (multiple low-energy optical phonons combine to form a high-energy phonon) and downconversion (a high-energy phonon splits into multiple low-energy phonons). Typically, up- and down-conversion is dominated by the third-order interaction. Thus, the perturbation Hamiltonian used in Fermi's golden rule for phonon up- and downconversion is
H p − p ′ = Ψ κ p κ p ′ κ p ″ , α α ′ α ″ Ψ κ p κ p ′ κ p ″ , α α ′ α ″ = ∑ i j k ∑ α β γ Ψ i j k , α β γ s i α κ p s j α ′ κ p ′ s k α ″ κ p ″ ( ⟨ m ⟩ ω κ p , α ) 1 / 2 e x p [ i ( κ p ⋅ r i + κ p ′ ⋅ r j + κ p ⋅ r k ) ] , {\displaystyle {\begin{aligned}\mathrm {H} _{p-p}'&=\Psi _{{\boldsymbol {\kappa }}_{p}{\boldsymbol {\kappa }}_{p}'{\boldsymbol {\kappa }}_{p}'',\alpha \alpha '\alpha ''}\\\Psi _{{\boldsymbol {\kappa }}_{p}{\boldsymbol {\kappa }}_{p}'{\boldsymbol {\kappa }}_{p}'',\alpha \alpha '\alpha ''}&=\sum _{ijk}\sum _{\alpha \beta \gamma }\Psi _{ijk,\alpha \beta \gamma }{\frac {\mathbf {s} _{i\alpha {\boldsymbol {\kappa }}_{p}}\mathbf {s} _{j\alpha '{\boldsymbol {\kappa }}_{p}'}\mathbf {s} _{k\alpha ''{\boldsymbol {\kappa }}_{p}''}}{(\langle m\rangle \omega _{\kappa _{p},\alpha })^{1/2}}}\mathrm {exp} [i({\boldsymbol {\kappa }}_{p}\cdot r_{i}+{\boldsymbol {\kappa }}_{p}'\cdot r_{j}+{\boldsymbol {\kappa }}_{p}\cdot r_{k})],\end{aligned}}}
where s i α κ p {\displaystyle \mathbf {s} _{i\alpha \mathbf {\kappa _{p}} }} is the direction of displacement for atom i due to the phonon. The resulting downconversion rate, from Fermi's golden rule, is
γ ˙ p − p = ℏ π 16 | Ψ κ p κ p ′ κ p ″ , α α ′ α ″ | 2 × δ κ p , κ p ′ + κ p ″ δ ( ω κ p , α − ω κ p ′ , α ′ − ω κ ″ , α ″ ) × ( f p ′ + f p ″ + 1 ) , {\displaystyle {\dot {\gamma }}_{p-p}={\frac {\hbar \pi }{16}}|\Psi _{{\boldsymbol {\kappa }}_{p}{\boldsymbol {\kappa }}_{p}'{\boldsymbol {\kappa }}_{p}'',\alpha \alpha '\alpha ''}|^{2}\times \delta _{\kappa _{p},\kappa _{p}'+\kappa _{p}''}\delta (\omega _{\kappa _{p},\alpha }-\omega _{\kappa _{p}',\alpha '}-\omega _{\kappa '',\alpha ''})\times