Reasoning about real-time repetitions: terminating and nonterminating

It is common for a real-time system to contain a nonterminating process monitoring an input and controlling an output. Hence, a real-time program development method needs to support nonterminating repetitions. In this paper we develop a general proof rule for reasoning about possibly nonterminating repetitions. The rule makes use of a Floyd-Hoare-style loop invariant that is maintained by each iteration of the repetition, a Jones-style relation between the pre- and post-states on each iteration, and a deadline specifying an upper bound on the starting time of each iteration. The general rule is proved correct with respect to a predicative semantics. In the case of a terminating repetition the rule reduces to the standard rule extended to handle real time. Other special cases include repetitions whose bodies are guaranteed to terminate, nonterminating repetitions with the constant true as a guard, and repetitions whose termination is guaranteed by the inclusion of a fixed deadline. (C) 2002 Elsevier Science B.V. All rights reserved.

Field quantization, photons and non-Hermitean modes

Field quantization in unstable optical systems is treated by expanding the vector potential in terms of non-Hermitean (Fox-Li) modes. We define non-Hermitean modes and their adjoints in both the cavity and external regions and make use of the important bi-orthogonality relationships that exist within each mode set. We employ a standard canonical quantization procedure involving the introduction of generalized coordinates and momenta for the electromagnetic (EM) field. Three-dimensional systems are treated, making use of the paraxial and monochromaticity approximations for the cavity non-Hermitean modes. We show that the quantum EM field is equivalent to a set of quantum harmonic oscillators (QHOs), associated with either the cavity or the external region non-Hermitean modes, and thus confirming the validity of the photon model in unstable optical systems. Unlike in the conventional (Hermitean mode) case, the annihilation and creation operators we define for each QHO are not Hermitean adjoints. It is shown that the quantum Hamiltonian for the EM field is the sum of non-commuting cavity and external region contributions, each of which can be expressed as a sum of independent QHO Hamiltonians for each non-Hermitean mode, except that the external field Hamiltonian also includes a coupling term responsible for external non-Hermitean mode photon exchange processes. The non-commutativity of certain cavity and external region annihilation and creation operators is associated with cavity energy gain and loss processes, and may be described in terms of surface integrals involving cavity and external region non-Hermitean mode functions on the cavity-external region boundary. Using the essential states approach and the rotating wave approximation, our results are applied to the spontaneous decay of a two-level atom inside an unstable cavity. We find that atomic transitions leading to cavity non-Hermitean mode photon absorption are associated with a different coupling constant to that for transitions leading to photon emission, a feature consequent on the use of non-Hermitean mode functions. We show that under certain conditions the spontaneous decay rate is enhanced by the Petermann factor.

Automatic construction of explicit R matrices for the one-parameter families of irreducible typical highest weight (0(m)vertical bar alpha(n)) representations of U-q[gl(m vertical bar n)]

We detail the automatic construction of R matrices corresponding to (the tensor products of) the (O-m\alpha(n)) families of highest-weight representations of the quantum superalgebras Uq[gl(m\n)]. These representations are irreducible, contain a free complex parameter a, and are 2(mn)-dimensional. Our R matrices are actually (sparse) rank 4 tensors, containing a total of 2(4mn) components, each of which is in general an algebraic expression in the two complex variables q and a. Although the constructions are straightforward, we describe them in full here, to fill a perceived gap in the literature. As the algorithms are generally impracticable for manual calculation, we have implemented the entire process in MATHEMATICA; illustrating our results with U-q [gl(3\1)]. (C) 2002 Published by Elsevier Science B.V.

Branching processes and computational collapse of discretized unimodal mappings

In computer simulations of smooth dynamical systems, the original phase space is replaced by machine arithmetic, which is a finite set. The resulting spatially discretized dynamical systems do not inherit all functional properties of the original systems, such as surjectivity and existence of absolutely continuous invariant measures. This can lead to computational collapse to fixed points or short cycles. The paper studies loss of such properties in spatial discretizations of dynamical systems induced by unimodal mappings of the unit interval. The problem reduces to studying set-valued negative semitrajectories of the discretized system. As the grid is refined, the asymptotic behavior of the cardinality structure of the semitrajectories follows probabilistic laws corresponding to a branching process. The transition probabilities of this process are explicitly calculated. These results are illustrated by the example of the discretized logistic mapping.

Multiple axon guidance cues establish the olfactory topographic map: how do these cues interact?

Each primary olfactory neuron stochastically expresses one of similar to1000 odorant receptors. The total population of these neurons therefore consists of similar to1,000 distinct subpopulations, each of which are mosaically dispersed throughout one of four semi-annular zones in the nasal cavity. The axons of these different subpopulations are initially intermingled within the olfactory nerve. However, upon reaching the olfactory bulb, they sort out and converge so that axons expressing the same odorant receptor typically target one or two glomeruli. The spatial location of each of these 1800 glomeruli are topographically-fixed in the olfactory bulb and are invariant from animal to animal. Thus, while odorant receptors are expressed mosaically by neurons throughout the olfactory neuroepithelium their axons sort out, converge and target the same glomerulus within the olfactory bulb. How is such precise and reproducible topographic targeting generated? While some of the mechanisms governing the growth cone guidance of olfactory sensory neurons are understood, the cues responsible for homing axons to their target site remain elusive.

