We present a study of the nonlinear intersubband ͑ISB͒ response of conduction electrons in a GaAs/Al 0.3 Ga 0.7 As quantum well to strong terahertz ͑THz͒ radiation, using a density-matrix approach combined with time-dependent density-functional theory. We demonstrate coherent control of ISB optical bistability, using THz control pulses to induce picosecond switching between the bistable states. The switching speed is determined by the ISB relaxation and decoherence times, T 1 and T 2 . © 2004 American Institute of Physics. ͓DOI: 10.1063/1.1751611͔Intersubband ͑ISB͒ transitions in semiconductor quantum wells take place on a meV energy scale and are therefore attractive for terahertz ͑THz͒ device applications.1 Many ISB effects of practical interest occur in the nonlinear regime, such as second-and third-harmonic generation, 2 intensitydependent saturation of photo-absorption, 3,4 directional control over photocurrents, 5 generation of ultrashort THz pulses, 6 plasma instability, 7 or optical bistability. 8,9 Inspired by the photoabsorption experiments by Craig et al., 4 this letter presents a theoretical study of the optical bistability region in a strongly driven, modulation n-doped GaAs/Al 0.3 Ga 0.7 As quantum well. We demonstrate that ISB bistability can be manipulated on a picosecond time scale by short THz control pulses. This opens up opportunities for experimental study of optical bistability, which in the long run may lead to THz applications such as high-speed alloptical modulators and switches.Most previous theoretical studies of nonlinear ISB dynamics were based on density-matrix approaches, in Hartree approximation 10-12 or using exchange-only semiconductor Bloch equations. [13][14][15] These studies showed that the collective ISB electron dynamics is strongly influenced by depolarization and exchange-correlation ͑xc͒ many-body effects. 15 We will account for these effects using timedependent density-functional theory, which has the advantage of formal and computational simplicity.We describe the conduction subbands in effective-mass approximation, where m*ϭ0.067m and e*ϭe/ͱ⑀, ⑀ϭ13, are the effective mass and charge for GaAs. Initially, the conduction electrons are in the ground state, characterized by single-particle states of the form ⌿ jq ʈ 0 (r) ϭA Ϫ1/2 e iq ʈ r ʈ j 0 (z), with r ʈ and q ʈ in the x -y plane. The envelope function for the jth subband j 0 (z) follows selfconsistently from a one-dimensional Kohn-Sham ͑KS͒ equation, 16 with the ground-state densityHere, E jq ʈ ϭ j ϩប 2 q ʈ 2 /2m*, and j and F are the subband and Fermi energy levels. We choose the electronic sheet density N s such that only the lowest subband is occupied, in which case F ϭប 2 N s /m*ϩ 1 . Under the influence of THz driving fields, linearly polarized along z, the time-dependent states have the form ⌿ jq ʈ (r,t)ϭA Ϫ1/2 e iq ʈ r ʈ j (z,t), with initial condition ⌿ jq ʈ (r,t 0 )ϭ⌿ jq ʈ 0 (r). In the absence of disorder and phonons, the time evolution of the envelope functions follows from a time-dependent KS equationHer...