Parallel Computation of 3D Morse-Smale Complexes

The Morse-Smale complex is a topological structure that captures the behavior of the gradient of a scalar function on a manifold. This paper discusses scalable techniques to compute the Morse-Smale complex of scalar functions defined on large three-dimensional structured grids. Computing the Morse-Smale complex of three-dimensional domains is challenging as compared to two-dimensional domains because of the non-trivial structure introduced by the two types of saddle criticalities. We present a parallel shared-memory algorithm to compute the Morse-Smale complex based on Forman's discrete Morse theory. The algorithm achieves scalability via synergistic use of the CPU and the GPU. We first prove that the discrete gradient on the domain can be computed independently for each cell and hence can be implemented on the GPU. Second, we describe a two-step graph traversal algorithm to compute the 1-saddle-2-saddle connections efficiently and in parallel on the CPU. Simultaneously, the extremasaddle connections are computed using a tree traversal algorithm on the GPU.

Cascaded walks in protein sequence space: use of artificial sequences in remote homology detection between natural proteins

Over the past two decades, many ingenious efforts have been made in protein remote homology detection. Because homologous proteins often diversify extensively in sequence, it is challenging to demonstrate such relatedness through entirely sequence-driven searches. Here, we describe a computational method for the generation of `protein-like' sequences that serves to bridge gaps in protein sequence space. Sequence profile information, as embodied in a position-specific scoring matrix of multiply aligned sequences of bona fide family members, serves as the starting point in this algorithm. The observed amino acid propensity and the selection of a random number dictate the selection of a residue for each position in the sequence. In a systematic manner, and by applying a `roulette-wheel' selection approach at each position, we generate parent family-like sequences and thus facilitate an enlargement of sequence space around the family. When generated for a large number of families, we demonstrate that they expand the utility of natural intermediately related sequences in linking distant proteins. In 91% of the assessed examples, inclusion of designed sequences improved fold coverage by 5-10% over searches made in their absence. Furthermore, with several examples from proteins adopting folds such as TIM, globin, lipocalin and others, we demonstrate that the success of including designed sequences in a database positively sensitized methods such as PSI-BLAST and Cascade PSI-BLAST and is a promising opportunity for enormously improved remote homology recognition using sequence information alone.

Natural Convection in a Vertical Channel with Arrays of Flush-Mounted Heaters on Opposite Conductive Walls

Three-dimensional numerical study of natural convection in a vertical channel with flush-mounted discrete heaters on opposite conductive substrate walls is carried out in the present work. Detailed flow and heat transfer characteristics are presented for various Grashof numbers. The heat transfer effects on one wall by the presence of heaters on its opposite wall is examined. It is found that heat transfer rates on one wall are increased by the presence of heaters on its opposite wall. The thermal boundary layers on the opposite walls complement each other for enhanced heat transfer. The effects of spacing between the heated walls, spacings between heaters and substrate conductivity on flow and heat transfer are examined. Existence of optimum spacings between the heated walls for maximum heat transfer and mass flow are observed. It is found that the heat transfer and fluid flow do not follow the same optimum spacings. Mass flow rate reaches maximum value at a wall spacing greater than the spacing for maximum heat transfer. This is because the interaction of thermal boundary layers on individual walls ceases at a lower spacing before the velocity boundary layers separate each other. It is found that increased spacings between heaters reduce individual heater temperatures provided the heaters close to exit on both substrates avail sufficient substrate potions on the exit side. Insufficient substrate portions between the exit heaters and the exit cause abnormal local temperature rise in the exit heaters which are the hottest ones among all the heaters. Optimal heater spacings exist for minimum hottest heater temperature rise. Correlations are presented for dimensionless mass flow rate, temperature maximum, and average Nusselt number.

Synthetic studies toward geranylated PPAP natural products oblongifolin A, oblongifolin D, and enervosanone

