基于精密光程差调制的时域干涉信号闭合相位检测方法研究

闭合相位法是实现长基线恒星光干涉高分辨成像的重要技术手段之一,获得精确的闭合相位信息是进行光干涉图像重构的先决条件.提出一种基于精密光程差调制的时域干涉信号闭合相位检测方法,在3路干涉臂上进行非冗余精密光程调制,并通过多次干涉测量结合数据拟合的方法消除光程差调制中存在的正弦误差,使得光程调制的精度达到20 nm以内.引入高速探测器件提升时域干涉信号的采样频率,对探测器上获得的时域干涉信号进行傅立叶变换处理,获得3路干涉臂精确的闭合相位信息.室内实验结果表明,基于精密光程调制的时域信号闭合相位计算精度可以达到1/50波长以内.     

Imaging the heart of astrophysical objects with optical long-baseline interferometry

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: (1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; (2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long-baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image-reconstruction process was ??blind?? in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image-reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with six to ten telescopes, should be well suited to provide images of complex sceneries.     

Estimating Solar Flux Density at Low Radio Frequencies Using a Sky Brightness Model

Sky models have been used in the past to calibrate individual low radio frequency telescopes. In this article we generalize this approach from a single antenna to a two element interferometer, and formulate the problem in a way that allows us to estimate the flux density of the Sun using the normalized cross-correlations (visibilities) measured on a low resolution interferometric baseline. For wide field-of-view instruments, typically the case at low radio frequencies, this approach can provide robust absolute solar flux calibration for well characterized antennas and receiver systems. It can provide a reliable and computationally lean method for extracting parameters of physical interest using a small fraction of the voluminous interferometric data, which can be computationally prohibitively expensive to calibrate and image using conventional approaches. We demonstrate this technique by applying it to data from the Murchison Widefield Array and assess its reliability.     

Solar flares in the EUV observed from OSO-5

Solar flares in three broad EUV spectral bands have been observed from OSO-5 with a grating spectrophotometer. Results are given for three large flares of March 12, March 21 and April 21, 1969. In general the time dependence of flare intensity in each band is characterized by a slowly varying component with impulsive bursts superimposed. Bands 2 (465–630 Å) and 3 (760–1030 Å) are quite similar in their time variations, but band 1 (280–370 Å) shows less impulsive structure, and declines more slowly. Absolute EUV intensities for the flares are estimated, and a comparison made with the 2800 mc s^–1 radio emission. A flare model is proposed to account for the EUV time variations during a large flare.     

Why chromatic imaging matters

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u ? v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction.     

Compressed sensing imaging techniques for radio interferometry

Radio interferometry probes astrophysical signals through incomplete and noisy Fourier measurements. The theory of compressed sensing demonstrates that such measurements may actually suffice for accurate reconstruction of sparse or compressible signals. We propose new generic imaging techniques based on convex optimization for global minimization problems defined in this context. The versatility of the framework notably allows introduction of specific prior information on the signals, which offers the possibility of significant improvements of reconstruction relative to the standard local matching pursuit algorithm CLEAN used in radio astronomy. We illustrate the potential of the approach by studying reconstruction performances on simulations of two different kinds of signals observed with very generic interferometric configurations. The first kind is an intensity field of compact astrophysical objects. The second kind is the imprint of cosmic strings in the temperature field of the cosmic microwave background radiation, of particular interest for cosmology.     

Speckle interferometry of nearby multiple stars. V. 2002–2006 positional measurements

We present the measurements of positional parameters of 194 nearby binary stars performed at the 6-m BTA telescope of the SAO RAS from 2002 through 2006 applying the speckle interferometric technique. The observations were conducted at central filter wavelengths ranging from 545 to 800 nm using a speckle interferometer equipped with a fast CCD and a three-stage image intensifier. A significant part of the observed systems (80 stars) are pairs, the duality of which was discovered or suspected from the Hipparcos satellite observations. The remaining stars are visual binaries and interferometric binary systems with orbital motion period estimates from several to tens of years, as well as the pairs with slow relative motion of the components, used for gaging the image scales and positional angle calibrations.     

Point source detection performance of Hard X-ray Modulation Telescope imaging observation

The Hard X-ray Modulation Telescope(HXMT) will perform an all-sky survey in the hard X-ray band as well as deep imaging of a series of small sky regions.We expect various compact objects to be detected in these imaging observations. Point source detection performance of HXMT imaging observation depends not only on the instrument but also on the data analysis method that is applied since images are reconstructed from HXMT observed data with numerical methods. The denoising technique used plays an important part in the HXMT imaging data analysis pipeline along with demodulation and source detection. In this paper we have implemented several methods for denoising HXMT data and evaluated the point source detection performances in terms of sensitivities and location accuracies. The results show that direct demodulation with 1-fold cross-correlation should be the default reconstruction and regularization method, although both sensitivity and location accuracy could be further improved by selecting and tuning numerical methods in data analysis used for HXMT imaging observations.     

A laboratory discharge lamp source for the O2 atmospheric bands and OI 5577 Å and 6300 Å lines of a known line width

A laboratory discharge lamp is described which strongly emits the forbidden OI 5577 and 6300 Å lines and the O₂ (0?0) atmospheric band. Experimental measurements confirm that these atomic and molecular species are in thermal equilibrium with one another, so that a rotational temperature measurement of the O₂ atmospheric band allows one to deduce the line widths of the 5577 and 6300 Å emissions. This thus provides a useful calibration source for interferometric measurements of these emissions.     

Exploring three faint source detections methods for aperture synthesis radio images

Wide-field radio interferometric images often contain a large population of faint compact sources. Due to their low intensity/noise ratio, these objects can be easily missed by automated detection methods, which have been classically based on thresholding techniques after local noise estimation. The aim of this paper is to present and analyse the performance of several alternative or complementary techniques to thresholding. We compare three different algorithms to increase the detection rate of faint objects. The first technique consists of combining wavelet decomposition with local thresholding. The second technique is based on the structural behaviour of the neighbourhood of each pixel. Finally, the third algorithm uses local features extracted from a bank of filters and a boosting classifier to perform the detections. The methods’ performances are evaluated using simulations and radio mosaics from the Giant Metrewave Radio Telescope and the Australia Telescope Compact Array. We show that the new methods perform better than well-known state of the art methods such as SExtractor, SAD and DUCHAMP at detecting faint sources of radio interferometric images.     

The accurate determination of stellar angular diameters using broad-band stellar interferometry

The angular diameter of a star can be estimated from interferometric observations by fitting the data with the visibility function for a uniformly illuminated disc and then using published correction factors to convert the uniform-disc angular diameter to the limb-darkened angular diameter. The correction factors are strictly valid only for monochromatic light. We investigate the effect of using a broad bandwidth, and present a simple method for calculating broad-band correction factors from the monochromatic factors.
The technique of fitting the data with a uniform-disc visibility function is only useful for stars with compact atmospheres and 'typical' limb-darkening profiles. It should not be applied to stars with extended atmospheres or that show extreme limb darkening. These stars have visibility functions that are qualitatively different from a uniform-disc visibility function, so they can be distinguished observationally from compact-atmosphere stars.     

Two‐dimensional speckle spectroscopy of Hα features

In May 2002, the solar chromosphere was observed with a two‐dimensional spectrometer which is mounted in the German Vacuum Tower Telescope (VTT) at the Observatorio del Teide/Tenerife. The aim of this observation was to investigate the fine structure of the solar chromosphere seen in Hα. We took narrow‐band filtergrams (Δλ ≈ 72 mÅ) by scanning through this line. Broad‐band images taken strictly simultaneously with the narrow‐band filtergrams were restored by speckle methods. The instantaneous optical transfer function from this restoration procedure was used for the reconstruction of the narrow‐band images. Some results of this high spatial resolution observation are presented below.     

Few projections astrotomography: radio astronomical approach to 3D reconstruction

We develop a radio astronomical approach to 3D‐reconstruction in few projections tomography. It is based on the 2‐CLEAN DSA method which consists of two clean algorithms by using a synthesized beam. In complex cases two extreme solutions are used for the analysis of the image structure. These solutions determine the limits of permissible energy redistribution on the image among the components of small and large scales. Two variants of 3D‐reconstruction proceeding from a set of two‐dimensional projections (3D_2D) and from a set of one‐dimensional ones (3D_1D) are considered. It is shown that the quality of 3D_2D‐reconstruction should be similar to the quality of 2D_1D‐reconstruction if the same number of equally spaced scans is used. But a doubled number of projections is required for 3D_1D‐reconstruction. We have simulated 3D‐reconstruction of an optically thin emitting object. The present technique is a component of astrotomography and it has good prospects for a wide range of remote sensing. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Gaussianity of the cosmic microwave background: smooth goodness-of-fit tests applied to interferometric data

We adapt the smooth tests of goodness-of-fit developed by Rayner and Best to the study of the non-Gaussianity of interferometric observations of the cosmic microwave background (CMB). The interferometric measurements (visibilities) are transformed into signal-to-noise ratio eigenmodes, and then the method is applied directly in Fourier space. This transformation allows us to perform the analysis in different subsets of eigenmodes according to their signal-to-noise ratio level. The method can also deal with non-uniform or incomplete coverage of the UV plane. We explore here two possibilities: we analyse either the real and imaginary parts of the complex visibilities (Gaussianly distributed under the Gaussianity hypothesis) or their phases (uniformly distributed under the Gaussianity hypothesis). The power of the method in discriminating between Gaussian and non-Gaussian distributions is studied by using several kinds of non-Gaussian simulations. On the one hand, we introduce a certain degree of non-Gaussianity directly into the Fourier space using the Edgeworth expansion, and afterwards the desired correlation is introduced. On the other hand, we consider interferometric observations of a map with topological defects (cosmic strings). To these previous non-Gaussian simulations we add different noise levels and quantify the required signal-to-noise ratio necessary to achieve a detection of these non-Gaussian features. Finally, we have also studied the ability of the method to constrain the so-called non-linear coupling constant f _NL using χ² simulations. The whole method is illustrated here by application to simulated data from the Very Small Array interferometer.     

Grad–Shafranov Reconstruction of Magnetic Clouds: Overview and Improvements

The Grad–Shafranov reconstruction is a method of estimating the orientation (invariant axis) and cross section of magnetic flux ropes using the data from a single spacecraft. It can be applied to various magnetic structures such as magnetic clouds (MCs) and flux ropes embedded in the magnetopause and in the solar wind. We develop a number of improvements of this technique and show some examples of the reconstruction procedure of interplanetary coronal mass ejections (ICMEs) observed at 1 AU by the STEREO, Wind, and ACE spacecraft during the minimum following Solar Cycle 23. The analysis is conducted not only for ideal localized ICME events but also for non-trivial cases of magnetic clouds in fast solar wind. The Grad–Shafranov reconstruction gives reasonable results for the sample events, although it possesses certain limitations, which need to be taken into account during the interpretation of the model results.     

Ultraviolet, hard X-ray, and gamma-ray emission of solar flares recorded by VUSS-L and SONG instruments in 2001–2003

The results of simultaneous measurements of variations of UV radiation (in a band near the hydrogen L_α line, 121.6 nm) and hard X-ray and gamma-ray radiation (50 keV-200 MeV) performed by the VUSS-L and SONG instruments, respectively, onboard the CORONAS-F spacecraft are presented for periods of solar flares. Variations in the L_α ultraviolet radiation during the impulsive phase of a flare are shown to be synchronous with those of hard X-ray radiation. Temporal variations of UV and X-ray fluxes correspond to the progressive heating of higher and higher regions of the solar atmosphere and the energy transfer from the lower layers of the solar atmosphere to the coronal areas of flare regions. The energy of electrons in beams arising during the impulsive phase of flares can be as high as 500 keV. The velocity of the energy propagation from the regions of its release to the upper layers of the solar atmosphere can reach several tens of kilometers per second.     

Smear fitting: a new image-deconvolution method for interferometric data

A new technique is presented for producing images from interferometric data. The method, 'smear fitting', makes the constraints necessary for interferometric imaging double as a model, with uncertainties, of the sky brightness distribution. It does this by modelling the sky with a set of functions and then convolving each component with its own elliptical Gaussian to account for the uncertainty in its shape and location that arises from noise. This yields much sharper resolution than clean for significantly detected features, without sacrificing any sensitivity. Using appropriate functional forms for the components provides both a scientifically interesting model and imaging constraints that tend to be better than those used by traditional deconvolution methods. This allows it to avoid the most serious problems that limit the imaging quality of those methods. Comparisons of smear fitting to clean and maximum entropy are given, using both real and simulated observations. It is also shown that the famous Rayleigh criterion (resolution = wavelength/baseline) is inappropriate for interferometers as it does not consider the reliability of the measurements.     

Radio and X-ray observations of a multiple impulsive solar burst with high time resolution

A well-developed multiple impulsive microwave burst occurred on February 17, 1979 simultaneously with a hard X-ray burst and a large group of type III bursts at metric wavelengths. The whole event is composed of several subgroups of elementary spike bursts. Detailed comparisons between these three classes of emissions with high time resolution of 0.5 s reveal that individual type III bursts coincide in time with corresponding elementary X-ray and microwave spike bursts. It suggests that a non-thermal electron pulse generating a type III spike burst is produced simultaneously with those responsible for the corresponding hard X-ray and microwave spike bursts. The rise and decay characteristic time scales of the elementary spike burst are 1 s, 1 s and 3 s for type III, hard X-ray and microwave emissions respectively. Radio interferometric observations made at 17 GHz reveal that the spatial structure varies from one subgroup to others while it remains unchanged in a subgroup. Spectral evolution of the microwave burst seems to be closely related to the spatial evolution. The spatial evolution together with the spectral evolution suggests that the electron-accelerating region shifts to a different location after it stays at one location for several tens of seconds, duration of a subgroup of elementary spike bursts. We discuss several requirements for a model of the impulsive burst which come out from these observational results, and propose a migrating double-source model.     

On 3D Reconstruction of Coronal Mass Ejections: I. Method Description and Application to SECCHI-COR Data

The data from SECCHI-COR1 and SECCHI-COR2 coronagraphs onboard the STEREO mission, which was launched in October 2006, provide us with the first-ever stereoscopic images of the Sun’s corona. These observations were found to be useful in inferring the three-dimensional structure of coronal mass ejections (CMEs) and their propagation direction in space. We apply four methods for reconstructing CMEs: i) Forward modeling technique; ii) Local correlation tracking (to identify the same feature in COR Ahead and COR Behind images) plus tie-point reconstruction technique; iii) Center of mass of the structures in a given epipolar plane plus tie-point reconstruction technique; iv) Polarization ratio technique. The four techniques are applied to three structured CMEs observed by COR1 and COR2 instruments, respectively, on 15 May 2007, 31 August 2007, and 25 March 2008. A comparison of the results obtained from the application of the four reconstruction algorithms is presented and discussed.     

