Modelling the neutral gas environment of comets with special application to P/Halley

A technique has been developed which allows relatively accurate modelling of cometary gas production from nothing more than a visible light curve. Application to P/Halley suggests the production rate of parent molecules will be about 2.6 × 10²⁹ per second on March 10, 1986, for example. The uncertainties and intrinsic limitations in this approach are outlined. The theory is then extended to predictions of abundance of other gaseous species, and a photometric model of these gases provided. Combined with the dust model of N. Divine, preliminary predictions of the luminance of P/Halley as seen in any direction from inside the coma or outside can be provided for λλ3000–7000.     

Numerical simulations of the radiance from the comet 46P/Wirtanen in the various configurations of the measurements during “Rosetta” mission

The work we present deals with the spectrometric measurements of VIRTIS instrument of the Comet P/Wirtanen planned for the Rosetta mission. This spectrometer can monitor (VIRTIS M channel: 0.250μm – 0.980μm; Δκ=20cm⁻¹; 0.980 – 5.0 μm; Δκ = 5cm⁻¹; VIRTIS H channel: 2.0 μm – 5.0 μm; Δκ=2cm⁻¹) the nucleus and the coma in order to provide a general picture of coma's composition, the production of gas and dust, the relationship of coma production to surface composition and the structure and variation of mineralogy of the nucleus surface. During the mission the observation conditions of the spectroscopic investigation change due to different relative positions spacecraft/comet, and to the different illumination conditions of the surface at various distances of the comet to the Sun. The nucleus surface is continuously modified by the ice sublimation accompanied by gas and dust emission. Consequently the surface also its spectrophotometric properties changes and their monitoring can give a new insight. The important role of simulations is to predict the results of measurements in various experimental condition what, in the future, can help in interpretation of the measured data.

In this paper the first results of our simulation the radiance from the comet in the 0.25–5.0μm spectral range at two distances from the Sun (1AU and 3AU) are shown. The distance between the Rosetta orbiter and the nucleus surface as well as the sun zenith angles are taken into account according to the Rosetta mission phases. In fact the surface and coma properties vary along the comet orbit, and should be taken into account in our calculations. The optical parameters of the dust on the surface (e.g. reflectance) and in the coma (e.g. Q_ext) were calculated from optical constants of possible comet analogues. The thermodynamic parameters of the comet are taken from the models of comet evolution. Through this kind of modelling it is possible to identify the surface characteristics in spectra of the radiation from the surface of nucleus transmitted through the coma loaded with dust and gases.

Even if the “Rosetta mission” is postponed, with the consequence of a target change, we think that our idea and the method used for the simulations can be useful also for the new Rosetta target - the comet 67P/Churyumov Gerasimenko.     

Interplanetary scintillations of PKS 2314+03 during occultation by comet Halley

Strong interplanetary scintillations (IPS) of the quasar 2314+03 were recorded at 103 MHz at Thaltej-Ahmedabad, India with a transit type correlation interferometer on 18, 19 and 20 December 1985, as the radio source was predicted to be occulted by the ion tail of the comet Halley.

On 18th through 20th very strong scintillations, with periodicities of 1 sec average were observed, their amplitude progressively decreasing as the source approached the tail-end. The rms variations of scintillating flux of the source on 18, 19 & 20 were about 18, 11 & 4.7 Jy, as against 3.3 Jy on control days 17 and 21 December for solar elongation of 87°.

Assuming Gaussian irregularities with weak scattering, the rms density variations, ΔN, of 10, 6, 3 and 1 elec./cm³ on 18 through 21 December, from the comet nucleus towards its tail-end, varied as (ΔN) ∝ r^−3.3 as against (ΔN) ∝ r⁻² in the solar plasma.

Quasi-periodic modulations of the enhanced scintillating flux possibly imply 10⁴ km scale-size ion condensations and mean electron density of 10³ to 10⁴ electrons/cm³ in the Halley's plasma tail.     

Cometary atmosphere (gas and dust)

