Intrinsic high-frequency oscillations (^2.5 THz) in current and corresponding quantum-well density have been simulated for the first time for a fixed-bias voltage in the negative differential resistance (NDR) region of the current-voltage (I-V) characteristics of a resonant tunneling diode. Scattering and self-consistency are included. Hysteresis and "plateaulike" behavior of the time-averaged I-V curve are simulated in the NDR region. Intrinsic bistability is manifested by the phenomenon of unstable electron charge buildup and ejection from the quantum well.PACS numbers: 73.40.Gk The resonant tunneling of electrons through doublebarrier structures is known to produce a region of negative differential resistance (NDR) in the I-V relationship. This has been well established experimentally, 1 " 5 and a simple understanding of the phenomena may be obtained through analyzing the transmission coefficient of the one-dimensional resonant tunneling diode (RTD) profile. 6 " 9 Transmission-coefficient-based calculations can include self-consistency and duplicate gross features of RTD I-V curves, but they are impractical for simulating time-dependent nonlinear effects, such as occur when the bias is suddenly switched (transient response or frequency-dependent behavior) and other time-dependent nonlinear phenomena such as the high-frequency oscillations at fixed bias in the NDR region simulated in this paper.As the full characterization of RTD devices requires simulating a time-evolving system, the Wigner distribution function 10 (WDF) approach has emerged as the most promising method. 11 " 15 Scattering within the device and self-consistency 13,15 " 17 may be included, and the WDF approach has the added advantage of lending itself to a "particle" interpretation. 15,16 ' 18 We present here a more accurate and fully time-dependent and self-consistent quantum transport simulation 15 which includes scattering as a function of temperature and which demonstrates current oscillations at a fixed bias. We note, however, that other time-dependent and selfconsistent scattering calculations have been reported prior to the current work, 12 " 14 but they differ from ours in some key respects. 19 The method we have chosen to incorporate self-consistency and scattering allows us to address more accurately the controversial "bistability" issue reported 4,20 and challenged 21 previously.The experimentally reported I-V curves exhibit a characteristic "plateaulike" behavior 5,22 and hysteresis. 20,21 The first view proposed by Goldman, Tsui, and Cunningham 4 attributed this entire behavior to "intrinsic bistability," that is, charge is dynamically built up and ejected from the well, and through an electrostatic feedback mechanism, two current states exist. Sollner, 21 who represents the second view, disagrees, claiming that instead of Goldman's hysteresis being due to a charging of the well, it is due to external circuit-induced oscillations of the current in the NDR region of the device. We have observed intrinsic oscillations (which have a f...