A perturbed martingale approach to global optimization

A new global stochastic search, guided mainly through derivative-free directional information computable from the sample statistical moments of the design variables within a Monte Carlo setup, is proposed. The search is aided by imparting to the directional update term additional layers of random perturbations referred to as `coalescence' and `scrambling'. A selection step, constituting yet another avenue for random perturbation, completes the global search. The direction-driven nature of the search is manifest in the local extremization and coalescence components, which are posed as martingale problems that yield gain-like update terms upon discretization. As anticipated and numerically demonstrated, to a limited extent, against the problem of parameter recovery given the chaotic response histories of a couple of nonlinear oscillators, the proposed method appears to offer a more rational, more accurate and faster alternative to most available evolutionary schemes, prominently the particle swarm optimization. (C) 2014 Elsevier B.V. All rights reserved.

Salmonella methylglyoxal detoxification by STM3117-encoded lactoylglutathione lyase affects virulence in coordination with Salmonella pathogenicity island 2 and phagosomal acidification

Intracellular pathogens such as Salmonella enterica serovar Typhimurium (S. Typhimurium) manipulate their host cells through the interplay of various virulence factors. A multitude of such virulence factors are encoded on the genome of S. Typhimurium and are usually organized in pathogenicity islands. The virulence-associated genomic stretch of STM3117-3120 has structural features of pathogenicity islands and is present exclusively in non-typhoidal serovars of Salmonella. It encodes metabolic enzymes predicted to be involved in methylglyoxal metabolism. STM3117-encoded lactoylglutathione lyase significantly impacts the proliferation of intracellular Salmonella. The deletion mutant of STM3117 (Delta lgl) fails to grow in epithelial cells but hyper-replicates in macrophages. This difference in proliferation outcome was the consequence of failure to detoxify methylglyoxal by Delta lgl, which was also reflected in the form of oxidative DNA damage and upregulation of kefB in the mutant. Within macrophages, the toxicity of methylglyoxal adducts elicits the potassium efflux channel (KefB) in the mutant which subsequently modulates the acidification of mutant-containing vacuoles (MCVs). The perturbation in the pH of the MCV milieu and bacterial cytosol enhances the Salmonella pathogenicity island 2 translocation in Delta lgl, increasing its net growth within macrophages. In epithelial cells, however, the maturation of Delta lgl-containing vacuoles were affected as these non-phagocytic cells maintain less acidic vacuoles compared to those in macrophages. Remarkably, ectopic expression of Toll-like receptors 2 and 4 on epithelial cells partially restored the survival of Delta lgl. This study identified a novel metabolic enzyme in S. Typhimurium whose activity during intracellular infection within a given host cell type differentially affected the virulence of the bacteria.

Two timescale convergent Q-learning for sleep-scheduling in wireless sensor networks

In this paper, we consider an intrusion detection application for Wireless Sensor Networks. We study the problem of scheduling the sleep times of the individual sensors, where the objective is to maximize the network lifetime while keeping the tracking error to a minimum. We formulate this problem as a partially-observable Markov decision process (POMDP) with continuous stateaction spaces, in a manner similar to Fuemmeler and Veeravalli (IEEE Trans Signal Process 56(5), 2091-2101, 2008). However, unlike their formulation, we consider infinite horizon discounted and average cost objectives as performance criteria. For each criterion, we propose a convergent on-policy Q-learning algorithm that operates on two timescales, while employing function approximation. Feature-based representations and function approximation is necessary to handle the curse of dimensionality associated with the underlying POMDP. Our proposed algorithm incorporates a policy gradient update using a one-simulation simultaneous perturbation stochastic approximation estimate on the faster timescale, while the Q-value parameter (arising from a linear function approximation architecture for the Q-values) is updated in an on-policy temporal difference algorithm-like fashion on the slower timescale. The feature selection scheme employed in each of our algorithms manages the energy and tracking components in a manner that assists the search for the optimal sleep-scheduling policy. For the sake of comparison, in both discounted and average settings, we also develop a function approximation analogue of the Q-learning algorithm. This algorithm, unlike the two-timescale variant, does not possess theoretical convergence guarantees. Finally, we also adapt our algorithms to include a stochastic iterative estimation scheme for the intruder's mobility model and this is useful in settings where the latter is not known. Our simulation results on a synthetic 2-dimensional network setting suggest that our algorithms result in better tracking accuracy at the cost of only a few additional sensors, in comparison to a recent prior work.

Newton-based stochastic optimization using q-Gaussian smoothed functional algorithms

