Performance metrics in a hybrid MPI-OpenMP based molecular dynamics simulation with short-range interactions

We discuss the computational bottlenecks in molecular dynamics (MD) and describe the challenges in parallelizing the computation-intensive tasks. We present a hybrid algorithm using MPI (Message Passing Interface) with OpenMP threads for parallelizing a generalized MD computation scheme for systems with short range interatomic interactions. The algorithm is discussed in the context of nano-indentation of Chromium films with carbon indenters using the Embedded Atom Method potential for Cr-Cr interaction and the Morse potential for Cr-C interactions. We study the performance of our algorithm for a range of MPI-thread combinations and find the performance to depend strongly on the computational task and load sharing in the multi-core processor. The algorithm scaled poorly with MPI and our hybrid schemes were observed to outperform the pure message passing scheme, despite utilizing the same number of processors or cores in the cluster. Speed-up achieved by our algorithm compared favorably with that achieved by standard MD packages. (C) 2013 Elsevier Inc. All rights reserved.

Phylogeography unplugged: comparative surveys in the genomic era

In March 2012, the authors met at the National Evolutionary Synthesis Center (NESCent) in Durham, North Carolina, USA, to discuss approaches and cooperative ventures in Indo-Pacific phylogeography. The group emerged with a series of findings: (1) Marine population structure is complex, but single locus mtDNA studies continue to provide powerful first assessment of phylogeographic patterns. (2) These patterns gain greater significance/power when resolved in a diversity of taxa. New analytical tools are emerging to address these analyses with multi-taxon approaches. (3) Genome-wide analyses are warranted if selection is indicated by surveys of standard markers. Such indicators can include discordance between genetic loci, or between genetic loci and morphology. Phylogeographic information provides a valuable context for studies of selection and adaptation. (4) Phylogeographic inferences are greatly enhanced by an understanding of the biology and ecology of study organisms. (5) Thorough, range-wide sampling of taxa is the foundation for robust phylogeographic inference. (6) Congruent geographic and taxonomic sampling by the Indo-Pacific community of scientists would facilitate better comparative analyses. The group concluded that at this stage of technology and software development, judicious rather than wholesale application of genomics appears to be the most robust course for marine phylogeographic studies. Therefore, our group intends to affirm the value of traditional (''unplugged'') approaches, such as those based on mtDNA sequencing and microsatellites, along with essential field studies, in an era with increasing emphasis on genomic approaches.

A Harmonic Suppression Scheme for Open-End Winding Split-Phase IM Drive Using Capacitive Filters for the Full Speed Range

Voltage source inverter (VSI)-fed six-phase induction motor (IM) drives have high 6n +/- 1, n = odd-order harmonic currents. This is because these currents, driven by the corresponding harmonic voltages in the inverter output, are limited only by the stator leakage impedance, as these harmonics are absent in the back electromotive force of the motor. To suppress the harmonic currents, either bulky inductive harmonic filters or complex pulsewidth modulation (PWM) techniques have to be used. This paper proposes a harmonic elimination scheme using switched capacitor filters for a VSI-fed split-phase IM drive. Two 3-phase inverters fed from capacitors are used on the open-end side of the motor to suppress 6n +/- 1, n = odd-order harmonics. A PWM scheme that can suppress the harmonics as well as balance the capacitor voltage is also proposed. The capacitor fed inverters are switched so that the fundamental voltage is not affected, and the fundamental power is always drawn from the main inverters. The proposed scheme is verified with a detailed experimental study. The effectiveness of the scheme is demonstrated by comparing the results with those obtained by disabling the capacitor fed inverters.

Micro-mechanical stages with enhanced range

We present concepts and an optimization-based methodology for the design of micro-mechanical stages that have not only high precision but also an enhanced range. Joint-free distributed compliant designs provide high precision and easy manufacturability at macro and micro scales. The range of motion is enhanced by using displacement-amplifying compliant mechanisms (DaCMs). The main issue addressed in this paper is how to retain the decoupling between the X and Y motions in the stage when it is equipped with DaCMs. The natural frequency of the stage is also not compromised in enhancing the range. The optimized design has 2.5 times more range than the designs reported in the literature. Furthermore, the sensitivity improved by a factor of two when the stage is optimized for an accelerometer.

Variability of cloud liquid water and ice over South Asia from TMI estimates

