A class of self-stabilizing MCA learning algorithms

In this letter, we propose a class of self-stabilizing learning algorithms for minor component analysis (MCA), which includes a few well-known MCA learning algorithms. Self-stabilizing means that the sign of the weight vector length change is independent of the presented input vector. For these algorithms, rigorous global convergence proof is given and the convergence rate is also discussed. By combining the positive properties of these algorithms, a new learning algorithm is proposed which can improve the performance. Simulations are employed to confirm our theoretical results.

A modified discrete random pore model allowing for different initial surface reactivity

We investigate here a modification of the discrete random pore model [Bhatia SK, Vartak BJ, Carbon 1996;34:1383], by including an additional rate constant which takes into account the different reactivity of the initial pore surface having attached functional groups and hydrogens, relative to the subsequently exposed surface. It is observed that the relative initial reactivity has a significant effect on the conversion and structural evolution, underscoring the importance of initial surface chemistry. The model is tested against experimental data on chemically controlled char oxidation and steam gasification at various temperatures. It is seen that the variations of the reaction rate and surface area with conversion are better represented by the present approach than earlier random pore models. The results clearly indicate the improvement of model predictions in the low conversion region, where the effect of the initially attached functional groups and hydrogens is more significant, particularly for char oxidation. It is also seen that, for the data examined, the initial surface chemistry is less important for steam gasification as compared to the oxidation reaction. Further development of the approach must also incorporate the dynamics of surface complexation, which is not considered here.

SPAnDX: a process-based population dynamics model to explore management and climate change impacts on an invasive alien plant, Acacia nilotica

This paper describes a process-based metapopulation dynamics and phenology model of prickly acacia, Acacia nilotica, an invasive alien species in Australia. The model, SPAnDX, describes the interactions between riparian and upland sub-populations of A. nilotica within livestock paddocks, including the effects of extrinsic factors such as temperature, soil moisture availability and atmospheric concentrations of carbon dioxide. The model includes the effects of management events such as changing the livestock species or stocking rate, applying fire, and herbicide application. The predicted population behaviour of A. nilotica was sensitive to climate. Using 35 years daily weather datasets for five representative sites spanning the range of conditions that A. nilotica is found in Australia, the model predicted biomass levels that closely accord with expected values at each site. SPAnDX can be used as a decision-support tool in integrated weed management, and to explore the sensitivity of cultural management practices to climate change throughout the range of A. nilotica. The cohort-based DYMEX modelling package used to build and run SPAnDX provided several advantages over more traditional population modelling approaches (e.g. an appropriate specific formalism (discrete time, cohort-based, process-oriented), user-friendly graphical environment, extensible library of reusable components, and useful and flexible input/output support framework). (C) 2003 Published by Elsevier Science B.V.

A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics

The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to Moore's law, i.e., to double every 18-24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures.

Neuromuscular-skeletal constraints on the acquisition of skill in a discrete torque production task

The organisation of the human neuromuscular-skeletal system allows an extremely wide variety of actions to be performed, often with great dexterity. Adaptations associated with skill acquisition occur at all levels of the neuromuscular-skeletal system although all neural adaptations are inevitably constrained by the organisation of the actuating apparatus (muscles and bones). We quantified the extent to which skill acquisition in an isometric task set is influenced by the mechanical properties of the muscles used to produce the required actions. Initial performance was greatly dependent upon the specific combination of torques required in each variant of the experimental task. Five consecutive days of practice improved the performance to a similar degree across eight actions despite differences in the torques required about the elbow and forearm. The proportional improvement in performance was also similar when the actions were performed at either 20 or 40% of participants' maximum voluntary torque capacity. The skill acquired during practice was successfully extrapolated to variants of the task requiring more torque than that required during practice. We conclude that while the extent to which skill can be acquired in isometric actions is independent of the specific combination of joint torques required for target acquisition, the nature of the kinetic adaptations leading to the performance improvement in isometric actions is influenced by the neural and mechanical properties of the actuating muscles.

Oscillatory regulation of hes1: Discrete stochastic delay modelling and simulation

Discrete stochastic simulations are a powerful tool for understanding the dynamics of chemical kinetics when there are small-to-moderate numbers of certain molecular species. In this paper we introduce delays into the stochastic simulation algorithm, thus mimicking delays associated with transcription and translation. We then show that this process may well explain more faithfully than continuous deterministic models the observed sustained oscillations in expression levels of hes1 mRNA and Hes1 protein.

Scaling and non-scaling of muscle activity, kinematics, and dynamics in sit-ups with different degrees of difficulty

The purpose of this study was to investigate how the CNS adjusts motor patterns for variants of a complex axial movement-the situp. Adjustments were induced by changing the support surface contact and mass distribution of the body. Healthy adults performed straight-legged sit-ups, 3 s in duration, with support added to or removed from the lumbar trunk, or with mass added to the head or to the legs. Each of these interventions either increased or decreased the difficulty of the task. The study addressed the extent to which changes in sit-up difficulty are compensated by scaling of muscle activity, kinematics, and dynamics versus the extent to which they are compensated by changing discretely the motor pattern. The analysis of muscle activity, kinematics, and dynamics focused on the first 30-40% of the sit-up-the trunk flexion phase-since this is the most critical part of the movement. Our results demonstrate that, in some respects, sit-up kinematics and dynamics scaled with difficulty, but in other respects, they did not. Muscle activity also scaled, in many respects, but in more difficult sit-ups, abdominal flexor activity decreased instead of increased. Non-scaling changes in these parameters suggest that complex movements, such as the sit-up, may require discrete changes in motor pattern in order to deal with large loads, which challenge the available leverage. (C) 2005 Elsevier Ltd. All rights reserved.