The embryonic development of the bodywall and nervous system of the cestode flatworm Hymenolepis diminuta

Cestodes (tapeworms) are a derived, parasitic clade of the phylum Platyhelminthes (flatworms). The cestode body wall represents an adaptation to its endoparasitic lifestyle. The epidermis forms a nonciliated syncytium, and both muscular and nervous system are reduced. Morphological differences between cestodes and free-living flatworms become apparent already during early embryogenesis. Cestodes have a complex life cycle that begins with an infectious larva, called the oncosphere. In regard to cell number, cestode oncospheres are among the simplest multicellular organisms, containing in the order of 50-100 cells. As part of our continuing effort to analyze embryonic development in flatworms, we describe here the staining pattern obtained with acTub in embryos and larvae of the cestode Hymenolepis diminuta and, briefly, the monogenean Neoheterocotyle rhinobatidis. In addition, we labeled the embryonic musculature of Hymenolepis with phalloidin. In Hymenolepis embryos, two different cell types that we interpret as neurons and epidermal gland cells express acTub. There exist only two neurons that develop close to the midline at the anterior pole of the embryo. The axons of these two neurons project posteriorly into the center of the oncosphere, where they innervate the complex of muscles that is attached to the booklets. In addition to neurons, acTub labels a small and invariant set of epidermal gland cells that develop at superficial positions, anteriorly adjacent to the neurons, in the dorsal midline, and around the posteriorly located hooklets. During late stages of embryogenesis they spread and form a complete covering of the embryo. We discuss these data in the broader context of platyhelminth embryology.

New Zealand climate in the Neogene and implications for global atmospheric circulation

New Zealand has a good Neogene plant fossil record. During the Miocene it was without high topography and it was highly maritime, meaning that its climate, and the resulting vegetation, would be controlled dominantly by zonal climate conditions. Its vegetation record during this time suggests the climate passed from an ever-wet and cool but frostless phase in the Early Miocene in which Nothofagus subgenus Brassospora was prominent. Then it became seasonally dry, with vegetation in which palms and Eucalyptus were prominent and fires were frequent, and in the mid-Miocene, it developed a dry-climate vegetation dominated by Casuarinaceae. These changes are reflected in a sedimentological change from acidic to alkaline chemistry and the appearance of regular charcoal in the record. The vegetation then changed again to include a prominent herb component including Chenopodiaceae and Asteraceae. Sphagnum became prominent, and Nothofagus returned, but mainly as the subgenus Fuscospora (presently restricted to temperate climates). This is interpreted as a return to a generally wet, but now cold climate, in which outbreaks of cold polar air and frost were frequent. The transient drying out of a small maritime island and the accompanying vegetation/climate sequence could be explained by a higher frequency of the Sub-Tropical High Pressure (STHP) cells (the descending limbs of the Hadley cells) over New Zealand during the Miocene. This may have resulted from an increased frequency of 'blocking', a synoptic situation which occurs in the region today. An alternative hypothesis, that the global STHP belt lay at a significantly higher latitude in the early Neogene (perhaps 55degreesS) than today (about 30degreesS), is considered less likely because of physical constraints on STHP belt latitude. In either case, the difference between the early Neogene and present situation may have been a response to an increased polar-equatorial temperature gradient. This contrasts with current climate models for the geological past in which the latitude of the High Pressure belt impact is held invariant though geological time. (C) 2003 Elsevier Science B.V. All rights reserved.

Electron exchange coupling for single-donor solid-state spin qubits

Intervalley interference between degenerate conduction band minima has been shown to lead to oscillations in the exchange energy between neighboring phosphorus donor electron states in silicon [B. Koiller, X. Hu, and S. Das Sarma, Phys. Rev. Lett. 88, 027903 (2002); Phys. Rev. B 66, 115201 (2002)]. These same effects lead to an extreme sensitivity of the exchange energy on the relative orientation of the donor atoms, an issue of crucial importance in the construction of silicon-based spin quantum computers. In this article we calculate the donor electron exchange coupling as a function of donor position incorporating the full Bloch structure of the Kohn-Luttinger electron wave functions. It is found that due to the rapidly oscillating nature of the terms they produce, the periodic part of the Bloch functions can be safely ignored in the Heitler-London integrals as was done by Koiller, Hu, and Das Sarma, significantly reducing the complexity of calculations. We address issues of fabrication and calculate the expected exchange coupling between neighboring donors that have been implanted into the silicon substrate using an 15 keV ion beam in the so-called top down fabrication scheme for a Kane solid-state quantum computer. In addition, we calculate the exchange coupling as a function of the voltage bias on control gates used to manipulate the electron wave functions and implement quantum logic operations in the Kane proposal, and find that these gate biases can be used to both increase and decrease the magnitude of the exchange coupling between neighboring donor electrons. The zero-bias results reconfirm those previously obtained by Koiller, Hu, and Das Sarma.

