Exact quantum Monte Carlo methods for correlated electron systems

Thema dieser Arbeit ist die Entwicklung und Kombination verschiedener numerischer Methoden, sowie deren Anwendung auf Probleme stark korrelierter Elektronensysteme. Solche Materialien zeigen viele interessante physikalische Eigenschaften, wie z.B. Supraleitung und magnetische Ordnung und spielen eine bedeutende Rolle in technischen Anwendungen. Es werden zwei verschiedene Modelle behandelt: das Hubbard-Modell und das Kondo-Gitter-Modell (KLM). In den letzten Jahrzehnten konnten bereits viele Erkenntnisse durch die numerische Lösung dieser Modelle gewonnen werden. Dennoch bleibt der physikalische Ursprung vieler Effekte verborgen. Grund dafür ist die Beschränkung aktueller Methoden auf bestimmte Parameterbereiche. Eine der stärksten Einschränkungen ist das Fehlen effizienter Algorithmen für tiefe Temperaturen.rnrnBasierend auf dem Blankenbecler-Scalapino-Sugar Quanten-Monte-Carlo (BSS-QMC) Algorithmus präsentieren wir eine numerisch exakte Methode, die das Hubbard-Modell und das KLM effizient bei sehr tiefen Temperaturen löst. Diese Methode wird auf den Mott-Übergang im zweidimensionalen Hubbard-Modell angewendet. Im Gegensatz zu früheren Studien können wir einen Mott-Übergang bei endlichen Temperaturen und endlichen Wechselwirkungen klar ausschließen.rnrnAuf der Basis dieses exakten BSS-QMC Algorithmus, haben wir einen Störstellenlöser für die dynamische Molekularfeld Theorie (DMFT) sowie ihre Cluster Erweiterungen (CDMFT) entwickelt. Die DMFT ist die vorherrschende Theorie stark korrelierter Systeme, bei denen übliche Bandstrukturrechnungen versagen. Eine Hauptlimitation ist dabei die Verfügbarkeit effizienter Störstellenlöser für das intrinsische Quantenproblem. Der in dieser Arbeit entwickelte Algorithmus hat das gleiche überlegene Skalierungsverhalten mit der inversen Temperatur wie BSS-QMC. Wir untersuchen den Mott-Übergang im Rahmen der DMFT und analysieren den Einfluss von systematischen Fehlern auf diesen Übergang.rnrnEin weiteres prominentes Thema ist die Vernachlässigung von nicht-lokalen Wechselwirkungen in der DMFT. Hierzu kombinieren wir direkte BSS-QMC Gitterrechnungen mit CDMFT für das halb gefüllte zweidimensionale anisotrope Hubbard Modell, das dotierte Hubbard Modell und das KLM. Die Ergebnisse für die verschiedenen Modelle unterscheiden sich stark: während nicht-lokale Korrelationen eine wichtige Rolle im zweidimensionalen (anisotropen) Modell spielen, ist in der paramagnetischen Phase die Impulsabhängigkeit der Selbstenergie für stark dotierte Systeme und für das KLM deutlich schwächer. Eine bemerkenswerte Erkenntnis ist, dass die Selbstenergie sich durch die nicht-wechselwirkende Dispersion parametrisieren lässt. Die spezielle Struktur der Selbstenergie im Impulsraum kann sehr nützlich für die Klassifizierung von elektronischen Korrelationseffekten sein und öffnet den Weg für die Entwicklung neuer Schemata über die Grenzen der DMFT hinaus.

Lamellar body ultrastructure revisited: high-pressure freezing and cryo-electron microscopy of vitreous sections

Lamellar bodies are the storage sites for lung surfactant within type II alveolar epithelial cells. The structure-function models of lamellar bodies are based on microscopic analyses of chemically fixed tissue. Despite available alternative fixation methods that are less prone to artifacts, such as cryofixation by high-pressure freezing, the nature of the lung, being mostly air filled, makes it difficult to take advantage of these improved methods. In this paper, we propose a new approach and show for the first time the ultrastructure of intracellular lamellar bodies based on cryo-electron microscopy of vitreous sections in the range of nanometer resolution. Thus, unspoiled by chemical fixation, dehydration and contrasting agents, a close to native structure is revealed. Our approach uses perfluorocarbon to substitute the air in the alveoli. Lung tissue was subsequently high-pressure frozen, cryosectioned and observed in a cryo-electron microscope. The lamellar bodies clearly show a tight lamellar morphology. The periodicity of these lamellae was 7.3 nm. Lamellar bifurcations were observed in our cryosections. The technical approach described in this paper allows the examination of the native cellular ultrastructure of the surfactant system under near in vivo conditions, and therefore opens up prospectives for scrutinizing various theories of lamellar body biogenesis, exocytosis and recycling.

