Diffusion Anomaly in Silicalite and VPI-5 from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations on rigid and flexible framework models of silicalite and a rigid framework model of the aluminophosphate VPI-5 for different sorbate diameters are reported. The sorbate-host interactions are modeled in terms of simple atom-atom Lennard-Jones interactions. The results suggest that the diffusion coefficient exhibits an anomaly as gamma approaches unity. The MD results confirm the existence of a linear regime for sorbate diameters significantly smaller than the channel diameter and an anomalous regime observed for sorbate diameters comparable to the channel diameter. The power spectra obtained by Fourier transformation of the velocity autocorrelation function indicate that there is an increase in the intensity of the low-frequency component for the velocity component parallel to the direction of motion for the sorbate diameter in the anomalous regime. The present results suggest that the diffusion anomaly is observed irrespective of (1) the geometry and topology of the pore structure and (2) the nature of the host material. The results are compared with the work of Derouane and co-workers, who have suggested the existence of ''floating molecules'' on the basis of earlier theoretical and computational approaches.

From the Icosahedron to Natural Triangulations of CP(2) and S(2) x S(2)

We present two constructions in this paper: (a) a 10-vertex triangulation CP(10)(2) of the complex projective plane CP(2) as a subcomplex of the join of the standard sphere (S(4)(2)) and the standard real projective plane (RP(6)(2), the decahedron), its automorphism group is A(4); (b) a 12-vertex triangulation (S(2) x S(2))(12) of S(2) x S(2) with automorphism group 2S(5), the Schur double cover of the symmetric group S(5). It is obtained by generalized bistellar moves from a simplicial subdivision of the standard cell structure of S(2) x S(2). Both constructions have surprising and intimate relationships with the icosahedron. It is well known that CP(2) has S(2) x S(2) as a two-fold branched cover; we construct the triangulation CP(10)(2) of CP(2) by presenting a simplicial realization of this covering map S(2) x S(2) -> CP(2). The domain of this simplicial map is a simplicial subdivision of the standard cell structure of S(2) x S(2), different from the triangulation alluded to in (b). This gives a new proof that Kuhnel's CP(9)(2) triangulates CP(2). It is also shown that CP(10)(2) and (S(2) x S(2))(12) induce the standard piecewise linear structure on CP(2) and S(2) x S(2) respectively.

Computational algorithms inspired by biological processes and evolution

In recent times computational algorithms inspired by biological processes and evolution are gaining much popularity for solving science and engineering problems. These algorithms are broadly classified into evolutionary computation and swarm intelligence algorithms, which are derived based on the analogy of natural evolution and biological activities. These include genetic algorithms, genetic programming, differential evolution, particle swarm optimization, ant colony optimization, artificial neural networks, etc. The algorithms being random-search techniques, use some heuristics to guide the search towards optimal solution and speed-up the convergence to obtain the global optimal solutions. The bio-inspired methods have several attractive features and advantages compared to conventional optimization solvers. They also facilitate the advantage of simulation and optimization environment simultaneously to solve hard-to-define (in simple expressions), real-world problems. These biologically inspired methods have provided novel ways of problem-solving for practical problems in traffic routing, networking, games, industry, robotics, economics, mechanical, chemical, electrical, civil, water resources and others fields. This article discusses the key features and development of bio-inspired computational algorithms, and their scope for application in science and engineering fields.

A study on the usefulness of Support Vector Machines for the realtime computational simulation of soft biological organs

Realistic and realtime computational simulation of soft biological organs (e.g., liver, kidney) is necessary when one tries to build a quality surgical simulator that can simulate surgical procedures involving these organs. Since the realistic simulation of these soft biological organs should account for both nonlinear material behavior and large deformation, achieving realistic simulations in realtime using continuum mechanics based numerical techniques necessitates the use of a supercomputer or a high end computer cluster which are costly. Hence there is a need to employ soft computing techniques like Support Vector Machines (SVMs) which can do function approximation, and hence could achieve physically realistic simulations in realtime by making use of just a desktop computer. Present work tries to simulate a pig liver in realtime. Liver is assumed to be homogeneous, isotropic, and hyperelastic. Hyperelastic material constants are taken from the literature. An SVM is employed to achieve realistic simulations in realtime, using just a desktop computer. The code for the SVM is obtained from [1]. The SVM is trained using the dataset generated by performing hyperelastic analyses on the liver geometry, using the commercial finite element software package ANSYS. The methodology followed in the present work closely follows the one followed in [2] except that [2] uses Artificial Neural Networks (ANNs) while the present work uses SVMs to achieve realistic simulations in realtime. Results indicate the speed and accuracy that is obtained by employing the SVM for the targeted realistic and realtime simulation of the liver.

Efficient method of moving shadow detection and vehicle classification

