Phonovoltaic

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.^[1] The pV cell is an application of heat transfer physics^[2] 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.^[1]^[3] 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.^[1] For a PV cell, many materials are available with a band gap ( $\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,\mathrm {O} }$ ). Thus, novel materials are required with both energetic optical phonon modes ( $E_{p,\mathrm {O} }\geq 100$ meV, e.g., graphene, diamond, or boron nitride) and a small band gap ( $\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.^[1]^[4] 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,^[5]^[6] 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 ( $\eta _{\phi }$ ) is quantified by the ratio of the band gap and optical phonon energy, that is^[1]

{\begin{aligned}\eta _{\phi }=\Delta E_{e,g}/E_{p,\mathrm {O} }.\end{aligned}}

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 ( $\eta _{QE}$ ) is quantified by the tendency of a hot optical phonon to downconvert rather than generate an electron-hole pair, that is^[1]

{\begin{aligned}\eta _{QE}={\dot {\gamma }}_{e-p}^{*}={\frac {{\dot {\gamma }}_{e-p}}{{\dot {\gamma }}_{e-p}+{\dot {\gamma }}_{p-p}}},\end{aligned}}

where ${\dot {\gamma }}_{e-p}$ is the rate of generation and ${\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 ( $\eta _{pn}$ ) is limited by the Carnot efficiency given by^[1]

{\begin{aligned}\eta _{\mathrm {C} }=1-{\frac {T_{pV}}{T_{p,\mathrm {O} }}},\end{aligned}}

where $T_{\mathrm {pV} }$ is the temperature of the pV cell and $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_{\mathrm {B} }T$ , where $k_{\mathrm {B} }$ is the Boltzmann constant and $T$ is the temperature). Indeed, the p-n junction efficiency is approximately^[1]

{\begin{aligned}\eta _{pn}\simeq \eta _{\mathrm {C} }[1-0.75\exp(-0.1\Delta E_{e,g}/k_{\mathrm {B} }T_{\mathrm {pV} })].\end{aligned}}

Thus, the overall efficiency ( $\eta _{\mathrm {pV} }$ ) is^[1]

{\begin{aligned}\eta _{\mathrm {pV} }=\eta _{\mathrm {C} }{\dot {\gamma }}_{e-p}^{*}{\frac {\Delta E_{e,g}}{E_{p,\mathrm {O} }}}[1-0.75\exp(-0.1\Delta E_{e,g}/k_{\mathrm {B} }T_{pV})],\end{aligned}}

where the temperature independent terms become the material figure of merit ( $Z_{\mathrm {pV} }$ ),^[1]

{\begin{aligned}Z_{\mathrm {pV} }={\dot {\gamma }}_{e-p}^{*}{\frac {\Delta E_{e,g}}{E_{p,\mathrm {O} }}}.\end{aligned}}

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_{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, $\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 ^[7]

${\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 $\langle m\rangle$ is the average atomic mass, $\omega _{\kappa _{p},\alpha }$ and $\mathbf {d} _{{\boldsymbol {\kappa }}_{p},\alpha }$ are the frequency and atomic displacement due to a phonon with polarization $\alpha$ and momentum ${\boldsymbol {\kappa }}_{p}$ , and $|{\boldsymbol {\kappa }}_{e},i\rangle$ is the electron wavefunction for an electron with momentum ${\boldsymbol {\kappa }}_{e}$ in band i. From Fermi's golden rule

${\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,\mathbf {\kappa } _{e}}$ is the energy of an electron in band i and momentum $\mathbf {\kappa } _{e}$ , $f_{e,i,\mathbf {\kappa } _{e}}$ is the corresponding electron occupation, and $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\alpha }$ ) from the equilibrium position ( $r_{i}$ ) of atom i in direction $\alpha$ in direction, i.e.,^[2]^[8]

${\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 $\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 $\Gamma _{ij,\alpha \beta }$ and $\Psi _{ijk,\alpha \beta \gamma }$ are the second- and third-order force constants between atoms i, j, and k when moved in along coordinate $\alpha$ , $\beta$ , and $\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 ^[9]

${\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 $\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

${\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 (f_{p}'+f_{p}''+1),$

where two phonons are produced with polarization $\alpha '$ and $\alpha ''$ and momentum $\mathbf {\kappa } _{p}'$ and $\mathbf {\kappa } _{p}''$ .

The suitability of graphene as a phonovoltaic material

As outlined above, an efficient pV cell requires a material with an optical phonon mode more energetic than the bandgap, which in turn is much more energetic than the thermal energy at the intended operating temperature $(E_{p,\mathrm {O} }\simeq \Delta E_{e,g}\gg k_{\mathrm {B} }T_{\mathrm {pV} })$ . Furthermore, the pV cell requires a material wherein a hot optical phonon prefers to produce an electron rather than multiple low energy acoustic phonons ( ${\dot {\gamma }}_{e-p}^{*}\rightarrow 1$ ).

Very few materials offer this combination of properties. Indeed, the vast majority of crystals have optical phonon energies limited to below 50 meV, and those with more energetic optical phonons tend to have much more energetic band gaps. In general, a material with a first-row element (periodic table) is required to have a highly energetic optical phonon. However, the high electronegativity of a first-row elements tends to create a very large band gap, as in diamond and the boron nitride allotropes.^[10] Graphene is one of the few materials which diverges from this trend, with no bandgap and an exceptionally energetic optical phonon mode near 200 meV.^[11]^[12] Thus, graphene has been the initial target for development of a phonovoltaic material through the opening and tuning of its bandgap.^[13]^[3]

Opening and tuning the bandgap of graphene has received substantial attention, and numerous strategies have been suggested and investigated. These include the use of uniaxial strain,^[14] electric fields,^[15] and chemical doping and functionalization.^[16] In general, these mechanisms work by either changing the symmetry of graphene (both Carbon atoms in the unit cell are identical) or hybridization ( $sp^{2}$ ).

In the first phonovoltaic material investigations, it has been suggested that the latter technique destroys the electron-phonon coupling while the former preserves it.^[13]^[3] In particular, these investigations predict that hydrogenating graphene, to produce graphane, reduces the electron-phonon coupling so substantially that the material figure of merit vanishes;^[13] and that doping graphene with boron nitride maintains the strong electron-phonon coupling in graphene, such that its figure of merit is predicted to reach 0.65 and enable heat harvesting with twice the efficiency of a typical thermoelectric generator.^[3]

References

^ ^a ^b ^c ^d ^e ^f ^g ^h ⁱ ^j ^k ^l ^m Melnick, C.; Kaviany, M. (March 2016). "Phonovoltaic I. Harvesting hot optical phonons in a nanoscale p-n junction". Physical Review B. 93 (9): 094302. Bibcode:2016PhRvB..93i4302M. doi:10.1103/PhysRevB.93.094302.
^ ^a ^b Kaviany, M. (2014). Heat transfer physics (2nd ed.). New York: Cambridge University Press. ISBN 978-1-107041783.
^ ^a ^b ^c ^d Melnick, C.; Kaviany, M. (December 2016). "Phonovoltaic. III. Electron-phonon coupling and figure of merit of graphene:BN". Physical Review B. 94 (24): 245412. Bibcode:2016PhRvB..94x5412M. doi:10.1103/PhysRevB.94.245412.
^ da Silva, L. (2004). "Micro-thermoelectric cooler: interfacial effects on thermal and electrical transport". Int. J. Heat Mass Trans. 47 (10–11): 2417. doi:10.1016/j.ijheatmasstransfer.2003.11.024.
^ Span, G. (2007). "Miniaturized TEG with thermal generation of free carriers". Physica Status Solidi RRL. 1 (6): 241. Bibcode:2007PSSRR...1..241S. doi:10.1002/pssr.200701171. S2CID 44740940.
^ Wagner, M. (2007). "Thermoelectric power generation using large-area Si/SiGe pn-junctions with carying Ge content". Semiconductor Science and Technology. 22 (1): S173. Bibcode:2007SeScT..22S.173W. doi:10.1088/0268-1242/22/1/S41. S2CID 27885891.
^ Baroni, S.; de Gironcoli, S.; Corso, A. D.; Giannozzi, P. (2001). "Phonons and related crystal properties from density-functional perturbation theory". Rev. Mod. Phys. 73 (2): 515–562. arXiv:cond-mat/0012092. Bibcode:2001RvMP...73..515B. doi:10.1103/RevModPhys.73.515. S2CID 19021637.
^ Srivistava, G. (1990). The Physics of Phonons. Bristol: Adam Hilger.
^ Hellman, O.; Abrikosov, I. A. (2013). "Temperature-dependent effective third-order interatomic force constants from first principles". Phys. Rev. B. 88 (14): 144301. arXiv:1308.5436. Bibcode:2013PhRvB..88n4301H. doi:10.1103/PhysRevB.88.144301. S2CID 39340446.
^ Madelung, O. (2004). Semiconductors: Data Handbook. Berlin: Springer.
^ Adachi, S. (2005). Properties of Group-IV, III-V, and II-IV Semiconductors. Sussex: Wiley.
^ Maultzsch, J.; et, al. (2004). "Phonon Dispersion in Graphite". Physical Review Letters. 92 (7): 075501. Bibcode:2004PhRvL..92g5501M. doi:10.1103/PhysRevLett.92.075501. PMID 14995866.
^ ^a ^b ^c Melnick, C.; Kaviany, M. (2016). "Phonovoltaic. II. Tuning band gap to optical phonon in graphite". Physical Review B. 93 (12): 125203. Bibcode:2016PhRvB..93l5203M. doi:10.1103/PhysRevB.93.125203.
^ Periera, V.M. (2009). "Tight-binding approach to uniaxial strain in graphene". Physical Review B. 80 (4): 045401. arXiv:0811.4396. Bibcode:2009PhRvB..80d5401P. doi:10.1103/PhysRevB.80.045401. S2CID 119255776.
^ Park, C. (2015). "Electronic Properties of Bilayer Graphene Strongly Coupled to Interlayer Stacking and an External Electric Field". Physical Review Letters. 115 (1): 015502. arXiv:1505.01959. Bibcode:2015PhRvL.115a5502P. doi:10.1103/PhysRevLett.115.015502. PMID 26182105. S2CID 9975742.
^ Georgakilas, V. (2012). "Functionalization of Graphene: Covalent and Non-Covalent Approaches, Derivatives and Applications". Chemical Reviews. 112 (11): 6156–6214. doi:10.1021/cr3000412. PMID 23009634.

[pv1-1] ^ ^a ^b ^c ^d ^e ^f ^g ^h ⁱ ^j ^k ^l ^m Melnick, C.; Kaviany, M. (March 2016). "Phonovoltaic I. Harvesting hot optical phonons in a nanoscale p-n junction". Physical Review B. 93 (9): 094302. Bibcode:2016PhRvB..93i4302M. doi:10.1103/PhysRevB.93.094302.

[HTPbook-2] Kaviany, M. (2014). Heat transfer physics (2nd ed.). New York: Cambridge University Press. ISBN 978-1-107041783.

[pv3-3] Melnick, C.; Kaviany, M. (December 2016). "Phonovoltaic. III. Electron-phonon coupling and figure of merit of graphene:BN". Physical Review B. 94 (24): 245412. Bibcode:2016PhRvB..94x5412M. doi:10.1103/PhysRevB.94.245412.

[4] Silva, L. (2004). "Micro-thermoelectric cooler: interfacial effects on thermal and electrical transport". Int. J. Heat Mass Trans. 47 (10–11): 2417. doi:10.1016/j.ijheatmasstransfer.2003.11.024.

[5] Span, G. (2007). "Miniaturized TEG with thermal generation of free carriers". Physica Status Solidi RRL. 1 (6): 241. Bibcode:2007PSSRR...1..241S. doi:10.1002/pssr.200701171. S2CID 44740940.

[6] Wagner, M. (2007). "Thermoelectric power generation using large-area Si/SiGe pn-junctions with carying Ge content". Semiconductor Science and Technology. 22 (1): S173. Bibcode:2007SeScT..22S.173W. doi:10.1088/0268-1242/22/1/S41. S2CID 27885891.

[7] Baroni, S.; de Gironcoli, S.; Corso, A. D.; Giannozzi, P. (2001). "Phonons and related crystal properties from density-functional perturbation theory". Rev. Mod. Phys. 73 (2): 515–562. arXiv:cond-mat/0012092. Bibcode:2001RvMP...73..515B. doi:10.1103/RevModPhys.73.515. S2CID 19021637.

[8] Srivistava, G. (1990). The Physics of Phonons. Bristol: Adam Hilger.

[9] Hellman, O.; Abrikosov, I. A. (2013). "Temperature-dependent effective third-order interatomic force constants from first principles". Phys. Rev. B. 88 (14): 144301. arXiv:1308.5436. Bibcode:2013PhRvB..88n4301H. doi:10.1103/PhysRevB.88.144301. S2CID 39340446.

[madelung-10] Madelung, O. (2004). Semiconductors: Data Handbook. Berlin: Springer.

[adachi-11] Adachi, S. (2005). Properties of Group-IV, III-V, and II-IV Semiconductors. Sussex: Wiley.

[Maultzsch-12] Maultzsch, J.; et, al. (2004). "Phonon Dispersion in Graphite". Physical Review Letters. 92 (7): 075501. Bibcode:2004PhRvL..92g5501M. doi:10.1103/PhysRevLett.92.075501. PMID 14995866.

[pv2-13] Melnick, C.; Kaviany, M. (2016). "Phonovoltaic. II. Tuning band gap to optical phonon in graphite". Physical Review B. 93 (12): 125203. Bibcode:2016PhRvB..93l5203M. doi:10.1103/PhysRevB.93.125203.

[pereira-14] Periera, V.M. (2009). "Tight-binding approach to uniaxial strain in graphene". Physical Review B. 80 (4): 045401. arXiv:0811.4396. Bibcode:2009PhRvB..80d5401P. doi:10.1103/PhysRevB.80.045401. S2CID 119255776.

[park-15] Park, C. (2015). "Electronic Properties of Bilayer Graphene Strongly Coupled to Interlayer Stacking and an External Electric Field". Physical Review Letters. 115 (1): 015502. arXiv:1505.01959. Bibcode:2015PhRvL.115a5502P. doi:10.1103/PhysRevLett.115.015502. PMID 26182105. S2CID 9975742.

[georgakilas-16] Georgakilas, V. (2012). "Functionalization of Graphene: Covalent and Non-Covalent Approaches, Derivatives and Applications". Chemical Reviews. 112 (11): 6156–6214. doi:10.1021/cr3000412. PMID 23009634.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]