We present the first q-Gaussian smoothed functional (SF) estimator of the Hessian and the first Newton-based stochastic optimization algorithm that estimates both the Hessian and the gradient of the objective function using q-Gaussian perturbations. Our algorithm requires only two system simulations (regardless of the parameter dimension) and estimates both the gradient and the Hessian at each update epoch using these. We also present a proof of convergence of the proposed algorithm. In a related recent work (Ghoshdastidar, Dukkipati, & Bhatnagar, 2014), we presented gradient SF algorithms based on the q-Gaussian perturbations. Our work extends prior work on SF algorithms by generalizing the class of perturbation distributions as most distributions reported in the literature for which SF algorithms are known to work turn out to be special cases of the q-Gaussian distribution. Besides studying the convergence properties of our algorithm analytically, we also show the results of numerical simulations on a model of a queuing network, that illustrate the significance of the proposed method. In particular, we observe that our algorithm performs better in most cases, over a wide range of q-values, in comparison to Newton SF algorithms with the Gaussian and Cauchy perturbations, as well as the gradient q-Gaussian SF algorithms. (C) 2014 Elsevier Ltd. All rights reserved.

Instabilities in a fluid overlying an inclined anisotropic and inhomogeneous porous layer

In this paper, linear stability analysis on a Newtonian fluid film flowing under the effect of gravity over an inclined porous medium saturated with the same fluid in isothermal condition is carried out. The focus is placed on the effect of the anisotropic and inhomogeneous variations in the permeability of the porous medium on the shear mode and surface mode instabilities. The fluid-porous system is modelled by a coupled two-dimensional Navier-Stokes/Darcy problem. The perturbation equations are solved numerically using the Chebyshev collocation method. Detailed stability characteristics as a function of the depth ratio (the ratio of the depth of the fluid layer to that of the porous layer), the anisotropic parameter (the ratio of the permeability in the direction of the basic flow to that in the direction transverse to the basic flow) and the inhomogeneity functions are presented.

Noise Detected NMR Spectroscopy

Spin noise phenomenon was predicted way back in 1946. However, experimental investigations regarding spin noise became possible only recently with major technological improvements in NMR hardware. These experiments have several potential novel applications and also demand refinements in the existing theoretical framework to explain the phenomenon. Elegance of noise spectroscopy in gathering information about the properties of a system lies in the fact that it does not require external perturbation, and the system remains in thermal equilibrium. Spin noise is intrinsic magnetic fluctuations, and both longitudinal and transverse components have been detected independently in many systems. Detection of fluctuating longitudinal magnetization leads to field of Magnetic Resonance Force Microscopy (MRFM) that can efficiently probe very few spins even down to the level of single spin utilizing ultrasensitive cantilevers. Transverse component of spin noise, which can simultaneously monitor different resonances over a given frequency range enabling one to distinguish between different chemical environments, has also received considerable attention, and found many novel applications. These experiments demand a detailed understanding of the underlying spin noise phenomenon in order to perform perturbation-free magnetic resonance and widen the highly promising application area. Detailed investigations of noise magnetization have been performed recently using force microscopy on equilibrium ensemble of paramagnetic alkali atoms. It was observed that random fluctuations generate spontaneous spin coherences which has similar characteristics as generated by macroscopic magnetization of polarized ensemble in terms of precession and relaxation properties. Several other intrinsic properties like g-factors, isotope-abundance ratios, hyperfine splitting, spin coherence lifetimes etc. also have been achieved without having to excite the sample. In contrast to MRFM-approaches, detection of transverse spin noise also offers novel applications, attracting considerable attention. This has unique advantage as different resonances over a given frequency range enable one to distinguish between different chemical environments. Since these noise signatures scale inversely with sample size, these approaches lead to the possibility of non-perturbative magnetic resonance of small systems down to nano-scale. In this review, these different approaches will be highlighted with main emphasis on transverse spin noise investigations.

Adaptive mesh generation for singularly perturbed fourth-order ordinary differential equations

In this paper, we consider a singularly perturbed boundary-value problem for fourth-order ordinary differential equation (ODE) whose highest-order derivative is multiplied by a small perturbation parameter. To solve this ODE, we transform the differential equation into a coupled system of two singularly perturbed ODEs. The classical central difference scheme is used to discretize the system of ODEs on a nonuniform mesh which is generated by equidistribution of a positive monitor function. We have shown that the proposed technique provides first-order accuracy independent of the perturbation parameter. Numerical experiments are provided to validate the theoretical results.

A robust layer adapted difference method for singularly perturbed two-parameter parabolic problems

A finite difference method for a time-dependent singularly perturbed convection-diffusion-reaction problem involving two small parameters in one space dimension is considered. We use the classical implicit Euler method for time discretization and upwind scheme on the Shishkin-Bakhvalov mesh for spatial discretization. The method is analysed for convergence and is shown to be uniform with respect to both the perturbation parameters. The use of the Shishkin-Bakhvalov mesh gives first-order convergence unlike the Shishkin mesh where convergence is deteriorated due to the presence of a logarithmic factor. Numerical results are presented to validate the theoretical estimates obtained.

