Thermal Fluctuation Spectroscopy of DNA Thermal Denaturation

We have developed the technique of thermal fluctuation spectroscopy to measure the thermal fluctuations in a system. This technique is particularly useful to study the denaturation dynamics of biomolecules like DNA. Here we present a study of the thermal fluctuations during the thermal denaturation (or melting) of double-stranded DNA. We find that the thermal denaturation of heteropolymeric DNA is accompanied by large, non-Gaussian thermal fluctuations. The thermal fluctuations show a two-peak structure as a function of temperature. Calculations of enthalpy exchanged show that the first peak comes from the denaturation of AT rich regions and the second peak from denaturation of GC rich regions. The large fluctuations are almost absent in homopolymeric DNA. We suggest that bubble formation and cooperative opening and closing dynamics of basepairs causes the additional fluctuation at the first peak and a large cooperative transition from a partially molten DNA to a completely denatured state causes the additional fluctuation at the second peak.

Characterization of unitary processes with independent and stationary increments

This is a continuation of the earlier work (Publ. Res. Inst. Math. Sci. 45 (2009) 745-785) to characterize unitary stationary independent increment Gaussian processes. The earlier assumption of uniform continuity is replaced by weak continuity and with technical assumptions on the domain of the generator, unitary equivalence of the process to the solution of an appropriate Hudson-Parthasarathy equation is proved.

Narrowband random excitation of a limit cycle system

The response of the Van der Pol oscillator to stationary narrowband Gaussian excitation is considered. The central frequency of excitation is taken to be in the neighborhood of the system limit cycle frequency. The solution is obtained using a non-Gaussian closure approximation on the probability density function of the response. The validity of the solution is examined with the help of a stochastic stability analysis. Solution based on Stratonovich''s quasistatic averaging technique is also obtained. The comparison of the theoretical solutions with the digital simulations shows that the theoretical estimates are reasonably good.

Detection of edges from projections

In a number of applications of computerized tomography, the ultimate goal is to detect and characterize objects within a cross section. Detection of edges of different contrast regions yields the required information. The problem of detecting edges from projection data is addressed. It is shown that the class of linear edge detection operators used on images can be used for detection of edges directly from projection data. This not only reduces the computational burden but also avoids the difficulties of postprocessing a reconstructed image. This is accomplished by a convolution backprojection operation. For example, with the Marr-Hildreth edge detection operator, the filtering function that is to be used on the projection data is the Radon transform of the Laplacian of the 2-D Gaussian function which is combined with the reconstruction filter. Simulation results showing the efficacy of the proposed method and a comparison with edges detected from the reconstructed image are presented

A Geostatistical Study pf Alumina-Content in Ores of An Indian Iron-Ore Deposit

Mining and blending operations in the high grade iron ore deposit under study are performed to optimize recovery with minimal alumina content while maintaining required levels of other chemical component and a proper mix of ore types. In the present work the regionalisation of alumina in the ores has been studied independently and its effects on global and local recoverable tonnage as well as on alternatives of mining operations have been evaluated. The global tonnage recovery curves for blocks (20m x 20m x 12m) obtained by simulation closely approximated the curves obtained theoretically using a change of support under the discretised gaussian model. Variations in block size up to 80m x 20m x 12m did not affect the recovery as the horizontal dimensions of the blocks are small in relation to the range of the variogram. A comparison of the local tonnage recovery curves obtained through multiple conditional simulations made with that obtained by the method of uniform conditioning of block grades on an estimate of panel 100m x 100m x 12m panel grade reveals comparable results only in panels which have been well conditioned and possesing an ensemble simulation mean close to the ordinary kriged value for the panel. Study of simple alternative sequence of mining on the conditionally simulated deposit shows that concentration of mining operations simultaneously on a single bench enhances the fluctuation in alumina values of ore mined.

Path integral description of polymers using fractional Brownian walks

The statistical properties of fractional Brownian walks are used to construct a path integral representation of the conformations of polymers with different degrees of bond correlation. We specifically derive an expression for the distribution function of the chains’ end‐to‐end distance, and evaluate it by several independent methods, including direct evaluation of the discrete limit of the path integral, decomposition into normal modes, and solution of a partial differential equation. The distribution function is found to be Gaussian in the spatial coordinates of the monomer positions, as in the random walk description of the chain, but the contour variables, which specify the location of the monomer along the chain backbone, now depend on an index h, the degree of correlation of the fractional Brownian walk. The special case of h=1/2 corresponds to the random walk. In constructing the normal mode picture of the chain, we conjecture the existence of a theorem regarding the zeros of the Bessel function.

Solvation dynamics in liquid water. A novel interplay between librational and diffusive modes

