Beyond fractional derivatives: local approximation of other convolution integrals

Dynamic systems involving convolution integrals with decaying kernels, of which fractionally damped systems form a special case, are non-local in time and hence infinite dimensional. Straightforward numerical solution of such systems up to time t needs O(t(2)) computations owing to the repeated evaluation of integrals over intervals that grow like t. Finite-dimensional and local approximations are thus desirable. We present here an approximation method which first rewrites the evolution equation as a coupled in finite-dimensional system with no convolution, and then uses Galerkin approximation with finite elements to obtain linear, finite-dimensional, constant coefficient approximations for the convolution. This paper is a broad generalization, based on a new insight, of our prior work with fractional order derivatives (Singh & Chatterjee 2006 Nonlinear Dyn. 45, 183-206). In particular, the decaying kernels we can address are now generalized to the Laplace transforms of known functions; of these, the power law kernel of fractional order differentiation is a special case. The approximation can be refined easily. The local nature of the approximation allows numerical solution up to time t with O(t) computations. Examples with several different kernels show excellent performance. A key feature of our approach is that the dynamic system in which the convolution integral appears is itself approximated using another system, as distinct from numerically approximating just the solution for the given initial values; this allows non-standard uses of the approximation, e. g. in stability analyses.

Random network behaviour of protein structures

Geometric and structural constraints greatly restrict the selection of folds adapted by protein backbones, and yet, folded proteins show an astounding diversity in functionality. For structure to have any bearing on function, it is thus imperative that, apart from the protein backbone, other tunable degrees of freedom be accountable. Here, we focus on side-chain interactions, which non-covalently link amino acids in folded proteins to form a network structure. At a coarse-grained level, we show that the network conforms remarkably well to realizations of random graphs and displays associated percolation behavior. Thus, within the rigid framework of the protein backbone that restricts the structure space, the side-chain interactions exhibit an element of randomness, which account for the functional flexibility and diversity shown by proteins. However, at a finer level, the network exhibits deviations from these random graphs which, as we demonstrate for a few specific examples, reflect the intrinsic uniqueness in the structure and stability, and perhaps specificity in the functioning of biological proteins.

Growth inhibition of human promonocytic leukaemic U937 cells by interferon gamma is irreversible and not cell cycle phase-specific.

Growth of human promonocytic leukaemic U937 cells was found arrested within 24 h upon exposure to interferon gamma (IFN-gamma). Removal of the interferon did not result in the resumption of growth, as is evident from the absence of doubling of viable cell count and(3)H-thymidine incorporation. 5-Bromo-2'-deoxyuridine-based flow cytometric analysis of the growth-arrested cells, 24 h subsequent to the removal of IFN-gamma, showed absence of DNA synthesis, confirming the irreversible nature of the growth inhibition. Propidium iodide-based flow cytometric analysis of the growth-arrested cells showed a distribution which is typical of a growth inhibition without resulting in the accumulation of cells in any specific phase of the cell cycle. These results indicated that IFN-gamma arrested growth of U937 cells in an irreversible and cell cycle phase-independent manner. These observations were in contrast to our earlier report on the reversible and cell cycle phase-specific growth inhibition of human amniotic (fetal epithelial) WISH cells by the interferon. Copyright 1999 Academic Press.

Wigner distributions for finite-state systems without redundant phase-point operators

