Hybrid Modelling of Titan's Interaction with the Magnetosphere of Saturn

In this dissertation we study the interaction between Saturn's moon Titan and the magnetospheric plasma and magnetic field. The method of research is a three-dimensional computer simulation model, that is used to simulate this interaction. The simulation model used is a hybrid model. Hybrid models enable individual tracking or tracing of ions and also take into account the particle motion in the propagation of the electromagnetic fields. The hybrid model has been developed at the Finnish Meteorological Institute. This thesis gives a general description of the effects that the solar wind has on Earth and other planets of our solar system. Planetary satellites can also have similar interactions with the solar wind but also with the plasma flows of planetary magnetospheres. Titan is clearly the largest among the satellites of Saturn and also the only known satellite with a dense atmosphere. It is the atmosphere that makes Titan's plasma interaction with the magnetosphere of Saturn so unique. Nevertheless, comparisons with the plasma interactions of other solar system bodies are valuable. Detecting charged plasma particles requires in situ measurements obtainable through scientific spacecraft. The Cassini mission has been one of the most remarkable international efforts in space science. Since 2004 the measurements and images obtained from instruments onboard the Cassini spacecraft have increased the scientific knowledge of Saturn as well as its satellites and magnetosphere in a way no one was probably able to predict. The current level of science on Titan is practically unthinkable without the Cassini mission. Many of the observations by Cassini instrument teams have influenced this research both the direct measurements of Titan as well as observations of its plasma environment. The theoretical principles of the hybrid modelling approach are presented in connection to the broader context of plasma simulations. The developed hybrid model is described in detail: e.g. the way the equations of the hybrid model are solved is shown explicitly. Several simulation techniques, such as the grid structure and various boundary conditions, are discussed in detail as well. The testing and monitoring of simulation runs is presented as an essential routine when running sophisticated and complex models. Several significant improvements of the model, that are in preparation, are also discussed. A main part of this dissertation are four scientific articles based on the results of the Titan model. The Titan model developed during the course of the Ph.D. research has been shown to be an important tool to understand Titan's plasma interaction. One reason for this is that the structures of the magnetic field around Titan are very much three-dimensional. The simulation results give a general picture of the magnetic fields in the vicinity of Titan. The magnetic fine structure of Titan's wake as seen in the simulations seems connected to Alfvén waves an important wave mode in space plasmas. The particle escape from Titan is also a major part of these studies. Our simulations show a bending or turning of Titan's ionotail that we have shown to be a direct result of the basic principles in plasma physics. Furthermore, the ion flux from the magnetosphere of Saturn into Titan's upper atmosphere has been studied. The modelled ion flux has asymmetries that would likely have a large impact in the heating in different parts of Titan's upper atmosphere.

The watershed simulation and forecasting system in the National Board of and the Environment

Yhteenveto: Vesistömalleihin perustuva vesistöjen seuranta- ja ennustejärjestelmä vesi- ja ympäristöhallinnossa

Components of fractionated stallion seminal plasma and the effects of seminal plasma on sperm longevity

