Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
36works
0followers
23topics
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

36 published item(s)

preprint2016arXiv

Stellar Dynamics around a Massive Black Hole I: Secular Collisionless Theory

We present a theory in three parts, of the secular dynamics of a (Keplerian) stellar system of mass $M$ orbiting a black hole of mass $M_\bullet \gg M$. Here we describe the collisionless dynamics; Papers II and III are on the (collisional) theory of Resonant Relaxation. The mass ratio, $\varepsilon = M/M_\bullet \ll 1$, is a natural small parameter implying a separation of time scales between the short Kepler orbital periods and the longer orbital precessional periods. The collisionless Boltzmann equation (CBE) for the stellar distribution function (DF) is averaged over the fast Kepler orbital phase using the method of multiple scales. The orbit-averaged system is described by a secular DF, $F$, in a reduced phase space. $F$ obeys a secular CBE that includes stellar self-gravity, general relativistic corrections up to 1.5 post-Newtonian order, and external sources varying over secular times. Secular dynamics, even with general time dependence, conserves the semi-major axis of every star. This additional integral of motion promotes extra regularity of the stellar orbits, and enables the construction of equilibria, $F_0$, through a secular Jeans theorem. A linearized secular CBE determines the response and stability of $F_0$. Spherical, non-rotating equilibria may support long-lived, warp-like distortions. We also prove that an axisymmetric, zero-thickness, flat disc is secularly stable to all in-plane perturbations, when its DF, $F_0$, is a monotonic function of the angular momentum at fixed energy.

preprint2016arXiv

Stellar Dynamics around a Massive Black Hole II: Resonant Relaxation

We present a first-principles theory of Resonant Relaxation (RR) of a low mass stellar system orbiting a more massive black hole (MBH). We first extend the kinetic theory of Gilbert (1968) to include the Keplerian field of a black hole of mass $M_\bullet$. Specializing to a Keplerian stellar system of mass $M \ll M_\bullet$, we use the orbit-averaging method of Sridhar & Touma (2015; Paper I) to derive a kinetic equation for RR. This describes the collisional evolution of a system of $N \gg 1$ Gaussian Rings in a reduced 5-dim space, under the combined actions of self-gravity, 1 PN and 1.5 PN relativistic effects of the MBH and an arbitrary external potential. In general geometries RR is driven by both apsidal and nodal resonances, so the distinction between scalar-RR and vector-RR disappears. The system passes through a sequence of quasi-steady secular collisionless equilibria, driven by irreversible 2-Ring correlations that accrue through gravitational interactions, both direct and collective. This correlation function is related to a `wake function', which is the linear response of the system to the perturbation of a chosen Ring. The wake function is easier to appreciate, and satisfies a simpler equation, than the correlation function. We discuss general implications for the interplay of secular dynamics and non-equilibrium statistical mechanics in the evolution of Keplerian stellar systems toward secular thermodynamic equilibria, and set the stage for applications to the RR of axisymmetric discs in Paper III.

preprint2016arXiv

Using skewness and the first-digit phenomenon to identify dynamical transitions in cardiac models

Disruptions in the normal rhythmic functioning of the heart, termed as arrhythmia, often result from qualitative changes in the excitation dynamics of the organ. The transitions between different types of arrhythmia are accompanied by alterations in the spatiotemporal pattern of electrical activity that can be measured by observing the time-intervals between successive excitations of different regions of the cardiac tissue. Using biophysically detailed models of cardiac activity we show that the distribution of these time-intervals exhibit a systematic change in their skewness during such dynamical transitions. Further, the leading digits of the normalized intervals appear to fit Benford's law better at these transition points. This raises the possibility of using these observations to design a clinical indicator for identifying changes in the nature of arrhythmia. More importantly, our results reveal an intriguing relation between the changing skewness of a distribution and its agreement with Benford's law, both of which have been independently proposed earlier as indicators of regime shift in dynamical systems.

preprint2016arXiv

When the clock strikes: Modeling the relation between circadian rhythms and cardiac arrhythmias

It has recently been observed that the occurrence of sudden cardiac death has a close statistical relationship with the time of day, viz., ventricular fibrillation is most likely to occur between 12 am-6 am, with 6 pm-12 am being the next most likely period. Consequently there has been significant interest in understanding how cardiac activity is influenced by the circadian clock, i.e., temporal oscillations in physiological activity with a period close to 24 hours and synchronized with the day-night cycle. Although studies have identified the genetic basis of circadian rhythms at the intracellular level, the mechanisms by which they influence cardiac pathologies are not yet fully understood. Evidence has suggested that diurnal variations in the conductance properties of ion channel proteins that govern the excitation dynamics of cardiac cells may provide the crucial link. In this paper, we investigate the relationship between the circadian rhythm as manifested in modulations of ion channel properties and the susceptibility to cardiac arrhythmias by using a mathematical model that describes electrical activity in ventricular tissue. We show that changes in the channel conductance that lead to extreme values for the duration of action potentials in cardiac cells can result either in abnormally high-frequency reentrant activity or spontaneous conduction block of excitation waves. Both phenomena increase the likelihood of wavebreaks that are known to initiate potentially life-threatening arrhythmias. Thus, disruptive cardiac excitation dynamics are most likely to occur in time-intervals of the day-night cycle during which the channel properties are closest to these extreme values, providing an intriguing relation between circadian rhythms and cardiac pathologies.

preprint2015arXiv

Beyond Thomas--Fermi analysis of the density profiles of a miscible two-component Bose--Einstein condensate

We investigate a harmonically trapped two-component Bose--Einstein condensate within the miscible regime, close to its boundaries, for different ratios of effective intra- and inter-species interactions. We derive analytically a universal equation for the density around the different boundaries in one, two and three dimensions, for both the coexisting and spatially separated regimes. We also present a general procedure to solve the Thomas--Fermi approximation in all three spatial dimensionalities, reducing the complexity of the Thomas--Fermi problem for the spatially separated case in one and three dimensions to a single numerical inversion. Finally, we analytically determine the frontier between the two different regimes of the system.

preprint2015arXiv

Designer dynamics through chaotic traps: Controlling complex behavior in driven nonlinear systems

Control schemes for dynamical systems typically involve stabilizing unstable periodic orbits. In this paper we introduce a new paradigm of control that involves `trapping' the dynamics arbitrarily close to any desired trajectory. This is achieved by a state-dependent dynamical selection of the input signal applied to the driven nonlinear system. An emergent property of the trapping process is that the signal changes in a chaotic sequence: a manifestation of chaos-induced order. The simplicity of the control scheme makes it easily implementable in experimental systems.