In this study, the Tropical Rainfall Measurement Mission based Microwave Imager estimates (2A12) have been used to compare and contrast the characteristics of cloud liquid water and ice over the Indian land region and the ocean surrounding it, during the premonsoon (May) and monsoon (June-September) seasons. Based on the spatial homogeneity of rainfall, we have selected five regions for our study (three over ocean, two over land). Comparison across three ocean regions suggests that the cloud liquid water (CLW) over the orographically influenced Arabian Sea (close to the Indian west coast) behaves differently from the CLW over a trapped ocean (Bay of Bengal) or an open ocean (equatorial Indian Ocean). Specifically, the Arabian Sea region shows higher liquid water for a lower range of rainfall, whereas the Bay of Bengal and the equatorial Indian Ocean show higher liquid water for a higher range of rainfall. Apart from geographic differences, we also documented seasonal differences by comparing CLW profiles between monsoon and premonsoon periods, as well as between early and peak phases of the monsoon. We find that the CLW during the lean periods of rainfall (May or June) is higher than during the peak and late monsoon season (July-September) for raining clouds. As active and break phases are important signatures of the monsoon progression, we also analysed the differences in CLW during various phases of the monsoon, namely, active, break, active-to-break and break-to-active transition phases. We find that the cloud liquid water content during the break-to-active transition phase is significantly higher than during the active-to-break transition phase over central India. We speculate that this could be attributed to higher amount of aerosol loading over this region during the break phase. We lend credence to this aerosol-CLW/rain association by comparing the central Indian CLW with that over southeast Asia (where the aerosol loading is significantly smaller) and find that in the latter region, there are no significant differences in CLW during the different phases of the monsoon. While our hypothesis needs to be further investigated with numerical models, the results presented in this study can potentially serve as a good benchmark in evaluating the performance of cloud resolving models over the Indian region.

12 V / kilo-Farad Range Lead-Carbon Hybrid Ultracapacitors and Their Envisaged Applications

12 V / kilo-Farad (kF) range substrate-integrated lead-carbon hybrid ultracapacitors (HUCs) wherein the conventional positive plates of lead-acid batteries are replaced with substrate-integrated PbO2 positive plates and the negative plates are replaced with carbon-coated graphitic electrodes, providing totally non-faradaic and corrosion-free electrodes, are developed and performance tested. Constant-current discharge data at varying load-currents, constant-power discharge data at varying power values, and the capacitance data at different temperature for a 12 V / kF range substrate-integrated lead-carbon HUC are described along with its resistance, leakage current, self-discharge and cycle-life characteristics.

A Stochastic Chemical Dynamic Approach to Correlate Autoimmunity and Optimal Vitamin-D Range

Motivated by several recent experimental observations that vitamin-D could interact with antigen presenting cells (APCs) and T-lymphocyte cells (T-cells) to promote and to regulate different stages of immune response, we developed a coarse grained but general kinetic model in an attempt to capture the role of vitamin-D in immunomodulatory responses. Our kinetic model, developed using the ideas of chemical network theory, leads to a system of nine coupled equations that we solve both by direct and by stochastic (Gillespie) methods. Both the analyses consistently provide detail information on the dependence of immune response to the variation of critical rate parameters. We find that although vitamin-D plays a negligible role in the initial immune response, it exerts a profound influence in the long term, especially in helping the system to achieve a new, stable steady state. The study explores the role of vitamin-D in preserving an observed bistability in the phase diagram (spanned by system parameters) of immune regulation, thus allowing the response to tolerate a wide range of pathogenic stimulation which could help in resisting autoimmune diseases. We also study how vitamin-D affects the time dependent population of dendritic cells that connect between innate and adaptive immune responses. Variations in dose dependent response of anti-inflammatory and pro-inflammatory T-cell populations to vitamin-D correlate well with recent experimental results. Our kinetic model allows for an estimation of the range of optimum level of vitamin-D required for smooth functioning of the immune system and for control of both hyper-regulation and inflammation. Most importantly, the present study reveals that an overdose or toxic level of vitamin-D or any steroid analogue could give rise to too large a tolerant response, leading to an inefficacy in adaptive immune function.

Cryptic species and Miocene diversification of Palaearctic naked-toed geckos (Squamata: Gekkonidae) in the Indian dry zone

