Metastability of Liquid Water on Mars

A simple model of local heat transport on Mars demonstrates that transient melting of ice may occur in depressions and gullies nearly anywhere on the planet where thin ice is illuminated by normal-incidence insolation. An experiment has been performed to confirm the model of evaporation rate at low pressure. Reduction of radiative cooling due to gully geometry is shown to be important. Since appropriate meteorological, topographic, and optical conditions may occur on slopes nearly anywhere on the planet, hydrological features such as gullies would likely form only where such ice might accumulate, notably in sheltered locations at high latitudes. It is suggested that cold-trapping of winter condensation could concentrate a sufficient amount of ice to allow seasonal melting in gullies.     

Accretion disc dynamics in extreme mass ratio compact binaries

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios   q = M ₂/ M ₁ < 0.1  , applicable to ultracompact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For   q ∼ 0.1  , the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q . At   q = 0.01  , the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com, which has   q = 0.02  . The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultracompact binaries.     

Investigating thermal evolution of the self-gravitating one dimensional molecular cloud by smoothed particle hydrodynamics

The heating of the ion-neutral (or ambipolar) diffusion may affect the thermal phases of the molecular clouds. We present an investigation on the effect of this heating mechanism in the thermal instability of the molecular clouds. A weakly ionized one-dimensional slab geometry, which is allowed for self-gravity and ambipolar diffusion, is chosen to study its thermal phases. We use the thermodynamic evolution of the slab to obtain the regions where slab cloud becomes thermally unstable. We investigate this evolution using the model of ambipolar diffusion with two-fluid smoothed particle hydrodynamics, as outlined by Hosking and Whitworth. Firstly, some parts of the technique are improved to test the pioneer works on behavior of the ambipolar diffusion in an isothermal self-gravitating slab. Afterwards, the improved two-fluid technique is used for thermal evolution of the slab. The results show that the thermal instability may persist inhomogeneities with a large density contrast at the intermediate parts of the cloud. We suggest that this feature may be responsible for the planet formation in the intermediate regions of a collapsing molecular cloud and/or may also be relevant to the formation of star forming dense cores in the clumps.     

A new photometric analysis of the eclipsing binary system ST Carinae

ST Carinae is an eclipsing binary with a period of 0 _. ^d 90165 which is believed to consist of an A0V primary and a secondary of type F5-8IV. About 900 observations inUBV, made by Somerville in 1963 but unreduced at that time, are presented. The Wood model is used to obtain orbital elements, and four different solutions of the light curves are presented. These are also computed with the solutions obtained by previous investigators of this system. The solutions indicate a reasonably consistent geometry, but there is still substantial uncertainty with regard to the mass-ratio and dynamical status of ST Car. The possibility exists that ST Car is in the initial and rapid stage of mass exchange in which the primary component fills its Roche lobe and is losing mass to its companion. The two components of ST Car appear to be of normal dimensions, but additional work is needed to clarify the exact status of this system.     

Probing shock geometry via the charge to mass ratio dependence of heavy ion spectra from multiple spacecraft observations of the 2013 November 4 event

In large Solar Energetic Particle(SEP) events, ions can be accelerated at coronal mass ejection(CME)-driven shocks to very high energies. The spectra of heavy ions in many large SEP events show features such as roll-overs or spectral breaks. In some events when the spectra are plotted in terms of energy/nucleon, they can be shifted relative to each other to make the spectral breaks align. The amount of shift is charge to mass ratio(Q/A) dependent and varies from event to event. This can be understood if the spectra of heavy ions are organized by the diffusion coefficients(Cohen et al. 2005). In the work of Li et al.(2009), the Q/A dependence of the scaling is related to shock geometry when the CME-driven shock is close to the Sun. For events where multiple in-situ spacecraft observations exist, one may expect that different spacecraft are connected to different portions of the CME-driven shock that have different shock geometries, therefore yielding different Q/A dependence. In this work, we examine one SEP event which occurred on 2013 November 4. We study the Q/A dependence of the energy scaling for heavy ion spectra using helium, oxygen and iron ions. Observations from STEREO-A, STEREO-B and ACE are examined. We find that the scalings are different for different spacecraft. We suggest that this is because ACE, STEREOA and STEREO-B are connected to different parts of the shock that have different shock geometries. Our analysis indicates that studying the Q/A scaling of in-situ particle spectra can serve as a powerful tool to remotely examine the shock geometry for large SEP events.     

Properties of the hermean regolith: iii. disk-resolved vis-NIR reflectance spectra and implications for the abundance of iron

