GNU Archimedes

Last updated
GNU Archimedes
Original author(s) Jean Michel Sellier
Developer(s) GNU project
Stable release
2.0.1 [1]   OOjs UI icon edit-ltr-progressive.svg / 30 April 2013
Operating system Linux, UNIX
Type TCAD
License 2007: GPL-3.0-or-later [lower-alpha 1]
2004: GPL-2.0-or-later [lower-alpha 2]
Website gnu.org/software/archimedes/

Archimedes is a TCAD package for use by engineers to design and simulate submicron and mesoscopic semiconductor devices. Archimedes is free software and thus it can be copied, modified and redistributed under GPL. Archimedes uses the Ensemble Monte Carlo method and is able to simulate physics effects and transport for electrons and heavy holes in Silicon, Germanium, GaAs, InSb, AlSb, AlAs, AlxInxSb, AlxIn(1-x)Sb, AlP, AlSb, GaP, GaSb, InP and their compounds (III-V semiconductor materials), along with Silicon Oxide. Applied and/or self-consistent electrostatic and magnetic fields are handled with the Poisson and Faraday equations.

Contents

The GNU project has announced in May, 2012 that the software package Aeneas [2] will be substituted by Archimedes, making this one the GNU package for Monte Carlo semiconductor devices simulations. [3]

Introduction

Archimedes is the GNU package for semiconductor device simulations that has been released for the first time on 2005 under GPL. It has been created by Jean Michel Sellier who is, since then, the leader of the project and the main developer. It is a free software and thus it can be copied, modified and redistributed under GPL.

Archimedes belongs to the well-known family of TCAD software, i.e. tools utilized to assist the development of technologically relevant products. In particular, this package assists engineers in designing and simulating submicron and mesoscopic semiconductor devices. In a next-future[ timeframe? ] version Archimedes will also be able to simulate nanodevices, using the Wigner Monte Carlo formalism [4] (an experimental release can be found at [5] ). Today Archimedes is used in several big companies for simulation and production purposes.

Archimedes is also useful for teaching purposes since everybody can access the sources, modify and test them. Today, it is used for teaching courses in several hundreds universities all around the world. Furthermore, a simplified version, developed for students, is available on nanoHUB.org.

The Ensemble Monte Carlo method is the method that Archimedes uses to simulate and predict the behavior of a devices. Being the Monte Carlo very stable and reliable, Archimedes can be used to know the characteristics of a device even before this last is built.

The physics and geometry of a device is described simply by a script, which makes, in this sense, Archimedes a powerful tool for the simulation of quite general semiconductor devices. [6]

Archimedes is able to simulate a plenty of physics effects and transport for electrons and heavy holes in Silicon, Germanium, GaAs, InSb, AlSb, AlAs, AlxInxSb, AlxIn(1-x)Sb, AlP, AlSb, GaP, GaSb, InP and their compounds (III-V semiconductor materials), along with Silicon Oxide, the applied and/or self-consistent electrostatic and magnetic fields by means of Poisson and Faraday equation. It is also able to deal with heterostructures.

Boltzmann transport equation

The Boltzmann transport equation model has been the main tool used in the analysis of transport in semiconductors. The BTE equation is given by:

The distribution function, f, is a dimensionless function which is used to extract all observables of interest and gives a full depiction of electron distribution in both real and k-space. Further, it physically represents the probability of particle occupation of energy k at position r and time t. In addition, due to being a seven-dimensional integro-differential equation (six dimensions in the phase space and one in time) the solution to the BTE is cumbersome and can be solved in closed analytical form under very special restrictions. Numerically, solution to the BTE is employed using either a deterministic method or a stochastic method. The deterministic method solution is based on a grid-based numerical method such as the spherical harmonics approach, whereas the Monte Carlo is the stochastic approach used to solve the BTE.

Monte Carlo method

The semiclassical Monte Carlo method is a statistical method used to yield exact solution to the Boltzmann transport equation which includes complex band structure and scattering processes. This approach is semiclassical for the reason that scattering mechanisms are treated quantum mechanically using the Fermi's Golden Rule, whereas the transport between scattering events is treated using the classical particle notion. The Monte Carlo model in essence tracks the particle trajectory at each free flight and chooses a corresponding scattering mechanism stochastically. Two of the great advantages of semiclassical Monte Carlo are its capability to provide accurate quantum mechanical treatment of various distinct scattering mechanisms within the scattering terms, and the absence of assumption about the form of carrier distribution in energy or k-space. The semiclassical equation describing the motion of an electron is

