Surface composition and physical properties of several trans-neptunian objects from the Hapke scattering theory and Shkuratov model

Several different trans-neptunian objects have been studied in order to investigate their physical and chemical properties. New observations in the 1.1-1.4 μm range, obtained with the ISAAC instrument, are presented in order to complete previous observations carried out with FORS1 in the visible and SINFONI in the near infrared. All of the observations have been performed at the ESO/Very Large Telescope. We analyze the spectra of six different objects (2003 AZ₈₃, Echeclus, Ixion, 2002 AW₁₉₇, 1999 DE₉ and 2003 FY₁₂₈) in the 0.45-2.3 μm range with the model of Hapke (Hapke, B. [1981]. J. Geophys. Res. 86, 4571-4586) and the method of Shkuratov et al. (Shkuratov, Y., Starukhina, L., Hoffmann, H., Arnold, G. [1999]. Icarus 137, 235-246). Water ice is found on two objects, and in particular it is confirmed in its amorphous and crystalline states on 2003 AZ₈₄ surface. Upper limits on the water ice content are given for the other four TNOs investigated, confirming previous results (Barkume, K.M., Brown, M.E., Schaller, E.L. [2008]. Astron. J. 135, 55-67; Guilbert, A., Alvarez-Candal, A., Merlin, F., Barucci, M.A., Dumas, C., de Bergh, C., Delsanti, A. [2009]. Icarus 201, 272-283). Whatever the absorption features in the near infrared, all objects but one exhibit a moderate red slope in the visible, as most TNOs and Centaurs. We discuss the implications of the presence of water ice and the probable sources of the red slope.     

康普顿化的计算-推广的Kompaneets方程与Monte Carlo方法比较

在X射线天文学以及辐射物理学中,当硬X射线穿过" 冷"的等离子体时所发生的Compton软化是一种重要的辐射转移过程.简要介绍推广的Kompaneets方程,该方程在hv〈mec2及kTe〈mec2较宽松的条件下广泛成立,因此不仅能处理Compton硬化过程,而且也适用于Compton软化过程,后者在目前快速发展的X射线和γ射线天文学中十分重要.基于此方程,我们对天体物理中4种常见辐射谱(Gauss型发射谱线、黑体辐射谱、幂律辐射谱和热轫致辐射谱)在Compton软化情况下的谱演化进行了数值求解,并同Monte Carlo模拟结果做比较,证实了推广的Kompaneets方程的正确性和有效性.最后指出此方程在X射线天文学和γ射线天文学中的重要意义和潜在应用.     

Studies of Planetary Atmospheres by Stellar Occultations: Application to Mars and Venus

We consider the application of the stellar occultation method to the studies of planetary atmospheres and its history and briefly describe the instruments designed for such measurements (SPICAM/Mars-96, GOMOS/ENVISAT). In comparison with solar occultations, this method allows the profiles to be measured almost at any time of the day and at any location of the planet, irrespective of the orbit of the spacecraft from which observations are carried out. Based on the measuring characteristics of the SPICAM-Light UV spectrometer for the spectral range 118–320 nm with a resolution of 0.9 nm (for the ESA Mars Express Mission; launched in June 2003), we simulate the capabilities of the method to study the Martian atmosphere. In stellar occultation measurements, the stellar spectrum changes because of the absorption by CO₂ and O₃, other gases, and aerosols. The profiles of the CO₂ and O₃ density (and, hence, the temperature) and the aerosol content can be restored by solving the inverse problem. Observations of bright stars (no fewer than 30) three to five times in a turn allow us to measure the atmospheric density at altitudes 10–150 km with an accuracy of about 2% and the temperature at altitudes 20–130 km with an accuracy of 3 K. Ozone is measured with an accuracy of several percent at altitudes 25–40 km or lower, depending on the conditions. Optically thin clouds and hazes, particularly on the nightside where no measurements are possible in reflected light, can be studied. The SPICAV experiment, which is similar to SPICAM-Light, is part of the Venus Express (to be launched in 2005) scientific payload. On Venus, stellar occultations can be used to measure the atmospheric temperature and density above clouds at altitudes up to 130–150 km and to study the SO₂ profile. The results of our simulations can be easily extended to instruments with different measuring characteristics.     

Dynamics of Two Planets in the 3/2 Mean-motion Resonance: Application to the Planetary System of the Pulsar PSR B1257+12

