Modelling coupled physical-biogeochemical processes in ice-covered oceans

The last decades have seen a large effort of the scientific community to study and understand the physics of sea ice. We currently have a wide - even though still not exhaustive - knowledge of the sea ice dynamics and thermodynamics and of their temporal and spatial variability. Sea ice biogeochemistry is instead largely unknown. Sea ice algae production may account for up to 25% of overall primary production in ice-covered waters of the Southern Ocean. However, the influence of physical factors, such as the location of ice formation, the role of snow cover and light availability on sea ice primary production is poorly understood. There are only sparse localized observations and little knowledge of the functioning of sea ice biogeochemistry at larger scales. Modelling becomes then an auxiliary tool to help qualifying and quantifying the role of sea ice biogeochemistry in the ocean dynamics. In this thesis, a novel approach is used for the modelling and coupling of sea ice biogeochemistry - and in particular its primary production - to sea ice physics. Previous attempts were based on the coupling of rather complex sea ice physical models to empirical or relatively simple biological or biogeochemical models. The focus is moved here to a more biologically-oriented point of view. A simple, however comprehensive, physical model of the sea ice thermodynamics (ESIM) was developed and coupled to a novel sea ice implementation (BFM-SI) of the Biogeochemical Flux Model (BFM). The BFM is a comprehensive model, largely used and validated in the open ocean environment and in regional seas. The physical model has been developed having in mind the biogeochemical properties of sea ice and the physical inputs required to model sea ice biogeochemistry. The central concept of the coupling is the modelling of the Biologically-Active-Layer (BAL), which is the time-varying fraction of sea ice that is continuously connected to the ocean via brines pockets and channels and it acts as rich habitat for many microorganisms. The physical model provides the key physical properties of the BAL (e.g., brines volume, temperature and salinity), and the BFM-SI simulates the physiological and ecological response of the biological community to the physical enviroment. The new biogeochemical model is also coupled to the pelagic BFM through the exchange of organic and inorganic matter at the boundaries between the two systems . This is done by computing the entrapment of matter and gases when sea ice grows and release to the ocean when sea ice melts to ensure mass conservation. The model was tested in different ice-covered regions of the world ocean to test the generality of the parameterizations. The focus was particularly on the regions of landfast ice, where primary production is generally large. The implementation of the BFM in sea ice and the coupling structure in General Circulation Models will add a new component to the latters (and in general to Earth System Models), which will be able to provide adequate estimate of the role and importance of sea ice biogeochemistry in the global carbon cycle.

Development of Clinical Decision Support Systems based on Mathematical Models of Physiological Systems