Link invariants associated with gauge equivalent solutions of the Yang-Baxter equation: The one-parameter family of minimal typical representations of Uq[gl(2|1)]

In this paper we investigate the construction of state models for link invariants using representations of the braid group obtained from various gauge choices for a solution of the trigonometric Yang-Baxter equation. Our results show that it is possible to obtain invariants of regular isotopy (as defined by Kauffman) which may not be ambient isotopic. We illustrate our results with explicit computations using solutions of the trigonometric Yang-Baxter equation associated with the one-parameter family of minimal typical representations of the quantum superalgebra U-q,[gl(2/1)]. We have implemented MATHEMATICA code to evaluate the invariants for all prime knots up to 10 crossings.

A new wavelet-based adaptive method for solving population balance equations

A new wavelet-based adaptive framework for solving population balance equations (PBEs) is proposed in this work. The technique is general, powerful and efficient without the need for prior assumptions about the characteristics of the processes. Because there are steeply varying number densities across a size range, a new strategy is developed to select the optimal order of resolution and the collocation points based on an interpolating wavelet transform (IWT). The proposed technique has been tested for size-independent agglomeration, agglomeration with a linear summation kernel and agglomeration with a nonlinear kernel. In all cases, the predicted and analytical particle size distributions (PSDs) are in excellent agreement. Further work on the solution of the general population balance equations with nucleation, growth and agglomeration and the solution of steady-state population balance equations will be presented in this framework. (C) 2002 Elsevier Science B.V. All rights reserved.

Implementação de formulações do método dos elementos de contorno para associação de placas no espaço

Tese (doutorado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Mecânica, 2016.

Measuring and Controlling the Chaotic Motion of Profits

The study of economic systems has generated deep interest in exploring the complexity of chaotic motions in economy. Due to important developments in nonlinear dynamics, the last two decades have witnessed strong revival of interest in nonlinear endogenous business chaotic models. The inability to predict the behavior of dynamical systems in the presence of chaos suggests the application of chaos control methods, when we are more interested in obtaining regular behavior. In the present article, we study a specific economic model from the literature. More precisely, a system of three ordinary differential equations gather the variables of profits, reinvestments and financial flow of borrowings in the structure of a firm. Firstly, using results of symbolic dynamics, we characterize the topological entropy and the parameter space ordering of kneading sequences, associated with one-dimensional maps that reproduce significant aspects of the model dynamics. The analysis of the variation of this numerical invariant, in some realistic system parameter region, allows us to quantify and to distinguish different chaotic regimes. Finally, we show that complicated behavior arising from the chaotic firm model can be controlled without changing its original properties and the dynamics can be turned into the desired attracting time periodic motion (a stable steady state or into a regular cycle). The orbit stabilization is illustrated by the application of a feedback control technique initially developed by Romeiras et al. [1992]. This work provides another illustration of how our understanding of economic models can be enhanced by the theoretical and numerical investigation of nonlinear dynamical systems modeled by ordinary differential equations.

Lanczos algorithm on the grassmann manifold

5th. European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2008) 8th. World Congress on Computational Mechanics (WCCM8)

Genus and Braid Index Associated to Sequences of Renormalizable Lorenz Maps

We describe the Lorenz links generated by renormalizable Lorenz maps with reducible kneading invariant (K(f)(-), = K(f)(+)) = (X, Y) * (S, W) in terms of the links corresponding to each factor. This gives one new kind of operation that permits us to generate new knots and links from the ones corresponding to the factors of the *-product. Using this result we obtain explicit formulas for the genus and the braid index of this renormalizable Lorenz knots and links. Then we obtain explicit formulas for sequences of these invariants, associated to sequences of renormalizable Lorenz maps with kneading invariant (X, Y) * (S,W)*(n), concluding that both grow exponentially. This is specially relevant, since it is known that topological entropy is constant on the archipelagoes of renormalization.

New physics contributions to A(FB)(t(t)over-bar) at the Tevatron

The Tevatron has measured a discrepancy relative to the standard model prediction in the forward-backward asymmetry in top quark pair production. This asymmetry grows with the rapidity difference of the two top quarks. It also increases with the invariant mass of the t (t) over bar pair, reaching, for high invariant masses, 3.4 standard deviations above the next-to-leading order prediction for the charge asymmetry of QCD. However, perfect agreement between experiment and the standard model was found in both total and differential cross section of top quark pair production. As this result could be a sign of new physics we have parametrized this new physics in terms of a complete set of dimension six operators involving the top quark. We have then used a Markov chain Monte Carlo approach in order to find the best set of parameters that fits the data, using all available data regarding top quark pair production at the Tevatron. We have found that just a very small number of operators are able to fit the data better than the standard model.