Competition between single and double electron transfer in collisions of doubly charged molecular pyrrole ions with neutral pyrrole molecules

Cross-sections have been determined for one- and two-electron transfer channels in the collisions of keV gas-phase doubly charged pyrrole ions with pyrrole molecules. Measured single and double electron transfer total cross-sections approximate 45 Å2 and 15 Å2, respectively. A combination of symmetric resonance charge exchange and multistate curve-crossing models has been invoked to describe these reactions.

Search for a heavy particle decaying into an electron and a muon with the ATLAS detector in sqrt[s] = 7 TeV pp collisions at the LHC

This Letter presents the first search for a heavy particle decaying into an e ± μ(-/+) final state in sqrt[s] = 7 TeV pp collisions at the LHC. The data were recorded by the ATLAS detector during 2010 and correspond to a total integrated luminosity of 35 pb(-1). No excess above the standard model background expectation is observed. Exclusions at 95% confidence level are placed on two representative models. In an R-parity violating supersymmetric model, tau sneutrinos with a mass below 0.75 TeV are excluded, assuming all R-parity violating couplings are zero except λ(311)' = 0.11 and λ312 = 0.07. In a lepton flavor violating model, a Z'-like vector boson with masses of 0.70-1.00 TeV and corresponding cross sections times branching ratios of 0.175-0.183 pb is excluded. These results extend to higher mass R-parity violating sneutrinos and lepton flavor violating Z's than previous constraints from the Tevatron.

Photon-beam subsource sensitivity to the initial electron-beam parameters

One limitation to the widespread implementation of Monte Carlo (MC) patient dose-calculation algorithms for radiotherapy is the lack of a general and accurate source model of the accelerator radiation source. Our aim in this work is to investigate the sensitivity of the photon-beam subsource distributions in a MC source model (with target, primary collimator, and flattening filter photon subsources and an electron subsource) for 6- and 18-MV photon beams when the energy and radial distributions of initial electrons striking a linac target change. For this purpose, phase-space data (PSD) was calculated for various mean electron energies striking the target, various normally distributed electron energy spread, and various normally distributed electron radial intensity distributions. All PSD was analyzed in terms of energy, fluence, and energy fluence distributions, which were compared between the different parameter sets. The energy spread was found to have a negligible influence on the subsource distributions. The mean energy and radial intensity significantly changed the target subsource distribution shapes and intensities. For the primary collimator and flattening filter subsources, the distribution shapes of the fluence and energy fluence changed little for different mean electron energies striking the target, however, their relative intensity compared with the target subsource change, which can be accounted for by a scaling factor. This study indicates that adjustments to MC source models can likely be limited to adjusting the target subsource in conjunction with scaling the relative intensity and energy spectrum of the primary collimator, flattening filter, and electron subsources when the energy and radial distributions of the initial electron-beam change.

Electron microscopy of high pressure frozen samples: bridging the gap between cellular ultrastructure and atomic resolution

Transmission electron microscopy has provided most of what is known about the ultrastructural organization of tissues, cells, and organelles. Due to tremendous advances in crystallography and magnetic resonance imaging, almost any protein can now be modeled at atomic resolution. To fully understand the workings of biological "nanomachines" it is necessary to obtain images of intact macromolecular assemblies in situ. Although the resolution power of electron microscopes is on the atomic scale, in biological samples artifacts introduced by aldehyde fixation, dehydration and staining, but also section thickness reduces it to some nanometers. Cryofixation by high pressure freezing circumvents many of the artifacts since it allows vitrifying biological samples of about 200 mum in thickness and immobilizes complex macromolecular assemblies in their native state in situ. To exploit the perfect structural preservation of frozen hydrated sections, sophisticated instruments are needed, e.g., high voltage electron microscopes equipped with precise goniometers that work at low temperature and digital cameras of high sensitivity and pixel number. With them, it is possible to generate high resolution tomograms, i.e., 3D views of subcellular structures. This review describes theory and applications of the high pressure cryofixation methodology and compares its results with those of conventional procedures. Moreover, recent findings will be discussed showing that molecular models of proteins can be fitted into depicted organellar ultrastructure of images of frozen hydrated sections. High pressure freezing of tissue is the base which may lead to precise models of macromolecular assemblies in situ, and thus to a better understanding of the function of complex cellular structures.

