Observation of a power-law memory kernel for fluctuations within a single protein molecule

The fluctuation of the distance between a fluorescein-tyrosine pair within a single protein complex was directly monitored in real time by photoinduced electron transfer and found to be a stationary, time-reversible, and non-Markovian Gaussian process. Within the generalized Langevin equation formalism, we experimentally determine the memory kernel K(t), which is proportional to the autocorrelation function of the random fluctuating force. K(t) is a power-law decay, t(-0.51 +/- 0.07) in a broad range of time scales (10(-3)-10 s). Such a long-time memory effect could have implications for protein functions.

Primordial Perturbations from a Self-interacting Curvaton

Inflation is a period of accelerated expansion in the very early universe, which has the appealing aspect that it can create primordial perturbations via quantum fluctuations. These primordial perturbations have been observed in the cosmic microwave background, and these perturbations also function as the seeds of all large-scale structure in the universe. Curvaton models are simple modifications of the standard inflationary paradigm, where inflation is driven by the energy density of the inflaton, but another field, the curvaton, is responsible for producing the primordial perturbations. The curvaton decays after inflation as ended, where the isocurvature perturbations of the curvaton are converted into adiabatic perturbations. Since the curvaton must decay, it must have some interactions. Additionally realistic curvaton models typically have some self-interactions. In this work we consider self-interacting curvaton models, where the self-interaction is a monomial in the potential, suppressed by the Planck scale, and thus the self-interaction is very weak. Nevertheless, since the self-interaction makes the equations of motion non-linear, it can modify the behaviour of the model very drastically. The most intriguing aspect of this behaviour is that the final properties of the perturbations become highly dependent on the initial values. Departures of Gaussian distribution are important observables of the primordial perturbations. Due to the non-linearity of the self-interacting curvaton model and its sensitivity to initial conditions, it can produce significant non-Gaussianity of the primordial perturbations. In this work we investigate the non-Gaussianity produced by the self-interacting curvaton, and demonstrate that the non-Gaussianity parameters do not obey the analytically derived approximate relations often cited in the literature. Furthermore we also consider a self-interacting curvaton with a mass in the TeV-scale. Motivated by realistic particle physics models such as the Minimally Supersymmetric Standard Model, we demonstrate that a curvaton model within the mass range can be responsible for the observed perturbations if it can decay late enough.

Observation of the Ωb- baryon and measurement of the properties of the Ξb- and Ωb- baryons

We report the observation of the bottom, doubly-strange baryon Omega^-_b through the decay chain Omega^-_b -> J/psi Omega^-, where J/psi -> mu^+ mu^-, Omega^- -> Lambda K^-, and Lambda -> p pi^-, using 4.2 fb^{-1} of data from p\bar p collisions at sqrt{s}=1.96 TeV, and recorded with the Collider Detector at Fermilab. A signal is observed whose probability of arising from a background fluctuation is 4.0 * 10^{-8}, or 5.5 Gaussian standard deviations. The Omega^-_b mass is measured to be 6054.4 +/- 6.8 (stat.) +/- 0.9 (syst.) MeV/c^2. The lifetime of the Omega^-_b baryon is measured to be 1.13^{+0.53}_{-0.40}(stat.) +/- 0.02(syst.)$ ps. In addition, for the \Xi^-_b baryon we measure a mass of 5790.9 +/- 2.6(stat.) +/- 0.8(syst.) MeV/c^2 and a lifetime of 1.56^{+0.27}_{-0.25}(stat.) +/-0.02(syst.) ps. Under the assumption that the \Xi_b^- and \Omega_b^- are produced with similar kinematic distributions to the \Lambda^0_b baryon, we find sigma(Xi_b^-) B(Xi_b^- -> J/psi Xi^-)}/ sigma(Lambda^0_b) B(Lambda^0_b -> J/psi Lambda)} = 0.167^{+0.037}_{-0.025}(stat.) +/-0.012(syst.) and sigma(Omega_b^-) B(Omega_b^- -> J/psi Omega^-)/ sigma(Lambda^0_b) B(Lambda^0_b -> J/psi Lambda)} = 0.045^{+0.017}_{-0.012}(stat.) +/- 0.004(syst.) for baryons produced with transverse momentum in the range of 6-20 GeV/c.

Compliance of rough cylinders in compression

Compression of a rough turned cylinder between two hard, smooth, flat plates has been analysed with the aid of a mathematical model based on statistical analysis. It is assumed that the asperity peak heights follow Gaussian or normal and beta distribution functions and that the loaded asperities comply as though they are completely isolated from the neighbouring ones. Equations have been developed for the loadcompliance relation of the real surface using a simplified relation of the form W0 = K1δn for the load-compliance of a single asperity. Parameters K1 and n have considerable influence on the load-compliance curve and they depend on the material, tip angle of the asperity, standard deviation of the asperity peak height distribution and the density of the asperities.