Seminal plasma (SP) is the fluid portion of semen, secreted by the epididymides and the accessory glands before and during ejaculation. The stallion s ejaculate is a series of jets that differ in sperm concentration, semen volume and biochemical composition. Before the actual ejaculation, a clear and watery pre-sperm fluid is secreted. The first three jets form the sperm-rich fractions, and contain ¾ of the total number of sperm. The semen volume and sperm concentration in each of the jets decrease towards the end of the ejaculation, and the last jets are sperm-poor fractions with a low sperm concentration. The aims of these studies were to examine the effects of the different SP fractions, and the presence of SP, on sperm survival during storage. Pre-sperm fluid, and semen fractions with a high (sperm-rich) and low (sperm-poor) sperm concentration were collected in five experiments. The levels of selected enzymes, electrolytes and proteins in different SP fractions were determined. These studies also aimed at assessing the individual variation in the levels of the selected SP components and in the effects of SP on spermatozoa. The association between the components of SP and semen quality, sperm longevity, and fertility was examined with a stepwise linear regression analysis. Compared to samples containing SP during storage, centrifugation and the subsequent removal of SP reduced sperm motility parameters during 24 h of cooled storage in all SP fractions, but sperm membrane integrity was not affected. Some of the measured post-thaw motility parameters were also higher in samples containing SP compared to samples stored without SP. In contrast, the proportion of DNA-damaged spermatozoa was greater in the samples stored with SP than those without SP, and this effect was seen in both sperm-rich and sperm-poor fractions. There were no differences in DNA integrity between fractions stored with SP, but the sperm-rich fraction showed less DNA damage than the sperm-poor fraction after SP removal. The differences between fractions in sperm motility after cooled storage were non-significant. The levels of alkaline phosphatase, acid phosphatase and β-glucuronidase were higher in the sperm-rich fractions compared to the sperm-poor fractions, while the concentrations of calcium and magnesium were higher in sperm-poor fractions than in sperm-rich fractions. The concentrations of sodium and chloride were highest in pre-sperm fluid. In the sperm-poor fraction, the level of potassium was associated with the maintenance of sperm motility during storage. The levels of alkaline and acid phosphatase were associated with sperm concentration and the total number of spermatozoa in the ejaculates. None of the measured SP components were correlated to the first cycle pregnancy rate. In summary, the removal of SP improved DNA integrity after cooled storage compared with samples containing SP. There were no differences in the maintenance of sperm motility between the sperm-rich and sperm-poor fractions and whole ejaculates during cooled storage, irrespective of the presence of SP. The lowest rate of DNA damage was found in the sperm-rich fractions stored without SP. In practice, the results presented in this thesis support the use of individual modifications of semen processing techniques for cooled transported semen from subfertile stallions.

Simulation and graph mining tools for improving gene mapping efficiency

Gene mapping is a systematic search for genes that affect observable characteristics of an organism. In this thesis we offer computational tools to improve the efficiency of (disease) gene-mapping efforts. In the first part of the thesis we propose an efficient simulation procedure for generating realistic genetical data from isolated populations. Simulated data is useful for evaluating hypothesised gene-mapping study designs and computational analysis tools. As an example of such evaluation, we demonstrate how a population-based study design can be a powerful alternative to traditional family-based designs in association-based gene-mapping projects. In the second part of the thesis we consider a prioritisation of a (typically large) set of putative disease-associated genes acquired from an initial gene-mapping analysis. Prioritisation is necessary to be able to focus on the most promising candidates. We show how to harness the current biomedical knowledge for the prioritisation task by integrating various publicly available biological databases into a weighted biological graph. We then demonstrate how to find and evaluate connections between entities, such as genes and diseases, from this unified schema by graph mining techniques. Finally, in the last part of the thesis, we define the concept of reliable subgraph and the corresponding subgraph extraction problem. Reliable subgraphs concisely describe strong and independent connections between two given vertices in a random graph, and hence they are especially useful for visualising such connections. We propose novel algorithms for extracting reliable subgraphs from large random graphs. The efficiency and scalability of the proposed graph mining methods are backed by extensive experiments on real data. While our application focus is in genetics, the concepts and algorithms can be applied to other domains as well. We demonstrate this generality by considering coauthor graphs in addition to biological graphs in the experiments.

Towards Understanding Powder Behavior Via Simulation

The aim of this study was to investigate powder and tablet behavior at the level of mechanical interactions between single particles. Various aspects of powder packing, mixing, compression, and bond formation were examined with the aid of computer simulations. The packing and mixing simulations were based on spring forces interacting between particles. Packing and breakage simulations included systems in which permanent bonds were formed and broken between particles, based on their interaction strengths. During the process, a new simulation environment based on Newtonian mechanics and elementary interactions between the particles was created, and a new method for evaluating mixing was developed. Powder behavior is a complicated process, and many of its aspects are still unclear. Powders as a whole exhibit some aspects of solids and others of liquids. Therefore, their physics is far from clear. However, using relatively simple models based on particle-particle interaction, many powder properties could be replicated during this work. Simulated packing densities were similar to values reported in the literature. The method developed for describing powder mixing correlated well with previous methods. The new method can be applied to determine mixing in completely homogeneous materials, without dividing them into different components. As such, it can describe the efficiency of the mixing method, regardless of the powder's initial setup. The mixing efficiency at different vibrations was examined, and we found that certain combinations of amplitude, direction, and frequencies resulted in better mixing while using less energy. Simulations using exponential force potentials between particles were able to explain the elementary compression behavior of tablets, and create force distributions that were similar to the pressure distributions reported in the literature. Tablet-breaking simulations resulted in breaking strengths that were similar to measured tablet breaking strengths. In general, many aspects of powder behavior can be explained with mechanical interactions at the particle level, and single particle properties can be reliably linked to powder behavior with accurate simulations.