Integration of computational models and experimental characterization to study internal frost damage in cementitious materials

The objective of this doctoral research is to investigate the internal frost damage due to crystallization pore pressure in porous cement-based materials by developing computational and experimental characterization tools. As an essential component of the U.S. infrastructure system, the durability of concrete has significant impact on maintenance costs. In cold climates, freeze-thaw damage is a major issue affecting the durability of concrete. The deleterious effects of the freeze-thaw cycle depend on the microscale characteristics of concrete such as the pore sizes and the pore distribution, as well as the environmental conditions. Recent theories attribute internal frost damage of concrete is caused by crystallization pore pressure in the cold environment. The pore structures have significant impact on freeze-thaw durability of cement/concrete samples. The scanning electron microscope (SEM) and transmission X-ray microscopy (TXM) techniques were applied to characterize freeze-thaw damage within pore structure. In the microscale pore system, the crystallization pressures at sub-cooling temperatures were calculated using interface energy balance with thermodynamic analysis. The multi-phase Extended Finite Element Modeling (XFEM) and bilinear Cohesive Zone Modeling (CZM) were developed to simulate the internal frost damage of heterogeneous cement-based material samples. The fracture simulation with these two techniques were validated by comparing the predicted fracture behavior with the captured damage from compact tension (CT) and single-edge notched beam (SEB) bending tests. The study applied the developed computational tools to simulate the internal frost damage caused by ice crystallization with the two dimensional (2-D) SEM and three dimensional (3-D) reconstructed SEM and TXM digital samples. The pore pressure calculated from thermodynamic analysis was input for model simulation. The 2-D and 3-D bilinear CZM predicted the crack initiation and propagation within cement paste microstructure. The favorably predicted crack paths in concrete/cement samples indicate the developed bilinear CZM techniques have the ability to capture crack nucleation and propagation in cement-based material samples with multiphase and associated interface. By comparing the computational prediction with the actual damaged samples, it also indicates that the ice crystallization pressure is the main mechanism for the internal frost damage in cementitious materials.

MULTISCALE EXAMINATION AND MODELING OF ELECTRON TRANSPORT IN NANOSCALE MATERIALS AND DEVICES

For half a century the integrated circuits (ICs) that make up the heart of electronic devices have been steadily improving by shrinking at an exponential rate. However, as the current crop of ICs get smaller and the insulating layers involved become thinner, electrons leak through due to quantum mechanical tunneling. This is one of several issues which will bring an end to this incredible streak of exponential improvement of this type of transistor device, after which future improvements will have to come from employing fundamentally different transistor architecture rather than fine tuning and miniaturizing the metal-oxide-semiconductor field effect transistors (MOSFETs) in use today. Several new transistor designs, some designed and built here at Michigan Tech, involve electrons tunneling their way through arrays of nanoparticles. We use a multi-scale approach to model these devices and study their behavior. For investigating the tunneling characteristics of the individual junctions, we use a first-principles approach to model conduction between sub-nanometer gold particles. To estimate the change in energy due to the movement of individual electrons, we use the finite element method to calculate electrostatic capacitances. The kinetic Monte Carlo method allows us to use our knowledge of these details to simulate the dynamics of an entire device— sometimes consisting of hundreds of individual particles—and watch as a device ‘turns on’ and starts conducting an electric current. Scanning tunneling microscopy (STM) and the closely related scanning tunneling spectroscopy (STS) are a family of powerful experimental techniques that allow for the probing and imaging of surfaces and molecules at atomic resolution. However, interpretation of the results often requires comparison with theoretical and computational models. We have developed a new method for calculating STM topographs and STS spectra. This method combines an established method for approximating the geometric variation of the electronic density of states, with a modern method for calculating spin-dependent tunneling currents, offering a unique balance between accuracy and accessibility.

