The fully implicit stochastic-alpha method for stiff stochastic differential equations

A fully implicit integration method for stochastic differential equations with significant multiplicative noise and stiffness in both the drift and diffusion coefficients has been constructed, analyzed and illustrated with numerical examples in this work. The method has strong order 1.0 consistency and has user-selectable parameters that allow the user to expand the stability region of the method to cover almost the entire drift-diffusion stability plane. The large stability region enables the method to take computationally efficient time steps. A system of chemical Langevin equations simulated with the method illustrates its computational efficiency.

Spatial correlation in grain misorientation distribution

Grain misorientation was studied in relation to the nearest neighbor's mutual distance using electron back-scattered diffraction measurements. The misorientation correlation function was defined as the probability density for the occurrence of a certain misorientation between pairs of grains separated by a certain distance. Scale-invariant spatial correlation between neighbor grains was manifested by a power law dependence of the preferred misorientation vs. inter-granular distance in various materials after diverse strain paths. The obtained negative scaling exponents were in the range of -2 +/- 0.3 for high-angle grain boundaries. The exponent decreased in the presence of low-angle grain boundaries or dynamic recrystallization, indicating faster decay of correlations. The correlations vanished in annealed materials. The results were interpreted in terms of lattice incompatibility and continuity conditions at the interface between neighboring grains. Grain-size effects on texture development, as well as the implications of such spatial correlations on texture modeling, were discussed.

Influence of surface chemistry of mesoporous alumina with wide pore distribution on controlled drug release

The crucial role of the drug carrier surface chemical moeities on the uptake and in vitro release of drug is discussed here in a systematic manner. Mesoporous alumina with a wide pore size distribution (2-7 nm) functionalized with various hydrophilic and hydrophobic surface chemical groups was employed as the carrier for delivery of the model drug ibuprofen. Surface functionalization with hydrophobic groups resulted in low degree of drug loading (approximately 20%) and fast rate of release (85% over a period of 5 h) whereas hydrophilic groups resulted in a significantly higher drug payloads (21%-45%) and slower rate of release (12%-40% over a period of 5 h). Depending on the chemical moiety, the diffusion controlled (proportional to time(-0.5)) drug release was additionally observed to be dependent on the mode of arrangement of the functional groups on the alumina surface as well as on the pore characteristics of the matrix. For all mesoporous alumina systems the drug dosages were far lower than the maximum recommended therapeutic dosages (MRTD) for oral delivery. We envisage that the present study would aid in the design of delivery systems capable of sustained release of multiple drugs.

Vertical distribution of aerosols over the east coast of India inferred from airborne LIDAR measurements

The information on altitude distribution of aerosols in the atmosphere is essential in assessing the impact of aerosol warming on thermal structure and stability of the atmosphere.In addition, aerosol altitude distribution is needed to address complex problems such as the radiative interaction of aerosols in the presence of clouds. With this objective,an extensive, multi-institutional and multi-platform field experiment (ICARB-Integrated Campaign for Aerosols, gases and Radiation Budget) was carried out under the Geosphere Biosphere Programme of the Indian Space Research Organization (ISRO-GBP) over continental India and adjoining oceans during March to May 2006. Here, we present airborne LIDAR measurements carried out over the east Coast of the India during the ICARB field campaign. An increase in aerosol extinction (scattering + absorption) was observed from the surface upwards with a maximum around 2 to 4 km. Aerosol extinction at higher atmospheric layers (>2 km) was two to three times larger compared to that of the surface. A large fraction (75-85%) of aerosol column optical depth was contributed by aerosols located above 1 km. The aerosol layer heights (defined in this paper as the height at which the gradient in extinction coefficient changes sign) showed a gradual decrease with an increase in the offshore distance. A large fraction (60-75%) of aerosol was found located above clouds indicating enhanced aerosol absorption above clouds. Our study implies that a detailed statistical evaluation of the temporal frequency and spatial extent of elevated aerosol layers is necessary to assess their significance to the climate. This is feasible using data from space-borne lidars such as CALIPSO,which fly in formation with other satellites like MODIS AQUA and MISR, as part of the A-Train constellation.

Continuous Distribution Kinetics for Photopolymerization of Alkyl Methacrylates

The photopolymerization of methyl,ethyl,butyl, and hexyl methacrylates in solution was studied. The effect of initial initiator and monomer concentrations on the time evolution of polymer concentration (M) over bar (n) and PDI was examined. The reversible chain addition and beta-scission, and primary radical termination steps were included in the mechanism along with the classical steps. The rate equations were derived using continuous distribution kinetics and solved numerically to fit the experimental data. The regressed rate coefficients compared well with the literature data. The model predicted the instantaneous increase in (M) over bar (n) and PDI to steady state values. The rate coefficients exhibited a linear increase with the size of alkyl chain of the alkyl methacrylates.

