A novel filtering framework through Girsanov correction for the identification of nonlinear dynamical systems

Using a Girsanov change of measures, we propose novel variations within a particle-filtering algorithm, as applied to the inverse problem of state and parameter estimations of nonlinear dynamical systems of engineering interest, toward weakly correcting for the linearization or integration errors that almost invariably occur whilst numerically propagating the process dynamics, typically governed by nonlinear stochastic differential equations (SDEs). Specifically, the correction for linearization, provided by the likelihood or the Radon-Nikodym derivative, is incorporated within the evolving flow in two steps. Once the likelihood, an exponential martingale, is split into a product of two factors, correction owing to the first factor is implemented via rejection sampling in the first step. The second factor, which is directly computable, is accounted for via two different schemes, one employing resampling and the other using a gain-weighted innovation term added to the drift field of the process dynamics thereby overcoming the problem of sample dispersion posed by resampling. The proposed strategies, employed as add-ons to existing particle filters, the bootstrap and auxiliary SIR filters in this work, are found to non-trivially improve the convergence and accuracy of the estimates and also yield reduced mean square errors of such estimates vis-a-vis those obtained through the parent-filtering schemes.

Effect of temperature and strain rate on the deformation behavior and microstructure of a homogenized AZ31 magnesium alloy

The plastic deformation behavior and dynamic recrystallization (DRX) in homogenized AZ31 Mg alloy was investigated in uniaxial compression in the temperature range between 150 and 400 degrees C with strain rates ranging from 10(-3) to 10(2) s(-1). Twinning was found to contribute significantly during the early stages of deformation. The onset of twinning was examined in detail by recourse to the examination of the appearance of first local maxima before peak strain in the stress-strain responses and the second derivative of stress with strain. High strain hardening rate was observed immediately after the onset of twinning and was found to increase with the Zener-Hollomon parameter. DRX was observed at temperatures above 250 degrees C whereas deformation at lower temperatures (< 250 degrees C) leads to extensive twinning at all the strain rates. At intermediate temperatures of 250-300 degrees C, plastic strains tend to localize near grain/twin boundaries, confining DRX only to these regions. Increase in the temperature promotes non-basal slip, which, in turn, leads to uniform deformation; DRX too becomes uniform. Deformation behavior in three different regimes of temperature is discussed. The dependence of critical stress for the onset of DRX and peak flow stress on temperature and strain rate are also described. (C) 2013 Elsevier B.V. All rights reserved.

The first teresomatan caecilian (Amphibia: Gymnophiona) from the Eastern Ghats of India-a new species of Gegeneophis Peters, 1880

A new species of caecilian amphibian, Gegeneophis orientalis sp. nov., is described based on a series of nine specimens from high elevation (ca. 1,200 m) habitats in the Eastern Ghats in the states of Andhra Pradesh and Odisha, India. This species differs from all other congeners in having only bicuspid teeth in the outer as well as inner rows. The new species is the first caecilian reported from the state of Odisha, the first teresomatan caecilian from the Eastern Ghats, and is the only Indian indotyphlid known from outside the Western Ghats region.

Cyclic AMP in Mycobacteria: the second messenger comes first

Cyclic AMP (cAMP) has emerged as a pivotal molecule for signalling in all life forms. Mycobacterial genomes have been found to encode for numerous proteins that are involved in cAMP generation, degradation and utilization. Many of these proteins have domain organizations unique to mycobacteria. This review summarizes recent advances in mechanisms of cAMP synthesis and degradation, focusing on the processes by which cAMP modulates mycobacterial signalling. We explore its impact on the physiology of the organism and on the discourse between M. tuberculosis and its host.

Tb3+ sensitization in a deoxycholate organogel matrix, and selective quenching of luminescence by an aromatic nitro derivative

In this article, we present the discovery of a metallo-organogel derived from a Tb3+ salt and sodium deoxycholate (NaDCh) in methanol. The gel was made luminescent through sensitization of Tb3+ by doping with 2,3-dihydroxynaphthalene (DHN) in micromolar concentrations. Rheological measurements of the mechanical properties of the organogel confirmed the characteristics of a true gel. Significant quenching of Tb3+ luminescence was observed in the deoxycholate gel matrix by 2,4,7-trinitrofluorenone (TNF), but not by several other polynitro aromatics. Microscopic studies (AFM, TEM and SEM) revealed a highly entangled fibrous network. The xerogels retained luminescent properties suggesting the possibility for application in coatings, etc.