Three centered hydrogen bonds of the type C=O...H(N)...X-C in diphenyloxamide derivatives involving halogens and a rotating CF3 group: NMR, QTAIM, NCI and NBO studies

The existence of three centered C=O...H(N)...X-C hydrogen bonds (H-bonds) involving organic fluorine and other halogens in diphenyloxamide derivatives has been explored by NMR spectroscopy and quantum theoretical studies. The three centered H-bond with the participation of a rotating CF3 group and the F...H-N intramolecular hydrogen bonds, a rare observation of its kind in organofluorine compounds, has been detected. It is also unambiguously established by a number of one and two dimensional NMR experiments, such as temperature perturbation, solvent titration, N-15-H-1 HSQC, and F-19-H-1 HOESY, and is also confirmed by theoretical calculations, such as quantum theory of atoms in molecules (QTAIM), natural bond orbital (NBO) and non-covalent interaction (NCI).

THREE-LAYER FLUID FLOW OVER A SMALL OBSTRUCTION ON THE BOTTOM OF A CHANNEL

Many boundary value problems occur in a natural way while studying fluid flow problems in a channel. The solutions of two such boundary value problems are obtained and analysed in the context of flow problems involving three layers of fluids of different constant densities in a channel, associated with an impermeable bottom that has a small undulation. The top surface of the channel is either bounded by a rigid lid or free to the atmosphere. The fluid in each layer is assumed to be inviscid and incompressible, and the flow is irrotational and two-dimensional. Only waves that are stationary with respect to the bottom profile are considered in this paper. The effect of surface tension is neglected. In the process of obtaining solutions for both the problems, regular perturbation analysis along with a Fourier transform technique is employed to derive the first-order corrections of some important physical quantities. Two types of bottom topography, such as concave and convex, are considered to derive the profiles of the interfaces. We observe that the profiles are oscillatory in nature, representing waves of variable amplitude with distinct wave numbers propagating downstream and with no wave upstream. The observations are presented in tabular and graphical forms.

The tarani mutation alters surface curvature in Arabidopsis leaves by perturbing the patterns of surface expansion and cell division

The leaf surface usually stays flat, maintained by coordinated growth. Growth perturbation can introduce overall surface curvature, which can be negative, giving a saddle-shaped leaf, or positive, giving a cup-like leaf. Little is known about the molecular mechanisms that underlie leaf flatness, primarily because only a few mutants with altered surface curvature have been isolated and studied. Characterization of mutants of the CINCINNATA-like TCP genes in Antirrhinum and Arabidopsis have revealed that their products help maintain flatness by balancing the pattern of cell proliferation and surface expansion between the margin and the central zone during leaf morphogenesis. On the other hand, deletion of two homologous PEAPOD genes causes cup-shaped leaves in Arabidopsis due to excess division of dispersed meristemoid cells. Here, we report the isolation and characterization of an Arabidopsis mutant, tarani (tni), with enlarged, cup-shaped leaves. Morphometric analyses showed that the positive curvature of the tni leaf is linked to excess growth at the centre compared to the margin. By monitoring the dynamic pattern of CYCLIN D3;2 expression, we show that the shape of the primary arrest front is strongly convex in growing tni leaves, leading to excess mitotic expansion synchronized with excess cell proliferation at the centre. Reduction of cell proliferation and of endogenous gibberellic acid levels rescued the tni phenotype. Genetic interactions demonstrated that TNI maintains leaf flatness independent of TCPs and PEAPODs.

Maximum group velocity in a one-dimensional model with a sinusoidally varying staggered potential

We use Floquet theory to study the maximum value of the stroboscopic group velocity in a one-dimensional tight-binding model subjected to an on-site staggered potential varying sinusoidally in time. The results obtained by numerically diagonalizing the Floquet operator are analyzed using a variety of analytical schemes. In the low-frequency limit we use adiabatic theory, while in the high-frequency limit the Magnus expansion of the Floquet Hamiltonian turns out to be appropriate. When the magnitude of the staggered potential is much greater or much less than the hopping, we use degenerate Floquet perturbation theory; we find that dynamical localization occurs in the former case when the maximum group velocity vanishes. Finally, starting from an ``engineered'' initial state where the particles (taken to be hard-core bosons) are localized in one part of the chain, we demonstrate that the existence of a maximum stroboscopic group velocity manifests in a light-cone-like spreading of the particles in real space.

Phase diagram of the half-filled ionic Hubbard model