Disk-resolved reflectance spectra of the surface of Mercury (longitudes 240-300°), obtained in the visual (vis) and near-infrared (NIR) spectral region, are presented and analyzed. The observations were made at the 2.6-m Nordic Optical Telescope with the ALFOSC low-resolution spectrograph on 20 and 22 June 1999 in the wavelength range 520-970 nm with a footprint size of 700 km on the mid-disk of Mercury. A method which enables more accurate correction for telluric line absorptions and atmospheric extinction than that applied on previously published vis-NIR spectra of Mercury is introduced. The resulting reflectance spectra are remarkably linear, lack significant absorption features, and have optical slopes comparable to remotely sensed lunar pure anorthosites. The relation between spectral slope and photometric geometry found by Warell (2002, Icarus 156, 313-317) is confirmed and is explained as caused by strongly backscattering particles with embedded submicroscopic metallic iron in a mature regolith. With the theoretical maturation model of Hapke (2001, J. Geophys. Res. 106 (E5), 10039-10073) an abundance of 0.05-0.3 wt% submicroscopic metallic iron in the regolith for silicate grain sizes in the range 10-80 μm is determined, implying a ferrous iron content in mafic minerals intrinsically lower than that of the lunar highlands. A binary crustal composition model with anorthite linearly mixed with pyroxene provides better spectral fits than a pure anorthitic composition. Comparison with mature lunar pure anorthosite spectra yields a confident upper limit to the FeO content of 3 wt% under the assumption that the surfaces are similarly matured, but this figure probably represents a considerable overestimate. The average mercurian regolith does not seem to be substantially more weathered than the most mature lunar highland soils in terms of abundance of submicroscopic metallic iron, indicating that a steady-state maturation level has been reached. However, the strong relation between optical spectral slope and photometric geometry may imply that the majority of regolith particles are more fine-grained than their lunar counterparts and that the regolith is admixed with complex agglutinate weathering products which are more abundant and more transparent than those of the lunar highlands. This is consistent with more energetic impacts and a higher rate of impact melt production in an iron-poor regolith. An observed relation between the spectral slope and latitude provides evidence that the Ostwald ripening process may be operating at equatorial latitudes on Mercury.     

Light Curve Analysis of Early-Type Overcontact Systems

A method is presented to include irradiation effects in the modeling of over-contact binaries. The Roche potential is numerically modified in a way to account for mechanical effects of the mutual irradiation of the two binary components. The efficiency of radiative interaction is parametrized by means of the fraction of radiative over gravitational forces for each component of the binary. The modified Roche geometry is implemented as binary model in an eclipse light curve solution code, which is based on the general Wilson–Devinney scheme. As an example the method is applied to the early-type over-contact system V606 Cen.     

Bayesian calibrated significance levels applied to the spectral tilt and hemispherical asymmetry

Bayesian model selection provides a formal method of determining the level of support for new parameters in a model. However, if there is not a specific enough underlying physical motivation for the new parameters it can be hard to assign them meaningful priors, an essential ingredient of Bayesian model selection. Here we look at methods maximizing the prior so as to work out what is the maximum support the data could give for the new parameters. If the maximum support is not high enough then one can confidently conclude that the new parameters are unnecessary without needing to worry that some other prior may make them significant. We discuss a computationally efficient means of doing this which involves mapping p-values on to upper bounds of the Bayes factor (or odds) for the new parameters. A p-value of 0.05 (1.96σ) corresponds to odds less than or equal to 5:2, which is below the 'weak' support at best threshold. A p-value of 0.0003 (3.6σ) corresponds to odds of less than or equal to 150:1, which is the 'strong' support at best threshold. Applying this method we find that the odds on the scalar spectral index being different from one are 49:1 at best. We also find that the odds that there is primordial hemispherical asymmetry in the cosmic microwave background are 9:1 at best.     

Diffusion of magnetic flux elements on a fractal geometry

Recent observations have indicated that magnetic field elements are distributed on the Sun in fractal patterns with dimension D < 2. We suggest that the transport of magnetic field elements across the solar surface should be treated as diffusion on a fractal geometry. We review a semi-analytical, theoretical treatment of fractal diffusion. Comparison with observations of small-scale motions of solar magnetic flux concentrations indicates that fractal diffusion may be taking place with dimension in the range 1.3 to 1.8. It is shown that, compared to the predictions that would be made for two-dimensional diffusion, fractal diffusion in this range would lead to an increased level of in situ flux cancellation in decaying active regions by 7% to 35%. Other work in specialities outside of solar physics may be useful in explaining solar magnetic phenomena.     