A microscopic calculation of the solvation dynamics of an ion in liquid water is presented. The calculated solvation time correlation function shows an ultrafast Gaussian decay which carries about 70%–90% of the strength followed by a biexponential decay with time constants equal to 250 fs and 1 ps. These results are in excellent agreement with the computer simulations of Maroncelli and Fleming and also with the experimental findings of Barbara and Jarzeba. In addition, we find that both the rotational librations and the intermolecular translational vibrational modes of water contribute significantly to the initial Gaussian decay.

Some new results regarding spikes and a heuristic for spike construction

his paper addresses the problem of minimizing the number of columns with superdiagonal nonzeroes (viz., spiked columns) in a square, nonsingular linear system of equations which is to be solved by Gaussian elimination. The exact focus is on a class of min-spike heuristics in which the rows and columns of the coefficient matrix are first permuted to block lower-triangular form. Subsequently, the number of spiked columns in each irreducible block and their heights above the diagonal are minimized heuristically. We show that ifevery column in an irreducible block has exactly two nonzeroes, i.e., is a doubleton, then there is exactly one spiked column. Further, if there is at least one non-doubleton column, there isalways an optimal permutation of rows and columns under whichnone of the doubleton columns are spiked. An analysis of a few benchmark linear programs suggests that singleton and doubleton columns can abound in practice. Hence, it appears that the results of this paper can be practically useful. In the rest of the paper, we develop a polynomial-time min-spike heuristic based on the above results and on a graph-theoretic interpretation of doubleton columns.

E. C. G. Sudarshan and Symmetry in Classical Dynamics, Optics and Quantum Mechanics

We review work initiated and inspired by Sudarshan in relativistic dynamics, beam optics, partial coherence theory, Wigner distribution methods, multimode quantum optical squeezing, and geometric phases. The 1963 No Interaction Theorem using Dirac's instant form and particle World Line Conditions is recalled. Later attempts to overcome this result exploiting constrained Hamiltonian theory, reformulation of the World Line Conditions and extending Dirac's formalism, are reviewed. Dirac's front form leads to a formulation of Fourier Optics for the Maxwell field, determining the actions of First Order Systems (corresponding to matrices of Sp(2,R) and Sp(4,R)) on polarization in a consistent manner. These groups also help characterize properties and propagation of partially coherent Gaussian Schell Model beams, leading to invariant quality parameters and the new Twist phase. The higher dimensional groups Sp(2n,R) appear in the theory of Wigner distributions and in quantum optics. Elegant criteria for a Gaussian phase space function to be a Wigner distribution, expressions for multimode uncertainty principles and squeezing are described. In geometric phase theory we highlight the use of invariance properties that lead to a kinematical formulation and the important role of Bargmann invariants. Special features of these phases arising from unitary Lie group representations, and a new formulation based on the idea of Null Phase Curves, are presented.

Symmetry Properties of the Large-Deviation Function of the Velocity of a Self-Propelled Polar Particle

A geometrically polar granular rod confined in 2D geometry, subjected to a sinusoidal vertical oscillation, undergoes noisy self-propulsion in a direction determined by its polarity. When surrounded by a medium of crystalline spherical beads, it displays substantial negative fluctuations in its velocity. We find that the large-deviation function (LDF) for the normalized velocity is strongly non-Gaussian with a kink at zero velocity, and that the antisymmetric part of the LDF is linear, resembling the fluctuation relation known for entropy production, even when the velocity distribution is clearly non-Gaussian. We extract an analogue of the phase-space contraction rate and find that it compares well with an independent estimate based on the persistence of forward and reverse velocities.

Optimal pulse spacing for dynamical decoupling in the presence of a purely dephasing spin bath

Maintaining quantum coherence is a crucial requirement for quantum computation; hence protecting quantum systems against their irreversible corruption due to environmental noise is an important open problem. Dynamical decoupling (DD) is an effective method for reducing decoherence with a low control overhead. It also plays an important role in quantum metrology, where, for instance, it is employed in multiparameter estimation. While a sequence of equidistant control pulses the Carr-Purcell-Meiboom-Gill (CPMG) sequence] has been ubiquitously used for decoupling, Uhrig recently proposed that a nonequidistant pulse sequence the Uhrig dynamic decoupling (UDD) sequence] may enhance DD performance, especially for systems where the spectral density of the environment has a sharp frequency cutoff. On the other hand, equidistant sequences outperform UDD for soft cutoffs. The relative advantage provided by UDD for intermediate regimes is not clear. In this paper, we analyze the relative DD performance in this regime experimentally, using solid-state nuclear magnetic resonance. Our system qubits are C-13 nuclear spins and the environment consists of a H-1 nuclear spin bath whose spectral density is close to a normal (Gaussian) distribution. We find that in the presence of such a bath, the CPMG sequence outperforms the UDD sequence. An analogy between dynamical decoupling and interference effects in optics provides an intuitive explanation as to why the CPMG sequence performs better than any nonequidistant DD sequence in the presence of this kind of environmental noise.

