Scalable Visualization of Galaxies, Oceans, and Brains

One of the challenges in scientific visualization is to generate software libraries suitable for the large-scale data emerging from tera-scale simulations and instruments. We describe the efforts currently under way at SDSC and NPACI to address these challenges. The scope of the SDSC project spans data handling, graphics, visualization, and scientific application domains. Components of the research focus on the following areas: intelligent data storage, layout and handling, using an associated “Floor-Plan” (meta data); performance optimization on parallel architectures; extension of SDSC’s scalable, parallel, direct volume renderer to allow perspective viewing; and interactive rendering of fractional images (“imagelets”), which facilitates the examination of large datasets. These concepts are coordinated within a data-visualization pipeline, which operates on component data blocks sized to fit within the available computing resources. A key feature of the scheme is that the meta data, which tag the data blocks, can be propagated and applied consistently. This is possible at the disk level, in distributing the computations across parallel processors; in “imagelet” composition; and in feature tagging. The work reflects the emerging challenges and opportunities presented by the ongoing progress in high-performance computing (HPC) and the deployment of the data, computational, and visualization Grids.

Twisted parafermions

A new type of nonlocal currents (quasi-particles), which we call twisted parafermions, and its corresponding twisted Z-algebra are found. The system consists of one spin-1 bosonic field and six nonlocal fields of fractional spins. Jacobi-type identities for the twisted parafermions are derived, and a new conformal field theory is constructed from these currents. As an application, a parafermionic representation of the twisted affine current algebra A(2)((2)) is given.

The relationship between plasma lactate parameters, W-peak and 1-h cycling performance in women

Purpose: The relationship between six descriptors of lactate increase, peak (V) over dot O-2,W-peak, and 1-h cycling performance were compared in 24 trained, female cyclists (peak (V) over dot O-2 = 48.11 +/- 6.32 mL . kg(-1) . min(-1)). Methods: The six descriptors of lactate increase were: 1) lactate threshold (LT; the power output at which plasma lactate concentration begins to increase above the resting level during an incremental exercise test), 2) LT1 (the power output at which plasma lactate increases by 1 mM or more), 3) LTD (the lactate threshold calculated by the D-max method), 4) LTMOD (the lactate threshold calculated by a modified D-max method), 5) L4 (the power output at which plasma lactate reaches a concentration of 4 mmol-L-1), and 6) LTLOG (the power output at which plasma lactate concentration begins to increase when the log([La-]) is plotted against the log (power output)). Subjects first completed a peak (V) over dot O-2 test on a cycle ergometer. Finger-tip capillary blood was sampled within 30 s of the end of each 3-min stage for analysis of plasma lactate. Endurance performance was assessed 7 d later using a 1-h cycle test (OHT) in which subjects were directed to achieve the highest possible average power output. Results: The mean power output (W) for the OHT (+/- SD) was 183.01 +/- 18.88, and for each lactate variable was: LT (138.54 +/- 46.61), LT1 (179.17 +/- 27.25), LTLOG (143.97 +/- 45.74), L4 (198.09 +/- 33.84), LTD (178.79 +/- 24.07), LTMOD (212.28 +/- 31.75). Average power output during the OHT was more strongly correlated with all plasma lactate parameters (0.61 < r < 0.84) and W-peak (r = 0.81) than with peak (V) over dot O-2 (r = 0.55). The six lactate parameters were strongly correlated with each other (0.54 < r < 0.91) and of the six lactate parameters, LTD correlated best with endurance performance (r = 0.84). Conclusions: It was concluded that plasma lactate parameters and W-peak provide better indices of endurance performance than peak (V) over dot O-2 and that, of the six descriptors of lactate increase measured in this study, LTD is most strongly related to 1-h cycling performance in trained, female cyclists.

An evaluation of the time dependence of the anisotropy of the water diffusion tensor in acute human ischemia

We have performed MRI examinations to determine the water diffusion tensor in the brain of six patients who were admitted to the hospital within 12 h after the onset of cerebral ischemic symptoms. The examinations have been carried out immediately after admission, and thereafter at varying intervals up to 90 days post admission. Maps of the trace of the diffusion tensor, the fractional anisotropy and the lattice index, as well as maps of cerebral blood perfusion parameters, were generated to quantitatively assess the character of the water diffusion tensor in the infarcted area. In patients with significant perfusion deficits and substantial lesion volume changes, four of six cases, our measurements show a monotonic and significant decrease in the diffusion anisotropy within the ischemic lesion as a function of time. We propose that retrospective analysis of this quantity, in combination with brain tissue segmentation and cerebral perfusion maps, may be used in future studies to assess the severity of the ischemic event. (C) 1999 Elsevier Science Inc.