This paper considers the dynamics of two planets, as the planets B and C of the pulsar PSR B1257+12, near a 3/2 mean-motion resonance. A two-degrees-of-freedom model, in the framework of the general three-body planar problem, is used and the solutions are analyzed through surfaces of section and Fourier techniques in the full phase space of the system.On the possibility of existence of a fourth planet in distant orbit, see Wolszczan et al., 2000     

Models for the X‐ray spectra and variability of luminous accreting black holes

The X‐ray spectra of luminous Seyfert 1 galaxies often appear to be reflection dominated. In a number of Narrow Line Seyfert 1 (NLS1) galaxies and galactic black holes in the very high state, the variability of the continuum and of the iron line are decoupled, the reflected component being often much less variable than the continuum. These properties have been interpreted as effects of gravitational light bending. In this framework, we present detailed Monte‐Carlo simulations of the reflection continuum in the Kerr metric. These calculations confirm that the spectra and variability behaviour of these sources can be reproduced by the light bending model. As an alternative to the light bending model, we show that similar observational properties are expected from radiation pressure dominated discs subject to violent clumping instabilities and, as a result, have a highly inhomogeneous two‐phase structure. In this model, most of the observed spectral and variability features originate from the complex geometrical structure of the inner regions of near‐Eddington accretion flows and are therefore a signature of accretion physics rather than general relativity. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

行星大气定理及其天灾预测应用

在界定行星大气运动范围后,提出行星大气运动三定理。从行星磁场角度研究不完全电离磁流体(incompleteionospheregaseousfluid)得到行星磁力线绝热压缩反弹致冷效应,给出大气运动出现的冷池(coolpool)现象的物理机制。最后从天文学角度讨论和展望大气灾害的可预报性。     

Hapke mixture modeling applied to VNIR spectra of mafic mineral mixtures and shergottites: Implications for quantitative analysis of satellite data

The mineralogy of Mars is well understood on a qualitative level at a global scale due to satellite data. Quantitative analysis of visible and near-infrared (VNIR) satellite data is a desirable but nontrivial task, due partly to the nonlinearity of VNIR reflectance spectra from the mineral mixtures of the Martian surface. In this study, we investigated the use of the Hapke radiative transfer model to generate linearly mixed single scattering albedo data from nonlinearly mixed VNIR reflectance data and then quantitatively analyzed them using the linear spectral mixture model. Simplifications to the Hapke equation were tested accounting for variables that would be unknown when using satellite data. Mineral mixture spectra from the RELAB spectral library were degraded to test the robustness of the unmixing technique in the face of data that mimic some of the complexities of satellite spectral data collected at Mars. A final test was performed on spectra from shergottite meteorites to assess the technique against real Martian mineral mixtures. The simplified Hapke routine produced robust abundance estimates within 5–10% accuracy when applied to laboratory standard spectra from the synthetic mixtures of igneous minerals in agreement with previous studies. The results of tests involving degraded data to mimic the low spectral contrast of the Martian surface and the lack of a priori knowledge of the constituent mineral spectral endmembers, however, were less encouraging, with errors in abundance estimation greater than 25%. These results cast doubt on the utility of Hapke unmixing for the quantitative analysis of VNIR data of the surface of Mars.     

The total and relative contribution of the relevant absorption processes to the opacity of DB white dwarf atmospheres in the UV and VUV regions

The main aim of this work is to estimate the total contribution of the processes of     molecular ion photodissociation and     collisional absorption charge exchange to the opacity of DB white dwarf atmospheres, and compare this with the contribution of     and other relevant radiative absorption processes included in standard models.
The method for the calculations of the molecular ion     photodissociation cross-sections is based on the dipole approximation and quantum-mechanical treatment of the internuclear motion, while the quasi-classical method for describing absorption processes in     collisions is based on the quasi-static approximation.
Absorption coefficients are calculated in the region  50 nm ≤λ≤ 850 nm  and compared with the corresponding coefficients of other relevant absorption processes; the calculations of the optical depth of the atmosphere layers considered are performed in the far-UV and VUV regions; the contribution of the relevant absorption processes to the opacity of DB white dwarf atmospheres is examined.
We examined the spectral ranges in which the total     and     absorption processes dominate in particular layers of DB white dwarf atmospheres. In addition, we show that in the region of  λ≲ 70 nm  the process of     atom photoionization is also important, in spite of the fact that the ratio of hydrogen and helium abundances in the DB white dwarf atmosphere considered is  1:10⁵  .     