Optimizing simulation model on forest policy.

XVIII IUFRO World Congress, Ljubljana 1986.

Forest sector simulation models as methodological tools in forest policy analysis.

XVIII IUFRO World Congress, Ljubljana 1986.

Application of an ecological simulation model to Lake Päijänne

Tiivistelmä: Simulointimallin soveltaminen Pohjois-Päijänteellä

Effects of gender, and SLCO1B1 and CYP7A1 polymorphisms on plasma bile acids in humans and the pharmacokinetics of therapeutic ursodeoxycholic acid

Bile acids are important steroid-derived molecules essential for fat absorption in the small intestine. They are produced in the liver and secreted into the bile. Bile acids are transported by bile flow to the small intestine, where they aid the digestion of lipids. Most bile acids are reabsorbed in the small intestine and return to the liver through the portal vein. The whole recycling process is referred to as the enterohepatic circulation, during which only a small amount of bile acids are removed from the body via faeces. The enterohepatic circulation of bile acids involves the delicate coordination of a number of bile acid transporters expressed in the liver and the small intestine. Organic anion transporting polypeptide 1B1 (OATP1B1), encoded by the solute carrier organic anion transporter family, member 1B1 (SLCO1B1) gene, mediates the sodium independent hepatocellular uptake of bile acids. Two common SNPs in the SLCO1B1 gene are well known to affect the transport activity of OATP1B1. Moreover, bile acid synthesis is an important elimination route for cholesterol. Cholesterol 7α-hydroxylase (CYP7A1) is the rate-limiting enzyme of bile acid production. The aim of this thesis was to investigate the effects of SLCO1B1 polymorphism on the fasting plasma levels of individual endogenous bile acids and a bile acid synthesis marker, and the pharmacokinetics of exogenously administered ursodeoxycholic acid (UDCA). Furthermore, the effects of CYP7A1 genetic polymorphism and gender on the fasting plasma concentrations of individual endogenous bile acids and the bile acid synthesis marker were evaluated. Firstly, a high performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) method for the determination of bile acids was developed (Study I). A retrospective study examined the effects of SLCO1B1 genetic polymorphism on the fasting plasma concentrations of individual bile acids and a bile acid synthesis marker in 65 healthy subjects (Study II). In another retrospective study with 143 healthy individuals, the effects of CYP7A1 genetic polymorphism and gender as well as SLCO1B1 polymorphism on the fasting plasma levels of individual bile acids and the bile acid synthesis marker were investigated (Study III). The effects of SLCO1B1 polymorphism on the pharmacokinetics of exogenously administered UDCA were evaluated in a prospective genotype panel study including 27 healthy volunteers (Study IV). A robust, sensitive and simple HPLC-MS/MS method was developed for the simultaneous determination of 16 individual bile acids in human plasma. The method validation parameters for all the analytes met the requirements of the FDA (Food and Drug Administration) bioanalytical guidelines. This HPLC-MS/MS method was applied in Studies II-IV. In Study II, the fasting plasma concentrations of several bile acids and the bile acid synthesis marker seemed to be affected by SLCO1B1 genetic polymorphism, but these findings were not replicated in Study III with a larger sample size. Moreover, SLCO1B1 polymorphism had no effect on the pharmacokinetic parameters of exogenously administered UDCA. Furthermore, no consistent association was observed between CYP7A1 genetic polymorphism and the fasting plasma concentrations of individual bile acids or the bile acid synthesis marker. In contrast, gender had a major effect on the fasting plasma concentrations of several bile acids and also total bile acids. In conclusion, gender, but not SLCO1B1 or CYP7A1 polymorphisms, has a major effect on the fasting plasma concentrations of individual bile acids. Moreover, the common genetic polymorphism of CYP7A1 is unlikely to influence the activity of CYP7A1 under normal physiological conditions. OATP1B1 does not play an important role in the in vivo disposition of exogenously administered UDCA.

Modelling vacuum arcs : from plasma initiation to surface interactions

