BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

L. L. Bonilla

L. L. Bonilla contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechcond-mat.mes-hallBiological Physicsnlin.PScond-mat.softmath-phmath.MPcond-mat.mtrl-sciTissues and OrgansBiomoleculescond-mat.dis-nnnlin.AOnlin.CDphysics.med-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
29works
0followers
14topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2023arXiv

Scale free chaos in the confined Vicsek flocking model

The Vicsek model encompasses the paradigm of active dry matter. Motivated by collective behavior of insects in swarms, we have studied finite size effects and criticality in the three dimensional, harmonically confined Vicsek model. We have discovered a phase transition that exists for appropriate noise and small confinement strength. On the critical line of confinement versus noise, swarms are in a state of scale-free chaos characterized by minimal correlation time, correlation length proportional to swarm size and topological data analysis. The critical line separates dispersed single clusters from confined multicluster swarms. Scale-free chaotic swarms occupy a compact region of space and comprise a recognizable `condensed' nucleus and particles leaving and entering it. Susceptibility, correlation length, dynamic correlation function and largest Lyapunov exponent obey power laws. The critical line and a narrow criticality region close to it move simultaneously to zero confinement strength for infinitely many particles. At the end of the first chaotic window of confinement, there is another phase transition to infinitely dense clusters of finite size that may be termed flocking black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bifurcation theory captures band formation in the Vicsek model of flock formation

Collective behavior occurs ubiquitously in nature and it plays a key role in bacterial colonies, mammalian cells or flocks of birds. Here, we examine the average density and velocity of self-propelled particles, which are described by a system of partial differential equations near the flocking transition of the Vicsek model. This agent-based model illustrates the trend towards flock formation of animals that align their velocities to an average of those of their neighbors. Near the flocking transition, particle density and velocity obey partial differential equations that include a parameter $ε$ measuring the distance to the bifurcation point. We have obtained analytically the Riemann invariants in one and two spatial dimensions for the hyperbolic ($ε=0$) and parabolic ($ε\neq 0$) system and, under periodic initial-boundary value conditions, we show that the solutions include wave trains. Additionally, we have found wave trains having oscillation frequencies that agree with those predicted by a linearization approximation and that may propagate at angles depending on the initial condition. The wave amplitudes increase with time for the hyperbolic system but are stabilized to finite values for the parabolic system. To integrate the partial differential equations, we design a basic numerical scheme which is first order in time and space. To mitigate numerical dissipation and ensure good resolution of the wave features, we also use a high order accurate WENO5 reconstruction procedure in space and a third order accurate Runge-Kutta scheme in time. Comparisons with direct simulations of the Vicsek model confirm these predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two dimensional soliton in tumor induced angiogenesis

Ensemble averages of a stochastic model show that, after a formation stage, the tips of active blood vessels in an angiogenic network form a moving two dimensional stable diffusive soliton, which advances toward sources of growth factor. Here we use methods of multiple scales to find the diffusive soliton as a solution of a deterministic equation for the mean density of active endothelial cells tips. We characterize the diffusive soliton shape in a general geometry, and find that its vector velocity and the trajectory of its center of mass along curvilinear coordinates solve appropriate collective coordinate equations. The vessel tip density predicted by the soliton compares well with that obtained by ensemble averages of simulations of the stochastic model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic Model of Tumor-induced Angiogenesis: Ensemble Averages and Deterministic Equations

A recent conceptual model of tumor-driven angiogenesis including branching, elongation, and anastomosis of blood vessels captures some of the intrinsic multiscale structures of this complex system, yet allowing to extract a deterministic integro-partial differential description of the vessel tip density [Phys. Rev. E 90, 062716 (2014)]. Here we solve the stochastic model, show that ensemble averages over many realizations correspond to the deterministic equations, and fit the anastomosis rate coefficient so that the total number of vessel tips evolves similarly in the deterministic and ensemble averaged stochastic descriptions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Wavelength selection of rippling patterns in myxobacteria

Rippling patterns of myxobacteria appear in starving colonies before they aggregate to form fruiting bodies. These periodic traveling cell density waves arise from the coordination of individual cell reversals, resulting from an internal clock regulating them, and from contact signaling during bacterial collisions. Here we revisit a mathematical model of rippling in myxobacteria due to Igoshin et al.\ [Proc. Natl. Acad. Sci. USA {\bf 98}, 14913 (2001) and Phys. Rev. E {\bf 70}, 041911 (2004)]. Bacteria in this model are phase oscillators with an extra internal phase through which they are coupled to a mean-field of oppositely moving bacteria. Previously, patterns for this model were obtained only by numerical methods and it was not possible to find their wavenumber analytically. We derive an evolution equation for the reversal point density that selects the pattern wavenumber in the weak signaling limit, show the validity of the selection rule by solving numerically the model equations and describe other stable patterns in the strong signaling limit. The nonlocal mean-field coupling tends to decohere and confine patterns. Under appropriate circumstances, it can annihilate the patterns leaving a constant density state via a nonequilibrium phase transition reminiscent of destruction of synchronization in the Kuramoto model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measuring strain and rotation fields at the dislocation core in graphene

Strain fields, dislocations and defects may be used to control electronic properties of graphene. By using advanced imaging techniques with high-resolution transmission electron microscopes, we have measured the strain and rotation fields about dislocations in monolayer graphene with single-atom sensitivity. These fields differ qualitatively from those given by conventional linear elasticity. However, atom positions calculated from two dimensional (2D) discrete elasticity and three dimensional discrete periodized Föppl-von Kármán equations (dpFvKEs) yield fields close to experiments when determined by geometric phase analysis. 2D theories produce symmetric fields whereas those from experiments exhibit asymmetries. Numerical solutions of dpFvKEs provide strain and rotation fields of dislocation dipoles and pairs that also exhibit asymmetries and, compared with experiments, may yield information on out-of-plane displacements of atoms. While discrete theories need to be solved numerically, analytical formulas for strains and rotation about dislocations can be obtained from 2D Mindlin's hyperstress theory. These formulas are very useful for fitting experimental data and provide a template to ascertain the importance of nonlinear and nonplanar effects. Measuring the parameters of this theory, we find two characteristic lengths between three and four times the lattice spacings that control dilatation and rotation about a dislocation. At larger distances from the dislocation core, the elastic fields decay to those of conventional elasticity. Our results may be relevant for strain engineering in graphene and other 2D materials of current interest.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of optically injected currents in carbon nanotubes

We consider theoretically the dynamics of electric currents optically injected in carbon nanotubes. Although the plasma oscillations are not seen in these systems, the main effect on the carrier's motion is due to strongly nonuniform space-charge Coulomb forces produced by time-dependent separation of injected electron and hole densities. We calculate evolution of the dipole moment characterizing the time- and coordinate-dependent charge density distributions and analyze different regimes of the dynamics. The developed time-dependent dipole moment leads to a dipole radiation in the THz frequency range for typical parameters of injected currents.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hybrid modeling of tumor-induced angiogenesis

When modeling of tumor-driven angiogenesis, a major source of analytical and computational complexity is the strong coupling between the kinetic parameters of the relevant stochastic branching-and-growth of the capillary network, and the family of interacting underlying fields. To reduce this complexity, we take advantage of the system intrinsic multiscale structure: we describe the stochastic dynamics of the cells at the vessel tip at their natural mesoscale, whereas we describe the deterministic dynamics of the underlying fields at a larger macroscale. Here, we set up a conceptual stochastic model including branching, elongation, and anastomosis of vessels and derive a mean field approximation for their densities. This leads to a deterministic integro-partial differential system that describes the formation of the stochastic vessel network. We discuss the proper capillary injecting boundary conditions and include the results of relevant numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Noise enhanced spontaneous chaos in semiconductor superlattices at room temperature

Physical systems exhibiting fast spontaneous chaotic oscillations are used to generate high-quality true random sequences in random number generators. The concept of using fast practical entropy sources to produce true random sequences is crucial to make storage and transfer of data more secure at very high speeds. While the first high-speed devices were chaotic semiconductor lasers, the discovery of spontaneous chaos in semiconductor superlattices at room temperature provides a valuable nanotechnology alternative. Spontaneous chaos was observed in 1996 experiments at temperatures below liquid nitrogen. Here we show spontaneous chaos at room temperature appears in idealized superlattices for voltage ranges where sharp transitions between different oscillation modes occur. Internal and external noises broaden these voltage ranges and enhance the sensitivity to initial conditions in the superlattice snail-shaped chaotic attractor thereby rendering spontaneous chaos more robust.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sawtooth patterns in biomolecules force-extension curves: an equilibrium-statistical-mechanics theory

We analyze the force-extension curve for a general class of systems, which are described at the mesoscopic level by a free energy depending on the extension of its components. Similarly to what is done in real experiments, the total length of the system is the controlled parameter. This imposes a global constraint in the minimization procedure leading to the equilibrium values of the extensions. As a consequence, the force-extension curve has multiple branches in a certain range of forces. The stability of these branches is governed by the free energy: there are a series of first-order phase transitions at certain values of the total length, in which the free energy itself is continuous but its first derivative, the force, has a finite jump. This behavior is completely similar to the one observed in real experiments with biomolecules like proteins, and other complex systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Model of ripples in graphene

We propose a model of ripples in suspended graphene sheets based on plate equations that are made discrete with the periodicity of the honeycomb lattice and then periodized. In addition, the equation for the displacements with respect to the planar configuration contains a double-well site potential, a nonlinear friction and a multiplicative white noise term satisfying the fluctuation-dissipation theorem. The nonlinear friction terms agree with those proposed by Eichler et al [Nature Nanotech. {\bf 6}, 339 (2011)] to explain their experiments with a graphene resonator. The site double-well potential indicates that the carbon atoms at each lattice point have equal probability to move upward or downward off-plane. For the considered parameter values, the relaxation time due to friction is much larger than the periods of membrane vibrations and the noise is quite small. Then ripples with no preferred orientation appear as long-lived metastable states for any temperature. Numerical solutions confirm this picture.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical method for hydrodynamic modulation equations describing Bloch oscillations in semiconductor superlattices

We present a finite difference method to solve a new type of nonlocal hydrodynamic equations that arise in the theory of spatially inhomogeneous Bloch oscillations in semiconductor superlattices. The hydrodynamic equations describe the evolution of the electron density, electric field and the complex amplitude of the Bloch oscillations for the electron current density and the mean energy density. These equations contain averages over the Bloch phase which are integrals of the unknown electric field and are derived by singular perturbation methods. Among the solutions of the hydrodynamic equations, at a 70 K lattice temperature, there are spatially inhomogeneous Bloch oscillations coexisting with moving electric field domains and Gunn-type oscillations of the current. At higher temperature (300 K) only Bloch oscillations remain. These novel solutions are found for restitution coefficients in a narrow interval below their critical values and disappear for larger values. We use an efficient numerical method based on an implicit second-order finite difference scheme for both the electric field equation (of drift-diffusion type) and the parabolic equation for the complex amplitude. Double integrals appearing in the nonlocal hydrodynamic equations are calculated by means of expansions in modified Bessel functions. We use numerical simulations to ascertain the convergence of the method. If the complex amplitude equation is solved using a first order scheme for restitution coefficients near their critical values, a spurious convection arises that annihilates the complex amplitude in the part of the superlattice that is closer to the cathode. This numerical artifact disappears if the space step is appropriately reduced or we use the second-order numerical scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ripples in a graphene membrane coupled to Glauber spins

We propose a theory of ripples in suspended graphene sheets based on two-dimensional elasticity equations that are made discrete on the honeycomb lattice and then periodized. At each point carbon atoms are coupled to Ising spins whose values indicate the atoms local trend to move vertically off-plane. The Ising spins are in contact with a thermal bath and evolve according to Glauber dynamics. In the limit of slow spin flip compared to membrane vibrations, ripples with no preferred orientation appear as long-lived metastable states for any temperature. Numerical solutions confirm this picture.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spin-oscillator model for DNA/RNA unzipping by mechanical force

We model unzipping of DNA/RNA molecules subject to an external force by a spin-oscillator system. The system comprises a macroscopic degree of freedom, represented by a one-dimensional oscillator, and internal degrees of freedom, represented by Glauber spins with nearest-neighbor interaction and a coupling constant proportional to the oscillator position. At a critical value $F_c$ of an applied external force $F$, the oscillator rest position (order parameter) changes abruptly and the system undergoes a first-order phase transition. When the external force is cycled at different rates, the extension given by the oscillator position exhibits a hysteresis cycle at high loading rates whereas it moves reversibly over the equilibrium force-extension curve at very low loading rates. Under constant force, the logarithm of the residence time at the stable and metastable oscillator rest position is proportional to $(F-F_c)$ as in an Arrhenius law.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Strain and rotation fields of dislocations in graphene

Strain and rotation fields of dislocations in monolayer graphene have been mapped in a recent experiment. These fields are finite everywhere and differ from those given by linear elasticity which does not consider rotation explicitly and predicts infinite rotation and strains at the dislocation point. A derivative regularization of two-dimensional linear elasticity fixes these shortcomings. The theory adds rotation, dislocation and residual strain energies to the usual elastic energy. There are two extra material constants that determine the size of the dislocation core and are determined from experimental data. These findings are useful for studies of dislocations in graphene and for analyzing incipient plasticity in two dimensional crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spatially confined Bloch oscillations in semiconductor superlattices

In a semiconductor superlattice with long scattering times, damping of Bloch oscillations due to scattering is so small that convective nonlinearities may compensate it and Bloch oscillations persist even in the hydrodynamic regime. In this case, numerical solutions show that there are stable Bloch oscillations confined to a region near the collector with inhomogeneous field, charge, current density and energy density profiles. These Bloch oscillations disappear when damping due to inelastic collisions becomes sufficiently strong.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statics and dynamics of a harmonic oscillator coupled to a one-dimensional Ising system

We investigate an oscillator linearly coupled with a one-dimensional Ising system. The coupling gives rise to drastic changes both in the oscillator statics and dynamics. Firstly, there appears a second order phase transition, with the oscillator stable rest position as its order parameter. Secondly, for fast spins, the oscillator dynamics is described by an effective equation with a nonlinear friction term that drives the oscillator towards the stable equilibrium state.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theory of defect dynamics in graphene: defect groupings and their stability

We use our theory of periodized discrete elasticity to characterize defects in graphene as the cores of dislocations or groups of dislocations. Earlier numerical implementations of the theory predicted some of the simpler defect groupings observed in subsequent Transmission Electron Microscope experiments. Here we derive the more complicated defect groupings of three or four defect pairs from our theory, show that they correspond to the cores of two pairs of dislocation dipoles and ascertain their stability.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theory of spatially inhomogneous Bloch oscillations in semiconductor superlattices

In a semiconductor superlattice with long scattering times, damping of Bloch oscillations due to scattering is so small that nonlinearities may compensate it and Bloch oscillations persist even in the hydrodynamic regime. To demonstrate this, a Boltzmann-Poisson transport model of miniband superlattices with inelastic collisions is proposed and hydrodynamic equations for electron density, electric field and the complex amplitude of the Bloch oscillations are derived by singular perturbation methods. For appropriate parameter ranges, numerical solutions of these equations show stable Bloch oscillations with spatially inhomogeneous field, charge, current density and energy density profiles. These Bloch oscillations disappear as scattering times become sufficiently short. For sufficiently low lattice temperatures, Bloch and Gunn type oscillations mediated by electric field, current and energy domains coexist for a range of voltages. For larger lattice temperatures (300 K), there are only Bloch oscillations with stationary amplitude and electric field profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Wave fronts, pulses and wave trains in photoexcited superlattices behaving as excitable or oscillatory media

Undoped and strongly photoexcited semiconductor superlattices with field-dependent recombination behave as excitable or oscillatory media with spatially discrete nonlinear convection and diffusion. Infinitely long, dc-current-biased superlattices behaving as excitable media exhibit wave fronts with increasing or decreasing profiles, whose velocities can be calculated by means of asymptotic methods. These superlattices can also support pulses of the electric field. Pulses moving downstream with the flux of electrons can be constructed from their component wave fronts, whereas pulses advancing upstream do so slowly and experience saltatory motion: they change slowly in long intervals of time separated by fast transitions during which the pulses jump to the previous superlattice period. Photoexcited superlattices can also behave as oscillatory media and exhibit wave trains.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nonequilibrium dynamics of a fast oscillator coupled to Glauber spins

A fast harmonic oscillator is linearly coupled with a system of Ising spins that are in contact with a thermal bath, and evolve under a slow Glauber dynamics at dimensionless temperature $θ$. The spins have a coupling constant proportional to the oscillator position. The oscillator-spin interaction produces a second order phase transition at $θ=1$ with the oscillator position as its order parameter: the equilibrium position is zero for $θ>1$ and non-zero for $θ< 1$. For $θ<1$, the dynamics of this system is quite different from relaxation to equilibrium. For most initial conditions, the oscillator position performs modulated oscillations about one of the stable equilibrium positions with a long relaxation time. For random initial conditions and a sufficiently large spin system, the unstable zero position of the oscillator is stabilized after a relaxation time proportional to $θ$. If the spin system is smaller, the situation is the same until the oscillator position is close to zero, then it crosses over to a neighborhood of a stable equilibrium position about which keeps oscillating for an exponentially long relaxation time. These results of stochastic simulations are predicted by modulation equations obtained from a multiple scale analysis of macroscopic equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase transitions in a mechanical system coupled to Glauber spins

A harmonic oscillator linearly coupled with a linear chain of Ising spins is investigated. The $N$ spins in the chain interact with their nearest neighbours with a coupling constant proportional to the oscillator position and to $N^{-1/2}$, are in contact with a thermal bath at temperature $T$, and evolve under Glauber dynamics. The oscillator position is a stochastic process due to the oscillator-spin interaction which produces drastic changes in the equilibrium behaviour and the dynamics of the oscillator. Firstly, there is a second order phase transition at a critical temperature $T_c$ whose order parameter is the oscillator stable rest position: this position is zero above $T_c$ and different from zero below $T_c$. This transition appears because the oscillator moves in an effective potential equal to the harmonic term plus the free energy of the spin system at fixed oscillator position. Secondly, assuming fast spin relaxation (compared to the oscillator natural period), the oscillator dynamical behaviour is described by an effective equation containing a nonlinear friction term that drives the oscillator towards the stable equilibrium state of the effective potential. The analytical results are compared with numerical simulation throughout the paper.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Two mini-band model for self-sustained oscillations of the current through resonant tunneling semiconductor superlattices

A two miniband model for electron transport in semiconductor superlattices that includes scattering and interminiband tunnelling is proposed. The model is formulated in terms of Wigner functions in a basis spanned by Pauli matrices, includes electron-electron scattering in the Hartree approximation and modified Bhatnagar-Gross-Krook collision tems. For strong applied fields, balance equations for the electric field and the miniband populations are derived using a Chapman-Enskog perturbation technique. These equations are then solved numerically for a dc voltage biased superlattice. Results include self-sustained current oscillations due to repeated nucleation of electric field pulses at the injecting contact region and their motion towards the collector. Numerical reconstruction of the Wigner functions shows that the miniband with higher energy is empty during most of the oscillation period: it becomes populated only when the local electric field (corresponding to the passing pulse) is sufficiently large to trigger resonant tunneling.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Photo-excited semiconductor superlattices as constrained excitable media: Motion of dipole domains and current self-oscillations

A model for charge transport in undoped, photo-excited semiconductor superlattices, which includes the dependence of the electron-hole recombination on the electric field and on the photo-excitation intensity through the field-dependent recombination coefficient, is proposed and analyzed. Under dc voltage bias and high photo-excitation intensities, there appear self-sustained oscillations of the current due to a repeated homogeneous nucleation of a number of charge dipole waves inside the superlattice. In contrast to the case of a constant recombination coefficient, nucleated dipole waves can split for a field-dependent recombination coefficient in two oppositely moving dipoles. The key for understanding these unusual properties is that these superlattices have a unique static electric-field domain. At the same time, their dynamical behavior is akin to the one of an extended excitable system: an appropriate finite disturbance of the unique stable fixed point may cause a large excursion in phase space before returning to the stable state and trigger pulses and wave trains. The voltage bias constraint causes new waves to be nucleated when old ones reach the contact.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Kinetics of helium bubble formation in nuclear materials

The formation and growth of helium bubbles due to self-irradiation in plutonium has been modelled by a discrete kinetic equations for the number densities of bubbles having $k$ atoms. Analysis of these equations shows that the bubble size distribution function can be approximated by a composite of: (i) the solution of partial differential equations describing the continuum limit of the theory but corrected to take into account the effects of discreteness, and (ii) a local expansion about the advancing leading edge of the distribution function in size space. Both approximations contribute to the memory term in a close integrodifferential equation for the monomer concentration of single helium atoms. The present boundary layer theory for discrete equations is compared to the numerical solution of the full kinetic model and to previous approximation of Schaldach and Wolfer involving a truncated system of moment equations.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

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&#39;&#39; 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.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Aging in the Linear Harmonic Oscillator

The low temperature Monte Carlo dynamics of an ensemble of linear harmonic oscillators shows some entropic barriers related to the difficulty of finding the directions in configurational space which decrease the energy. This mechanism is enough to observe some typical non-equilibrium features of glassy systems like activated-type behavior and aging in the correlation function and in the response function. Due to the absence of interactions the model only displays a one-step relaxation process.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Asymptotic analysis of the Gunn effect with realistic boundary conditions

A general asymptotic analysis of the Gunn effect in n-GaAs under general boundary conditions for metal-semiconductor contacts is presented. Depending on the parameter values in the boundary condition of the injecting contact, different types of waves mediate the Gunn effect. The periodic current oscillation typical of the Gunn effect may be caused by moving charge-monopole accumulation or depletion layers, or by low or high-field charge-dipole solitary waves. A new instability caused by multiple shedding of (low field) dipole waves is found. In all cases the shape of the current oscillation is described in detail: we show the direct relationship between its major features (maxima, minima, plateau&#39;s, ...) and several critical currents (which depend on the values of the contact parameters). Our results open the possibility of measuring contact parameters from the analysis of the shape of the current oscillation.

0 citations0 reviewsNo code linkedOpen access