CCD photometry of variable stars in the Magellanic Clouds — VI. The eclipsing binary HV 982 in the Large Magellanic Cloud

Using improved techniques, high-quality CCD uvbyVI photometry has been obtained for the eclipsing binary HV 982 in the Large Magellanic Cloud (LMC). International Ultraviolet Explorer ultraviolet spectrophotometry was also obtained. These data have been analysed using the Wilson–Devinney synthetic light-curve code and Kurucz low-metallicity model atmospheres as well as the EBOP code. The system is detached and the orbit is eccentric. Apsidal motion is detected with apsidal period 205 ± 7 yr. The effective temperatures of the components are found via flux fitting to be T _eff,1 = 28 000 ± 5000 K and T _eff,2 = 27 200 ± 5000 K. The large errors result from uncertainties over the appropriate interstellar extinction correction. The system plausibly comprises two ∼ 8 M stars of radius 6–7 R separated by ∼ 30 R. For pedagogical and historical interest, the near simultaneity of the eclipse minima at different wavelengths is used to constrain the constancy of the speed of light with wavelength and the mass of the photon, yielding m _γ < 10⁻⁴¹ kg. Because of the great distance to HV 982, this limit is some 10² times smaller than previously achieved with eclipse timings, but it is nevertheless 10 orders of magnitude less stringent than that which is provided by satellite measurements of planetary magnetic fields.     

Influence of the YORP effect on rotation rates of near‐Earth asteroids

The distribution of near‐Earth asteroid (NEA) rotation rates differs considerably from the similar distribution of Main Belt asteroids (MBAs) by the presence of excesses of fast and slow rotators, which are not observed or not so prominent in the distribution for MBAs. Among possible reasons for the difference, there can be influence of solar radiation on spin rate of small NEAs, the so‐called “YORP effect,” which appears due to reflection, absorption, and IR re‐emission of the sunlight by an irregularly shaped rotating asteroid. It is known that the YORP‐effect action strongly depends on the amount of solar energy obtained by the body (insolation), its size, and albedo. The analysis of observation data has shown that: (1) the mean diameter of NEAs decreases from the middle of the distribution to its ends, that is, the excesses of slow rotators (ω ≤ 2 rev day^?1) and fast rotators (ω ≥ 8 rev day^?1) are composed of smaller NEAs than in the middle of the distribution; (2) NEAs of both excesses are in the orbits where their insolation is about 8–10% larger than that of NEAs in the middle of the distribution; and (3) the objects in both excesses have a little lower albedo on average than that of objects in the middle of the distribution. All these results qualitatively agree well with the YORP‐effect action and may be considered as independent arguments in favor of it.     

The orbital evolution of NEA 30825 1900 TG1

The orbital evolution of the near-Earth asteroid (NEA) 30825 1990 TG1 has been studied by numerical integration of the equations of its motion over the 100 000-year time interval with allowance for perturbations from eight major planets and Pluto, and the variations in its osculating orbit over this time interval were determined. The numerical integrations were performed using two methods: the Bulirsch-Stoer method and the Everhart method. The comparative analysis of the two resulting orbital evolutions of motion is presented for the time interval examined. The evolution of the asteroid motion is qualitatively the same for both variants, but the rate of evolution of the orbital elements is different. Our research confirms the known fact that the application of different integrators to the study of the long-term evolution of the NEA orbit may lead to different evolution tracks.     

The Recovery as an Important Part of NEA Astrometric Follow-up

