Estimation of heavy rainfall during cyclonic storms from microwave observations using nonlinear approach over Indian Ocean

In this study, an effort has been made to study heavy rainfall events during cyclonic storms over Indian Ocean. This estimate is based on microwave observations from tropical rainfall measuring mission (TRMM) Microwave Imager (TMI). Regional scattering index (SI) developed for Indian region based on measurements at 19-, 21- and 85-GHz brightness temperature and polarization corrected temperature (PCT) at 85 GHz have been utilized in this study. These PCT and SI are collocated against Precipitation Radar (PR) onboard TRMM to establish a relationship between rainfall rate, PCT and SI. The retrieval technique using both linear and nonlinear regressions has been developed utilizing SI, PCT and the combination of SI and PCT. The results have been compared with the observations from PR. It was found that a nonlinear algorithm using combination of SI and PCT is more accurate than linear algorithm or nonlinear algorithm using either SI or PCT. Statistical comparison with PR exhibits the correlation coefficients (CC) of 0.68, 0.66 and 0.70, and root mean square error (RMSE) of 1.78, 1.96 and 1.68 mm/h from the observations of SI, PCT and combination of SI and PCT respectively using linear regressions. When nonlinear regression is used, the CC of 0.73, 0.71, 0.79 and RMSE of 1.64, 1.95, 1.54 mm/h are observed from the observations of SI, PCT and combination of SI and PCT, respectively. The error statistics for high rain events (above 10 mm/h) shows the CC of 0.58, 0.59, 0.60 and RMSE of 5.07, 5.47, 5.03 mm/h from the observations of SI, PCT and combination of SI and PCT, respectively, using linear regression, and on the other hand, use of nonlinear regression yields the CC of 0.66, 0.64, 0.71 and RMSE of 4.68, 5.78 and 4.02 mm/h from the observations of SI, PCT and combined SI and PCT, respectively.

Subdiffusion in hair bundle dynamics: The role of protein conformational fluctuations

The detection of sound signals in vertebrates involves a complex network of different mechano-sensory elements in the inner ear. An especially important element in this network is the hair bundle, an antenna-like array of stereocilia containing gated ion channels that operate under the control of one or more adaptation motors. Deflections of the hair bundle by sound vibrations or thermal fluctuations transiently open the ion channels, allowing the flow of ions through them, and producing an electrical signal in the process, eventually causing the sensation of hearing. Recent high frequency (0.1-10 kHz) measurements by Kozlov et al. Proc. Natl. Acad. Sci. U. S. A. 109, 2896 (2012)] of the power spectrum and the mean square displacement of the thermal fluctuations of the hair bundle suggest that in this regime the dynamics of the hair bundle are subdiffusive. This finding has been explained in terms of the simple Brownian motion of a filament connecting neighboring stereocilia (the tip link), which is modeled as a viscoelastic spring. In the present paper, the diffusive anomalies of the hair bundle are ascribed to tip link fluctuations that evolve by fractional Brownian motion, which originates in fractional Gaussian noise and is characterized by a power law memory. The predictions of this model for the power spectrum of the hair bundle and its mean square displacement are consistent with the experimental data and the known properties of the tip link. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4768902]

The Coil-to-Helix Transition in IlvN Regulates the Allosteric Control of Escherichia coli Acetohydroxyacid Synthase I

The solution structure of IlvN, the regulatory subunit of Escherichia coil acetohydroxyacid synthase I, in the valine-bound form has been determined using high-resolution multidimensional, multinuclear nuclear magnetic resonance (NMR) methods. IlvN in the presence or absence of the effector molecule is present as a 22.5 kDa dimeric molecule. The ensemble of 20 low-energy structures shows a backbone root-mean-square deviation of 0.73 +/- 0.13 angstrom and a root-mean-square deviation of 1.16 +/- 0.13 angstrom for all heavy atoms. Furthermore, more than 98% of the backbone phi and psi dihedral angles occupy the allowed and additionally allowed regions of the Ramachandran map, which is indicative of the fact that the structures are of high stereochemical quality. Each protomer exhibits a beta alpha beta beta alpha beta alpha topology that is a characteristic feature of the ACT domain seen in metabolic enzymes. In the valine-bound form, IlvN exists apparently as a single conformer. In the free form, IlvN exists as a mixture of conformational states that are in intermediate exchange on the NMR time scale. Thus, a large shift in the conformational equilibrium is observed upon going from the free form to the bound form. The structure of the valine-bound form of IlvN was found to be similar to that of the ACT domain of the unliganded form of IlvH. Comparisons of the structures of the unliganded forms of these proteins suggest significant differences. The structural and conformational properties of IlvN determined here have allowed a better understanding of the mechanism of regulation of branched chain amino acid biosynthesis.

