Pattern formation with a conservation law

Pattern formation in systems with a conserved quantity is considered by studying the appropriate amplitude equations. The conservation law leads to a large-scale neutral mode that must be included in the asymptotic analysis for pattern formation near onset. Near a stationary bifurcation, the usual Ginzburg--Landau equation for the amplitude of the pattern is then coupled to an equation for the large-scale mode. These amplitude equations show that for certain parameters all roll-type solutions are unstable. This new instability differs from the Eckhaus instability in that it is amplitude-driven and is supercritical. Beyond the stability boundary, there exist stable stationary solutions in the form of strongly modulated patterns. The envelope of these modulations is calculated in terms of Jacobi elliptic functions and, away from the onset of modulation, is closely approximated by a sech profile. Numerical simulations indicate that as the modulation becomes more pronounced, the envelope broadens. A number of applications are considered, including convection with fixed-flux boundaries and convection in a magnetic field, resulting in new instabilities for these systems.

Exotic dynamics in a firing rate model of neural tissue with threshold accommodation

Many of the equations describing the dynamics of neural systems are written in terms of firing rate functions, which themselves are often taken to be threshold functions of synaptic activity. Dating back to work by Hill in 1936 it has been recognized that more realistic models of neural tissue can be obtained with the introduction of state-dependent dynamic thresholds. In this paper we treat a specific phenomenological model of threshold accommodation that mimics many of the properties originally described by Hill. Importantly we explore the consequences of this dynamic threshold at the tissue level, by modifying a standard neural field model of Wilson-Cowan type. As in the case without threshold accommodation classical Mexican-Hat connectivity is shown to allow for the existence of spatially localized states (bumps) in both one and two dimensions. Importantly an analysis of bump stability in one dimension, using recent Evans function techniques, shows that bumps may undergo instabilities leading to the emergence of both breathers and traveling waves. Moreover, a similar analysis for traveling pulses leads to the conditions necessary to observe a stable traveling breather. In the regime where a bump solution does not exist direct numerical simulations show the possibility of self-replicating bumps via a form of bump splitting. Simulations in two space dimensions show analogous localized and traveling solutions to those seen in one dimension. Indeed dynamical behavior in this neural model appears reminiscent of that seen in other dissipative systems that support localized structures, and in particular those of coupled cubic complex Ginzburg-Landau equations. Further numerical explorations illustrate that the traveling pulses in this model exhibit particle like properties, similar to those of dispersive solitons observed in some three component reaction-diffusion systems. A preliminary account of this work first appeared in S Coombes and M R Owen, Bumps, breathers, and waves in a neural network with spike frequency adaptation, Physical Review Letters 94 (2005), 148102(1-4).

Dissipative solitons in fiber lasers

Dissipative solitons (also known as auto-solitons) are stable, nonlinear, time-or space-localized solitary waves that occur due to the balance between energy excitation and dissipation. We review the theory of dissipative solitons applied to fiber laser systems. The discussion context includes the classical Ginzburg-Landau and Maxwell-Bloch equations and their modifications that allow describing laser-cavity-produced waves. Practical examples of laser systems generating dissipative solitons are discussed.

Inertial mass of a vortex in cuprate superconductors

We present here a calculation of the inertial mass of a moving vortex in cuprate superconductors. This is a poorly known basic quantity of obvious interest in vortex dynamics. The motion of a vortex causes a dipolar density distortion and an associated electric field which is screened. The energy cost of the density distortion as well as the related screened electric field contributes to the vortex mass, which is small because of efficient screening. As a preliminary, we present a discussion and calculation of the vortex mass using a microscopically derivable phase-only action functional for the far region which shows that the contribution from the far region is negligible and that most of it arises from the (small) core region of the vortex. A calculation based on a phenomenological Ginzburg-Landau functional is performed in the core region. Unfortunately such a calculation is unreliable; the reasons for it are discussed. A credible calculation of the vortex mass thus requires a fully microscopic non-coarse-grained theory. This is developed, and results are presented for an s-wave BCS-like gap, with parameters appropriate to the cuprates. The mass, about 0.5m(e) per layer, for a magnetic field along the c axis arises from deformation of quasiparticle states bound in the core and screening effects mentioned above. We discuss earlier results, possible extensions to d-wave symmetry, and observability of effects dependent on the inertial mass. [S0163-1829(97)05534-3].

Theoretical studies on magnetic superconductors