The number of known near-Earth asteroids (NEAs) has increased rapidly in recent years due to large surveys. This discovery process has to be followed by follow-up observations to obtain a sufficient number of precise astrometric data needed for an accurate orbit determination of newly discovered bodies.Accurate orbit determination requires observations from at least two oppositions. If asteroids are not found in the next apparition, different from the discovery apparition, then they can be considered lost. This is particularly embarrassing for NEAs. If data for different apparitions are not found in the course of precovery surveys or in other archive data, then it is necessary to prepare targeted observations of a particular NEA in the second convenient apparition. Therefore NEA recovery is a very important part of NEA follow-up.We discuss here methods, techniques, and results of planned recoveries at the Klet' Observatory using a 0.57-m telescope equipped with a CCD detector. The Klet' NEA recovery subprogram has brought 21 planned NEA recoveries since 1997, including seven NEAs belonging to the potentially hazardous asteroid category.We briefly mention the overall work on NEA recoveries provided by several NEO follow-up programs as well as the need for communication resources supporting astrometric observers. Finally we present here a planned extension of the Klet' NEA recovery subprogram to fainter objects by means of a new 1.06-m reflector.     

The effect of non-gravitational forces on the median inclination of short-period comets

The numbered Jupiter family comets (orbital periods   P < 20 yr  ) have a median orbital inclination of about     . In this paper, we integrate the orbits of these comets into the future, under the influence of both typical non-gravitational forces and planetary perturbation, using a Bulirsch–Stoer integrator. In the case where non-gravitational forces were not acting, the median inclination of those comets that remained on   P < 20 yr  orbits increased at the rate of  (1.92 ± 0.12) × 10⁻³ deg yr⁻¹  for the first 3600 yr of the integration. During this time the population of the original family decreases, such that the half-life is about 13 200 ± 800 yr. The introduction of non-gravitational forces slows down the rate of increase in inclination to a value of around  (1.23 ± 0.16) × 10⁻³ deg yr⁻¹  . This rate of increase in inclination was found to be only weakly dependent on the non-gravitational parameters used during the integration. After a few thousand years, the rate of change in inclination decreases, and after 20 000 yr the inclinations of those initial Jupiter family members that still have orbits with   P < 20 yr  become constant at about     , independent of whether non-gravitational forces are acting or not. The presently known Jupiter family of comets is losing members at the rate of one in every 67 yr. To maintain the family in equilibrium, Jupiter has to capture comets at a similar rate, and these captured comets have to be of low inclination to compensate for the pumping up of inclinations by gravitational perturbation.     

A Chandra X-ray observation of the globular cluster Terzan 1

We present a ∼19-ks Chandra Advanced CCD Imaging Spectrometer (ACIS)-S observation of the globular cluster Terzan 1. 14 sources are detected within 1.4 arcmin of the cluster centre with two of these sources predicted to be not associated with the cluster (background active galactic nuclei or foreground objects). The neutron star X-ray transient, X1732−304, has previously been observed in outburst within this globular cluster with the outburst seen to last for at least 12 yr. Here, we find four sources that are consistent with the ROSAT position for this transient, but none of the sources are fully consistent with the position of a radio source detected with the Very Large Array that is likely associated with the transient. The most likely candidate for the quiescent counterpart of the transient has a relatively soft spectrum and an unabsorbed 0.5–10 keV luminosity of  2.6 × 10³² erg s⁻¹  , quite typical of other quiescent neutron stars. Assuming standard core cooling, from the quiescent flux of this source we predict long (>400 yr) quiescent episodes to allow the neutron star to cool. Alternatively, enhanced core cooling processes are needed to cool down the core. However, if we do not detect the quiescent counterpart of the transient this gives an unabsorbed 0.5–10 keV luminosity upper limit of  8 × 10³¹ erg s⁻¹  . We also discuss other X-ray sources within the globular cluster. From the estimated stellar encounter rate of this cluster we find that the number of sources we detect is significantly higher than expected by the relationship of Pooley et al.     

Circinus X-1: survivor of a highly asymmetric supernova