First passage probability during random vibration

The probability that a random process crosses an arbitrary level for the first time is expressed as a Gram—Charlier series, the leading term of which is the Poisson approximation. The coefficients of this series are related to the moments of the number of level crossings. The results are applicable to both stationary and non-stationary processes. Some numerical results are presented for the response process of a linear single-degree-of-freedom oscillator under Gaussian white noise excitation.

Moment-Matched Lognormal Modeling of Uplink Interference with Power Control and Cell Selection

We develop an alternate characterization of the statistical distribution of the inter-cell interference power observed in the uplink of CDMA systems. We show that the lognormal distribution better matches the cumulative distribution and complementary cumulative distribution functions of the uplink interference than the conventionally assumed Gaussian distribution and variants based on it. This is in spite of the fact that many users together contribute to uplink interference, with the number of users and their locations both being random. Our observations hold even in the presence of power control and cell selection, which have hitherto been used to justify the Gaussian distribution approximation. The parameters of the lognormal are obtained by matching moments, for which detailed analytical expressions that incorporate wireless propagation, cellular layout, power control, and cell selection parameters are developed. The moment-matched lognormal model, while not perfect, is an order of magnitude better in modeling the interference power distribution.

Nonlinear vibration analysis of composite laminated and sandwich plates with random material properties

Nonlinear vibration analysis is performed using a C-0 assumed strain interpolated finite element plate model based on Reddy's third order theory. An earlier model is modified to include the effect of transverse shear variation along the plate thickness and Von-Karman nonlinear strain terms. Monte Carlo Simulation with Latin Hypercube Sampling technique is used to obtain the variance of linear and nonlinear natural frequencies of the plate due to randomness in its material properties. Numerical results are obtained for composite plates with different aspect ratio, stacking sequence and oscillation amplitude ratio. The numerical results are validated with the available literature. It is found that the nonlinear frequencies show increasing non-Gaussian probability density function with increasing amplitude of vibration and show dual peaks at high amplitude ratios. This chaotic nature of the dispersion of nonlinear eigenvalues is also r

Simple kinetic sensor to structural transitions

Driven nonequilibrium structural phase transformation has been probed using time-varying resistance fluctuations or noise. We demonstrate that the non-Gaussian component (NGC) of noise obtained by evaluating the higher-order statistics of fluctuations, serves as a simple kinetic detector of these phase transitions. Using the Martensite transformation in free-standing wires of nickel-titanium binary alloys as a prototype, we observe clear deviations from the Gaussian background in the transformation zone, indicative of the long-range correlations in the system as the phase transforms. The viability of non-Gaussian statistics as a robust probe to structural phase transition was also confirmed by comparing the results from differential scanning calorimetry measurements. We further studied the response of the NGC to the modifications in the microstructure on repeated thermal cycling, as well as the variations in the temperature-drive rate, and explained the results using established simplistic models based on the different competing time scales. Our experiments (i) suggest an alternative method to estimate the transformation temperature scales with high accuracy and (ii) establish a connection between the material-specific evolution of microstructure to the statistics of its linear response. Since the method depends on an in-built long-range correlation during transformation, it could be portable to other structural transitions, as well as to materials of different physical origin and size.

Essays on the Modeling and Prediction of Volatility and Higher Moments of Stock Returns (summary section only)

One of the most fundamental and widely accepted ideas in finance is that investors are compensated through higher returns for taking on non-diversifiable risk. Hence the quantification, modeling and prediction of risk have been, and still are one of the most prolific research areas in financial economics. It was recognized early on that there are predictable patterns in the variance of speculative prices. Later research has shown that there may also be systematic variation in the skewness and kurtosis of financial returns. Lacking in the literature so far, is an out-of-sample forecast evaluation of the potential benefits of these new more complicated models with time-varying higher moments. Such an evaluation is the topic of this dissertation. Essay 1 investigates the forecast performance of the GARCH (1,1) model when estimated with 9 different error distributions on Standard and Poor’s 500 Index Future returns. By utilizing the theory of realized variance to construct an appropriate ex post measure of variance from intra-day data it is shown that allowing for a leptokurtic error distribution leads to significant improvements in variance forecasts compared to using the normal distribution. This result holds for daily, weekly as well as monthly forecast horizons. It is also found that allowing for skewness and time variation in the higher moments of the distribution does not further improve forecasts. In Essay 2, by using 20 years of daily Standard and Poor 500 index returns, it is found that density forecasts are much improved by allowing for constant excess kurtosis but not improved by allowing for skewness. By allowing the kurtosis and skewness to be time varying the density forecasts are not further improved but on the contrary made slightly worse. In Essay 3 a new model incorporating conditional variance, skewness and kurtosis based on the Normal Inverse Gaussian (NIG) distribution is proposed. The new model and two previously used NIG models are evaluated by their Value at Risk (VaR) forecasts on a long series of daily Standard and Poor’s 500 returns. The results show that only the new model produces satisfactory VaR forecasts for both 1% and 5% VaR Taken together the results of the thesis show that kurtosis appears not to exhibit predictable time variation, whereas there is found some predictability in the skewness. However, the dynamic properties of the skewness are not completely captured by any of the models.