Cramer-Rao-Type Bounds for Sparse Bayesian Learning

In this paper, we derive Hybrid, Bayesian and Marginalized Cramer-Rao lower bounds (HCRB, BCRB and MCRB) for the single and multiple measurement vector Sparse Bayesian Learning (SBL) problem of estimating compressible vectors and their prior distribution parameters. We assume the unknown vector to be drawn from a compressible Student-prior distribution. We derive CRBs that encompass the deterministic or random nature of the unknown parameters of the prior distribution and the regression noise variance. We extend the MCRB to the case where the compressible vector is distributed according to a general compressible prior distribution, of which the generalized Pareto distribution is a special case. We use the derived bounds to uncover the relationship between the compressibility and Mean Square Error (MSE) in the estimates. Further, we illustrate the tightness and utility of the bounds through simulations, by comparing them with the MSE performance of two popular SBL-based estimators. We find that the MCRB is generally the tightest among the bounds derived and that the MSE performance of the Expectation-Maximization (EM) algorithm coincides with the MCRB for the compressible vector. We also illustrate the dependence of the MSE performance of SBL based estimators on the compressibility of the vector for several values of the number of observations and at different signal powers.

Joint data detection and dominant singular mode estimation in time varying reciprocal MIMO systems

This paper proposes an algorithm for joint data detection and tracking of the dominant singular mode of a time varying channel at the transmitter and receiver of a time division duplex multiple input multiple output beamforming system. The method proposed is a modified expectation maximization algorithm which utilizes an initial estimate to track the dominant modes of the channel at the transmitter and the receiver blindly; and simultaneously detects the un known data. Furthermore, the estimates are constrained to be within a confidence interval of the previous estimate in order to improve the tracking performance and mitigate the effect of error propagation. Monte-Carlo simulation results of the symbol error rate and the mean square inner product between the estimated and the true singular vector are plotted to show the performance benefits offered by the proposed method compared to existing techniques.

A risk-estimation-based formulation for speech enhancement and its relation to Wiener filtering

The goal of speech enhancement algorithms is to provide an estimate of clean speech starting from noisy observations. The often-employed cost function is the mean square error (MSE). However, the MSE can never be computed in practice. Therefore, it becomes necessary to find practical alternatives to the MSE. In image denoising problems, the cost function (also referred to as risk) is often replaced by an unbiased estimator. Motivated by this approach, we reformulate the problem of speech enhancement from the perspective of risk minimization. Some recent contributions in risk estimation have employed Stein's unbiased risk estimator (SURE) together with a parametric denoising function, which is a linear expansion of threshold/bases (LET). We show that the first-order case of SURE-LET results in a Wiener-filter type solution if the denoising function is made frequency-dependent. We also provide enhancement results obtained with both techniques and characterize the improvement by means of local as well as global SNR calculations.

Structure, molecular dynamics, and stress in a linear polymer

We present a study correlating uniaxial stress in a polymer with its underlying structure when it is strained. The uniaxial stress is significantly influenced by the mean-square bond length and mean bond angle. In contrast, the size and shape of the polymer, typically represented by the end-to-end length, mass ratio, and radius of gyration, contribute negligibly. Among externally set control variables, density and polymer chain length play a critical role in influencing the anisotropic uniaxial stress. Short chain polymers more or less behave like rigid molecules. Temperature and rate of loading, in the range considered, have a very mild effect on the uniaxial stress.

Composition, structure and electrical properties of DC reactive magnetron sputtered Al2O3 thin films

Thin films of alumina (Al2O3) were deposited over Si < 1 0 0 > substrates at room temperature at an oxygen gas pressure of 0.03 Pa and sputtering power of 60 W using DC reactive magnetron sputtering. The composition of the as-deposited film was analyzed by X-ray photoelectron spectroscopy and the O/Al atomic ratio was found to be 1.72. The films were then annealed in vacuum to 350, 550 and 750 degrees C and X-ray diffraction results revealed that both as-deposited and post deposition annealed films were amorphous. The surface morphology and topography of the films was studied using scanning electron microscopy and atomic force microscopy, respectively. A progressive decrease in the root mean square (RMS) roughness of the films from 1.53 nm to 0.7 nm was observed with increase in the annealing temperature. Al-Al2O3-Al thin film capacitors were then fabricated on p-type Si < 1 0 0 > substrate to study the effect of temperature and frequency on the dielectric property of the films and the results are discussed.

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.

