Universal simulation of Hamiltonian dynamics for quantum systems with finite-dimensional state spaces

What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum system? Dodd [Phys. Rev. A 65, 040301(R) (2002)] provided a partial solution to this problem in the form of an efficient algorithm to simulate any desired two-body Hamiltonian evolution using any fixed two-body entangling N-qubit Hamiltonian, and local unitaries. We extend this result to the case where the component systems are qudits, that is, have D dimensions. As a consequence we explain how universal quantum computation can be performed with any fixed two-body entangling N-qudit Hamiltonian, and local unitaries.

Factors affecting the rate of phosphocreatine resynthesis following intense exercise

Within the skeletal muscle cell at the onset of muscular contraction, phosphocreatine (PCr) represents the most immediate reserve for the rephosphorylation of adenosine triphosphate (ATP). As a result, its concentration can be reduced to less than 30% of resting levels during intense exercise. As a fall in the level of PCr appears to adversely affect muscle contraction, and therefore power output in a subsequent bout, maximising the rate of PCr resynthesis during a brief recovery period will be of benefit to an athlete involved in activities which demand intermittent exercise. Although this resynthesis process simply involves the rephosphorylation of creatine by aerobically produced ATP (with the release of protons), it has both a fast and slow component, each proceeding at a rate that is controlled by different components of the creatine kinase equilibrium. The initial fast phase appears to proceed at a rate independent of muscle pH. Instead, its rate appears to be controlled by adenosine diphosphate (ADP) levels; either directly through its free cytosolic concentration, or indirectly, through its effect on the free energy of ATP hydrolysis. Once this fast phase of recovery is complete, there is a secondary slower phase that appears almost certainly rate-dependant on the return of the muscle cell to homeostatic intracellular pH. Given the importance of oxidative phosphorylation in this resynthesis process, those individuals with an elevated aerobic power should be able to resynthesise PCr at a more rapid rate than their sedentary counterparts. However, results from studies that have used phosphorus nuclear magnetic resonance (P-31-NMR) spectroscopy, have been somewhat inconsistent with respect to the relationship between aerobic power and PCr recovery following intense exercise. Because of the methodological constraints that appear to have limited a number of these studies, further research in this area is warranted.

Are transitions in human gait determined by mechanical, kinetic or energetic factors?

It is currently unclear whether it is the need to maintain metabolic efficiency, the need to keep skeletal loading below critical force levels, or simple mechanical factors that drive the walk-to-run (W R) and run-to-walk (R-W) transitions in human gait. Eighteen adults (9 males and 9 females) locomoted on an instrumented treadmill using their preferred gait. Each completed 2 ascending (W-R) and 2 descending (R-W) series of trials under three levels of loading (0%, 15% and 30% body weight). For each trial, participants locomoted for 60 s at each of 9 different speeds -4 speeds both above and below their preferred transition speed (PTS) plus their PTS. Evidence was sought for critical levels of key kinetic (maximum vertical force, impulse, first peak force, time to first peak force and maximum loading rate), energetic (oxygen consumption, transport cost) and mechanical variables (limb lengths, strength) predictive of the gait transition. Analyses suggested the kinetic variables of time to first peak force and loading rate as the most likely determinants of the W-R and R-W transitions. (C) 2003 Elsevier Science B.V. All rights reserved.

Neural compensation for compliant loads during rhythmic movement

An experiment was performed to characterise the movement kinematics and the electromyogram (EMG) during rhythmic voluntary flexion and extension of the wrist against different compliant (elastic-viscous-inertial) loads. Three levels of each type of load, and an unloaded condition, were employed. The movements were paced at a frequency of I Hz by an auditory metronome, and visual feedback of wrist displacement in relation to a target amplitude of 100degrees was provided. Electro-myographic recordings were obtained from flexor carpi radialis (FCR) and extensor carpi radialis brevis (ECR). The movement profiles generated in the ten experimental conditions were indistinguishable, indicating that the CNS was able to compensate completely for the imposed changes in the task dynamics. When the level of viscous load was elevated, this compensation took the form of an increase in the rate of initial rise of the flexor and the extensor EMG burst. In response to increases in inertial load, the flexor and extensor EMG bursts commenced and terminated earlier in the movement cycle, and tended to be of greater duration. When the movements were performed in opposition to an elastic load, both the onset and offset of EMG activity occurred later than in the unloaded condition. There was also a net reduction in extensor burst duration with increases in elastic load, and an increase in the rate of initial rise of the extensor burst. Less pronounced alterations in the rate of initial rise of the flexor EMG burst were also observed. In all instances, increases in the magnitude of the external load led to elevations in the overall level of muscle activation. These data reveal that the elements of the central command that are modified in response to the imposition of a compliant load are contingent, not only upon the magnitude, but also upon the character of the load.

On the validity of thermogravimetric determination of carbon gasification kinetics