preprint2015arXiv

Migration into a Companion's Trap: Disruption of Multiplanet Systems in Binaries

Most exoplanetary systems in binary stars are of S--type, and consist of one or more planets orbiting a primary star with a wide binary stellar companion. Gravitational forcing of a single planet by a sufficiently inclined binary orbit can induce large amplitude oscillations of the planet's eccentricity and inclination through the Kozai-Lidov (KL) instability. KL cycling was invoked to explain: the large eccentricities of planetary orbits; the family of close--in hot Jupiters; and the retrograde planetary orbits in eccentric binary systems. However, several kinds of perturbations can quench the KL instability, by inducing fast periapse precessions which stabilize circular orbits of all inclinations: these could be a Jupiter--mass planet, a massive remnant disc or general relativistic precession. Indeed, mutual gravitational perturbations in multiplanet S--type systems can be strong enough to lend a certain dynamical rigidity to their orbital planes. Here we present a new and faster process that is driven by this very agent inhibiting KL cycling. Planetary perturbations enable secular oscillations of planetary eccentricities and inclinations, also called Laplace--Lagrange (LL) eigenmodes. Interactions with a remnant disc of planetesimals can make planets migrate, causing a drift of LL mode periods which can bring one or more LL modes into resonance with binary orbital motion. The results can be dramatic, ranging from excitation of large eccentricities and mutual inclinations to total disruption. Not requiring special physical or initial conditions, binary resonant driving is generic and could have profoundly altered the architecture of many S--type multiplanet systems. It can also weaken the multiplanet occurrence rate in wide binaries, and affect planet formation in close binaries.

preprint2014arXiv

Large-scale dynamo action due to $α$ fluctuations in a linear shear flow

We present a model of large-scale dynamo action in a shear flow that has stochastic, zero-mean fluctuations of the $α$ parameter. This is based on a minimal extension of the Kraichnan-Moffatt model, to include a background linear shear and Galilean-invariant $α$-statistics. Using the first order smoothing approximation we derive a linear integro-differential equation for the large-scale magnetic field, which is non perturbative in the shearing rate $S\,$, and the $α$-correlation time $τ_α\,$. The white-noise case, $τ_α= 0\,$, is solved exactly, and it is concluded that the necessary condition for dynamo action is identical to the Kraichnan-Moffatt model without shear; this is because white-noise does not allow for memory effects, whereas shear needs time to act. To explore memory effects we reduce the integro-differential equation to a partial differential equation, valid for slowly varying fields when $τ_α$ is small but non zero. Seeking exponential modal solutions, we solve the modal dispersion relation and obtain an explicit expression for the growth rate as a function of the six independent parameters of the problem. A non zero $τ_α$ gives rise to new physical scales, and dynamo action is completely different from the white-noise case; e.g. even weak $α$ fluctuations can give rise to a dynamo. We argue that, at any wavenumber, both Moffatt drift and Shear always contribute to increasing the growth rate. Two examples are presented: (a) a Moffatt drift dynamo in the absence of shear; (b) a Shear dynamo in the absence of Moffatt drift.

