Computational chemistry on Fujitsu vector-parallel processors: Development and performance of applications software

In this and a preceding paper, we provide an introduction to the Fujitsu VPP range of vector-parallel supercomputers and to some of the computational chemistry software available for the VPP. Here, we consider the implementation and performance of seven popular chemistry application packages. The codes discussed range from classical molecular dynamics to semiempirical and ab initio quantum chemistry. All have evolved from sequential codes, and have typically been parallelised using a replicated data approach. As such they are well suited to the large-memory/fast-processor architecture of the VPP. For one code, CASTEP, a distributed-memory data-driven parallelisation scheme is presented. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

Inability of adult Limnodynastes peronii (Amphibia : Anura) to thermally acclimate locomotor performance

Despite several studies on adult amphibians, only larvae of the striped marsh frog (Limnodynastes peronii) have been reported to possess the ability to compensate for the effects of cool temperature on locomotor performance by thermal acclimation. In this study, we investigated whether this thermal acclimatory ability is shared by adult L. peronii. We exposed adult L. peronii to either 18 or 30 degrees C for 8 weeks and tested their swimming and jumping performance at six temperatures between 8 and 35 degrees C. Acute changes in temperature affected both maximum swimming and jumping performance, however there was no difference between the two treatment groups in locomotor performance between 8 and 30 degrees C. Maximum swimming velocity of both groups increased from 0.62 +/- 0.02 at 8 degrees C to 1.02 +/- 0.03 m s(-1) at 30 degrees C, while maximum jump distance increased from similar to 20 to > 60 cm over the same temperature range. Although adult L. peronii acclimated to 18 degrees C failed to produce a locomotor response at 35 degrees C, this most likely reflected a change in thermal tolerance limits with acclimation rather than modifications in the locomotor system. As all adult amphibians studied to date are incapable of thermally acclimating locomotor performance, including adults of L. peronii, this acclimatory capacity appears to be absent from the adult stage of development. (C) 2000 Elsevier Science Inc. All rights reserved.

Sensitivity to measurement perturbation of single-atom dynamics in cavity QED

We consider continuous observation of the nonlinear dynamics of single atom trapped in an optical cavity by a standing wave with intensity modulation. The motion of the atom changes the phase of the field which is then monitored by homodyne detection of the output field. We show that the conditional Hilbert space dynamics of this system, subject to measurement-induced perturbations, depends strongly on whether the corresponding classical dynamics is regular or chaotic. If the classical dynamics is chaotic, the distribution of conditional Hilbert space vectors corresponding to different observation records tends to be orthogonal. This is a characteristic feature of hypersensitivity to perturbation for quantum chaotic systems.

Quantification of free mycophenolic acid by high-performance liquid chromatography-atmospheric pressure chemical ionisation tandem mass spectrometry

To facilitate the investigation of free mycophenolic acid concentrations we developed a high-performance liquid chromatography tandem mass spectrometry method using indomethacin as an internal standard. Free drug was isolated from plasma samples (500 mul) using ultrafiltration, The analytes were extracted from the ultrafiltrate (200 mul) using C-18 solid-phase extraction. Detection was by selected reactant monitoring of mycophenolic acid (m/z 318.9-->190.9) and the internal standard (m/z 356.0-->297.1) with an atmospheric pressure chemical ionisation interface. The total chromatographic analysis time was 12 min. The method was found to be linear over the range investigated, 2.5-200 mug/l (r>0.990, n=6). The relative recovery of the method for the control samples studied (7.5, 40.0 and 150 mug/l) ranged from 95 to 104%. The imprecision of the method, expressed in terms of intra- and inter-day coefficients of variation, was

Atoms in an amplitude-modulated standing wave - dynamics and pathways to quantum chaos

Cold rubidium atoms are subjected to an amplitude-modulated far-detuned standing wave of light to form a quantum-driven pendulum. Here we discuss the dynamics of these atoms. Phase space resonances and chaotic transients of the system exhibit dynamics which can be useful in many atom optics applications as they can be utilized as means for phase space state preparation. We explain the occurrence of distinct peaks in the atomic momentum distribution, analyse them in detail and give evidence for the importance of the system for quantum chaos and decoherence studies.

