基于小波方法的平滑算法在引力透镜中的应用

基于离散小波变换(DWT)方法,提出了一种可用于计算三维数值模拟样本面密度的平滑算法.为检验方法的有效性,利用该算法研究了两组不同质量解析度的引力透镜数值模拟样本,样本采用了暗物质晕的等温椭球模型,使用蒙特卡罗方法生成.计算结果表明此算法能够在很高的精度上构建引力透镜模拟样本的面密度分布轮廓,由面密度计算出来的透镜的临界曲线和焦散曲线也能较好地和理论曲线吻合,结果是令人满意的.同时比较了三组不同的小波基的计算结果,包括Daub4,Daub6和B-spline 3th,给出了最优的选择.在不损失平滑效果的同时,此算法具有非常高的速度,非常适合于处理以后更高精度的N体数值模拟.     

Spectral analyses of solar photospheric fluctuations

The application of the fast-Fourier-transform (FFT) algorithm to calculating one-dimensional and bi-dimensional (temporal and spatial), power and cross-power (coherence and phase) spectra is examined for solar photospheric fluctuations. Alternative methods for smoothing raw spectra, direct averaging (employing various weights) and indirect truncation of the correlation function, are compared, and indirect smoothing is compared with spectra calculated by mean-lagged-product (MLP) methods. Besides providing the raw spectrum, FFT techniques easily allow computing a series of spectra with varying amounts of smoothing. From these spectra a range of satisfactory compromise between resolution and stability can be determined which helps in the interpretation of spectral trends, and in identifying more clearly the existence and significance of spectral features. For bi-dimensional spectra presented as contour plots, this range of satisfactory smoothing can be restricted, particularly when spectral trends must be represented by small-scale contours. Equivalent spectra (i.e. comparable equivalent degrees of freedom) computed or smoothed by different methods have minor, but not negligible, differences. Examination of these differences favors computing of FFT spectra smoothed by averaging for photospheric fluctuations.     

An optimal method for the power spectrum measurement

An aliasing effect brought up by mass assignment onto Fast Fourier Transformation (FFT) grids may bias measurement of the power spectrum of large scale structures. In this paper, based on the Beylkin's unequally spaced FFT technique, we propose a new precise method to extract the true power spectrum of a large discrete data set. We compare the traditional mass assignment schemes with the new method using the Daub6 and the 3rd-order B-spline scaling functions. Our measurement of Poisson samples and samples of N-body simulations shows that the B-spline scaling function is an optimal choice for mass assignment in the sense that (1) it has a compact support in real space and thus yields an efficient algorithm (2) without any extra corrections. The Fourier space behavior of the 3rd-order B-spline scaling function enables it to be able to accurately recover the true power spectrum with errors less than 5% up to k ＜ kN. It is expected that such a method can be applied to higher order statistics in Fourier space and will enable us to have a precision capture of the non-Gaussian features in the large scale structure of the universe.     

Retrieval of atmospheric mass densities in lower thermosphere below 200 km from precise orbit of re-entry object CZ-3B R/B by analytical and numerical methods

Taking the re-entry object CZ-3B R/B (COSPAR identifier 2012-018D, NORAD catalog number 38253) as an example, retrieval of atmospheric mass densities in lower thermosphere below 200 km from its rebuilt precise orbit is studied in this paper. Two methodologies, i.e. analytical and numerical methods, are adopted in the retrieval. Basic principles of these two methodologies are briefly introduced. Based on the short-arc sparse observational data accumulated in the high accuracy re-entry prediction, orbit determinations of re-entry object CZ-3B R/B are performed sectionally, and then its precise orbit is rebuilt. According to the orbit theory, the variation of orbital semi-major axis of re-entry object CZ-3B R/B induced by atmospheric drag perturbation only is derived from the rebuilt precise orbit. In the derivation of secular change of the orbital semi-major axis of re-entry object CZ-3B R/B induced by atmospheric drag perturbation only, the time-span is set as one minute tentatively. And then retrieval results of atmospheric mass densities in lower thermosphere below 200 km by analytical and numerical methods are presented, as well as their bias deviations from the calculated results of the NRLMSISE-00 empirical model of the atmosphere. Setting bias deviation bands, the corresponding ‘confidence coefficients’ of the retrieved atmospheric mass densities with respect to the model values are given. Average bias deviations of the retrieved atmospheric mass densities by analytical and numerical methods from the model values are also calculated respectively. On the whole, the retrieved atmospheric mass densities by numerical method approach to the model values more closely; the differences between the retrieved results and the model values are relatively smaller at the peaks of atmospheric mass densities than the other places.     