There is important progress now in the identifications and measurements of primary (parent) molecules in the inner coma of Comet Halley. H₂O, CO₂ and CO are definitely in the list, CH and some complicate organic molecules are suspected. Gas production rate for water vapor is Q_H2O 10³⁰ s⁻¹. The bulk of data doesn't contradict to the Whipple model of nucleus (with clathrate modification). Pronounced spatial structure of gaseous flow in the coma was observed, but in general measured properties of neutral gas in the coma of Comet Halley are not very different from predicted. Situation for dust is different. In situ dust measurements show that size spectrum and optical properties of particles in coma are substantively declining from predicted on the base of groundbased photometry. However there are discrepancies between Vega and Giotto dust counter data. Dust in the inner coma didn't prevent the succesful imaging of nucleus by TV on Vega 1 and 2.     

Near infrared observations of comet Halley from Vega 2

Spatial distribution of the continuum radiation in the range of 0.95–1.9 μm presumes total dust production rate of the comet of 10ρ tonne s⁻¹ (ρ is the dust material density) and its angular distribution proportional cos . Observations of the water vapor band at 1.38 μ m reveal strong jets, their time shift from the dust jet measured in situ is consistent with gas velocity of 0.82±0.1 km s⁻¹ and dust velocity of 0.55±0.08 km s⁻¹. The OH vibrational-rotational bands observed are excided directly via photolysis of water vapor. Water vapor production rate deduced from the H₂O band and OH band intensities is 8×10²⁹ s⁻¹. Intensity of the CN(0,0) band result in the CN column density of 9×10¹² cm⁻², i.e. larger by a factor of 3 than given by the violet band.     

The GPS measured SITEC caused by the very intense solar flare on July 14, 2000

This work studies the sudden increases in total electron content of the ionosphere caused by the very intense solar flare on July 14, 2000. Total electron content (TEC) data observed from a Global Positioning System (GPS) network are used to calculate the flare-induced TEC increment, δTEC_f, and variation rate, dTEC_f/dt. It is found that both dTEC_f/dt and δTEC_f are closely related with the solar zenith angles. To explain the observation results, we derived a simple relationship between the partial derivative of the flare-induced TEC, ∂TEC_f/∂t, which is a good approximation for dTEC_f/dt, and the solar zenith angle χ, as well as the effective flare radiation flux I_f, according to the well-known Chapman theory of ionization. The derived formula predicted that ∂TEC_f/∂t is proportional to I_f and inverse proportional to Chapman function ch(χ). This theoretical prediction not only explains the correlation of dTEC_f/dt and δTEC_f with χ as shown in our TEC observation, but also gives a way to deduce I_f from TEC observation of GPS network. Thus, the present work shows that GPS observation is a powerful tool in the observation and investigation of solar flare effects on the ionosphere, i.e., the sudden ionospheric disturbances, which is a significant phenomenon of space weather.     

Modeling of the very high velocity impact process with respect to in-situ ionization measurements

This overview deals with very high impact velocities, where complete vaporization of an impacting cosmic dust particle is to be expected upon expansion from the high pressure high temperature state behind the stopping shock (v > 15 km/s). The topics discussed are the mechanics and thermodynamics of compression, adiabatic release, equation of state and nonequilibrium states upon expansion. The case of very high particle porosity (ρ 1 g/cm³) and the case of very small dust masses (m < 10⁻¹⁷ g) are discussed from what one presently knows. The possibility of three body collisions in the expanding gas phase is discussed briefly. The effect of oblique impact is discussed with respect to its relevance to the ionization process. The numbers communicated are up to the highest “experimental” impact velocities (80 km/s, Halley mission). As one goes to lower impact velocities (20 < v < 30 km/s) there is still complete vaporization of the dust particle but ionization out of the bulk of the particle becomes low. Other than thermal processes may become important. Ideas are outlined to understand their physical nature.     

A Far Infrared Photometer (FIRP) for the Infrared Telescope in Space (IRTS)

We describe the design and calibration of the Far-Infrared Photometer (FIRP), one of four focal plane instruments on the Infrared Telescope in Space (IRTS). The FIRP will provide absolute photometry in four bands centered at 150, 250, 400, and 700 μm with spectral resolution λ/Δλ ≈ 3 and spatial resolution ΔΘ = 0.5 degrees. High sensitivity is achieved by using bolometric detectors operated at 300 mK in an AC bridge circuit. The closed-cycle ³He refrigerator can be recycled in orbit. A 2 K shutter provides a zero reference for each field of view. More than 10% of the sky will be surveyed during the ≈3 week mission lifetime with a sensitivity of <10⁻¹³ W·cm⁻²·sr⁻¹ per 0.5 degree pixel.     