Moving shadow detection and removal from the extracted foreground regions of video frames, aim to limit the risk of misconsideration of moving shadows as a part of moving objects. This operation thus enhances the rate of accuracy in detection and classification of moving objects. With a similar reasoning, the present paper proposes an efficient method for the discrimination of moving object and moving shadow regions in a video sequence, with no human intervention. Also, it requires less computational burden and works effectively under dynamic traffic road conditions on highways (with and without marking lines), street ways (with and without marking lines). Further, we have used scale-invariant feature transform-based features for the classification of moving vehicles (with and without shadow regions), which enhances the effectiveness of the proposed method. The potentiality of the method is tested with various data sets collected from different road traffic scenarios, and its superiority is compared with the existing methods. (C) 2013 Elsevier GmbH. All rights reserved.

On piercing (pseudo)lines and boxes

We say a family of geometric objects C has (l;k)-property if every subfamily C0C of cardinality at most lisk- piercable. In this paper we investigate the existence of g(k;d)such that if any family of objects C in Rd has the (g(k;d);k)-property, then C is k-piercable. Danzer and Gr̈ unbaum showed that g(k;d)is infinite for fami-lies of boxes and translates of centrally symmetric convex hexagons. In this paper we show that any family of pseudo-lines(lines) with (k2+k+ 1;k)-property is k-piercable and extend this result to certain families of objects with discrete intersections. This is the first positive result for arbitrary k for a general family of objects. We also pose a relaxed ver-sion of the above question and show that any family of boxes in Rd with (k2d;k)-property is 2dk- piercable.

Effect of Inlet Jet injection angle on a calandria based reactor-An investigation with numerical analysis of heat and mass transfer for an optimum reactor design

Heat and mass transfer studies in a calandria based reactor is quite complex both due to geometry and due to the complex mixing flow. It is challenging to devise optimum operating conditions with efficient but safe working range for such a complex configuration. Numerical study known to be very effective is taken up for investigation. In the present study a 3D RANS code with turbulence model has been used to compute the flow fields and to get the heat transfer characteristics to understand certain design parameters of engineering importance. The angle of injection and of the coolant liquid has a large effect on the heat transfer within the reactor.

THE LEGACY OF GN RAMACHANDRAN AND THE DEVELOPMENT OF STRUCTURAL BIOLOGY IN INDIA

G.N. Ramachandran is among the founding fathers of structural molecular biology. He made pioneering contributions in computational biology, modelling and what we now call bioinformatics. The triple helical coiled coil structure of collagen proposed by him forms the basis of much of collagen research at the molecular level. The Ramachandran map remains the simplest descriptor and tool for validation of protein structures. He has left his imprint on almost all aspects of biomolecular conformation. His contributions in the area of theoretical crystallography have been outstanding. His legacy has provided inspiration for the further development of structural biology in India. After a pause, computational biology and bioinformatics are in a resurgent phase. One of the two schools established by Ramachandran pioneered the development of macromolecular crystallography, which has now grown into an important component of modern biological research in India. Macromolecular NMR studies in the country are presently gathering momentum. Structural biology in India is now poised to again approach heights of the kind that Ramachandran conquered more than a generation ago.

Monte Carlo simulation of light transport in turbid medium with embedded object-spherical, cylindrical, ellipsoidal, or cuboidal objects embedded within multilayered tissues

Monte Carlo modeling of light transport in multilayered tissue (MCML) is modified to incorporate objects of various shapes (sphere, ellipsoid, cylinder, or cuboid) with a refractive-index mismatched boundary. These geometries would be useful for modeling lymph nodes, tumors, blood vessels, capillaries, bones, the head, and other body parts. Mesh-based Monte Carlo (MMC) has also been used to compare the results from the MCML with embedded objects (MCML-EO). Our simulation assumes a realistic tissue model and can also handle the transmission/reflection at the object-tissue boundary due to the mismatch of the refractive index. Simulation of MCML-EO takes a few seconds, whereas MMC takes nearly an hour for the same geometry and optical properties. Contour plots of fluence distribution from MCML-EO and MMC correlate well. This study assists one to decide on the tool to use for modeling light propagation in biological tissue with objects of regular shapes embedded in it. For irregular inhomogeneity in the model (tissue), MMC has to be used. If the embedded objects (inhomogeneity) are of regular geometry (shapes), then MCML-EO is a better option, as simulations like Raman scattering, fluorescent imaging, and optical coherence tomography are currently possible only with MCML. (C) 2014 Society of Photo-Optical Instrumentation Engineers (SPIE)

Small strong epsilon nets

Let P be a set of n points in R-d. A point x is said to be a centerpoint of P if x is contained in every convex object that contains more than dn/d+1 points of P. We call a point x a strong centerpoint for a family of objects C if x is an element of P is contained in every object C is an element of C that contains more than a constant fraction of points of P. A strong centerpoint does not exist even for halfspaces in R-2. We prove that a strong centerpoint exists for axis-parallel boxes in Rd and give exact bounds. We then extend this to small strong epsilon-nets in the plane. Let epsilon(S)(i) represent the smallest real number in 0, 1] such that there exists an epsilon(S)(i)-net of size i with respect to S. We prove upper and lower bounds for epsilon(S)(i) where S is the family of axis-parallel rectangles, halfspaces and disks. (C) 2014 Elsevier B.V. All rights reserved.

The beta hairpin structure within ribosomal protein S5 mediates interplay between domains II and IV and regulates HCV IRES function