Amorphous carbon around carbon stars

In this paper we present the results of a simplified model to determine the flux emerging from dust envelopes around cool stars. The model proposed holds under the hypotheses of negligible scattering effects and spherical geometry of the dust cloud.The aim of this work is to compare the effects of a graphitic or amorphous composition of the carbon grains in the envelopes. To do this we have used, for the first time, experimental extinction data obtained in the laboratory for submicron amorphous carbon particles.The model has been used to fit the FIR spectral trend of 78 optically thin sources and to reproduce the full spectra of two of the most IR luminous optically-thick sources: CIT 6 and IRC+10216.Our calculations indicate clearly that solid carbon particles around these sources may be amorphous rather than crystalline.     

Steady dynamos in finite domains: an integral equation approach

The paper deals with the integral equation approach to steady kinematic dynamo models in finite domains based on Biot‐Savart's law. The role of the electric potential at the boundary is worked out explicitly. As an example, a modified version of the simple spherical α‐effect dynamo model proposed by Krause and Steenbeck is considered in which the α‐coeffcient is no longer constant but may vary with the radial coordinate. In particular, the results for the original model are re‐derived. Possible applications of this integral equation approach for numerical simulations of dynamos in arbitrary geometry and for an “inverse dynamo theory” are sketched.     

The evolved B[e] star HD 87643: observations and a radiation-driven disc wind model for B[e] stars

New high-resolution spectroscopic and medium-resolution spectropolarimetric data of the B[e] star HD 87643 are presented, complemented with optical broad- and narrow-band imaging. The spectrum of HD 87643 exhibits the hybrid characteristics well known to be representative of the group of B[e] stars; a fast wind with an expansion velocity in excess of 1000 km s⁻¹ is measured in the hydrogen and helium lines, while a slower component is traced by lower excitation lines and forbidden lines. Clues to the geometry of the rapidly expanding circumstellar shell are provided by the startling polarization changes across Hα. Comparison with published schematic calculations indicates that the polarizing material is located in a slowly rotating, expanding disc structure. A hydrodynamical model is then presented, the results of which are consistent with the original two-wind concept for B[e] stars, and which exhibits kinematic properties that may well explain the observed spectral features in HD 87643. The model calculations use as input a B star undergoing mass loss, surrounded by an optically thick disc. The resulting configuration consists of a fast polar wind from the star and a slowly expanding disc wind. The model also predicts that the stellar wind at intermediate latitudes is slower and denser than in the polar region.     

Naturally occurring continued fractions in the variation of Kepler's equation

A family of functions involving integrals of universal functions is introduced. These functions have some interesting mathematical properties including the fact that they may be expressed as Gaussian continued fractions. An unique method of performing the integration is demonstrated which indicates why these functions may be important in the variation of Kepler's equation.This work was supported at the Charles Stark Draper Laboratory, Inc. by the National Aeronautics and Space Administration under Contract NAS9-17560.     

Real-time Solution of Atomic Clock Operating Parameters Based on Adaptive Kalman Filtering

The real-tine and accurate knowledge of the operation of an atomic clock is of great significance to satellite navigation and time keeping. In the light of the problems which exist in the real-time solution of the atomic clock operating parameters by means of the standard Kalman filtering, the method for applying the adaptive Kalman filtering model based on the prediction residual to the real-time solution of the clock parameters is proposed. An adaptive factor is constructed by the model based on the standard Kalman filtering model so that the accurate and reliable estimates of the operating parameters may be given rather fast after the frequency modulation or phase modulation of the atomic clock. Finally, in combination with the 3 computational examples of the work on the navigation and time keeping, the effectiveness of the method is verified.     

Polarization evolution accompanying the very early sharp decline of gamma-ray burst X-ray afterglows

In the synchrotron radiation model, the polarization property depends on both the configuration of the magnetic field and the geometry of the visible emitting region. Some peculiar behaviours in the X-ray afterglows of gamma-ray bursts (GRBs) observed with Swift , such as energetic flares and a plateau followed by a sharp drop, might be highly linearly polarized because the outflows powering these behaviours may be dominated by Poynting flux. The breakdown of the symmetry of the visible emitting region may also be well hidden in the peculiar X-ray data and may give rise to interesting polarization signatures. In this paper, we focus on the polarization accompanying the very early sharp decline of GRB X-ray afterglows. We show that strong polarization evolution is possible in both the high latitude emission model and the dying central engine model, which are used to interpret this sharp X-ray decline. It is thus not easy to efficiently probe the physical origin of the very early X-ray sharp decline with future polarimetry. Strong polarization evolution is also possible in the decline phase of X-ray flares and in the shallow decline phase of X-ray light curves characterized by chromatic X-ray versus optical breaks. A detector such as the X-ray Telescope (XRT), but with polarization capability, on board a satellite like Swift would be suitable for testing our predictions.     