Thermal and radiation curing of phenylethynyl terminated macromers

The thermal and gamma-irradiation induced curing of two phenylethynyl terminated composite resin systems, DFB/BPF and PETI5A, was investigated. Thermal curing of these matrix resin samples was performed at a temperature of 360 degrees C, gamma irradiation of the samples was conducted at 300 degrees C at a dose rate of 2.2 kGy h(-1). The reaction and subsequent loss of ethynyl groups in the resins for both cure methods was demonstrated by observing the decrease of the 2215 cm(-1) peak in the Raman spectra of the resins. Fully cured resin samples were found to have glass transition temperatures of 244-246 degrees C and 278-280 degrees C for DFB/BPF and PETI5A respectively. Similar relationships between T-g and fractional conversion were observed in both resins. The apparent polymerization rate, R-p, for thermal cure at 360 degrees C, was found to be 4.79 x 10(-2)% s(-1) in PETI5A and 3.22 x 10(-2)% s(-1) in DFB/BPF. Catastrophic degradation under nitrogen was observed to commence near 450 degrees C and 530 degrees C, with 5% weight losses occurring at 455 degrees C and 540 degrees C for DFB/BPF and PETI5A respectively. Gamma radiation induced cure at 300 degrees C was shown to be feasible, with full cure being reached with doses of 40 kGy for DFB/BPF and 100 kGy for PETI5A.

Determining forest structural attributes using an inverted geometric-optical model in mixed eucalypt forests, Southeast Queensland, Australia

The Montreal Process indicators are intended to provide a common framework for assessing and reviewing progress toward sustainable forest management. The potential of a combined geometrical-optical/spectral mixture analysis model was assessed for mapping the Montreal Process age class and successional age indicators at a regional scale using Landsat Thematic data. The project location is an area of eucalyptus forest in Emu Creek State Forest, Southeast Queensland, Australia. A quantitative model relating the spectral reflectance of a forest to the illumination geometry, slope, and aspect of the terrain surface and the size, shape, and density, and canopy size. Inversion of this model necessitated the use of spectral mixture analysis to recover subpixel information on the fractional extent of ground scene elements (such as sunlit canopy, shaded canopy, sunlit background, and shaded background). Results obtained fron a sensitivity analysis allowed improved allocation of resources to maximize the predictive accuracy of the model. It was found that modeled estimates of crown cover projection, canopy size, and tree densities had significant agreement with field and air photo-interpreted estimates. However, the accuracy of the successional stage classification was limited. The results obtained highlight the potential for future integration of high and moderate spatial resolution-imaging sensors for monitoring forest structure and condition. (C) Elsevier Science Inc., 2000.

Tubular system volume changes in twitch fibres from toad and rat skeletal muscle assessed by confocal microscopy

The volume of the extracellular compartment (tubular system) within intact muscle fibres from cane toad and rat was measured under various conditions using confocal microscopy. Under physiological conditions at rest, the fractional volume of the tubular system (t-sys(Vol)) was 1.38 +/- 0.09% (n = 17),1.41 +/- 0.09% (n = 12) and 0.83 +/- 0.07% (n = 12) of the total fibre volume in the twitch fibres from toad iliofibularis muscle, rat extensor digitorum longus muscle and rat soleus muscle, respectively. In toad muscle fibres, the t-sys(Vol) decreased by 30% when the tubular system was fully depolarized and decreased by 15% when membrane cholesterol was depleted from the tubular system with methyl-beta-cyclodextrin but did not change as the sarcomere length was changed from 1.93 to 3.30 mum. There was also an increase by 30% and a decrease by 25% in t-sys(Vol) when toad fibres were equilibrated in solutions that were 2.5-fold hypertonic and 50% hypotonic, respectively. When the changes in total fibre volume were taken into consideration, the t-sys(Vol) expressed as a percentage of the isotonic fibre volume did actually decrease as tonicity increased, revealing that the tubular system in intact fibres cannot be compressed below 0.9% of the isotonic fibre volume. The results can be explained in terms of forces acting at the level of the tubular wall. These observations have important physiological implications showing that the tubular system is a dynamic membrane structure capable of changing its volume in response to the membrane potential, cholesterol depletion and osmotic stress but not when the sarcomere length is changed in resting muscle.

The kinetics of baculovirus adsorption to insect cells in suspension culture

