The Durham/ESO SLODAR optical turbulence profiler

An instrument for monitoring of the vertical profile of atmospheric optical turbulence strength, employing the Slope Detection and Ranging (SLODAR) double star technique applied to a small telescope, has been developed by Durham University and the European South Observatory. The system has been deployed at the Cerro Paranal observatory in Chile for statistical characterization of the site. The instrument is configured to sample the turbulence at altitudes below 1.5 km with a vertical resolution of approximately 170 m. The system also functions as a general-purpose seeing monitor, measuring the integrated optical turbulence strength for the whole atmosphere, and hence the seeing width. We give technical details of the prototype and present data to characterize its performance. Comparisons with contemporaneous measurements from a differential image motion monitor (DIMM) and a multi-aperture scintillation sensor (MASS) are discussed. Statistical results for the optical turbulence profile at the Paranal site are presented. We find that, in the median case, 49 per cent of the total optical turbulence strength is associated with the surface layer (below 100 m), 35 per cent with the 'free atmosphere' (above 1500 m) and 16 per cent with the intermediate altitudes (100–1500 m).     

Large-amplitude periodic variations in the angular diameter of R Leonis

We report the first direct detection of long-term periodic diameter variations in a Mira variable. Angular diameter measurements of the 313-d period variable R Leonis at 833 nm and 940 nm obtained between 1996 February and 1997 June using the Cambridge Optical Aperture Synthesis Telescope (COAST) and the William Herschel Telescope (WHT) show a cyclic modulation of the apparent stellar diameter by approximately 35 per cent. The agreement between these new data and archival measurements from 1992 January suggests coherence in the modulation over a 5-yr period. Our data are consistent with recent models which suggest that, in photometric bands with only weak to moderate molecular contamination, periodic variations in stellar diameter of order 50 per cent can be maintained. The measurements indicate that the apparent stellar diameter was largest at visual phase 0.5 and that any phase shifts between the visual light curve and those at 833 and 940 nm were at most 0.05. The large offset (∼ 0.25) between the phase of the observed diameter maximum and that predicted for the photospheric continuum diameter variations suggests that our observations are more sensitive to the changing temperature structure of the outer atmosphere than to the deeper continuum-forming layers.     

High-altitude wind velocity at Oukaimeden observatory

This paper aims at studying the wind at 200 mbar over the Moroccan observatory Oukaimeden, as high-altitude winds have been adopted as a useful parameter for site characterization in terms of the suitability of a site for the development of some adaptive optics techniques. The data used come from the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis data base, which is widely acknowledged as being reliable. Statistical analyses of 200-mbar wind speed since 1983 are performed. Comparison with some of the main observatory sites worldwide qualifies Oukaimeden as one of the best observatory sites in terms of 200-mbar wind statistics. Our analysis of a record of seeing measurements during the years 2003 and 2004 concludes that while 200-mbar wind speed can be used as a parameter for ranking astronomical sites in term of their suitability for adaptive optics, it cannot be used for the whole atmospheric seeing prediction. A comparison of monthly values of the seeing parameter at Oukaimeden, La Silla and Paranal demonstrates the high seeing quality of Oukaimeden, as the seeing values measured were lower than those of La Silla and Paranal for most of the time during the comparison period. Furthermore, a statistical analysis of atmospheric stratified seeing, wavefront coherence time and isoplanatic angle measured with a Multi-Aperture Scintillation Sensor instrument over Paranal from 2004 to 2007 have been performed. We found good correlations between 200-mbar wind velocity and levels 4, 5 and 6 seeing, wavefront coherence time and isoplanatic angle, with corresponding correlation coefficients of 0.74, 0.79, 0.70, 0.97 and 0.78.     

Monitoring Secular Orbital Period Variations of Some Eclipsing Binaries at the Ankara University Observatory

This study concerns the long-term monitoring of the secular variation character in the orbital period of some short-period eclipsing binaries observed at the Ankara University Observatory. Among the systems of our observing list are CK Boo, V502 Oph and V836 Cyg that show long-term secular variations in their orbital periods. We use classical O-C diagram analysis technique as a tool to reveal the character of the period variations of these binary systems.     

中国最早的观象台发掘

在山西襄汾县陶寺镇附近发掘出距今约4000年的最早观象台兼祭祀台遗址,它是由观测点、夯土柱和柱间狭缝组成,用于观测日出方位变化,确定回归年长度来制订历法。模拟观测结果表明,位于东南和东北的二狭缝可准确测定冬至和夏至日期,而此二狭缝之间有10个土柱,应象征着视太阳向北和向南每位移一个土柱为一个节气。由此推之,帝尧时的历法特征是将一岁分成20个节气的阳历,上古时的阴阳五行历即十月太阳历是源于它奠定的基础。     

