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.

Conotoxin TVIIA, a novel peptide from the venom of Conus tulipa - 2. Three-dimensional solution structure

The three-dimensional solution structure of conotoxin TVIIA, a 30-residue polypeptide from the venom of the piscivorous cone snail Conus tulipa, has been determined using 2D H-1 NMR spectroscopy. TVIIA contains six cysteine residues which form a 'four-loop' structural framework common to many peptides from Conus venoms including the omega-, delta-, kappa-, and mu O-conotoxins. However, TVIIA does not belong to these well-characterized pharmacological classes of conotoxins, but displays high sequence identity with conotoxin GS, a muscle sodium channel blocker from Conus geographus. Structure calculations were based on 562 interproton distance restraints inferred from NOE data, together with 18 backbone and nine side-chain torsion angle restraints derived from spin-spin coupling constants. The final family of 20 structures had mean pairwise rms differences over residues 2-27 of 0.18 +/- 0.05 Angstrom for the backbone atoms and 1.39 +/- 0.33 Angstrom for all heavy atoms. The structure consists of a triple-stranded, antiparallel beta sheet with +2x, -1 topology (residues 7-9, 16-20 and 23-27) and several beta turns. The core of the molecule is formed by three disulfide bonds which form a cystine knot motif common to many toxic and inhibitory polypeptides. The global fold, molecular shape and distribution of amino-acid sidechains in TVIIA is similar to that previously reported for conotoxin GS, and comparison with other four-loop conotoxin structures provides further indication that TVIIA and GS represent a new and distinct subgroup of this structural family. The structure of TVIIA determined in this study provides the basis for determining a structure-activity relationship for these molecules and their interaction with target receptors.

Cyclotron effective masses in layered metals

Many layered metals such as quasi-two-dimensional organic molecular crystals show properties consistent with a Fermi-liquid description at low temperatures. The effective masses extracted from the temperature dependence of the magnetic oscillations observed in these materials are in the range, m(c)*/m(e) similar to 1 - 7, suggesting that these systems are strongly correlated. However, the ratio m(c)*/m(e) contains both the renormalization due to the electron-electron interaction and the periodic potential of the lattice. We show that for any quasi-two-dimensional band structure, the cyclotron mass is proportional to the density-of-states at the Fermi energy. Due to Luttinger's theorem, this result is also valid in the presence of interactions. We then evaluate m(c) for several model band structures for the beta, kappa, and theta families of (BEDT-TTF)(2)X, where BEDT-TTF is bis-(ethylenedithia-tetrathiafulvalene) and X is an anion. We find that for kappa-(BEDT-TTF)(2)X, the cyclotron mass of the beta orbit, m(c)*(beta) is close to 2 m(c)*(alpha), where m(c)*(alpha) is the effective mass of the alpha orbit. This result is fairly insensitive to the band-structure details. For a wide range of materials we compare values of the cyclotron mass deduced from band-structure calculations to values deduced from measurements of magnetic oscillations and the specific-heat coefficient gamma.

Mesoionic pyridopyrimidinylium and pyridooxazinylium olates and non-mesoionic pyridopyrimidinones. Structures in the solid state, solution, and matrices

X-Ray crystal structures, C-13 NMR spectra and theoretical calculations (B3LYP/6-31G*) are reported for the mesoionic (zwitterionic) pyridopyrimidinylium- and pyridooxazinyliumolates 2a, 3a and 5a,b as well as the enol ether 11b and the enamine 11c. The 1-NH compounds like 1a, 2a and 3a exist in the mesoionic form in the crystal and in solution, but the OH tautomers such as 1b and 2b dominate in the gas phase as revealed by the Ar matrix IR spectra in conjunction with DFT calculations. All data indicate that the mesoionic compounds can be regarded as intramolecular pyridine-ketene zwitterions (cf. 16 --> 17) with a high degree of positive charge on the pyridinium nitrogen, a long pyridinium N-CO bond (ca. 1.44-1.49 Angstrom), and normal C=O double bonds (ca. 1.22 Angstrom). All mesoionic compounds exhibit a pronounced tilting of the olate C=O groups (the C=O groups formally derived from a ketene) towards the pyridinium nitrogen, giving NCO angles of 110-118 degrees. Calculations reveal a hydrogen bond with 6-CH, analogous to what is found in ketene-pyridine zwitterions and the C3O2-pyridine complex. The 2-OH tautomers of type 1b, 2b, and 11 also show a high degree of zwitterionic character as indicated by the canonical structures 11 12.