Solution nuclear magnetic resonance structure of the GATase subunit and structural Basis of the interaction between GATase and ATPPase subunits in a two-subunit-type GMPS from methanocaldococcus jannaschii

The solution structure of the monomeric glutamine amidotransferase (GATase) subunit of the Methanocaldococcus janaschii (Mj) guanosine monophosphate synthetase (GMPS) has been determined using high-resolution nuclear magnetic resonance methods. Gel filtration chromatography and N-15 backbone relaxation studies have shown that the Mj GATase subunit is present in solution as a 21 kDa (188-residue) monomer. The ensemble of 20 lowest-energy structures showed root-mean-square deviations of 0.35 +/- 0.06 angstrom for backbone atoms and 0.8 +/- 0.06 angstrom for all heavy atoms. Furthermore, 99.4% of the backbone dihedral angles are present in the allowed region of the Ramachandran map, indicating the stereochemical quality of the structure. The core of the tertiary structure of the GATase is composed of a seven-stranded mixed beta-sheet that is fenced by five alpha-helices. The Mj GATase is similar in structure to the Pyrococcus horikoshi (Ph) GATase subunit. Nuclear magnetic resonance (NMR) chemical shift perturbations and changes in line width were monitored to identify residues on GATase that were responsible for interaction with magnesium and the ATPPase subunit, respectively. These interaction studies showed that a common surface exists for the metal ion binding as well as for the protein-protein interaction. The dissociation constant for the GATase-Mg2+ interaction has been found to be similar to 1 mM, which implies that interaction is very weak and falls in the fast chemical exchange regime. The GATase-ATPPase interaction, on the other hand, falls in the intermediate chemical exchange regime on the NMR time scale. The implication of this interaction in terms of the regulation of the GATase activity of holo GMPS is discussed.

Feature Search in the Grassmanian in Online Reinforcement Learning

We consider the problem of finding the best features for value function approximation in reinforcement learning and develop an online algorithm to optimize the mean square Bellman error objective. For any given feature value, our algorithm performs gradient search in the parameter space via a residual gradient scheme and, on a slower timescale, also performs gradient search in the Grassman manifold of features. We present a proof of convergence of our algorithm. We show empirical results using our algorithm as well as a similar algorithm that uses temporal difference learning in place of the residual gradient scheme for the faster timescale updates.

Cooperative Particle Swarm Optimization Based Receiver for Large-Dimension MIMO-ZPSC Systems

In this paper, we propose a cooperative particle swarm optimization (CPSO) based channel estimation/equalization scheme for multiple-input multiple-output zero-padded single-carrier (MIMO-ZPSC) systems with large dimensions in frequency selective channels. We estimate the channel state information at the receiver in time domain using a PSO based algorithm during training phase. Using the estimated channel, we perform information symbol detection in the frequency domain using FFT based processing. For this detection, we use a low complexity OLA (OverLap Add) likelihood ascent search equalizer which uses minimum mean square (MMSE) equalizer solution as the initial solution. Multiple iterations between channel estimation and data detection are carried out which significantly improves the mean square error and bit error rate performance of the receiver.

Channel estimation at the transmitter in a reciprocal MIMO spatial multiplexing system

This paper considers the problem of channel estimation at the transmitter in a spatial multiplexing-based Time Division Duplex (TDD) Multiple Input Multiple Output (MIMO) system with perfect CSIR. A novel channel-dependent Reverse Channel Training (RCT) sequence is proposed, using which the transmitter estimates the beamforming vectors for forward link data transmission. This training sequence is designed based on the following two metrics: (i) a capacity lower bound, and (ii) the mean square error in the estimate. The performance of the proposed training scheme is analyzed and is shown to significantly outperform the conventional orthogonal RCT sequence. Also, in the case where the transmitter uses water-filling power allocation for data transmission, a novel RCT sequence is proposed and optimized with respect to the MSE in estimating the transmit covariance matrix.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Channel estimation relying on the minimum bit-error-ratio criterion for BPSK and QPSK signals

The authors consider the channel estimation problem in the context of a linear equaliser designed for a frequency selective channel, which relies on the minimum bit-error-ratio (MBER) optimisation framework. Previous literature has shown that the MBER-based signal detection may outperform its minimum-mean-square-error (MMSE) counterpart in the bit-error-ratio performance sense. In this study, they develop a framework for channel estimation by first discretising the parameter space and then posing it as a detection problem. Explicitly, the MBER cost function (CF) is derived and its performance studied, when transmitting binary phase shift keying (BPSK) and quadrature phase shift keying (QPSK) signals. It is demonstrated that the MBER based CF aided scheme is capable of outperforming existing MMSE, least square-based solutions.