Numerical simulation of comet nuclei I. Water-ice comets

A one-dimensional simulation of pure water-ice cometary nuclei is presented, and the effect of the nucleus as a heat reservoir is considered. The phase transition from amorphous to crystalline ice is studied for two cases: (1) where the released latent heat goes entirely into heating adjacent layers and (2) where the released latent heat goes entirely into sublimation. For a Halley-like orbit it was found that for case 1 the phase boundary penetrates about 15 m on the first orbit and does not advance until sublimation brings the surface to some 10 m from the phase boundary. For case 2 the phase boundary penetrates about 1 m below the surface and remains at this depth as the surface sublimates. For an orbit like that of Schwassmann-Wachmann 1 the phase boundary penetrates about 50 m initially for case 1 and about 1 m for case 2. There is no further transformation until the entire comet is heated slowly to near the transition temperature, after which the entire nucleus is converted to crystalline ice. For an Encke-type orbit case 1 gives a nearly continuous transition of the entire nucleus to crystalline ice, while for case 2 the initial penetration is about 8 m and remains at this depth relative to the surface as sublimation decreases the cometary radius. Thus the entire comet is converted to crystalline ice just before it is completely dissipated.     

X-ray synchrotron emission from the oblique shock in the jet of the powerful radio galaxy 3C 346

We report the first detection, with Chandra , of X-ray emission from the jet of the powerful narrow-line radio galaxy 3C 346. X-rays are detected from the bright radio and optical knot at which the jet apparently bends by approximately 70°. The Chandra observation also reveals a bright galaxy-scale atmosphere within the previously known cluster and provides a good X-ray spectrum for the bright core of 3C 346. The X-ray emission from the knot is synchrotron radiation, as seen in lower-power sources. In common with these sources, there is evidence of morphological differences between the radio/optical and X-ray structures, and the spectrum is inconsistent with a one-component continuous-injection model. We suggest that the X-ray-bright knot is associated with a strong oblique shock in a moderately relativistic, light jet, at ∼ 20° to the line of sight, and that this shock is caused by the jet interacting with the wake in the cluster medium behind the companion galaxy of 3C 346. The general jet curvature can result from pressure gradients in the cluster atmosphere.     

Einsteins' gravity with massive gravitons

Graviton may, in principle, have a small non-zero mass. In this paper the relevant theory of the massive graviton with six polarisations is developed. The drastic impact of a non-zero mass of the graviton on cosmology is also illustrated.     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

The 550-au Mission: a critical discussion

We have studied the science rationale, goals and requirements for a mission aimed at using the gravitational lensing from the Sun as a way of achieving high angular resolution and high signal amplification. We find that such a mission concept is compromised by several practical problems. Most severe are the effects due to the plasma in the solar atmosphere which cause refraction and scattering of the propagating rays. These effects either limit the frequencies that can be observed to those above ∼1 THz, or they move the optical point outwards beyond the vacuum value of ≥550 au. (Thus for observing frequency of 300 GHz the optical point is moved outwards to ∼ 680 au.) Density fluctuations in the inner solar atmosphere will further cause random pathlength differences for different rays. The corrections for the radiation from the Sun itself will also be a major challenge at any wavelength used, but could be mitigated with coronographic techniques. Given reasonable constraints on the spacecraft (particularly in terms of size and propulsion), source selection as well as severe navigational constraints further add to the difficulties for a potential mission. Nevertheless, unbiased surveys of small-scale structure on the sky at short wavelengths might be the most promising application of such a mission.     

A one-dimensional model for the interaction between continental-scale ice sheets and atmospheric stationary waves

 The great continental ice sheets of the Pleistocene represented a significant topographic obstacle to the westerly winds in northern midlatitudes. This work explores how consequent changes in the atmospheric stationary wave pattern might have affected the shape and growth of the ice sheets themselves. A one dimensional (1-D) model is developed which permits an examination of the types and magnitudes of the feedbacks that might be expected. When plausible temperature perturbations are introduced at the ice-sheet margin which are proportional to the stationary wave amplitude, the equilibrium shape of the ice sheet is significantly altered, and depends on the sign of the perturbation. The proposed feedback also affects the response of the ice sheet to time-varying climate forcing. The results suggest that the evolution of a continental-scale ice sheet with a land-based margin may be significantly determined by the changes it induces in the atmospheric circulation. Received: 1 October 1999 / Accepted: 17 July 2000     