Three-dimensional structure of RK-1: A novel alpha-defensin peptide

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of RK-1, an antimicrobial peptide from rabbit kidney recently discovered from homology screening based on the distinctive physicochemical properties of the corticostatins/defensins. RK-1 consists of 32 residues, including six cysteines arranged into three disulfide bonds. It exhibits antimicrobial activity against Escherichia coli and activates Ca2+ channels in vitro. Through its physicochemical similarity, identical cysteine spacing, and linkage to the corticostatins/defensins, it was presumed to be a member of this family. However, RK-1 lacks both a large number of arginines in the primary sequence and a high overall positive charge, which are characteristic of this family of peptides. The three-dimensional solution structure, determined by NMR, consists of a triple-stranded antiparallel beta -sheet and a series of turns and is similar to the known structures of other alpha -defensins. This has enabled the definitive classification of RK-1 as a member of this family of antimicrobial peptides. Ultracentrifuge measurements confirmed that like rabbit neutrophil defensins, RK-1 is monomeric in solution, in contrast to human neutrophil defensins, which are dimeric.

The burden of mental disorders in Victoria, 1996

Background: Between 1998 and 1999, a burden of disease assessment was carried out in Victoria, Australia applying and improving on the methods of the Global Burden of Disease Study. This paper describes the methods and results of the calculations of the burden due to 22 mental disorders, adding 14 conditions not included in previous burden of disease estimates, Methods: The National Survey of Mental Health and Wellbeing provided recent data on the occurrence of the major adult mental disorders in Australia. Data from international studies and expert advice further contributed to the construction of disease models, describing each condition in terms of incidence, average duration and level of severity, with adjustments for comorbidity with other mental disorders. Disability weights for the time spent in different states of mental ill health were borrowed mainly from a study in the Netherlands, supplemented by weights derived in a local extrapolation exercise. Results: Mental disorders were the third largest group of conditions contributing to the burden of disease in Victoria, ranking behind cancers and cardiovascular diseases. Depression was the greatest cause of disability in both men and women. Eight other mental disorders in men and seven in women ranked among the top twenty causes of disability. Conclusions: Insufficient information on the natural history of many of the mental disorders, the limited information on the validity of mental disorder diagnoses in community surveys and considerable differences between ICD-10 and DSM-IV defined diagnoses were the main concerns about the accuracy of the estimates. Similar and often greater concerns have been raised in relation to the estimation of the burden from common non-fatal physical conditions such as asthma, diabetes and osteoarthritis. In comparison, psychiatric epidemiology can boast greater scientific rigour in setting standards for population surveys.

Determination of the stability constants of mercury(II) complexes of mixed donor macrobicyclic encapsulating ligands