preprint2012arXiv

Counter-rotating stellar discs around a massive black hole: self-consistent, time-dependent dynamics

We formulate the collisionless Boltzmann equation (CBE) for dense star clusters that lie within the radius of influence of a massive black hole in galactic nuclei. Our approach to these nearly Keplerian systems follows that of Sridhar and Touma (1999): Delaunay canonical variables are used to describe stellar orbits and we average over the fast Keplerian orbital phases. The stellar distribution function (DF) evolves on the longer time scale of precessional motions, whose dynamics is governed by a Hamiltonian, given by the orbit-averaged self-gravitational potential of the cluster. We specialize to razor-thin, planar discs and consider two counter-rotating ("$\pm$") populations of stars. To describe discs of small eccentricities, we expand the $\pm$ Hamiltonian to fourth order in the eccentricities, with coefficients that depend self-consistently on the $\pm$ DFs. We construct approximate $\pm$ dynamical invariants and use Jeans' theorem to construct time-dependent $\pm$ DFs, which are completely described by their centroid coordinates and shape matrices. When the centroid eccentricities are larger than the dispersion in eccentricities, the $\pm$ centroids obey a set of 4 autonomous ordinary differential equations. We show that these can be cast as a two-degree of freedom Hamiltonian system which is nonlinear, yet integrable. We study the linear instability of initially circular discs and derive a criterion for the counter-rotating instability. We then explore the rich nonlinear dynamics of counter-rotating discs, with focus on the variety of steadily precessing eccentric configurations that are allowed. The stability and properties of these configurations are studied as functions of parameters such as the disc mass ratios and angular momentum.

preprint2012arXiv

Estimation of Effort in Software Cost Analysis for Heterogenous Dataset using Fuzzy Analogy

One of the significant objectives of software engineering community is to use effective and useful models for precise calculation of effort in software cost estimation. The existing techniques cannot handle the dataset having categorical variables efficiently including the commonly used analogy method. Also, the project attributes of cost estimation are measured in terms of linguistic values whose imprecision leads to confusion and ambiguity while explaining the process. There are no definite set of models which can efficiently handle the dataset having categorical variables and endure the major hindrances such as imprecision and uncertainty without taking the classical intervals and numeric value approaches. In this paper, a new approach based on fuzzy logic, linguistic quantifiers and analogy based reasoning is proposed to enhance the performance of the effort estimation in software projects dealing with numerical and categorical data. The performance of this proposed method illustrates that there is a realistic validation of the results while using historical heterogeneous dataset. The results were analyzed using the Mean Magnitude Relative Error (MMRE) and indicates that the proposed method can produce more explicable results than the methods which are in vogue.

preprint2012arXiv

Light beams with general direction and polarization: global description and geometric phase

We construct the manifold describing the family of plane monochromatic light waves with all directions, polarizations, phases and intensities. A smooth description of polarization, valid over the entire sphere S^2 of directions, is given through the construction of an orthogonal basis pair of complex polarization vectors for each direction; any light beam is then uniquely and smoothly specified by giving its direction and two complex amplitudes. This implies that the space of all light beams is the six dimensional manifold S^2 X C^2, the Cartesian product of a sphere and a two dimensional complex vector space. A Hopf map (i.e mapping the two complex amplitudes to the Stokes parameters) then leads to the four dimensional manifold S^2 X S^2 which describes beams with all directions and polarization states. This product of two spheres can be viewed as an ordered pair of two points on a single sphere, in contrast to earlier work in which the same system was represented using Majorana's mapping of the states of a spin one quantum system to an unordered pair of points on a sphere. This is a different manifold, CP^2, two dimensional complex projective space, which does not faithfully represent the full space of all directions and polarizations. Following the now-standard framework, we exhibit the fibre bundle whose total space is the set of all light beams of non-zero intensity, and base space S^2 X S^2. We give the U(1) connection which determines the geometric phase as the line integral of a one-form along a closed curve in the total space. Bases are classified as globally smooth, global but singular, and local, with the last type of basis being defined only when the curve traversed by the system is given. Existing as well as new formulae for the geometric phase are presented in this overall framework.