The correlation between the Calern solar diameter measurements and the solar irradiance

The objective of this paper is to present some results deduced from the analysis of (space-based) solar irradiance observations used jointly with (ground-based) solar diameter variations. The idea which is explored consists in searching a possible influence of the variability of the Sun's whole shape on the luminosity. It is shown that such an effect, albeit small, may occur. Thus, the global geometry of the Sun – which is not a perfect static ellipsoid – would have to be taken into account when attempting to model the irradiance. Our very preliminary results may help to construct empirical models that can be used, in turn to force any model of the thermal structure of the ocean and atmosphere to deduce climate variations, if any.     

Propagation of two dimensional cylindrical fast magnetoacoustic solitary waves in a warm dust plasma

A set of multi-fluid equations and Maxwell’s equations are carried out to investigate the properties of nonlinear fast magnetoacoustic solitary waves with the combined effects of dusty plasma pressure and transverse perturbation in the bounded cylindrical geometry. The reductive perturbation method has been applied to the dynamical system causeway and the derived two dimensional cylindrical Kadomtsev–Petviashvili equation (CKP) predicts different natures of solitons in complex plasma. Under a suitable coordinate transformation the CKP equation can be solved analytically. The change in the soliton structure due to mass of dust, ion temperature, ion density, and dust temperature is studied by numerical calculation of the CKP equation. It is noted that the dust cylindrical fast magnetoacoustic solitary waves in warm plasmas may disappear slowly because of an increase in dust mass. The present analysis could be helpful for understanding the nonlinear ion-acoustic solitary waves propagating in interstellar medium and pulsar wind,which contain an excess of superthermal particles.     

On the depths and shapes of the freshest kilometer-scale simple craters on the lunar maria: A new crater shape model

Recent work on the shapes of small, simple impact craters on the Moon has shown that the parabolic ideal does not well represent the vast majority of these craters. They are hyperbolic in shape and usually resemble a cone more than a parabola. A parabolic shape also does not fit the most commonly held archetype for simple craters in general (Linné), which is also hyperbolic. In addition, Linné itself may not be the best model for fresh simple craters, in terms of cross-sectional shape, although shape data to compare it to have heretofore been lacking. Here, the “free shadowfront method” for determining the shapes of simple craters is used to measure 64 fresh simple craters on five lunar maria to test both assumptions. Laser altimetry cross sections, available for many of the craters measured herein, are used to complement and spot-check the shadow measurement results, and thereby demonstrate the efficacy of the free shadowfront method. A new shape model is established, and two craters that better fit this model than Linné are identified. These are located at 24.45° N/328.12° E and 31.35° N/296.46° E and have diameters of 1.40 and 2.73 km, respectively. An apparent dichotomy between fresh simple craters smaller than 2.5 km and those larger than this is observed. Flat floors are found to be ubiquitous among the larger craters, but rare and small in extent in smaller ones. A slide in one crater which appears to be an incipient flat floor suggests a major mode of formation for these flat floors.     

Higher dimensional inhomogeneous cosmological models in Lyra geometry

We have presented in homogeneous cosmological models within the framework of Lyra geometry. We have considered an inhomogeneous spherically symmetric higher dimensional model in presence of a mass less scalar field whose potential has a flat part. The scalar field is considered to be inhomogeneous. Also an inhomogeneous cosmological model is derived in a Kaluza-Klein type of space time. The matter field is taken as an inhomogeneous distribution of fluid. It is observed that there is no singularity at finite past in our model and the desirable feature of dimensional reduction is also possible for the extra space. This revised version was published online in July 2006 with corrections to the Cover Date.     

On disk accretion

A simple but sufficiently accurate method for calculating an accretion disk structure is developed. The detailed analysis of the accretion disk is fulfilled by using this method. The effect of turbulent heat transfer on the disk structure is taken into account along with the effect of radiative transfer. The turbulent heat transfer is shown to play an important role, and may be even predominant in the inner disk region. The dependence of temperatureT and density on thez-coordinate is found. Simple analytic expressions are proposed for the run of the density in all the disk zones. It is shown that the inner disk region is convectively stable. The main parameters of all the zones are derived. The geometry of regular motions is studied; the regular hydrodynamical flows are found to appear directed both toward and outside the central object. A detailed comparison is made with the results of other authors.