We used HARPS to measure oscillations in the low-mass star tau Cet. Although the data were compromised by instrumental noise, we have been able to extract the main features of the oscillations. We found tau Cet to oscillate with an amplitude that is about half that of the Sun, and with a mode lifetime that is slightly shorter than solar. The large frequency separation is 169 muHz, and we have identified modes with degrees 0, 1, 2, and 3. We used the frequencies to estimate the mean density of the star to an accuracy of 0.45% which, combined with the interferometric radius, gives a mass of 0.783 +/- 0.012 M_sun (1.6%).

Recent revisions of the determination of the solar composition have resulted in solar models in marked disagreement with helioseismic inferences. The effect of the composition change on the model is largely caused by the change in the opacity. Thus we wish to determine an intrinsic opacity change that would compensate for the revision of the composition. By comparing models computed with the old and revised composition we determine the required opacity change. Models are computed with the opacity thus modified and used as reference in helioseismic inversions to determine the difference between the solar and model sound speed. An opacity increase varying from around 30 per cent near the base of the convection zone to a few percent in the solar core results in a sound-speed profile, with the revised composition, which is essentially indistinguishable from the original solar model. As a function of the logarithm of temperature this is well represented by a simple cubic fit. The physical realism of such a change remains debatable, however.

Asteroseismology provides us with a unique opportunity to improve our understanding of stellar structure and evolution. Recent developments, including the first systematic studies of solar-like pulsators, have boosted the impact of this field of research within Astrophysics and have led to a significant increase in the size of the research community. In the present paper we start by reviewing the basic observational and theoretical properties of classical and solar-like pulsators and present results from some of the most recent and outstanding studies of these stars. We centre our review on those classes of pulsators for which interferometric studies are expected to provide a significant input. We discuss current limitations to asteroseismic studies, including difficulties in mode identification and in the accurate determination of global parameters of pulsating stars, and, after a brief review of those aspects of interferometry that are most relevant in this context, anticipate how interferometric observations may contribute to overcome these limitations. Moreover, we present results of recent pilot studies of pulsating stars involving both asteroseismic and interferometric constraints and look into the future, summarizing ongoing efforts concerning the development of future instruments and satellite missions which are expected to have an impact in this field of research.

We have observed oscillations in the nearby G2 subgiant star beta Hyi using high-precision velocity observations obtained over more than a week with the HARPS and UCLES spectrographs. The oscillation frequencies show a regular comb structure, as expected for solar-like oscillations, but with several l=1 modes being strongly affected by avoided crossings. The data, combined with those we obtained five years earlier, allow us to identify 28 oscillation modes. By scaling the large frequency separation from the Sun, we measure the mean density of beta Hyi to an accuracy of 0.6%. The amplitudes of the oscillations are about 2.5 times solar and the mode lifetime is 2.3 d. A detailed comparison of the mixed l=1 modes with theoretical models should allow a precise estimate of the age of the star.

Scientific output varies between research fields and between disciplines within a field such as astrophysics. Even in fields where publication is the primary output, there is considerable variation in publication and hence in citation rates. Data from the Smithsonian/NASA Astrophysics Data System is used to illustrate this problem and argue against a "one size fits all" approach to performance metrics, especially over the short time-span covered by the Research Assessment Exercise (soon underway in the UK).

A central step in time-distance local helioseismology techniques is to obtain travel times of packets of wave signals between points or sets of points on the visible surface. Standard ways of determining group or phase travel times involve cross-correlating the signal between locations at the solar photosphere and determining the shift of the envelope of this cross correlation function, or a zero crossing, using a standard wavelet or a reference wave packet. Here a novel method is described which makes use of triple correlations, i.e. cross-correlating signal between three locations. By using an average triple correlation as reference, differential travel times can be extracted in a straightforward manner.