Radiation from a Superbolide

Numerical simulation of the destruction, evaporation, deceleration, and emission of the Chelyabinsk superbolide has been carried out. The model assumes that the main energy is radiated in the stage when the asteroid is already completely destroyed and does not have solidity (quasi-liquid approximation). The radiation transfer during the motion is taken into account in the approximation of radiative heat conductivity and volumetric emission. The distributions of temperatures and densities are obtained at the moments when the bolide is at different altitudes. The intensity of radiation at the Earth’s surface is calculated at certain times by solving the radiative transfer equation along the rays passing through the luminous region using the air and LL-chondrite vapor absorption coefficients. The features of superbolide radiation, the contribution of air and vapor to radiation, the size of the luminous region, and the radiation spectrum have been considered. The calculated efficiency of radiation—17% of the kinetic energy of a cosmic body—agrees with the results of observations. It is shown that due to anisotropy of the superbolide radiation, the determination of luminous efficiency from measurements can depend on the observation point. For estimations, the pointsource approximation can be used, but in general, the source luminous efficiency is unknown, and its location is determined with some error; therefore, numerical simulation is required to reliably estimate the consequences of space body falls.     

'Stokes imaging' of the accretion region in magnetic cataclysmic variables — I. Conception and realization

Modelling the polarized cyclotron emission from magnetic cataclysmic variables has been a pivotal technique for determining the structure of the accretion zones on the white dwarf. To date, model solutions have been obtained from trial fits to the intensity and polarization data, which have been constructed from emission regions (for example arcs and spots) put in by hand. These models were all inferred indirectly from arguments based on the polarization and X-ray light curves.   We present a more analytical and objective technique using optimization by a genetic algorithm, Tikhonov regularization and Powell's method that robustly models the details of polarized emission.   To demonstrate the success of this technique, we show the results of several simulations in which we calculated the intensity and polarization curves from arbitrarily shaped emission regions on the surface of a sphere and then applied our code to these curves to recover the original test data. We also show how adding artificial noise affects the outcome of the optimization technique.     

Chromaticity of gravitational microlensing events

In this paper, we investigate the colour changes of gravitational microlensing events caused by the two different mechanisms of differential amplification for a limb-darkened extended source and blending. From this investigation, we find that the colour changes of limb-darkened extended source events (colour curves) have dramatically different characteristics depending on whether the lens transits the source star or not. We show that for a source transit event, the lens proper motion can be determined by simply measuring the turning time of the colour curve instead of fitting the overall colour or light curves. We also find that even for a very small fraction of blended light, the colour changes induced by blending are equivalent to those induced by limb darkening, causing serious distortion in the observed colour curve. Therefore, to obtain useful information about the lens and source star from the colour curve of an event, it will be essential to correct for blending. We discuss various methods of blending correction .     

433 Eros Global Basemap from NEAR Shoemaker MSI Images

During its year-long orbital survey of Asteroid 433 Eros (February 14, 2000 to February 12, 2001) the NEAR Shoemaker spacecraft returned over 160,000 images of the surface, obtained under a wide variety of viewing conditions and resolutions. To handle this large volume of images of Eros, we have utilized specialized techniques for projecting, mosaicking, and photometrically analyzing images of irregular objects. Local incidence and emission angles are calculated on a pixel-by-pixel basis from the global shape model to produce a normalized reflectance map of the entire asteroid. This map provides not only a medium resolution (∼10 m/pixel) scientific product in its own right but also a regional context for high-resolution (<1 m/pixel) image sequences.     

