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.

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.

Characterization of pore size distributions of mesoporous materials from adsorption isotherms

In this article, a new hybrid model for estimating the pore size distribution of micro- and mesoporous materials is developed, and tested with the adsorption data of nitrogen, oxygen, and argon on ordered mesoporous materials reported in the literature. For the micropore region, the model uses the Dubinin-Rudushkevich (DR) isotherm with the Chen-Yang modification. A recent isotherm model of the authors for nonporous materials, which uses a continuum-mechanical model for the multilayer region and the Unilan model for the submonolayer region, has been extended for adsorption in mesopores. The experimental data is inverted using regularization to obtain the pore size distribution. The present model was found to be successful in predicting the pore size distribution of pure as well as binary physical mixtures of MCM-41 synthesized with different templates, with results in agreement with those from the XRD method and nonlocal density functional theory. It was found that various other recent methods, as well as the classical Broekhoff and de Beer method, underpredict the pore diameter of MCM-41. The present model has been successfully applied to MCM-48, SBA's, CMK, KIT, HMS, FSM, MTS, mesoporous fly ash, and a large number of other regular mesoporous materials.

Quantitative analysis of vascular to cardiac selectivity of L- and T-type voltage-operated calcium channel antagonists in human tissues

1. Classical L-type voltage-operated calcium channel (VOCC) antagonists dilate blood vessels, depress myocardial contractility and slow cardiac conduction. 2. We compared four L-type VOCC antagonists and a novel tetralol derivative, mibefradil, reportedly 10-fold more selective for T- (transient) over L-type VOCC in two in vitro assays of human tissue, namely isolated small arteries from the aortic vasa vasorum in a myograph and right atrial trabeculae muscle under isometric force conditions. 3. In arteries contracted with K+ (62 mmol/L), the relaxation pIC(50) values for the VOCC antagonists felodipine, nifedipine, amlodipine, verapamil and mibefradil were 8.30, 7.78, 6.64, 6.26 and 6.22, respectively. In atrial trabeculae, the pIC(50) values to inhibit the inotropic response to a submaximal concentration of isoprenaline (6 nmol/L) for felodipine, nifedipine, verapamil, amlodipine and mibefradil were 7.21, 6.95, 6.91, 5.94 and 4.61, respectively. 4. Taking the anti-log (pIC(50) vessel - pIC(50) atrium) the vascular relaxation to cardiac depression potency ratios for mibefradil, felodipine, nifedipine, amlodipine and verapamil were 41, 12, 7, 5 and 0.22, respectively. 5. We conclude that, in human tissue assays, perhaps T- over L-type VOCC selectivity confers the most favourable vascular selectivity on mibefradil. Alternatively, splice variants of L-type VOCC in the vasculature (CaV1.2b) may be more sensitive to mibefradil than the splice variants in the heart (CaV1.2a).

Vortices and solitons in Bose-condensates

The new science of nonlinear atom optics and atom lasers is evolving rapidly. There are similarities between many related areas in modern photonic and atom optics, particularly at the mean-field level. In both cases we can often use classical nonlinear wave equations to describe classical solitons, vortices, and other nonlinear structure. Atom-molecular coupling can be used to play the role of second-harmonic generation. This leads to novel types of soliton. In addition, quantum effects at low densities are likely to be readily observable.

Quantum chaos in atom optics

I shall discuss the quantum and classical dynamics of a class of nonlinear Hamiltonian systems. The discussion will be restricted to systems with one degree of freedom. Such systems cannot exhibit chaos, unless the Hamiltonians are time dependent. Thus we shall consider systems with a potential function that has a higher than quadratic dependence on the position and, furthermore, we shall allow the potential function to be a periodic function of time. This is the simplest class of Hamiltonian system that can exhibit chaotic dynamics. I shall show how such systems can be realized in atom optics, where very cord atoms interact with optical dipole potentials of a far-off resonance laser. Such systems are ideal for quantum chaos studies as (i) the energy of the atom is small and action scales are of the order of Planck's constant, (ii) the systems are almost perfectly isolated from the decohering effects of the environment and (iii) optical methods enable exquisite time dependent control of the mechanical potentials seen by the atoms.

Quantum measurement and stochastic processes in mesoscopic conductors