Improved cardiac performance after ischemia in aged rats supplemented with vitamin E and alpha-lipoic acid

The purpose of these experiments was to examine the effects of dietary antioxidant supplementation with vitamin E (VE) and alpha -lipoic acid (alpha -LA) on biochemical and physiological responses to in vivo myocardial ischemia-reperfusion (I-R) in aged rats. Male Fischer-334 rats (18 mo old) were assigned to either 1) a control diet (CON) or 2) a VE and alpha -LA supplemented diet (ANTIOX). After a 14-wk feeding period, animals in each group underwent an in vivo I-R protocol (25 min of myocardial ischemia and 15 min of reperfusion). During reperfusion, peak arterial pressure was significantly higher (P < 0.05) in ANTIOX animals compared with CON diet animals. I-R resulted in a significant increase (P < 0.05) in myocardial lipid peroxidation in CON diet animals but not in ANTIOX animals. Compared with ANTIOX animals, heart homogenates from CON animals experienced significantly less (P < 0.05) oxidative damage when exposed to five different in vitro radical producing systems. These data indicate that dietary supplementation with VE and -LA protects the aged rat heart from I-R-induced lipid peroxidation by scavenging numerous reactive oxygen species. Importantly, this protection is associated with improved cardiac performance during reperfusion.

Context-free and context-sensitive dynamics in recurrent neural networks

Continuous-valued recurrent neural networks can learn mechanisms for processing context-free languages. The dynamics of such networks is usually based on damped oscillation around fixed points in state space and requires that the dynamical components are arranged in certain ways. It is shown that qualitatively similar dynamics with similar constraints hold for a(n)b(n)c(n), a context-sensitive language. The additional difficulty with a(n)b(n)c(n), compared with the context-free language a(n)b(n), consists of 'counting up' and 'counting down' letters simultaneously. The network solution is to oscillate in two principal dimensions, one for counting up and one for counting down. This study focuses on the dynamics employed by the sequential cascaded network, in contrast to the simple recurrent network, and the use of backpropagation through time. Found solutions generalize well beyond training data, however, learning is not reliable. The contribution of this study lies in demonstrating how the dynamics in recurrent neural networks that process context-free languages can also be employed in processing some context-sensitive languages (traditionally thought of as requiring additional computation resources). This continuity of mechanism between language classes contributes to our understanding of neural networks in modelling language learning and processing.

Simultaneous steady state and dynamic design of process systems using structural indicators

The simultaneous design of the steady-state and dynamic performance of a process has the ability to satisfy much more demanding dynamic performance criteria than the design of dynamics only by the connection of a control system. A method for designing process dynamics based on the use of a linearised systems' eigenvalues has been developed. The eigenvalues are associated with system states using the unit perturbation spectral resolution (UPSR), characterising the dynamics of each state. The design method uses a homotopy approach to determine a final design which satisfies both steady-state and dynamic performance criteria. A highly interacting single stage forced circulation evaporator system, including control loops, was designed by this method with the goal of reducing the time taken for the liquid composition to reach steady-state. Initially the system was successfully redesigned to speed up the eigenvalue associated with the liquid composition state, but this did not result in an improved startup performance. Further analysis showed that the integral action of the composition controller was the source of the limiting eigenvalue. Design changes made to speed up this eigenvalue did result in an improved startup performance. The proposed approach provides a structured way to address the design-control interface, giving significant insight into the dynamic behaviour of the system such that a systematic design or redesign of an existing system can be undertaken with confidence.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Controlling chaos in a Lorenz-like system using feedback

We demonstrate that the dynamics of an autonomous chaotic laser can be controlled to a periodic or steady state under self-synchronization. In general, past the chaos threshold the dependence of the laser output on feedback applied to the pump is submerged in the Lorenz-like chaotic pulsation. However there exist specific feedback delays that stabilize the chaos to periodic behavior or even steady state. The range of control depends critically on the feedback delay time and amplitude. Our experimental results are compared with the complex Lorenz equations which show good agreement.