VEGA three-channel spectrometer experiment: Spectroscopy of comet Halley in the visible and near-ultraviolet ranges

A total of 3600 spectra of Comet Halley in the 275–710 nm were obtained on March, 8, 9, 10 and 11, 1986, from the VEGA 2 spacecraft. The emissions of OH, NH, CN, C₃, CH, C₂, NH₂ and H₂O⁺ are identified. From the OH intensity in the (0,0) band: 1.1 Megarayleigh at 5400 km from the nucleus, it can be inferred that the OH production rate was (1.4 ± 0.5)×10³⁰ molecules s⁻¹. The NH, C₃, CH and NH₂ bands became comparatively more intense at distances from the nucleus shorter than 3000km. At 06:40 U.T. when the instrument field of view was 6000×4500 km, two jets were observed. Spectra from the jets show significant differences with other spectra. Inside a jet NH, C₃ and NH₂ are comparatively more intense and the rotational distributions of OH, CN and C₂ are strongly distorted. This shows that part of the observed emissions probably comes from radicals directly produced in the excited state during the initial process of photolysis of the parent molecules.     

Cometary probe of the Venera-Halley mission

Venera-Halley mission is to be launched to Venus in Dec. 1984. It will fly by Venus in June 1985. Separation of the cometary probe and Venera descend module will take place at that time. The gravitational swing-by at Venus will provide the encounter with the Halley comet in March 1986. The remote sensing of the inner coma (TV-imagery, spectrometry in the region from 1200 A to 12 μm, polarimetry) and of the nucleus, direct measurements of dust fluxes, dust composition, plasma and magnetic field are planned in the framework of multinational cooperation.     

Ion formation by high- and medium-velocities dust impacts from laboratory measurements and Halley results

The ion formation processes by dust impacts have been studied qualitatively as well as quantitatively by dust accelerator laboratory measurements. Iron, carbon and metallized glass particles in the femto- to nano-gram mass range had been impacted on various metal targets in a velocity regime of v = 2 - 64 km/s. In the high velocity regime as relevant for the (retrograde) Halley encounter more than 99% of the ions produced are singly charged atomic, the rest molecular ones. The ion/atom ratios are apparently modified SIMS yields, the modification parameter being impact velocity dependent. A semiempirical formula was deduced for the determination of mass and density of the impacting particle from target and projectile ion yields. When evaluating the Halley encounter results, the elemental distribution of p/Halley dust appeared nearly to be solar; the organic fraction (CHON) could be characterized in a rough manner as fairly unsaturated. Oligomers of the monomers C₂H₂ (65%), CH₂O (25%), and HCN (10%) are probable.

With medium velocities (for prograde comet encounter), i.e. v = 15-30 km/s molecular ion types govern the mass spectra. Consequently, more chemical information of the projectile can be expected in this case, additional to the elemental distribution. Mass and density of the impinging dust particles can be determined as well.     

Physical properties of cometary dust: Relation of DIDSY data to grain properties

The Dust Impact Detection System (DIDSY) for the Giotto Halley Mission consists of two types of sensors for the detection of cometary dust particles: two impact plasma sensors and five piezo-electric momentum sensors. One sensor of each type is covered by a penetration film. A 1 μm thick aluminum film covers an impact plasma sensor. One momentum sensor is mounted onto the rear shield behind the 1 mm front shield made from aluminum. The parameters measured are the total charge released upon impact and the amplitude of the acoustic signal generated by the impact. Both quantities depend on the mass and speed of the impacting particles. At the impact speed of 68 km/sec the mass of cometary dust particles can be determined in the mass range from 10⁻¹⁷ g to 10⁻³ g. From the difference in the countrates measured by the sensors with and without penetration film the average bulk density of dust particles of masses 10⁻¹⁴ g and 10⁻⁶ g can be determined. With appropriate calibration an accuracy of a factor of 2 for both the mass and density determination can be obtained.     

Effect of track structure and radioprotectors on the induction of oncogenic transformation in murine fibroblasts by heavy ions