Porous interstellar grains

Recent observations of stellar composition suggest that elements in the Sun are significantly more abundant than in other stars. The reduction in the available element budget implies a drastic revision in current models of interstellar dust. Theoretical models are therefore exploring fluffy, porous physical structure for the grain material. Since a detailed exact treatment of extinction cross-sections is mandatory for a correct understanding of the nature of interstellar dust, we present a technique based on the multipole expansions of the electromagnetic field, which has proven to be general, flexible and powerful in treating scattering of light by porous, composite, arbitrarily shaped particles. The results of this study speak in favour of core–mantle structures characterized by the presence of porosities.     

Theoretical accuracy of synthetic aperture sonar micronavigation using a displaced phase-center antenna

The Cramer-Rao lower bounds on the cross-track translation and rotation of a displaced phase-center antenna (DPCA) in the slant range plane between two successive pings (known as DPCA sway and yaw in what follows) are computed, assuming statistically homogeneous backscatter. These bounds are validated using experimental data from a 118-182-kHz sonar, showing an accuracy of the order of 20 microns on the ping-to-ping cross-track displacements. Next, the accuracy required on the DPCA sway and yaw in order to achieve a given synthetic aperture sonar (SAS) beampattern specification, specified by the expected SAS array gain, is computed as a function of the number P of pings in the SAS. Higher accuracy is required when P increases to counter the accumulation of errors during the integration of the elementary ping-to-ping estimates: the standard deviation must decrease as P/sup -1/2/ for the DPCA sway and P/sup -3/2/ for the yaw. Finally, by combining the above results, the lower bounds on DPCA micronavigation accuracy are established. These bounds set an upper limit to the SAS length achievable in practice. The maximum gain Q in cross-range resolution achievable by a DPCA micronavigated SAS is computed as a function of the key SAS parameters. These theoretical predictions are compared with simulations and experimental results.     

Metamorphism in Archaean greenstone belts: calculated fluid compositions and implications for gold mineralization

An inescapable consequence of the metamorphism of greenstone belt sequences is the release of a large volume of metamorphic fluid of low salinity with chemical characteristics controlled by the mineral assemblages involved in the devolatilization reactions. For mafic and ultramafic sequences, the composition of fluids released at upper greenschist to lower amphibolite facies conditions for the necessary relatively hot geotherm corresponds to those inferred for greenstone gold deposits (X_CO2= 0.2–0.3). This result follows from the calculation of mineral equilibria in the model system CaO–MgO–FeO–Al₂O₃–SiO₂–H₂O–CO₂, using a new, expanded, internally consistent dataset. Greenstone metamorphism cannot have involved much crustal over-thickening, because very shallow levels of greenstone belts are preserved. Such orogeny can be accounted for if compressive deformation of the crust is accompanied by thinning of the mantle lithosphere. In this case, the observed metamorphism, which was contemporaneous with deformation, is of the low-P high-T type. For this type of metamorphism, the metamorphic peak should have occurred earlier at deeper levels in the crust; i.e. the piezothermal array should be of the ‘deeper-earlier’type. However, at shallow crustal levels, the piezothermal array is likely to have been of ‘deeper-later’type, as a consequence of erosion. Thus, while the lower crust reached maximum temperatures, and partially melted to produce the observed granites, mid-crustal levels were releasing fluids prograde into shallow crustal levels that were already retrograde. We propose that these fluids are responsible for the gold mineralization. Thus, the contemporaneity of igneous activity and gold mineralization is a natural consequence of the thermal evolution, and does not mean that the mineralization has to be a consequence of igneous processes. Upward migration of metamorphic fluid, via appropriate structurally controlled pathways, will bring the fluid into contact with mineral assemblages that have equilibrated with a fluid with significantly lower X_CO2. These assemblages are therefore grossly out of equilibrium with the fluid. In the case of infiltrated metabasic rocks, intense carbonation and sulphidation is predicted. If, as seems reasonable, gold is mobilized by the fluid generated by devolatilization, then the combination of processes proposed, most of which are an inevitable consequence of the metamorphism, leads to the formation of greenstone gold deposits predominantly from metamorphic fluids.