In the last years of research, I focused my studies on different physiological problems. Together with my supervisors, I developed/improved different mathematical models in order to create valid tools useful for a better understanding of important clinical issues. The aim of all this work is to develop tools for learning and understanding cardiac and cerebrovascular physiology as well as pathology, generating research questions and developing clinical decision support systems useful for intensive care unit patients. I. ICP-model Designed for Medical Education We developed a comprehensive cerebral blood flow and intracranial pressure model to simulate and study the complex interactions in cerebrovascular dynamics caused by multiple simultaneous alterations, including normal and abnormal functional states of auto-regulation of the brain. Individual published equations (derived from prior animal and human studies) were implemented into a comprehensive simulation program. Included in the normal physiological modelling was: intracranial pressure, cerebral blood flow, blood pressure, and carbon dioxide (CO2) partial pressure. We also added external and pathological perturbations, such as head up position and intracranial haemorrhage. The model performed clinically realistically given inputs of published traumatized patients, and cases encountered by clinicians. The pulsatile nature of the output graphics was easy for clinicians to interpret. The manoeuvres simulated include changes of basic physiological inputs (e.g. blood pressure, central venous pressure, CO2 tension, head up position, and respiratory effects on vascular pressures) as well as pathological inputs (e.g. acute intracranial bleeding, and obstruction of cerebrospinal outflow). Based on the results, we believe the model would be useful to teach complex relationships of brain haemodynamics and study clinical research questions such as the optimal head-up position, the effects of intracranial haemorrhage on cerebral haemodynamics, as well as the best CO2 concentration to reach the optimal compromise between intracranial pressure and perfusion. We believe this model would be useful for both beginners and advanced learners. It could be used by practicing clinicians to model individual patients (entering the effects of needed clinical manipulations, and then running the model to test for optimal combinations of therapeutic manoeuvres). II. A Heterogeneous Cerebrovascular Mathematical Model Cerebrovascular pathologies are extremely complex, due to the multitude of factors acting simultaneously on cerebral haemodynamics. In this work, the mathematical model of cerebral haemodynamics and intracranial pressure dynamics, described in the point I, is extended to account for heterogeneity in cerebral blood flow. The model includes the Circle of Willis, six regional districts independently regulated by autoregulation and CO2 reactivity, distal cortical anastomoses, venous circulation, the cerebrospinal fluid circulation, and the intracranial pressure-volume relationship. Results agree with data in the literature and highlight the existence of a monotonic relationship between transient hyperemic response and the autoregulation gain. During unilateral internal carotid artery stenosis, local blood flow regulation is progressively lost in the ipsilateral territory with the presence of a steal phenomenon, while the anterior communicating artery plays the major role to redistribute the available blood flow. Conversely, distal collateral circulation plays a major role during unilateral occlusion of the middle cerebral artery. In conclusion, the model is able to reproduce several different pathological conditions characterized by heterogeneity in cerebrovascular haemodynamics and can not only explain generalized results in terms of physiological mechanisms involved, but also, by individualizing parameters, may represent a valuable tool to help with difficult clinical decisions. III. Effect of Cushing Response on Systemic Arterial Pressure. During cerebral hypoxic conditions, the sympathetic system causes an increase in arterial pressure (Cushing response), creating a link between the cerebral and the systemic circulation. This work investigates the complex relationships among cerebrovascular dynamics, intracranial pressure, Cushing response, and short-term systemic regulation, during plateau waves, by means of an original mathematical model. The model incorporates the pulsating heart, the pulmonary circulation and the systemic circulation, with an accurate description of the cerebral circulation and the intracranial pressure dynamics (same model as in the first paragraph). Various regulatory mechanisms are included: cerebral autoregulation, local blood flow control by oxygen (O2) and/or CO2 changes, sympathetic and vagal regulation of cardiovascular parameters by several reflex mechanisms (chemoreceptors, lung-stretch receptors, baroreceptors). The Cushing response has been described assuming a dramatic increase in sympathetic activity to vessels during a fall in brain O2 delivery. With this assumption, the model is able to simulate the cardiovascular effects experimentally observed when intracranial pressure is artificially elevated and maintained at constant level (arterial pressure increase and bradicardia). According to the model, these effects arise from the interaction between the Cushing response and the baroreflex response (secondary to arterial pressure increase). Then, patients with severe head injury have been simulated by reducing intracranial compliance and cerebrospinal fluid reabsorption. With these changes, oscillations with plateau waves developed. In these conditions, model results indicate that the Cushing response may have both positive effects, reducing the duration of the plateau phase via an increase in cerebral perfusion pressure, and negative effects, increasing the intracranial pressure plateau level, with a risk of greater compression of the cerebral vessels. This model may be of value to assist clinicians in finding the balance between clinical benefits of the Cushing response and its shortcomings. IV. Comprehensive Cardiopulmonary Simulation Model for the Analysis of Hypercapnic Respiratory Failure We developed a new comprehensive cardiopulmonary model that takes into account the mutual interactions between the cardiovascular and the respiratory systems along with their short-term regulatory mechanisms. The model includes the heart, systemic and pulmonary circulations, lung mechanics, gas exchange and transport equations, and cardio-ventilatory control. Results show good agreement with published patient data in case of normoxic and hyperoxic hypercapnia simulations. In particular, simulations predict a moderate increase in mean systemic arterial pressure and heart rate, with almost no change in cardiac output, paralleled by a relevant increase in minute ventilation, tidal volume and respiratory rate. The model can represent a valid tool for clinical practice and medical research, providing an alternative way to experience-based clinical decisions. In conclusion, models are not only capable of summarizing current knowledge, but also identifying missing knowledge. In the former case they can serve as training aids for teaching the operation of complex systems, especially if the model can be used to demonstrate the outcome of experiments. In the latter case they generate experiments to be performed to gather the missing data.

Genesis of metasomatic sapphirine-corundum-spinel-bearing granulites in Sri Lanka