where F is the electric field, E(k) is the energy dispersion relation, and k is the momentum wave vector. To solve the above equation, one needs strong knowledge of the band structure (E(k)). The E(k) relation describes how the particle moves inside the device, in addition to depicting useful information necessary for transport such as the density of states (DOS) and the particle velocity. A Full-band E(K) relation can be obtained using the semi-empirical pseudopotential method. [7]

Screenshots

Notes

  1. GPL-3.0-or-later since 2007-08-25 (0.0.8).
  2. GPL-2.0-or-later from 2004-10-12 (0.0.1) until 2007-01-10 (0.0.7).

Related Research Articles

In solid state physics, a particle's effective mass is the mass that it seems to have when responding to forces, or the mass that it seems to have when interacting with other identical particles in a thermal distribution. One of the results from the band theory of solids is that the movement of particles in a periodic potential, over long distances larger than the lattice spacing, can be very different from their motion in a vacuum. The effective mass is a quantity that is used to simplify band structures by modeling the behavior of a free particle with that mass. For some purposes and some materials, the effective mass can be considered to be a simple constant of a material. In general, however, the value of effective mass depends on the purpose for which it is used, and can vary depending on a number of factors.

<span class="mw-page-title-main">Quantum well</span> Concept in quantum mechanics

A quantum well is a potential well with only discrete energy values.

A heterojunction is an interface between two layers or regions of dissimilar semiconductors. These semiconducting materials have unequal band gaps as opposed to a homojunction. It is often advantageous to engineer the electronic energy bands in many solid-state device applications, including semiconductor lasers, solar cells and transistors. The combination of multiple heterojunctions together in a device is called a heterostructure, although the two terms are commonly used interchangeably. The requirement that each material be a semiconductor with unequal band gaps is somewhat loose, especially on small length scales, where electronic properties depend on spatial properties. A more modern definition of heterojunction is the interface between any two solid-state materials, including crystalline and amorphous structures of metallic, insulating, fast ion conductor and semiconducting materials.

In solid-state physics, the electron mobility characterises how quickly an electron can move through a metal or semiconductor when pushed or pulled by an electric field. There is an analogous quantity for holes, called hole mobility. The term carrier mobility refers in general to both electron and hole mobility.

A superlattice is a periodic structure of layers of two materials. Typically, the thickness of one layer is several nanometers. It can also refer to a lower-dimensional structure such as an array of quantum dots or quantum wells.

The Compton wavelength is a quantum mechanical property of a particle, defined as the wavelength of a photon the energy of which is the same as the rest energy of that particle. It was introduced by Arthur Compton in 1923 in his explanation of the scattering of photons by electrons.

Quantum-cascade lasers (QCLs) are semiconductor lasers that emit in the mid- to far-infrared portion of the electromagnetic spectrum and were first demonstrated by Jérôme Faist, Federico Capasso, Deborah Sivco, Carlo Sirtori, Albert Hutchinson, and Alfred Cho at Bell Laboratories in 1994.

The Gross–Pitaevskii equation describes the ground state of a quantum system of identical bosons using the Hartree–Fock approximation and the pseudopotential interaction model.

The beam propagation method (BPM) is an approximation technique for simulating the propagation of light in slowly varying optical waveguides. It is essentially the same as the so-called parabolic equation (PE) method in underwater acoustics. Both BPM and the PE were first introduced in the 1970s. When a wave propagates along a waveguide for a large distance, rigorous numerical simulation is difficult. The BPM relies on approximate differential equations which are also called the one-way models. These one-way models involve only a first order derivative in the variable z and they can be solved as "initial" value problem. The "initial" value problem does not involve time, rather it is for the spatial variable z.

In theoretical physics, the eikonal approximation is an approximative method useful in wave scattering equations, which occur in optics, seismology, quantum mechanics, quantum electrodynamics, and partial wave expansion.

<span class="mw-page-title-main">Radiative transfer equation and diffusion theory for photon transport in biological tissue</span>

Photon transport in biological tissue can be equivalently modeled numerically with Monte Carlo simulations or analytically by the radiative transfer equation (RTE). However, the RTE is difficult to solve without introducing approximations. A common approximation summarized here is the diffusion approximation. Overall, solutions to the diffusion equation for photon transport are more computationally efficient, but less accurate than Monte Carlo simulations.

Defect types include atom vacancies, adatoms, steps, and kinks that occur most frequently at surfaces due to the finite material size causing crystal discontinuity. What all types of defects have in common, whether surface or bulk defects, is that they produce dangling bonds that have specific electron energy levels different from those of the bulk. This difference occurs because these states cannot be described with periodic Bloch waves due to the change in electron potential energy caused by the missing ion cores just outside the surface. Hence, these are localized states that require separate solutions to the Schrödinger equation so that electron energies can be properly described. The break in periodicity results in a decrease in conductivity due to defect scattering.