Variable effects of larval size on post-metamorphic performance in the field

Larval quality may be capable of explaining much of the variation in the recruitment and subsequent population dynamics of benthic marine invertebrates. Whilst the effects of larval nutritional condition on adult performance have received the most attention, recent work has shown that larval size may also be an important and ubiquitous source of variation in larval quality. We examined the effects of variation in larval size on the post-metamorphic survival and growth of Watersipora subtorquata in 2 very different habitats - experimental substrata and pier pilings. We found strong effects of larval size on colony performance, although these varied among experiments. For colonies on experimental substrata, larval size positively affected adult survival and, initially, growth. However, after 3 wk in the field, there was no relationship between larval size and colony size, possibly because colonies were completely surrounded by newly settled organisms. Larval size also positively affected post-metamorphic growth of colonies on pier pilings, but, surprisingly, colonies that came from larger larvae had lower survival than colonies from smaller larvae. Overall, variation in larval size will strongly affect the recruitment and subsequent performance of adults in this species, although this may vary among different habitats. This study highlights the importance of examining the effects of larval quality on adult performance in as realistic conditions as possible, because of the strong interaction between larval size effects and the environment.

PREDBALB/c: a system for the prediction of peptide binding to H2(d) molecules, a haplotype of the BALB/c mouse

PREDBALB/c is a computational system that predicts peptides binding to the major histocompatibility complex-2 (H2(d)) of the BALB/c mouse, an important laboratory model organism. The predictions include the complete set of H2(d) class I ( H2-K-d, H2-L-d and H2-D-d) and class II (I-E-d and I-A(d)) molecules. The prediction system utilizes quantitative matrices, which were rigorously validated using experimentally determined binders and non-binders and also by in vivo studies using viral proteins. The prediction performance of PREDBALB/c is of very high accuracy. To our knowledge, this is the first online server for the prediction of peptides binding to a complete set of major histocompatibility complex molecules in a model organism (H2(d) haplotype). PREDBALB/c is available at http://antigen.i2r.a-star.edu.sg/predBalbc/.

Analysis of dynamic models for structural insight and model reduction .2. A multi-stage compressor shutdown case-study

In this second paper, the three structural measures which have been developed are used in the modelling of a three stage centrifugal synthesis gas compressor. The goal of this case study is to determine the essential mathematical structure which must be incorporated into the compressor model to accurately model the shutdown of this system. A simple, accurate and functional model of the system is created via three structural measures. It was found that the model can be correctly reduced into its basic modes and that the order of the differential system can be reduced from 51(st) to 20(th). Of the 31 differential equational 21 reduce to algebraic relations, 8 become constants and 2 can be deleted thereby increasing the algebraic set from 70 to 91 equations. An interpretation is also obtained as to which physical phenomena are dominating the dynamics of the compressor add whether the compressor will enter surge during the shutdown. Comparisons of the reduced model performance against the full model are given, showing the accuracy and applicability of the approach. Copyright (C) 1996 Elsevier Science Ltd

Quantum dynamics of an atomic Bose-Einstein condensate in a double-well potential

We consider the quantum dynamics of a neutral atom Bose-Einstein condensate in a double-well potential, including many-body hard-sphere interactions. Using a mean-field factorization we show that the coherent oscillations due to tunneling are suppressed when the number of atoms exceeds a critical value. An exact quantum solution, in a two-mode approximation, shows that the mean-field solution is modulated by a quantum collapse and revival sequence.

Structure optimization combining soft-core interaction functions, the diffusion equation method, and molecular dynamics

Smoothing the potential energy surface for structure optimization is a general and commonly applied strategy. We propose a combination of soft-core potential energy functions and a variation of the diffusion equation method to smooth potential energy surfaces, which is applicable to complex systems such as protein structures; The performance of the method was demonstrated by comparison with simulated annealing using the refinement of the undecapeptide Cyclosporin A as a test case. Simulations were repeated many times using different initial conditions and structures since the methods are heuristic and results are only meaningful in a statistical sense.