We sampled Palaearctic naked-toed geckos from across their range in India and used two mitochondrial and two nuclear genes to reconstruct relationships within a global phylogeny. Published sequences of Peninsular Indian Hemidactylus allow us to contrast these two groups in dating analyses - providing insights into the history of the Indian dry zone. Palaearctic naked-toed geckos first moved onto the Indian Plate in the Oligocene, with higher-level diversification probably linked to collision of the Indian and Eurasian plates, and subsequent dispersal into-India and diversification with increasing Miocene aridity. An independent gekkonid radiation with species in the dry zone, Hemidactylus diversified during the same period in Peninsular India. Our results demonstrate that dry zone taxa across India may date back to at least the Miocene, with a potential historical climatic barrier between the Indus and Peninsular Indian Divisions. `Cyrtopodion' aravallense is revealed to be a complex with seven genetically and environmentally divergent lineages that began diversifying in the late Miocene, congruent with increased aridity in north-western India. This discovery of cryptic diversity in the Indian dry zone represents the first terrestrial vertebrate radiation from north-western central India and highlights how little we understand of the regions' biodiversity, emphasizing the need for systematic geographic sampling and multiline evidence to reveal true patterns of diversity. The ancestor of `Cyrtopodion' aravallense came into the nascent Indian dry zone in the Miocene and has since diversified, potentially in the absence of any sympatric scansorial rupicolous geckos. Cyrtopodion scabrum represents a unique case of commensalism and shows phylogeographic structure in its presumed native range. The taxonomic implications of our study include a number of undescribed species, recognition of `Cyrtopodion' as a distinct lineage and the non-monophyly of Altiphylax.

Long range emission enhancement and anisotropy in coupled quantum dots induced by aligned gold nanoantenna

Quantum dot arrays have been projected as the material of choice for next generation displays and photodetectors. Extensive ongoing research aims at improving optical and electrical efficiencies of such devices. We report experimental results on non-local long range emission intensity enhancement and anisotropy in quantum dot assemblies induced by isolated and partially aligned gold nanoantennas. Spatially resolved photoluminescence clearly demonstrate that the effect is maximum, when the longitudinal surface plasmon resonance of the nanoantenna is resonant with the emission maxima of the quantum dots. We estimated the decay length of this enhancement to be similar to 2.6 mu m, which is considerably larger than the range of near field interaction of metal nanoantenna. Numerical simulations qualitatively capture the near field behavior of the nanorods but fail to match the experimentally observed non-local effects. We have suggested how strong interactions of quantum dots in the close packed assemblies, mediated by the nanoantennas, could lead to such observed behavior. (C) 2014 AIP Publishing LLC.

Orientational order as the origin of the long-range hydrophobic effect

The long range attractive force between two hydrophobic surfaces immersed in water is observed to decrease exponentially with their separation-this distance-dependence of effective force is known as the hydrophobic force law (HFL). We explore the microscopic origin of HFL by studying distance-dependent attraction between two parallel rods immersed in 2D Mercedes Benz model of water. This model is found to exhibit a well-defined HFL. Although the phenomenon is conventionally explained by density-dependent theories, we identify orientation, rather than density, as the relevant order parameter. The range of density variation is noticeably shorter than that of orientational heterogeneity. The latter is comparable to the observed distances of hydrophobic force. At large separation, attraction between the rods arises primarily from a destructive interference among the inwardly propagating oppositely oriented heterogeneity generated in water by the two rods. As the rods are brought closer, the interference increases leading to a decrease in heterogeneity and concomitant decrease in free energy of the system, giving rise to the effective attraction. We notice formation of hexagonal ice-like structures at the onset of attractive region which suggests that metastable free energy minimum may play a role in the origin of HFL. (C) 2015 AIP Publishing LLC.

HCN channels enhance spike phase coherence and regulate the phase of spikes and LFPs in the theta-frequency range

What are the implications for the existence of subthreshold ion channels, their localization profiles, and plasticity on local field potentials (LFPs)? Here, we assessed the role of hyperpolarization-activated cyclic-nucleotide-gated (HCN) channels in altering hippocampal theta-frequency LFPs and the associated spike phase. We presented spatiotemporally randomized, balanced theta-modulated excitatory and inhibitory inputs to somatically aligned, morphologically realistic pyramidal neuron models spread across a cylindrical neuropil. We computed LFPs from seven electrode sites and found that the insertion of an experimentally constrained HCN-conductance gradient into these neurons introduced a location- dependent lead in the LFP phase without significantly altering its amplitude. Further, neurons fired action potentials at a specific theta phase of the LFP, and the insertion of HCN channels introduced large lags in this spike phase and a striking enhancement in neuronal spike-phase coherence. Importantly, graded changes in either HCN conductance or its half-maximal activation voltage resulted in graded changes in LFP and spike phases. Our conclusions on the impact of HCN channels on LFPs and spike phase were invariant to changes in neuropil size, to morphological heterogeneity, to excitatory or inhibitory synaptic scaling, and to shifts in the onset phase of inhibitory inputs. Finally, we selectively abolished the inductive lead in the impedance phase introduced by HCN channels without altering neuronal excitability and found that this inductive phase lead contributed significantly to changes in LFP and spike phase. Our results uncover specific roles for HCN channels and their plasticity in phase-coding schemas and in the formation and dynamic reconfiguration of neuronal cell assemblies.

