Diffraction tomography from intensity measurements: an evolutionary stochastic search to invert experimental data

We develop iterative diffraction tomography algorithms, which are similar to the distorted Born algorithms, for inverting scattered intensity data. Within the Born approximation, the unknown scattered field is expressed as a multiplicative perturbation to the incident field. With this, the forward equation becomes stable, which helps us compute nearly oscillation-free solutions that have immediate bearing on the accuracy of the Jacobian computed for use in a deterministic Gauss-Newton (GN) reconstruction. However, since the data are inherently noisy and the sensitivity of measurement to refractive index away from the detectors is poor, we report a derivative-free evolutionary stochastic scheme, providing strictly additive updates in order to bridge the measurement-prediction misfit, to arrive at the refractive index distribution from intensity transport data. The superiority of the stochastic algorithm over the GN scheme for similar settings is demonstrated by the reconstruction of the refractive index profile from simulated and experimentally acquired intensity data. (C) 2014 Optical Society of America

EFFECTIVE Q CRITERION FOR DISK STABILITY IN AN EXTERNAL POTENTIAL

standard Q criterion (with Q > 1) describes the stability against local, axisymmetric perturbations in a disk supported by rotation and random motion. Most astrophysical disks, however, are under the influence of an external gravitational potential, which can significantly affect their stability. A typical example is a galactic disk embedded in a dark matter halo. Here, we do a linear perturbation analysis for a disk in an external potential and obtain a generalized dispersion relation and the effective stability criterion. An external potential, such as that due to the dark matter halo concentric with the disk, contributes to the unperturbed rotational field and significantly increases its stability. We obtain the values for the effective Q parameter for the Milky Way and for a low surface brightness galaxy, UGC 7321. We find that in each case the stellar disk by itself is barely stable and it is the dark matter halo that stabilizes the disk against local, axisymmetric gravitational instabilities. Thus, the dark matter halo is necessary to ensure local disk stability. This result has been largely missed so far because in practice the Q parameter for a galactic disk is obtained using the observed rotational field that already includes the effect of the halo.

Insights into the AIEE of 1,8-Naphthalimides (NPIs): Inverse Effects of Intermolecular Interactions in Solution and Aggregates

Systematic structural perturbation has been used to fine-tune and understand the luminescence properties of three new 1,8-naphthalimides (NPIs) in solution and aggregates. The NPIs show blue emission in the solution state and their fluorescence quantum yields are dependent upon their molecular rigidity. In concentrated solutions of the NPIs, intermolecular interactions were found to quench the fluorescence due to the formation of excimers. In contrast, upon aggregation (in THF/H2O mixtures), the NPIs show aggregation-induced emission enhancement (AIEE). The NPIs also show moderately high solid-state emission quantum yields (ca. 10-12.7 %). The AIEE behaviour of the NPIs depends on their molecular rigidity and the nature of their intermolecular interactions. The NPIs 1-3 show different extents of intermolecular (pi-pi and C-H center dot center dot center dot O) interactions in their solid-state crystal structures depending on their substituents. Detailed photophysical, computational and structural investigations suggest that an optimal balance of structural flexibility and intermolecular communication is necessary for achieving AIEE characteristics in these NPIs.

Dynamics of the reaction between the free end of a tethered self-avoiding polymer and a flat penetrable surface: A renormalization group study

The average time tau(r) for one end of a long, self-avoiding polymer to interact for the first time with a flat penetrable surface to which it is attached at the other end is shown here to scale essentially as the square of the chain's contour length N. This result is obtained within the framework of the Wilemski-Fixman approximation to diffusion-limited reactions, in which the reaction time is expressed as a time correlation function of a ``sink'' term. In the present work, this sink-sink correlation function is calculated using perturbation expansions in the excluded volume and the polymer-surface interactions, with renormalization group methods being used to resum the expansion into a power law form. The quadratic dependence of tau(r) on N mirrors the behavior of the average time tau(c) of a free random walk to cyclize, but contrasts with the cyclization time of a free self-avoiding walk (SAW), for which tau(r) similar to N-2.2. A simulation study by Cheng and Makarov J. Phys. Chem. B 114, 3321 (2010)] of the chain-end reaction time of an SAW on a flat impenetrable surface leads to the same N-2.2 behavior, which is surprising given the reduced conformational space a tethered polymer has to explore in order to react. (C) 2014 AIP Publishing LLC.

Boundary perturbations and the Helmholtz equation in three dimensions

We propose an analytic perturbative scheme in the spirit of Lord Rayleigh's work for determining the eigenvalues of the Helmholtz equation in three dimensions inside an arbitrary boundary where the eigenfunction satisfies either the Dirichlet boundary condition or the Neumann boundary condition. Although numerous works are available in the literature for arbitrary boundaries in two dimensions, to the best of our knowledge the formulation in three dimensions is proposed for the first time. In this novel prescription, we have expanded the arbitrary boundary in terms of spherical harmonics about an equivalent sphere and obtained perturbative closed-form solutions at each order for the problem in terms of corrections to the equivalent spherical boundary for both the boundary conditions. This formulation is in parallel with the standard time-independent Rayleigh-Schrodinger perturbation theory. The efficacy of the method is tested by comparing the perturbative values against the numerically calculated eigenvalues for spheroidal, superegg and superquadric shaped boundaries. It is shown that this perturbation works quite well even for wide departure from spherical shape and for higher excited states too. We believe this formulation would find applications in the field of quantum dots and acoustical cavities.

Universal power law in crossover from integrability to quantum chaos

We study models of interacting fermions in one dimension to investigate the crossover from integrability to nonintegrability, i.e., quantum chaos, as a function of system size. Using exact diagonalization of finite-sized systems, we study this crossover by obtaining the energy level statistics and Drude weight associated with transport. Our results reinforce the idea that for system size L -> infinity nonintegrability sets in for an arbitrarily small integrability-breaking perturbation. The crossover value of the perturbation scales as a power law similar to L-eta when the integrable system is gapless. The exponent eta approximate to 3 appears to be robust to microscopic details and the precise form of the perturbation. We conjecture that the exponent in the power law is characteristic of the random matrix ensemble describing the nonintegrable system. For systems with a gap, the crossover scaling appears to be faster than a power law.

Finite size scaling in crossover among different random matrix ensembles in microscopic lattice models

Using numerical diagonalization we study the crossover among different random matrix ensembles (Poissonian, Gaussian orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE) and Gaussian symplectic ensemble (GSE)) realized in two different microscopic models. The specific diagnostic tool used to study the crossovers is the level spacing distribution. The first model is a one-dimensional lattice model of interacting hard-core bosons (or equivalently spin 1/2 objects) and the other a higher dimensional model of non-interacting particles with disorder and spin-orbit coupling. We find that the perturbation causing the crossover among the different ensembles scales to zero with system size as a power law with an exponent that depends on the ensembles between which the crossover takes place. This exponent is independent of microscopic details of the perturbation. We also find that the crossover from the Poissonian ensemble to the other three is dominated by the Poissonian to GOE crossover which introduces level repulsion while the crossover from GOE to GUE or GOE to GSE associated with symmetry breaking introduces a subdominant contribution. We also conjecture that the exponent is dependent on whether the system contains interactions among the elementary degrees of freedom or not and is independent of the dimensionality of the system.

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.