Response of nonlinear systems to narrow-band excitation

The hardening cubic spring oscillator is studied under narrow-band gaussian excitation. Equivalent linearization leads to multiple steady states. The realizability of the solution is discussed through stochastic stability analysis. Theoretical results are supported by numerical simulation.

A Wavelet based Teager Energy Operator for Spike Detection in Microelectrode Array Recordings

Spike detection in neural recordings is the initial step in the creation of brain machine interfaces. The Teager energy operator (TEO) treats a spike as an increase in the `local' energy and detects this increase. The performance of TEO in detecting action potential spikes suffers due to its sensitivity to the frequency of spikes in the presence of noise which is present in microelectrode array (MEA) recordings. The multiresolution TEO (mTEO) method overcomes this shortcoming of the TEO by tuning the parameter k to an optimal value m so as to match to frequency of the spike. In this paper, we present an algorithm for the mTEO using the multiresolution structure of wavelets along with inbuilt lowpass filtering of the subband signals. The algorithm is efficient and can be implemented for real-time processing of neural signals for spike detection. The performance of the algorithm is tested on a simulated neural signal with 10 spike templates obtained from [14]. The background noise is modeled as a colored Gaussian random process. Using the noise standard deviation and autocorrelation functions obtained from recorded data, background noise was simulated by an autoregressive (AR(5)) filter. The simulations show a spike detection accuracy of 90%and above with less than 5% false positives at an SNR of 2.35 dB as compared to 80% accuracy and 10% false positives reported [6] on simulated neural signals.

Performance analysis of a suprathreshold, stochastic resonance based nonlinear detector

In this paper we introduce a nonlinear detector based on the phenomenon of suprathreshold stochastic resonance (SSR). We first present a model (an array of 1-bit quantizers) that demonstrates the SSR phenomenon. We then use this as a pre-processor to the conventional matched filter. We employ the Neyman-Pearson(NP) detection strategy and compare the performances of the matched filter, the SSR-based detector and the optimal detector. Although the proposed detector is non-optimal, for non-Gaussian noises with heavy tails (leptokurtic) it shows better performance than the matched filter. In situations where the noise is known to be leptokurtic without the availability of the exact knowledge of its distribution, the proposed detector turns out to be a better choice than the matched filter.

Cluster based core vector machine

Core Vector Machine(CVM) is suitable for efficient large-scale pattern classification. In this paper, a method for improving the performance of CVM with Gaussian kernel function irrespective of the orderings of patterns belonging to different classes within the data set is proposed. This method employs a selective sampling based training of CVM using a novel kernel based scalable hierarchical clustering algorithm. Empirical studies made on synthetic and real world data sets show that the proposed strategy performs well on large data sets.

A suprathreshold stochastic resonance based pre-processor for direction-of-arrival estimation

In this paper we propose a nonlinear preprocessor for enhancing the performance of processors used for direction-of-arrival (DOA) estimation in heavy-tailed non-Gaussian noise. The preprocessor based on the phenomenon of suprathreshold stochastic resonance (SSR), provides SNR gain. The preprocessed data is used for DOA estimation by the MUSIC algorithm. Simulation results are presented to show that the SSR preprocessor provides a significant improvement in the performance of MUSIC in heavy-tailed noise at low SNR.

Orbital evolution and orbital phase resolved spectroscopy of the HMXB pulsar 4U 1538-52 with RXTE-PCA and BeppoSAX

We report here results from detailed timing and spectral studies of the high mass X-ray binary pulsar 4U 1538-52 over several binary periods using observations made with the Rossi X-ray Timing Explorer (RXTE) and BeppoSAX satellites. Pulse timing analysis with the 2003 RXTE data over two binary orbits confirms an eccentric orbit of the system. Combining the orbitial parameters determined from this observation with the earlier measurements we did not find any evidence of orbital decay in this X-ray binary. We have carried out orbital phase resolved spectroscopy to measure changes in the spectral parameters with orbital phase, particularly the absorption column density and the iron line flux. The RXTE-PCA spectra in the 3-20 keV energy range were fitted with a power law and a high energy cut-off along with a Gaussian line at similar to 6.4 keV, whereas the BeppoSAX spectra needed only a power law and Gaussian emission line at similar to 6.4keV in the restricted energy range of 0.3-10.0 keV. An absorption along the line of sight was included for both the RXTE and BeppoSAX data. The variation of the free spectral parameters over the binary orbit was investigated and we found that the variation of the column density of absorbing material in the line of sight with orbital phase is in reasonable agreement with a simple model of a spherically symmetric stellar wind from the companion star.