preprint2012arXiv

Mean field dynamo action in renovating shearing flows

We study mean field dynamo action in renovating flows with finite and non zero correlation time ($τ$) in the presence of shear. Previous results obtained when shear was absent are generalized to the case with shear. The question of whether the mean magnetic field can grow in the presence of shear and non helical turbulence, as seen in numerical simulations, is examined. We show in a general manner that, if the motions are strictly non helical, then such mean field dynamo action is not possible. This result is not limited to low (fluid or magnetic) Reynolds numbers nor does it use any closure approximation; it only assumes that the flow renovates itself after each time interval $τ$. Specifying to a particular form of the renovating flow with helicity, we recover the standard dispersion relation of the $α^2 Ω$ dynamo, in the small $τ$ or large wavelength limit. Thus mean fields grow even in the presence of rapidly growing fluctuations, surprisingly, in a manner predicted by the standard quasilinear closure, even though such a closure is not strictly justified. Our work also suggests the possibility of obtaining mean field dynamo growth in the presence of helicity fluctuations, although having a coherent helicity will be more efficient.

preprint2012arXiv

Optimization of fuzzy analogy in software cost estimation using linguistic variables

One of the most important objectives of software engineering community has been the increase of useful models that beneficially explain the development of life cycle and precisely calculate the effort of software cost estimation. In analogy concept, there is deficiency in handling the datasets containing categorical variables though there are innumerable methods to estimate the cost. Due to the nature of software engineering domain, generally project attributes are often measured in terms of linguistic values such as very low, low, high and very high. The imprecise nature of such value represents the uncertainty and vagueness in their elucidation. However, there is no efficient method that can directly deal with the categorical variables and tolerate such imprecision and uncertainty without taking the classical intervals and numeric value approaches. In this paper, a new approach for optimization based on fuzzy logic, linguistic quantifiers and analogy based reasoning is proposed to improve the performance of the effort in software project when they are described in either numerical or categorical data. The performance of this proposed method exemplifies a pragmatic validation based on the historical NASA dataset. The results were analyzed using the prediction criterion and indicates that the proposed method can produce more explainable results than other machine learning methods.

preprint2012arXiv

Unstable $m=1$ modes of counter-rotating Keplerian discs

We study the linear $m=1$ counter-rotating instability in a two-component, nearly Keplerian disc. Our goal is to understand these \emph{slow} modes in discs orbiting massive black holes in galactic nuclei. They are of interest not only because they are of large spatial scale--and can hence dominate observations--but also because they can be growing modes that are readily excited by accretion events. Self-gravity being nonlocal, the eigenvalue problem results in a pair of coupled integral equations, which we derive for a two-component softened gravity disc. We solve this integral eigenvalue problem numerically for various values of mass fraction in the counter-rotating component. The eigenvalues are in general complex, being real only in the absence of the counter-rotating component, or imaginary when both components have identical surface density profiles. Our main results are as follows: (i) the pattern speed appears to be non negative, with the growth (or damping) rate being larger for larger values of the pattern speed; (ii) for a given value of the pattern speed, the growth (or damping) rate increases as the mass in the counter-rotating component increases; (iii) the number of nodes of the eigenfunctions decreases with increasing pattern speed and growth rate. Observations of lopsided brightness distributions would then be dominated by modes with the least number of nodes, which also possess the largest pattern speeds and growth rates.

preprint2011arXiv

A Classical Fuzzy Approach for Software Effort Estimation on Machine Learning Technique

Software Cost Estimation with resounding reliability,productivity and development effort is a challenging and onerous task. This has incited the software community to give much needed thrust and delve into extensive research in software effort estimation for evolving sophisticated methods. Estimation by analogy is one of the expedient techniques in software effort estimation field. However, the methodology utilized for the estimation of software effort by analogy is not able to handle the categorical data in an explicit and precise manner. A new approach has been developed in this paper to estimate software effort for projects represented by categorical or numerical data using reasoning by analogy and fuzzy approach. The existing historical data sets, analyzed with fuzzy logic, produce accurate results in comparison to the data set analyzed with the earlier methodologies.

preprint2011arXiv

Critical role of pinning defects in scroll-wave breakup in active media

The breakup of rotating scroll waves in three-dimensional excitable media has been linked to important biological processes. The known mechanisms for this transition almost exclusively involve the dynamics of the scroll filament, i.e., the line connecting the phase singularities. In this paper, we describe a novel defect-induced route to breakup of a scroll wave pinned by an inexcitable obstacle partially extending through the bulk of the medium. The wave is helically wound around the defect inducing sudden changes in velocity components of the wavefront at the obstacle boundary. This results in breakup far from the filament, eventually giving rise to spatiotemporal chaos. Our results suggest a potentially critical role of pinning obstacles in the onset of life-threatening disturbances of cardiac activity.