The influence of various culture parameters on the attachment of a recombinant baculovirus to suspended insect cells was examined under normal culture conditions. These parameters included cell density, multiplicity of infection, and composition of the cell growth medium. It was found that the fractional rate of virus attachment was independent of the multiplicity of infection but dependent on the cell density. A first order mathematical model was used to simulate the adsorption kinetics and predict the efficiency of virus attachment under the various culture conditions. This calculated efficiency of virus attachment was observed to decrease at high cell densities, which was attributed to cell clumping. It was also observed that virus attachment was more efficient in Sf900II serum free medium than it was in IPL-41 serum-supplemented medium. This effect was attributed to the protein in serum which may coat the cells and so inhibit adsorption. A general discussion relating the observations made in-these experiments to the kinetics of recombinant baculovirus adsorption to suspended insect cells is presented.

Recovery of poly-3-hydroxybutyrate from recombinant Escherichia coli by homogenization and centrifugation

Poly-3-hydroxybutyrate from recombinant E. coli was recovered using homogenization and continuous centrifugation with a purity of 94%. Final protein and DNA concentrations were 1.0% w/w and 1.9% w/w, respectively, when a hypochlorite treatment was employed prior to centrifugation. High fractional cell debris removal (94%) was achieved with two centrifugation steps.

Reversal of cardiovascular remodelling with candesart

Cardiovascular remodelling, defined as ventricular and vascular hypertrophy together with fibrosis, characterises hypertension following inhibition of the production of the endogenous vasodilator, nitric oxide (NO). This study has determined whether the cardiovascular remodelling following chronic NO synthase inhibition can e reversed by administration of the selective angiotensin II AT(1)-receptor antagonist, candesartan. Male Wistar rats were treated with L-nitroarginine methyl ester (L-NAME, 400 mg/l in drinking water) for eight weeks and with candesartan cilexetil (2 mg/kg/day by oral gavage) for the last four weeks. L-NAME-treated rats became hypertensive with systolic blood pressure increasing from 110 +/- 4 mmHg (control) to 170 +/- 10 mmHg. Rats developed left ventricular hypertrophy (control 1.70 +/- 0.06; L-NAME 2.10 +/- 0.04 mg/kg body wt) with markedly increased deposition of perivascular and interstitial collagen. Candesartan returned blood pressure, left ventricular weights and collagen deposition to control values. Echo cardiographic assessment showed concentric hypertrophy with an increased fractional shortening; this was reversed by candesartan treatment. Heart failure was not evident. In the isolated Langendorff heart, diastolic stiffness increased in L-NAME-treated rats while the rate of increase in pressure (+dP/dt) increased after eight weeks only; candesartan reduced collagen deposition and normalised +dP/dt. In isolated left ventricular papillary muscles, the potency (negative log EC50) of noradrenaline as a positive inotropic compound was unchanged, (control 6.56 +/- 0.14); maximal increase in force before ectopic beats was reduced from 5.0 +/- 0.4 mN to 2.0 +/- 0.2 mN. Noradrenaline potency as a vasoconstrictor in thoracic aortic rings was unchanged, but maximal contraction was markedly reduced from 25.2 +/- 2.0 mN to 3.0 +/- 0.3 mN; this was partially reversed by candesartan treatment. Thus, chronic inhibition of NO production with L-NAME induces hypertension, hypertrophy and fibrosis with increased toxicity and significant decreases in vascular responses to noradrenaline. These changes were at least partially reversible by treatment with candesartan, implying a significant role of AT(1)-receptors in L-NAME-induced cardiovascular changes.

Set-valued Markov chains and negative semitrajectories of discretized dynamical systems

Computer simulation of dynamical systems involves a phase space which is the finite set of machine arithmetic. Rounding state values of the continuous system to this grid yields a spatially discrete dynamical system, often with different dynamical behaviour. Discretization of an invertible smooth system gives a system with set-valued negative semitrajectories. As the grid is refined, asymptotic behaviour of the semitrajectories follows probabilistic laws which correspond to a set-valued Markov chain, whose transition probabilities can be explicitly calculated. The results are illustrated for two-dimensional dynamical systems obtained by discretization of fractional linear transformations of the unit disc in the complex plane.

Petrology and geochemistry of the Bandas del Sur Formation, Las Canadas Edifice, Tenerife (Canary Islands)