The exceptional X-ray variability of the dwarf Seyfert nucleus NGC 4395

An analysis of the X-ray variability of the low-luminosity Seyfert nucleus NGC 4395, based on a long XMM–Newton observation, is presented. The power spectrum shows a clear break from a flat spectrum  (α≈ 1)  to a steeper spectrum  (α≈ 2)  at a frequency   f _br= 0.5–3.0 × 10⁻³ Hz  , comparable to the highest characteristic frequency found previously in a Seyfert galaxy. This extends the measured   M _BH− f _br  values to lower M _BH than previous studies of Seyfert galaxies, and is consistent with an inverse scaling of variability frequency with black hole mass. The variations observed are among the most violent seen in an active galactic nuclei to date, with the fractional rms amplitude  ( F _var)  exceeding 100 per cent in the softest band. The amplitude of the variations seems intrinsically higher in NGC 4395 than most other Seyfert galaxies, even after accounting for the differences in characteristic frequencies. The origin of this difference is not clear, but it is unlikely to be a high accretion rate (   L / L _Edd≲ 20  per cent for NGC 4395). The variations clearly follow the linear rms–flux relation, further supporting the idea that this is a ubiquitous characteristics of accreting black holes. The variations are highly coherent between different energy bands with any frequency-dependent time delay limited to ≲1 per cent.     

"陶寺观象台"发掘与古代文献的印证

古代文献记载帝尧时的历法存在两种不同说法。羲和观测四仲中星定季节当是不实之辞,并以闰月成岁推论为使用农历也没依据。羲和观测太阳出没方位变化而"生十日"、即制订产生十月太阳历才符合事实。在文[1]中的推论,也与后一种说法可互相印证。本文对文[1]的推论作出了论据补充,与文[1]成姐妹篇。     

Power spectra of cosmic ray intensity variations at the ground and atmospheric pressure oscillations

Power spectrum estimates are computed for cosmic ray and pressure variations in the frequency range of 1.6 × 10⁻⁶ to 4.15 × 10⁻³ c/sec for three data sets each of 27 days length (2 min interval) recorded at Chacaltaya (Bolivia) during 1965–1966. The general trend of these spectra showy ƒ⁻ fit having exponent values from −1.5 to −1.9 for L.F. side. From the semi-diurnal peaks both in cosmic rays and pressure, the average value of pressure coefficient is found to be 0.3 per cent/mm of Hg. For the same three sets of data, a detailed analysis of cross-spectrum, coherence and phase relationship between cosmic rays and pressure is carried out in the frequency range of 1–12 c/day. Besides semi-diurnal peak, variations of 3, 4, 7 and 8 c/day are found to be common both in cosmic rays and pressure. The value of pressure coefficient and residual amplitude of cosmic rays for these particular frequencies are calculated.     

Variability and spectral energy distributions of low-luminosity active galactic nuclei: a simultaneous X-ray/UV look with Swift

We have observed four low-luminosity active galactic nuclei (AGNs) classified as type 1 Low-Ionization Nuclear Emission-Line Regions (LINERs) with the X-Ray Telescope (XRT) and the Ultraviolet–Optical Telescope (UVOT) onboard Swift , in an attempt to clarify the main powering mechanism of this class of nearby sources. Among our targets, we detect X-ray variability in NGC 3998 for the first time. The light curves of this object reveal variations of up to 30 per cent amplitude in half a day, with no significant spectral variability on this time-scale. We also observe a decrease of ∼30 per cent over 9 d, with significant spectral softening. Moreover, the X-ray flux is ∼40 per cent lower than observed in previous years. Variability is detected in M81 as well, at levels comparable to those reported previously: a flux increase in the hard X-rays (1–10 keV) of 30 per cent in ∼3 h and variations by up to a factor of 2 within a few years. This X-ray behaviour is similar to that of higher luminosity, Seyfert-type objects. Using previous high-angular-resolution imaging data from the Hubble Space Telescope ( HST ), we evaluate the diffuse UV emission due to the host galaxy and isolate the nuclear flux in our UVOT observations. All sources are detected in the UV band, at levels similar to those of the previous observations with HST . The XRT (0.2–10 keV) spectra are well described by single power laws and the UV-to-X-ray flux ratios are again consistent with those of Seyferts and radio-loud AGNs of higher luminosity. The similarity in X-ray variability and broad-band energy distributions suggests the presence of similar accretion and radiation processes in low- and high-luminosity AGNs.     

Measurements of Atmospheric Extinction at a Ground Level Observatory