The goal of the present study is to understand the mechanism of mass transfer, the composition and the role of fluids during crustal metasomatism in high-temperature metamorphic terranes. A well constrained case study, a locality at Rupaha, Sri Lanka was selected. It is located in the Highland Complex of Sri Lanka, which represents a small, but important fragment of the super-continent Gondwana. Excellent exposures of ultramafic rocks, which are embedded in granulites, were found at 10 localities. These provide a unique background for understanding the metasomatic processes. The boundary between the ultramafic and the granulite rocks are lined with metasomatic reaction zones up to 50cm in width. Progressing from the ultramafics to the granulite host rock, three distinct zones with the following mineral assemblages can be distinguished: (1). phlogopite + spinel + sapphirine, (2). spinel + sapphirine and (3). corundum + biotite + plagioclase. In order to assess the P-T-t path, the peak metamorphism and the exhumation history were constrained using different thermobarometers, as well as a diffusion model of garnet zoning. A maximum temperature of 875 ± 20oC (Opx-Cpx thermometer) and at the peak pressure of 9.0 ± 0.1 kbar (Grt-Cpx-Pl-Qtz) was calculated for the silicic granulite. The ultramafic rocks recorded a peak temperature of 840 ± 70oC (Opx-Cpx thermometer) at 9 kbar. Coexisting spinel and sapphirine from the reaction zone yield a temperature of 820 ± 40oC. This is in agreement with the peak-temperatures recorded in the adjacent granulites and ultramafics rocks. The structural concordance of the ultramafic rocks with the siliceous granulite host rock further support the suggestion, that all units have experienced the same peak metamorphism. Diffusion modeling of retrograde zoning in garnets from mafic granulites suggests a three-step cooling history. A maximum cooling rate of 1oC/Ma is estimated during the initial stage of cooling, followed by a cooling rate of ~30oC/Ma. The outermost rims of garnet indicate a slightly slower cooling rate at about 10-15oC/Ma. The sequences of mineral zones, containing a variety of Al-rich, silica undersaturated minerals in the reaction zones separating the ultramafic rocks from the silica-rich rocks can be explained by a diffusion model. This involves the diffusion of Mg from ultramafic rocks across the layers, and K and Si diffuse in opposite direction. Chemical potential of Mg and Si generated continuous monotonic gradient, allowing steady state diffusional transport across the profile. The strong enrichment in Al, and the considerable loss of Si, during the formation of reaction bands can be inferred from isocon diagrams. Some Al was probably added to the reaction zones, while Si was lost. This is most likely due to fluids percolating parallel to the zones at the boundary of the rock units. This study has shown that not only pressure and temperature conditions but most importantly PH2O and the concentration of the chlorine and fluorine in aqueous fluids also control the mass transport in different geological environments.

2 H-NMR und Viskositätsuntersuchungen zum Studium der molekularen Dynamik unterkühlter Flüssigkeiten unter Hochdruck

In der vorliegenden Arbeit wurde die Druckabhängigkeit der molekularen Dynamik mittels 2H-NMR und Viskositätsmessungen untersucht. Für die Messungen wurde der niedermolekulare organische Glasbildner ortho-Terphenyl (OTP) ausgewählt, da dieser aufgrund einer Vielzahl vorliegender Arbeiten als Modellsubstanz angesehen werden kann. Daneben wurden auch Messungen an Salol durchgeführt.Die Untersuchungen erstreckten sich über einen weiten Druck- und Temperaturbereich ausgehend von der Schmelze bis weit in die unterkühlte Flüssigkeit. Dieser Bereich wurde aufgrund experimenteller Voraussetzungen immer durch eine Druckerhöhung erreicht.Beide Substanzen zeigten druckabhängig ein Verhalten, das dem der Temperaturvariation bei Normaldruck sehr ähnelt. Auf einer Zeitskala der molekularen Dynamik von 10E-9 s bis zu 10E+2 s wurde daher am Beispiel von OTP ein Druck-Temperatur-Zeit-Superpositionsprinzip diskutiert. Zudem konnte eine Temperatur-Dichte-Skalierung mit rho T-1/4 erfolgreich durchgeführt werden. Dies entspricht einem rein repulsiven Potentialverlauf mit rho -12±3 .Zur Entscheidung, ob die Verteilungsbreiten der mittleren Rotationskorrelationszeiten durch Druckvariation beeinflußt werden, wurden auch Ergebnisse anderer experimenteller Methoden herangezogen. Unter Hinzuziehung aller Meßergebnisse kann sowohl eine Temperatur- als auch Druckabhängigkeit der Verteilungsbreite bestätigt werden. Zur Auswertung von Viskositätsdaten wurde ein Verfahren vorgestellt, das eine quantitative Aussage über den Fragilitätsindex von unterkühlten Flüssigkeiten auch dann zuläßt, wenn die Messungen nicht bis zur Glasübergangstemperatur Tg durchgeführt werden. Die Auswertung der druckabhängigen Viskositätsdaten von OTP und Salol zeigt einen sehr differenzierten druckabhängigen Verlauf des Fragilitätsindexes für beide Glasbildner. OTP zeigt zunächst eine leichte Abnahme und danach wieder eine Zunahme des Fragilitätsindexes, dieses Ergebnis wird auch von Simulationsdaten, die der Literatur entnommen wurden, unterstützt. Salol hingegen zeigt zunächst eine deutliche Zunahme und danach eine Abnahme des Fragilitätsindexes. Das unterschiedliche Verhalten der beiden Glasbildner mit ähnlichem Fragilitätsindex bei Normaldruck wird auf die Wasserstoffbrückenbindungen innerhalb von Salol zurückgeführt.

