User Tools
qcl:working_principle
Differences
This shows you the differences between two versions of the page.
| Next revision | Previous revision | ||
|
qcl:working_principle [2017/07/24 12:51] thomas.grange created |
qcl:working_principle [2018/07/16 11:36] (current) stefan.birner |
||
|---|---|---|---|
| Line 1: | Line 1: | ||
| + | ====== Working principle ====== | ||
| The code is based on the non-equilibrium Green's functions (NEGF) formalism (also known as the Keldysh, or Kadanoff-Baym formalism). This formalism allows to account for both quantum transport effects (i.e. coherent transport effects, such as resonant tunneling), as well as scattering mechanisms. | The code is based on the non-equilibrium Green's functions (NEGF) formalism (also known as the Keldysh, or Kadanoff-Baym formalism). This formalism allows to account for both quantum transport effects (i.e. coherent transport effects, such as resonant tunneling), as well as scattering mechanisms. | ||
| Line 12: | Line 13: | ||
| In the beginning of the calculation, the single-band effective mass Schrödinger equation is solved in real space. | In the beginning of the calculation, the single-band effective mass Schrödinger equation is solved in real space. | ||
| - | The calculated energy levels and wave functions are then used as input to the NEGF algorithm. | + | The energy levels (i.e. minibands) are selected up to an energy cut-off which is specified in the input file by <code><Energy_Range_Axial unit="meV">150</Energy_Range_Axial> </code> |
| - | The wave functions are termed **modes** and the NEGF algorithm is written in terms of **mode space** and not **real space** to make it computationally more efficient. | + | The minibands are transformed into a localized basis of modes which we name "Reduced Real Space" basis. These wave functions are then used as a basis in the NEGF algorithm. |
| - | The number of QCL periods that are input to this Schrödinger equation are specified in ''<Number_of_lateral_periods_for_band_structure>''. | + | |
| - | In general, the core of the NEGF algorithm should be rather independent of this number, e.g. | + | As a second step, the scattering coupling terms are calculated for each of the accounted mechanism (optical and acoustic phonons, charged impurities, interface roughness, alloy disorder, electron-electron scattering). |
| - | <code> | + | |
| - | <Number_of_lateral_periods_for_band_structure> 4 | + | |
| - | </Number_of_lateral_periods_for_band_structure> | + | |
| - | </code> | + | |
| - | should lead to very similar results compared to a value of ''5'' but the numerical values might differ slightly. | + | |
| {{ :qcl:ldos_emin_emax_okay.jpg?200 |}} | {{ :qcl:ldos_emin_emax_okay.jpg?200 |}} | ||
| - | |||
| - | As a second step, the scattering coupling terms are calculated for each of the accounted mechanism (optical and acoustic phonons, charged impurities, interface roughness, alloy disorder). | ||
| Then, the main part of the calculation consists in the self-consistent NEGF solver. Starting from an initial guess of the Green's functions, the self-energies are calculated. The Green's functions are then calculated iteratively. | Then, the main part of the calculation consists in the self-consistent NEGF solver. Starting from an initial guess of the Green's functions, the self-energies are calculated. The Green's functions are then calculated iteratively. | ||
| Simultaneously, the mean-field electrostatic potential is calculted self-consistenly (Poisson's equation). Such iterations are made until convergence is reached for the Green's functions as well as for the calculated current. | Simultaneously, the mean-field electrostatic potential is calculted self-consistenly (Poisson's equation). Such iterations are made until convergence is reached for the Green's functions as well as for the calculated current. | ||
| - | From the Green's functions solution, the gain is then calculated (if requested in the input file). | + | From the Green's functions solution, the current density is directly calculated, as well as as the carrier density. To gain more insights into analysis, the populations, the density matrix, the oscillator strengths are displayed in different basis. |
| + | |||
| + | The gain is then calculated if requested in the input file. To this purpose, for each photon energy specified in the input file, an ac electromagnetic perturbation is considered. An additional self-consistent routine is used is to calculate the linear response of the Green's functions to this perturbation. The gain (or absorption) spectrum is then obtained from this linear response. | ||
qcl/working_principle.1500900689.txt.gz · Last modified: 2017/07/24 12:51 by thomas.grange