The oncogenic potential of high-energy ⁵⁶Fe particles (1 GeV/nucleon) accelerated with the Alternating Gradient Synchrotron at the Brookhaven National Laboratory was examined utilizing the mouse C3H 10T1/2 cell model. The dose-averaged LET for high-energy ⁵⁶Fe is estimated to be 143 keV/μm with the exposure conditions used in this study. For ⁵⁶Fe ions, the maximum relative biological effectiveness (RBE_max) values for cell survival and oncogenic transformation were 7.71 and 16.5 respectively. Compared to 150 keV/μm ⁴He nuclei, high-energy ⁵⁶Fe nuclei were significantly less effective in cell killing and oncogenic induction. The prostaglandin E₁ analog misoprostol, an effective oncoprotector of C3H 10T1/2 cells exposed to X rays, was evaluated for its potential as a radioprotector of oncogenic transformation with high-energy ⁵⁶Fe. Exposure of cells to misoprostol did not alter ⁵⁶Fe cytotoxicity or the rate of ⁵⁶Fe-induced oncogenic transformation.     

Fe XVII X-ray lines in solar coronal plasmas

Calculated intensities of the Fe X-ray lines due to transitions 2p⁶ − 2p⁵3d lines (near 15 Å) and 2p⁶ − 2p⁵3s lines (near 17 Å) are compared with measured line intensities in solar and tokamak spectra. For the solar spectra, temperature T_e is obtained from the ratio of the Fe 16.776 Å line to a nearby Fe line. We find excellent agreement for all the major Fe line features in the 15–17 Å region except the Fe 15.015 Å line, the observed flux of which is less than the theoretical by a factor f. We find that f strongly depends on the heliocentric angle θ of the emitting region, being smallest (0.2) when the region is nearest Sun centre, but nearly 1 near the limb. Attributing this to resonance scattering, we are able to deduce the path length and electron density from the observations. Possible application to stellar active regions is given.     

Submillimeter spectroscopy from space — Missions planned between 1995 and 2005

The submillimeter wavelength range (100 μm ≤ λ ≤ 1 mm) is exceptionally rich in atomic and molecular transitions that provide direct measures of the chemical composition, temperature, density, and cooling of the interstellar medium. In part, this is due to the fact that the most abundant oxygen- and carbon-bearing species as well as the most common hydrides have their lowest lying transitions within this wavelength region. However, because of the great abundance of many of these same species in our own atmosphere (e.g., H₂O and O₂), terrestrial absorption at submillimeter wavelengths is strong, rendering much of the sky inaccessible to ground- and even airborne observing platforms. In recognition of the scientific importance of the submillimeter range to astronomy and the limitations placed on its pursuit by the atmosphere, both NASA and ESA have embarked on ambitious programs to explore these wavelengths from space.     

Analysis of the nucleus and circumnuclear area of comet Halley with the “I.K.S.” infrared sounder from the “Vega” flyby probes

The “Vega” Soviet flyby probes to comet Halley will carry a French infrared sounder, called “I.K.S.”. In order to assess its observing capabilities, a theoretical model of the comet infrared emission was constructed. We show how the experiment results will be used to derive the nucleus size and radiative properties, and to study the distribution of gas and dust in the inner coma and circumnuclear area. A preliminary discussion is made of the relevance of the data in instances where the cometary phenomena would be more complex than assumed in the model.     

Optical observations of ion tail structures and velocity flows in comet Halley during the period of the five spacecraft encounters

Observations of the distribution and evolution of a number of the major constituents of the neutral coma (CN, C₂, CH, O, H, Na) of Comet Halley were made during two observing periods, each of 3 weeks duration, from the Table Mountain Observatory, California. The first period was pre-perihelion, in late November/December 1985. The second period, from Feb 28 to March 22 1986, covered the five close spacecraft encounters with Halley, and when ICE flew some 20 M Km upstream of Halley. Sodium emission was recorded in early Dec 1985 from the near-nuclear region at a heliocentric distance of 1.4 AU, an observation confirmed with the UCL Doppler Imaging system. The CN coma could be detected to an outer diameter of more than 4M Km in Dec 1985, and 5 – 6M Km in early March 1986, allowing the production of heavy cometary pick-up ions to be estimated. Observations of the cometary ion coma (H₂O⁺ and CO⁺ ions) showed considerable variability from day to day, particularly during the period of the spacecraft encounters. These observations have been used, in conjuction with the neutral coma data, to map the flow field of cometary ions. In early Dec. 1985, Halley developed a traditional “type I” ion tail, which persisted until late April 1986. It has also been possible to evaluate the ion flow fields within the narrow core of the ion tail, and in the surrounding diffuse, low density, regions populated by pick-up and extracted cometary ions, and by slowed solar wind ions. Tail disconnection events were observed on several occasions, particularly between the VEGA 2 and GIOTTO encounters, and with a highly spectacular event on March 19 1986.     