Translation initiation in Hepatitis C Virus (HCV) is mediated by Internal Ribosome Entry Site (IRES), which is independent of cap-structure and uses a limited number of canonical initiation factors. During translation initiation IRES-40S complex formation depends on high affinity interaction of IRES with ribosomal proteins. Earlier, it has been shown that ribosomal protein S5 (RPS5) interacts with HCV IRES. Here, we have extensively characterized the HCV IRES-RPS5 interaction and demonstrated its role in IRES function. Computational modelling and RNA-protein interaction studies demonstrated that the beta hairpin structure within RPS5 is critically required for the binding with domains II and IV. Mutations disrupting IRES-RPS5 interaction drastically reduced the 80S complex formation and the corresponding IRES activity. Computational analysis and UV cross-linking experiments using various IRES-mutants revealed interplay between domains II and IV mediated by RPS5. In addition, present study demonstrated that RPS5 interaction is unique to HCV IRES and is not involved in 40S-3 ` UTR interaction. Further, partial silencing of RPS5 resulted in preferential inhibition of HCV RNA translation. However, global translation was marginally affected by partial silencing of RPS5. Taken together, results provide novel molecular insights into IRES-RPS5 interaction and unravel its functional significance in mediating internal initiation of translation.

The legacy of G. N. Ramachandran and the development of structural biology in India

G. N. Ramachandran is among the founding fathers of structural molecular biology. He made pioneering contributions in computational biology, modelling and what we now call bioinformatics. The triple helical coiled coil structure of collagen proposed by him forms the basis of much of collagen research at the molecular level. The Ramachandran map remains the simplest descriptor and tool for validation of protein structures. He has left his imprint on almost all aspects of biomolecular conformation. His contributions in the area of theoretical crystallography have been outstanding. His legacy has provided inspiration for the further development of structural biology in India. After a pause, computational biology and bioinformatics are in a resurgent phase. One of the two schools established by Ramachandran pioneered the development of macromolecular crystallography, which has now grown into an important component of modern biological research in India. Macromolecular NMR studies in the country are presently gathering momentum. Structural biology in India is now poised to again approach heights of the kind that Ramachandran conquered more than a generation ago.

Efficient simulation and rendering of realistic motion of one-dimensional flexible objects

In gross motion of flexible one-dimensional (1D) objects such as cables, ropes, chains, ribbons and hair, the assumption of constant length is realistic and reasonable. The motion of the object also appears more natural if the motion or disturbance given at one end attenuates along the length of the object. In an earlier work, variational calculus was used to derive natural and length-preserving transformation of planar and spatial curves and implemented for flexible 1D objects discretized with a large number of straight segments. This paper proposes a novel idea to reduce computational effort and enable real-time and realistic simulation of the motion of flexible 1D objects. The key idea is to represent the flexible 1D object as a spline and move the underlying control polygon with much smaller number of segments. To preserve the length of the curve to within a prescribed tolerance as the control polygon is moved, the control polygon is adaptively modified by subdivision and merging. New theoretical results relating the length of the curve and the angle between the adjacent segments of the control polygon are derived for quadratic and cubic splines. Depending on the prescribed tolerance on length error, the theoretical results are used to obtain threshold angles for subdivision and merging. Simulation results for arbitrarily chosen planar and spatial curves whose one end is subjected to generic input motions are provided to illustrate the approach. (C) 2016 Elsevier Ltd. All rights reserved.

Structure and interactions of aspirin-like drugs

A pre-requisite for the elucidation of the mechanism of action of aspirin-like drugs, which are believed to exert their pharmacological effects through the inhibition of prostaglandin biosynthesis, is an understanding of their molecular geometry, the non-covalent interactions they are likely to be involved in, and the geometrical and the electronic consequences of such interactions. This has been sought to be achieved through the x-ray analysis of these drug molecules and their crystalline complexes with other suitable molecules. The results obtained from such studies have been discussed in terms of specific typical examples. For instance, antipyrine can form metal and hydrogen-bonded complexes; phenylbutazone can form ionic complexes with basic molecules. Complex formation is accompanied by characteristic changes in the molecular geometry and the electronic structure in both the cases. The results obtained so far appear to indicate that the important common invariant structural features of the fenamates, deduced from crystal structures, are retained even when complexation takes place.

Efficient algorithms for computing Reeb graphs

The Reeb graph tracks topology changes in level sets of a scalar function and finds applications in scientific visualization and geometric modeling. We describe an algorithm that constructs the Reeb graph of a Morse function defined on a 3-manifold. Our algorithm maintains connected components of the two dimensional levels sets as a dynamic graph and constructs the Reeb graph in O(nlogn+nlogg(loglogg)3) time, where n is the number of triangles in the tetrahedral mesh representing the 3-manifold and g is the maximum genus over all level sets of the function. We extend this algorithm to construct Reeb graphs of d-manifolds in O(nlogn(loglogn)3) time, where n is the number of triangles in the simplicial complex that represents the d-manifold. Our result is a significant improvement over the previously known O(n2) algorithm. Finally, we present experimental results of our implementation and demonstrate that our algorithm for 3-manifolds performs efficiently in practice.