CarrierCapture.jl

A set of codes to compute carrier capture and recombination rates in semiconducting compounds. This topic has a rich history starting from the work by Huang and Rhys. Our implementation was inspired by the approach (and FORTRAN code) employed by Alkauskas and coworkers, but has been adapted to also describe anharmonic potential energy surfaces.

Please cite: Kim, S.; Hood, S. N.; Gerwen, P. van; Whalley, L. D.; Walsh, A. CarrierCapture.Jl: Anharmonic Carrier Capture. Journal of Open Source Software 2020, 5 (47), 2102. https://doi.org/10.21105/joss.02102

Installation

The codes are written in Julia, while the scripts and Jupyter Notebooks also contain Python and use pymatgen and pawpyseed (This has been tested on Scientific Linux 7 and Linux Mint 18), which are assumed to be installed. The Brooglie package is used to solve the time-independent Schrödinger equation.

Install the package by:

julia> using Pkg

julia> Pkg.add(PackageSpec(url="https://github.com/WMD-group/CarrierCapture.jl.git"))

To run the unit tests for the package, use the Pkg.test function.

julia> Pkg.test("CarrierCapture")

Development

The project is hosted on Github. Please use the issue tracker for feature requests, bug reports and more general questions. If you would like to contribute, please do so via a pull request.

Usage

A typical workflow will consist of several steps, implemented in a series of short programs.

1. Prepare a sequence of atomic structure models with displacements that interpolate between two defect configurations (e.g. a site vacancy in charge states q=0 and q=+1). Run single-point energy calculations on these structures, and extract the total energies. Scripts for preprocessing may be found in script.

2. Find a best fit for the energy calculations of the deformed structures (potential) to generate potential energy surfaces (PES). Solve the 1D Schrödinger equation for each PES to obtain their phonon (nuclear) wavefunctions.

3. Construct configuration coordinate (conf_coord) to calculate the wavefunction overlap between each PES, which forms part of the temperature-dependent capture coefficient.

The command-line interface (GetPotential.jl and GetRate.jl) is depreciated. Use Jupyter Notebook examples as a template. A detailed description provided in the Usage section.

Examples

The following examples are provided to illustrate some of the applications of these codes. The input data has been generated from density functional theory (DFT) using VASP, but the framework can easily be adapted to accept output from other electronic structure calculators.

• Sn_Zn in CZTS: Harmonic approximation

• DX-center in GaAs: Anharmonic fitting

• Electron-phonon coupling: Electron-phonon coupling matrix element

Theory

The electronic matrix element frequently causes feelings of discomfort (Stoneham, 1981)

The capture of electrons or holes by point defects in a crystalline materials requires the consideration of a number of factors including the coupling between electronic and vibrational degrees of freeedom. Many theories and approximations have been developed to describe the reaction kinetics.

The capture coefficient between an initial and final state for this computational set up is given by (eq. 22 in Alkauskas and coworkers):

\[C_p = V \frac{2\pi}{\hbar} g W_{if}^2 \sum_m w_m \sum_n |\langle \xi_{im}| Q - Q_0 | \xi_{fn}\rangle|^2 \delta (\Delta E + m\hbar\Omega_i -n\hbar\Omega_f )\]

Here, $V$ is the volume of the supercell, $W_{if}$ is the electron-phonon overlap and $ξ_{im}$ and $ξ_{fn}$ describe the wavefunctions of the $m^{th}$ and $n^{th}$ phonons in the initial $i$ and final $f$ states. The final delta-function term serves to conserve energy and in practice is replaced by a smearing Gaussian of finite width $σ$.

User Warning

The values produced by this type of analysis procedure are sensitive to the quality of the input. We expect that most input data will have been generated by DFT where the basis set, k-points, and ionic forces have been carefully converged. In addition, the alignment of energy surfaces for defects in different charge states requires appropriate finite-size corrections (e.g. see Freysoldt and coworkers).

Extended Reading List

Theory Development

• Heny and Lang, Nonradiative capture and recombination by multiphonon emission in GaAs and GaP (1977)

Seminal contribution that introduces many important concepts

• Huang, Adiabatic approximation theory and static coupling theory of nonradiative transition (1981)

Context for the static approximation that we employ

• Stoneham, Non-radiative transitions in semiconductors (1981)

Review on theory and various models of recombination

• Markvart, Determination of potential surfaces from multiphonon transition rates (1981)

Discussion and treatment of anharmonicity

• Markvart, Semiclassical theory of non-radiative transitions (1981)

Semiclassical treatment of matrix elements following Landau and Holstein

Applications

• Kim et al, Anharmonic Lattice Relaxation during Non-radiative Carrier Capture (2019)

• Kim et al, Lone-pair effect on carrier capture in Cu2ZnSnS4 solar cells (2019)

• Kim et al, Identification of killer defects in kesterite thin-film solar cells (2018)

• Wickramaratne et al, Iron as a source of efficient Shockley-Read-Hall recombination in GaN (2016)