A set of codes to compute carrier capture and non-radiative recombination rates associated with point defects in semiconducting compounds.
Multiphonon process involving impurities 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.

Installation

The codes are written in Julia, while the scripts and Jupyter Notebooks also contain Python and use pymatgen and pawpyseed, 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 consists of several steps, implemented in a series of programs, which may be run from the command line. Input for the calculations is provided in input.yaml.

0. 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.

1. 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.

2. 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.

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 and consider using the doped package).

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.

• SnZn in Cu2ZnSnS4: 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 freedom. 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):

Here, V is the volume of the supercell, Wif is the electron-phonon overlap and ξim and ξfn describe the wavefunctions of the mth and nth 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 σ.

Citation

@article{kim2020carriercapture,
  title={Carriercapture. jl: Anharmonic carrier capture},
  author={Kim, Sunghyun and Hood, Samantha N and van Gerwen, Puck and Whalley, Lucy D and Walsh, Aron},
  journal={Journal of Open Source Software},
  volume={5},
  number={47},
  pages={2102},
  year={2020},
  doi={10.21105/joss.02102},
  url={https://joss.theoj.org/papers/10.21105/joss.02102},
}

Extended Reading List

Theory Development

• Henry 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 of CarrierCapture

• Wang et al, Upper efficiency limit of Sb2Se3 solar cells (2024)

• Kavanagh et al, Impact of metastable defect structures on carrier recombination in solar cells (2022)

• Kavanagh et al, Rapid recombination by cadmium vacancies in CdTe (2021)

• Whalley et al, Giant Huang–Rhys factor for electron capture by the iodine intersitial in perovskite solar cells (2021)

• Kim and Walsh, Ab initio calculation of the detailed balance limit to the photovoltaic efficiency of single p-n junction kesterite solar cells (2021)

• Dahliah et al, High-throughput computational search for high carrier lifetime, defect-tolerant solar absorbers (2021)

• Kim et al, Upper limit to the photovoltaic efficiency of imperfect crystals (2020)

• 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)