Electric field induced short range to long range structural ordering and its influence on the Eu+3 photoluminescence in the lead-free ferroelectric Na1/2Bi1/2TiO3

Eu+3 was incorporated into the lattice of a lead-free ferroelectric Na1/2Bi1/2TiO3 (NBT) as per the nominal formula Na0.5Bi0.5-xEuxTiO3. This system was investigated with regard to the Eu+3 photoluminescence (PL) and structural behaviour as a function of composition and electric field. Electric field was found to irreversibly change the features in the PL spectra and also in the x-ray diffraction patterns below the critical composition x = 0.025. Detailed analysis revealed that below the critical composition, electric field irreversibly suppresses the structural heterogeneity inherent of the host matrix NBT and brings about a long range ferroelectric state with rhombohedral (R3c) distortion. It is shown that the structural disorder on the nano-scale opens a new channel for radiative transition which manifests as a new emission line branching off from the main D-5(0)-> F-7(0) line along with a concomitant change in the relative intensity of the other crystal field induced Stark lines with different J values. The study suggests that Eu+3 luminescence can be used to probe the relative degree of field induced structural ordering in relaxor ferroelectrics and also in high performance piezoelectric alloys where electric field couples very strongly with the lattice and structural degrees of freedom. (C) 2015 AIP Publishing LLC.

Active dendrites mediate stratified gamma-range coincidence detection in hippocampal model neurons

Hippocampal pyramidal neurons exhibit gamma-phase preference in their spikes, selectively route inputs through gamma frequency multiplexing and are considered part of gamma-bound cell assemblies. How do these neurons exhibit gamma-frequency coincidence detection capabilities, a feature that is essential for the expression of these physiological observations, despite their slow membrane time constant? In this conductance-based modelling study, we developed quantitative metrics for the temporal window of integration/coincidence detection based on the spike-triggered average (STA) of the neuronal compartment. We employed these metrics in conjunction with quantitative measures for spike initiation dynamics to assess the emergence and dependence of coincidence detection and STA spectral selectivity on various ion channel combinations. We found that the presence of resonating conductances (hyperpolarization-activated cyclic nucleotide-gated or T-type calcium), either independently or synergistically when expressed together, led to the emergence of spectral selectivity in the spike initiation dynamics and a significant reduction in the coincidence detection window (CDW). The presence of A-type potassium channels, along with resonating conductances, reduced the STA characteristic frequency and broadened the CDW, but persistent sodium channels sharpened the CDW by strengthening the spectral selectivity in the STA. Finally, in a morphologically precise model endowed with experimentally constrained channel gradients, we found that somatodendritic compartments expressed functional maps of strong theta-frequency selectivity in spike initiation dynamics and gamma-range CDW. Our results reveal the heavy expression of resonating and spike-generating conductances as the mechanism underlying the robust emergence of stratified gamma-range coincidence detection in the dendrites of hippocampal and cortical pyramidal neurons.

Solid-state optical absorption from optimally tuned time-dependent range-separated hybrid density functional theory

We present a framework for obtaining reliable solid-state charge and optical excitations and spectra from optimally tuned range-separated hybrid density functional theory. The approach, which is fully couched within the formal framework of generalized Kohn-Sham theory, allows for the accurate prediction of exciton binding energies. We demonstrate our approach through first principles calculations of one- and two-particle excitations in pentacene, a molecular semiconducting crystal, where our work is in excellent agreement with experiments and prior computations. We further show that with one adjustable parameter, set to produce the known band gap, this method accurately predicts band structures and optical spectra of silicon and lithium fluoride, prototypical covalent and ionic solids. Our findings indicate that for a broad range of extended bulk systems, this method may provide a computationally inexpensive alternative to many-body perturbation theory, opening the door to studies of materials of increasing size and complexity.