The alpha method for solving differential algebraic inequality (DAI) systems

This paper describes an algorithm for ``direct numerical integration'' of the initial value Differential-Algebraic Inequalities (DAI) in a time stepping fashion using a sequential quadratic programming (SQP) method solver for detecting and satisfying active path constraints at each time step. The activation of a path constraint generally increases the condition number of the active discretized differential algebraic equation's (DAE) Jacobian and this difficulty is addressed by a regularization property of the alpha method. The algorithm is locally stable when index 1 and index 2 active path constraints and bounds are active. Subject to available regularization it is seen to be stable for active index 3 active path constraints in the numerical examples. For the high index active path constraints, the algorithm uses a user-selectable parameter to perturb the smaller singular values of the Jacobian with a view to reducing the condition number so that the simulation can proceed. The algorithm can be used as a relatively cheaper estimation tool for trajectory and control planning and in the context of model predictive control solutions. It can also be used to generate initial guess values of optimization variables used as input to inequality path constrained dynamic optimization problems. The method is illustrated with examples from space vehicle trajectory and robot path planning.

A lower bound to the survival probability and an approximate first passage time distribution for Markovian and non-Markovian dynamics in phase space

We derive a very general expression of the survival probability and the first passage time distribution for a particle executing Brownian motion in full phase space with an absorbing boundary condition at a point in the position space, which is valid irrespective of the statistical nature of the dynamics. The expression, together with the Jensen's inequality, naturally leads to a lower bound to the actual survival probability and an approximate first passage time distribution. These are expressed in terms of the position-position, velocity-velocity, and position-velocity variances. Knowledge of these variances enables one to compute a lower bound to the survival probability and consequently the first passage distribution function. As examples, we compute these for a Gaussian Markovian process and, in the case of non-Markovian process, with an exponentially decaying friction kernel and also with a power law friction kernel. Our analysis shows that the survival probability decays exponentially at the long time irrespective of the nature of the dynamics with an exponent equal to the transition state rate constant.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.

A study on the variations of optical and physical properties of aerosols over a tropical semi-arid station during grassland fire

The present paper records the results of a case study on the impact of an extensive grassland fire on the physical and optical properties of aerosols at a semi-arid station in southern India for the first time from ground based measurements using a MICROTOPS-II sunphotometer, an aethalometer and a quartz crystal microbalance impactor (QCM). Observations revealed a substantial increase in aerosol optical depth (AOD) at all wavelengths during burning days compared to normal days. High AOD values observed at shorter wavelengths suggest the dominance of accumulation mode particle loading over the study area. Daily mean aerosol size spectra shows, most of the time, power-law distribution. To characterize AOD, the Angstrom parameters (i.e., alpha and beta) were used. Wavelength exponent (1.38) and turbidity coefficient (0.21) are high during burning days compared to normal days, thereby suggesting an increase in accumulation mode particle loading. Aerosol size distribution suggested dominance of accumulation mode particle loading during burning days compared to normal days. A significant positive correlation was observed between AOD at 500 mn and water vapour and negative correlation between AOD at 500 nm and wind speed for burning and non-burning days. Diurnal variations of black carbon (BC) aerosol mass concentrations increased by a factor of similar to 2 in the morning and afternoon hours during burning period compared to normal days.

High Resolution Surface Imaging Using a Carbon Nanotube Array with Pointed Height distribution

In this paper we discuss a new technique to image the surfaces of metallic substrates using field emission from a pointed array of carbon nanotubes (CNTs). We consider a pointed height distribution of the CNT array under a diode configuration with two side gates maintained at a negative potential to obtain a highly intense beam of electrons localized at the center of the array. The CNT array on a metallic substrate is considered as the cathode and the test substrate as the anode. Scanning the test Substrate with the cathode reveals that the field emission current is highly sensitive to the surface features with nanometer resolution. Surface features of semi-circular, triangular and rectangular geometries (projections and grooves) are considered for simulation. This surface scanning/mapping technique can be applied for surface roughness measurements with nanoscale accuracy. micro/nano damage detection, high precision displacement sensors, vibrometers and accelerometers. among other applications.

Anomalous reaction-diffusion as a model of nonexponential DNA escape kinetics