A better understanding of vacuum arcs is desirable in many of today's 'big science' projects including linear colliders, fusion devices, and satellite systems. For the Compact Linear Collider (CLIC) design, radio-frequency (RF) breakdowns occurring in accelerating cavities influence efficiency optimisation and cost reduction issues. Studying vacuum arcs both theoretically as well as experimentally under well-defined and reproducible direct-current (DC) conditions is the first step towards exploring RF breakdowns. In this thesis, we have studied Cu DC vacuum arcs with a combination of experiments, a particle-in-cell (PIC) model of the arc plasma, and molecular dynamics (MD) simulations of the subsequent surface damaging mechanism. We have also developed the 2D Arc-PIC code and the physics model incorporated in it, especially for the purpose of modelling the plasma initiation in vacuum arcs. Assuming the presence of a field emitter at the cathode initially, we have identified the conditions for plasma formation and have studied the transitions from field emission stage to a fully developed arc. The 'footing' of the plasma is the cathode spot that supplies the arc continuously with particles; the high-density core of the plasma is located above this cathode spot. Our results have shown that once an arc plasma is initiated, and as long as energy is available, the arc is self-maintaining due to the plasma sheath that ensures enhanced field emission and sputtering. The plasma model can already give an estimate on how the time-to-breakdown changes with the neutral evaporation rate, which is yet to be determined by atomistic simulations. Due to the non-linearity of the problem, we have also performed a code-to-code comparison. The reproducibility of plasma behaviour and time-to-breakdown with independent codes increased confidence in the results presented here. Our MD simulations identified high-flux, high-energy ion bombardment as a possible mechanism forming the early-stage surface damage in vacuum arcs. In this mechanism, sputtering occurs mostly in clusters, as a consequence of overlapping heat spikes. Different-sized experimental and simulated craters were found to be self-similar with a crater depth-to-width ratio of about 0.23 (sim) - 0.26 (exp). Experiments, which we carried out to investigate the energy dependence of DC breakdown properties, point at an intrinsic connection between DC and RF scaling laws and suggest the possibility of accumulative effects influencing the field enhancement factor.

Literature Review: Investigating Drug Transport at the Blood-Brain Barrier and the Effects of P-glycoprotein on the Brain Distribution of Drugs.

The blood-brain barrier (BBB) is a unique barrier that strictly regulates the entry of endogenous substrates and xenobiotics into the brain. This is due to its tight junctions and the array of transporters and metabolic enzymes that are expressed. The determination of brain concentrations in vivo is difficult, laborious and expensive which means that there is interest in developing predictive tools of brain distribution. Predicting brain concentrations is important even in early drug development to ensure efficacy of central nervous system (CNS) targeted drugs and safety of non-CNS drugs. The literature review covers the most common current in vitro, in vivo and in silico methods of studying transport into the brain, concentrating on transporter effects. The consequences of efflux mediated by p-glycoprotein, the most widely characterized transporter expressed at the BBB, is also discussed. The aim of the experimental study was to build a pharmacokinetic (PK) model to describe p-glycoprotein substrate drug concentrations in the brain using commonly measured in vivo parameters of brain distribution. The possibility of replacing in vivo parameter values with their in vitro counterparts was also studied. All data for the study was taken from the literature. A simple 2-compartment PK model was built using the Stella™ software. Brain concentrations of morphine, loperamide and quinidine were simulated and compared with published studies. Correlation of in vitro measured efflux ratio (ER) from different studies was evaluated in addition to studying correlation between in vitro and in vivo measured ER. A Stella™ model was also constructed to simulate an in vitro transcellular monolayer experiment, to study the sensitivity of measured ER to changes in passive permeability and Michaelis-Menten kinetic parameter values. Interspecies differences in rats and mice were investigated with regards to brain permeability and drug binding in brain tissue. Although the PK brain model was able to capture the concentration-time profiles for all 3 compounds in both brain and plasma and performed fairly well for morphine, for quinidine it underestimated and for loperamide it overestimated brain concentrations. Because the ratio of concentrations in brain and blood is dependent on the ER, it is suggested that the variable values cited for this parameter and its inaccuracy could be one explanation for the failure of predictions. Validation of the model with more compounds is needed to draw further conclusions. In vitro ER showed variable correlation between studies, indicating variability due to experimental factors such as test concentration, but overall differences were small. Good correlation between in vitro and in vivo ER at low concentrations supports the possibility of using of in vitro ER in the PK model. The in vitro simulation illustrated that in the simulation setting, efflux is significant only with low passive permeability, which highlights the fact that the cell model used to measure ER must have low enough paracellular permeability to correctly mimic the in vivo situation.