We have analysed the kinematical parameters of Cir X-1 to constrain the nature of its companion star, the eccentricity of the binary and the pre-supernova parameter space. We argue that the companion is most likely to be a low-mass (≲2.0 M_⊙) unevolved star and that the eccentricity of the orbit is 0.94±0.04. We have evaluated the dynamical effects of the supernova explosion and we find it must have been asymmetric. On average , we find that a kick of ∼740 km s⁻¹ is needed to account for the recently measured radial velocity of +430 km s⁻¹ (Johnston, Fender & Wu) for this extreme system. The corresponding minimum kick velocity is ∼500 km s⁻¹. This is the largest kick needed to explain the motion of any observed binary system. If Cir X-1 is associated with the supernova remnant G321.9-0.3 then we find a limiting minimum age of this remnant of ∼60 000 yr. Furthermore, we predict that the companion star has lost ∼10 per cent of its mass as a result of stripping and ablation from the impact of the supernova shell shortly after the explosion.     

Simulations for the Original Distribution of KBOs

Using the N-body dynamical model that includes the sun, the 8 planets, Pluto, UB313 and massless particles, we simulate the orbital evolution of 551 Kuiper Belt Objects (KBOs) with known parameters. The initial conditions of the simulations are the currently observed orbital parameters. The integration backtracks from now to -10×10⁸ yr. The results show that about 10×10⁸ years ago, more than 1/3 of the presently observed KBOs resided in the region of the present Kuiper main belt, a few were located inside the Neptune orbit, and the rest were beyond 50AU; and that about 4.5×10⁸ years ago, all the objects in the Kuiper main belt exhibited a rather good normal distribution, without so many objects concentrated in the Neptune's 3:2 resonance region, as at present time.     

SGP4/SDP4模型精度分析

本文基于最新发布的SGP4/SDP4(Simplified General Perturbation Version 4/Simplified Deep-space Perturbation Version 4)模型设计了一套定轨方案,从空间目标库中挑选出不同类型和轨道参数的1120个目标进行计算,定量给出了SGP4/SDP4模型处理不同类型空间目标的定轨预报精度.结果表明:近地目标定轨精度为百米量级;半同步和同步轨道定轨精度平均为0.7和1.9km.椭圆轨道目标的定轨精度与偏心率有关,除少数e>0.8的椭圆轨道目标,绝大多数椭圆轨道目标定轨误差均小于10km.用SGP4/SDP4模型对近地目标预报3天,半同步轨道预报30天,同步轨道预报15天,椭圆轨道预报1天,预报误差一般不超过40km.     

A 408-MHz variability study of two fields in the region 6h30m, +45°

Variability in 408-MHz flux density, over a 1-yr interval between 1995 and 1996, has been investigated using the DRAO Synthesis Telescope for a complete sample of 322 sources with S _408 MHz > 80 mJy at galactic latitude b  = 15°; four sources are found to show significant flux density variations. These results, taken with those from a similar study at b  = 58°, indicate that 408-MHz variability over a time-scale of 1 yr is more uncommon than has previously been thought; variations at a level >10 per cent could have been detected in a statistically complete sample of about 150 of these sources and none was found. Variability over a period of 18 yr has also been investigated by comparing the DRAO data with those from the B3 catalogue for 182 sources with S _408 MHz > 200 mJy; five variables are found. Large flux density variations over a period of 18 yr are much more common than over 1 yr, with roughly one source in 15 showing variations at the 15 per cent level or greater. The statistics are not good enough to show any conclusive evidence for a dependence of variability on galactic latitude. Roughly half of the variable sources have spectra that are steep at high frequencies and flatten at low frequencies; the other half have flat radio spectra over the entire frequency range. Simultaneous observations at 1420 MHz, also with the DRAO Synthesis Telescope, have been used to investigate variability at this frequency as well, and indicate that over a period of 1 yr roughly half of the flat-spectrum sources vary at the 5 per cent level or greater.     

Eclipsing phenomena of the symbiotic star CH Cygni