Phasenverhalten von Polyethylen in Mischsystemen

Die vorliegende Arbeit beschäftigt sich mit dem Phasenverhalten von Polyethylen (PE) in nicht-reaktiven und in reaktiven Systemen. Von drei eng verteilten Polyethylenen (Mw = 6,4, 82 bzw. 380 kg/mol) in n-Hexan sowie für das System 2,2-Dimethylbutan / PE 82 wurde die Entmischung in Abhängigkeit von der Zusammensetzung, dem Druck und der Temperatur experimentell bestimmt. Die Modellierung der Trübungskurven erfolgte nach der Theorie von Sanchez und Lacombe. Dieser Ansatz beschreibt die Ergebnisse qualitativ und kann in einem engen Temperatur- und Druckbereich für gegebenes Molekulargewicht die kritische Temperatur und den kritischen Druck quantitativ vorhersagen. Durch Extrapolation der kritischen Temperatur der verschiedenen Lösungen von PE in n-Hexan auf unendliches Molekulargewicht nach Shultz-Flory wurde im Druckbereich von 20 bis 100 bar und im Temperaturbereich von 130 bis 200 °C eine Grenzlinie bestimmt. Diese Linie trennt unabhängig vom Molekulargewicht des Polymers und der Zusammensetzung der Mischung das Zweiphasengebiet vom homogenen Bereich. Im Fall des Mischlösungsmittels n-Hexan / 2,2-Dimethylbutan wurde für eine annähernd kritische Polymerkonzentration die Abhängigkeit der Entmischungsbedingungen von der Zusammensetzung untersucht. Durch einfache Erweiterung der Sanchez-Lacombe-Theorie und Einführen eines Fitparameters konnte das ternäre System beschrieben werden. An einer breit verteilten PE-Probe wurden Experimente zur Fraktionierung von PE in n-Hexan durchgeführt. Die Analyse der in den koexistenten Phasen enthaltenen Polymere lieferte Informationen über die Konzentration und die Molekulargewichtsverteilung des PE in diesen Phasen sowie die kritische Zusammensetzung der Mischung. Von verschiedenen PE-Lösungen (Mw = 0,5 kg/mol) wurde die polymerisationsinduzierte Phasenseparation in Isobornylmethacrylat mit und ohne Vernetzer untersucht. Mit 15 Gew.-% PE und in Abwesenheit von Vernetzer findet die Entmischung erst bei hohen Umsätzen statt. Die Charakterisierung der resultierenden Proben zeigte, dass sich etwas mehr als 5 Gew.-% PE im Polyisobornylmethacrylat lösen. Die Glasübergangstemperaturen der Polymermischungen steigen mit steigender Vernetzer- und sinkender Polyethylenkonzentration. Bei Proben mit 15 Gew.-% PE zeigte sich folgendes: 5 Gew.-% Vernetzer führen zu großen PE-Bereichen (150 - 200 nm) in der Matrix und der Kristallinitätsgrad ist gering. Bei der Polymermischung mit 10 Gew.-% Vernetzer bilden sich sehr kleine Polyethylenkristalle (< 80 nm) und der Kristallinitätsgrad ist hoch. Ohne Vernetzer hängt der Kristallinitätsgrad - wie bei reinem PE - von der Abkühlrate ab, mit Vernetzer ist er von ihr unabhängig.

Investigations on maar-diatreme volcanoes by inversion of magnetic and gravity data from the Eifel area, Germany