A method is presented for investigating the periodic signal content of time series in which a number of signals is present, such as arising from the observation of multiperiodic oscillating stars in observational asteroseismology. Standard Fourier analysis tends only to be effective in cases when the data are perfectly regularly sampled. During normal telescope operation it is often the case that there are large, diurnal, gaps in the data, that data are missing, or that the data are not regularly sampled at all. For this reason it is advantageous to perform the analysis as much as possible in the time domain. Furthermore, for quantitative analyses of the frequency content and power of all real signals, it is of importance to have good estimates of the errors on these parameters. This is easiest to perform if one can use linear combinations of the measurements. Here such a linear method is described. The method is based in part on well-known techniques in radio technology used in every FM radio receiver, and in part on the SOLA inverse method

A statistical analysis is presented of the turbulent velocity structure in the Orion Molecular Cloud at scales ranging from 70 AU to 30000 AU. Results are based on IR Fabry-Perot interferometric observations of shock and photon-excited H2 in the K-band S(1) v=1-0 line at 2.121micron and refer to the dynamical characteristics of warm perturbed gas. Observations establish that the Larson size-linewidth relation is obeyed to the smallest scales studied here extending the range of validity of this relationship by nearly 2 orders of magnitude. The velocity probability distribution function (PDF) is constructed showing extended exponential wings, providing evidence of intermittency, further supported by the skewness and kurtosis of the velocity distribution. Variance and kurtosis of the PDF of velocity differences are constructed as a function of lag. The variance shows an approximate power law dependence on lag, with exponent significantly lower than the Kolmogorov value, and with deviations below 2000AU which are attributed to outflows and possibly disk structures associated with low mass star formation within OMC1. The kurtosis shows strong deviation from a gaussian velocity field, providing evidence of velocity correlations at small lags. Results agree accurately with semi-empirical simulations in Eggers & Wang (1998). In addition, 170 individual H2 emitting clumps have been analysed with sizes between 500 and 2200 AU. These show considerable diversity with regard to PDFs and variance functions. Our analysis constitutes the first characterization of the turbulent velocity field at the scale of star formation and provide a dataset which models of star-forming regions should aim to reproduce.

We report here on serendipitous observations of the intense gamma-ray flare from SGR 1806-20 that occured on 27 December 2004. Unique data from the Cluster and Double Star-2 satellites, designed to study the Earth's magnetosphere, provide the first observational evidence of three separate timescales within the early (first 100ms) phases of this class of events. These observations reveal that, in addition to the initial very steep (<0.25ms) X-ray onset, there is firstly a 4.9ms exponential rise timescale followed by a continued exponential rise in intensity on a timescale of 70ms. These three timescales are a prominent feature of current theoretical models including the timescale (several ms) for fracture propagation in the crust of the neutron star.

Various dedicated satellite projects are underway or in advanced stages of planning to perform high-precision, long duration time series photometry of stars, with the purpose of using the frequencies of stellar oscillations to put new constraints on the internal structure of stars. It is known (cf. Brown, et al. 1994) that the effectiveness of oscillation frequencies in constraining stellar model parameters is significantly higher if classical parameters such as effective temperature, and luminosity are known with high precision. In order to optimize asteroseismic campaigns it is therefore useful to select targets from among candidates for which good spectroscopic and astrometric data already exists. This paper presents selection criteria, as well as redeterminations of stellar luminosity and reddening for stars satisfying these criteria.

We report the firm discovery of solar-like oscillations in a giant star. We monitored the star xi Hya (G7III) continuously during one month with the CORALIE spectrograph attached to the 1.2m Swiss Euler telescope. The 433 high-precision radial-velocity measurements clearly reveal multiple oscillation frequencies in the range 50 - 130 uHz, corresponding to periods between 2.0 and 5.5 hours. The amplitudes of the strongest modes are slightly smaller than 2 m/s. Current model calculations are compatible with the detected modes.

We present optical lightcurves of the gravitationally lensed components A (=A1+A2+A3) and B of the quadruple quasar RX J0911.4+0551 (z = 2.80). The observations were primarily obtained at the Nordic Optical Telescope between 1997 March and 2001 April and consist of 74 I-band data points for each component. The data allow the measurement of a time delay of 146 +- 8 days (2 sigma) between A and B, with B as the leading component. This value is significantly shorter than that predicted from simple models and indicates a very large external shear. Mass models including the main lens galaxy and the surrounding massive cluster of galaxies at z = 0.77, responsible for the external shear, yield H_0 = 71 +- 4 (random, 2 sigma) +- 8 (systematic) km/s/Mpc. The systematic model uncertainty is governed by the surface-mass density (convergence) at the location of the multiple images.