The Bandas del Sur Formation preserves a Quaternary extra-caldera record of central phonolitic explosive volcanism of the Las Canadas volcano at Tenerife. Volcanic rocks are bimodal in composition, being predominantly phonolitic pyroclastic deposits, several eruptions of which resulted in summit caldera collapse, alkali basaltic lavas erupted from many fissures around the flanks. For the pyroclastic deposits, there is a broad range of pumice glass compositions from phonotephrite to phonolite. The phonolite pyroclastic deposits are also characterized by a diverse, 7-8-phase phenocryst assemblage (alkali feldspar + biotite + sodian diopside + titanomagnetite + ilmenite + nosean-hauyne + titanite + apatite) with alkali feldspar dominant, in contrast to interbedded phonolite lavas that typically have lower phenocryst contents and lack hydrous phases. Petrological and geochemical data are consistent with fractional crystallization (involving the observed phenocryst assemblages) as the dominant process in the development of phonolite magmas. New stratigraphically constrained data indicate that petrological and geochemical differences exist between pyroclastic deposits of the last two explosive cycles of phonolitic volcanism. Cycle 2 (0.85-0.57 Ma) pyroclastic fall deposits commonly show a cryptic compositional zonation indicating that several eruptions tapped chemically, and probably thermally stratified magma systems. Evidence for magma mixing is most widespread in the pyroclastic deposits of Cycle 3 (0.37-0.17 Ma), which includes the presence of reversely and normally zoned phenocrysts, quenched mafic glass blebs in pumice, banded pumice, and bimodal to polymodal phenocryst compositional populations. Syn-eruptive mixing events involved mostly phonolite and tephriphonolite magmas, whereas a pre-eruptive mixing event involving basaltic magma is recorded in several banded pumice-bearing ignimbrites of Cycle 3. The periodic addition and mixing of basaltic magma ultimately may have triggered several eruptions. Recharge and underplating by basaltic magma is interpreted to have elevated sulphur contents (occurring as an exsolved gas phase) in the capping phonolitic magma reservoir. This promoted nosean-hauyne crystallization over nepheline, elevated SO3 contents in apatite, and possibly resulted in large, climatologically important SO2 emissions.

Fluid-mobile trace element constraints on the role of slab melting and implications for Archaean crustal growth models

We use published and new trace element data to identify element ratios which discriminate between arc magmas from the supra-subduction zone mantle wedge and those formed by direct melting of subducted crust (i.e. adakites). The clearest distinction is obtained with those element ratios which are strongly fractionated during refertilisation of the depleted mantle wedge, ultimately reflecting slab dehydration. Hence, adakites have significantly lower Pb/Nd and B/Be but higher Nb/Ta than typical arc magmas and continental crust as a whole. Although Li and Be are also overenriched in continental crust, behaviour of Li/Yb and Be/Nd is more complex and these ratios do not provide unique signatures of slab melting. Archaean tonalite-trondhjemite-granodiorites (TTGs) strongly resemble ordinary mantle wedge-derived arc magmas in terms of fluid-mobile trace element content, implying that they-did not form by slab melting but that they originated from mantle which was hydrated and enriched in elements lost from slabs during prograde dehydration. We suggest that Archaean TTGs formed by extensive fractional crystallisation from a mafic precursor. It is widely claimed that the time between the creation and subduction of oceanic lithosphere was significantly shorter in the Archaean (i.e. 20 Ma) than it is today. This difference was seen as an attractive explanation for the presumed preponderance of adakitic magmas during the first half of Earth's history. However, when we consider the effects of a higher potential mantle temperature on the thickness of oceanic crust, it follows that the mean age of oceanic lithosphere has remained virtually constant. Formation of adakites has therefore always depended on local plate geometry and not on potential mantle temperature.

MR image-based measurement of rates of change in volumes of brain structures. Part 1: Method and validation

A detailed analysis procedure is described for evaluating rates of volumetric change in brain structures based on structural magnetic resonance (MR) images. In this procedure, a series of image processing tools have been employed to address the problems encountered in measuring rates of change based on structural MR images. These tools include an algorithm for intensity non-uniforniity correction, a robust algorithm for three-dimensional image registration with sub-voxel precision and an algorithm for brain tissue segmentation. However, a unique feature in the procedure is the use of a fractional volume model that has been developed to provide a quantitative measure for the partial volume effect. With this model, the fractional constituent tissue volumes are evaluated for voxels at the tissue boundary that manifest partial volume effect, thus allowing tissue boundaries be defined at a sub-voxel level and in an automated fashion. Validation studies are presented on key algorithms including segmentation and registration. An overall assessment of the method is provided through the evaluation of the rates of brain atrophy in a group of normal elderly subjects for which the rate of brain atrophy due to normal aging is predictably small. An application of the method is given in Part 11 where the rates of brain atrophy in various brain regions are studied in relation to normal aging and Alzheimer's disease. (C) 2002 Elsevier Science Inc. All rights reserved.