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.

Parallel 3D Simulation of a Fault Gouge using the Lattice Solid Model

Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.

Timing and duration of the last interglacial inferred from high resolution U-series chronology of stalagmite growth in Southern Hemisphere

High-precision Th-230-U-238 ages for a stalagmite from Newdegate Cave in southern Tasmania, Australia define a rare record of precipitation between 100 and 155 ka before the present. The fastest stalagmite growth occurred between 129.2 +/- 1.6 and 122.1 +/- 2.0 ka (similar to 61.5 mm/ka), coinciding with a time of prolific coral growth from Western Australia (128-122 ka). This is the first high-resolution continental record in the Southern Hemisphere that can be compared and correlated with the marine record. Such correlation shows that in southern Australia the onset of full interglacial sea level and the initiation of highest precipitation on land were synchronous. The stalagmite growth rate between 129.2 and 142.2 ka (similar to 5.9 mm/ka) was lower than that between 142.2 and 154.5 ka (similar to 18.7 mm/ka), implying drier conditions during the Penultimate Deglaciation, despite rising temperature and sea level. This asymmetrical precipitation pattern is caused by latitudinal movement of subtropical highs and an associated Westerly circulation, in response to a changing Equator-to-Pole temperature gradient. Both marine and continental records in Australia strongly suggest that the insolation maximum between 126 and 128 ka at 65 degreesN was directly responsible for the maintenance of full Last Interglacial conditions, although the triggers that initiated Penultimate Deglaciation (at similar to 142 ka) remain unsolved. (C) 2001 Elsevier Science B.V. All rights reserved.

Hepatic disposition and metabolite kinetics of a homologous series of diflunisal esters

The hepatic disposition and metabolite kinetics of a homologous series of diflunisal O-acyl esters (acetyl, butanoyl, pentanoyl, anti hexanoyl) were determined using a single-pass perfused in situ rat liver preparation. The experiments were conducted using 2% BSA Krebs-Henseleit buffer (pH 7.4), and perfusions were performed at 30 mL/min in each liver. O-Acyl esters of diflunisal and pregenerated diflunisal were injected separately into the portal vein. The venous outflow samples containing the esters and metabolite diflunisal were analyzed by high performance liquid chromatography (HPLC). The normalized outflow concentration-time profiles for each parent ester and the formed metabolite, diflunisal, were analyzed using statistical moments analysis and the two-compartment dispersion model. Data (presented as mean +/- standard error for triplicate experiments) was compared using ANOVA repeated measures, significance level P < 0.05. The hepatic availability (AUC'), the fraction of the injected dose recovered in the outflowing perfusate, for O-acetyldiflunisal (C2D = 0.21 +/- 0.03) was significantly lower than the other esters (0.34-0.38). However, R-N/f(u), the removal efficiency number R-N divided by the unbound fraction in perfusate f(u), which represents the removal efficiency of unbound ester by the liver, was significantly higher for the most lipophilic ester (O-hexanoyldiflunisal, C6D = 16.50 +/- 0.22) compared to the other members of the series (9.57 to 11.17). The most lipophilic ester, C6D, had the largest permeability surface area (PS) product (94.52 +/- 38.20 mt min-l g-l liver) and tissue distribution value VT (35.62 +/- 11.33 mL g(-1) liver) in this series. The MTT of these O-acyl esters of diflunisal were not significantly different from one another. However, the metabolite diflunisal MTTs tended to increase with the increase in the parent ester lipophilicity (11.41 +/- 2.19 s for C2D to 38.63 +/- 9.81 s for C6D). The two-compartment dispersion model equations adequately described the outflow profiles for the parent esters and the metabolite diflunisal formed from the O-acyl esters of diflunisal in the liver.

Novel composite membranes for gas separation: Preparation and performance

High performance composite membranes based on molecular sieving silica (MSS) were synthesized using sols containing silicon co-polymers (methyltriethoxysilane and tetraethylorthosilicate). Alpha alumina supports were treated with hydrochloric acid prior to sol deposition. Permselectivity of CO2 over CH4 as high as 16.68 was achieved whilst permeability of CO2 up to 36.7 GPU (10(-6) cm(3) (STP) cm(-2) . s(-1) . cm Hg-1) was measured. The best membrane's permeability was finger printed during various stages of the synthesis process showing an increase in CO2/CH4 permselectivity by over 25 times from initial support condition (no membrane film) to the completion of pore structure tailoring. Transport measurement results indicate that the membrane pretreated with HCl has highest permselectivity and permeation rate. In particular, there is a definite cut-off pore size between 3.3 and 3.4 angstroms which is just below the kinetic diameters of Ar and CH4. This demonstrates that the mechanism for the separation in the prepared composite membrane is molecular sieving (activated diffusion), rather than Knudsen diffusion.

High level of aspartic acid-bond isomerization during the synthesis of an N-linked tau glycopeptide

An increased degree of utilization of the potential N-glycosylation site In the fourth repeat unit of the human tau protein may be involved in the inability of tau to bind to the corresponding tubulin sequence(s) and in the subsequent development of the paired helical filaments of Alzheimer's disease. To model these processes, we synthesized the octadecapeptide spanning this region without sugar, and with the addition of an N-acetyl-glucosamine moiety. The carbohydrate-protected, glycosylated asparagine was incorporated as a building block during conventional Fmoc-solid phase peptide synthesis. While the crude non-glycosylated analog was obtained as a single peptide, two peptides with, the identical, expected masses, in approximately equal amounts, were detected after the cleavage of the peracetylated glycopeptide. Surprisingly, the two glycopeptides switched positions on the reversed-phase high performance liquid chromatogram after removal of the sugar-protecting acetyl groups. Nuclear magnetic resonance spectroscopy and peptide sequencing identified the more hydrophobic deprotected peak as the target peptide, and the more hydrophilic deprotected peak as a peptide analog in which the aspartic acid-bond just preceding the glycosylated asparagine residue was isomerized resulting in the formation of a beta-peptide. The anomalous chromatographic behavior of the acetylated beta-isomer could be explained on the basis of the generation of an extended hydrophobic surface which is not present in any of the other three glycopeptide variants. Repetition of the syntheses, with altered conditions and reagents, revealed reproducibly high levels of aspartic acid-bond isomerization of the glycopeptide as well as lack of isomerization for the non-glycosylated parent analog. If similar increased aspartic acid-bond isomerization occurs in vivo, a protein modification well known to take place for both the amyloid deposits and the neurofibrillary tangles in Alzheimer's disease, this process may explain the aggregation of glycosylated tau into the paired helical filaments in the affected brains. Copyright (C) 1999 European Peptide Society and John Wiley & Sons, Ltd.