Miniaturizing components has always been a primary goal in the semiconductor industry because it cuts production cost and lets companies build smaller computers and other devices. Miniaturization, however, has increased dissipated power per unit area and made it a key limiting factor in integrated circuit performance. Temperature increase becomes relevant for relatively small-cross-sections wires, where it may affect normal semiconductor behavior. Besides, since the generation of heat is proportional to the frequency of operation for switching circuits, fast computers have larger heat generation than slow ones, an undesired effect for chips manufacturers. This article summaries physical concepts that describe the generation and conduction of heat in an integrated circuit, and presents numerical methods that model heat transfer from a macroscopic point of view.

In condensed matter physics, biexcitons are created from two free excitons.

The transport of heat in solids involves both electrons and vibrations of the atoms (phonons). When the solid is perfectly ordered over hundreds of thousands of atoms, this transport obeys established physics. However, when the size of the ordered regions decreases new physics can arise, thermal transport in nanostructures. In some cases heat transport is more effective, in others it is not.

The Monte Carlo method for electron transport is a semiclassical Monte Carlo (MC) approach of modeling semiconductor transport. Assuming the carrier motion consists of free flights interrupted by scattering mechanisms, a computer is utilized to simulate the trajectories of particles as they move across the device under the influence of an electric field using classical mechanics. The scattering events and the duration of particle flight is determined through the use of random numbers.

Biology Monte Carlo methods (BioMOCA) have been developed at the University of Illinois at Urbana-Champaign to simulate ion transport in an electrolyte environment through ion channels or nano-pores embedded in membranes. It is a 3-D particle-based Monte Carlo simulator for analyzing and studying the ion transport problem in ion channel systems or similar nanopores in wet/biological environments. The system simulated consists of a protein forming an ion channel (or an artificial nanopores like a Carbon Nano Tube, CNT), with a membrane (i.e. lipid bilayer) that separates two ion baths on either side. BioMOCA is based on two methodologies, namely the Boltzmann transport Monte Carlo (BTMC) and particle-particle-particle-mesh (P3M). The first one uses Monte Carlo method to solve the Boltzmann equation, while the later splits the electrostatic forces into short-range and long-range components.

McXtrace is an open source software package for performing Monte Carlo simulations of X-ray scattering experiments. While its chief objective is to aid in the optimization of beamlines at e.g. synchrotrons, it may also be used for data analysis and at laboratory sources and beamlines. McXtrace is free software released under the GNU GPL.

Heat transfer physics describes the kinetics of energy storage, transport, and energy transformation by principal energy carriers: phonons, electrons, fluid particles, and photons. Heat is thermal energy stored in temperature-dependent motion of particles including electrons, atomic nuclei, individual atoms, and molecules. Heat is transferred to and from matter by the principal energy carriers. The state of energy stored within matter, or transported by the carriers, is described by a combination of classical and quantum statistical mechanics. The energy is different made (converted) among various carriers. The heat transfer processes are governed by the rates at which various related physical phenomena occur, such as the rate of particle collisions in classical mechanics. These various states and kinetics determine the heat transfer, i.e., the net rate of energy storage or transport. Governing these process from the atomic level to macroscale are the laws of thermodynamics, including conservation of energy.

The semiconductor Bloch equations describe the optical response of semiconductors excited by coherent classical light sources, such as lasers. They are based on a full quantum theory, and form a closed set of integro-differential equations for the quantum dynamics of microscopic polarization and charge carrier distribution. The SBEs are named after the structural analogy to the optical Bloch equations that describe the excitation dynamics in a two-level atom interacting with a classical electromagnetic field. As the major complication beyond the atomic approach, the SBEs must address the many-body interactions resulting from Coulomb force among charges and the coupling among lattice vibrations and electrons.

References

  1. "Archimedes 2.0.1 has been released". 30 April 2013.
  2. « Aeneas », gnu.org, May 2012.
  3. Sellier, Jean Michel (2012-05-13). "Aeneas new release" (Mailing list). info-gnu. Retrieved 2012-05-13.
  4. E. Wigner, On the Quantum Correction for Thermodynamic Equilibrium (1932)
  5. J.M. Sellier, http://www.nano-archimedes.com
  6. "ALT Linux - p9 - archimedes-2.0.1-alt1_11 - 2D Quantum Monte Carlo simulator for semiconductor devices". packages.altlinux.org. Retrieved 2023-06-20.
  7. K. Hess, Monte Carlo Device Simulation: Full Band and Beyond, Technology (1991)