preprint2011arXiv

Robust method to determine the resolution of a superlens by analyzing the near-field image of a two-slit object

In the last decade, metamaterials-based superlenses, with a resolution below Abbe's diffraction limit, have emerged. To obtain a rough estimate of the resolution of such superlenses, imaging of two subwavelength slits, separated by a subwavelength gap \textit{d} is typically performed. The resolution $Δ$ of the lens corresponds to the minimum possible gap $d_{min}$ for which a distinct image of the two slits can be resolved ($Δ\sim d_{min}$). In this letter, we present a more quantitative estimate of the resolution of manufactured lenses by fitting analytical near-field image profiles, obtained from imaging a two-slit object with a theoretical negative-index lens of known resolution, to experimental data. We conclude the discussion by applying our analytical method to 3 case examples of superlensing from the literature. As shown, this method is particularly attractive for rapidly assessing the performance of fabricated superresolution lenses.

preprint2011arXiv

Transport coefficients for the shear dynamo problem at small Reynolds numbers

We build on the formulation developed in Sridhar & Singh (JFM, 664, 265, 2010), and present a theory of the \emph{shear dynamo problem} for small magnetic and fluid Reynolds numbers, but for arbitrary values of the shear parameter. Specializing to the case of a mean magnetic field that is slowly varying in time, explicit expressions for the transport coefficients, $α_{il}$ and $η_{iml}$, are derived. We prove that, when the velocity field is non helical, the transport coefficient $α_{il}$ vanishes. We then consider forced, stochastic dynamics for the incompressible velocity field at low Reynolds number. An exact, explicit solution for the velocity field is derived, and the velocity spectrum tensor is calculated in terms of the Galilean--invariant forcing statistics. We consider forcing statistics that is non helical, isotropic and delta-correlated-in-time, and specialize to the case when the mean-field is a function only of the spatial coordinate $X_3$ and time $τ\,$; this reduction is necessary for comparison with the numerical experiments of Brandenburg, R{ä}dler, Rheinhardt & Käpylä (ApJ, 676, 740, 2008). Explicit expressions are derived for all four components of the magnetic diffusivity tensor, $η_{ij}(τ)\,$. These are used to prove that the shear-current effect cannot be responsible for dynamo action at small $\re$ and $\rem$, but for all values of the shear parameter.

preprint2010arXiv

Nonperturbative quasilinear approach to the shear dynamo problem