The reactions of mercury(II) with the mixed donor encapsulating ligands 3,6,16-trithia-6,11,19-triazabicyclo[6.6.6]icosane (AMN(3)S(3)sar) and 1-amino-8-methyl-6,19-dithia-3,10,13,16-tetraazabicyclo[6.6.6]icosane (AMN(4)S(2)sar) have been studied. NMR ligand-ligand competition experiments with the ligands 1,4,8,11-tetraazaeyclotetradecane ([14]aneN(4)), 1-thia-4,7,10-triazacyclododecane ([12]aneN(3)S) and ethylenediaminetetraacetic acid (EDTA) with AMN(3)S(3)sar and Hg(II) indicated that [14]aneN(4) would be an appropriate competing ligand for the, determination of the Hg(II) stability constant. Calculations indicated the ratio of concentrations of AMN3S3sar, [14]aneN(4) and Hg(II) required for the determination of the stability constant ranged from 1:1:1 to 1:5:1. Refinement of the titration curves yielded log(10)K[Hg(AMN(3)S(3)sar)](2+) = 17.7. A similar competition titration resulted in the determination of the stability constant for the AMN(4)S(2)sar system as log(10)K[Hg(AMN(4)S(2)sar)](2+) = 19.5. The observed binding constants for the mixed N/S donor systems and the hexaaza analogues sar (3,6,10,13,16,19-hexaazabicyclo [6.6.6]icosane) and diamsar (1,8-diamino-3,6,10,13,16,19 -hexazabicyclo [6.6.6] icosane (log(10)K-[Hg(diamsar)](2+) = 26.4; log(10)K[Hg(sar)](2+) = 28.1) differ by approximately ten orders of magnitude. The difference is ascribed not to a cryptate effect but to a mismatch in the Hg-N and Hg-S bond lengths in the N/S systems.

Cost-effectiveness of cognitive behavioural therapy and selective serotonin reuptake inhibitors for major depression in children and adolescents

Objective: To assess from a health sector perspective the incremental cost-effectiveness of cognitive behavioural therapy (CBT) and selective serotonin reuptake inhibitors (SSRIs) for the treatment of major depressive disorder (MDD) in children and adolescents, compared to 'current practice'. Method: The health benefit is measured as a reduction in disability-adjusted life years (DALYs), based on effect size calculations from meta-analysis of randomised controlled trials. An assessment on second stage filter criteria ('equity'; 'strength of evidence', 'feasibility' and 'acceptability to stakeholders') is also undertaken to incorporate additional factors that impact on resource allocation decisions. Costs and benefits are tracked for the duration of a new episode of MDD arising in eligible children (age 6-17 years) in the Australian population in the year 2000. Simulation-modelling techniques are used to present a 95% uncertainty interval (UI) around the cost-effectiveness ratios. Results: Compared to current practice, CBT by public psychologists is the most cost-effective intervention for MDD in children and adolescents at A$9000 per DALY saved (95% UI A$3900 to A$24 000). SSRIs and CBT by other providers are less cost-effective but likely to be less than A$50 000 per DALY saved (> 80% chance). CBT is more effective than SSRIs in children and adolescents, resulting in a greater total health benefit (DALYs saved) than could be achieved with SSRIs. Issues that require attention for the CBT intervention include equity concerns, ensuring an adequate workforce, funding arrangements and acceptability to various stakeholders. Conclusions: Cognitive behavioural therapy provided by a public psychologist is the most effective and cost-effective option for the first-line treatment of MDD in children and adolescents. However, this option is not currently accessible by all patients and will require change in policy to allow more widespread uptake. It will also require 'start-up' costs and attention to ensuring an adequate workforce.

Cost-effectiveness of psychological and pharmacological interventions for generalized anxiety disorder and panic disorder

Objective: To assess from a health sector perspective the incremental cost-effectiveness of interventions for generalized anxiety disorder (cognitive behavioural therapy [CBT] and serotonin and noradrenaline reuptake inhibitors [SNRIs]) and panic disorder (CBT, selective serotonin reuptake inhibitors [SSRIs] and tricyclic antidepressants [TCAs]). Method: The health benefit is measured as a reduction in disability-adjusted life years (DALYs), based on effect size calculations from meta-analyses of randomised controlled trials. An assessment on second stage filters ('equity', 'strength of evidence', 'feasibility' and 'acceptability to stakeholders') is also undertaken to incorporate additional factors that impact on resource allocation decisions. Costs and benefits are calculated for a period of one year for the eligible population (prevalent cases of generalized anxiety disorder/panic disorder identified in the National Survey of Mental Health and Wellbeing, extrapolated to the Australian population in the year 2000 for those aged 18 years and older). Simulation modelling techniques are used to present 95% uncertainty intervals (UI) around the incremental cost-effectiveness ratios (ICERs). Results: Compared to current practice, CBT by a psychologist on a public salary is the most cost-effective intervention for both generalized anxiety disorder (A$6900/DALY saved; 95% UI A$4000 to A$12 000) and panic disorder (A$6800/DALY saved; 95% UI A$2900 to A$15 000). Cognitive behavioural therapy results in a greater total health benefit than the drug interventions for both anxiety disorders, although equity and feasibility concerns for CBT interventions are also greater. Conclusions: Cognitive behavioural therapy is the most effective and cost-effective intervention for generalized anxiety disorder and panic disorder. However, its implementation would require policy change to enable more widespread access to a sufficient number of trained therapists for the treatment of anxiety disorders.

Thermal analysis, nuclear magnetic resonance spectroscopy, and impedance spectroscopy of N,N-dimethyl-pyrrolidinium iodide: An ionic solid exhibiting rotator phases

N,N-dimethyl-pyrrolidinium iodide has been investigated using differential scanning calorimetry, nuclear magnetic resonance (NMR) spectroscopy, second moment calculations, and impedance spectroscopy. This pyrrolidinium salt exhibits two solid-solid phase transitions, one at 373 K having an entropy change, Delta S, of 38 J mol(-1) K-1 and one at 478 K having Delta S of 5.7 J mol(-1) K-1. The second moment calculations relate the lower temperature transition to a homogenization of the sample in terms of the mobility of the cations, while the high temperature phase transition is within the temperature region of isotropic tumbling of the cations. At higher temperatures a further decrease in the H-1 NMR linewidth is observed which is suggested to be due to diffusion of the cations. (C) 2005 American Institute of Physics.

Axial Young's modulus prediction of single-walled carbon nanotube arrays with diameters from nanometer to meter scales

Based on a self-similar array model, we systematically investigated the axial Young's modulus (Y-axis) of single-walled carbon nanotube (SWNT) arrays with diameters from nanometer to meter scales by an analytical approach. The results show that the Y-axis of SWNT arrays decreases dramatically with the increases of their hierarchy number (s) and is not sensitive to the specific size and constitution when s is the same, and the specific Young's modulus Y-axis(s) is independent of the packing configuration of SWNTs. Our calculations also show that the Y-axis of SWNT arrays with diameters of several micrometers is close to that of commercial high performance carbon fibers (CFs), but the Y-axis(s) of SWNT arrays is much better than that of high performance CFs. (C) 2005 American Institute of Physics.

Ferromagnetism, paramagnetism, and a Curie-Weiss metal in an electron-doped Hubbard model on a triangular lattice

Motivated by the unconventional properties and rich phase diagram of NaxCoO2 we consider the electronic and magnetic properties of a two-dimensional Hubbard model on an isotropic triangular lattice doped with electrons away from half-filling. Dynamical mean-field theory (DMFT) calculations predict that for negative intersite hopping amplitudes (t < 0) and an on-site Coulomb repulsion, U, comparable to the bandwidth, the system displays properties typical of a weakly correlated metal. In contrast, for t > 0 a large enhancement of the effective mass, itinerant ferromagnetism, and a metallic phase with a Curie-Weiss magnetic susceptibility are found in a broad electron doping range. The different behavior encountered is a consequence of the larger noninteracting density of states (DOS) at the Fermi level for t > 0 than for t < 0, which effectively enhances the mass and the scattering amplitude of the quasiparticles. The shape of the DOS is crucial for the occurrence of ferromagnetism as for t > 0 the energy cost of polarizing the system is much smaller than for t < 0. Our observation of Nagaoka ferromagnetism is consistent with the A-type antiferromagnetism (i.e., ferromagnetic layers stacked antiferromagnetically) observed in neutron scattering experiments on NaxCoO2. The transport and magnetic properties measured in NaxCoO2 are consistent with DMFT predictions of a metal close to the Mott insulator and we discuss the role of Na ordering in driving the system towards the Mott transition. We propose that the Curie-Weiss metal phase observed in NaxCoO2 is a consequence of the crossover from a bad metal with incoherent quasiparticles at temperatures T > T-* and Fermi liquid behavior with enhanced parameters below T-*, where T-* is a low energy coherence scale induced by strong local Coulomb electron correlations. Our analysis also shows that the one band Hubbard model on a triangular lattice is not enough to describe the unusual properties of NaxCoO2 and is used to identify the simplest relevant model that captures the essential physics in NaxCoO2. We propose a model which allows for the Na ordering phenomena observed in the system which, we propose, drives the system close to the Mott insulating phase even at large dopings.

A new algorithm and refined bounds for extended gcd computation

Extended gcd computation is interesting itself. It also plays a fundamental role in other calculations. We present a new algorithm for solving the extended gcd problem. This algorithm has a particularly simple description and is practical. It also provides refined bounds on the size of the multipliers obtained.

Paramagnetism of caesium titanium alum

Caesium titanium alum, CsTi(SO4)(2) . 12H(2)O, is a beta alum and exhibits a large trigonal field and a dynamic Jahn-Teller effect. Exact calculations of the linear (2)T(2)xe Jahn-Teller coupling show that in the strict S-6 Site symmetry the ground multiplet consists of a Kramers doublet 2 Gamma(6) with magnetic splitting factors g(parallel to)=1.1 and g perpendicular to=0, a Gamma(4) Gamma(5) doublet at similar to 60 cm(-1) with g(parallel to)=2.51 and g(perpendicular to)=0.06 and another Gamma(4) Gamma(5) doublet at similar to 270 cm(-1) with g(parallel to)=1.67 and g(perpendicular to)=1.83. The controversial g values observed below 4.2 K, g(parallel to)=1.25 and g(perpendicular to)=1.14, are shown to arise from low symmetry distortions. These distortions couple the vibronic levels and induce into the ground state the off-diagonal axial Zeeman interaction that exists between the first excited and the ground vibronic levels. (C) 1997 American Institute of Physics.

Population versus clinical view of case fatality from acute coronary heart disease - Results from the WHO MONICA Project 1985-1990

Background The clinical view of case fatality (CF) from acute myocardial infarction (AMI) in those reaching the hospital alive is different from the population view. Registration of both hospitalized AMI cases and out-of-hospital coronary heart disease (CHD) deaths in the WHO MONICA Project allows both views to be reconciled. The WHO MONICA Project provides the largest data set worldwide to explore the relationship between CHD CF and age, sex, coronary event rate, and first versus recurrent event. Methods and Results All 79 669 events of definite AMI or possible coronary death, occurring from 1985 to 90 among 5 725 762 people, 35 to 64 years of age, in 29 MONICA populations are the basis for CF calculations. Age-adjusted CF (percentage of CHD events that were fatal) was calculated across populations, stratified for different time periods, and related to age, sex, and CHD event rate. Median 28-day population CF was 49% (range, 35% to 60%) in men and 51% (range, 34% to 70%) in women and was particularly higher in women than men in populations in which CHD event rates were low. Median 28-day CF for hospitalized events was much lower: in men 22% (range, 15% to 36%) and in women 27% (range, 19% to 46%). Among hospitalized events CF was twice as high for recurrent as for first events. Conclusions Overall 28-day CF is halved for hospitalized events compared with all events and again nearly halved for hospitalized 24-hour survivors. Because approximately two thirds of 28-day CHD deaths in men and women occurred before reaching the hospital, opportunities for reducing CF through improved care in the acute event are limited. Major emphasis should be on primary and secondary prevention.