Thermogravimetric analysis has been widely applied in kinetic studies of carbon gasification, with the associated temporal weight change profiles being used to extract kinetic information and to validate gasification models. However the weight change profiles are not always governed by the intrinsic gasification activity because of the effect of chemisorption and its dynamics. In the present work we theoretically determine the criteria under which weight change profiles can be used to determine intrinsic kinetics for CO2 and O2 gasification by examining the region in which the chemisorption dynamics can be assumed pseudo-steady. It is found that the validity of the pseudo-steady assumption depends on the experimental conditions as well as on the initial surface area of carbon. Based on known mechanisms and rate constants an active surface area region is identified within which the steady state assumption is valid and the effect of chemisorption dynamics is negligible. The size of the permissible region is sensitive to the reaction temperature and gas pressure. The results indicate that in some cases the thermogravimetric data should be used with caution in kinetic studies. A large amount of literature on thermogravimetric analyzer determined char gasification kinetics is examined and the importance of chemisorption dynamics for the data assessed.

Predicting plant leaf area production: shoot assimilate accumulation and partitioning, and leaf area ratio, are stable for a wide range of sorghum population densities

Predicting plant leaf area production is required for modelling carbon balance and tiller dynamics in plant canopies. Plant leaf area production can be studied using a framework based on radiation intercepted, radiation use efficiency (RUE) and leaf area ratio (LAR) (ratio of leaf area to net above-ground biomass). The objective of this study was to test this framework for predicting leaf area production of sorghum during vegetative development by examining the stability of the contributing components over a large range of plant density. Four densities, varying from 2 to 16 plants m(-2), were implemented in a field experiment. Plants were either allowed to tiller or were maintained as uniculm by systematic tiller removal. In all cases, intercepted radiation was recorded daily and leaf area and shoot dry matter partitioning were quantified weekly at individual culm level. Up to anthesis, a unique relationship applied between fraction of intercepted radiation and leaf area index, and between shoot dry weight accumulation and amount of intercepted radiation, regardless of plant density. Partitioning of shoot assimilate between leaf, stem and head was also common across treatments up to anthesis, at both plant and culm levels. The relationship with thermal time (TT) from emergence of specific leaf area (SLA) and LAR of tillering plants did not change with plant density. In contrast, SLA of uniculm plants was appreciably lower under low-density conditions at any given TT from emergence. This was interpreted as a consequence of assimilate surplus arising from the inability of the plant to compensate by increasing the leaf area a culm could produce. It is argued that the stability of the extinction coefficient, RUE and plant LAR of tillering plants observed in these conditions provides a reliable way to predict leaf area production regardless of plant density. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.

Development of a generic crop model template in the cropping system model APSIM

The Agricultural Production Systems slMulator, APSIM, is a cropping system modelling environment that simulates the dynamics of soil-plant-management interactions within a single crop or a cropping system. Adaptation of previously developed crop models has resulted in multiple crop modules in APSIM, which have low scientific transparency and code efficiency. A generic crop model template (GCROP) has been developed to capture unifying physiological principles across crops (plant types) and to provide modular and efficient code for crop modelling. It comprises a standard crop interface to the APSIM engine, a generic crop model structure, a crop process library, and well-structured crop parameter files. The process library contains the major science underpinning the crop models and incorporates generic routines based on physiological principles for growth and development processes that are common across crops. It allows APSIM to simulate different crops using the same set of computer code. The generic model structure and parameter files provide an easy way to test, modify, exchange and compare modelling approaches at process level without necessitating changes in the code. The standard interface generalises the model inputs and outputs, and utilises a standard protocol to communicate with other APSIM modules through the APSIM engine. The crop template serves as a convenient means to test new insights and compare approaches to component modelling, while maintaining a focus on predictive capability. This paper describes and discusses the scientific basis, the design, implementation and future development of the crop template in APSIM. On this basis, we argue that the combination of good software engineering with sound crop science can enhance the rate of advance in crop modelling. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.

Surface diffusion of strong adsorbing vapours on porous carbon

Surface diffusion of strongly adsorbing hydrocarbon vapours on activated carbon was measured by using a constant molar flow method (D.D. Do, Dynamics of a semi-batch adsorber with constant molar supply rate: a method for studying adsorption rate of pure gas, Chem. Eng. Sci. 50 (1995) 549), where pure adsorbate is introduced into a semi-batch adsorber at a constant molar flow rate. The surface diffusivity was determined from the analysis of pressure response versus time, using a linear mathematical model developed earlier. To apply the linear theory over the non-linear range of the adsorption isotherm, we implement a differential increment method on the system which is initially equilibrated with some pre-determined loading. By conducting the experiments at different initial loadings, the surface diffusivity can be extracted as a function of loading. Propane, n-butane, n-hexane, benzene, and ethanol were used as diffusing adsorbate on a commercial activated carbon. It is found that the surface diffusivity of these strongly adsorbing vapours increases rapidly with loading, and the surface diffusion flux contributes significantly to the total flux and cannot be ignored. The surface diffusivity increases with temperature according to the Arrhenius law, and for the paraffins tested it decreases with the molecular weight of the adsorbate. (C) 2002 Elsevier Science Ltd. All rights reserved.

Control of ovulation with a GnRH agonist after superstimulation of follicular growth in buffalo: fertilization and embryo recovery