A concise approach of general utility toward mono- and di-geranylated PPAP frameworks employing `reconstructive aldol cyclization' as the key step is delineated. (C) 2012 Elsevier Ltd. All rights reserved.

Enantioselective syntheses of bioactive epoxyquinone natural products (+)-harveynone and (-)-asperpentyn

Stereo- and enantioselective syntheses of (+)-harveynone and (-)-asperpentyn are reported. (C) 2012 Elsevier Ltd. All rights reserved.

Singlet-state creation and universal quantum computation in NMR using a genetic algorithm

The experimental implementation of a quantum algorithm requires the decomposition of unitary operators. Here we treat unitary-operator decomposition as an optimization problem, and use a genetic algorithm-a global-optimization method inspired by nature's evolutionary process-for operator decomposition. We apply this method to NMR quantum information processing, and find a probabilistic way of performing universal quantum computation using global hard pulses. We also demonstrate the efficient creation of the singlet state (a special type of Bell state) directly from thermal equilibrium, using an optimum sequence of pulses. © 2012 American Physical Society.

Singlet-state creation and universal quantum computation in NMR using a genetic algorithm

The experimental implementation of a quantum algorithm requires the decomposition of unitary operators. Here we treat unitary-operator decomposition as an optimization problem, and use a genetic algorithm-a global-optimization method inspired by nature's evolutionary process-for operator decomposition. We apply this method to NMR quantum information processing, and find a probabilistic way of performing universal quantum computation using global hard pulses. We also demonstrate the efficient creation of the singlet state (a special type of Bell state) directly from thermal equilibrium, using an optimum sequence of pulses.

Model synthetic studies toward jiadifenin and majucin type seco-prezizaane natural products via a stereo- and enantioselective approach

A concise stereo- and enantioselective approach to seco-prezizaane sesquiterpenoids, leading to the acquisition of two bicyclic fragments, is delineated. (c) 2012 Elsevier Ltd. All rights reserved.

Study of conjugate natural convection between vertical coaxial rectangular cylinders

Laminar natural convection between two coaxial vertical rectangular cylinders is numerically studied in this work. The outer cylinder is connected with vertical rectangular inlet and outlet pipes. The inner cylinder dissipates volumetric heat. The fluid flow and heat transfer characteristics between the cylinders are analyzed in detail for various Grashof numbers. The heat transfer rates on the individual faces of the inner cylinder are reported. The bottom face of the inner cylinder is found to associate with much higher heat rates than those of the other faces. The average Nusselt number on bottom face is more than 2.5 times of the Nusselt number averaged on all the faces. At a given elevation, local Nusselt number on the inner cylinder faces increases towards cylinder edges. The effect of thermal condition of the walls of outer cylinder, inlet and outlet on the natural convection is analyzed. The thermal condition shows strong qualitative and quantitative impact on the fluid flow and heat transfer. The variation of induced flow rate, dimensionless maximum temperature and average Nusselt numbers with Grashof number is studied. Correlations for dimensionless buoyancy-induced mass flow rate and temperature maximum are presented. (c) 2012 Elsevier Ltd. All rights reserved.

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.

Use of Fick's law and Maxwell–Stefan equations in computation of multicomponent diffusion

This article does not have an abstract.

A finite population model of mlecular evolution: theory and computation

This article is concerned with the evolution of haploid organisms that reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed the quasispecies model in an attempt to understand such an evolutionary process. Their work has impacted antiviral treatment and vaccine design strategies. Yet, predictions of the quasispecies model are at best viewed as a guideline, primarily because it assumes an infinite population size, whereas realistic population sizes can be quite small. In this paper we consider a population genetics-based model aimed at understanding the evolution of such organisms with finite population sizes and present a rigorous study of the convergence and computational issues that arise therein. Our first result is structural and shows that, at any time during the evolution, as the population size tends to infinity, the distribution of genomes predicted by our model converges to that predicted by the quasispecies model. This justifies the continued use of the quasispecies model to derive guidelines for intervention. While the stationary state in the quasispecies model is readily obtained, due to the explosion of the state space in our model, exact computations are prohibitive. Our second set of results are computational in nature and address this issue. We derive conditions on the parameters of evolution under which our stochastic model mixes rapidly. Further, for a class of widely used fitness landscapes we give a fast deterministic algorithm which computes the stationary distribution of our model. These computational tools are expected to serve as a framework for the modeling of strategies for the deployment of mutagenic drugs.

Natural convection in a series of thermally interacting cavities with wall mounted heaters

Laminar natural convection in a series of thermally interacting cavities is numerically studied. Each cavity consists of a conducting bottom wall with a surface mounted heater. The side walls of the cavities are isothermally cooled. Each cavity thermally interacts with its adjacent cavities through the conducting walls. Flow and heat transfer characteristics are studied in detail for various Rayleigh numbers. The convection characteristics in multiple cavities are compared with those in single independent cavity. The thermal interaction between the cavities results in lower temperatures compared with those in independent cavities. While heat is rejected into the adjacent upper cavity through some portion of the conducting wall, heat is received from the adjacent cavity through the remaining portion of the wall. The influence of substrate conductivity on heat exchange between adjacent cavities are examined. Substrate conductivity shows strong effect on temperature distribution. When cooling at both vertical sides is changed to one side cooling, the heat transfer characteristics are changed drastically and many interesting flow features are observed. Effects of cavity aspect ratio is studied and higher heat transfer rates are observed at higher aspect ratios. Correlations for dimensionless temperature maximum and average Nusselt number are presented in terms of Rayleigh number.

Stability of Adult Emergence and Activity/Rest Rhythms in Fruit Flies Drosophila melanogaster under Semi-Natural Condition

Here we report the results of a study aimed at examining stability of adult emergence and activity/rest rhythms under seminatural conditions (henceforth SN), in four large outbred fruit fly Drosophila melanogaster populations, selected for emergence in a narrow window of time under laboratory (henceforth LAB) light/dark (LD) cycles. When assessed under LAB, selected flies display enhanced stability in terms of higher amplitude, synchrony and accuracy in emergence and activity rhythms compared to controls. The present study was conducted to assess whether such differences in stability between selected and control populations, persist under SN where several gradually changing time-cues are present in their strongest form. The study revealed that under SN, emergence waveform of selected flies was modified, with even more enhanced peak and narrower gate-width compared to those observed in the LAB and compared to control populations in SN. Furthermore, flies from selected populations continued to exhibit enhanced synchrony and accuracy in their emergence and activity rhythms under SN compared to controls. Further analysis of zeitgeber effects revealed that enhanced stability in the rhythmicity of selected flies under SN was primarily due to increased sensitivity to light because emergence and activity rhythms of selected flies were as stable as controls under temperature cycles. These results thus suggest that stability of circadian rhythms in fruit flies D. melanogaster, which evolved as a consequence of selection for emergence in a narrow window of time under weak zeitgeber condition of LAB, persists robustly in the face of day-to-day variations in cycling environmental factors of nature.

Secure Computation in a Bidirectional Relay

Bidirectional relaying, where a relay helps two user nodes to exchange equal length binary messages, has been an active area of recent research. A popular strategy involves a modified Gaussian MAC, where the relay decodes the XOR of the two messages using the naturally-occurring sum of symbols simultaneously transmitted by user nodes. In this work, we consider the Gaussian MAC in bidirectional relaying with an additional secrecy constraint for protection against a honest but curious relay. The constraint is that, while the relay should decode the XOR, it should be fully ignorant of the individual messages of the users. We exploit the symbol addition that occurs in a Gaussian MAC to design explicit strategies that achieve perfect independence between the received symbols and individual transmitted messages. Our results actually hold for a more general scenario where the messages at the two user nodes come from a finite Abelian group G, and the relay must decode the sum within G of the two messages. We provide a lattice coding strategy and study optimal rate versus average power trade-offs for asymptotically large dimensions.