Nonadiabatic charge pumping by oscillating potentials in one dimension: Results for infinite system and finite ring

We study charge pumping when a combination of static potentials and potentials oscillating with a time period T is applied in a one-dimensional system of noninteracting electrons. We consider both an infinite system using the Dirac equation in the continuum approximation and a periodic ring with a finite number of sites using the tight-binding model. The infinite system is taken to be coupled to reservoirs on the two sides which are at the same chemical potential and temperature. We consider a model in which oscillating potentials help the electrons to access a transmission resonance produced by the static potentials and show that nonadiabatic pumping violates the simple sin phi rule which is obeyed by adiabatic two-site pumping. For the ring, we do not introduce any reservoirs, and we present a method for calculating the current averaged over an infinite time using the time evolution operator U(T) assuming a purely Hamiltonian evolution. We analytically show that the averaged current is zero if the Hamiltonian is real and time-reversal invariant. Numerical studies indicate another interesting result, namely, that the integrated current is zero for any time dependence of the potential if it is applied to only one site. Finally we study the effects of pumping at two sites on a ring at resonant and nonresonant frequencies, and show that the pumped current has different dependences on the pumping amplitude in the two cases.

Vision based indoor positioning in a retail environment

Modern smart phones often come with a significant amount of computational power and an integrated digital camera making them an ideal platform for intelligents assistants. This work is restricted to retail environments, where users could be provided with for example navigational in- structions to desired products or information about special offers within their close proximity. This kind of applications usually require information about the user's current location in the domain environment, which in our case corresponds to a retail store. We propose a vision based positioning approach that recognizes products the user's mobile phone's camera is currently pointing at. The products are related to locations within the store, which enables us to locate the user by pointing the mobile phone's camera to a group of products. The first step of our method is to extract meaningful features from digital images. We use the Scale- Invariant Feature Transform SIFT algorithm, which extracts features that are highly distinctive in the sense that they can be correctly matched against a large database of features from many images. We collect a comprehensive set of images from all meaningful locations within our domain and extract the SIFT features from each of these images. As the SIFT features are of high dimensionality and thus comparing individual features is infeasible, we apply the Bags of Keypoints method which creates a generic representation, visual category, from all features extracted from images taken from a specific location. A category for an unseen image can be deduced by extracting the corresponding SIFT features and by choosing the category that best fits the extracted features. We have applied the proposed method within a Finnish supermarket. We consider grocery shelves as categories which is a sufficient level of accuracy to help users navigate or to provide useful information about nearby products. We achieve a 40% accuracy which is quite low for commercial applications while significantly outperforming the random guess baseline. Our results suggest that the accuracy of the classification could be increased with a deeper analysis on the domain and by combining existing positioning methods with ours.

An accurate numerical integration scheme for finite rotations using rotation vector parametrization