Planetary microlensing signals from the orbital motion of the source star around the common barycentre

With several detections, the technique of gravitational microlensing has proven useful for studying planets that orbit stars at Galactic distances, and it can even be applied to detect planets in neighbouring galaxies. So far, planet detections by microlensing have been considered to result from a change in the bending of light and the resulting magnification caused by a planet around the foreground lens star. However, in complete analogy to the annual parallax effect caused by the revolution of the Earth around the Sun, the motion of the source star around the common barycentre with an orbiting planet can also lead to observable deviations in microlensing light curves that can provide evidence for the unseen companion. We discuss this effect in some detail and study the prospects of microlensing observations for revealing planets through this alternative detection channel. Given that small distances between lens and source star are favoured, and that the effect becomes nearly independent of the source distance, planets would remain detectable even if their host star is located outside the Milky Way with a sufficiently good photometry (exceeding present-day technology) being possible. From synthetic light curves arising from a Monte Carlo simulation, we find that the chances for such detections are not overwhelming and appear practically limited to the most massive planets (at least with current observational set-ups), but they are large enough for leaving the possibility that one or the other signal has already been observed. However, it may remain undetermined whether the planet actually orbits the source star or rather the lens star, which leaves us with an ambiguity not only with respect to its location, but also to its properties.     

亚毫米波射电天文用超导SIS接收机的研制(英文)

为了研制亚毫米波射电天文用超导SIS(超导 -绝缘体 -超导 )接收机 ,我们重点开展了如下研究 ,1 )Nb超导SIS结在其能隙频率附近的量子混频特性 ,及其结合高能隙超导薄膜 (NbTiN)和高电导率金属薄膜 (Al)分布结阵在 780 - 950GHz频率区间的量子混频特性 ;2 )亚毫米波超导混频器嵌入阻抗的数值和实验表征 ;3)高电流密度小面积Nb超导SIS结的制备和特性表征 ;4)一个 60 0 - 72 0GHz超导SIS混频器的研制和特性表征。本文详细介绍了相关的数值分析和实验测量结果。     

四种综合脉冲星时算法比较

由单颗脉冲星定义的脉冲星时受多种噪声源的影响,其短期和长期稳定度都不够好.为了削弱这些噪声源对单脉冲星时的影响,可以采取合适的算法对多个单脉冲星时进行综合得到综合脉冲星时,从而提高综合脉冲星时的长期稳定度.文中介绍4种综合脉冲星时算法:经典加权算法、小波分析算法、维纳滤波算法和小波域中的维纳滤波算法,将这4种算法分别应用于Arecibo天文台对两颗毫秒脉冲星PSR B1855+09和PSRB1937+21观测得到的计时残差并作出比较.     

Simulations of multigroup relativistic radiative transfer for shock waves in supernovae

The results of numerical simulations of light curves and spectra and the influence of special relativity effects at the epoch of supernova shock breakout are considered. The algorithm of the RADA code used in the simulations is described. Prospects for using the numerical simulation results to analyze and interpret available and future data from the SWIFT and LOBSTER space observatories are discussed.     

Assessment of Kinetic Energy of Meteoroids Detected by Satellite-Based Light Sensors