We present optical I-band light curves of the gravitationally lensed double QSO B1600+434 from observations obtained at the Nordic Optical Telescope (NOT) between April 1998 and November 1999. The photometry has been performed by simultaneous deconvolution of all the data frames, involving a numerical lens galaxy model. Four methods have been applied to determine the time delay between the two QSO components, giving a mean estimate of \Delta_t = 51+/-4 days (95% confidence level). This is the fourth optical time delay ever measured. Adopting a Omega=0.3, Lambda=0 Universe and using the mass model of Maller et al. (2000), this time-delay estimate yields a Hubble parameter of H_0=52 (+14, -8) km s^-1 Mpc^-1 (95% confidence level) where the errors include time-delay as well as model uncertainties. There are time-dependent offsets between the two (appropriately shifted) light curves that indicate the presence of external variations due to microlensing.

An unbiased method for improving the resolution of astronomical images is presented. The strategy at the core of this method is to establish a linear transformation between the recorded image and an improved image at some desirable resolution. In order to establish this transformation only the actual point spread function and a desired point spread function need be known. Any image actually recorded is not used in establishing the linear transformation between the recorded and improved image. This method has a number of advantages over other methods currently in use. It is not iterative which means it is not necessary to impose any criteria, objective or otherwise, to stop the iterations. The method does not require an artificial separation of the image into ``smooth'' and ``point-like'' components, and thus is unbiased with respect to the character of structures present in the image. The method produces a linear transformation between the recorded image and the deconvolved image and therefore the propagation of pixel-by-pixel flux error estimates into the deconvolved image is trivial. It is explicitly constrained to preserve photometry.

One of the most well-known tests of General Relativity (GR) results from combining measurements of the anomalous precession of the orbit of Mercury with a determination of the gravitational quadrupole moment of the Sun J_2. The latter can be done by inference from an integral relation between J_2 and the solar internal rotation. New observational data of high quality obtained from the Solar Heliospheric Satellite (SoHO) and from the Global Oscillations Network Group (GONG), allow the determination of the internal rotation velocity of the Sun as a function of radius and latitude with unprecedented spatial resolution and accuracy. As a consequence, a number of global properties of the Sun can also be determined with much higher accuracy, notably the gravitational quadrupole moment of the Sun. The anomalous precession of the orbit of Mercury is primarily due to GR effects but there are classical corrections the largest of which is that due to J_2. It is shown here that the data are currently consistent with the predictions of GR.