A Nearly Exact Reformulation of the Girsanov Linearization for Stochastically Driven Nonlinear Oscillators

The Girsanov linearization method (GLM), proposed earlier in Saha, N., and Roy, D., 2007, ``The Girsanov Linearisation Method for Stochastically Driven Nonlinear Oscillators,'' J. Appl. Mech., 74, pp. 885-897, is reformulated to arrive at a nearly exact, semianalytical, weak and explicit scheme for nonlinear mechanical oscillators under additive stochastic excitations. At the heart of the reformulated linearization is a temporally localized rejection sampling strategy that, combined with a resampling scheme, enables selecting from and appropriately modifying an ensemble of locally linearized trajectories while weakly applying the Girsanov correction (the Radon-Nikodym derivative) for the linearization errors. The semianalyticity is due to an explicit linearization of the nonlinear drift terms and it plays a crucial role in keeping the Radon-Nikodym derivative ``nearly bounded'' above by the inverse of the linearization time step (which means that only a subset of linearized trajectories with low, yet finite, probability exceeds this bound). Drift linearization is conveniently accomplished via the first few (lower order) terms in the associated stochastic (Ito) Taylor expansion to exclude (multiple) stochastic integrals from the numerical treatment. Similarly, the Radon-Nikodym derivative, which is a strictly positive, exponential (super-) martingale, is converted to a canonical form and evaluated over each time step without directly computing the stochastic integrals appearing in its argument. Through their numeric implementations for a few low-dimensional nonlinear oscillators, the proposed variants of the scheme, presently referred to as the Girsanov corrected linearization method (GCLM), are shown to exhibit remarkably higher numerical accuracy over a much larger range of the time step size than is possible with the local drift-linearization schemes on their own.

Investigation on non-exchange spring behaviour and exchange spring behaviour: A first order reversal curve analysis

The First Order Reversal Curve (FORC) method has been utilised to understand the magnetization reversal and the extent of the irreversible magnetization of the soft CoFe2O4-hard SrFe12O19 nanocomposite in the nonexchange spring and the exchange spring regime. The single peak switching behaviour in the FORC distribution of the exchange spring composite confirms the coherent reversal of the soft and hard phases. The onset of the nucleation field and the magnetization reversal by domain wall movement are also evident from the FORC measurements. (C) 2013 AIP Publishing LLC.

A prototype of a first-responder indoor localization system

In this paper we present an approach to build a prototype. model of a first-responder localization system intended for disaster relief operations. This system is useful to monitor and track the positions of the first-responders in an indoor environment, where GPS is not available. Each member of the first responder team is equipped with two zero-velocity-update-aided inertial navigation systems, one on each foot, a camera mounted on a helmet, and a processing platform strapped around the waist of the first responder, which fuses the data from the different sensors. The fusion algorithm runs real-time on the processing platform. The video is also processed using the DSP core of the computing machine. The processed data consisting of position, velocity, heading information along with video streams is transmitted to the command and control system via a local infrastructure WiFi network. A centralized cooperative localization algorithm, utilizing the information from Ultra Wideband based inter-agent ranging devices combined with the position estimates and uncertainties of each first responder, has also been implemented.

Base Triggered Enhancement of First Hyperpolarizability of a Keto-Enol Tautomer

This work describes the base triggered enhancement of first hyperpolarizability of a tautomeric organic molecule, namely, benzoylacetanilide (BA). We have used the hyper-Rayleigh scattering technique to measure the first hyperpolarizability (beta) of BA which exists in the pure keto form in water and as a keto-enol tautomer in ethanol. Its anion exists in equilibrium with the keto and enol forms at pH 11 in aqueous solution. The beta value of the anion form is 709 X 10(-30) esu, whereas that of the enol is 232 x 10(-3) esu and of the keto is 88 X 10(-30) esu. There is an enhancement of beta by similar to 8 times for the anion and similar to 3 times for the enol compared to the keto form. All these are achieved by altering the equilibrium between the three forms of BA by simple means. MP2 calculations reproduce the experimental trend, but the computed beta values are much lower than the measured values. DFT calculations with the standard B3LYP functional could not predict the right order in the beta values. The difference between experimental and calculated values is, perhaps, due to the fact that electron correlation effects are important in computing optical nonlinearities of large organic molecules and MP2 and B3LYP calculations done here for different forms of BA could not account for such effects adequately.