First definitive observations of meteor shower particles using a high-power large-aperture radar

We present the first clear observations of meteor shower activity from meteor-head echoes detected by a high-power large-aperture radar (HPLAR). Such observations have been performed at the Jicamarca VHF radar using its interferometric capabilities allowing the discrimination of meteor shower echoes from the much more frequent sporadic meteors. Until now, HPLARs were unable to distinguish meteor shower from the much more common sporadic meteor ones. In this work we have been able to detect and characterize the η-Aquariids (ETA) as well as the Perseids (PER) showers. The shower activity is more conspicuous for the ETA than for the PER shower due to the more favorable geometry. Namely, PER meteors come from low elevation angles, experiencing more filtering due to the combined Earth-atmosphere-radar instrument. In both cases, there is an excellent agreement between the measured mean velocity of the shower echoes and their expected velocity, within a fraction of 1 km s⁻¹. Besides the good agreement with expected visual results, HPLARs observe meteors with a variety of particles sizes and masses, not observed by any other technique. Taking into account the different viewing volumes, compare to optical observations Jicamarca observes more than 1000 times more ETA meteors. Our results indicate that Jicamarca and other HPLARs are able to detect the echoes from meteor showers, but without interferometric capabilities such populations are difficult to identify just from their velocity distributions, particularly if their velocity distributions are expected to be similar to the more dominant distributions of sporadic meteors.