A study of maar-diatreme volcanoes has been perfomed by inversion of gravity and magnetic data. The geophysical inverse problem has been solved by means of the damped nonlinear least-squares method. To ensure stability and convergence of the solution of the inverse problem, a mathematical tool, consisting in data weighting and model scaling, has been worked out. Theoretical gravity and magnetic modeling of maar-diatreme volcanoes has been conducted in order to get information, which is used for a simple rough qualitative and/or quantitative interpretation. The information also serves as a priori information to design models for the inversion and/or to assist the interpretation of inversion results. The results of theoretical modeling have been used to roughly estimate the heights and the dip angles of the walls of eight Eifel maar-diatremes — each taken as a whole. Inversemodeling has been conducted for the Schönfeld Maar (magnetics) and the Hausten-Morswiesen Maar (gravity and magnetics). The geometrical parameters of these maars, as well as the density and magnetic properties of the rocks filling them, have been estimated. For a reliable interpretation of the inversion results, beside the knowledge from theoretical modeling, it was resorted to other tools such like field transformations and spectral analysis for complementary information. Geologic models, based on thesynthesis of the respective interpretation results, are presented for the two maars mentioned above. The results gave more insight into the genesis, physics and posteruptive development of the maar-diatreme volcanoes. A classification of the maar-diatreme volcanoes into three main types has been elaborated. Relatively high magnetic anomalies are indicative of scoria cones embeded within maar-diatremes if they are not caused by a strong remanent component of the magnetization. Smaller (weaker) secondary gravity and magnetic anomalies on the background of the main anomaly of a maar-diatreme — especially in the boundary areas — are indicative for subsidence processes, which probably occurred in the late sedimentation phase of the posteruptive development. Contrary to postulates referring to kimberlite pipes, there exists no generalized systematics between diameter and height nor between geophysical anomaly and the dimensions of the maar-diatreme volcanoes. Although both maar-diatreme volcanoes and kimberlite pipes are products of phreatomagmatism, they probably formed in different thermodynamic and hydrogeological environments. In the case of kimberlite pipes, large amounts of magma and groundwater, certainly supplied by deep and large reservoirs, interacted under high pressure and temperature conditions. This led to a long period phreatomagmatic process and hence to the formation of large structures. Concerning the maar-diatreme and tuff-ring-diatreme volcanoes, the phreatomagmatic process takes place due to an interaction between magma from small and shallow magma chambers (probably segregated magmas) and small amounts of near-surface groundwater under low pressure and temperature conditions. This leads to shorter time eruptions and consequently to structures of smaller size in comparison with kimberlite pipes. Nevertheless, the results show that the diameter to height ratio for 50% of the studied maar-diatremes is around 1, whereby the dip angle of the diatreme walls is similar to that of the kimberlite pipes and lies between 70 and 85°. Note that these numerical characteristics, especially the dip angle, hold for the maars the diatremes of which — estimated by modeling — have the shape of a truncated cone. This indicates that the diatreme can not be completely resolved by inversion.

An investigation into complex inorganic materials with Mössbauer spectroscopy

Mössbauer Spektroskopie ist ein unverzichtbares Instrument für die Bestimmung von Oxidationszuständen und für die Analyse von lokalen Ordnungsphänomenen von Mössbauer aktiven Atomen. Weil es sich um eine lokale Methode handelt können sowohl kristalline als auch amorphe Materialien untersucht werden. Die Kombination von lokaler Prüfung mit Mössbauer Spektroskopie und globaler Untersuchung z.B. mit Röntgendiffraktometrie ermöglicht die Studie von Ordnungseffekten von statistisch besetzten Positionen in einer geordneten Matrix. Das wurde hier eingesetzt um die lokale Umgebung in zwei Serien von Heuslerverbindungen, Co2-xFe1+xSi and Co2Mn1-xFexAl zu untersuchen. Für die Co2Mn1-xFexAl Serie wurde eine L21 geordnete Phase in einer insgesamt B2 geordneten Probe detektiert. Ein Wechsel von der AlCu2Mn zu der CuHg2Ti Struktur wurde für die Co2-xFe1+xSi Proben gefunden. Die Transformation von einem Glas zu einem keramischen Material wurde mit 119Sn Mössbauer Spektroskopie untersucht. Die höhere Ordnung in der Keramik wurde von einer kleiner werdenden Mössbauerlinienbreite begleitet. Demzufolge geben die Modifikationen der Sn Umgebungen klar die Transformation des gesamten Materials wieder. Ist die lokale Umgebung von unregelmäßig auftretenden Atomen in einer amorphen Matrix von Interesse, sind lokal prüfende Methoden die zuverlässigsten Methoden die zur Verfügung stehen. In dieser Arbeit wurde 119Sn Mössbauer Spektroskopie eingesetzt um die Oxidationszustände, die lokalen Umgebungen und relativen Intensitäten von Zinn Atomen in einer Silikatmatrix zu bestimmen. Modifikationen dieser Parameter als Funktion von Prozess bestimmenden Parametern wie der Sauerstoffpartialdruck, die Temperatur, die Behandlungsdauer und der Abkühlprozess genauso wie der SnO2 Gehalt sind von Interesse, weil durch Reduktions- und Diffusionsprozesse Änderungen des Koordinations- und des Oxidationszustands der Zinnatome auftreten. Da diese Änderungen in der Glasmatrix verursachen, die das fertige Produkt im industriellen Fertigungsprozess ruinieren können sind diese feinen Veränderungen sehr wichtig. Wenigstens zwei Mössbauerlinien korrespondierend mit zwei verschiedenen Umgebungen für Sn2+ und Sn4+ sind für eine Analyse mit ausreichender Qualität notwendig. Durch Vergleich von den bestimmten Hyperfein Parametern mit den Parametern von Modelsubstanzen werden lokale Umgebungen der Zinnatome entworfen. Für Sn2+ werden zwei auf einer trigonalen Pyramide basierende Umgebungen mit variierender Anzahl von bindenden und nicht-bindenden Sauerstoffatomen formuliert. Für Sn4+ wurde eine tetraedrische und eine oktaedrische Umgebung postuliert. Die relativen Intensitäten der vier Mössbauerlinien wurden um ein Diffusions- und Reaktionsmodell zu entwickeln und um einen Satz von Diffusions- und Transferkoeffizienten zu bestimmen eingesetzt. Die bestimmten Diffusionskoeffizienten stimmen mit den Literaturdaten überein. Der Massentransferkoeffizient ist kleiner als der bestimmte Wert, aber immer noch in der gleichen Größenordnung. Im Gegensatz zu den Erwartungen ist der präsentierte Diffusionskoeffizient für Sn4+ bestimmt als der von Sn2+. Das wiederum kann durch Berücksichtigung von Elektronhoppingprozessen erklärt werden.

