Theory of polymer breaking under tension

We consider the breaking of a polymer molecule which is fixed at one end and is acted upon by a force at the other. The polymer is assumed to be a linear chain joined together by bonds which satisfy the Morse potential. The applied force is found to modify the Morse potential so that the minimum becomes metastable. Breaking is just the decay of this metastable bond, by causing it to go over the barrier. Increasing the force causes the potential to become more and more distorted and eventually leads to the disappearance of the barrier. The limiting force at which the barrier disappears is D(e)a/2,D-e with a the parameters characterizing the Morse potential. The rate of breaking is first calculated using multidimensional quantum transition state theory. We use the harmonic approximation to account for vibrations of all the units. It includes tunneling contributions to the rate, but is valid only above a certain critical temperature. It is possible to get an analytical expression for the rate of breaking. We have calculated the rate of breaking for a model, which mimics polyethylene. First we calculate the rate of breaking of a single bond, without worrying about the other bonds. Inclusion of other bonds under the harmonic approximation is found to lower this rate by at the most one order of magnitude. Quantum effects are found to increase the rate of breaking and are significant only at temperatures less than 150 K. At 300 K, the calculations predict a bond in polyethylene to have a lifetime of only seconds at a force which is only half the limiting force. Calculations were also done using the Lennard-Jones potential. The results for Lennard-Jones and Morse potentials were rather different, due to the different long-range behaviors of the two potentials. A calculation including friction was carried out, at the classical level, by assuming that each atom of the chain is coupled to its own collection of harmonic oscillators. Comparison of the results with the simulations of Oliveira and Taylor [J. Chem. Phys. 101, 10 118 (1994)] showed the rate to be two to three orders of magnitude higher. As a possible explanation of discrepancy, we consider the translational motion of the ends of the broken chains. Using a continuum approximation for the chain, we find that in the absence of friction, the rate of the process can be limited by the rate at which the two broken ends separate from one another and the lowering of the rate is at the most a factor of 2, for the parameters used in the simulation (for polyethylene). In the presence of friction, we find that the rate can be lowered by one to two orders of magnitude, making our results to be in reasonable agreement with the simulations.

Oneand twodimensional single quantumand multiplequantumNMR methodologies: tools for chiral analyses

It is well known that enantiomers cannot be distinguished by NMR spectroscopy unless diastereomorphic interactions are imposed. Several chiral aligning media have therefore been reported for their visualization, although extensive studies are carried out using the liquid crystal made of polypeptide poly-γ-benzyl-L-glutamate (PBLG) in organic solvent. In PBLG medium the spin systems are weakly coupled and the first order analyses of the spectra are generally possible. But due to large number of pair wise interactions of nuclear spins resulting in many degenerate transitions the 1H NMR spectra are not only complex but also broad and featureless, in addition to an indistinguishable overlap of the spectra of enantiomers. This enormous loss of resolution severely hinders the analyses of proton spectra, even for spin systems with 5–6 interacting protons, thereby restricting itsroutine application. In this review we discuss our recently developed several one and multidimensional NMR experiments to circumvent these difficulties taking specific examples of the molecules containing a single chiral centre.

Development of novel multiple quantummethodologies for the analyses of complex protonNMR spectra of scalar coupled spins

One of the significant advancements in Nuclear Magnetic Resonance spectroscopy (NMR) in combating the problem of spectral complexity for deriving the structure and conformational information is the incorporation of additional dimension and to spread the information content in a two dimensional space. This approach together with the manipulation of the dynamics of nuclear spins permitted the designing of appropriate pulse sequences leading to the evolution of diverse multidimensional NMR experiments. The desired spectral information can now be extracted in a simplified and an orchestrated manner. The indirect detection of multiple quantum (MQ) NMR frequencies is a step in this direction. The MQ technique has been extensively used in the study of molecules aligned in liquid crystalline media to reduce spectral complexity and to determine molecular geometries. Unlike in dipolar coupled systems, the size of the network of scalar coupled spins is not big in isotropic solutions and the MQ 1H detection is not routinely employed,although there are specific examples of spin topology filtering. In this brief review, we discuss our recent studies on the development and application of multiple quantum correlation and resolved techniques for the analyses of proton NMR spectra of scalar coupled spins.