Data compression by linear prediction for storage and transmission of EEG signals

The EEG time series has been subjected to various formalisms of analysis to extract meaningful information regarding the underlying neural events. In this paper the linear prediction (LP) method has been used for analysis and presentation of spectral array data for the better visualisation of background EEG activity. It has also been used for signal generation, efficient data storage and transmission of EEG. The LP method is compared with the standard Fourier method of compressed spectral array (CSA) of the multichannel EEG data. The autocorrelation autoregressive (AR) technique is used for obtaining the LP coefficients with a model order of 15. While the Fourier method reduces the data only by half, the LP method just requires the storage of signal variance and LP coefficients. The signal generated using white Gaussian noise as the input to the LP filter has a high correlation coefficient of 0.97 with that of original signal, thus making LP as a useful tool for storage and transmission of EEG. The biological significance of Fourier method and the LP method in respect to the microstructure of neuronal events in the generation of EEG is discussed.

Effect of surface roughness on diffusion-limited charge transfer

A theory is developed for diffusion-limited charge transfer on a non-fractally rough electrode. The perturbation expressions are obtained for concentration, current density and measured diffusion-limited current for arbitrary one- and two-dimensional surface profiles. The random surface model is employed for a rough electrode\electrolyte interface. In this model the gross geometrical property of an electrochemically active rough surface - the surface structure factor-is related to the average electrode current, current density and concentration. Under short and long time regimes, various morphological features of the rough electrodes, i.e. excess area (related to roughness slope), curvature, correlation length, etc. are related to the (average) current transients. A two-point Pade approximant is used to develop an all time average current expression in terms of partial morphological features of the rough surface. The inverse problem of predicting the surface structure factor from the observed transients is also described. Finally, the effect of surface roughness is studied for specific surface statistics, namely a Gaussian correlation function. It is shown how the surface roughness enhances the overall diffusion-limited charge transfer current.

Effects of Ultrafast Solvation on the Rate of Adiabatic Outer-Sphere Electron Transfer Reactions

Recent studies have demonstrated that solvation dynamics in many common dipolar liquids contain an initial, ultrafast Gaussian component which may contribute even more than 60% to the total solvation energy. It is also known that adiabatic electron transfer reactions often probe the high-frequency components of the relevant solvent friction (Hynes, J. T. J. Phys. Chem. 1986, 90, 3701). In this paper, we present a theoretical study of the effects of the ultrafast solvent polar modes on the adiabatic electron transfer reactions by using the formalism of Hynes. Calculations have been carried out for a model system and also for water and acetonitrile. It is found that, in general, the ultrafast modes can greatly enhance the rate of electron transfer, even by more than an order of magnitude, over the rate obtained by using only the slow overdamped modes usually considered. For water, this acceleration of the rate can be attributed to the high-frequency intermolecular vibrational and librational modes. For a weakly adiabatic reaction, the rate is virtually indistinguishable from the rate predicted by the Marcus transition state theory. Another important result is that even in this case of ultrafast underdamped solvation, energy diffusion appears to be efficient so that electron transfer reaction in water is controlled essentially by the barrier crossing dynamics. This is because the reactant well frequency is-directly proportional to the rate of the initial Gaussian decay of the solvation time correlation function. As a result, the value of the friction at the reactant well frequency rarely falls below the value required for the Kramers turnover except when the polarizability of the water molecules may be neglected. On the other hand, in acetonitrile, the rate of electron transfer reaction is found to be controlled by the energy diffusion dynamics, although a significant contribution to the rate comes also from the barrier crossing rate. Therefore, the present study calls for a need to understand the relaxation of the high-frequency modes in dipolar liquids.

Measurement of complex optical flow with use of an augmented generalized gradient scheme

We present a method for measuring the local velocities and first-order variations in velocities in a timevarying image. The scheme is an extension of the generalized gradient model that encompasses the local variation of velocity within a local patch of the image. Motion within a patch is analyzed in parallel by 42 different spatiotemporal filters derived from 6 linearly independent spatiotemporal kernels. No constraints are imposed on the image structure, and there is no need for smoothness constraints on the velocity field. The aperture problem does not arise so long as there is some two-dimensional structure in the patch being analyzed. Among the advantages of the scheme is that there is no requirement to calculate second or higher derivatives of the image function. This makes the scheme robust in the presence of noise. The spatiotemporal kernels are of simple form, involving Gaussian functions, and are biologically plausible receptive fields. The validity of the scheme is demonstrated by application to both synthetic and real video images sequences and by direct comparison with another recently published scheme Biol. Cybern. 63, 185 (1990)] for the measurement of complex optical flow.