In order to determine the atmospheric extinction at Physics Department Astronomical Observatory of the University of Extremadura, located at Badajoz (Spain), several stars were observed during some clear nights of atmospheric stability in the period 1998–2000, at optical wavelengths corresponding to the filters of the Kron-Cousins and Strömgren photometric systems. The determination of the extinction coefficients was made by assuming the Bouguer's law, which was shown to be a good approximation for this study. The results exhibited temporal variations and can be considered to be associated with clean atmospheres at locations of low altitude.     

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.     

The spectrum of electron content fluctuations in the ionosphere

Continuous records of the electron content of the ionosphere, from 1965 to 1970, are used to obtain power spectra covering periods from 30 sec to 2 yr at latitudes of 34°S and 42°S. At periods up to 5 min, amplitudes were less than 0.2 per cent of the total electron content. Variations produced by gravity waves were very common in the range 20–80 min, with no preferred periods. The r.m.s. amplitude per octave A₀ was about 10¹⁵ electrons/m², or 0.6 per cent of the mean electron content. The amplitude increased during the day, particularly in winter when periodic components predominated. The cut-off at about 17 min was sharply defined, giving a mean scale height for the neutral atmosphere (at 300 km) of about 43 km in summer, 47 km on winter days and 42 km on winter nights.

From 12 hr to 1 month A₀ was about 12 per cent of the mean electron content in both summer and winter at 34°S, and 10 per cent at 42°S. The 24 hr and 27 day peaks were largest just before sunspot maximum, and almost disappeared near sunspot minimum. Variations between 1 and 27 days reflect the random occurrence of ionospheric storms and show no consistent peaks. Day to day and night to night variations were both about 10 per cent of the background content for periods from 2 days to 2 yr, apart from a slight decrease between 1 and 6 months.     

Adaptive optics parameters connection to wind speed at the Teide Observatory

Current projects for large telescopes demand a proper knowledge of atmospheric turbulence to design efficient adaptive optics systems in order to reach large Strehl ratios. However, the proper characterization of the turbulence above a particular site requires long-term monitoring. Because of the lack of long-term information on turbulence, high-altitude winds (in particular winds at the 200 mbar pressure level) were proposed as a parameter for estimating the total turbulence at a particular site, with the advantage of records of winds going back several decades. We present the first complete study of atmospheric adaptive optics parameters above the Teide Observatory (Canary Islands, Spain) in relation to wind speed. On-site measurements of   C ² _N ( h )  profiles (more than 20 200 turbulence profiles) from G-SCIDAR (Generalized Scintillation Detection and Ranging) observations and wind vertical profiles from balloons have been used to calculate the seeing, the isoplanatic angle and the coherence time. The connection of these parameters to wind speeds at ground and at 200 mbar pressure level are shown and discussed. Our results confirm the well-known high quality of the Canary Islands astronomical observatories.     

The new intermediate polar RX J1238 − 38: a system below the period gap?

We present optical observations of the recently discovered ROSAT source RX J1238 − 38, which is a new member of the intermediate polar class of asynchronous magnetic cataclysmic variables (CVs). Optical photometry reveals two coherent periodicities at 1860 and 2147 s respectively, with similar amplitudes of ∼ 8 per cent. Infrared ( J -band) intensity variations are detected only at the 1860-s period, at an amplitude of ∼ 15 per cent. The initial hypothesis, that these two periods were the spin and synodic (i.e., beat) period respectively, appears not to be supported by the spectroscopic data. The emission lines vary on the longer photometric period, and radial velocity variations are detected at this period and at a longer period of ∼ 5300 s, which we identify as the spin and orbital periods respectively. The most likely explanation for the 1860-s period is that it is the first harmonic of the ω − Ω sideband, leading to an improved determination of the orbital period as 5077 s (= 84 min). If this interpretation is correct, RX J1238 − 38 joins EX Hya as the only other intermediate polar below the 2–3 h period gap, and with an orbital period close to the minimum for CVs with non-degenerate secondaries. The spin-modulated emission-line radial velocities and widths appear to be anticorrelated, with maximum width occurring at maximum blueshift. Such an anticorrelation is expected for aspect changes of accretion curtains. Polarimetric observations of RX J1238 − 38 were inconclusive, although we can put a limit of 0.4 per cent on any variability on the circular polarization, and certainly there is no indication of variations at the photometric or spectroscopic periods.     

Be star candidates in the Large Magellanic Cloud: the catalogue and comparison with the Small Magellanic Cloud sample