Radiation energies of bright flashes caused by disintegration of large meteoroids in the atmosphere have been measured using optical sensors on board geostationary satellites. Light curves versus time are available for some of the events. We have worked out several numerical techniques to derive the kinetic energy of the meteoroids that produced the flashes. Spectral opacities of vapor of various types of meteoroids were calculated for a wide range of possible temperatures and densities. Coefficients of conversion of kinetic energy to radiation energy were computed for chondritic and iron meteoroids 10 cm to 10 m in size using radiation–hydrodynamics numerical simulations. Luminous efficiency increases with body size and initial velocity. Some analytical approximations are presented for average conversion coefficients for irons and H-chondrites. A mean value of this coefficient for large meteoroids (1–10 m in size) is about 5–10%. The theory was tested by analyzing the light curves of several events in detail.Kinetic energies of impactors and energy–frequency distribution of 51 bolides, detected during 22 months of systematic observations in 1994–1996, are determined using theoretical values of luminous efficiencies and heat-transfer coefficients. The number of impacts in the energy range from 0.25 to 4 kt TNT is 25 per year and per total surface of the Earth.The energy–frequency distribution is in a rather good agreement with that derived from acoustic observations and the lunar crater record. Acoustic systems have registered one 1 Mt event in 12 years of observation. Optical systems have not detected such an event as yet due to a shorter time of observation. The probability of a 1 Mt impact was estimated by extrapolation of the observational data.     

Atlas of lunar radar maps at 70-cm wavelength

The intensity distribution of lunar radar echoes has been mapped for two-thirds of the earth-visible lunar surface at a wavelength of 70 cm. The depolarizing effects of the lunar surface were observed by simultaneously receiving the radar echoes in opposite polarizations. These echoes were mapped with areal resolutions of 25–100 km². Mappings with this resolution confirmed that the young craters have enhanced returns. A few craters were found to have enhanced echoes only from their rims. Backscattering differences were also observed between various areas within a mare, between different highland areas, and between maria and adjacent highlands. These scattering differences were interpreted with a simple model, which assumed that the surface backscattered with varying amounts of quasi-specular and diffuse power. Only an increase in the diffuse power was needed to give the numerical values of the enhancements.     

Mineralogical and spectroscopic investigation of the Tagish Lake carbonaceous chondrite by X‐ray diffraction and infrared reflectance spectroscopy

Abstract– We have carried out a sample‐correlated spectroscopic and mineralogical investigation of samples from seven different collection sites of the Tagish Lake C2 chondrite. Rietveld refinement of high‐resolution powder X‐ray diffraction (XRD) data was used to determine quantitative major mineral abundances. Thermal infrared (400–4500 cm⁻¹, 2.2–25.0 μm) spectra of the same samples were obtained using diffuse (biconical) reflectance infrared Fourier transform spectroscopy (DRIFTS). Our results are in good agreement with previous studies of the mineralogy of the Tagish Lake meteorite; we find however that Tagish Lake is more varied in major mineralogy than has previously been reported. In particular, we observed two new distinct lithologies, an inclusion‐poor magnetite‐ and sulfide‐rich lithology, and a carbonate‐rich, siderite‐dominated lithology in addition to the previously documented carbonate‐rich and carbonate‐poor lithologies. Grain density for each Tagish Lake sample was calculated from the measured mineral modal abundances and known mineral densities. For powders from three originally intact inclusion‐rich samples, the calculated grain density is 2.77 ± 0.05 g cm⁻³, in excellent agreement with those reported in the literature for other intact inclusion‐rich Tagish Lake samples. Tagish Lake disaggregated samples have a significantly higher calculated grain density due to their lower saponite‐serpentine content, likely a result of mineral separation in the meltwater holes from which they were collected; the disaggregated samples may not therefore adequately represent bulk samples of the Tagish Lake meteorite. The predominance of very fine‐grained material in the Tagish Lake samples investigated in this study is expected to produce infrared spectra representative of asteroidal regolith. Gypsum and talc have been found by XRD in powders from the inclusion‐rich, intact Tagish Lake samples in this study, and may have been present in the parent body; if present, these hydrous sulfates would complicate the interpretation of possible hydrated mineral features in asteroid infrared spectra.     

Bounded relative motion under zonal harmonics perturbations