Unique structures in a tumor herpesvirus revealed by cryo-electron tomography and microscopy.

Gammaherpesviruses, including the human pathogens Epstein-Barr virus and Kaposi's sarcoma-associated herpesvirus, are causative agents of lymphomas and other malignancies. The structural characterization of these viruses has been limited due to difficulties in obtaining adequate amount of virion particles. Here we report the first three-dimensional structural characterization of a whole gammaherpesvirus virion by an emerging integrated approach of cryo-electron tomography combined with single-particle cryo-electron microscopy, using murine gammaherpesvirus-68 (MHV-68) as a model system. We found that the MHV-68 virion consists of distinctive envelope and tegument compartments, and a highly conserved nucleocapsid. Two layers of tegument are identified: an inner tegument layer tethered to the underlying capsid and an outer, flexible tegument layer conforming to the overlying, pleomorphic envelope, consistent with the sequential viral tegumentation process inside host cells. Surprisingly, comparison of the MHV-68 virion and capsid reconstructions shows that the interactions between the capsid and inner tegument proteins are completely different from those observed in alpha and betaherpesviruses. These observations support the notion that the inner layer tegument across different subfamilies of herpesviruses has evolved significantly to confer specific characteristics related to viral-host interactions, in contrast to a highly conserved capsid for genome encapsidation and protection.

3.88 A structure of cytoplasmic polyhedrosis virus by cryo-electron microscopy.

Cytoplasmic polyhedrosis virus (CPV) is unique within the Reoviridae family in having a turreted single-layer capsid contained within polyhedrin inclusion bodies, yet being fully capable of cell entry and endogenous RNA transcription. Biochemical data have shown that the amino-terminal 79 residues of the CPV turret protein (TP) is sufficient to bring CPV or engineered proteins into the polyhedrin matrix for micro-encapsulation. Here we report the three-dimensional structure of CPV at 3.88 A resolution using single-particle cryo-electron microscopy. Our map clearly shows the turns and deep grooves of alpha-helices, the strand separation in beta-sheets, and densities for loops and many bulky side chains; thus permitting atomic model-building effort from cryo-electron microscopy maps. We observed a helix-to-beta-hairpin conformational change between the two conformational states of the capsid shell protein in the region directly interacting with genomic RNA. We have also discovered a messenger RNA release hole coupled with the mRNA capping machinery unique to CPV. Furthermore, we have identified the polyhedrin-binding domain, a structure that has potential in nanobiotechnology applications.

Three-dimensional electron beam dose calculations

The MDAH pencil-beam algorithm developed by Hogstrom et al (1981) has been widely used in clinics for electron beam dose calculations for radiotherapy treatment planning. The primary objective of this research was to address several deficiencies of that algorithm and to develop an enhanced version. Two enhancements have been incorporated into the pencil-beam algorithm; one models fluence rather than planar fluence, and the other models the bremsstrahlung dose using measured beam data. Comparisons of the resulting calculated dose distributions with measured dose distributions for several test phantoms have been made. From these results it is concluded (1) that the fluence-based algorithm is more accurate to use for the dose calculation in an inhomogeneous slab phantom, and (2) the fluence-based calculation provides only a limited improvement to the accuracy the calculated dose in the region just downstream of the lateral edge of an inhomogeneity. The source of the latter inaccuracy is believed primarily due to assumptions made in the pencil beam's modeling of the complex phantom or patient geometry.^ A pencil-beam redefinition model was developed for the calculation of electron beam dose distributions in three dimensions. The primary aim of this redefinition model was to solve the dosimetry problem presented by deep inhomogeneities, which was the major deficiency of the enhanced version of the MDAH pencil-beam algorithm. The pencil-beam redefinition model is based on the theory of electron transport by redefining the pencil beams at each layer of the medium. The unique approach of this model is that all the physical parameters of a given pencil beam are characterized for multiple energy bins. Comparisons of the calculated dose distributions with measured dose distributions for a homogeneous water phantom and for phantoms with deep inhomogeneities have been made. From these results it is concluded that the redefinition algorithm is superior to the conventional, fluence-based, pencil-beam algorithm, especially in predicting the dose distribution downstream of a local inhomogeneity. The accuracy of this algorithm appears sufficient for clinical use, and the algorithm is structured for future expansion of the physical model if required for site specific treatment planning problems. ^