In quantum measurement theory it is necessary to show how a, quantum source conditions a classical stochastic record of measured results. We discuss mesoscopic conductance using quantum stochastic calculus to elucidate the quantum nature of the measurement taking place in these systems. To illustrate the method we derive the current fluctuations in a two terminal mesoscopic circuit with two tunnel barriers containing a single quasi bound state on the well. The method enables us to focus on either the incoming/ outgoing Fermi fields in the leads, or on the irreversible dynamics of the well state itself. We show an equivalence between the approach of Buttiker and the Fermi quantum stochastic calculus for mesoscopic systems.

Pathway-specific targeting of GAB(A) a receptor subtypes to somatic and dendritic synapses in the central amygdala

Neurons in the central amygdala express two distinct types of ionotropic GABA receptor. One is the classical GABA(A) receptor that is blocked by low concentrations of bicuculline and positively modulated by benzodiazepines. The other is a novel type of ionotropic GABA receptor that is less sensitive to bicuculline but blocked by the GABA(C) receptor antagonist (1,2,5,6-tetrohydropyridine-4-yl) methylphosphinic acid (TPMPA) and by benzodiazepines. In this study, we examine the distribution of these two receptor types. Recordings of GABAergic miniature inhibitory postsynaptic currents (mIPSCs) showed a wide variation in amplitude. Most events had amplitudes of 100 pA. Large-amplitude events also had rise times faster than small-amplitude events. Large-amplitude events were fully blocked by 10 muM bicuculline but unaffected by TPMPA. Small amplitude events were partially blocked by both bicuculline and TPMPA. Focal application of hypertonic sucrose to the soma evoked large-amplitude mIPSCs, whereas focal dendritic application of sucrose evoked small-amplitude mIPSCs. Thus inhibitory synapses on the dendrites of neurons in the central amygdala express both types of GABA receptor, but somatic synapses expressed purely GABA(A) receptors. Minimal stimulation revealed that inhibitory inputs arising from the laterally located intercalated cells innervate dendritic synapses, whereas inhibitory inputs of medial origin innervated somatic inhibitory synapses. These results show that different types of ionotropic GABA receptors are targeted to spatially and functionally distinct synapses. Thus benzodiazepines will have different modulatory effects on different inhibitory pathways in the central amygdala.

Transfer matrix eigenvalues of the anisotropic multiparametric U model

A multiparametric extension of the anisotropic U model is discussed which maintains integrability. The R-matrix solving the Yang-Baxter equation is obtained through a twisting construction applied to the underlying U-q(sl (2/1)) superalgebraic structure which introduces the additional free parameters that arise in the model. Three forms of Bethe ansatz solution for the transfer matrix eigenvalues are given which we show to be equivalent.

The advantages of theft over toil: The design inference and arguing from ignorance

Intelligent design theorist William Dembski has proposed an explanatory filter for distinguishing between events due to chance, lawful regularity or design. We show that if Dembski's filter were adopted as a scientific heuristic, some classical developments in science would not be rational, and that Dembski's assertion that the filter reliably identifies rarefied design requires ignoring the state of background knowledge. If background information changes even slightly, the filter's conclusion will vary wildly. Dembski fails to overcome Hume's objections to arguments from design.

Group theory and quasiprobability integrals of Wigner functions

The integral of the Wigner function of a quantum-mechanical system over a region or its boundary in the classical phase plane, is called a quasiprobability integral. Unlike a true probability integral, its value may lie outside the interval [0, 1]. It is characterized by a corresponding selfadjoint operator, to be called a region or contour operator as appropriate, which is determined by the characteristic function of that region or contour. The spectral problem is studied for commuting families of region and contour operators associated with concentric discs and circles of given radius a. Their respective eigenvalues are determined as functions of a, in terms of the Gauss-Laguerre polynomials. These polynomials provide a basis of vectors in a Hilbert space carrying the positive discrete series representation of the algebra su(1, 1) approximate to so(2, 1). The explicit relation between the spectra of operators associated with discs and circles with proportional radii, is given in terms of the discrete variable Meixner polynomials.

Characterisation of the three-dimensional structure of earthworm burrow systems using image analysis and mathematical morphology

