Condensation and aggregation of solar corundum and corundum‐hibonite grains

Abstract— Forty‐three corundum grains (1–11 μm in size) and 5 corundum‐hibonite grains with corundum overgrown by hibonite (4–7 μm in size), were found in the matrix of the mineralogically pristine, ungrouped carbonaceous chondrite Acfer 094 by using cathodoluminescence imaging. Some of the corundum and corundum‐hibonite grains occur as aggregates of 2 to 6 grains having similar sizes. The oxygen isotopic compositions of some of the corundum‐bearing grains suggest their solar nebula origin. ²⁶Al‐²⁶Mg systematics of one corundum grain showed the canonical initial ²⁶Al/²⁷Al ratio, also suggesting a solar nebula origin. Quantitative evaluation of condensation and accretion processes made based on the homogeneous nucleation of corundum, diffusion‐controlled hibonite formation, collisions of grains in the nebula, and critical velocity for sticking, indicates that, in contrast to the hibonite‐bearing aggregates of corundum grains, the hibonite‐free corundum aggregates could not have formed in the slowly cooling nebular region with solar composition. We suggest instead that such aggregates formed near the protosun, either in a region that stayed above the condensation temperature of hibonite for a long time or in a chemically fractionated, Ca‐depleted region, and were subsequently physically removed from this hot region, e.g., by disk wind.     

Recipes for making synthetic CAIs,refractory residues,and minerals for rim‐forming experiments

Abstract— Recipes are presented for synthesizing various type A and type B Ca‐Al‐rich inclusions (CAIs), refractory volatilization residues, and the minerals forsterite and melilite that are required for experiments. These experiments (described in other works) aim to make two determinations: 1) the conditions under which the surfaces of CAIs were either “flash‐heated” or “volatilized subsolidus” to form a temporary ultra‐refractory residue, and 2) the conditions under which the residue was then metasomatized to form the mineral layers making up Wark‐Lovering (WL) rims on CAIs.     

The stability of hibonite,melilite and other aluminous phases in silicate melts: Implications for the origin of hibonite-bearing inclusions from carbonaceous chondrites

Abstract— Phase fields in which hibonite and silicate melt coexist with spinel, CaAl₄O₇, gehlenitic melilite, anorthite or corundum at 1 bar in the system CaO-MgO-Al₂O₃-SiO₂-TiO₂ were determined. The hibonites contain up to 1.7 wt% SiO₂. For TiO₂, the experimentally determined partition coefficients between hibonite and coexisting melt, D^Hib/L_i, vary from 0.8 to 2.1 and generally decrease with increasing TiO₂ in the liquid. Based on Ti partitioning between hibonite and melt, bulk inclusion compositions and hibonite-saturated liquidus phase diagrams, the hibonite in hibonite-poor fluffy Type A inclusions from Allende and at least some hibonite from hibonite-rich inclusions is relict, although much of the hibonite from hibonite-glass spherules probably crystallized metastably from a melt Bulk compositions for all of these CAIs are consistent with an origin as melilite + hibonite + spinel + perovskite phase assemblages that were partially altered and in some cases partially or completely melted The duration of the melting event was sufficient to remove any Na introduced by the alteration process but frequently insufficient to dissolve all of the original hibonite. Simple thermochemical models developed for meteoritic melilite and hibonite solid solutions were used to obtain equilibration temperatures of hibonite-bearing phase assemblages with vapor. Referenced to 10^?3 atm, hibonite + corundum + vapor equilibrated at ～1260 °C and hibonite + spinel ± melilite + vapor at 1215 ± 10 °C. If these temperatures reflect condensation in a cooling gas of solar composition, then hibonite ± corundum condensed first, followed by spinel and then melilite. The position of perovskite within this sequence is uncertain, but it probably began to condense before spinel. This sequence of phase appearances and relative temperatures is generally consistent with observed textures but differs from expectations based on classical condensation calculations in that equilibration temperatures are generally lower than predicted and melilite initially condenses with or even after spinel. Simple thermochemical models for the substitution of trace elements into the Ca site of meteoritic hibonites suggest that virtually all Eu is divalent in early condensate hibonites but that Eu²⁺/Eu³⁺ decreases by a factor of 20 or more during the course of condensation primarily because the ratio is proportional to the partial pressure of Al, which decreases dramatically as aluminous phases condense. The relative sizes of Eu and Yb anomalies in meteoritic hibonites and inclusions may be partly due to this effect     

Observation of chromospheric evaporation during the solar maximum mission

