Interstellar titanium-oxides in interplanetary dust

Interstellar gas abundances (Clayton et al., 1986) suggest that titanium may be bound up in dust and indeed, excess titanium in carbonaceous chondrites is attributed to mixing of interstellar and Solar System materials (Morton, 1974). Fine-grained chondritic interplanetary dust particles (lOPs) of cometary origin are relatively pristine early Solar System materials (Mackinnon and Rietmeijer, 1987; Rietmeijer, 1987) and show chemical and mineralogical signatures related to a pre-solar or nebular origin. For example, large OtH ratios suggest a presolar or interstellar dust component in some chondritic lOPs(Mackinnon and Rietmeijer, 1987). Ti/Si ratios (normalized to bulk CI) in lOPs and carbonaceous chondrite matrices exceed solar abundances but are similar to dust from comet Halley (Jessberger et al., 1987). The Ti-distribution in chondritic lOPs shows major, small-scale « 0.1 urn) variations (Flynn et al., 1978) consistent with heterogeneously distributed Ti-bearingphases. Analytical electron microscope (AEM) studies, in fact, have identified platey grains of Ti-metal, Ti407 and Ti s09 in two different lOPs (Mackinnon and Rietmeijer, 1987). The occurrence of Ti407 was related in situ low-temperature aqueous alteration and therefore implied the presence of BaTi03 (Rietmeijer and Mackinnon, 1984). Yet, the presence ofTis09 in an lOp which shows no evidence of aqueous alteration (Rietmeijer.and McKay, 1986) requires a different interpretation. The distribution of Ti-oxides in chondritic lOPs were investigated with ultra-microtomed thin sections of fluffy chondri tic lOP U2011*B (lSC allocation U2011C2) using a lEOL 2000FX AEM operating at an accelerating voltage of 200kV and with an attached Tracor Northern TN5500 energy dispersive spectrometer.

Grain size distributions of Magneli phases and metallic titanium in chondritic porous interplanetary dust particles

A mineralogical survey of chondritic interplanetary dust particles (IDPs)showed that these micrometeorites differ significantly in form and texture from components of carbonaceous chondrites and contain some mineral assemblages which do not occur in any meteorite class1. Models of chondritic IDP mineral evolution generally ignore the typical (ultra-) fine grain size of consituent minerals which range between 0.002-0.1µm in size2. The chondritic porous (CP) subset of chondritic IDPs is probably debris from short period comets although evidence for a cometary origin is still circumstantial3. If CP IDPs represent dust from regions of the Solar System in which comet accretion occurred, it can be argued that pervasive mineralogical evolution of IDP dust has been arrested due to cryogenic storage in comet nuclei. Thus, preservation in CP IDPs of "unusual meteorite minerals", such as oxides of tin, bismuth and titanium4, should not be dismissed casually. These minerals may contain specific information about processes that occurred in regions of the solar nebula, and early Solar System, which spawned the IDP parent bodies such as comets and C, P and D asteroids6. It is not fully appreciated that the apparent disparity between the mineralogy of CP IDPs and carbonaceous chondrite matrix may also be caused by the choice of electron-beam techniques with different analytical resolution. For example, Mg-Si-Fe distributions of Cl matrix obtained by "defocussed beam" microprobe analyses are displaced towards lower Fe-values when using analytical electron microscope (AEM)data which resolve individual mineral grains of various layer silicates and magnetite in the same matrix6,7. In general, "unusual meteorite minerals" in chondritic IDPs, such as metallic titanium, Tin01-n(Magneli phases) and anatase8 add to the mineral data base of fine-grained Solar System materials and provide constraints on processes that occurred in the early Solar System.

Mineralogy of chondritic porous aggregates : current status