We study the phase diagram of the ionic Hubbard model (IHM) at half filling on a Bethe lattice of infinite connectivity using dynamical mean-field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered ionic potential Delta and the on-site Hubbard U. We find that for a finite Delta and at zero temperature, long-range antiferromagnetic (AFM) order sets in beyond a threshold U = U-AF via a first-order phase transition. For U smaller than U-AF the system is a correlated band insulator. Both methods show a clear evidence for a quantum transition to a half-metal (HM) phase just after the AFM order is turned on, followed by the formation of an AFM insulator on further increasing U. We show that the results obtained within both methods have good qualitative and quantitative consistency in the intermediate-to-strong-coupling regime at zero temperature as well as at finite temperature. On increasing the temperature, the AFM order is lost via a first-order phase transition at a transition temperature T-AF(U,Delta) or, equivalently, on decreasing U below U-AF(T,Delta)], within both methods, for weak to intermediate values of U/t. In the strongly correlated regime, where the effective low-energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) transition from the AFM phase to the paramagnetic phase, but the CTQMC does. At a finite temperature T, DMFT + CTQMC shows a second phase transition (not seen within DMFT + IPT) on increasing U beyond U-AF. At U-N > U-AF, when the Neel temperature T-N for the effective Heisenberg model becomes lower than T, the AFM order is lost via a second-order transition. For U >> Delta, T-N similar to t(2)/U(1 - x(2)), where x = 2 Delta/U and thus T-N increases with increase in Delta/U. In the three-dimensional parameter space of (U/t, T/t, and Delta/t), as T increases, the surface of first-order transition at U-AF(T,Delta) and that of the second-order transition at U-N(T,Delta) approach each other, shrinking the range over which the AFM order is stable. There is a line of tricritical points that separates the surfaces of first- and second-order phase transitions.

Mycobacterium tuberculosis has diminished capacity to counteract redox stress induced by elevated levels of endogenous superoxide

Mycobacterium tuberculosis (Mtb) has evolved protective and detoxification mechanisms to maintain cytoplasmic redox balance in response to exogenous oxidative stress encountered inside host phagocytes. In contrast, little is known about the dynamic response of this pathogen to endogenous oxidative stress generated within Mtb. Using a noninvasive and specific biosensor of cytoplasmic redox state of Mtb, we for first time discovered a surprisingly high sensitivity of this pathogen to perturbation in redox homeostasis induced by elevated endogenous reactive oxygen species (ROS). We synthesized a series of hydroquinone-based small molecule ROS generators and found that ATD-3169 permeated mycobacteria to reliably enhance endogenous ROS including superoxide radicals. When Mtb strains including multidrug-resistant (MDR) and extensively drug-resistant (XDR) patient isolates were exposed to this compound, a dose-dependent, long-lasting, and irreversible oxidative shift in intramycobacterial redox potential was detected. Dynamic redox potential measurements revealed that Mtb had diminished capacity to restore cytoplasmic redox balance in comparison with Mycobacterium smegmatis (Msm), a fast growing nonpathogenic mycobacterial species. Accordingly, Mtb strains were extremely susceptible to inhibition by ATD-3169 but not Msm, suggesting a functional linkage between dynamic redox changes and survival. Microarray analysis showed major realignment of pathways involved in redox homeostasis, central metabolism, DNA repair, and cell wall lipid biosynthesis in response to ATD-3169, all consistent with enhanced endogenous ROS contributing to lethality induced by this compound. This work provides empirical evidence that the cytoplasmic redox poise of Mtb is uniquely sensitive to manipulation in steady-state endogenous ROS levels, thus revealing the importance of targeting intramycobacterial redox metabolism for controlling TB infection. (C) 2015 The Authors. Published by Elsevier Inc.

Origin of nonlinearity and plausible turbulence by hydromagnetic transient growth in accretion disks: Faster growth rate than magnetorotational instability

We investigate the evolution of hydromagnetic perturbations in a small section of accretion disks. It is known that molecular viscosity is negligible in accretion disks. Hence, it has been argued that a mechanism, known as magnetorotational instability (MRI), is responsible for transporting matter in the presence of a weak magnetic field. However, there are some shortcomings, which question the effectiveness of MRI. Now the question arises, whether other hydromagnetic effects, e.g., transient growth (TG), can play an important role in bringing nonlinearity into the system, even at weak magnetic fields. In addition, it should be determined whether MRI or TG is primarily responsible for revealing nonlinearity in order to make the flow turbulent. Our results prove explicitly that the flows with a high Reynolds number (Re), which is the case for realistic astrophysical accretion disks, exhibit nonlinearity via TG of perturbation modes faster than that by modes producing MRI. For a fixed wave vector, MRI dominates over transient effects only at low Re, lower than the value expected to be in astrophysical accretion disks, and low magnetic fields. This calls into serious question the (overall) persuasiveness of MRI in astrophysical accretion disks.