A sample of flares detected in 1980 with the Bent Crystal Spectrometer and the Hard X-Ray Burst Spectrometer on the Solar Maximum Mission satellite has been analysed to study the upward motions of part of the soft X-ray emitting plasma. These motions are inferred from the presence of secondary blue-shifted lines in the Ca XIX and Fe XXV spectral regions during the impulsive phase of disk flares. Limb flares do not show such blue-shifted lines indicating that the direction of the plasma motion is mainly radial and outward. The temporal association of these upward motions with the rise of the thermal phase and with the impulsive hard X-ray burst, as well as considerations of the plasma energetics, favour the interpretation of this phenomenon in terms of chromospheric evaporation. The two measureable parameters of the evaporating plasma, emission measure and velocity, depend on parameters related to the energy deposition and to the thermal phase. The evaporation velocity is found to be correlated with the spectral index of the hard X-ray flux and with the rise time of the thermal emission measure of the coronal plasma. The emission measure of the rising plasma is found to be correlated with the total energy deposited by the fast electrons in the chromosphere by collisions during the impulsive phase and with the maximum emission measure of the coronal plasma.     

On the hawking evaporation of Dirac particles in Kerr-Newman spacetime

We have obtained a solution of the Dirac equation near the horizon of a Kerr-Newman black hole and compared it with the solution of the Klein-Gordon equation. We first generalized the work of LIU Liao and XU Dian-yan for a quasi-extreme Kerr black hole to a quasi-extreme Kerr-Newman black hole, then further generalized to a general Kerr-Newman black hole, and thereby verified that the Dirac particles emitted by a general Kerr-Newman black hole do have a black-body energy spectrum.     

Substances of cometary grains estimated from evaporation and radiation pressure mechanisms

The shape and intensity distribution of tails for several large comets are estimated on the basis of grain properties in the solar radiation field. The following results are obtained: (1) The ratio of the maximum radiation pressure force to the gravitaional force acting on dust grains in cometary tails is found to be less than 2.5. This means that grains such as graphite particles in the size range 0.02–0.2 μm do not exist in them, because such particles would allow forces greater than 2.5 (2) Tail substances supplied near the time of perihelion passage for the Sun-grazing comet Ikeya-Seki (1965 VIII) and Comet Seki-Lines (1962 III) were composed of particular grains which had values of radiation pressure ratio less than 1.0. Therefore, it is concluded that the material was composed of silicate grains only, since iron grains had sublimated and there were no graphite particles.     

A simple model for the evaporation of black holes at final stages

We present a simple model for the evaporation of primordial black holes at final stages with the formation of a relic remnant with a mass of 1–10³ m _P1. The model takes into account the conservation of energy and the impossibility of passing through the state with the minimum possible mass. These relic remnants may account for a substantial fraction of dark matter in the Universe.     

Exploration of quantitative kinetic models for the evaporation of silicate melts in vacuum and in hydrogen

Abstract— Two basic approaches (pure component reference (PCR) and equilibrium reference (EQR)) to modeling silicate melt evaporation are explored. The PCR model calculates the maximum possible evaporation rates of the pure oxides from their equilibrium vapor pressures and rescales these rates according to the activities of the oxides in the silicate melts and the melt densities. The EQR model calculates the maximum possible evaporation rates based on the equilibrium vapor pressures of the melts. Differences between the calculated and experimentally determined evaporation rates are accounted for with evaporation (α^evap) coefficients that are only dependent on temperature. Two versions of the PCR model, Cases 1 and 2, are explored to try to resolve apparently contradictory conclusions about the composition of the evaporating species based on Mg and Si isotope fractionation during evaporation (species are not in thermodynamic equilibrium proportions) and direct measurements of gas species in Langmuir experiments (species are in roughly equilibrium proportions). The Case 2 and EQR models cannot explain the observed isotope fractionations unless evaporation occurred under non‐Rayleigh conditions, either because there was significant recondensation during the experiments or because diffusion was playing a limiting role. Whether or not the role of diffusion is included, the PCR and EQR models are able to reproduce the elemental results of evaporation experiments of “chondritic” melts from temperatures of 1700 to 2000 °C, and up to mass losses of about 95%. However, the models underestimate absolute evaporation rates in very Ca‐ and Al‐rich melts. This may reflect errors in the model used to estimate oxide activities. The EQR model can only reproduce the observed evaporation behavior of Na if, unlike the other oxides, its α^evap coefficient is close to unity. Based on available diffusion data, diffusion is not slow enough in “chondritic” or forsteritic melts to explain the isotopic fractionations of Mg and O in the evaporation experiments, but it may play a role in limiting Si isotope fractionation. Provided recondensation was not a significant factor in the experiments, at present PCR Case 1 appears to be the best model if both the Langmuir and the isotopic fractionation experiments are to be explained.     