On the solution of bivariate population balance equations for aggregation: X-discretization of space for expansion and contraction of computational domain

A new structured discretization of 2D space, named X-discretization, is proposed to solve bivariate population balance equations using the framework of minimal internal consistency of discretization of Chakraborty and Kumar [2007, A new framework for solution of multidimensional population balance equations. Chem. Eng. Sci. 62, 4112-4125] for breakup and aggregation of particles. The 2D space of particle constituents (internal attributes) is discretized into bins by using arbitrarily spaced constant composition radial lines and constant mass lines of slope -1. The quadrilaterals are triangulated by using straight lines pointing towards the mean composition line. The monotonicity of the new discretization makes is quite easy to implement, like a rectangular grid but with significantly reduced numerical dispersion. We use the new discretization of space to automate the expansion and contraction of the computational domain for the aggregation process, corresponding to the formation of larger particles and the disappearance of smaller particles by adding and removing the constant mass lines at the boundaries. The results show that the predictions of particle size distribution on fixed X-grid are in better agreement with the analytical solution than those obtained with the earlier techniques. The simulations carried out with expansion and/or contraction of the computational domain as population evolves show that the proposed strategy of evolving the computational domain with the aggregation process brings down the computational effort quite substantially; larger the extent of evolution, greater is the reduction in computational effort. (C) 2011 Elsevier Ltd. All rights reserved.

Five-stage Milstein methods for SDEs

In this paper, we consider the problem of computing numerical solutions for Ito stochastic differential equations (SDEs). The five-stage Milstein (FSM) methods are constructed for solving SDEs driven by an m-dimensional Wiener process. The FSM methods are fully explicit methods. It is proved that the FSM methods are convergent with strong order 1 for SDEs driven by an m-dimensional Wiener process. The analysis of stability (with multidimensional Wiener process) shows that the mean-square stable regions of the FSM methods are unbounded. The analysis of stability shows that the mean-square stable regions of the methods proposed in this paper are larger than the Milstein method and three-stage Milstein methods.

Structural Transformations, Composition Anomalies and a Dramatic Collapse of Linear Polymer Chains in Dilute Ethanol-Water Mixtures

Water-ethanol mixtures exhibit many interesting anomalies, such as negative excess partial molar volume of ethanol, excess sound absorption coefficient at low concentrations, and positive deviation from Raoult's law for vapor pressure, to mention a few. These anomalies have been attributed to different, often contradictory origins, but a quantitative understanding is still lacking. We show by computer simulation and theoretical analyses that these anomalies arise from the sudden emergence of a bicontinuous phase that occurs at a relatively low ethanol concentration of x(eth) approximate to 0.06-0.10 (that amounts to a volume fraction of 0.17-0.26, which is a significant range!). The bicontinuous phase is formed by aggregation of ethanol molecules, resulting in a weak phase transition whose nature is elucidated. We find that the microheterogeneous structure of the mixture gives rise to a pronounced nonmonotonic composition dependence of local compressibility and nonmonotonic dependence in the peak value of the radial distribution function of ethyl groups. A multidimensional free energy surface of pair association is shown to provide a molecular explanation of the known negative excess partial volume of ethanol in terms of parallel orientation and hence better packing of the ethyl groups in the mixture due to hydrophobic interactions. The energy distribution of the ethanol molecules indicates additional energy decay channels that explain the excess sound attenuation coefficient in aqueous alcohol mixtures. We studied the dependence of the solvation of a linear polymer chain on the composition of the water-ethanol solvent. We find that there is a sudden collapse of the polymer at x(eth) approximate to 0.05-a phenomenon which we attribute to the formation of the microheterogeneous structures in the binary mixture at low ethanol concentrations. Together with recent single molecule pulling experiments, these results provide new insight into the behavior of polymer chain and foreign solutes, such as enzymes, in aqueous binary mixtures.

Thermal instability and the feedback regulation of hot haloes in clusters, groups and galaxies

