### Abstract

Results from the application of our electrothermal simulator to n-type 0.15 µm gate In_{0.15}Ga_{0.85}As-Al_{0.28}Ga_{0.72}As HEMT structures are presented. The simulator involves an iterative procedure which alternately solves the Heat Diffusion Equation (HDE) and executes a Monte Carlo electronic transport algorithm. The net thermal flux generated during each Monte Carlo stage, calculated from the net rate of phonon emission, is fed into the thermal solution; the resulting temperature map is then used in the following Monte Carlo iteration. The HDE is solved through application of a novel analytical thermal resistance matrix technique which allows calculation of temperatures solely within the region of interest while including the large-scale boundary conditions. A novel charge injection scheme is applied for the treatment of side ohmic contacts, which avoids anomalous generation of thermal flux in adjacent regions. The characteristic 'thermal droop' is found in the I-V characteristics of the simulated device. Associated temperature distributions are shown to be spatially non-uniform with peak values and spatial locations dependent upon bias and the length of the containing die. Electron drift velocities and energies along the HEMT channel exhibit the largest shift on the inclusion of thermal self-consistency below the drain end of the gate, not at the location of the temperature peak.

Original language | English |
---|---|

Pages (from-to) | 207-211 |

Number of pages | 5 |

Journal | Journal of Computational Electronics |

Volume | 2 |

Issue number | 2-4 |

DOIs | |

Publication status | Published - Dec 2003 |

### Fingerprint

### Cite this

*Journal of Computational Electronics*,

*2*(2-4), 207-211. https://doi.org/10.1023/B:JCEL.0000011426.11111.64

**Electrothermal Monte Carlo Simulations of InGaAs/AlGaAs HEMTs.** / Pilgrim, N.J.; Batty, W.; Kelsall, R.W.

Research output: Contribution to journal › Article

*Journal of Computational Electronics*, vol. 2, no. 2-4, pp. 207-211. https://doi.org/10.1023/B:JCEL.0000011426.11111.64

}

TY - JOUR

T1 - Electrothermal Monte Carlo Simulations of InGaAs/AlGaAs HEMTs

AU - Pilgrim, N.J.

AU - Batty, W.

AU - Kelsall, R.W.

PY - 2003/12

Y1 - 2003/12

N2 - Results from the application of our electrothermal simulator to n-type 0.15 µm gate In0.15Ga0.85As-Al0.28Ga0.72As HEMT structures are presented. The simulator involves an iterative procedure which alternately solves the Heat Diffusion Equation (HDE) and executes a Monte Carlo electronic transport algorithm. The net thermal flux generated during each Monte Carlo stage, calculated from the net rate of phonon emission, is fed into the thermal solution; the resulting temperature map is then used in the following Monte Carlo iteration. The HDE is solved through application of a novel analytical thermal resistance matrix technique which allows calculation of temperatures solely within the region of interest while including the large-scale boundary conditions. A novel charge injection scheme is applied for the treatment of side ohmic contacts, which avoids anomalous generation of thermal flux in adjacent regions. The characteristic 'thermal droop' is found in the I-V characteristics of the simulated device. Associated temperature distributions are shown to be spatially non-uniform with peak values and spatial locations dependent upon bias and the length of the containing die. Electron drift velocities and energies along the HEMT channel exhibit the largest shift on the inclusion of thermal self-consistency below the drain end of the gate, not at the location of the temperature peak.

AB - Results from the application of our electrothermal simulator to n-type 0.15 µm gate In0.15Ga0.85As-Al0.28Ga0.72As HEMT structures are presented. The simulator involves an iterative procedure which alternately solves the Heat Diffusion Equation (HDE) and executes a Monte Carlo electronic transport algorithm. The net thermal flux generated during each Monte Carlo stage, calculated from the net rate of phonon emission, is fed into the thermal solution; the resulting temperature map is then used in the following Monte Carlo iteration. The HDE is solved through application of a novel analytical thermal resistance matrix technique which allows calculation of temperatures solely within the region of interest while including the large-scale boundary conditions. A novel charge injection scheme is applied for the treatment of side ohmic contacts, which avoids anomalous generation of thermal flux in adjacent regions. The characteristic 'thermal droop' is found in the I-V characteristics of the simulated device. Associated temperature distributions are shown to be spatially non-uniform with peak values and spatial locations dependent upon bias and the length of the containing die. Electron drift velocities and energies along the HEMT channel exhibit the largest shift on the inclusion of thermal self-consistency below the drain end of the gate, not at the location of the temperature peak.

U2 - 10.1023/B:JCEL.0000011426.11111.64

DO - 10.1023/B:JCEL.0000011426.11111.64

M3 - Article

VL - 2

SP - 207

EP - 211

JO - Journal of Computational Electronics

JF - Journal of Computational Electronics

SN - 1569-8025

IS - 2-4

ER -