Iron isotopes in chondrules: Implications for the role of evaporation during chondrule formation

Abstract— We have measured the δ⁵⁷Fe of olivines in nine Chainpur chondrules. All are within error of normal (typically 2σ ≤ 1–2%0). Most of the chondules could not have lost more than ～20% of their FeO by Rayleigh evaporation and none can have lost more than ～61%. Yet, the range of Fo contents in these chondrules is Fo_78–99.9. The isotopic compositions of the chondrules clearly demonstrate that, for instance, type I chondrules cannot form from type II chondrules by evaporation of FeO under Rayleigh conditions. The isotopic compositions also place constraints on the minimum cooling rates these chondrules could have experienced. These cooling rates must also be equal to or slower than those required to produce the chondrule textures. Assuming flash heating and evaporation rates like those measured in vacuum, the minimum cooling rates necessary to prevent detectable Fe isotopic fractionation via Rayleigh evaporation approach those needed to produce barred and porphyritic textures. The presence of hydrogen in the nebula, non‐linear cooling and other effects will all tend to increase the cooling rates required to prevent δ⁵⁷Fe > 1–2%0, perhaps by as much as 1–2 orders of magnitude. The two most likely ways that the cooling rates required to prevent δ⁵⁷Fe >1–2%0 can be kept below those needed to produce barred and porphyritic textures are (1) the pH₂ in the nebula was low enough to keep evaporation rates close to those in vacuum, or (2) back reaction of chondrules with Fe in the gas suppressed isotopic fractionation.     

The flux of secondary anti-deuterons and antihelium produced in the interstellar medium

Several measurements have been performed to find antiprotons in the primary cosmic radiation. Because it is difficult to get completely separated secondary produced antiprotons primary ones, calculations based on accelerator results have been performed for the flux of secondary produced antideuterons and antihelium.     

Anisotropy of Mg isotopic fractionation during evaporation and Mg self-diffusion of forsterite in vacuum

Evaporation of solid materials under low-pressure conditions could play important roles in chemical and isotopic fractionations in the early solar system. We have studied anisotropy of isotopic fractionation of ²⁶Mg and ²⁵Mg during kinetic evaporation of forsterite (Mg₂SiO₄), which is potentially a powerful tool to understand thermal histories of crystals in the early solar system. Ion-microprobe depth profiling revealed that the Mg isotopic zoning profiles of forsterite evaporated at 1500-1700 °C are notably differing along the a-, b-, and c-axes, which can be attributed to anisotropy in self-diffusion coefficient of Mg (D) and an isotopic fractionation factor for evaporation of Mg (α). The D and α were obtained from zoning profiles by applying the diffusion-controlled isotopic fractionation model of Wang et al. [1999. Evaporation of single crystal forsterite: Evaporation kinetics, magnesium isotope fractionation, and implications of mass-dependent isotopic fractionation of a diffusion-controlled reservoir. Geochim. Cosmochim. Acta 63(6), 953-966.].The D is largest and smallest along the a- and c-axes, respectively. The activation energy of 560-670 kJ/mol indicates that Mg diffusion at 1500-1700 °C occurred in the intrinsic diffusion regime.The α seems to be larger along the a- or c-axes than along the b-axis. The α along the a- or c-axes show weak temperature dependence. The α along all the crystallographic orientations is closer to unity than that expected from the kinetic theory of gases. These lines of evidence suggest that surface processes such as breaking of bonds and surface diffusion are responsible for the isotopic fractionation.     

Stability of artificial and natural satellites

A quantitative measure of stability based on Hill's definition is evaluated for direct and retrograde satellite orbits. These orbits are known as Poincaré's first kind in the restricted problem of three bodies. Onsets of possible instabilities and captures are established. A critical (maximum) value of the satellite's orbital radius is found for stability as a remarkably simple function of the massparameter. The results are applied to the natural satellites of the solar system.     

Determination of the electromagnetic field produced by a magnetic oblique-rotator