Brightness fluctuations of the Sun and p‐mode oscillations: The inverse problem and nonadiabatic waves in the photosphere

Using photometric observations of the Sun as a star (DIFOS, SoHO) we were able to solve the inverse heloiseismic problem and determine the global time‐dependent relative temperature fluctuations as functions of the geometric height. This was done under the adiabatic assumption. A mathematical tool was developed to solve the inverse problem, which is ill‐posed. The calculations were done using the numerical software Matlab 7. The adiabatic solution shows signs of temperature waves in the lower photosphere, which agrees with calculations done by Rodríguez Hidalgo et al. (2001) and Stodilka (2011). (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stokes imaging of AM Her systems using 3D inhomogeneous models – I. Description of the code and an application to V834 Cen

Polars (or AM Her systems) are cataclysmic variables without a disc due to the strong magnetic field of the white dwarf. Most of their emission comes from the region where the accretion column impacts the white dwarf and cools through cyclotron and bremsstrahlung processes. We present a new code, cyclops , to model the optical emission from these systems including the four Stokes parameters. It considers a three-dimensional region with the electronic density and temperature varying following a shock-like profile and a dipolar magnetic field. The radiative transfer is solved in steps considering the solution with non-null input radiation. The footprint of the column in the white dwarf surface is determined by the threading region in the equatorial plane, i.e. the region from where the flow follows the magnetic lines. The extinction caused by the Thomson scattering above the emitting region is optionally included. The search for the model parameters that best fit an observational data set is carried out using a hybrid approach: a genetic algorithm is used to seek for the regions of the parameter space having the best models and then an amoeba code refines the search. An example of the application to multi-wavelength data of V834 Cen is presented. The fit found is consistent with previous parameter estimates and is able to reproduce the features of V834 data in three wavebands.     

Massive star formation in galaxies: radiative transfer models of the UV to millimetre emission of starburst galaxies

We present illustrative models for the UV to millimetre emission of starburst galaxies which are treated as an ensemble of optically thick giant molecular clouds (GMCs) centrally illuminated by recently formed stars. The models follow the evolution of the GMCs owing to the ionization-induced expansion of the H  ii regions and the evolution of the stellar population within the GMC according to the Bruzual & Charlot stellar population synthesis models. The effect of transiently heated dust grains/PAHs on the radiative transfer, as well as multiple scattering, is taken into account.
The expansion of the H  ii regions and the formation of a narrow neutral shell naturally explain why the emission from PAHs dominates over that from hot dust in the near- to mid-infrared, an emerging characteristic of the infrared spectra of starburst galaxies.
The models allow us to relate the observed properties of a galaxy to its age and star formation history. We find that exponentially decaying 10⁷–10⁸ yr old bursts can explain the IRAS colours of starburst galaxies. The models are also shown to account satisfactorily for the multiwavelength data on the prototypical starburst galaxy M82 and NGC 6090, a starburst galaxy recently observed by ISO . In M82 we find evidence for two bursts separated by 10⁷ yr. In NGC 6090 we find that at least part of the far-infrared excess may be due to the age of the burst (6.4×10⁷ yr). We also make predictions about the evolution of the luminosity of starbursts at different wavelengths which indicate that far-infrared surveys may preferentially detect older starbursts than mid-infrared surveys.     

Recombination lines and free‐free continua formed in asymptotic ionized winds: Analytic solution for the radiative transfer

In dense hot star winds, the infrared and radio continua are dominated by free‐free opacity and recombination emission line spectra. In the case of a spherically symmetric outflow that is isothermal and expanding at constant radial speed, the radiative transfer for the continuum emission from a dense wind is analytic. Even the emission profile shape for a recombination line can be derived. Key to these derivations is that the opacity scales with only the square of the density. These results are well‐known. Here an extension of the derivation is developed that also allows for line blends and the inclusion of an additional power‐law dependence beyond just the density dependence. The additional power‐law is promoted as a representation of a radius dependent clumping factor. It is shown that differences in the line widths and equivalent widths of the emission lines depend on the steepness of the clumping power‐law. Assuming relative level populations in LTE in the upper levels of He II, an illustrative application of the model to Spitzer/IRS spectral data of the carbon‐rich star WR 90 is given (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effects of scattering and dust grain size on the temperature structure of protoplanetary discs: a three-layer approach

The temperature in the optically thick interior of protoplanetary discs is essential for the interpretation of millimetre observations of the discs, for the vertical structure of the discs, for models of the disc evolution and the planet formation, and for the chemistry in the discs. Since large icy grains have a large albedo even in the infrared, the effect of scattering of the diffuse radiation in the discs on the interior temperature should be examined. We have performed a series of numerical radiation transfer simulations, including isotropic scattering by grains with various typical sizes for the diffuse radiation as well as for the incident stellar radiation. We also have developed an analytic model including isotropic scattering to understand the physics concealed in the numerical results. With the analytic model, we have shown that the standard two-layer approach is valid only for grey opacity (i.e. grain size ≳10 μm) even without scattering. A three-layer interpretation is required for grain size ≲10 μm. When the grain size is 0.1–10 μm, the numerical simulations show that the isotropic scattering reduces the temperature of the disc interior. This reduction is nicely explained by the analytic three-layer model as a result of the energy loss by scatterings of the incident stellar radiation and of the warm diffuse radiation in the disc atmosphere. For grain size ≳10 μm (i.e. grey scattering), the numerical simulations show that the isotropic scattering does not affect the interior temperature. This is nicely explained by the analytic two-layer model; the energy loss by scattering in the disc atmosphere is exactly offset by the 'green-house effect' due to the scattering of the cold diffuse radiation in the interior.     

Surface slope probabilities from the spectra of weak radar echoes: Application to Mars

Slope probability densities were derived from the power spectra of radar echoes from Mars using integral inversion. The inverse problem is ill-posed; that is, small changes in the data can lead to large changes in the solution. We describe a method of stabilizing the inversion, which was necessary for echoes with signal-to-noise power spectral densities on the order of unity, and for those with broad spectral distributions. The resulting slope probabilities usually consisted of a component due to quasi-specular reflection which decreased rapidly with tilt, plus a broad, slowly decreasing, “diffuse” component due to scattering from (1) surface scales small compared with a radar wavelength, or (2) larger features with high slopes. In the absence of more complete polarization measurements, we are unable to In the absence of more complete polarization measurements, we are unable to distinguish between these possibilities. Root mean square tilts have been determined separately for the two cases. For case (1), values of rms tilt associated with surface features responsible for the quasi-specular echo are normally less than 3°; for case (2), values greater than 8° are common. Knowledge of the depolarization of radar returns would help distinguish between these possibilities.     

Computer Simulation of Sputtering of Lunar Regolith by Solar Wind Protons: Contribution to Change of Surface Composition and to Hydrogen Flux at the Lunar Poles

A computer simulation of the sputtering of lunar soil by solar wind protons was performed with the TRIM program. The rate of the sputtering-induced erosion of regolith particles was shown to be less than 0.2 Å per year. A preferential sputtering of Ca, Mg, and O was found along with a less intense sputtering of Fe, Si, and Ti. However, with no other selection mechanisms, surface concentrations of the atoms would differ from the volume ones by no more than 6 %. The enrichment of rims of regolith particles with iron occurs as a result of selective removal of lighter atoms from the lunar surface because of different energies of escape from the Moon's gravity. The energy distributions proved to be the same for all sorts of the sputtered atoms, except for implanted hydrogen; thus, a greater fraction of the atoms left on the lunar surface corresponds to heavier elements. According to simulation results, the concentration of reduced iron observed in the mature regolith could be attained during the time of regolith particle exposure to the present flux of solar wind (10⁵ years). Thus, sputtering can provide the concentration of Fe₀ observed in regolith. On periphery of a cloud of impact vapor the temperature is too low for an irreversible selective removal of evaporation products; thus, a meteoritic bombardment contributes to the formation of composition of the rims of regolith particles mainly through enrichment of the rims with elements from the bulk of the particles. The estimates of fluxes of backscattered solar wind protons and of sputtered protons, earlier implanted to the regolith, demonstrated that their contribution to the proton flux near the poles is only 10⁴ cm^–2 s^–1. This is by two orders of magnitude smaller than the proton flux from the Earth's magnetosphere which is, therefore, the main source of protons for permanently shaded polar craters of the Moon.