Eclipsing phenomena of the inner binary with a period of 756 d in the triple symbiotic system CH Cyg may have been detected in detailed spectrophotometric observations. The eclipse of the hot component by the red giant started on 1994 October 11 and finished between 1995 January 8 and 18. The ingress duration was less than one day. The radius of the red giant is estimated to be 288 ± 15 R⊙ from the duration of the eclipse. Assuming the bolometric correction of the red giant (M7 III) as 4, the distance to this object is estimated to be 307 ± 32 pc, which agrees well with that obtained in the observations by Hipparcos . The interstellar extinction in the direction of this object may be much lower than that in the nearby areas. It has been suggested that the outer binary system with an orbital period of about 15 yr is an eclipsing one. It seems unlikely, however, that the variation of the activity of this object with a time-scale of more than 10 yr was a result of eclipses.     

Monte Carlo simulations of star clusters — I. First Results

A revision of Stodoíkiewicz's Monte Carlo code is used to simulate evolution of star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The first calculations for isolated, equal-mass N -body systems with three-body energy generation according to Spitzer's formulae show good agreement with direct N -body calculations for N  = 2000, 4096 and 10 000 particles. The density, velocity, mass distributions, energy generation, number of binaries, etc., follow the N -body results. Only the number of escapers is slightly too high compared with N -body results, and there is no level-off anisotropy for advanced post-collapse evolution of Monte Carlo models as is seen in N -body simulations for N  ≤ 2000. For simulations with N  > 10 000 gravothermal oscillations are clearly visible. The calculations of N   2000, 4096, 10 000, 32 000 and 100 000 models take about 2, 6, 20, 130 and 2500 h, respectively. The Monte Carlo code is at least 10⁵ times faster than the N -body one for N  = 32 768 with special-purpose hardware. Thus it becomes possible to run several different models to improve statistical quality of the data and run individual models with N as large as 100 000. The Monte Carlo scheme can be regarded as a method which lies in the middle between direct N -body and Fokker–Planck models and combines most advantages of both methods.     

Orbital and thermal evolutions of four potential targets for a sample return space mission to a primitive near-Earth asteroid

In this paper, we present our study of the orbital and thermal evolutions, due to solar radiative heating, of four near-Earth asteroids (NEAs) considered as potential target candidates for sample return space missions to primitive asteroids. We used a dynamical model of the NEA population to estimate the most likely source region and orbital history of these objects. Then, for each asteroid, we integrated numerically over their entire lifetime a set of 14 initially indistinguishable orbit (clones), obtained by small variations of the nominal initial conditions. Using a thermal model, we then computed surface and sub-surface temperatures of these bodies during their dynamical history. Our aim is to determine whether these bodies are likely to have experienced high temperature level, and whether great temperature changes can be expected due to the orbital changes as well as their maximum and minimum values. Such information is important in the framework of sample return space missions whose goal is to bring back pristine materials. The knowledge of the temperature range of materials at different depth over the orbital evolution of potential targets can help defining sampling strategies that ensure the likelihood that unaltered material will be brought back. Our results suggest that for all the considered potential targets, the surface has experienced for some time temperatures greater than 400 K and at most 500 K with 50% probability. This probability drops rapidly with increasing temperature. Sub-surface materials at a depth of only 3 cm are much more protected from high temperature and generally do not reach temperatures exceeding 450 K (with 50% probability). They should thus be unaltered at this depth at least from a Sun-driven heating point of view. On the other hand, surface material for some of the considered objects can have a range of temperature which can make them less reliable as pristine materials. However, it is assumed here that the same material is constantly exposed to solar heat, while regolith turnover may occur. The latter can be caused by different processes such as seismic shaking and/or impact cratering. This would reduce the total time that materials are exposed to a certain temperature. Thus, it is very likely that a sample collected from any of the four considered targets, or any primitive NEA with similar dynamical properties, will have components that will be thermally unaltered as long as some of it comes from only 3 to 5 cm depth. Such a depth is not considered difficult to reach with some of the current designs of sampling devices.     

Variability of Seyfert galaxies in the near-infrared — II. The unusual starburst Seyfert NGC 7469