Atmospheric heterogeneous reactions of N 2 O 5 and NO 3 radicals with mineral dust particles

Die heterogenen Reaktionen von N2O5 bzw. NO3 auf mineralischen Staubpartikeln wurden untersucht, um deren Einfluss auf den Abbau atmosphärischer Stickoxide (NOx) sowie auf die chemische Veränderung der Staubpartikel während ihres Transportes durch die Atmosphäre besser verstehen zu können. Die experimentellen Studien wurden bei Atmosphärendruck, Raumtemperatur und unterschiedlichen relativen Luftfeuchten durchgeführt. Der Aufnahmekoeffizient γ(N2O5) von N2O5 auf dispergiertem Staub aus der Sahara wurde zu 0,020 ± 0,002 (1σ) bestimmt, unabhängig von der relativen Feuchte (0 - 67 %) sowie der N2O5-Konzentration (5x1011 - 3x1013 Moleküle cm-3).rnDie Analyse der Reaktionsprodukte in der Gasphase sowie auf der Partikeloberfläche führt zu der Annahme, dass N2O5 auf der Staubpartikeloberfläche zu Nitrat hydrolysiert wird. Es konnte kein Einfluss der relativen Feuchte auf den Aufnahmekoeffizienten ermittelt werden, was durch das vorhandene interlamellare Wasser, welches bis zu 10 % der Partikelmasse betragen kann, erklärbar ist. Der gemessene Wert des Aufnahmekoeffizienten ist unabhängig von der Eingangs-N2O5-Konzentration, was sich über die sehr große innere Oberfläche der Partikel erklären lässt. Dennoch ließ sich durch eine vorherige Konditionierung der Partikel mit gasförmigem HNO3, was eine Nitratanreicherung an der Oberfläche bewirkt, die Effizienz der N2O5-Aufnahme auf die Staubpartikel reduzieren. Zusätzliche Studien befassten sich mit der Bestimmung des Aufnahmekoeffizienten von N2O5 auf Illit-Partikeln und auf Teststaub aus Arizona. Bei einer relativen Luftfeuchte von 0 % wurden für γ(N2O5) Werte von 0,084 ± 0,019 (1σ) für Illit und von 0,010 ± 0,001 (1σ) für Arizona Teststaub ermittelt.rnUnter Anwendung einer neuartigen Messmethode, die auf der zeitgleichen Messung der Konzentrationsabnahme von NO3 und N2O5 relativ zueinander beruht, wurde das Verhältnis γ(NO3)/γ(N2O5) der Aufnahmekoeffizienten von NO3 und N2O5 auf Saharastaub zu 0,9 ± 0,4 (1σ) bestimmt. Dieser Wert war unabhängig von der relativen Feuchte, den NO3- und N2O5-Konzentrationen sowie der Reaktionszeit, obwohl eine Oberflächendeaktivierung für beide Spurenstoffe beobachtet wurde.

Prediction of new crystal structures under extreme conditions