Influence of gravity on the photomorphism of secondary moss protonemata

In dark-grown plantlets of the moss, Pottia intermedia, negatively gravitropic secondary protonemata differentiate from the superficial cells of leafy shoots. When transferred to the light, distal parts of the protonemata nearest to the apical cells begin to ramify and the apical cells of the side branches as well as of the main protonemal filaments often differentiate as buds. Dark-grown protonemata were oriented horizontally and illuminated from below with white light of different intensities. Only light with an intensity of 4.5 μmol·m⁻²·s⁻¹ was sufficient to induce: (a) phototropism in the apical cells, (b) light-directed initiation of branch primordia, and (c) directed growth of side branches and bud differentiation. Apical cells illuminated with light of lower (0.03–0.37 μmol·m⁻²·s⁻¹) intensity grew upwards (i.e., away from the light). It was shown that this upward growth was determined by the action of gravity. Although initiation of branch primordia was only slightly affected, their growth was strongly stimulated on the upper side of the protonemata.     

Collection of cometary dust

Rendezvous Missions to Comets lead to low velocities at the nucleus of the comet. The resulting impact velocity of the cometary dust on a target will range between 10 and 400 m/s. The dust particle which impacts on a target can be collected for a subsequent in-situ analysis.

The collection efficiency of a target depends in addition to obvious geometrical conditions upon the surface of the target. The surface characteristics can be divided into two groups:

• “dirty” surfaces, covered with silicate or hydrocarbon compounds (for example vacuum grease),

• “clean” surfaces, like gold (with additional sputtering).

This paper deals with the experimental and theoretical investigation of the collection efficiency of “clean” targets. Laboratory experiments are described which were conducted at the Technische Universität München, Lehrstuhl für Raumfahrttechnik, and the Max-Planck-Institut für Kernphysik, Heidelberg. In both experiments an electromagnetic accelerator is used to accelerate different types of dust in vacuum to velocities between 10 and 400 m/s.

The target is then examined under the microscope and a secondary ion mass spectrometer (which is a model of the laboratory carried on board of the spacecraft for “in situ” analysis). The adhesion of the dust grains at the target is evaluated experimentally in an ultracentrifuge.     

Comparative study of the dust emission of 19P/Borrelly (Deep space 1) and 1P/Halley

Images obtained by the Miniature Integrated Camera and Imaging Spectrometer (MICAS) experiment onboard the Deep Space 1 spacecraft which encountered comet 19P/Borrelly on September 22nd 2001 show a dust coma dominated by jets. In particular a major collimated dust jet on the sunward side of the nucleus was observed. Our approach to analyse these features is to integrate the observed intensity in concentric envelopes around the nucleus. The same procedures has been used on the Halley Multicolour Camera images of comet 1P/Halley acquired on March 14th 1986. We are able to show that at Borrelly the dust brightness dependence as a function of radial distance is different to that of Halley. At large distances both comets show constant values as the size of the concentric envelopes increases (as one would expect for force free radial outflow). For Halley the integral decreases as one gets closer to the nucleus. Borrelly shows opposite behaviour. The main cause for Halley's intensity distribution is either high optical thickness or particle fragmentation. For Borrelly, we have constructed a simple model of the brightness distribution near the nucleus. This indicates that the influence of deviations from point source geometry is insufficient to explain the observed steepening of the intensity profile close to the nucleus. Dust acceleration or fragmentation into submicron particles appear to be required. We also estimate the dust production rate of Borrelly with respect to Halley and compare their dust to gas ratios.