The aim of this work was to exemplify the specific contribution of both two- and three-dimensional (31)) X-ray computed tomography to characterise earthworm burrow systems. To achieve this purpose we used 3D mathematical morphology operators to characterise burrow systems resulting from the activity of an anecic (Aporrectodea noctunia), and an endogeic species (Allolobophora chlorotica), when both species were introduced either separately or together into artificial soil cores. Images of these soil cores were obtained using a medical X-ray tomography scanner. Three-dimensional reconstructions of burrow systems were obtained using a specifically developed segmentation algorithm. To study the differences between burrow systems, a set of classical tools of mathematical morphology (granulometries) were used. So-called granulometries based on different structuring elements clearly separated the different burrow systems. They enabled us to show that burrows made by the anecic species were fatter, longer, more vertical, more continuous but less sinuous than burrows of the endogeic species. The granulometry transform of the soil matrix showed that burrows made by A. nocturna were more evenly distributed than those of A. chlorotica. Although a good discrimination was possible when only one species was introduced into the soil cores, it was not possible to separate burrows of the two species from each other in cases where species were introduced into the same soil core. This limitation, partly due to the insufficient spatial resolution of the medical scanner, precluded the use of the morphological operators to study putative interactions between the two species.

Operational regimes for a simplified one-step artificial chaperone refolding method

The artificial chaperone method for protein refolding developed by Rozema et al. (Rozema, D.; Gellman, S. H. J. Am. Chem. Soc. 1995, 117 (8), 2373-2374) involves the sequential dilution of denatured protein into a buffer containing detergent (cetyltrimethylammonium bromide, CTAB) and then into a refolding buffer containing cyclodextrin WD). In this paper a simplified one-step artificial chaperone method is reported, whereby CTAB is added directly to the denatured solution, which is then diluted directly into a refolding buffer containing P-cyclodextrin (P-CD). This new method can be applied at high protein concentrations, resulting in smaller processing volumes and a more concentrated protein solution following refolding. The increase in achievable protein concentration results from the enhanced solubility of CTAB at elevated temperatures in concentrated denaturant. The refolding yields obtained for the new method were significantly higher than for control experiments lacking additives and were comparable to the yields obtained with the classical two-step approach. A study of the effect of beta-CD and CTAB concentrations on refolding yield suggested two operational regimes: slow stripping ( beta-CDXTABsimilar to1), most suited for higher protein concentrations, and fast stripping (beta-CD/CTABsimilar to2.7), best suited for lower protein concentrations. An increased chaotrope concentration resulted in higher refolding yields and an enlarged operational regime.

Interpretation of negative values of the interaction parameter in the adsorption equation through the effects of surface layer heterogeneity

The problem of the negative values of the interaction parameter in the equation of Frumkin has been analyzed with respect to the adsorption of nonionic molecules on energetically homogeneous surface. For this purpose, the adsorption states of a homologue series of ethoxylated nonionic surfactants on air/water interface have been determined using four different models and literature data (surface tension isotherms). The results obtained with the Frumkin adsorption isotherm imply repulsion between the adsorbed species (corresponding to negative values of the interaction parameter), while the classical lattice theory for energetically homogeneous surface (e.g., water/air) admits attraction alone. It appears that this serious contradiction can be overcome by assuming heterogeneity in the adsorption layer, that is, effects of partial condensation (formation of aggregates) on the surface. Such a phenomenon is suggested in the Fainerman-Lucassen-Reynders-Miller (FLM) 'Aggregation model'. Despite the limitations of the latter model (e.g., monodispersity of the aggregates), we have been able to estimate the sign and the order of magnitude of Frumkin's interaction parameter and the range of the aggregation numbers of the surface species. (C) 2004 Elsevier B.V All rights reserved.

Activation of NK cell cytotoxicity

Natural killer (NK) cells are innate effector lymphocytes necessary for defence against stressed, microbe-infected, or malignant cells. NK cells kill target cells by either of two major mechanisms that require direct contact between NK cells and target cells. In the first pathway, cytoplasmic granule toxins, predominantly a membrane-disrupting protein known as perforin, and a family of structurally related serine C, proteases (granzymes) with various substrate specificities, are secreted by exocytosis and together induce apoptosis of the target cell. The granule-exocytosis pathway potently activates cell-death mechanisms that operate through the activation of apoptotic cysteine proteases (caspases), but can also cause cell death in the absence of activated caspases. The second pathway involves the engagement of death receptors (e.g. Fas/CD95) on target cells by their cognate ligands (e.g. FasL on NK cells, resulting in classical caspase-dependent apoptosis. The comparative role of these pathways in the pathophysiology of many diseases is being dissected by analyses of gene-targeted mice that lack these molecules, and humans who have genetic mutations affecting these pathways. We are also now learning that the effector function of NK cells is controlled by interactions involving specific NK cell receptors and their cognate ligands, either on target cells, or other cells of the immune system. This review will discuss the functional importance of NK cell cytotoxicity and the receptor/ligand interactions that control these processes. (C) 2004 Elsevier Ltd. All rights reserved.