Optimal Guidance for Accurate Lunar Soft Landing with Minimum Fuel Consumption using Model Predictive Static Programming

In this paper the soft lunar landing with minimum fuel expenditure is formulated as a nonlinear optimal guidance problem. The realization of pinpoint soft landing with terminal velocity and position constraints is achieved using Model Predictive Static Programming (MPSP). The high accuracy of the terminal conditions is ensured as the formulation of the MPSP inherently poses final conditions as a set of hard constraints. The computational efficiency and fast convergence make the MPSP preferable for fixed final time onboard optimal guidance algorithm. It has also been observed that the minimum fuel requirement strongly depends on the choice of the final time (a critical point that is not given due importance in many literature). Hence, to optimally select the final time, a neural network is used to learn the mapping between various initial conditions in the domain of interest and the corresponding optimal flight time. To generate the training data set, the optimal final time is computed offline using a gradient based optimization technique. The effectiveness of the proposed method is demonstrated with rigorous simulation results.

A multi-level multi-scale approach to study essential genes in Mycobacterium tuberculosis

Background: The set of indispensable genes that are required by an organism to grow and sustain life are termed as essential genes. There is a strong interest in identification of the set of essential genes, particularly in pathogens, not only for a better understanding of the pathogen biology, but also for identifying drug targets and the minimal gene set for the organism. Essentiality is inherently a systems property and requires consideration of the system as a whole for their identification. The available experimental approaches capture some aspects but each method comes with its own limitations. Moreover, they do not explain the basis for essentiality in most cases. A powerful prediction method to recognize this gene pool including rationalization of the known essential genes in a given organism would be very useful. Here we describe a multi-level multi-scale approach to identify the essential gene pool in a deadly pathogen, Mycobacterium tuberculosis. Results: The multi-level workflow analyses the bacterial cell by studying (a) genome-wide gene expression profiles to identify the set of genes which show consistent and significant levels of expression in multiple samples of the same condition, (b) indispensability for growth by using gene expression integrated flux balance analysis of a genome-scale metabolic model, (c) importance for maintaining the integrity and flow in a protein-protein interaction network and (d) evolutionary conservation in a set of genomes of the same ecological niche. In the gene pool identified, the functional basis for essentiality has been addressed by studying residue level conservation and the sub-structure at the ligand binding pockets, from which essential amino acid residues in that pocket have also been identified. 283 genes were identified as essential genes with high-confidence. An agreement of about 73.5% is observed with that obtained from the experimental transposon mutagenesis technique. A large proportion of the identified genes belong to the class of intermediary metabolism and respiration. Conclusions: The multi-scale, multi-level approach described can be generally applied to other pathogens as well. The essential gene pool identified form a basis for designing experiments to probe their finer functional roles and also serve as a ready shortlist for identifying drug targets.

Community-managed forests and wildlife-friendly agriculture play a subsidiary but not substitutive role to protected areas for the endangered Asian elephant