One of the most important challenges in chemistry and material science is the connection between the contents of a compound and its chemical and physical properties. In solids, these are greatly influenced by the crystal structure.rnrnThe prediction of hitherto unknown crystal structures with regard to external conditions like pressure and temperature is therefore one of the most important goals to achieve in theoretical chemistry. The stable structure of a compound is the global minimum of the potential energy surface, which is the high dimensional representation of the enthalpy of the investigated system with respect to its structural parameters. The fact that the complexity of the problem grows exponentially with the system size is the reason why it can only be solved via heuristic strategies.rnrnImprovements to the artificial bee colony method, where the local exploration of the potential energy surface is done by a high number of independent walkers, are developed and implemented. This results in an improved communication scheme between these walkers. This directs the search towards the most promising areas of the potential energy surface.rnrnThe minima hopping method uses short molecular dynamics simulations at elevated temperatures to direct the structure search from one local minimum of the potential energy surface to the next. A modification, where the local information around each minimum is extracted and used in an optimization of the search direction, is developed and implemented. Our method uses this local information to increase the probability of finding new, lower local minima. This leads to an enhanced performance in the global optimization algorithm.rnrnHydrogen is a highly relevant system, due to the possibility of finding a metallic phase and even superconductor with a high critical temperature. An application of a structure prediction method on SiH12 finds stable crystal structures in this material. Additionally, it becomes metallic at relatively low pressures.

Origin of the high upconversion green luminescence efficiency in -NaYF4:2% Er3+,20% Yb3+

Site-selective spectroscopy in hexagonal beta-NaYF4:Er3+,Yb3+ has revealed different environments for Er3+ ions (multisite formation). The low-temperature S-4(3/2) -> (I15/2Er3+)-I-4 green emission depends on the excitation wavelength associated with the F-4(7/2) Er3+ level. We have studied the effect of hydrostatic pressure on the green, red, and blue Er3+ emission upon NIR excitation at similar to 980 nm, in order to establish the role played by energy resonance conditions and the multiple Er3+ sites due to the disordered structure for the upconversion (UC) process (energy tuning). The variation of photoluminescence spectra and lifetimes as a function of pressure and temperature reveals that the origin of the high green UC efficiency of the beta-NaYF4:Er3+,Yb3+ compound is mainly due to the multisite distribution, and the low phonon energy of the host lattice.

A novel interface for hybrid mock circulations to evaluate ventricular assist devices

This paper presents a novel mock circulation for the evaluation of ventricular assist devices (VADs), which is based on a hardware-in-the-loop concept. A numerical model of the human blood circulation runs in real time and computes instantaneous pressure, volume, and flow rate values. The VAD to be tested is connected to a numerical-hydraulic interface, which allows the interaction between the VAD and the numerical model of the circulation. The numerical-hydraulic interface consists of two pressure-controlled reservoirs, which apply the computed pressure values from the model to the VAD, and a flow probe to feed the resulting VAD flow rate back to the model. Experimental results are provided to show the proper interaction between a numerical model of the circulation and a mixed-flow blood pump.

Changes in respiratory mechanics and gas exchange during the acute respiratory distress syndrome

BACKGROUND: The time course of impairment of respiratory mechanics and gas exchange in the acute respiratory distress syndrome (ARDS) remains poorly defined. We assessed the changes in respiratory mechanics and gas exchange during ARDS. We hypothesized that due to the changes in respiratory mechanics over time, ventilatory strategies based on rigid volume or pressure limits might fail to prevent overdistension throughout the disease process. METHODS: Seventeen severe ARDS patients {PaO2/FiO2 10.1 (9.2-14.3) kPa; 76 (69-107) mmHg [median (25th-75th percentiles)] and bilateral infiltrates} were studied during the acute, intermediate, and late stages of ARDS (at 1-3, 4-6 and 7 days after diagnosis). Severity of lung injury, gas exchange, and hemodynamics were assessed. Pressure-volume (PV) curves of the respiratory system were obtained, and upper and lower inflection points (UIP, LIP) and recruitment were estimated. RESULTS: (1) UIP decreased from early to established (intermediate and late) ARDS [30 (28-30) cmH2O, 27 (25-30) cmH2O and 25 (23-28) cmH2O (P=0.014)]; (2) oxygenation improved in survivors and in patients with non-pulmonary etiology in late ARDS, whereas all patients developed hypercapnia from early to established ARDS; and (3) dead-space ventilation and pulmonary shunt were larger in patients with pulmonary etiology during late ARDS. CONCLUSION: We found a decrease in UIP from acute to established ARDS. If applied to our data, the inspiratory pressure limit advocated by the ARDSnet (30 cmH2O) would produce ventilation over the UIP, with a consequent increased risk of overdistension in 12%, 43% and 65% of our patients during the acute, intermediate and late phases of ARDS, respectively. Lung protective strategies based on fixed tidal volume or pressure limits may thus not fully avoid the risk of lung overdistension throughout ARDS.

ECG-triggered muscular counterpulsation for treatment of low cardiac output