From the observing campaigns of a number of helioseismic telescope networks such as the Global Oscillation Network Group (GONG) and also from the Solar Heliospheric Observatory satellite (SoHO), helioseismologists now have data on in excess of 20000 oscillation frequencies of the modes known collectively as the 5 minute oscillations. Depending on the order and degree of the modes these data are available either as individual frequencies of modes with known radial order n, spherical harmonic degree l, and order m, or as so-called `a-coefficients' which are obtained by fitting polynomials to ridges of power in the (l,\nu) diagram. Although there is a one-to-one relationship between these two types of data, the conversion between them is not straightforward which complicates comparison of helioseismic inference between the two types of data. Their relationship with each other and with the solar rotation rate is clarified in this paper.

The residence time of dust grains inside turbulent molecular clouds is calculated. Earlier estimates of this time by Boland & de Jong (1982) have been criticized by Prasad et al. (1987) for not taking into account properly the turbulent character of the velocity field in molecular clouds which implies that the trajectory of a dust grain should be treated as a random walk. Taking some minimal assumptions regarding the turbulent velocity field a time scale is derived that depends on the Reynolds number of the flow. The near proportionality with Reynolds number of this time scale results in a much longer time scale over which dust grains will remain inside a molecular cloud.

A common problem in astronomy is the determination of the time shift between two otherwise identical time series of measured flux from a variable source, in short the determination of a time lag. Two examples of where this problem occurs are in the determination of the Hubble constant from multiple images of gravitationally lensed variable quasars and also in the determination of distances to OH/IR stars. It is shown here that this problem can be seen as a restricted inversion problem. It therefore can be solved using the subtractive optimally localized averages (SOLA) method for inversion which has been described elsewhere (Pijpers & Thompson 1992, 1994 ; Pijpers & Wanders 1994).

We present an efficient method for two dimensional inversions for the solar rotation rate using the Subtractive Optimally Localized Averages (SOLA) method and a modification of the R1 X R1 technique proposed by Sekii (1993). The SOLA method is based on explicit construction of averaging kernels similar to the Backus-Gilbert method. The versatility and reliability of the SOLA method in reproducing a target form for the averaging kernel, in combination with the idea of the R1 X R1 decomposition, results in a computationally very efficient inversion algorithm. This is particularly important for full 2-D inversions of helioseismic data in which the number of modes runs into at least tens of thousands.

This paper was presented at the Institute for Mathematics and its Applications workshop "Inverse problems in wave propagation" and will appear in the series IMA volumes (Springer). A brief overview of applications of inversions within astronomy is presented and also an inventory of techniques commonly in use. Most of this paper is focussed on the method of Subtractive Optimally Localized Averages (SOLA) which is an adaptation of the Backus and Gilbert method. This method was originally developed for use in helioseismology where the Backus and Gilbert method is computationally too slow. Since then it has also been applied to the problem of reverberation mapping of active galactic nuclei and the differences between this inverse problem and the ones of helioseismology are also discussed.

It has been claimed that in stellar winds traversed by strong shocks the mechanism for driving the wind by sound wave pressure cannot operate because sound waves cannot propagate past the shocks. It is shown here that sound waves can propagate through shocks in one direction and that this is a sufficient condition for the sound wave pressure mechanism to work. A strong shock amplifies a sound wave passing through it and can drag the sound wave away from the star. It is immaterial for the sound wave pressure gradient that the sound wave vector points towards the star. Since the strong shocks drag the sound waves away, the star itself is the source for the sound waves propagating towards it.

In this paper a new method is presented to find the transfer function of the broad-line region in active galactic nuclei. The subtractive optimally localized averages (SOLA) method is a modified version of the Backus-Gilbert method and is presented as an alternative to the more often used maximum-entropy method. The SOLA method has been developed for use in helioseismology. It has been applied to the solar oscillation frequency splitting data currently available to deduce the internal rotation rate of the sun. The original SOLA method is reformulated in the present paper to cope with the slightly different problem of inverting time series. We use simulations to test the viability of the method and apply the SOLA method to the real data of the Seyfert-1 galaxy NGC 5548. We investigate the effects of measurement errors and how the resolution of the TF critically depends upon both the sampling rate and the photometric accuracy of the data. A uuencoded compressed postscript file of the paper which includes the figures is available by anonymous ftp at ftp://solaris.astro.uu.se/pub/articles/atmos/frank/PijWan.uue

We present the results of a short time scale monitoring of SiO maser emission (v=1 J=1-0 transition) in four known strong sources. These sources were monitored nightly for a period of about a month. The aim of these observations is to investigate the possible presence of variations in the maser lines on time scales of a few days to weeks, due to sound waves propagating out from the central star. If sound waves are responsible for the mass loss of certain cool giants, as suggested by Pijpers and Hearn (1989) and Pijpers and Habing (1989), local variations in density and relative velocity are expected just above the stellar photosphere. These could give rise to variations in any narrow spectral line formed in this region, and therefore in particular in the SiO maser lines. Our observations indicate that variations in the line shape (leading to relative changes in the intensity of about 20\%) occur in the SiO emission of Mira type stars, within short time scales of 10-20 days. The main component of the profile variability is consistent with a displacement of the velocity centroid of the dominant maser peaks, by about 1 km s$^{-1}$ in the average. Apparent variations in the total line flux were also found, but could be partially due to calibration uncertainties.