The potential to use a GnRH agonist bioimplant and injection of exogenous LH to control the time of ovulation in a multiple ovulation and embryo transfer (MOET) protocol was examined in buffalo. Mixed-parity buffalo (Bubalus bubalis; 4-15-year-old; 529 13 kg LW) were randomly assigned to one of five groups (n = 6): Group 1, conventional MOET protocol; Group 2, conventional MOET with 12 It delay in injection of PGF(2alpha); Group 3, implanted with GnRH agonist to block the pre-ovulatory surge release of LH; Group 4, implanted with GnRH agonist and injected with exogenous LH (Lutropin(R), 25 mg) 24 h after 4 days of superstimulation with FSH; Group 5, implanted with GnRH agonist and injected with LH 36 h after superstimulation with FSH. Ovarian follicular growth in all buffaloes was stimulated by treatment with FSH (Folltropin-V(R), 200 mg) administered over 4 days, and was monitored by ovarian ultrasonography. At the time of estrus, the number of follicles greater than or equal to8 mm. was greater (P < 0.05) for buffaloes in Group 2 (12.8) than for buffaloes in Groups 1 (8.5), 3 (7.3), 4 (6.1) and 5 (6.8), which did not differ. All buffaloes were mated by AI after spontaneous (Groups 1-3) or induced (Groups 4 and 5) ovulation. The respective number of buffalo that ovulated, number of corpora lutea, ovulation rate (%), and embryos + oocytes recovered were: Group 1 (2, 1.8 +/- 1.6, 18.0 +/- 13.6, 0.2 +/- 0.2); Group 2 (4, 6.1 +/- 2.9, 40.5 +/- 17.5, 3.7 +/- 2.1); Group 3 (0, 0, 0, 0); Group 4 (6, 4.3 +/- 1.2, 69.3 +/- 14.2, 2.0 +/- 0.9); and Group 5 (1, 2.5 +/- 2.5, 15.5 +/- 15.5, 2.1 +/- 2.1). All buffaloes in Group 4 ovulated after injection of LH and had a relatively high ovulation rate (69%) and embryo recovery (46%). It has been shown that the GnRH agonist-LH protocol can be used to improve the efficiency of MOET in buffalo. (C) 2002 Elsevier Science Inc. All rights reserved.

Mathematical modelling of schistosomiasis japonica: comparison of control strategies in the People's Republic of China

We present the first mathematical model on the transmission dynamics of Schistosoma japonicum. The work extends Barbour's classic model of schistosome transmission. It allows for the mammalian host heterogeneity characteristic of the S. japonicum life cycle, and solves the problem of under-specification of Barbour's model by the use of Chinese data we are collecting on human-bovine transmission in the Poyang Lake area of Jiangxi Province in China. The model predicts that in the lake/marshland areas of the Yangtze River basin: (1) once-early mass chemotherapy of humans is little better than twice-yearly mass chemotherapy in reducing human prevalence. Depending on the heterogeneity of prevalence within the population, targeted treatment of high prevalence groups, with lower overall coverage, can be more effective than mass treatment with higher overall coverage. Treatment confers a short term benefit only, with prevalence rising to endemic levels once chemotherapy programs are stopped (2) depending on the relative contributions of bovines and humans, bovine treatment can benefit humans almost as much as human treatment. Like human treatment, bovine treatment confers a short-term benefit. A combination of human and bovine treatment will dramatically reduce human prevalence and maintains the reduction for a longer period of time than treatment of a single host, although human prevalence rises once treatment ceases; (3) assuming 75% coverage of bovines, a bovine vaccine which acts on worm fecundity must have about 75% efficacy to reduce the reproduction rate below one and ensure mid-term reduction and long-term elimination of the parasite. Such a vaccination program should be accompanied by an initial period of human treatment to instigate a short-term reduction in prevalence, following which the reduction is enhanced by vaccine effects; (4) if the bovine vaccine is only 45% efficacious (the level of current prototype vaccines) it will lower the endemic prevalence, but will not result in elimination. If it is accompanied by an initial period of human treatment and by a 45% improvement in human sanitation or a 30% reduction in contaminated water contact by humans, elimination is then possible. (C) 2002 Elsevier Science B.V. All rights reserved.

Modelling activated sludge flocculation using population balances

Activated sludge flocculation was modelled using population balances. The model followed the dynamics of activated sludge flocculation providing a good approximation of the change in mean floe size with time. Increasing the average velocity gradient decreased the final floe size. The breakage rate coefficient and collision efficiency also varied with the average velocity gradient. A power law relationship was found for the increase in breakage rate coefficient with increasing average velocity gradient. Further investigation will be conducted to determine the relationship between the collision efficiency and particle size to provide a better approximation of dynamic changes in the floe size distribution during flocculation. (C) 2002 Published by Elsevier Science B.V.

Exchange rate sensitivity of Australian international equity funds

This paper analyses the exchange rate exposure displayed by a sample of Australian international equity trusts (IET). Exchange rate exposure is also examined in the context of differing economic climates with particular emphasis on the Asian crisis in mid-1997. It is found that there is evidence of exchange rate exposure particularly in the context of a multiple exchange rate model. Exposure varies substantially between three alternative time periods with different exposure apparent subsequent to the Asian crisis than prior to this event.