BACKGROUND: Skeletal muscular counterpulsation (MCP) has been used as a new noninvasive technique for treatment of low cardiac output. The MCP method is based on ECG-triggered skeletal muscle stimulation. The purpose of the present study was to evaluate acute hemodynamic changes induced by MCP in the experimental animal. METHODS: Eight anaesthetized pigs (43+/-4 kg) were studied at rest and after IV â-blockade (10 mg propranolol) before and after MCP. Muscular counterpulsation was performed on both thighs using trains (75 ms duration) of multiple biphasic electrical impulses with a width of 1 ms and a frequency of 200 Hz at low (10 V) and high (30 V) amplitude. ECG-triggering was used to synchronize stimulation to a given time point. LV pressure-volume relations were determined using the conductance catheter. After baseline measurements, MCP was carried out for 10 minutes at low and high stimulation amplitude. The optimal time point for MCP was determined from LV pressure-volume loops using different stimulation time points during systole and diastole. Best results were observed during end-systole and, therefore, this time point was used for stimulation. RESULTS: Under control conditions, MCP was associated with a significant decrease in pulmonary vascular resistance (-18%), a decrease in systemic vascular resistance (-11%) and stroke work index (-4%), whereas cardiac index (+2%) and ejection fraction (+6%) increased slightly. Pressure-volume loops showed a leftward shift with a decrease in end-systolic volume. After â-blockade, cardiac function decreased (HR, MAP, EF, dP/dt max), but it improved with skeletal muscle stimulation (HR +10% and CI +17%, EF +5%). There was a significant decrease in pulmonary (-19%) and systemic vascular resistance (-29%). CONCLUSIONS: In the animal model, ECG-triggered skeletal muscular counterpulsation is associated with a significant improvement in cardiac function at baseline and after IV â-blockade. Thus, MCP represents a new, non-invasive technique which improves cardiac function by diastolic compression of the peripheral arteries and veins, with a decrease in systemic vascular resistance and increase in cardiac output.

Transdermal nitroglycerin in patients with subarachnoid hemorrhage

Delayed ischemic neurological deficit (DIND) following cerebral vasospasm remains a cause for high morbidity and mortality in patients with subarachnoid hemorrhage (SAH). There is experimental and clinical evidence of positive effects of nitric oxide (NO) donors on cerebral vasospasm. We therefore analysed the effect of transdermal nitroglycerin in patients with SAH measuring transcranial Doppler velocities (TCD), cerebral blood flow (CBF) and DIND. Nitroglycerin was used in a target dose of 14 microg/kg/h. TCD assessment was performed daily. CBF measurements were done using the perfusion CT-technique. Blood pressure, volume intake and vasopressor administration, were registered. Nine patients were randomly assigned either to the nitroglycerin group (N-group) and eight patients in the control group (C-group). Mean TCD values in the extracranial portion of the internal carotid artery (ICA) were lower in the N-group (p<0.005). Mean TCD in the middle cerebral arteries (MCA) showed no difference. The Lindegaard ratio was higher in the N-group (p<0.04). CBF in the N-group was higher than in the C-group (p<0.03). Even though nitroglycerin reduces blood pressure and lowers ICA TCD-values and increases the Lindegaard ratio, a higher CBF was measured in the N-group. Thus, nitroglycerin influences the cerebral vascular tone and increases CBF. SAH therapy with nitroglycerin is possible without increasing the risk of DIND. The exact timing of onset, duration and reduction of nitroglycerin administration in respect to the appearance of vasospasm may have a strong impact on the success of such a therapy.

A cardiovascular mathematical model of graded head-up tilt

A lumped parameter model of the cardiovascular system has been developed and optimized using experimental data obtained from 13 healthy subjects during graded head-up tilt (HUT) from the supine position to [Formula: see text]. The model includes descriptions of the left and right heart, direct ventricular interaction through the septum and pericardium, the systemic and pulmonary circulations, nonlinear pressure volume relationship of the lower body compartment, arterial and cardiopulmonary baroreceptors, as well as autoregulatory mechanisms. A number of important features, including the separate effects of arterial and cardiopulmonary baroreflexes, and autoregulation in the lower body, as well as diastolic ventricular interaction through the pericardium have been included and tested for their significance. Furthermore, the individual effect of parameter associated with heart failure, including LV and RV contractility, baseline systemic vascular resistance, pulmonary vascular resistance, total blood volume, LV diastolic stiffness and reflex gain on HUT response have also been investigated. Our fitted model compares favorably with our experimental measurements and published literature at a range of tilt angles, in terms of both global and regional hemodynamic variables. Compared to the normal condition, a simulated congestive heart failure condition produced a blunted response to HUT with regards to the percentage changes in cardiac output, stroke volume, end diastolic volume and effector response (i.e., heart contractility, venous unstressed volume, systemic vascular resistance and heart rate) with progressive tilting.