We study large-scale dynamo action due to turbulence in the presence of a linear shear flow. Our treatment is quasilinear and equivalent to the standard `first order smoothing approximation'. However it is non perturbative in the shear strength. We first derive an integro-differential equation for the evolution of the mean magnetic field, by systematic use of the shearing coordinate transformation and the Galilean invariance of the linear shear flow. We show that, for non helical turbulence, the time evolution of the cross-shear components of the mean field do not depend on any other components excepting themselves; this is valid for any Galilean-invariant velocity field, independent of its dynamics. Hence, to all orders in the shear parameter, there is no shear-current type effect for non helical turbulence in a linear shear flow, in quasilinear theory in the limit of zero resistivity. We then develop a systematic approximation of the integro-differential equation for the case when the mean magnetic field varies slowly compared to the turbulence correlation time. For non-helical turbulence, the resulting partial differential equations can again be solved by making a shearing coordinate transformation in Fourier space. The resulting solutions are in the form of shearing waves, labeled by the wavenumber in the sheared coordinates. These shearing waves can grow at early and intermediate times but are expected to decay in the long time limit.

preprint2010arXiv

Response to sub-threshold stimulus is enhanced by spatially heterogeneous activity

Sub-threshold stimuli cannot initiate excitations in active media, but surprisingly as we show in this paper, they can alter the time-evolution of spatially heterogeneous activity by modifying the recovery dynamics. This results in significant reduction of waveback velocity which may lead to spatial coherence, terminating all activity in the medium including spatiotemporal chaos. We analytically derive model-independent conditions for which such behavior can be observed.

preprint2010arXiv

Slow light, open cavity formation, and large longitudinal electric field on slab waveguide made of indefinite permittivity materials

The optical properties of slab waveguides made of indefinite permittivity ($\vep$) materials (IEM) are considered. In this medium the transverse permittivity is negative while the longitudinal permittivity is positive. At any given frequency the waveguide supports an infinite number of transverse magnetic (TM) eigenmodes. For a slab waveguide with a fixed thickness, at most only one TM mode is forward-wave. The rest of them are backward waves which can have very large phase index. At a critical thickness, the waveguide supports degenerate forward- and backward-wave modes with zero group velocity. Above the critical thickness, the waveguide supports complex-conjugate decay modes instead of propagating modes. The presence of loss in IEMs will lift the degeneracy, resulting in modes with finite group velocity. Feasible realization is proposed. The performance of IEM waveguide is analyzed and possible applications are discussed which are supported by numerical calculations. These slab waveguides can be used to make optical delay lines in optical buffers to slow down and trap light, to form open cavities, to generate strong longitudinal electric fields, and as phase shifters in optical integrated circuits.

preprint2010arXiv

The shear dynamo problem for small magnetic Reynolds numbers

We study large-scale dynamo action due to turbulence in the presence of a linear shear flow, in the low conductivity limit. Our treatment is nonperturbative in the shear strength and makes systematic use of both the shearing coordinate transformation and the Galilean invariance of the linear shear flow. The velocity fluctuations are assumed to have low magnetic Reynolds number (Rm) but could have arbitrary fluid Reynolds number. The magnetic fluctuations are determined to lowest order in Rm by explicit calculation of the resistive Green's function for the linear shear flow. The mean electromotive force is calculated and an integro-differential equation is derived for the time evolution of the mean magnetic field. In this equation, velocity fluctuations contribute to two different kinds of terms, the C and D terms, in which first and second spatial derivatives of the mean magnetic field, respectively, appear inside the spacetime integrals. The contribution of the D terms is such that the time evolution of the cross-shear components of the mean field do not depend on any other components excepting themselves. Therefore, to lowest order in Rm but to all orders in the shear strength, the D terms cannot give rise to a shear-current assisted dynamo effect. Casting the integro-differential equation in Fourier space, we show that the normal modes of the theory are a set of shearing waves, labelled by their sheared wavevectors. The integral kernels are expressed in terms of the velocity spectrum tensor, which is the fundamental quantity that needs to be specified to complete the integro-differential equation description of the time evolution of the mean magnetic field.

preprint2009arXiv

Anomalous drift of spiral waves in heterogeneous excitable media

We study the drift of spiral waves in a simple model of heterogeneous excitable medium, having gradients in local excitability or cellular coupling. For the first time, we report the anomalous drift of spiral waves towards regions having higher excitability, in contrast to all earlier observations in reaction-diffusion models of excitable media. Such anomalous drift can promote the onset of complex spatio-temporal patterns, e.g., those responsible for life-threatening arrhythmias in the heart.

preprint2009arXiv

Magnetohydrodynamic turbulence in a strongly magnetised plasma

I present a review of incompressible magnetohydrodynamic (MHD) turbulence in a strongly magnetised plasma. The approach is phenomenological even where a more rigorous theory is available, so that a reader armed with paper, pencil and some determination may be able to work through most of the physics. The focus is on the inertial-range spectra for very large (fluid and magnetic) Reynolds numbers. These theories of the inertial-range are built on two important facts: (i) Kraichnan's insight that the turbulent cascades are a result of nonlinear interactions between oppositely directed wavepackets of Elsasser fields; (ii) these oppositely directed wavepackets do not exchange energy, but contribute only to changing each other's spatial structures. I begin with a description and critique of the Iroshnikov-Kraichnan theory, and explore the fundamental departures necessitated by the anisotropic nature of the turbulence. Derivations of the inertial-range spectra of four regimes of MHD turbulence -- the balanced weak, balanced strong, imbalanced weak and the imbalanced strong cascades -- are then presented. The need for studying the spectra of imbalanced turbulence when the waves on the outer scale have a short correlation time is briefly discussed.

preprint2009arXiv

Shear dynamo problem: Quasilinear kinematic theory

Large-scale dynamo action due to turbulence in the presence of a linear shear flow is studied. Our treatment is quasilinear and kinematic but is non perturbative in the shear strength. We derive the integro-differential equation for the evolution of the mean magnetic field, by systematic use of the shearing coordinate transformation and the Galilean invariance of the linear shear flow. For non helical turbulence the time evolution of the cross-shear components of the mean field do not depend on any other components excepting themselves. This is valid for any Galilean-invariant velocity field, independent of its dynamics. Hence the shear-current assisted dynamo is essentially absent, although large-scale non helical dynamo action is not ruled out.

preprint2009arXiv

Slow m=1 instabilities of softened gravity Keplerian discs

We present the simplest model that permits a largely analytical exploration of the m=1 counter-rotating instability in a "hot" nearly Keplerian disc of collisionless self-gravitating matter. The model consists of a two-component softened gravity disc, whose linear modes are analysed using WKB. The modes are slow in the sense that their (complex) frequency is smaller than the Keplerian orbital frequency by a factor which is of order the ratio of the disc mass to the mass of the central object. Very simple analytical expressions are derived for the precession frequencies and growth rates of local modes; it is shown that a nearly Keplerian disc must be unrealistically hot to avoid an overstability. Global modes are constructed for the case of zero net rotation.

preprint2007arXiv

Alternative approach to all-angle negative refraction in two-dimensional photonic crystals

We show that with an appropriate surface modification, a slab of photonic crystal can be made to allow wave transmission within the band gap. Furthermore, negative refraction and all-angle-negative-refraction (AANR) can be achieved by this surface modification in frequency windows that were not realized before in two-dimensional photonic crystals [C. Luo et al, Phys. Rev. B 65, 201104 (2002)]. This approach to AANR leads to new applications in flat lens imaging. Previous flat lens using photonic crystals requires object-image distance u+v less than or equal to the lens thickness d, u+v d. Our approach can be used to design flat lens with u+v=sd with s>>1, thus being able to image large and/or far away objects. Our results are confirmed by FDTD simulations.

preprint2007arXiv

Controlling spatiotemporal chaos and spiral turbulence in excitable media: A review

Excitable media are a generic class of models used to simulate a wide variety of natural systems including cardiac tissue. Propagation of excitation waves in this medium results in the formation of characteristic patterns such as rotating spiral waves. Instabilities in these structures may lead to spatiotemporal chaos through spiral turbulence, which has been linked to clinically diagnosed conditions such as cardiac fibrillation. Usual methods for controlling such phenomena involve very large amplitude perturbations and have several drawbacks. There have been several recent attempts to develop low-amplitude control procedures for spatiotemporal chaos in excitable media which are reviewed in this paper. The control schemes have been broadly classified by us into three types: (i) global, (ii) non-global spatially-extended and (iii) local, depending on the way the control signal is applied, and we discuss the merits and drawbacks for each.

preprint1997arXiv

Microwave properties of $(Pr_xY_{1-x})Ba_2Cu_3O_{7-δ}$ : Influence of magnetic scattering

We report measurements of the surface impedance $Z_s=R_s+iX_s$ of $(Pr_xY_{1-x})Ba_2Cu_3O_{7-δ}$, $(x=0,0.15,0.23,0.3,0.4,0.5)$. Increasing $Pr$ concentration leads to some striking results not observed in samples doped by non-magnetic constituents. The three principal features of the $R_s(T)$ data - multiple structure in the transition, a high residual resistance and, at high $Pr$ concentrations, an upturn of the low $T$ data, are all characteristic of the influence of magnetic scattering on superconductivity, and appear to be common to materials where magnetism and superconductivity coexist. The low $T$ behavior of $λ(T)$ appears to change from $T$ to $T^4$ at large $Pr$ doping, and provides evidence of the influence of magnetic pairbreaking of the $Pr$.

preprint1997arXiv

New features in the microwave response of YBCO crystals: Evidence for a multi-component order parameter

New features are reported in precision measurements of the complex microwave conductivity of high quality $YBa_{2}Cu_{3}O_{6.95}$ crystals grown in $BaZrO_{3}$ crucibles. A third peak in the normal conductivity, $σ_{1}(T)$, at around $80K$, and enhanced pair conductivity $σ_{2}(T)$ below $\sim 65K$ are observed. The data are inconsistent with a single order parameter, and instead are indicative of multi-component superconductivity. Overall, these results point to the presence of multiple pairing interactions in $YBa_{2}Cu_{3}O_{6.95}$ and also provide a natural explanation to account for the low temperature $35K$ conductivity peak observed in all $YBa_{2}Cu_{3}O_{6.95}$ crystals.

preprint1997arXiv

Three dimensional, axisymmetric cusps without chaos

We construct three dimensional axisymmetric, cuspy density distributions, whose potentials are of Stäckel form in parabolic coordinates. As in Sridhar and Touma (1997), a black hole of arbitrary mass may be added at the centre, without destroying the Stäckel form of the potentials. The construction uses a classic method, originally due to Kuzmin (1956), which is here extended to parabolic coordinates. The models are highly oblate, and the cusps are &#34;weak&#34;, with the density, $ρ\propto 1/r^k$, where $0<k<1$.

preprint1996arXiv

Cusps without chaos

We present cuspy, non-axisymmetric, scale-free mass models of discs, whose gravitational potentials are of Stäckel form in parabolic coordinates. A black hole may be added at the centre, without in any way affecting the Stäckel form; the dynamics in these potentials is, of course, fully integrable. The surface density, $Σ_{disc}\propto 1/r^γ$, where $0 <γ< 1$ corresponds to steep cusps for which the central force diverges. Thus cusps, black holes, and non-axisymmetry are not a sure recipe for chaos, as is generally assumed. A new family of orbits, lens orbits, emerges to replace the box orbits of models of elliptical galaxies that have constant-density cores. Loop orbits are conspicuous by their absence. Both lenses and boxlets (the other family of orbits), can be elongated in the direction of the density distribution, a property that is favourable for the construction of non-axisymmetric, self-consistent equilibrium models of elliptical galaxies.

preprint1996arXiv

MHD Turbulence Revisited

Kraichnan (1965) proposed that MHD turbulence occurs as a result of collisions between oppositely directed Alfvén wave packets. Recent work has generated some controversy over the nature of non linear couplings between colliding Alfvén waves. We find that the resolution to much of the confusion lies in the existence of a new type of turbulence, intermediate turbulence, in which the cascade of energy in the inertial range exhibits properties intermediate between those of weak and strong turbulent cascades. Some properties of intermediate MHD turbulence are: (i) in common with weak turbulent cascades, wave packets belonging to the inertial range are long lived; (ii) however, components of the strain tensor are so large that, similar to the situation in strong turbulence, perturbation theory is not applicable; (iii) the breakdown of perturbation theory results from the divergence of neighboring field lines due to wave packets whose perturbations in velocity and magnetic fields are localized, but whose perturbations in displacement are not; (iv) 3--wave interactions dominate individual collisions between wave packets, but interactions of all orders $n\geq 3$ make comparable contributions to the intermediate turbulent energy cascade; (v) successive collisions are correlated since wave packets are distorted as they follow diverging field lines; (vi) in common with the weak MHD cascade, there is no parallel cascade of energy, and the cascade to small perpendicular scales strengthens as it reaches higher wave numbers; (vii) For an appropriate weak excitation, there is a natural progression from a weak, through an intermediate, to a strong cascade.

preprint1996arXiv

Vortex Response and Critical Fields observed via rf penetration depth measurements on the superconductor YNi_2B_2C

Measurements of the rf penetration depth λ(T,H,θ) are used to study the superconducting order parameter, vortex dynamics in the mixed state and delineate critical fields in the borocarbide superconductor YNi_2B_2C. The lower critical field has an anomalous T dependence, H_{c1}(T)=1.12[1-(T/T_c)] kOe, which is however consistent with independent superfluid density measurements at microwave frequencies. The vortex response is dominated by viscous flux flow, indicative of extremely weak pinning, and is parametrized by a field scale H_{c2,eff}. The angular dependence of the vortex contribution λ(θ) is in good agreement with the Coffey-Clem model. Structure is seen in the depairing transition in the vicinity of the upper critical field, with the existence of well-defined critical fields H_{c2a}, H_{c2b} and H_{c2c}, with the vortex field scale H_{c2,eff} closest to H_{c2b}. Overall the measurements indicate that YNi_2B_2C has a rich and unusual field dependence of its transport parameters.

preprint1995arXiv

Experimental studies of Chaos and Localization in Quantum Wavefunctions

Wavefunctions in chaotic and disordered quantum billiards are studied experimentally using thin microwave cavities. The chaotic wavefunctions display universal density distributions and density auto-correlations in agreement with expressions derived fro m a 0-D nonlinear $σ$-model of supersymmetry, which coincides with Random Matrix Theory. In contrast, disordered wavefunctions show deviations from this universal behavior due to Anderson localization. A systematic behavior of the distribution function is studied as a function of the localization length, and can be understood in the framework of a 1-D version of the nonlinear $σ$-model.

preprint1995arXiv

Microwave Properties of Borocarbide Superconductors LnNi2B2C (Ln = Y, Er, Tm, Ho)

We report measurements of the microwave surface impedance of the borocarbide family of superconductors LnNi2B2C (Ln=Y, Er, Tm, Ho). The experiments enable direct measurements of the superfluid density, and are particularly sensitive to the influence of magnetic pairbreaking. In HoNi2B2C the antiferromagnetic transition is clearly observed at zero field, and leads to a drastic reduction of the superfluid density, which recovers at lower temperatures. In ErNi2B2C the antiferromagnetic transition is not seen in zero field data. Magnetic effects are responsible for anomalies in the low temperature surface impedance below approximately 4K in HoNi2B2C and TmNi2B2C. The temperature dependence of the microwave impedance disagrees with simple BCS calculations.