Optimum Electron Temperature and Density for Short-Wavelength Plasma-Lasing from Nickel-like ions

Soft X-ray lasing across a Ni-like plasma gain-medium requires optimum electron temperature and density for attaining to the Ni-like ion stage and for population inversion in the View the MathML source3d94d1(J=0)→3d94p1(J=1) laser transition. Various scaling laws, function of operating parameters, were compared with respect to their predictions for optimum temperatures and densities. It is shown that the widely adopted local thermodynamic equilibrium (LTE) model underestimates the optimum plasma-lasing conditions. On the other hand, non-LTE models, especially when complemented with dielectronic recombination, provided accurate prediction of the optimum plasma-lasing conditions. It is further shown that, for targets with Z equal or greater than the rare-earth elements (e.g. Sm), the optimum electron density for plasma-lasing is not accessible for pump-pulses at View the MathML sourceλ=1ω=1μm. This observation explains a fundamental difficulty in saturating the wavelength of plasma-based X-ray lasers below 6.8 nm, unless using 2ω2ω pumping.

Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Respiration rates and electron transport system activities of copepod species obtained during Maria S. Merian cruise MSM17/3

Respiration rates and electron transport system (ETS) activities were measured in dominant copepod species from the northern Benguela upwelling system in January-February 2011 to assess the accuracy of the ETS assay in predicting in vivo respiration rates. Individual respiration rates varied from 0.06 to 1.60 µL O2/h/ind, while ETS activities converted to oxygen consumption ranged from 0.14 to 4.46 µL O2/h/ind. ETS activities were significantly correlated with respiration rates (r**2 = 0.79, p = 0.0001). R:ETS ratios were lowest in slow-moving Eucalanidae (0.11) and highest in diapausing Calanoides carinatus copepodids CV (0.76) while fast-moving copepods showed intermediate R:ETS (0.23-0.37). 82% of the variance of respiration rates could be explained by differences in dry mass, temperature and the activity level of different copepod species. Three regression equations were derived to calculate respiration rates for diapausing, slow- and fast-moving copepods, respectively, based on parameters such as body mass and temperature. Thus, knowledge about the activity level and behavioral characteristics of copepod species can significantly increase the predictive accuracy of metabolic models, which will help to better understand and quantify the impact of copepods on nutrient and carbon fluxes in marine ecosystems.

(Table 7) Representative electron microprobe analyses for natroalunite and larger FeOx spheroids of experiment ADSU6

Acid-sulfate alteration of basalt by SO2-bearing volcanic vapors has been proposed as one possible origin for sulfate-rich deposits on Mars. To better define mineralogical signatures of acid-sulfate alteration, laboratory experiments were performed to investigate alteration pathways and geochemical processes during reaction of basalt with sulfuric acid. Pyroclastic cinders composed of phenocrysts including plagioclase, olivine, and augite embedded in glass were reacted with sulfuric acid at 145 °C for up to 137 days at a range of fluid : rock ratios. During the experiments, the phenocrysts reacted rapidly to form secondary products, while the glass was unreactive. Major products included amorphous silica, anhydrite, and Fe-rich natroalunite, along with minor iron oxides/oxyhydroxides (probably hematite) and trace levels of other sulfates. At the lowest fluid : rock ratio, hexahydrite and an unidentified Fe-silicate phase also occurred as major products. Reaction-path models indicated that formation of the products required both slow dissolution of glass and kinetic inhibitions to precipitation of a number of minerals including phyllosilicates and other aluminosilicates as well as Al- and Fe-oxides/oxyhydroxides. Similar models performed for Martian basalt compositions predict that the initial stages of acid-sulfate alteration of pyroclastic deposits on Mars should result in formation of amorphous silica, anhydrite, Fe-bearing natroalunite, and kieserite, along with relict basaltic glass. In addition, analysis of the experimental products indicates that Fe-bearing natroalunite produces a Mössbauer spectrum closely resembling that of jarosite, suggesting that it should be considered an alternative to the component in sulfate-rich bedrocks at Meridiani Planum that has previously been identified as jarosite.