We show that data from recent experiments carried out on the kinetics of DNA escape from alpha-hemolysin nanopores [M. Wiggin, C. Tropini, C. T. Cossa, N. N. Jetha, and A. Marziali, Biophys. J. 95, 5317 (2008)] may be rationalized by a model of chain dynamics based on the anomalous diffusion of a particle moving in a harmonic well in the presence of a delta function sink. The experiments of Wiggin found, among other things, that the occasional occurrence of unusually long escape times in the distribution of chain trapping events led to nonexponential decays in the survival probability, S(t), of the DNA molecules within the nanopore. Wiggin ascribed this nonexponentiality to the existence of a distribution of trapping potentials, which they suggested was theresult of stochastic interactions between the bases of the DNA and the amino acids located on the surface of the nanopore. Based on this idea, they showed that the experimentally determined S(t) could be well fit in both the short and long time regimes by a function of the form (1+t/tau)(-alpha) (the so called Becquerel function). In our model, S(t) is found to be given by a Mittag-Leffler function at short times and by a generalized Mittag-Leffler function at long times. By suitable choice of certain parameter values, these functions are found to fit the experimental S(t) even better than the Becquerel function. Anomalous diffusion of DNA within the trap prior to escape over a barrier of fixed height may therefore provide a second, plausible explanation of the data, and may offer fresh perspectives on similar trapping and escape problems.

H N.M.R. studies of protected alpha-aminoisobutyric acid containing peptides. Chemical shift nonequivalence of benzyloxycarbonyl methylene protons

The benzylic methylene protons in a large number of benzyloxycarbonyl alpha-aminoisobutyric acid (Z-Aib) containing peptides, show chemical shift nonequivalence. The magnitude of the geminal nonequivalence is correlated with the involvement of the urethane carbonyl group, in an intramolecular hydrogen bond. Studies of the model compounds Z-Aib-Aib-Ala-NHMe, and Z-Aib-Aib-Aib-Pro-OMe clearly establish the presence of intramolecular hydrogen bonds, involving the urethane CO group. In both compounds marked anisochrony of the benzylic methylene protons is demonstrated. In Z-Aib-Aib-Pro-OMe, where a 4 leads to 1 hydrogen bonded beta-turn is not possible, the benzylic-CH2-protons appear as a singlet in CDCl3 and have a very small chemical shift difference in (CD3)2SO. The observation of such nonequivalence is of value in establishing whether the amino terminal Aib-Pro beta-turn is retained in large peptide-fragments of alamethicin.

CD4+ alpha beta T cell and gamma delta T cell responses to Mycobacterium tuberculosis. Similarities and differences in Ag recognition, cytotoxic effector function, and cytokine production.

CD4+ and gamma delta T cells are activated readily by Mycobacterium tuberculosis. To examine their role in the human immune response to M. tuberculosis, CD4+ and gamma delta T cells from healthy tuberculin-positive donor were studied for patterns of Ag recognition, cytotoxicity, and cytokine production in response to M. tuberculosis-infected mononuclear phagocytes. Both T cell subsets responded to intact M. tuberculosis and its cytosolic Ags. However, CD4+ and gamma delta T cells differed in the range of cytosolic Ags recognized: reactivity to a wide m.w. range of Ags for CD4+ T cells, and a restricted pattern for gamma delta T cells, with dominance of Ags of 10 to 15 kDa. Both T cell subsets were equally cytotoxic for M. tuberculosis-infected monocytes. Furthermore, both CD4+ and gamma delta T cells produced large amounts of IFN-gamma: mean pg/ml of IFN-gamma in supernatants was 2458 +/- 213 for CD4+ and 2349 +/- 245 for gamma delta T cells. By filter-spot ELISA (ELISPOT), the frequency of IFN-gamma-secreting gamma delta T cells was one-half of that of CD4+ T cells in response to M. tuberculosis, suggesting that gamma delta T cells on a per cell basis were more efficient producers of IFN-gamma than CD4+ T cells. In contrast, CD4+ T cells produced more IL-2 than gamma delta T cells, which correlated with diminished T cell proliferation of gamma delta T cells compared with CD4+ T cells. These results indicate that CD4+ and gamma delta T cell subsets have similar effector functions (cytotoxicity, IFN-gamma production) in response to M. tuberculosis-infected macrophages, despite differences in the Ags recognized, IL-2 production, and efficiency of IFN-gamma production.

An HI study of three long-tailed irregular galaxies in the cluster Abell 1367

We present the results on the distribution and kinematics of HI gas with higher sensitivity and in one case of higher spectral resolution as well than reported earlier, of three irregular galaxies CGCG 097073, 097079 and 097087 (UGC 06697) in the cluster Abell 1367. These galaxies are known to exhibit long (50 - 75 kpc) tails of radio continuum and optical emission lines (H alpha) pointing away from the cluster centre and arcs of starformation on the opposite sides of the tails, These features as well as the HI properties, with two of the galaxies (CGCG 097073 and 097079) exhibiting sharper gradients in HI intensity on the side of the tails, are consistent with the HI gas being affected by the ram pressure of the intracluster medium. However the HI emission in all the three galaxies extends to much smaller distances than the radio-continuum and H alpha tails, and are possibly still bound to the parent galaxies. Approximately 20 - 30 per cent of the HI mass is seen to accumulate on the downstream side due to the effects of ram pressure.