The discovery of magnetic superconductors has posed the problem of the coexistence of two kinds of orders (magnetic and superconducting) in some temperature intervals in these systems. New microscopic mechanisms developed by us to explain the coexistence and reentrant behaviour are reported. The mechanism for antiferromagnetic superconductors which shows enhancement of superconductivity below the magnetic transition is found relevant for rare-earth systems having less than half-filled f-atomic shells. The theory will be compared with the experimental results of SmRh4B4 system. A phenomenological treatment based on a generalized Ginzburg-Landau approach will also be presented to explain the anomalous behaviour of the second critical field in some antiferromagnetic superconductors. These magnetic superconductors provide two kinds of Bose fields, namely, phonons and magnons which interact with each other and also with the conduction electrons. Theoretical studies of the effects of the excitations of these modes on superconducting pairing and magnetic ordering in these systems will be discussed.

Magnetic superconductors: Model theories and experimental properties of rare-earth compounds

Superconducting and magnetically long-range ordered states were believed to be mutually exclusive phenomena. The discovery of rare-earth compounds in recent years, which exhibit both superconductivity and magnetic ordering (ferromagnetic, antiferromagnetic or sinusoidal), has led to considerable theoretical and experimental work on such systems. In the present article, we give a review of various theoretical models and important experimental results. In the theoretical sections, we start with the Abrikosov-Gorkov pair breaking theory for dilute alloys and discuss its improvement in the work of Müller-Hartmann and Zittartz. Then, in the context of magnetic superconductors, various microscopic theories that have been advanced are presented. These predict re-entrant behaviour in some systems (ferromagnetic superconductors) and coexistence regions in others (particularly antiferromagnetic superconductors). Following this, phenomenological generalized Ginzburg-Landau theories for two kinds of orders (superconducting and magnetic) are presented. A section dealing with renormalization group analysis of phase diagrams in magnetic superconductors is given. In experimental sections, the properties of each rare-earth compounds (ternary as well as some tetranery) are reviewed. These involve susceptibility, heat capacity, resistivity, upper critical field, neutron scattering and magnetic resonance measurements. The anomalous behaviour of the upper critical field of antiferromagnetic superconductors near the Néel temperature is discussed both in theory sections and experimental section for various systems.

A theory of the upper critical field in antiferromagnetic superconductors

A generalized Ginzburg-Landau approach is used to study the nonmonotonic temperature dependence of the upper critical field H c 2(T) in antiferromagnetic superconductors RE(Mo)6S8; RE = Dy, Tb, Gd. It is found that electrodynamic effects incorporated through screening and indirect coupling between the staggered magnetization M Q (T) and superconducting order parameter psgr cannot explain the observed nonmonotonicity. This suggests that the direct coupling between the two order parameters should be considered to understand the experimental results, a finding which is consistent with recent microscopic calculations.

Study of magnetic field effect on the specific heat and entropy in DyBa2Cu3O7-x

The change in the specific heat by the application of magnetic field up to 161 for high temperature superconductor system for DyBa2Cu3O7-x by Revaz et al. [23] is examined through the phenomenological Ginzburg-Landau(G-L) theory of anisotropic Type-II superconductors. The observed specific heat anomaly near T-c with magnetic field is explained qualitatively through the expression <Delta C > = (B-a/T-c) t/(1 - t)(alpha Theta(gamma)lambda(2)(m)(0)), which is the anisotropic formulation of the G-L theory in the London limit developed by Kogan and coworkers; relating to the change in specific heat Delta C for the variation of applied magnetic field for different orientations with c-axis. The analysis of this equation explains satisfactorily the specific heat anomaly near T-c and determines the anisotropic ratio gamma as 5.608, which is close to the experimental value 5.3 +/- 0.5given in the paper of Revaz et al. for this system. (C) 2010 Elsevier B.V. All rights reserved.

Cell-dynamical simulation of magnetic hysteresis in the two-dimensional Ising system

We present results from numerical simulations using a ‘‘cell-dynamical system’’ to obtain solutions to the time-dependent Ginzburg-Landau equation for a scalar, two-dimensional (2D), (Φ2)2 model in the presence of a sinusoidal external magnetic field. Our results confirm a recent scaling law proposed by Rao, Krishnamurthy, and Pandit [Phys. Rev. B 42, 856 (1990)], and are also in excellent agreement with recent Monte Carlo simulations of hysteretic behavior of 2D Ising spins by Lo and Pelcovits [Phys. Rev. A 42, 7471 (1990)].