We set up Wigner distributions for N-state quantum systems following a Dirac-inspired approach. In contrast to much of the work in this study, requiring a 2N x 2N phase space, particularly when N is even, our approach is uniformly based on an N x N phase-space grid and thereby avoids the necessity of having to invoke a `quadrupled' phase space and hence the attendant redundance. Both N odd and even cases are analysed in detail and it is found that there are striking differences between the two. While the N odd case permits full implementation of the marginal property, the even case does so only in a restricted sense. This has the consequence that in the even case one is led to several equally good definitions of the Wigner distributions as opposed to the odd case where the choice turns out to be unique.

Random security scheme selection for mobile transactions

Security in a mobile communication environment is always a matter for concern, even after deploying many security techniques at device, network, and application levels. The end-to-end security for mobile applications can be made robust by developing dynamic schemes at application level which makes use of the existing security techniques varying in terms of space, time, and attacks complexities. In this paper we present a security techniques selection scheme for mobile transactions, called the Transactions-Based Security Scheme (TBSS). The TBSS uses intelligence to study, and analyzes the security implications of transactions under execution based on certain criterion such as user behaviors, transaction sensitivity levels, and credibility factors computed over the previous transactions by the users, network vulnerability, and device characteristics. The TBSS identifies a suitable level of security techniques from the repository, which consists of symmetric, and asymmetric types of security algorithms arranged in three complexity levels, covering various encryption/decryption techniques, digital signature schemes, andhashing techniques. From this identified level, one of the techniques is deployed randomly. The results shows that, there is a considerable reduction in security cost compared to static schemes, which employ pre-fixed security techniques to secure the transactions data.

Melting of an Anchored Bilayer: Phase Transitions in the Organic-Inorganic Hybrid Pervoskite (CH3NH3)(CH3(CH2)(n)NH3)(2)Pb2I7 (n=11, 13, 15, 17)

The conformation, organization, and phase transitions of alkyl chains in organic-inorganic hybrids based on the double pervoskite-slab lead iodides, (CH3NH3)(CH3(CH2)(n)NH3)(2)Pb2I7 (n = 11, 13, 15, 17) have been investigated by X-ray diffraction, calorimetry, and infrared vibrational spectroscopy. In these hybrid solids, double pervoskite (CH3NH3)Pb2I7 slabs are interleaved with alkyl ammonium chains with the anchored alkyl chains arranged as tilted bilayers and adopting a planar all-trans conformation at room temperature. The (CH3NH3)(CH3(CH2)(n)NH3)(2)Pb2I7 compounds exhibit a single reversible phase transition above room temperature with the associated enthalpy change varying linearly with alkyl chain length. This transition corresponds to the melting in two-dimensions of the alkyl chains of the anchored bilayer and is characterized by increased conformational disorder of the methylene units of the chain and loss of tilt angle coherence leading to an increase in the interslab spacing. By monitoring features in the infrared spectra that are characteristic of the global conformation of the alkyl chains, a quantitative relation between conformational disorder and melting of the anchored bilayer is established. It is found that, irrespective of the alkyl chain length, melting occurs when at least 60% of the chains in the anchored bilayer of (CH3NH3)(CH3(CH2)(n)NH3)(2)Pb2I7 have one or more gauche defects. This concentration is determined by the underlying lattice to which the alkyl chains are anchored.

Consistent quaternion interpolation for objective finite element approximation of geometrically exact beam

We explore an isoparametric interpolation of total quaternion for geometrically consistent, strain-objective and path-independent finite element solutions of the geometrically exact beam. This interpolation is a variant of the broader class known as slerp. The equivalence between the proposed interpolation and that of relative rotation is shown without any recourse to local bijection between quaternions and rotations. We show that, for a two-noded beam element, the use of relative rotation is not mandatory for attaining consistency cum objectivity and an appropriate interpolation of total rotation variables is sufficient. The interpolation of total quaternion, which is computationally more efficient than the one based on local rotations, converts nodal rotation vectors to quaternions and interpolates them in a manner consistent with the character of the rotation manifold. This interpolation, unlike the additive interpolation of total rotation, corresponds to a geodesic on the rotation manifold. For beam elements with more than two nodes, however, a consistent extension of the proposed quaternion interpolation is difficult. Alternatively, a quaternion-based procedure involving interpolation of relative rotations is proposed for such higher order elements. We also briefly discuss a strategy for the removal of possible singularity in the interpolation of quaternions, proposed in [I. Romero, The interpolation of rotations and its application to finite element models of geometrically exact rods, Comput. Mech. 34 (2004) 121–133]. The strain-objectivity and path-independence of solutions are justified theoretically and then demonstrated through numerical experiments. This study, being focused only on the interpolation of rotations, uses a standard finite element discretization, as adopted by Simo and Vu-Quoc [J.C. Simo, L. Vu-Quoc, A three-dimensional finite rod model part II: computational aspects, Comput. Methods Appl. Mech. Engrg. 58 (1986) 79–116]. The rotation update is achieved via quaternion multiplication followed by the extraction of the rotation vector. Nodal rotations are stored in terms of rotation vectors and no secondary storages are required.

Influence of magnetic clusters on electrical and magnetic properties of In1-xMnxSb/GaAs dilute magnetic semiconductor grown by liquid phase epitaxy

In1-xMnxSb films have been grown with different Mn doping concentrations (x = 0.0085, 0.018, 0.029 and 0.04) beyond the equilibrium 14 solubility limit by liquid phase epitaxy. We have studied temperature dependent resistivity, the Hall effect, magnetoresistance and magnetization for all compositions. Saturation in magnetization observed even at room temperature suggests the existence of ferromagnetic clusters in the film which has been verified by scanning electron microscopy studies. The anomalous Hall coefficient is found to be negative. Remnant field present on the surface of the clusters seems to affect the anomalous Hall effect at very low fields (below 350 Gauss). In the zero field resistivity, a variable-range hopping conduction mechanism dominates below 3.5 K for all samples above which activated behavior is predominant. The temperature dependence of the magnetization measurement shows a magnetic ordering below 10 K which is consistent with electrical measurements. (c) 2007 Elsevier Ltd. All rights reserved.

Localized reversible nanoscale phase separation in Pr0.63Ca0.37MnO3 single crystal using a scanning tunneling microscope tip

We report the destabilization of the charge ordered insulating (COI) state in a localized region of Pr0.63Ca0.37MnO3 single crystal by current injection using a scanning tunneling microscope tip. This leads to controlled phase separation and formation of localized metallic nanoislands in the COI matrix which have been detected by local tunneling conductance mapping. The metallic regions thus created persist even after reducing the injected current to lower values. The original conductance state can be restored by injecting a current of similar magnitude but of opposite polarity. We thus achieve reversible nanoscale phase separation that gives rise to the possibility to "write, read, and erase" nanosized conducting regions in an insulating matrix with high spatial resolution. (c) 2007 American Institute of Physics.

Bond-order wave phase, spin solitons, and thermodynamics of a frustrated linear spin-1/2 Heisenberg antiferromagnet

The linear spin-1/2 Heisenberg antiferromagnet with exchanges J(1) and J(2) between first and second neighbors has a bond-order wave (BOW) phase that starts at the fluid-dimer transition at J(2)/J(1)=0.2411 and is particularly simple at J(2)/J(1)=1/2. The BOW phase has a doubly degenerate singlet ground state, broken inversion symmetry, and a finite-energy gap E-m to the lowest-triplet state. The interval 0.4 < J(2)/J(1) < 1.0 has large E-m and small finite-size corrections. Exact solutions are presented up to N = 28 spins with either periodic or open boundary conditions and for thermodynamics up to N = 18. The elementary excitations of the BOW phase with large E-m are topological spin-1/2 solitons that separate BOWs with opposite phase in a regular array of spins. The molar spin susceptibility chi(M)(T) is exponentially small for T << E-m and increases nearly linearly with T to a broad maximum. J(1) and J(2) spin chains approximate the magnetic properties of the BOW phase of Hubbard-type models and provide a starting point for modeling alkali-tetracyanoquinodimethane salts.

Formation of ω-Al7Cu2Fe phase during laser processing of quasicrystal-forming Al-Cu-Fe alloy

The formation of an ω-Al7Cu2Fe phase during laser cladding of quasicrystal-forming Al65Cu23.3Fe11.7 alloy on a pure aluminium substrate is reported. This phase is found to nucleate at the periphery of primary icosahedral-phase particles. A large number of ω-phase particles form an envelope around the icosahedral phase. On the outer side, they form an interface with an agr-Al solid solution. Detailed transmission electron microscopic observations show that the ω phase exhibits an orientation relationship with the icosahedral phase. Analysis of experimental results suggests that the ω phase forms by precipitation on an icosahedral phase by heterogeneous nucleation and grows into the aluminium-rich melt until supersaturation is exhausted. The microstructural observations are explained in terms of available models of phase transformations.

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation-based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1-D as well as the 3-D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1-D problems.

Optimal Nonlinear Control and Estimation for a Reusable Launch Vehicle During Reentry Phase

In this paper a nonlinear optimal controller has been designed for aerodynamic control during the reentry phase of the Reusable Launch Vehicle (RLV). The controller has been designed based on a recently developed technique Optimal Dynamic Inversion (ODI). For full state feedback the controller has required full information about the system states. In this work an Extended Kalman filter (EKF) is developed to estimate the states. The vehicle (RLV) has been has been consider as a nonlinear Six-Degree-Of-Freedom (6-DOF) model. The simulation results shows that EKF gives a very good estimation of the states and it is working well with ODI. The resultant trajectories are very similar to those obtained by perfect state feedback using ODI only.

Network flow-control using asynchronous stochastic approximation

We propose several stochastic approximation implementations for related algorithms in flow-control of communication networks. First, a discrete-time implementation of Kelly's primal flow-control algorithm is proposed. Convergence with probability 1 is shown, even in the presence of communication delays and stochastic effects seen in link congestion indications. This ensues from an analysis of the flow-control algorithm using the asynchronous stochastic approximation (ASA) framework. Two relevant enhancements are then pursued: a) an implementation of the primal algorithm using second-order information, and b) an implementation where edge-routers rectify misbehaving flows. Next, discretetime implementations of Kelly's dual algorithm and primaldual algorithm are proposed. Simulation results a) verifying the proposed algorithms and, b) comparing the stability properties are presented.

New approximation algorithms for minimum cycle bases of graphs

We consider the problem of computing an approximate minimum cycle basis of an undirected edge-weighted graph G with m edges and n vertices; the extension to directed graphs is also discussed. In this problem, a {0,1} incidence vector is associated with each cycle and the vector space over F-2 generated by these vectors is the cycle space of G. A set of cycles is called a cycle basis of G if it forms a basis for its cycle space. A cycle basis where the sum of the weights of the cycles is minimum is called a minimum cycle basis of G. Cycle bases of low weight are useful in a number of contexts, e.g. the analysis of electrical networks, structural engineering, chemistry, and surface reconstruction. We present two new algorithms to compute an approximate minimum cycle basis. For any integer k >= 1, we give (2k - 1)-approximation algorithms with expected running time 0(kmn(1+2/k) + mn((1+1/k)(omega-1))) and deterministic running time 0(n(3+2/k)), respectively. Here omega is the best exponent of matrix multiplication. It is presently known that omega < 2.376. Both algorithms are o(m(omega)) for dense graphs. This is the first time that any algorithm which computes sparse cycle bases with a guarantee drops below the Theta(m(omega)) bound. We also present a 2-approximation algorithm with O(m(omega) root n log n) expected running time, a linear time 2-approximation algorithm for planar graphs and an O(n(3)) time 2.42-approximation algorithm for the complete Euclidean graph in the plane.