Evolving noisy oscillatory dynamics in genetic regulatory networks

We introduce a genetic programming (GP) approach for evolving genetic networks that demonstrate desired dynamics when simulated as a discrete stochastic process. Our representation of genetic networks is based on a biochemical reaction model including key elements such as transcription, translation and post-translational modifications. The stochastic, reaction-based GP system is similar but not identical with algorithmic chemistries. We evolved genetic networks with noisy oscillatory dynamics. The results show the practicality of evolving particular dynamics in gene regulatory networks when modelled with intrinsic noise.

Time-reversal test for stochastic quantum dynamics

The calculation of quantum dynamics is currently a central issue in theoretical physics, with diverse applications ranging from ultracold atomic Bose-Einstein condensates to condensed matter, biology, and even astrophysics. Here we demonstrate a conceptually simple method of determining the regime of validity of stochastic simulations of unitary quantum dynamics by employing a time-reversal test. We apply this test to a simulation of the evolution of a quantum anharmonic oscillator with up to 6.022×1023 (Avogadro's number) of particles. This system is realizable as a Bose-Einstein condensate in an optical lattice, for which the time-reversal procedure could be implemented experimentally.

Superchemistry: dynamics of coupled atomic and molecular Bose condensates

We analyze the dynamics of a dilute, trapped Bose-condensed atomic gas coupled to a diatomic molecular Bose gas by coherent Raman transitions. This system is shown to result in a new type of “superchemistry,” in which giant collective oscillations between the atomic and the molecular gas can occur. The phenomenon is caused by stimulated emission of bosonic atoms or molecules into their condensate phases.

Quantum dynamics of three coupled atomic Bose-Einstein condensates

The simplest model of three coupled Bose-Einstein condensates is investigated using a group theoretical method. The stationary solutions are determined using the SU(3) group under the mean-field approximation. This semiclassical analysis, using system symmetries, shows a transition in the dynamics of the system from self trapping to delocalization at a critical value for the coupling between the condensates. The global dynamics are investigated by examination of the stable points, and our analysis shows that the structure of the stable points depends on the ratio of the condensate coupling to the particle-particle interaction, and undergoes bifurcations as this ratio is varied. This semiclassical model is compared to a full quantum treatment, which also displays a dynamical transition. The quantum case has collapse and revival sequences superimposed on the semiclassical dynamics, reflecting the underlying discreteness of the spectrum. Nonzero circular current states are also demonstrated as one of the higher-dimensional effects displayed in this system.

Substructure and Dynamics of the Fornax Cluster

We present the first dynamical analysis of a galaxy cluster to include a large fraction of dwarf galaxies. Our sample of 108 Fornax Cluster members measured with the UK Schmidt Telescope FLAIR-II spectrograph contains 55 dwarf galaxies (15.5 > b(j) > 18.0 or -16 > M-B > -13.5). H alpha emission shows that of the dwarfs are star forming, twice the fraction implied by morphological classifications. The total sample has a mean velocity of 1493 +/- 36 kms s(-1) and a velocity dispersion of 374 +/- 26 km s(-1). The dwarf galaxies form a distinct population: their velocity dispersion (429 +/- 41 km s(-1)) is larger than that of the giants () at the 98% confidence level. This suggests that the dwarf population is dominated by infalling objects whereas the giants are virialized. The Fornax system has two components, the main Fornax Cluster centered on NGC 1399 with cz = 1478 km s(-1) and sigma (cz) = 370 km s(-1) and a subcluster centered 3 degrees to the southwest including NGC 1316 with cz = 1583 km s(-1) and sigma (cz) = 377 km s(-1). This partition is preferred over a single cluster at the 99% confidence level. The subcluster, a site of intense star formation, is bound to Fornax and probably infalling toward the cluster core for the first time. We discuss the implications of this substructure for distance estimates of the Fornax Cluster. We determine the cluster mass profile using the method of Diaferio, which does not assume a virialized sample. The mass within a projected radius of 1.4 Mpc is (7 +/- 2) x 10(13) M-., and the mass-to-light ratio is 300 +/- 100 M-./L-.. The mass is consistent with values derived from the projected mass virial estimator and X-ray measurements at smaller radii.

Evolution of Wolbachia pipientis transmission dynamics in insects

Wolbachia pipientis is an intracellular bacterial parasite of arthropods that enhances its transmission by manipulating host reproduction, most commonly by inducing cytoplasmic incompatibility. The discovery of isolates with modified cytoplasmic incompatibility phenotypes and others with novel virulence properties is an indication of the potential breadth of evolutionary strategies employed by Wolbachia.

Dependence of transient dynamics in a class-C laser upon variation of inversion with time

The transient statistics of a gain-switched coherently pumped class-C laser displays a linear correlation between the first passage time and subsequent peak intensity. Measurements are reported showing a positive or negative sign of this linear correlation, controlled through the switching time and the laser detuning. Further measurements of the small-signal laser gain combined with calculations involving a three-level laser model indicate that this sign fundamentally depends upon the way the laser inversion varies during the gain switching, despite the added dynamics of the laser polarization in the class-C laser. [S1050-2947(97)07112-6].