We present a catalogue with coordinates and photometric data of 2446 Be star candidates in the Large Magellanic Cloud (LMC), based on a search of the OGLE II data base. The I -band light curves of these stars show outbursts in 24 per cent of the sample (Type-1 stars), high and low states in 10 per cent, periodic variations in 6 per cent (Type-3 stars), and stochastic variations in 60 per cent of the cases. We report on the result of the statistical study of light curves of Type-1 and Type-3 stars in the LMC, and the comparison with the previously reported results of the Small Magellanic Cloud (SMC) sample. We find a statistically significant difference between amplitude, duration and asymmetry distributions of outbursts in both galaxies. Outbursts of SMC Type-1 stars are usually brighter, longer and with a slower decline. We find a bimodal distribution of periods of Type-3 stars in both galaxies, probably related to the recently discovered double periodic blue variables. We find also period and amplitude distributions of Type-3 LMC stars statistically different from those of the SMC stars. Our findings above suggest that the mechanisms causing the observed photometric variability of Type-1 and Type-3 stars could depend on metallicity. Moreover, they suggest that the outbursts are not primarily caused by stellar winds.     

Compact radio sources variable at 151 MHz – I. Catalogue of sources

The Cambridge Low-Frequency Synthesis Telescope has been used to produce a representative sample of low-frequency variable sources. 20 fields, each covering an area of approximately 9°×9° cosec  δ , have been observed at 151 MHz at between 2 and 10 epochs over the period from 1984 to 1996. At each epoch, maps were made with rms noise levels of typically 10–15 mJy beam⁻¹. From a total of ∼6000 sources detected on these maps, 207 are found for which the flux density variations between at least two epochs appear significant at greater than the 3 σ level. A numerical model is used to assess the true significance of the variability, given the analysis method adopted. This shows that for about half of the sources which appear to vary by >3 σ the variability is genuine. For the other half it is caused by random statistical fluctuations; most of the spurious variables vary by ≲3.5 σ between a single pair of epochs. A catalogue of the variable sources is presented, which includes an estimate of the probability that a given source is a genuine variable. Fractional flux density variations of between 5 and 100 per cent (typically 15–25 per cent) have been detected on a range of time-scales from 1 to 12 years.     

Numerical simulations of the Lyman α forest – a comparison of gadget-2 and enzo

We compare simulations of the Lyman α forest performed with two different hydrodynamical codes, gadget-2 and enzo . A comparison of the dark matter power spectrum for simulations run with identical initial conditions show differences of 1–3 per cent at the scales relevant for quantitative studies of the Lyman α forest. This allows a meaningful comparison of the effect of the different implementations of the hydrodynamic part of the two codes. Using the same cooling and heating algorithm in both codes, the differences in the temperature and the density probability distribution function are of the order of 10 per cent. The differences are comparable to the effects of box size and resolution on these statistics. When self-converged results for each code are taken into account, the differences in the flux power spectrum – the statistics most widely used for estimating the matter power spectrum and cosmological parameters from Lyman α forest data – are about 5 per cent. This is again comparable to the effects of box size and resolution. Numerical uncertainties due to a particular implementation of solving the hydrodynamic or gravitational equations appear therefore to contribute only moderately to the error budget in estimates of the flux power spectrum from numerical simulations. We further find that the differences in the flux power spectrum for enzo simulations run with and without adaptive mesh refinement are also of the order of 5 per cent or smaller. The latter require 10 times less CPU time making the CPU time requirement similar to that of a version of gadget-2 that is optimized for Lyman α forest simulations.     

Constraints on turbulent pressure in the X-ray haloes of giant elliptical galaxies from resonant scattering

The dense cores of X-ray emitting gaseous haloes of large elliptical galaxies with temperatures   kT ≲ 0.8  keV  show two prominent Fe  xvii emission features, which provide a sensitive diagnostic tool to measure the effects of resonant scattering. We present here high-resolution spectra of five bright nearby elliptical galaxies, obtained with the reflection grating spectrometers (RGS) on the XMM-Newton satellite. The spectra for the cores of four of the galaxies show the Fe  xvii line at 15.01 Å being suppressed by resonant scattering. The data for NGC 4636 in particular allow the effects of resonant scattering to be studied in detail and to prove that the 15.01 Å line is suppressed only in the dense core and not in the surrounding regions. Using deprojected density and temperature profiles for this galaxy obtained with the Chandra satellite, we model the radial intensity profiles of the strongest resonance lines, accounting for the effects of resonant scattering, for different values of the characteristic turbulent velocity. Comparing the model to the data, we find that the isotropic turbulent velocities on spatial scales smaller than ≈1 kpc are less than 100 km s⁻¹ and the turbulent pressure support in the galaxy core is smaller than 5 per cent of the thermal pressure at the 90 per cent confidence level, and less than 20 per cent at 95 per cent confidence. Neglecting the effects of resonant scattering in spectral fitting of the inner 2 kpc core of NGC 4636 will lead to underestimates of the chemical abundances of Fe and O by ∼10–20 per cent.     

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.