Global conservation policy is increasingly debating the feasibility of reconciling wildlife conservation and human resource requirements in land uses outside protected areas (PAs). However, there are few quantitative assessments of whether or to what extent these `wildlife-friendly' land uses fulfill a fundamental function of PAs-to separate biodiversity from anthropogenic threats. We distinguish the role of wildlife-friendly land uses as being (a) subsidiary, whereby they augment PAs with secondary habitat, or (b) substitutive, wherein they provide comparable habitat to PAs. We tested our hypotheses by investigating the influence of land use and human presence on space-use intensity of the endangered Asian elephant (Elephas maximus) in a fragmented landscape comprising PAs and wildlife-friendly land uses. We applied multistate occupancy models to spatial data on elephant occurrence to estimate and model the overall probability of elephants using a site, and the conditional probability of high-intensity use given that elephants use a site. The probability of elephants using a site regardless of intensity did not vary between PAs and wildlife-friendly land uses. However, high-intensity use declined with distance to PM, and this effect was accentuated by an increase in village density. Therefore, while wildlife-friendly land uses did play a subsidiary conservation role, their potential to substitute for PAs was offset by a strong human presence. Our findings demonstrate the need to evaluate the role of wildlife-friendly land uses in landscape-scale conservation; for species that have conflicting resource requirements with people, PAs are likely to provide crucial refuge from growing anthropogenic threats. (C) 2014 Elsevier Ltd. All rights reserved.

PocketMatch: A new algorithm to compare binding sites in protein structures

Recognizing similarities and deriving relationships among protein molecules is a fundamental requirement in present-day biology. Similarities can be present at various levels which can be detected through comparison of protein sequences or their structural folds. In some cases similarities obscure at these levels could be present merely in the substructures at their binding sites. Inferring functional similarities between protein molecules by comparing their binding sites is still largely exploratory and not as yet a routine protocol. One of the main reasons for this is the limitation in the choice of appropriate analytical tools that can compare binding sites with high sensitivity. To benefit from the enormous amount of structural data that is being rapidly accumulated, it is essential to have high throughput tools that enable large scale binding site comparison. Results: Here we present a new algorithm PocketMatch for comparison of binding sites in a frame invariant manner. Each binding site is represented by 90 lists of sorted distances capturing shape and chemical nature of the site. The sorted arrays are then aligned using an incremental alignment method and scored to obtain PMScores for pairs of sites. A comprehensive sensitivity analysis and an extensive validation of the algorithm have been carried out. A comparison with other site matching algorithms is also presented. Perturbation studies where the geometry of a given site was retained but the residue types were changed randomly, indicated that chance similarities were virtually non-existent. Our analysis also demonstrates that shape information alone is insufficient to discriminate between diverse binding sites, unless combined with chemical nature of amino acids. Conclusion: A new algorithm has been developed to compare binding sites in accurate, efficient and high-throughput manner. Though the representation used is conceptually simplistic, we demonstrate that along with the new alignment strategy used, it is sufficient to enable binding comparison with high sensitivity. Novel methodology has also been presented for validating the algorithm for accuracy and sensitivity with respect to geometry and chemical nature of the site. The method is also fast and takes about 1/250(th) second for one comparison on a single processor. A parallel version on BlueGene has also been implemented.

targetTB: A target identification pipeline for Mycobacterium tuberculosis through an interactome, reactome and genome-scale structural analysis

Background: Tuberculosis still remains one of the largest killer infectious diseases, warranting the identification of newer targets and drugs. Identification and validation of appropriate targets for designing drugs are critical steps in drug discovery, which are at present major bottle-necks. A majority of drugs in current clinical use for many diseases have been designed without the knowledge of the targets, perhaps because standard methodologies to identify such targets in a high-throughput fashion do not really exist. With different kinds of 'omics' data that are now available, computational approaches can be powerful means of obtaining short-lists of possible targets for further experimental validation. Results: We report a comprehensive in silico target identification pipeline, targetTB, for Mycobacterium tuberculosis. The pipeline incorporates a network analysis of the protein-protein interactome, a flux balance analysis of the reactome, experimentally derived phenotype essentiality data, sequence analyses and a structural assessment of targetability, using novel algorithms recently developed by us. Using flux balance analysis and network analysis, proteins critical for survival of M. tuberculosis are first identified, followed by comparative genomics with the host, finally incorporating a novel structural analysis of the binding sites to assess the feasibility of a protein as a target. Further analyses include correlation with expression data and non-similarity to gut flora proteins as well as 'anti-targets' in the host, leading to the identification of 451 high-confidence targets. Through phylogenetic profiling against 228 pathogen genomes, shortlisted targets have been further explored to identify broad-spectrum antibiotic targets, while also identifying those specific to tuberculosis. Targets that address mycobacterial persistence and drug resistance mechanisms are also analysed. Conclusion: The pipeline developed provides rational schema for drug target identification that are likely to have high rates of success, which is expected to save enormous amounts of money, resources and time in the drug discovery process. A thorough comparison with previously suggested targets in the literature demonstrates the usefulness of the integrated approach used in our study, highlighting the importance of systems-level analyses in particular. The method has the potential to be used as a general strategy for target identification and validation and hence significantly impact most drug discovery programmes.

The scale-invariance of spatial patterning in a developing system

Regulating systems, that is, those which exhibit scale-invariant patterns in the adult, are supposed, to do so on account of interactions between cells during development. The nature of these interactions has to be such that the system of positional information (ldquomaprdquo) in the embryo also regulates. To our knowledge, this supposition regarding a regulating map has not been subjected to a direct test in any embryonic system. Here we do so by means of a simple and novel criterion and use it to examine tip regeneration in the mulicellular stage (slug) ofDictyostelium discoideum. When anterior, tip-containing fragments of slugs are amputated, a new tip spontaneously regenerates at the cut surface of the (remaining) posterior fragment. The time needed for regeneration to occur depends on the relative size of the amputated fragment but is independent of the total size of the slug. We conclude from this finding that there is at least one system underlying positional information in the slug which regulates.

A conditional random field-based downscaling method for assessment of climate change impact on multisite daily precipitation in the Mahanadi basin

Downscaling to station-scale hydrologic variables from large-scale atmospheric variables simulated by general circulation models (GCMs) is usually necessary to assess the hydrologic impact of climate change. This work presents CRF-downscaling, a new probabilistic downscaling method that represents the daily precipitation sequence as a conditional random field (CRF). The conditional distribution of the precipitation sequence at a site, given the daily atmospheric (large-scale) variable sequence, is modeled as a linear chain CRF. CRFs do not make assumptions on independence of observations, which gives them flexibility in using high-dimensional feature vectors. Maximum likelihood parameter estimation for the model is performed using limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) optimization. Maximum a posteriori estimation is used to determine the most likely precipitation sequence for a given set of atmospheric input variables using the Viterbi algorithm. Direct classification of dry/wet days as well as precipitation amount is achieved within a single modeling framework. The model is used to project the future cumulative distribution function of precipitation. Uncertainty in precipitation prediction is addressed through a modified Viterbi algorithm that predicts the n most likely sequences. The model is applied for downscaling monsoon (June-September) daily precipitation at eight sites in the Mahanadi basin in Orissa, India, using the MIROC3.2 medium-resolution GCM. The predicted distributions at all sites show an increase in the number of wet days, and also an increase in wet day precipitation amounts. A comparison of current and future predicted probability density functions for daily precipitation shows a change in shape of the density function with decreasing probability of lower precipitation and increasing probability of higher precipitation.

Nonlocal scale effects on wave propagation in multi-walled carbon nanotubes

This paper represents the effect of nonlocal scale parameter on the wave propagation in multi-walled carbon nanotubes (MWCNTs). Each wall of the MWCNT is modeled as first order shear deformation beams and the van der Waals interactions between the walls are modeled as distributed springs. The studies shows that the scale parameter introduces certain band gap region in both flexural and shear wave mode where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite (or group speed tends to zero). The frequency at which this phenomenon occurs is called the ``Escape frequency''. The analysis shows that, for a given N-walled carbon nanotube (CNT). the nonlocal scaling parameter has a significant effect on the shear wave modes of the N - 1 walls. The escape frequencies of the flexural and shear wave modes of the N-walls are inversely proportionl to the nonlocal scaling parameter. It is also shown that the cut-off frequencies are independent of the nonlocal scale parameter. (C) 2009 Elsevier B.V. All rights reserved.

Diurnal-scale signatures of monsoon rainfall over the Indian region from TRMM satellite observations

One of the most important modes of summer season precipitation variability over the Indian region, the diurnal cycle, is studied using the Tropical Rainfall Measuring Mission 3-hourly, 0.25 degrees x 0.25 degrees 3B42 rainfall product for nine years (1999-2007). Most previous studies have provided an analysis of a single year or a few years of satellite-or station-based rainfall data. Our study aims to systematically analyze the statistical characteristics of the diurnal-scale signature of rainfall over the Indian and surrounding regions. Using harmonic analysis, we extract the signal corresponding to diurnal and subdiurnal variability. Subsequently, the 3-hourly time period or the octet of rainfall peak for this filtered signal, referred to as the ``peak octet,'' is estimated, with care taken to eliminate spurious peaks arising out of Gibbs oscillations. Our analysis suggests that over the Bay of Bengal, there are three distinct modes of the peak octet of diurnal rainfall corresponding to 1130, 1430, and 1730 Indian standard time (IST), from the north central to south bay. This finding could be seen to be consistent with southward propagation of the diurnal rainfall pattern reported by earlier studies. Over the Arabian Sea, there is a spatially coherent pattern in the mode of the peak octet (1430 IST), in a region where it rains for more than 30% of the time. In the equatorial Indian Ocean, while most of the western part shows a late night/early morning peak, the eastern part does not show a spatially coherent pattern in the mode of the peak octet owing to the occurrence of a ual maxima (early morng and early/late afternoon). The imalayan foothills were found to have a mode of peak octet corresponding to 0230 IST, whereas over the Burmese mountains and the Western Ghats (west coast of India) the rainfall peaks during late afternoon/early evening (1430-1730 IST). This implies that the phase of the diurnal cycle over inland orography (e. g., Himalayas) is significantly different from coastal orography (e. g., Western Ghats). We also find that over the Gangetic plains, the peak octet is around 1430 IST, a few hours earlier compared to the typical early evening maxima over land.

Benchmark scenarios for the NMSSM

We discuss constrained and semi--constrained versions of the next--to--minimal supersymmetric extension of the Standard Model (NMSSM) in which a singlet Higgs superfield is added to the two doublet superfields that are present in the minimal extension (MSSM). This leads to a richer Higgs and neutralino spectrum and allows for many interesting phenomena that are not present in the MSSM. In particular, light Higgs particles are still allowed by current constraints and could appear as decay products of the heavier Higgs states, rendering their search rather difficult at the LHC. We propose benchmark scenarios which address the new phenomenological features, consistent with present constraints from colliders and with the dark matter relic density, and with (semi--)universal soft terms at the GUT scale. We present the corresponding spectra for the Higgs particles, their couplings to gauge bosons and fermions and their most important decay branching ratios. A brief survey of the search strategies for these states at the LHC is given.

Why do birds participate in mixed-species foraging flocks? A large-scale synthesis

Mixed-species flocks of foraging birds have been documented from terrestrial habitats all over the world and are thought to form for either improved feeding efficiency or better protection from predators. Two kinds of flock participants are recognized: those that join other species ('followers') and are therefore likely to be the recipients of the benefits of flock participation and those that are joined ('leaders'). Through comparative analyses, using a large sample of flocks from around the world, we show that (1) 'followers' tend to be smaller, more insectivorous, and feed in higher strata than matched species that participate in flocks to a lesser extent and (2) 'leaders' tend to be cooperative breeders more often than matched species that are not known to lead flocks. Furthermore, meta-analyses of published results from across the world showed that bird species in terrestrial mixed-species flocks increase foraging rates and reduce vigilance compared to when they are solitary or in conspecific groups. Moreover, the increase in foraging rates is seen only with flock followers and not flock leaders. These findings suggest a role for predation in the evolution of mixed-species flocking. Species that are vulnerable to predation follow species whose vigilance they can exploit. By doing so, they are able to reduce their own vigilance and forage at higher rates. (C) 2009 The Association for the Study of Animal Behaviour. Published by Elsevier Ltd. All rights reserved.

Scale-up criteria of square tank surface aerator

Oxygen transfer rate and the corresponding power requirement to operate the rotor are vital for design and scale-up of surface aerators. Present study develops simulation or scale-up criterion correlating the oxygen transsimulation fer coefficient and power number along with a parameter governing theoretical power per unit volume (X, which is defined as equal to (FR1/3)-R-4/3, where F and R are impellers' Fronde and Reynolds number, respectively). Based on such scale-up criteria, design considerations are developed to save energy requirements while designing square tank surface aerators. It has been demonstrated that energy can be saved substantially if the aeration tanks are run at relatively higher input powers. It is also demonstrated that smaller sized tanks are more energy conservative and economical when compared to big sized tanks, while aerating the same volume of water, and at the same time by maintaining a constant input power in all the tanks irrespective of their size. An example illustrating how energy can be reduced while designing different sized aerators is given. The results presented have a wide application in biotechnology and bioengineering areas with a particular emphasis on the design of appropriate surface aeration systems.

Scale-up criteria of square tank surface aerator

Oxygen transfer rate and the corresponding power requirement to operate the rotor are vital for design and scale-up of surface aerators. Present study develops simulation or scale-up criterion correlating the oxygen transsimulation fer coefficient and power number along with a parameter governing theoretical power per unit volume (X, which is defined as equal to (FR1/3)-R-4/3, where F and R are impellers' Fronde and Reynolds number, respectively). Based on such scale-up criteria, design considerations are developed to save energy requirements while designing square tank surface aerators. It has been demonstrated that energy can be saved substantially if the aeration tanks are run at relatively higher input powers. It is also demonstrated that smaller sized tanks are more energy conservative and economical when compared to big sized tanks, while aerating the same volume of water, and at the same time by maintaining a constant input power in all the tanks irrespective of their size. An example illustrating how energy can be reduced while designing different sized aerators is given. The results presented have a wide application in biotechnology and bioengineering areas with a particular emphasis on the design of appropriate surface aeration systems.

Growth, self-assembly and dynamics of nano-scale films at fluid interfaces

Ultrathin films at fluid interfaces are important not only from a fundamental point of view as 2D complex fluids but have also become increasingly relevant in the development of novel functional materials. There has been an explosion in the synthesis work in this area over the last decade, giving rise to many exotic nanostructures at fluid interfaces. However, the factors controlling particle nucleation, growth and self-assembly at interfaces are poorly understood on a quantitative level. We will outline some of the recent attempts in this direction. Some of the selected investigations examining the macroscopic mechanical properties of molecular and particulate films at fluid interfaces will be reviewed. We conclude with a discussion of the electronic properties of these films that have potential technological and biological applications.

Photophysics of conjugated polymers: interplay between Forster energy migration and defect concentration in shaping a photochemical funnel in PPV

Recent single molecule experiments have suggested the existence of a photochemical funnel in the photophysics of conjugated polymers, like poly[2-methoxy-5-(2'-ethylhexyl)oxy-1,4-phenylenevinylene] (MEH-PPV). The funnel is believed to be a consequence of the presence of conformational or chemical defects along the polymer chain and efficient non-radiative energy transfer among different chromophore segments. Here we address the effect of the excitation energy dynamics on the photophysics of PPV. The PPV chain is modeled as a polymer with the length distribution of chromophores given either by a Gaussian or by a Poisson distribution. We observe that the Poisson distribution of the segment lengths explains the photophysics of PPV better than the Gaussian distribution. A recently proposed version of an extended particle-in-a-box' model is used to calculate the exciton energies and the transition dipole moments of the chromophores, and a master equation to describe the excitation energy transfer among different chromophores. The rate of energy transfer is assumed to be given here, as a first approximation, by the well-known Forster expression. The observed excitation population dynamics confirms the photochemical funneling of excitation energy from shorter to longer chromophores of the polymer chain. The time scale of spectral shift and energy transfer for our model polymer, with realistic values of optical parameters, is in the range of 200-300 ps. We find that the excitation energy may not always migrate towards the longest chromophore segments in the polymer chain as the efficiency of energy transfer between chromophores depends on the separation distance between the two and their relative orientation.