The bright Seyfert 1 and starburst galaxy NGC 7469 has been monitored at JHKL for 13 yr. Its variations have been small, amounting to less than 0.3 mag, except for an event in 1989 when the flux from its active nucleus declined to about one quarter of its normal level.   Using the technique of flux variation gradients (FVGs), the variations have been analysed, taking into account the known JHK fluxes from the star-forming ring surrounding the nucleus. By comparing the FVGs of NGC 7469 with those of other Seyfert galaxies, the nucleus of NGC 7469 is found to suffer very little reddening.   The active nucleus normally contributes more than 50 per cent of the flux measured in the K and L bands through an aperture of 12 arcsec diameter.   Infrared excesses (above the fluxes expected from stars) are observed from the circumnuclear ring at both K and L and it is suggested that they arise from warm, small particles that also emit the aromatic hydrocarbon features seen in this galaxy.     

A detailed abundance analysis of the hot post-AGB star ZNG-1 in M10

We present a model-atmosphere analysis for the bright ( V ∼13) star ZNG-1, in the globular cluster M10. From high-resolution ( R ∼40 000) optical spectra we confirm ZNG-1 to be a post-asymptotic giant branch (post-AGB) star. The derived atmospheric parameters are T _eff=26 500±1000 K and log  g =3.6±0.2 dex . A differential abundance analysis reveals a chemical composition typical of hot post-AGB objects, with ZNG-1 being generally metal poor, although helium is approximately solar. The most interesting feature is the large carbon underabundance of more than 1.3 dex. This carbon deficiency, along with an observed nitrogen enhancement relative to other elements, may suggest that ZNG-1 evolved off the AGB before the third dredge-up occurred. Also, iron depletions observed in other similar stars suggest that gas–dust fractionation in the AGB progenitor could be responsible for the observed composition of these objects. However, we need not invoke either scenario since the chemical composition of ZNG-1 is in good agreement with abundances found for a Population II star of the same metallicity.     

The long-term evolution and initial size of comets 46P/Wirtanen and 67P/Churyumov–Gerasimenko

We present a new method to study the long-term evolution of cometary nuclei in order to estimate their original size, and we consider the case of comets 46P/Wirtanen (hereafter 46P) and 67P/Churyumov–Gerasimenko (hereafter 67P). We calculate the past evolution of the orbital elements of both comets over 100 000 yr using a Bulirsch–Stoer integrator and over 450 000 yr using a Radau integrator, and we incorporate a realistic model of the erosion of their nucleus. Their long-term orbital evolution is prominently chaotic, resulting from several close encounters with planets, and this result is independent of the choice of the integrator and of the presence or not of non-gravitational forces. The dynamical lifetime of comet 46P is estimated at ∼133 000 yr and that of comet 67P at ∼105 000 yr. Our erosion model assumes a spherical nucleus composed of a macroscopic mixture of two thermally decoupled components, dust and pure water ice. Erosion strongly depends upon the active fraction and the density of the nucleus. It mainly takes place at heliocentric distances <4 au and lasts for only ∼7 per cent of the lifetime. Assuming a density of 300 kg m⁻³ and an average active fraction over time of 10 per cent, we find an initial radius of ∼1.3 km for 46P and ∼2.8 km for 67P. Upper limit are obtained assuming a density of 100 kg m⁻³ and an active fraction of 100 per cent, and amounts to 21 km for 46P and 25 km for 67P. Erosion acts as a rejuvenating process of the surface so that exposed materials on the surface may only contain very little quantities of primordial materials. However, materials located just under it (a few centimetres to metres) may still be much less evolved. We will apply this method to several other comets in the future.     

Dynamics of rotationally fissioned asteroids: Source of observed small asteroid systems