Nonlinear hydrodynamics of a hard-sphere fluid near the glass transition

We conduct a numerical study of the dynamic behavior of a dense hard-sphere fluid by deriving and integrating a set of Langevin equations. The statics of the system is described by a free-energy functional of the Ramakrishnan-Yussouff form. We find that the system exhibits glassy behavior as evidenced through a stretched exponential decay and a two-stage relaxation of the density correlation function. The characteristic times grow with increasing density according to the Vogel-Fulcher law. The wave-number dependence of the kinetics is extensively explored. The connection of our results with experiment, mode-coupling theory, and molecular-dynamics results is discussed.

Shear-induced enhancement of self-diffusion in interacting colloidal suspensions

We set up the generalized Langevin equations describing coupled single-particle and collective motion in a suspension of interacting colloidal particles in a shear how and use these to show that the measured self-diffusion coefficients in these systems should be strongly dependent on shear rate epsilon. Three regimes are found: (i) an initial const+epsilon(.2), followed by (ii) a large regime of epsilon(.1/2) behavior, crossing over to an asymptotic power-law approach (iii) D-o - const x epsilon(.-1/2) to the Stokes-Einstein value D-o. The shear dependence is isotropic up to very large shear rates and increases with the interparticle interaction strength. Our results provide a straightforward explanation of recent experiments and simulations on sheared colloids.

Dynamics of Fractional Brownian Walks

We investigate the dynamics of polymers whose solution configurations are represented by fractional Brownian walks. The calculation of the two dynamical quantities considered here, the longest relaxation time tau(r) and the intrinsic viscosity [eta], is formulated in terms of Langevin equations and is carried out within the continuum approach developed in an earlier paper. Our results for tau(r) and [eta] reproduce known scaling relations and provide reasonable numerical estimates of scaling amplitudes. The possible relevance of the work to the study of globular proteins and other compact polymeric phases is discussed.

Vortex dynamics and upper critical fields in ultrathin Bi films

Current-voltage (I-V) characteristics of quench condensed, superconducting, ultrathin Bi films in a magnetic field are reported. These I-V's show hysteresis for all films, grown both with and without thin Ge underlayers. Films on Ge underlayers, close to superconductor-insulator transition, show a peak in the critical current, indicating a structural transformation of the vortex solid. These underlayers, used to make the films more homogeneous, are found to be more effective in pinning the vortices. The upper critical fields (B-c2) of these films are determined from the resistive transitions in perpendicular magnetic field. The temperature dependence of the upper critical field is found to differ significantly from Ginzburg-Landau theory, after modifications for disorder.

Classical integrability in the BTZ black hole

Using the fact the BTZ black hole is a quotient of AdS(3) we show that classical string propagation in the BTZ background is integrable. We construct the flat connection and its monodromy matrix which generates the non-local charges. From examining the general behaviour of the eigen values of the monodromy matrix we determine the set of integral equations which constrain them. These equations imply that each classical solution is characterized by a density function in the complex plane. For classical solutions which correspond to geodesics and winding strings we solve for the eigen values of the monodromy matrix explicitly and show that geodesics correspond to zero density in the complex plane. We solve the integral equations for BMN and magnon like solutions and obtain their dispersion relation. We show that the set of integral equations which constrain the eigen values of the monodromy matrix can be identified with the continuum limit of the Bethe equations of a twisted SL(2, R) spin chain at one loop. The Landau-Lifshitz equations from the spin chain can also be identified with the sigma model equations of motion.

Physically consistent simulation of mesoscale chemical kinetics: The non-negative FIS-alpha method

Biochemical pathways involving chemical kinetics in medium concentrations (i.e., at mesoscale) of the reacting molecules can be approximated as chemical Langevin equations (CLE) systems. We address the physically consistent non-negative simulation of the CLE sample paths as well as the issue of non-Lipschitz diffusion coefficients when a species approaches depletion and any stiffness due to faster reactions. The non-negative Fully Implicit Stochastic alpha (FIS alpha) method in which stopped reaction channels due to depleted reactants are deleted until a reactant concentration rises again, for non-negativity preservation and in which a positive definite Jacobian is maintained to deal with possible stiffness, is proposed and analysed. The method is illustrated with the computation of active Protein Kinase C response in the Protein Kinase C pathway. (C) 2011 Elsevier Inc. All rights reserved.