We present global multidimensional numerical simulations of the plasma that pervades the dark matter haloes of clusters, groups and massive galaxies (the intracluster medium; ICM). Observations of clusters and groups imply that such haloes are roughly in global thermal equilibrium, with heating balancing cooling when averaged over sufficiently long time- and length-scales; the ICM is, however, very likely to be locally thermally unstable. Using simple observationally motivated heating prescriptions, we show that local thermal instability (TI) can produce a multiphase medium with similar to 104 K cold filaments condensing out of the hot ICM only when the ratio of the TI time-scale in the hot plasma (tTI) to the free-fall time-scale (tff) satisfies tTI/tff? 10. This criterion quantitatively explains why cold gas and star formation are preferentially observed in low-entropy clusters and groups. In addition, the interplay among heating, cooling and TI reduces the net cooling rate and the mass accretion rate at small radii by factors of similar to 100 relative to cooling-flow models. This dramatic reduction is in line with observations. The feedback efficiency required to prevent a cooling flow is similar to 10-3 for clusters and decreases for lower mass haloes; supernova heating may be energetically sufficient to balance cooling in galactic haloes. We further argue that the ICM self-adjusts so that tTI/tff? 10 at all radii. When this criterion is not satisfied, cold filaments condense out of the hot phase and reduce the density of the ICM. These cold filaments can power the black hole and/or stellar feedback required for global thermal balance, which drives tTI/tff? 10. In comparison to clusters, groups have central cores with lower densities and larger radii. This can account for the deviations from self-similarity in the X-ray luminositytemperature () relation. The high-velocity clouds observed in the Galactic halo can be due to local TI producing multiphase gas close to the virial radius if the density of the hot plasma in the Galactic halo is >rsim 10-5 cm-3 at large radii.

Measurement and Applications of Long-Range Heteronuclear Scalar Couplings: Recent Experimental and Theoretical Developments

The use of long-range heteronuclear couplings, in association with 1H1H scalar couplings and NOE restraints, has acquired growing importance for the determination of the relative stereochemistry, and structural and conformational information of organic and biological molecules. However, the routine use of such couplings is hindered by the inherent difficulties in their measurement. Prior to the advancement in experimental techniques, both long-range homo- and heteronuclear scalar couplings were not easily accessible, especially for very large molecules. The development of a large number of multidimensional NMR experimental methodologies has alleviated the complications associated with the measurement of couplings of smaller strengths. Subsequent application of these methods and the utilization of determined J-couplings for structure calculations have revolutionized this area of research. Problems in organic, inorganic and biophysical chemistry have also been solved by utilizing the short- and long-range heteronuclear couplings. In this minireview, we discuss the advantages and limitations of a number of experimental techniques reported in recent times for the measurement of long-range heteronuclear couplings and a few selected applications of such couplings. This includes the study of medium- to larger-sized molecules in a variety of applications, especially in the study of hydrogen bonding in biological systems. The utilization of these couplings in conjunction with theoretical calculations to arrive at conclusions on the hyperconjugation, configurational analysis and the effect of the electronegativity of the substituents is also discussed.

Partition function for hindered, one-and multi-dimensional rotors

Writing the hindered rotor (hr) partition function as the trace of (rho) over cap = e(-beta(H) over cap hr), we approximate it by the sum of contributions from a set of points in position space. The contribution of the density matrix from each point is approximated by performing a local harmonic expansion around it. The highlight of this method is that it can be easily extended to multidimensional systems. Local harmonic expansion leads to a breakdown of the method a low temperatures. In order to calculate the partition function at low temperatures, we suggest a matrix multiplication procedure. The results obtained using these methods closely agree with the exact partition function at all temperature ranges. Our method bypasses the evaluation of eigenvalues and eigenfunctions and evaluates the density matrix for internal rotation directly. We also suggest a procedure to account for the antisymmetry of the total wavefunction in the same. (C) 2012 Elsevier B.V. All rights reserved.

Directional Diffusion Regulator (DDR) for some numerical solvers of hyperbolic conservation laws

