Chaotic motion of space charge wavefronts in semiconductors under time-independent voltage bias

A standard drift-diffusion model of space charge wave propagation in semiconductors has been studied numerically and analytically under dc voltage bias. For sufficiently long samples, appropriate contact resistivity and applied voltage - such that the sample is biased in a regime of negative differential resistance - we find chaos in the propagation of nonlinear fronts (charge monopoles of alternating sign) of electric field. The chaos is always low-dimensional, but has a complex spatial structure; this behavior can be interpreted using a finite dimensional asymptotic model in which the front (charge monopole) positions and the electrical current are the only dynamical variables.Comment: 12 pages, 8 figure

Chapman-Enskog method and synchronization of globally coupled oscillators

The Chapman-Enskog method of kinetic theory is applied to two problems of synchronization of globally coupled phase oscillators. First, a modified Kuramoto model is obtained in the limit of small inertia from a more general model which includes ``inertial'' effects. Second, a modified Chapman-Enskog method is used to derive the amplitude equation for an O(2) Takens-Bogdanov bifurcation corresponding to the tricritical point of the Kuramoto model with a bimodal distribution of oscillator natural frequencies. This latter calculation shows that the Chapman-Enskog method is a convenient alternative to normal form calculations.Comment: 7 pages, 2-column Revtex, no figures, minor change

Dynamics of Electric Field Domains and Oscillations of the Photocurrent in a Simple Superlattice Model

A discrete model is introduced to account for the time-periodic oscillations of the photocurrent in a superlattice observed by Kwok et al, in an undoped 40 period AlAs/GaAs superlattice. Basic ingredients are an effective negative differential resistance due to the sequential resonant tunneling of the photoexcited carriers through the potential barriers, and a rate equation for the holes that incorporates photogeneration and recombination. The photoexciting laser acts as a damping factor ending the oscillations when its power is large enough. The model explains: (i) the known oscillatory static I-V characteristic curve through the formation of a domain wall connecting high and low electric field domains, and (ii) the photocurrent and photoluminescence time-dependent oscillations after the domain wall is formed. In our model, they arise from the combined motion of the wall and the shift of the values of the electric field at the domains. Up to a certain value of the photoexcitation, the non-uniform field profile with two domains turns out to be metastable: after the photocurrent oscillations have ceased, the field profile slowly relaxes toward the uniform stationary solution (which is reached on a much longer time scale). Multiple stability of stationary states and hysteresis are also found. An interpretation of the oscillations in the photoluminescence spectrum is also given.Comment: 34 pages, REVTeX 3.0, 10 figures upon request, MA/UC3M/07/9

Current-voltage characteristic and stability in resonant-tunneling n-doped semiconductor superlattices

We review the occurrence of electric-field domains in doped superlattices within a discrete drift model. A complete analysis of the construction and stability of stationary field profiles having two domains is carried out. As a consequence, we can provide a simple analytical estimation for the doping density above which stable stable domains occur. This bound may be useful for the design of superlattices exhibiting self-sustained current oscillations. Furthermore we explain why stable domains occur in superlattices in contrast to the usual Gunn diode.Comment: Tex file and 3 postscript figure

Stationary states and phase diagram for a model of the Gunn effect under realistic boundary conditions

A general formulation of boundary conditions for semiconductor-metal contacts follows from a phenomenological procedure sketched here. The resulting boundary conditions, which incorporate only physically well-defined parameters, are used to study the classical unipolar drift-diffusion model for the Gunn effect. The analysis of its stationary solutions reveals the presence of bistability and hysteresis for a certain range of contact parameters. Several types of Gunn effect are predicted to occur in the model, when no stable stationary solution exists, depending on the value of the parameters of the injecting contact appearing in the boundary condition. In this way, the critical role played by contacts in the Gunn effect is clearly stablished.Comment: 10 pages, 6 Post-Script figure

Chaotic dynamics of electric-field domains in periodically driven superlattices

Self-sustained time-dependent current oscillations under dc voltage bias have been observed in recent experiments on n-doped semiconductor superlattices with sequential resonant tunneling. The current oscillations are caused by the motion and recycling of the domain wall separating low- and high-electric- field regions of the superlattice, as the analysis of a discrete drift model shows and experimental evidence supports. Numerical simulation shows that different nonlinear dynamical regimes of the domain wall appear when an external microwave signal is superimposed on the dc bias and its driving frequency and driving amplitude vary. On the frequency - amplitude parameter plane, there are regions of entrainment and quasiperiodicity forming Arnol'd tongues. Chaos is demonstrated to appear at the boundaries of the tongues and in the regions where they overlap. Coexistence of up to four electric-field domains randomly nucleated in space is detected under ac+dc driving.Comment: 9 pages, LaTex, RevTex. 12 uuencoded figures (1.8M) should be requested by e-mail from the autho

Influence of primary particle density in the morphology of agglomerates

Agglomeration processes occur in many different realms of science such as colloid and aerosol formation or formation of bacterial colonies. We study the influence of primary particle density in agglomerate structure using diffusion-controlled Monte Carlo simulations with realistic space scales through different regimes (DLA and DLCA). The equivalence of Monte Carlo time steps to real time scales is given by Hirsch's hydrodynamical theory of Brownian motion. Agglomerate behavior at different time stages of the simulations suggests that three indices (fractal exponent, coordination number and eccentricity index) characterize agglomerate geometry. Using these indices, we have found that the initial density of primary particles greatly influences the final structure of the agglomerate as observed in recent experimental works.Comment: 11 pages, 13 figures, PRE, to appea