Chondritic porous aggregates (CPA's) belong to an important subset of small particles (usually between 5 and 50 micrometers) collected from the stratosphere by high flying aircraft. These aggregates are approximately chondritic in elemental abundance and are composed of many thousands of small­er, submicrometer particles. CPA particles have been the subject of intensive study during the past few years [1-3] and there is strong evidence that they are a new class of extraterrestrial material not represented in the meteorite collection [3,4]. However, CPA's may be related to carbonaceous chondrites and in fact, both may be part of a continuum of primitive extraterrestrial materials [5]. The importance of CPA's stems from suggestions that they are very primitive solar system material possibly derived from early formed proto­ planets, chondritic parent bodies, or comets [3, 6]. To better understand the origin and evolution of these particles, we have attempted to summarize all of the mineralogical data on identified CPA's published since about 1976.

Calcite precipitation induced by clay-bleach cation exchange in andesitic reservoir rocks

Experiments were carried out on the sodium hypochlorite bleach sensitivity of a deep subsurface andesitic reservoir in order to predict possible deleterious mineral transformations during a downhole clean-up job. Experiments involved examination of core samples from the reservoir using an Environmental Scanning Electron Microscope (ESEM) with an attached Energy Dispersive Spectrometer (EDS) before and after the samples were immersed in bleach. Bleach immersion of whole-rock samples resulted in rapid (less than 1 min) precipitation of abundant 3.0-10.0-μm-wide calcite rhombs within clay-associated micropores and on clay and feldspar grain surfaces. Abundant microporefilling calcite rhombs also formed in pure separates of constituent chlorite/corrensite, whereas no calcite formed in a pure separate of constituent zeolite. These experiments indicate that corrensite is the likely calcium source in this experimental fluid-rock system. Formation of calcite occurs via a cation exchange reaction in which calcium in the smectitic interlayers of corrensite exchanges for sodium in the bleach. Serious formation damage due to calcite precipitation would have occurred in the andesite reservoir had it been exposed to bleach. This finding gives credence to earlier suggestions that cation exchange reactions have the potential to cause calcite precipitation in some sandstone reservoirs when exposed to drilling, completion or stimulation fluids. © 1993.

Metastable carbon in two chondritic porous interplanetary dust particles

An understanding of carbonaceous matter in primitive extraterrestrial materials is an essential component of studies on dust evolution in the interstellar medium and the early history of the Solar System. We have suggested previously that a record of graphitization is preserved in chondritic porous (CP) aggregates and carbonaceous chondrites1,2 and that the detailed mineralogy of CP aggregates can place boundary conditions on the nature of both physical and chemical processes which occurred at the time of their formation2,3. Here, we report further analytical electron microscope (AEM) studies on carbonaceous material in two CP aggregates which suggest that a record of hydrocarbon carbonization may also be preserved in these materials. This suggestion is, based upon the presence of well-ordered carbon-2H (lonsdaleite) in CP aggregates W7029*A and W7010*A2. This carbon is a metastable phase resulting from hydrous pyrolysis below 300-350°C and may be a precursor to poorly graphitized carbons (PGCs) in primitive extraterrestrial materials2. © 1987 Nature Publishing Group.

New phyllosilicate types in a carbonaceous chondrite matrix

CARBONACEOUS chondrites provide valuable information as they are the least altered examples of early Solar System material1. The matrix constitutes a major proportion of carbonaceous chondrites. Despite many past attempts, unambiguous identification of the minerals in the matrix has not been totally successful2. This is mainly due to the extremely fine-grained nature of the matrix phases. Recently, progress in the characterisation of these phases has been made by electron diffraction studies3,4. We present here the direct observation, by high resolution imaging, of phases in carbonaceous chondrite matrices. We used ion-thinned sections from the Murchison C2(M) meteorite for transmission electron microscopy. The Murchison matrix contains both ordered and disordered inter-growths of serpentine-like and brucite-like layers. Such mixed-layer structures are new types of layer silicates. © 1979 Nature Publishing Group.

On gravitational fields of stars in f(R) gravity theories

Einstein's general relativity is a classical theory of gravitation: it is a postulate on the coupling between the four-dimensional, continuos spacetime and the matter fields in the universe, and it yields their dynamical evolution. It is believed that general relativity must be replaced by a quantum theory of gravity at least at extremely high energies of the early universe and at regions of strong curvature of spacetime, cf. black holes. Various attempts to quantize gravity, including conceptually new models such as string theory, have suggested that modification to general relativity might show up even at lower energy scales. On the other hand, also the late time acceleration of the expansion of the universe, known as the dark energy problem, might originate from new gravitational physics. Thus, although there has been no direct experimental evidence contradicting general relativity so far - on the contrary, it has passed a variety of observational tests - it is a question worth asking, why should the effective theory of gravity be of the exact form of general relativity? If general relativity is modified, how do the predictions of the theory change? Furthermore, how far can we go with the changes before we are face with contradictions with the experiments? Along with the changes, could there be new phenomena, which we could measure to find hints of the form of the quantum theory of gravity? This thesis is on a class of modified gravity theories called f(R) models, and in particular on the effects of changing the theory of gravity on stellar solutions. It is discussed how experimental constraints from the measurements in the Solar System restrict the form of f(R) theories. Moreover, it is shown that models, which do not differ from general relativity at the weak field scale of the Solar System, can produce very different predictions for dense stars like neutron stars. Due to the nature of f(R) models, the role of independent connection of the spacetime is emphasized throughout the thesis.

Lightcurve inversion for asteroid spins and shapes

Knowledge of the physical properties of asteroids is crucial in many branches of solar-system research. Knowledge of the spin states and shapes is needed, e.g., for accurate orbit determination and to study the history and evolution of the asteroids. In my thesis, I present new methods for using photometric lightcurves of asteroids in the determination of their spin states and shapes. The convex inversion method makes use of a general polyhedron shape model and provides us at best with an unambiguous spin solution and a convex shape solution that reproduces the main features of the original shape. Deriving information about the non-convex shape features is, in principle, also possible, but usually requires a priori information about the object. Alternatively, a distribution of non-convex solutions, describing the scale of the non-convexities, is also possible to be obtained. Due to insufficient number of absolute observations and inaccurately defined asteroid phase curves, the $c/b$-ratio, i.e., the flatness of the shape model is often somewhat ill-defined. However, especially in the case of elongated objects, the flatness seems to be quite well constrained, even in the case when only relative lightcurves are available. The results prove that it is, contrary to the earlier misbelief, possible to derive shape information from the lightcurve data if a sufficiently wide range of observing geometries is covered by the observations. Along with the more accurate shape models, also the rotational states, i.e., spin vectors and rotation periods, are defined with improved accuracy. The shape solutions obtained so far reveal a population of irregular objects whose most descriptive shape characteristics, however, can be expressed with only a few parameters. Preliminary statistical analyses for the shapes suggests that there are correlations between shape and other physical properties, such as the size, rotation period and taxonomic type of the asteroids. More shape data of, especially, the smallest and largest asteroids, as well as the fast and slow rotators is called for in order to be able to study the statistics more thoroughly.

Multiwavelength studies of regolith effects in planetary remote sensing

A large proportion of our knowledge about the surfaces of atmosphereless solar-system bodies is obtained through remote-sensing measurements. The measurements can be carried out either as ground-based telescopic observations or space-based observations from orbiting spacecraft. In both cases, the measurement geometry normally varies during the observations due to the orbital motion of the target body, the spacecraft, etc.. As a result, the data are acquired over a variety of viewing and illumination angles. Surfaces of planetary bodies are usually covered with a layer of loose, broken-up rock material called the regolith whose physical properties affect the directional dependence of remote-sensed measurements. It is of utmost importance for correct interpretation of the remote-sensed data to understand the processes behind this alteration. In the thesis, the multi-angular effects that the physical properties of the regolith have on remote-sensing measurements are studied in two regimes of electromagnetic radiation, visible to near infrared and soft X-rays. These effects are here termed generally the regolith effects in remote sensing. Although the physical mechanisms that are important in these regions are largely different, notable similarities arise in the methodology that is used in the study of the regolith effects, including the characterization of the regolith both in experimental studies and in numerical simulations. Several novel experimental setups have been constructed for the thesis. Alongside the experimental work, theoretical modelling has been carried out, and results from both approaches are presented. Modelling of the directional behaviour of light scattered from a regolith is utilized to obtain shape and spin-state information of several asteroids from telescopic observations and to assess the surface roughness and single-scattering properties of lunar maria from spacecraft observations. One of the main conclusions is that the azimuthal direction is an important factor in detailed studies of planetary surfaces. In addition, even a single parameter, such as porosity, can alter the light scattering properties of a regolith significantly. Surface roughness of the regolith is found to alter the elemental fluorescence line ratios of a surface obtained through planetary soft X-ray spectrometry. The results presented in the thesis are among the first to report this phenomenon. Regolith effects need to be taken into account in the analysis of remote-sensed data, providing opportunities for retrieving physical parameters of the surface through inverse methods.

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.

On ion escape from Venus

This doctoral thesis is about the solar wind influence on the atmosphere of the planet Venus. A numerical plasma simulation model was developed for the interaction between Venus and the solar wind to study the erosion of charged particles from the Venus upper atmosphere. The developed model is a hybrid simulation where ions are treated as particles and electrons are modelled as a fluid. The simulation was used to study the solar wind induced ion escape from Venus as observed by the European Space Agency's Venus Express and NASA's Pioneer Venus Orbiter spacecraft. Especially, observations made by the ASPERA-4 particle instrument onboard Venus Express were studied. The thesis consists of an introductory part and four peer-reviewed articles published in scientific journals. In the introduction Venus is presented as one of the terrestrial planets in the Solar System and the main findings of the work are discussed within the wider context of planetary physics. Venus is the closest neighbouring planet to the Earth and the most earthlike planet in its size and mass orbiting the Sun. Whereas the atmosphere of the Earth consists mainly of nitrogen and oxygen, Venus has a hot carbon dioxide atmosphere, which is dominated by the greenhouse effect. Venus has all of its water in the atmosphere, which is only a fraction of the Earth's total water supply. Since planets developed presumably in similar conditions in the young Solar System, why Venus and Earth became so different in many respects? One important feature of Venus is that the planet does not have an intrinsic magnetic field. This makes it possible for the solar wind, a continuous stream of charged particles from the Sun, to flow close to Venus and to pick up ions from the planet's upper atmosphere. The strong intrinsic magnetic field of the Earth dominates the terrestrial magnetosphere and deflects the solar wind flow far away from the atmosphere. The region around Venus where the planet's atmosphere interacts with the solar wind is called the plasma environment or the induced magnetosphere. Main findings of the work include new knowledge about the movement of escaping planetary ions in the Venusian induced magnetosphere. Further, the developed simulation model was used to study how the solar wind conditions affect the ion escape from Venus. Especially, the global three-dimensional structure of the Venusian particle and magnetic environment was studied. The results help to interpret spacecraft observations around the planet. Finally, several remaining questions were identified, which could potentially improve our knowledge of the Venus ion escape and guide the future development of planetary plasma simulations.

Embracing Integrability in Stellar and Planetary Dynamics

Hamiltonian systems in stellar and planetary dynamics are typically near integrable. For example, Solar System planets are almost in two-body orbits, and in simulations of the Galaxy, the orbits of stars seem regular. For such systems, sophisticated numerical methods can be developed through integrable approximations. Following this theme, we discuss three distinct problems. We start by considering numerical integration techniques for planetary systems. Perturbation methods (that utilize the integrability of the two-body motion) are preferred over conventional "blind" integration schemes. We introduce perturbation methods formulated with Cartesian variables. In our numerical comparisons, these are superior to their conventional counterparts, but, by definition, lack the energy-preserving properties of symplectic integrators. However, they are exceptionally well suited for relatively short-term integrations in which moderately high positional accuracy is required. The next exercise falls into the category of stability questions in solar systems. Traditionally, the interest has been on the orbital stability of planets, which have been quantified, e.g., by Liapunov exponents. We offer a complementary aspect by considering the protective effect that massive gas giants, like Jupiter, can offer to Earth-like planets inside the habitable zone of a planetary system. Our method produces a single quantity, called the escape rate, which characterizes the system of giant planets. We obtain some interesting results by computing escape rates for the Solar System. Galaxy modelling is our third and final topic. Because of the sheer number of stars (about 10^11 in Milky Way) galaxies are often modelled as smooth potentials hosting distributions of stars. Unfortunately, only a handful of suitable potentials are integrable (harmonic oscillator, isochrone and Stäckel potential). This severely limits the possibilities of finding an integrable approximation for an observed galaxy. A solution to this problem is torus construction; a method for numerically creating a foliation of invariant phase-space tori corresponding to a given target Hamiltonian. Canonically, the invariant tori are constructed by deforming the tori of some existing integrable toy Hamiltonian. Our contribution is to demonstrate how this can be accomplished by using a Stäckel toy Hamiltonian in ellipsoidal coordinates.

Light scattering in random media with wavelength-scale structures : astronomical and industrial applications

Light scattering, or scattering and absorption of electromagnetic waves, is an important tool in all remote-sensing observations. In astronomy, the light scattered or absorbed by a distant object can be the only source of information. In Solar-system studies, the light-scattering methods are employed when interpreting observations of atmosphereless bodies such as asteroids, atmospheres of planets, and cometary or interplanetary dust. Our Earth is constantly monitored from artificial satellites at different wavelengths. With remote sensing of Earth the light-scattering methods are not the only source of information: there is always the possibility to make in situ measurements. The satellite-based remote sensing is, however, superior in the sense of speed and coverage if only the scattered signal can be reliably interpreted. The optical properties of many industrial products play a key role in their quality. Especially for products such as paint and paper, the ability to obscure the background and to reflect light is of utmost importance. High-grade papers are evaluated based on their brightness, opacity, color, and gloss. In product development, there is a need for computer-based simulation methods that could predict the optical properties and, therefore, could be used in optimizing the quality while reducing the material costs. With paper, for instance, pilot experiments with an actual paper machine can be very time- and resource-consuming. The light-scattering methods presented in this thesis solve rigorously the interaction of light and material with wavelength-scale structures. These methods are computationally demanding, thus the speed and accuracy of the methods play a key role. Different implementations of the discrete-dipole approximation are compared in the thesis and the results provide practical guidelines in choosing a suitable code. In addition, a novel method is presented for the numerical computations of orientation-averaged light-scattering properties of a particle, and the method is compared against existing techniques. Simulation of light scattering for various targets and the possible problems arising from the finite size of the model target are discussed in the thesis. Scattering by single particles and small clusters is considered, as well as scattering in particulate media, and scattering in continuous media with porosity or surface roughness. Various techniques for modeling the scattering media are presented and the results are applied to optimizing the structure of paper. However, the same methods can be applied in light-scattering studies of Solar-system regoliths or cometary dust, or in any remote-sensing problem involving light scattering in random media with wavelength-scale structures.

A parity-violatingU4-gravitation theory

The general structure of a metric-torsion theory of gravitation allows a parity-violating contribution to the complete action which is linear in the curvature tensor and vanishes identically in the absence of torsion. The resulting action involves, apart from the constant ¯K E =8pgr/c4, a coupling (B) which governs the strength of the parity interaction mediated by torsion. In this model the Brans-Dicke scalar field generates the torsion field, even though it has zero spin. The interesting consequence of the theory is that its results for the solar-system differ very little from those obtained from Brans-Dicke (BD) theory. Therefore the theory is indistinguishable from BD theory in solar-system experiments.

A new theory of gravitation with parity-violating propagating torsion

A new theory of gravitation has been proposed in a more general space-time than Riemannian. It is a generalization of the ECSK and Brans-Dicke (BD) theory of gravitation. It is found that, in contrast to the standard the ECSK theory, a parity-violating propagating torsion is generated by the BD scalar field. The interesting consequence of the theory is that it can successfully predict solar system experimental results to desired accuracy.