A numerical integration procedure for rotational motion using a rotation vector parametrization is explored from an engineering perspective by using rudimentary vector analysis. The incremental rotation vector, angular velocity and acceleration correspond to different tangent spaces of the rotation manifold at different times and have a non-vectorial character. We rewrite the equation of motion in terms of vectors lying in the same tangent space, facilitating vector space operations consistent with the underlying geometric structure. While any integration algorithm (that works within a vector space setting) may be used, we presently employ a family of explicit Runge-Kutta algorithms to solve this equation. While this work is primarily motivated out of a need for highly accurate numerical solutions of dissipative rotational systems of engineering interest, we also compare the numerical performance of the present scheme with some of the invariant preserving schemes, namely ALGO-C1, STW, LIEMIDEA] and SUBCYC-M. Numerical results show better local accuracy via the present approach vis-a-vis the preserving algorithms. It is also noted that the preserving algorithms do not simultaneously preserve all constants of motion. We incorporate adaptive time-stepping within the present scheme and this in turn enables still higher accuracy and a `near preservation' of constants of motion over significantly longer intervals. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Euclideanization of Majorana and Weyl Fermions

A continuous procedure is presented for euclideanization of Majorana and Weyl fermions without doubling their degrees of freedom. The Euclidean theory so obtained is SO(4) invariant and Osterwalder-Schrader (OS) positive. This enables us to define a one-complex parameter family of the N=1 supersymmetric Yang-Mills (SSYM) theories which interpolate between the Minkowski and a Euclidean SSYM theory. The interpolating action, and hence the Euclidean action, manifests all the continous symmetries of the original Minkowski space theory.

Embracing Integrability in Stellar and Planetary Dynamics

Hamiltonian systems in stellar and planetary dynamics are typically near integrable. For example, Solar System planets are almost in two-body orbits, and in simulations of the Galaxy, the orbits of stars seem regular. For such systems, sophisticated numerical methods can be developed through integrable approximations. Following this theme, we discuss three distinct problems. We start by considering numerical integration techniques for planetary systems. Perturbation methods (that utilize the integrability of the two-body motion) are preferred over conventional "blind" integration schemes. We introduce perturbation methods formulated with Cartesian variables. In our numerical comparisons, these are superior to their conventional counterparts, but, by definition, lack the energy-preserving properties of symplectic integrators. However, they are exceptionally well suited for relatively short-term integrations in which moderately high positional accuracy is required. The next exercise falls into the category of stability questions in solar systems. Traditionally, the interest has been on the orbital stability of planets, which have been quantified, e.g., by Liapunov exponents. We offer a complementary aspect by considering the protective effect that massive gas giants, like Jupiter, can offer to Earth-like planets inside the habitable zone of a planetary system. Our method produces a single quantity, called the escape rate, which characterizes the system of giant planets. We obtain some interesting results by computing escape rates for the Solar System. Galaxy modelling is our third and final topic. Because of the sheer number of stars (about 10^11 in Milky Way) galaxies are often modelled as smooth potentials hosting distributions of stars. Unfortunately, only a handful of suitable potentials are integrable (harmonic oscillator, isochrone and Stäckel potential). This severely limits the possibilities of finding an integrable approximation for an observed galaxy. A solution to this problem is torus construction; a method for numerically creating a foliation of invariant phase-space tori corresponding to a given target Hamiltonian. Canonically, the invariant tori are constructed by deforming the tori of some existing integrable toy Hamiltonian. Our contribution is to demonstrate how this can be accomplished by using a Stäckel toy Hamiltonian in ellipsoidal coordinates.

Path-integral derivation of the superconformal anomaly for the N=1 supersymmetric Yang-Mills theory

Fujikawa's method of evaluating the supercurrent and the superconformal current anomalies, using the heat-kernel regularization scheme, is extended to theories with gauge invariance, in particular, to the off-shell N=1 supersymmetric Yang-Mills (SSYM) theory. The Jacobians of supersymmetry and superconformal transformations are finite. Although the gauge-fixing term is not supersymmetric and the regularization scheme is not manifestly supersymmetric, we find that the regularized Jacobians are gauge invariant and finite and they can be expressed in such a way that there is no one-loop supercurrent anomaly for the N=1 SSYM theory. The superconformal anomaly is nonzero and the anomaly agrees with a similar result obtained using other methods.

Complete Pick positivity and unitary invariance

The characteristic function for a contraction is a classical complete unitary invariant devised by Sz.-Nagy and Foias. Just as a contraction is related to the Szego kernel k(S) (z, w) = (1 - z (w) over tilde)(-1) for |z|, |w| < 1, by means of (1/k(S))(T,T*) >= 0, we consider an arbitrary open connected domain Omega in C-n, a complete Pick kernel k on Omega and a tuple T = (T-1, ..., T-n) of commuting bounded operators on a complex separable Hilbert space H such that (1/k)(T,T*) >= 0. For a complete Pick kernel the 1/k functional calculus makes sense in a beautiful way. It turns out that the model theory works very well and a characteristic function can be associated with T. Moreover, the characteristic function is then a complete unitary invariant for a suitable class of tuples T.

Thermal expansion of doped lanthanum gallates

Thermal expansion of several compositions of Sr and Mg-doped LaGaO3 including an A-site deficient composition (La0.9Sr0.1)(0.98)(Ga0.8Mg0.2)O-2.821 were measured in the temperature range from 298 to 1273 K. The effect of doping on thermal expansion was studied by varying the composition at one site of the perovskite structure (either A or B), while keeping the composition at the other site invariant. Thermal expansion varied nonlinearly with temperature and exhibited an inflexion between 550 and 620 K, probably related to the change in crystal structure from orthorhombic to rhombohedral. The dependence of average thermal expansion coefficient (alpha (av)) on the dopant concentration on either A or B site of the perovskite structure was found to be linear, when the composition at the other site was kept constant. Mg doping on the B-site had a greater effect on the average thermal expansion coefficient than Sr doping on the A-site. Cation deficiency at the A-site decreases thermal expansion when compositions at both sites are held constant.

Structural Relationships in 2,3-Bis-n-decyloxyanthracene and 12-Hydroxystearic Acid Molecular Gels and Aerogels Processed in Supercritical CO2

Supercritical carbon dioxide is used to prepare aerogels of two reference molecular organogelators, 2,3-bis-n-decyloxyanthracene (DDOA) (luminescent molecule) and 12-hydroxystearic acid (HSA). Electron microscopy reveals the fibrillar morphology of the aggregates generated by the protocol. SAXS and SANS measurements show that DDOA aerogels are crystalline materials exhibiting three morphs: (1) arrangements of the crystalline solid (2D p6m), (2) a second hexagonal morph slightly more compact, and (3) a packing specific of the fibers in the gel. Aggregates specific of the aerogel (volume fraction being typically phi approximate to 0.60) are developed over larger distances (similar to 1000 angstrom) and bear fewer defaults and residual strains than aggregates in the crystalline and gel phases. Porod, Scherrer and Debye-Bueche analyses of the scattering data have been performed. The first five diffraction peaks show small variations in position and intensity assigned to the variation of the number of fibers and their degree of vicinity within hexagonal bundles of the related SAFIN according to the Oster model. Conclusions are supported by the guidelines offered by the analysis of the situation in HSA aerogels for which the diffraction pattern can be described by two coexisting lamellar-like arrangements. The porosity of the aerogel, as measured by its specific surface extracted from the scattering invariant analysis, is only 1.8 times less than that of the swollen gel and is characteristic of a very porous material.

A structural principle for a class of adaptive systems

A feature common to many adaptive systems for identification and control is the adjustment.of gain parameters in a manner ensuring the stability of the overall system. This paper puts forward a principle which assures such a result for arbitrary systems which are linear and time invariant except for the adjustable parameters. The principle only demands that a transfer function be positive real. This transfer function dependent on the structure of the system with respect to the parameters. Several examples from adaptive identification, control and observer schemes are given as illustrations of the conceptual simplification provided by the structural principle.

A Bilinear Constitutive Model for Isotropic Bimodulus Materials

Proper formulation of stress-strain relations, particularly in tension-compression situations for isotropic biomodulus materials, is an unresolved problem. Ambartsumyan's model [8] and Jones' weighted compliance matrix model [9] do not satisfy the principle of coordinate invariance. Shapiro's first stress invariant model [10] is too simple a model to describe the behavior of real materials. In fact, Rigbi [13] has raised a question about the compatibility of bimodularity with isotropy in a solid. Medri [2] has opined that linear principal strain-principal stress relations are fictitious, and warned that the bilinear approximation of uniaxial stress-strain behavior leads to ill-working bimodulus material model under combined loading. In the present work, a general bilinear constitutive model has been presented and described in biaxial principal stress plane with zonewise linear principal strain-principal stress relations. Elastic coefficients in the model are characterized based on the signs of (i) principal stresses, (ii) principal strains, and (iii) on the value of strain energy component ratio ER greater than or less than unity. The last criterion is used in tension-compression and compression-tension situations to account for different shear moduli in pure shear stress and pure shear strain states as well as unequal cross compliances.

Adaptive optimal load balancing in a nonhomogeneous multiserver system with a central job scheduler

A model comprising several servers, each equipped with its own queue and with possibly different service speeds, is considered. Each server receives a dedicated arrival stream of jobs; there is also a stream of generic jobs that arrive to a job scheduler and can be individually allocated to any of the servers. It is shown that if the arrival streams are all Poisson and all jobs have the same exponentially distributed service requirements, the probabilistic splitting of the generic stream that minimizes the average job response time is such that it balances the server idle times in a weighted least-squares sense, where the weighting coefficients are related to the service speeds of the servers. The corresponding result holds for nonexponentially distributed service times if the service speeds are all equal. This result is used to develop adaptive quasi-static algorithms for allocating jobs in the generic arrival stream when the load parameters are unknown. The algorithms utilize server idle-time measurements which are sent periodically to the central job scheduler. A model is developed for these measurements, and the result mentioned is used to cast the problem into one of finding a projection of the root of an affine function, when only noisy values of the function can be observed

Control of Agent Swarms using Generalized Centroidal Cyclic Pursuit Laws

One of the major tasks in swarm intelligence is to design decentralized but homogenoeus strategies to enable controlling the behaviour of swarms of agents. It has been shown in the literature that the point of convergence and motion of a swarm of autonomous mobile agents can be controlled by using cyclic pursuit laws. In cyclic pursuit, there exists a predefined cyclic connection between agents and each agent pursues the next agent in the cycle. In this paper we generalize this idea to a case where an agent pursues a point which is the weighted average of the positions of the remaining agents. This point correspond to a particular pursuit sequence. Using this concept of centroidal cyclic pursuit, the behavior of the agents is analyzed such that, by suitably selecting the agents' gain, the rendezvous point of the agents can be controlled, directed linear motion of the agents can be achieved, and the trajectories of the agents can be changed by switching between the pursuit sequences keeping some of the behaviors of the agents invariant. Simulation experiments are given to support the analytical proofs.

Grad–Shafranov Reconstruction of Magnetic Clouds: Overview and Improvements

The Grad–Shafranov reconstruction is a method of estimating the orientation (invariant axis) and cross section of magnetic flux ropes using the data from a single spacecraft. It can be applied to various magnetic structures such as magnetic clouds (MCs) and flux ropes embedded in the magnetopause and in the solar wind. We develop a number of improvements of this technique and show some examples of the reconstruction procedure of interplanetary coronal mass ejections (ICMEs) observed at 1 AU by the STEREO, Wind, and ACE spacecraft during the minimum following Solar Cycle 23. The analysis is conducted not only for ideal localized ICME events but also for non-trivial cases of magnetic clouds in fast solar wind. The Grad–Shafranov reconstruction gives reasonable results for the sample events, although it possesses certain limitations, which need to be taken into account during the interpretation of the model results.

Analysis of sequence dependent variations in secondary and tertiary structure of tRNA molecules

The double helical regions of the five tRNA(Phe) and two tRNA(Asp) crystal structures have been analyzed using the local basepair step parameters. The sequence dependent effects in the mini double helices of tRNA are very similar to those observed in the crystal structures of oligonucleotides in the A-form, the purine-pyrimidine and purine-purine steps have small roll angles when compared to the fiber models of A-DNA as well as A-RNA, while the pyrimidine-purine doublet steps have large roll angles. The orientation of the basepairs in the D-stem is unusual and invariant i.e. they are different from the other three stems but are very similar in all the five tRNA(Phe) crystal structures, presumably due to tertiary interaction of the Watson-Crick basepairs with other bases, with all bases being highly conserved. The origin of the differences between the tertiary structures of tRNA(Phe) and tRNA(Asp) from yeast has also been investigated. It is found that even though the angle between the acceptor arm and the D-stem is very similar in the two structures, the angle subtended by the acceptor arm and the anticodon arm is smaller in the tRNA(Phe) structure (by more than 10 degrees). This is due to differences in the orientation of the two mini helices constituting the anticodon arm, which are inclined to each other by approximately 25 degrees in tRNA(Phe) and 16 degrees in tRNA(Asp). In addition, the acceptor arm, the D-stem and the anticodon stem are nearly coplanar in tRNA(Phe), while in tRNA(Asp) the anticodon stem projects out of the plane defined by the acceptor arm and the anticodon stem. These two features together lead to a larger separation between the acceptor and anticodon ends in tRNA(Asp) and indicate that the junction between the D-stem and the anticodon stem is quite variable, with features characteristic of a ball-and-socket type joint and determined for each tRNA molecule by the base sequence at the junction.