A computational tool called ``Directional Diffusion Regulator (DDR)'' is proposed to bring forth real multidimensional physics into the upwind discretization in some numerical schemes of hyperbolic conservation laws. The direction based regulator when used with dimension splitting solvers, is set to moderate the excess multidimensional diffusion and hence cause genuine multidimensional upwinding like effect. The basic idea of this regulator driven method is to retain a full upwind scheme across local discontinuities, with the upwind bias decreasing smoothly to a minimum in the farthest direction. The discontinuous solutions are quantified as gradients and the regulator parameter across a typical finite volume interface or a finite difference interpolation point is formulated based on fractional local maximum gradient in any of the weak solution flow variables (say density, pressure, temperature, Mach number or even wave velocity etc.). DDR is applied to both the non-convective as well as whole unsplit dissipative flux terms of some numerical schemes, mainly of Local Lax-Friedrichs, to solve some benchmark problems describing inviscid compressible flow, shallow water dynamics and magneto-hydrodynamics. The first order solutions consistently improved depending on the extent of grid non-alignment to discontinuities, with the major influence due to regulation of non-convective diffusion. The application is also experimented on schemes such as Roe, Jameson-Schmidt-Turkel and some second order accurate methods. The consistent improvement in accuracy either at moderate or marked levels, for a variety of problems and with increasing grid size, reasonably indicate a scope for DDR as a regular tool to impart genuine multidimensional upwinding effect in a simpler framework. (C) 2012 Elsevier Inc. All rights reserved.

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.

Loop transformations: convexity, pruning, and optimization

High-level loop transformations are a key instrument in mapping computational kernels to effectively exploit the resources in modern processor architectures. Nevertheless, selecting required compositions of loop transformations to achieve this remains a significantly challenging task; current compilers may be off by orders of magnitude in performance compared to hand-optimized programs. To address this fundamental challenge, we first present a convex characterization of all distinct, semantics-preserving, multidimensional affine transformations. We then bring together algebraic, algorithmic, and performance analysis results to design a tractable optimization algorithm over this highly expressive space. Our framework has been implemented and validated experimentally on a representative set of benchmarks running on state-of-the-art multi-core platforms.

Quadrature approximation properties of the spiral-phase quadrature transform

The notion of the 1-D analytic signal is well understood and has found many applications. At the heart of the analytic signal concept is the Hilbert transform. The problem in extending the concept of analytic signal to higher dimensions is that there is no unique multidimensional definition of the Hilbert transform. Also, the notion of analyticity is not so well under stood in higher dimensions. Of the several 2-D extensions of the Hilbert transform, the spiral-phase quadrature transform or the Riesz transform seems to be the natural extension and has attracted a lot of attention mainly due to its isotropic properties. From the Riesz transform, Larkin et al. constructed a vortex operator, which approximates the quadratures based on asymptotic stationary-phase analysis. In this paper, we show an alternative proof for the quadrature approximation property by invoking the quasi-eigenfunction property of linear, shift-invariant systems. We show that the vortex operator comes up as a natural consequence of applying this property. We also characterize the quadrature approximation error in terms of its energy as well as the peak spatial-domain error. Such results are available for 1-D signals, but their counter part for 2-D signals have not been provided. We also provide simulation results to supplement the analytical calculations.

IceCube: efficient targeted mining in data cubes

We address the problem of mining targeted association rules over multidimensional market-basket data. Here, each transaction has, in addition to the set of purchased items, ancillary dimension attributes associated with it. Based on these dimensions, transactions can be visualized as distributed over cells of an n-dimensional cube. In this framework, a targeted association rule is of the form {X -> Y} R, where R is a convex region in the cube and X. Y is a traditional association rule within region R. We first describe the TOARM algorithm, based on classical techniques, for identifying targeted association rules. Then, we discuss the concepts of bottom-up aggregation and cubing, leading to the CellUnion technique. This approach is further extended, using notions of cube-count interleaving and credit-based pruning, to derive the IceCube algorithm. Our experiments demonstrate that IceCube consistently provides the best execution time performance, especially for large and complex data cubes.

Measurement of in-plane and out-of-plane displacements and strains using digital holographic moire

The paper proposes a non-destructive method for simultaneous measurement of in-plane and out-of-plane displacements and strains undergone by a deformed specimen from a single moire fringe pattern obtained on the specimen in a dual beam digital holographic interferometry setup. The moire fringe pattern encodes multiple interference phases which carry the information on multidimensional deformation. The interference field is segmented in each column and is modeled as multicomponent quadratic/cubic frequency-modulated signal in each segment. Subsequently, the product form of modified cubic phase function is used for accurate estimation of phase parameters. The estimated phase parameters are further utilized for direct estimation of the unwrapped interference phases and phase derivatives. The simulation and experimental results are provided to validate the effectiveness of the proposed method.