The problem of finding natural bounded relative trajectories between the different units of a distributed space system is of great interest to the astrodynamics community. This is because most popular initialization methods still fail to establish long-term bounded relative motion when gravitational perturbations are involved. Recent numerical searches based on dynamical systems theory and ergodic maps have demonstrated that bounded relative trajectories not only exist but may extend up to hundreds of kilometers, i.e., well beyond the reach of currently available techniques. To remedy this, we introduce a novel approach that relies on neither linearized equations nor mean-to-osculating orbit element mappings. The proposed algorithm applies to rotationally symmetric bodies and is based on a numerical method for computing quasi-periodic invariant tori via stroboscopic maps, including extra constraints to fix the average of the nodal period and RAAN drift between two consecutive equatorial plane crossings of the quasi-periodic solutions. In this way, bounded relative trajectories of arbitrary size can be found with great accuracy as long as these are allowed by the natural dynamics and the physical constraints of the system (e.g., the surface of the gravitational attractor). This holds under any number of zonal harmonics perturbations and for arbitrary time intervals as demonstrated by numerical simulations about an Earth-like planet and the highly oblate primary of the binary asteroid (66391) 1999 KW4.     

Strong lensing optical depths in a ΛCDM universe

We investigate strong gravitational lensing in the concordance ΛCDM cosmology by carrying out ray tracing along past light cones through the Millennium Simulation, the largest simulation of cosmic structure formation ever carried out. We extend previous ray-tracing methods in order to take full advantage of the large volume and the excellent spatial and mass resolution of the simulation. As a function of source redshift we evaluate the probability that an image will be highly magnified, will be highly elongated or will be one of a set of multiple images. We show that such strong lensing events can almost always be traced to a single dominant lensing object and we study the mass and redshift distribution of these primary lenses. We fit analytic models to the simulated dark haloes in order to study how our optical depth measurements are affected by the limited resolution of the simulation and of the lensing planes that we construct from it. We conclude that such effects lead us to underestimate total strong lensing cross-sections by about 15 per cent. This is smaller than the effects expected from our neglect of the baryonic components of galaxies. Finally we investigate whether strong lensing is enhanced by material in front of or behind the primary lens. Although strong lensing lines of sight are indeed biased towards higher than average mean densities, this additional matter typically contributes only a few per cent of the total surface density.     

A genetic algorithm for the non-parametric inversion of strong lensing systems

We present a non-parametric technique to infer the projected mass distribution of a gravitational lens system with multiple strong-lensed images. The technique involves a dynamic grid in the lens plane on which the mass distribution of the lens is approximated by a sum of basis functions, one per grid cell. We used the projected mass densities of Plummer spheres as basis functions. A genetic algorithm then determines the mass distribution of the lens by forcing images of a single source, projected back on to the source plane, to coincide as well as possible. Averaging several tens of solutions removes the random fluctuations that are introduced by the reproduction process of genomes in the genetic algorithm and highlights those features common to all solutions. Given the positions of the images and the redshifts of the sources and the lens, we show that the mass of a gravitational lens can be retrieved with an accuracy of a few percent and that, if the sources sufficiently cover the caustics, the mass distribution of the gravitational lens can also be reliably retrieved. A major advantage of the algorithm is that it makes full use of the information contained in the radial images, unlike methods that minimize the residuals of the lens equation, and is thus able to accurately reconstruct also the inner parts of the lens.     

A new data structure for accelerating kinetic Monte Carlo method

The kinetic Monte Carlo simulation is a rigorous numerical approach to study the chemistry on dust grains in cold dense interstellar clouds. By tracking every single reaction in chemical networks step by step, this approach produces more precise results than other approaches but takes too much computing time. Here we present a method of a new data structure, which is applicable to any physical conditions and chemical networks, to save computing time for the Monte Carlo algorithm. Using the improved structure,the calculating time is reduced by 80 percent compared with the linear structure when applied to the osu-2008 chemical network at 10K. We investigate the effect of the encounter desorption in cold cores using the kinetic Monte Carlo model with an accelerating data structure. We found that the encounter desorption remarkably decreases the abundance of grain-surface H_2 but slightly influences the abundances of other species on the grain.