Aaron M Burger1, Radhe Agarwal2, Alexey Aprelev3, Edward Schruba1, Alejandro Gutierrez-Perez2, Vladimir M Fridkin2,4, Jonathan E Spanier1,2,3. 1. Department of Electrical and Computer Engineering, Drexel University, Philadelphia, PA 19104, USA. 2. Department of Materials Science and Engineering, Drexel University, Philadelphia, PA 19104, USA. 3. Department of Physics, Drexel University, Philadelphia, PA 19104, USA. 4. Shubnikov Institute of Crystallography, Russian Academy of Sciences, Leninsky Prospect 59, Moscow, 117333, Russian Federation.
Abstract
The quantum phenomenon of shift photovoltaic current was predicted decades ago, but this effect was never observed directly because shift and ballistic currents coexist. The atomic-scale relaxation time of shift, along with the absence of a photo-Hall behavior, has made decisive measurement of shift elusive. Here, we report a facile, direct-current, steady-state method for unambiguous determination of shift by means of the simultaneous measurements of linear and circular bulk photovoltaic currents under magnetic field, in a sillenite piezoelectric crystal. Comparison with theoretical predictions permits estimation of the signature length scale for shift. Remarkably, shift and ballistic photovoltaic currents under monochromatic illumination simultaneously flow in opposite directions. Disentangling the shift and ballistic contributions opens the way for quantitative, fundamental insight into and practical understanding of these radically different photovoltaic current mechanisms and their relationship.
The quantum phenomenon of shift photovoltaic current was predicted decades ago, but this effect was never observed directly because shift and ballistic currents coexist. The atomic-scale relaxation time of shift, along with the absence of a photo-Hall behavior, has made decisive measurement of shift elusive. Here, we report a facile, direct-current, steady-state method for unambiguous determination of shift by means of the simultaneous measurements of linear and circular bulk photovoltaic currents under magnetic field, in a sillenite piezoelectric crystal. Comparison with theoretical predictions permits estimation of the signature length scale for shift. Remarkably, shift and ballistic photovoltaic currents under monochromatic illumination simultaneously flow in opposite directions. Disentangling the shift and ballistic contributions opens the way for quantitative, fundamental insight into and practical understanding of these radically different photovoltaic current mechanisms and their relationship.
The bulk photovoltaic current in crystals among 20 symmetry point groups that lack a center of inversion consists of two parts: ballistic and shift photovoltaic currents. Shift photovoltaic current (), which contributes to the bulk photovoltaic effect (BPE) (–), has recently generated substantial interest in solid state science, physical chemistry, optoelectronics, and engineering, in part because this remarkable phenomenon involving inversion symmetry breaking has been proposed as a mechanism for photogenerated current, and because of its possible connections with electronic structure in topological insulators () and Weyl semimetals (), photovoltaics in the absence of a junction (–), and nonlinear light-matter interactions such as second-harmonic generation (). The BPE is a widely explored phenomenon in noncentrosymmetric materials; violation of the principle of detailed balance in the BPE and its application in nanoscale electrode geometry have yielded unity or higher quantum efficiency, and overcoming of the Shockley-Queisser limit on power conversion efficiency (), not attainable using traditional p-n junction photovoltaics. Although the quantum phenomenon of shift photovoltaic current was theoretically predicted decades ago (), the atomic-scale relaxation time of shift current (≈0.1 to 1 fs) has made it experimentally inaccessible.Decades following the experimental discovery of the BPE (, ) and subsequent development of theoretical models of ballistic and shift currents (, , ), expression of the shift mechanism of the bulk photovoltaic current in the form of a first principles density functional theoretical framework (, ) advancing the seminal theoretical description () has reinvigorated this research, opening the way for materials design and prediction of shift current excitations (, ). In a number of recent experimental works (–), it has been reported that the shift mechanism may be responsible for any photovoltaic current in a crystal lacking inversion symmetry and may account for enhanced BPE (, ) and even anomalously large carrier collection length arising from local excitation (, , ).Uniform illumination of a homogeneous asymmetric optically active crystal permits observation of the bulk photovoltaic current , where the first part is the linear bulk photovoltaic current, obtained with linearly polarized light, the second part is the circular photovoltaic current, obtained using circularly polarized light, I is the light intensity, e and are the components of the incident light polarization vector, i[ee*] determines the degree of circular polarization σ, and and are third and second rank tensors; , where q is the photon wave vector. This phenomenon was termed the BPE, and the first experimental observations of the linear BPE were in ferroelectrics SbSI (point group C2) () and LiNbO3 (C3) (), and of the circular BPE, in Te (D3) (, ). A comprehensive review of the theories of the linear and circular BPEs in a broad group of ferroelectrics and piezoelectrics has been presented ().The microscopic nature of the BPE is related to two mechanisms: ballistic and shift. The ballistic mechanism is associated with the excitation of nonthermalized (hot) carriers in a noncentrosymmetric crystal, leading to an asymmetric distribution of hot carrier momenta in the band and a violation of the Boltzmann principle of detailed balance (). In the ballistic mechanism, the photoexcited hot carriers lose their energy and relax to the bottom of the band over free path, that is, thermalization length (). The relaxation time of the ballistic current τb ≈ 10− 12 − 10− 14 s , where Ω(= eB/m*c) is the Larmor frequency and B is the classic magnetic field. This permits observation of the Hall effect under illumination for the ballistic linear and circular bulk photovoltaic currents and of the mobility of nonthermalized carriers (, , , ).Shift current is a separate, quantum mechanism whereby accounting for the nondiagonal elements of the density matrix, a consistent kinetic theory was developed () and later within the framework of density functional theory (, ). Shift itself does not have a transport mechanism, it does not involve inertia, and its short (atomic) relaxation time τs prevents observation of current kinetics or a photo-Hall effect, because for any classical magnetic field B, . These characteristics explain why shift photovoltaic current has never been observed experimentally despite a number of reports attributing measured photovoltaic currents to shift (–), including transient current measured indirectly from terahertz radiation (, ). Also, it was shown that impurity band transitions can contribute to both shift and ballistic currents and are strongly dependent on the nature of the impurity ().
RESULTS
The ballistic jb and shift js currents are expressed explicitly as jb = eαI(ℏω)−1l0ξϕ and js = eαI(ℏω)−1R (, ). Here, ξ is the photoexcitation asymmetry parameter, which is a function of wave vector but is represented here as a scalar. l0 is the free path of nonthermalized carriers, I is the light intensity, ℏω is the photon energy, α is the absorption coefficient, and ϕ is the quantum yield. The quantity R, called the shift vector, has the meaning of carrier displacement in real space per absorbed photon. Under a very rough approximation |R| ≈ aξ, where a is the lattice parameter (for cubic sillenites a ≈ 1 nm).The linear BPE current is the sum of the ballistic and shift contributions. It has been shown that the circular BPE current is purely ballistic (, , ). To experimentally separate the ballistic and shift currents, we consider here the simultaneous linear and circular currents in a cubic sillenite under magnetic field based on the following relations. The linear bulk photovoltaic current, for example, along [100], isThe mobility of nonthermalized electronsmay be obtained from the photo-Hall effect for current collected under circularly polarized light (jc), where is the Hall component under circularly polarized light, and B is the magnetic field applied along [010].From measurement of the Hall component of linear current , we obtain the expressions for ballistic jb and shift js currentsThese equations are based on two main theoretical properties (, ). First, shift current does not contribute to the Hall signal, and second, the current under circularly polarized light is purely ballistic.As shown in (, ), the magnetic field changes the asymmetry parameter in the Hall direction and could create additional current, which depends linearly on B. Under this condition, the expression for shift current js is thenwhere x (≳1) is a model-dependent parameter characterizing magneto-induced asymmetry [see p. 92 in ()].We investigated photovoltaic currents in single-crystal, optically active cubic piezoelectric Bi12GeO20 (point group T) having one linear component (in the compact Voigt notation) and one circular component . A circular and linear BPE in sillenites has been observed previously (). To observe the circular effect, we use an electro-optic modulation technique (, ), measuring the linear and circular photovoltaic currents and its Hall component under a magnetic field B ≲ 1 T. The energy gap of Bi12GeO20 is ≈3.2 eV (), and we performed the measurements in the extrinsic (impurity) region (450 to 650 nm). It was shown that sillenites have n-type photoconductivity (, ) arising from Bi+3 donors that produce a broad absorption shoulder (≈2.3 to 3.2 eV) (, ), and the mobility of thermalized electrons is extremely small (μth ≈ 10− 2 − 10− 6 cm2 V− 1 s− 1) and depends on the filled trap concentration in the extrinsic region.Measurements were performed at room temperature (≈300 K) under polarization-modulated monochromatic laser illumination intensities ranging from ≈300 to 2100 mW/cm2 using commercial single crystals. Photovoltaic currents jl and jc were collected through transparent electrodes on the (001) and () faces, and other electrodes were used to measure their corresponding Hall components and (Fig. 1, A and B, and Materials and Methods).
Fig. 1
Experimental scheme for probing bulk photovoltaic current in the extrinsic spectral.
(A) Schematic illustration of the sample geometry and electrode and field configuration. (B) Block diagram illustrating experimental configuration. (1) Bi12GeO20 crystal, (2) lock-in amplifier, (3) modulator driving electronics, (4) function generator, (5) Kr/Ar ion laser, (6) polarizer, (7) electromagnet, (8) electro-optic modulator, (9) λ/2 or λ/4 achromatic wave plate, and (10) Hall contacts. (C) Measured optical absorption in Bi12GeO20.
Experimental scheme for probing bulk photovoltaic current in the extrinsic spectral.
(A) Schematic illustration of the sample geometry and electrode and field configuration. (B) Block diagram illustrating experimental configuration. (1) Bi12GeO20 crystal, (2) lock-in amplifier, (3) modulator driving electronics, (4) function generator, (5) Kr/Ar ion laser, (6) polarizer, (7) electromagnet, (8) electro-optic modulator, (9) λ/2 or λ/4 achromatic wave plate, and (10) Hall contacts. (C) Measured optical absorption in Bi12GeO20.The measured absorption spectrum (Fig. 1C) quantifies the extrinsic character in the range studied, where a broad absorption shoulder has been attributed to an antisite defect (BiGe) (). The angular dependence of the linear photovoltaic current in [100] (β is the angle between the light polarization and [010]) is in agreement with tensorial properties of crystal (Fig. 2A). The obtained spectral dependences of jl and jc were collected for a series of wavelengths across the spectral region of λ ≈ 450 to 650 nm (Fig. 2B). Shown in Fig. 2C is the Hall component of the circular photovoltaic current as a function of magnetic field for right (σ+) and left (σ−) circularly polarized light for λ = 488.1 nm. The measurements of jc and its Hall components indicate that the nonthermalized mobility μnth approaches 600 cm2 V–1 s− 1, exceeding by at least 104 the thermalized mobility (, ). The high value for nonthermalized carriers has been obtained from photovoltaic Hall measurements for many piezoelectrics and ferroelectrics (, , ). The spectral dependence of mobility may be caused by thermalization, related to electron-defect center interaction. The ballistic current jb is proportional to the product of μnth and n, the concentration of nonthermalized electrons. We suppose that the minimum in μnth in the region of the impurity band is due to the stronger interaction of carriers in this spectral region and that n is larger in this spectral range because of the density of impurity centers in this range. These results confirm the extrinsic origin of photovoltaic current by excitation of donor levels and, significantly, that thermalization of nonequilibrium carriers is caused by scattering by ionized impurities. This can explain why the nonthermalized mobility is large outside the impurity band and low within it.
Fig. 2
Bulk photovoltaic currents under linearly and circularly polarized illumination in Bi12GeO20, and the corresponding Hall component.
(A) Measured variation in photocurrent jl, in Bi12GeO20, under linearly polarized λ = 488.1 nm light, as a function of β, the angle between the light polarization and [010]. (B) Measured spectral dependences of linear jl and circular jc photovoltaic currents, normalized to incident power. (C) Hall component of the circular photovoltaic current as a function of magnetic field, for right (σ+) and left (σ−) circularly polarized λ = 488.1 nm light.
Bulk photovoltaic currents under linearly and circularly polarized illumination in Bi12GeO20, and the corresponding Hall component.
(A) Measured variation in photocurrent jl, in Bi12GeO20, under linearly polarized λ = 488.1 nm light, as a function of β, the angle between the light polarization and [010]. (B) Measured spectral dependences of linear jl and circular jc photovoltaic currents, normalized to incident power. (C) Hall component of the circular photovoltaic current as a function of magnetic field, for right (σ+) and left (σ−) circularly polarized λ = 488.1 nm light.The spectral dependences of jb and js, each normalized to incident power, are presented in Fig. 3A, and corresponding spectral dependences of the tensor (Glass) coefficient for jb and js are presented in Fig. 3B. Our data and analysis reveal that shift exists in the extrinsic region and that the ballistic and shift currents are of the same order of magnitude, as predicted by theory (, , , ), but flow in opposite directions; the latter is a result not anticipated from the theory. The spectral distribution of μnth has a sharp minimum that is close to the maximum of js and jb (Fig. 3C). It shows, in accordance with (), that in this spectral interval, the impurity band transitions are responsible for ballistic and shift currents, permitting a rough estimate of R. On the basis of the data in Fig. 3A, for ℏω = 2.5 eV, R = 0.4 × 10− 2 nm, and correspondingly, ξ ≈ R/a ≈ 0.4 × 10− 2. This value is in good agreement with the asymmetry parameter obtained for the Coulomb impurity model (). For ℏω = 2.5 eV, ϕl0 ≈ 0.2 nm, and for ϕ < 1, l0 can reach a few nanometers.
Fig. 3
Separation of shift and ballistic contributions to the bulk photovoltaic current.
(A) Spectral dependences of shift js and ballistic jb photovoltaic currents, each normalized to incident power. (B) Corresponding spectral shift and ballistic Glass coefficients. (C) Spectral dependence of the nonthermalized carrier mobility in Bi12GeO20.
Separation of shift and ballistic contributions to the bulk photovoltaic current.
(A) Spectral dependences of shift js and ballistic jb photovoltaic currents, each normalized to incident power. (B) Corresponding spectral shift and ballistic Glass coefficients. (C) Spectral dependence of the nonthermalized carrier mobility in Bi12GeO20.
DISCUSSION
We note that under small classic magnetic field B, no detectable change was observed in jl or jc with B, and therefore, we conclude that the effect of magneto-induced asymmetry (Eq. 5) is small. The change of due to the Faraday effect can be neglected in comparison with the photo-Hall effect. In principle, the mobilities for linear and circular currents are different even in the absence of the magneto-photovoltaic effect. However, following the theory (, , , ), this difference is within a factor of two, and accordingly does not influence the value of js by more than an order of magnitude. With regard to the validity of Eqs. 2 to 5, it is well established that the photovoltaic current under circularly polarized light in the absence of magnetic field is purely ballistic (injection). While the mechanism for photovoltaic current for circularly polarized light in the presence of a magnetic field may be altered because of the change in the balance between the energy renormalization current and shift current, we have assumed here that in a magnetic field, the photovoltaic current in circularly polarized light remains ballistic and that the Hall mobility under circularly polarized light is ballistic. Last, we note that shift and ballistic currents in the extrinsic regime may be distinguished in other materials. Application of the same method to the sillenite Bi12SiO20 also permits their separation, observation of a different sign, but different spectral dependence (Supplementary Materials).Decades following its discovery and despite a well-developed theoretical framework for shift and ballistic currents, quantifying the extent of their contributions from these coexisting quantum and quasi-classical mechanisms has remained experimentally inaccessible. Analysis of the steady-state linear and circular photovoltaic currents for the separation of shift and ballistic contributions to the photovoltaic current in perovskites and other crystals that belong to 18 optically active point groups, including those of extrinsic origin, is now feasible. We anticipate that facile experimental disambiguation of shift should encourage new exploration of connections of first principles–based calculations with experiment. The experimental prescription holds great promise for exploiting these remarkable light-matter interactions and may be helpful in efforts to obtain higher power conversion efficiency in devices based on inversion symmetry breaking of the BPE (, ).
MATERIALS AND METHODS
Measurements were performed using commercially sourced (MetaLaser, China) Bi12GeO20 single crystals with dimensions of 5 by 5 by 5 mm3 along the x, y, and z axes, which coincide with < 100 >. Optical absorption of the identical Bi12GeO20 crystal was collected using ultraviolet-visible spectrophotometry (Shimadzu 1450). Transparent conducting electrodes (indium tin oxide) were deposited by radio frequency sputtering onto (001) and () faces to collect jl and jc, while other electrodes were used to measure their corresponding Hall components and . For electro-optic polarization modulation–assisted measurements of linearly and circularly polarized photovoltaic currents, the crystal was illuminated with selected lines from a Kr-Ar ion mixed-gas laser (Coherent Innova I90-6) through an electro-optic modulator permitting generation of positive (σ+) and negative (σ−) circular light polarization and, correspondingly, the changed sign of circular current jc. A linearly polarized beam from the laser is polarization modulated at 41 Hz, resulting in oscillation between two mutually perpendicular states before passing through a wave plate, λ/2 or λ/4 (Thorlabs), providing for measurement of jl and jc, respectively. Polarization modulation was achieved using a Pockels cell electro-optic modulator (ConOptics, KD*P 350-50) driven by the amplified output of a function generator using voltage magnitude selected for wavelength-dependent half-wave voltage, where a typical half wave voltage for this model at, for example, 500 nm, is 455 V. The resultant polarization states were studied using an analyzer and phototransistor prior to current measurements in the crystals. Current was recorded using a Stanford Research SR-830 digital lock-in amplifier. This phase-sensitive detector was referenced to the synchronous TTL output from the function generator. Photovoltaic current (Hall current) measurements were recorded with z-direction (x direction) electrodes under short circuit. Typical settings for the lock-in amplifier include low-pass filtering with 18 dB per decade roll-off, a 3-s integration time constant, and the input impedance from the lock-in amplifier is 1 kilohm (for current measurements). Detailed descriptions of similar modulation current measurement methods can be found elsewhere, for example, in () and ().