Approximate solutions for the electromagnetic fields produced by a uniformly magnetized oblique rotator in vacuum are derived systematically in various regions by using general relativistic considerations. The results, which are expressed in compact vector notation, are compared with the component expression of Deutsch (1955) and Backus (1956). In our method, the potentials (, A) and fields (E, B) due to the rotating magnet are represented in terms of the vector potential for a rest magnet, with the proper correction for rotation. In doing this, the transformation laws for various quantities between the rest and rotating frames play important roles. The meaning of rotating magnetic field lines is also considered in connection with the apparent paradox in the gedanken-experiment of a unipolar inductor.     

Determination of the electromagnetic field produced by a magnetic oblique-rotator

The inertial effect on the structure of the magnetosphere of a rotating star is investigated, in the corotation approximation for a surrounding quasi-neutral plasma. The equation of motion reduces to a usual static balance equation between the electromagnetic and the centrifugal forces, in the rotating frame. However the MHD condition, which can be regarded as a special form of the generalized Ohm's law, is modified by the inclusion of inertial effect, with a violation of the frozen-in condition in case of a general (i.e., not restricted to corotation) plasma motion. The inertial effect on the electromagnetic field is summarized in a partial scalar potential named the non-Backus potential, which is proportional to the centrifugal potential in the corotation approximation.An approximate solution of this corotation problem is given, in which another characteristic radiusr _M appears besides the light radiusr _L . This radius defines a distance beyond which the inertial effect becomes dominant over the electromagnetic one, and is useful in estimating the magnitude of the terminal velocity of a centrifugal wind. A few examples of the modification of dipole magnetic field due to the inertial effect are visualized. In an oblique-rotation case, it can be seen that such a warp of the neutral sheet (the surface ofB _r =0) is reproduced as observed in the Jovian magnetosphere.     

Determination of the electromagnetic field produced by a magnetic oblique-rotator

The structure of the corotating region, which forms an inner portion of a stellar magnetosphere, is reconsidered in a quasi-neutral case by taking into account the inertial effects of electrons as well as that of ions up to the first order in their mass ratio (δ=m_?/m₊). It is emphasized first that the magnetosphere is not globally equipotential even in the frame rotating with a central star (i.e. ?#0, where ? is the ‘non-Backus’ potential) due at least to the inertial effects of plasma particles. However, it is shown that the condition ?=0 is asymptotically recovered in the corotating region owing to the presence of the drift current which can be taken into account only when δ is not entirely neglected. This fact suggests that the deviation of the plasma motion in the outer magnetosphere from the corotation can be attributed to the non-zero ?. A globally self-consistent solution is obtained under this condition (?=0). In contrast with the solutions in the ‘force-free’ and the ‘mass-less-electron’ approximations, this solution has a disk structure in the corotation zone in which the plasma and the current density are concentrated to a thin disk near the magnetic equator. Owing to this sheet current in the disk the lines of force of the stellar magnetic field are modified to form a very elongated shape (the magnetodisk) if the plasma β-value is fairly large. Such a disk structure seems to be a common feature in the high β inner magnetospheres of various types of stars.     

Determination of the electromagnetic field produced by a magnetic oblique-rotator

The electromagnetic field produced by a magnetic dipole moment, , which is rotating obliquely surrounded by a corotating plasma sphere, is investigated. This corotating-plasma approximation has the same order of accuracy as the force-free one but has somewhat different physical implications. In the former the effect of non-electromagnetic forces such as the inertial force are included, though in somewhat artificial manner, as a departure from the strict MHD condition and this fact seems to guarantee the existence of physical solutions.Analogous to the relativistic force-free equation, a set of two differential equations (the corotation equation) are derived for the scalar functions associated with the electric and magnetic fields. A self-consistent solution of these equations is given and it is shown that this solution has no singularity, in spite of apparent divergence in the formal solution, on the light cylinder. It is concluded from this solution that, even in the extreme case of the largest possible corotation-radius (i.e.b=r _L , wherer _L is the light radius), the existence of a corotating plasma does not alter the field structure drastically from the vacuum case. It is also suggested through this treatment that inclusion of the inertial term in generalized Ohm's law might be essential in considering the centrifugal-wind problem.     

The cooling of flare produced plasmas in the solar corona

Solar flare X-rays, at energies less than 10 keV, are emitted by hot plasmas located in the corona. Three plasma cooling models are examined in detail. The cooling of the electrons by Coulomb collisions with ions at a lower temperature would require the observed material to occupy very large volumes. Cooling could take place by conduction or by radiation and observations are proposed which would allow the dominant cooling mechanism to be established.On leave during a portion of this work as University Research Fellow in Astronomy, University of Leicester, England.