Synthesis, structural studies and catalytic activity of dioxidomolybdenum(VI) complexes with aroylhydrazones of naphthol-derivative

Reaction of the salicylhydrazone of 2-hydroxy-1-naphthaldehyde (H2L1), anthranylhydrazone of 2hydroxy-l-naphthaldehyde (H2L2), benzoylhydrazone of 2-hydroxy-1-acetonaphthone (H2L3) and anthranylhydrazone of 2-hydroxy-1-acetonaphthone (H2L4; general abbreviation H2L) with MoO2(acac)21 afforded a series of 5- and 6- coordinate Mo(VI) complexes of the type MoO2L1-2(ROH)] where R = C2H5 (1) and CH3 (2)], and MoO2L3-4] (3 and 4). The substrate binding capacity of 1 has been demonstrated by the formation of one mononuclear mixed-ligand dioxidomolybdenum complex MoO2L1(Q)] (where Q= gamma-picoline (la)). Molecular structure of all the complexes (I, la, 2,3 and 4) is determined by X-ray crystallography, demonstrating the dibasic tridentate behavior of ligands. All the complexes show two irreversible reductive responses within the potential window -0.73 to -1.08 V, due to Movl/Mov and Mov/Mow processes. Catalytic potential of these complexes was tested for the oxidation of benzoin using 30% aqueous H2O2 as an oxidant in methanol. At least four reaction products, benzoic acid, benzaldehydedimethylacetal, methyl benzoate and benzil were obtained with the 95-99% conversion under optimized reaction conditions. Oxidative bromination of salicylaldehyde, a functional mimic of haloperoxidases, in aqueous 1-1202/KEr in the presence of HC1O4 at room temperature has also been carried out successfully. (C) 2013 Elsevier Ltd. All rights reserved.

First-principles DFT plus GW study of oxygen vacancies in rutile TiO2

We perform first-principles calculations of the quasiparticle defect states, charge transition levels, and formation energies of oxygen vacancies in rutile titanium dioxide. The calculations are done within the recently developed combined DFT + GW formalism, including the necessary electrostatic corrections for the supercells with charged defects. We find the oxygen vacancy to be a negative U defect, where U is the defect electron addition energy. For Fermi level values below similar to 2.8 eV (relative to the valence-band maximum), we find the +2 charge state of the vacancy to be the most stable, while above 2.8 eV we find that the neutral charge state is the most stable.

First-Principles Study of Structure, Vibrational, and Elastic Properties of Stoichiometric and Calcium-Deficient Hydroxyapatite

Hydroxyapatite (HAp), a primary constituent of human bone, is usually nonstoichiometric with varying Ca/P molar ratios, with the well-known fact that Ca deficiency can cause marked reductions in its mechanical properties. To gain insights into the mechanism of this degradation, we employ first-principles calculations based on density functional theory and determine the effects of Ca deficiency on structure, vibrational, and elastic properties of HAp. Our simulation results confirm a considerable reduction in the elastic constants of HAp due to Ca deficiency, which was experimentally reported earlier. Stress-induced transformation of the Ca-deficient defected structure into a metastable state upon the application of stress could be a reason for this. Local structural stability of HAp and Ca-deficient HAp structures is assessed with full phonon dispersion studies. Further, specific signatures in the computed vibrational spectra for Ca deficiency in HAp can be utilized in experimental characterization of different types of defected HAp.

Holographic stress tensor at finite coupling

We calculate one, two and three point functions of the holographic stress tensor for any bulk Lagrangian of the form L (g(ab), R-abcd, del(e) R-abcd). Using the first law of entanglement, a simple method has recently been proposed to compute the holographic stress tensor arising from a higher derivative gravity dual. The stress tensor is proportional to a dimension dependent factor which depends on the higher derivative couplings. In this paper, we identify this proportionality constant with a B-type trace anomaly in even dimensions for any bulk Lagrangian of the above form. This in turn relates to C-T, the coefficient appearing in the two point function of stress tensors. We use a background field method to compute the two and three point function of stress tensors for any bulk Lagrangian of the above form in arbitrary dimensions. As an application we consider general situations where eta/s for holographic plasmas is less than the KSS bound.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.