We present a model of near-Earth asteroid (NEA) rotational fission and ensuing dynamics that describes the creation of synchronous binaries and all other observed NEA systems including: doubly synchronous binaries, high-e binaries, ternary systems, and contact binaries. Our model only presupposes the Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) effect, “rubble pile” asteroid geophysics, and gravitational interactions. The YORP effect torques a “rubble pile” asteroid until the asteroid reaches its fission spin limit and the components enter orbit about each other (Scheeres, D.J. [2007]. Icarus 189, 370-385). Non-spherical gravitational potentials couple the spin states to the orbit state and chaotically drive the system towards the observed asteroid classes along two evolutionary tracks primarily distinguished by mass ratio. Related to this is a new binary process termed secondary fission - the secondary asteroid of the binary system is rotationally accelerated via gravitational torques until it fissions, thus creating a chaotic ternary system. The initially chaotic binary can be stabilized to create a synchronous binary by components of the fissioned secondary asteroid impacting the primary asteroid, solar gravitational perturbations, and mutual body tides. These results emphasize the importance of the initial component size distribution and configuration within the parent asteroid. NEAs may go through multiple binary cycles and many YORP-induced rotational fissions during their approximately 10 Myr lifetime in the inner Solar System. Rotational fission and the ensuing dynamics are responsible for all NEA systems including the most commonly observed synchronous binaries.     

A survey of orbits of co-orbitals of Mars

Many asteroids with a semimajor axis close to that of Mars have been discovered in the last several years. Potentially some of these could be in 1:1 resonance with Mars, much as are the classic Trojan asteroids with Jupiter, and its lesser-known horseshoe companions with Earth. In the 1990s, two Trojan companions of Mars, 5261 Eureka and 1998 VF₃₁, were discovered, librating about the L₅ Lagrange point, 60° behind Mars in its orbit. Although several other potential Mars Trojans have been identified, our orbital calculations show only one other known asteroid, 1999 UJ₇, to be a Trojan, associated with the L₄ Lagrange point, 60° ahead of Mars in its orbit. We further find that asteroid 36017 (1999 ND₄₃) is a horseshoe librator, alternating with periods of Trojan motion. This asteroid makes repeated close approaches to Earth and has a chaotic orbit whose behavior can be confidently predicted for less than 3000 years. We identify two objects, 2001 HW₁₅ and 2000 TG₂, within the resonant region capable of undergoing what we designate “circulation transition”, in which objects can pass between circulation outside the orbit of Mars and circulation inside it, or vice versa. The eccentricity of the orbit of Mars appears to play an important role in circulation transition and in horseshoe motion. Based on the orbits and on spectroscopic data, the Trojan asteroids of Mars may be primordial bodies, while some co-orbital bodies may be in a temporary state of motion.     

Searching for stars closely encountering with the solar system

Based on a new version of the Hipparcos catalog and currently available radial velocity data, we have searched for stars that either have encountered or will encounter the solar neighborhood within less than 3 pc in the time interval from −2 Myr to +2 Myr. Nine new candidates within 30 pc of the Sun have been found. To construct the stellar orbits relative to the solar orbit, we have used the epicyclic approximation. We show that, given the errors in the observational data, the probability that the well-known star HIP 89 825 (GL 710) encountering with the Sun most closely falls into the Oort cloud is 0.86 in the time interval 1.45 ± 0.06 Myr. This star also has a nonzero probability, 1 × 10⁻⁴, of falling into the region d < 1000 AU, where its influence on Kuiper Belt objects becomes possible.     

The X-ray spectrum and variability of the Seyfert 2 galaxy NGC 7172

We present evidence of flux variability, on both short (hours) and long (months) time-scales, of the Seyfert 2 galaxy NGC 7172. These results are based on the ASCA observation of NGC 7172 performed in 1996 May. The source was detected at a rather low flux level, about 3 times fainter than its usual state (including 1 yr before, when it was also observed by ASCA ).   The source also varied by about 30 per cent during the observation, confirming the presence of a type 1 nucleus in its centre. However, its spectrum appears to be flatter than the typical Seyfert 1 spectrum (in agreement with findings on other Seyfert 2s), posing problems for the unification model unless complex absorption is invoked.