首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
We outline a method for fitting binary-lens caustic-crossing microlensing events based on the alternative model parametrization proposed and detailed by Cassan. As an illustration of our methodology, we present an analysis of OGLE-2007-BLG-472, a double-peaked Galactic microlensing event with a source crossing the whole caustic structure in less than three days. In order to identify all possible models we conduct an extensive search of the parameter space, followed by a refinement of the parameters with a Markov Chain Monte Carlo algorithm. We find a number of low-  χ2  regions in the parameter space, which lead to several distinct competitive best models. We examine the parameters for each of them, and estimate their physical properties. We find that our fitting strategy locates several minima that are difficult to find with other modelling strategies and is therefore a more appropriate method to fit this type of event.  相似文献   

2.
We analyse a series of line profile observations of the He  i 6678 line in ζ Oph. A period analysis on these data using the mode and moments of the line profile confirms the two previously known periods. We describe a new method of mode identification for pulsating stars in which the calculated profiles are directly fitted to observed profiles. The method yields the full set of pulsational parameters including the spherical harmonic degree, ℓ, and azimuthal number, m . Application of the method to these data confirms the mode identifications previously suggested for the two periodicities. We find that the derived pulsational parameters are physically realistic and conclude that non-radial pulsation is the most likely explanation for the travelling subfeatures. However, a unique mode identification is still not possible – several non-sectorial modes fit the data as well as the usually adopted sectorial identifications. The predicted photometric amplitudes are in good accord with upper limits derived from photometric observations. We conclude that ζ Oph is a star in the β Cep instability strip in which two modes of high degree (probably ℓ=4 and ℓ=8) are excited. We present an interpretation of these findings in which the cause of the low-order line profile and light variations in periodic Be stars is corotating photospheric clouds, while the travelling subfeatures are incidental to the Be phenomenon and are a result of non-radial pulsation.  相似文献   

3.
The protoplanetary discs of T Tauri and Herbig Ae/Be stars have previously been studied using geometric disc models to fit their spectral energy distribution (SED). The simulations provide a means to reproduce the signatures of various circumstellar structures, which are related to different levels of infrared excess. With the aim of improving our previous model, which assumed a simple flat-disc configuration, we adopt here a reprocessing flared-disc model that assumes hydrostatic, radiative equilibrium. We have developed a method to optimize the parameter estimation based on genetic algorithms (GAs). This paper describes the implementation of the new code, which has been applied to Herbig stars from the Pico dos Dias Survey catalogue, in order to illustrate the quality of the fitting for a variety of SED shapes. The star AB Aur was used as a test of the GA parameter estimation, and demonstrates that the new code reproduces successfully a canonical example of the flared-disc model. The GA method gives a good quality of fit, but the range of input parameters must be chosen with caution, as unrealistic disc parameters can be derived. It is confirmed that the flared-disc model fits the flattened SEDs typical of Herbig stars; however, embedded objects (increasing SED slope) and debris discs (steeply decreasing SED slope) are not well fitted with this configuration. Even considering the limitation of the derived parameters, the automatic process of SED fitting provides an interesting tool for the statistical analysis of the circumstellar luminosity of large samples of young stars.  相似文献   

4.
The direct and indirect methods for predicting the ionospheric total electron content (TEC) by using the ARIMA (p, d, q) model in the theory of time-series analysis are proposed. The direct method predicts directly the time series composed of the TECs at the every grid point. The indirect method is: at first, to fit the ionospheric TECs with the spherical harmonics and obtain the time series of the fitting coefficients; then, to predict forward the fitting coefficients with the ARIMA (p, d, q) model in the theory of time-series analysis; and finally, to calculate the TEC for the specific time and grid point. Using the ionospheric data of the International GPS Service (IGS) in the period from 1 January 2004 to 31 January 2005, the testing and comparisons on the proposed methods are performed, the results indicate that within 12 days, the results of the two methods are basically consistent, and that for the predictions longer than 12 days, the accuracy of the indirect method is higher than the direct method. In the duration of 20 days, relative to the total ionospheric grid points, the percentage of the grid points with a prediction error less than 3 TECUs is about 80%; as the duration becomes longer, compared with the indirect method, the decline of this percentage is not apparent for the direct method. Obviously, the direct method suits the prediction of regional TECs, and the indirect method suits the prediction of global TECs.  相似文献   

5.
We study the mass distribution in six nearby  ( z < 0.06)  relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters, we approximate their total density distribution by the Navarro, Frenk & White (NFW) formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120–420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from the N -body simulations of dark matter haloes. We find that the best-fitting Galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters, we conclude that the trend of decreasing concentration for higher masses found in the cosmological N -body simulations is consistent with the data. By scaling the individual cluster data by mass, we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range  1.5 < c < 14  and the anisotropy parameter in the range  −1.1 < β < 0.5  at the 95 per cent confidence level.  相似文献   

6.
We present a new method for determining physical parameters of RRab variables exclusively from multicolour light curves. Our method is an inverse photometric Baade–Wesselink analysis which, using a non-linear least-squares algorithm, searches for the effective temperature ( T eff) and pulsational velocity ( V p) curves and other physical parameters that best fit the observed light curves, utilizing synthetic colours and bolometric corrections from static atmosphere models. The T eff and V p curves are initially derived from empirical relations then they are varied by the fitting algorithm. The method yields the variations and the absolute values of the radius, the effective temperature, the visual brightness and the luminosity of individual objects. Distance and mass are also determined. The method is tested on nine RRab stars subjected to Baade–Wesselink analyses earlier by several authors. The physical parameters derived by our method using only the light-curve data of these stars are well within their possible ranges defined by direct Baade–Wesselink and other techniques. A new empirical relation between the I C magnitude and the pulsational velocity is also presented, which allows to construct the V p curve of an RRab star purely from photometric observations to an accuracy of about 3.5 km s−1.  相似文献   

7.
The role of radial velocity (RV) jitter in extrasolar planet search surveys is discussed. Based on the maximum-likelihood principle, improved statistical algorithms for RV fitting and period search are developed. These algorithms incorporate a built-in jitter determination, so that resulting estimations of planetary parameters account for this jitter automatically. This approach is applied to RV data for several extrasolar planetary systems. It is shown that many RV planet search surveys suffer from periodic systematic errors which increase effective RV jitter and can lead to erroneous conclusions. For instance, the planet candidate HD 74156 d may be a false detection made due to annual systematic errors.  相似文献   

8.
赋形旋转抛物面天线母线拟合常用的方法包括多项式整体拟合和等间隔分段拟合,多项式整体拟合结果阶次高、计算量大,且高阶多项式在边缘处拟合的结果容易振荡;而等间隔分段拟合则是对母线数据按等间隔分段,分段方式具有一定的盲目性,容易造成拟合参数多、光滑性差等问题.针对上述问题,提出了一种基于母线拟合残差分布的自适应分段拟合方法,该方法包括初始整体拟合和分段拟合两步,初始整体拟合用以确定各离散点的拟合残差分布,分段拟合先参照残差分布情况对离散点数据分段,然后采用低阶多项式对各段数据进行拟合.经过实例拟合对比,该方法可避免高阶拟合的不稳定性,减少了分段数,更适用于赋形旋转抛物面天线母线的拟合.  相似文献   

9.
Spherical mass density models are used to fit the central surface brightness profiles of early-type galaxies which are generated from Nuker law parameters obtained from the literature. The mass density and the corresponding potential are in an analytical form. It is shown that only a few mass density components are necessary to obtain a good fit and that for all power-law galaxies and for the core galaxies that we consider, most or all of the mass density components must have cusps to provide good fits. The applied quadratic programming fitting allows for a method of deprojection, which is reliable and convenient. The results can be used directly for further dynamical modelling.  相似文献   

10.
The data reduction process for optical emission-line observations of galaxies using the TAURUS-2 Fabry–Perot interferometer mounted on the 3.9-m Anglo-Australian Telescope is described in detail. The initial steps (bias subtraction, flat-fielding, etc.) are the same as for calibration of CCD images, and the wavelength calibration is similar to that in optical spectroscopy. The final steps are specific to Fabry–Perot instruments, and include the fitting of several instrumental parameters and a phase correction to convert the raw ( x ,  y ,  z ) data cube into a useful position–velocity ( α ,  δ ,  v ) cube. Software has been written to assist with the latter steps of the data reduction. H α observations of NGC 1808, NGC 2442 and Circinus are used to demonstrate the reduction process.  相似文献   

11.
The large-scale matter distribution in three different simulations of CDM models is investigated and compared with corresponding results of the Zel'dovich theory of non-linear gravitational instability. We show that the basic characteristics of wall-like structure elements are well described by this theory, and that they can be expressed by the cosmological parameters and a few spectral moments of the perturbation spectrum. Therefore the characteristics of such elements provide reasonable estimates of these parameters. We show that the compressed matter is relaxed and gravitationally confined and manifests itself in the existence of walls as (quasi-)stationary structure elements with a lifetime restricted by their disruption into high-density clouds.
The matter distribution is investigated in both real and redshift spaces. In both cases almost the same particles form the walls, and we estimate differences in corresponding wall characteristics. The same methods are applied to several mock catalogues of 'galaxies', which allows us to characterize a large-scale bias between the spatial distribution of dark matter and of simulated 'galaxies'.  相似文献   

12.
We introduce a modified version of a standard power spectrum ‘peak‐bagging’ technique which is designed to gain some of the advantages that fitting the entire low‐degree p‐mode power spectrum simultaneously would bring, but without the problems involved in fitting a model incorporating many hundreds of parameters. Employing Monte‐Carlo simulations we show that by using this modified fitting code it is possible to determine the true background level in the vicinity of the p‐mode peaks. In addition to this we show how small biases in other mode parameters, which are related to inaccurate estimates of the true background, are also consequently removed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)  相似文献   

13.
During the past decade, several observational and theoretical works have provided evidence of the binary nature of η Carinae. Nevertheless, there is still no direct determination of the orbital parameters, and the different current models give contradictory results. The orbit is, in general, assumed to coincide with the Homunculus equator although the observations are not conclusive. Among all systems, η Car has the advantage that it is possible to observe both the direct emission of line transitions in the central source and its reflection by the Homunculus, which is dependent on the orbital inclination. In this work, we studied the orbital phase-dependent hydrogen Paschen spectra reflected by the south-east lobe of the Homunculus to constrain the orbital parameters of η Car and determine its inclination with respect to the Homunculus axis. Assuming that the emission excess originates in the wind–wind shock region, we were able to model the latitude dependence of the spectral line profiles. For the first time, we were able to estimate the orbital inclination of η Car with respect to the observer and to the Homunculus axis. The best fit occurs for an orbital inclination to the line of sight of   i ∼ 60°± 10°  , and   i *∼ 35°± 10°  with respect to the Homunculus axis, indicating that the angular momenta of the central object and the orbit are not aligned. We were also able to fix the phase angle of conjunction as  ∼−40°  , showing that periastron passage occurs shortly after conjunction.  相似文献   

14.
Using the near-infrared spectral stellar library of Cenarro et al., the behaviour of the Mg  i line at 8807 Å and nearby TiO bands is analyzed in terms of the effective temperature, surface gravity and metallicity of the library stars. New spectroscopic indices for both spectral features – namely MgI and sTiO – are defined, and their sensitivities to different signal-to-noise ratios, spectral resolutions, flux calibrations and sky emission-line residuals are characterized. The two new indices exhibit interesting properties. In particular, MgI is a good indicator of the Mg abundance, whereas sTiO is a powerful dwarf-to-giant discriminator for cold spectral types. Empirical fitting polynomials that reproduce the strength of the new indices as a function of the stellar atmospheric parameters are computed, and a fortran routine with the fitting function predictions is made available. A thorough study of several error sources, non-solar [Mg/Fe] ratios and their influence on the fitting function residuals is also presented. From this analysis, an [Mg/Fe] underabundance of  ∼−0.04  is derived for the Galactic open cluster M67.  相似文献   

15.
We study the distribution of dark matter in dwarf spheroidal galaxies by modelling the moments of their line-of-sight velocity distributions. We discuss different dark matter density profiles, both cuspy and possessing flat density cores. The predictions are made in the framework of standard dynamical theory of two-component (stars and dark matter) spherical systems with different velocity distributions. We compare the predicted velocity dispersion profiles to observations in the case of Fornax and Draco dwarfs. For isotropic models the dark haloes with cores are found to fit the data better than those with cusps. Anisotropic models are studied by fitting two parameters, dark mass and velocity anisotropy, to the data. In this case all profiles yield good fits, but the steeper the cusp of the profile, the more tangential is the velocity distribution required to fit the data. To resolve this well-known degeneracy of density profile versus velocity anisotropy, we obtain predictions for the kurtosis of the line-of-sight velocity distribution for models found to provide best fits to the velocity dispersion profiles. It turns out that profiles with cores typically yield higher values of kurtosis which decrease more steeply with distance than the cuspy profiles, which will allow us to discriminate between the profiles once the kurtosis measurements become available. We also show that with present quality of the data the alternative explanation of velocity dispersions in terms of Modified Newtonian Dynamics cannot yet be ruled out.  相似文献   

16.
If X , which follows a power-law distribution, is observed subject to Gaussian measurement error e , then   X + e   is distributed as the convolution of the power-law and Gaussian distributions. Maximum-likelihood estimation of the parameters of the two distributions is considered. Large-sample formulae are given for the covariance matrix of the estimated parameters, and implementation of a small-sample method (the jackknife) is also described. Other topics dealt with are tests for goodness of fit of the posited distribution, and tests whether special cases (no measurement errors or an infinite upper limit to the power-law distribution) may be preferred. The application of the methodology is illustrated by fitting convolved distributions to masses of giant molecular clouds in M33 and the Large Magellanic Cloud (LMC), and to H  i cloud masses in the LMC.  相似文献   

17.
Evidence of jet precession in many galactic and extragalactic sources has been reported in the literature. Much of this evidence is based on studies of the kinematics of the jet knots, which depends on the correct identification of the components to determine their respective proper motions and position angles on the plane of the sky. Identification problems related to fitting procedures, as well as observations poorly sampled in time, may influence the follow-up of the components in time, which consequently might contribute to a misinterpretation of the data. In order to deal with these limitations, we introduce a very powerful statistical tool to analyse jet precession: the cross-entropy method for continuous multi-extremal optimization. Only based on the raw data of the jet components (right ascension and declination offsets from the core), the cross-entropy method searches for the precession model parameters that better represent the data. In this work we present a large number of tests to validate this technique, using synthetic precessing jets built from a given set of precession parameters. With the aim of recovering these parameters, we applied the cross-entropy method to our precession model, varying exhaustively the quantities associated with the method. Our results have shown that even in the most challenging tests, the cross-entropy method was able to find the correct parameters within a 1 per cent level. Even for a non-precessing jet, our optimization method could point out successfully the lack of precession.  相似文献   

18.
The photometric light modulation of ζ Andromedae originates from the distorted geometry of the primary, and additionally, from spots of which parameters (temperature, size, location) are variable in time. We present spot modelling results for six two-colour light curves which show that spots preferably appear on the stellar surface towards the companion star and opposite to it, where the distortion also causes dimming. Therefore, simple fitting of the measured data for the ellipticity effect does not yield correct result. Instead, ellipticity calculated from exact stellar parameters should be removed from the data to get reliable spotted light curves.  相似文献   

19.
We describe a method of determining the system parameters in non-eclipsing interacting binaries. We find that the extent to which an observer sees the shape of the Roche lobe of the secondary star governs the amount of distortion of the absorption-line profiles. The width and degree of asymmetry of the phase-resolved absorption-line profiles show a characteristic shape, which depends primarily on the binary inclination and gravity darkening exponent. We show that, in principle, by obtaining high spectral and time resolution spectra of quiescent cataclysmic variables or low-mass X-ray binaries in which the mass-losing star is visible, fitting the shape of absorption-line profiles will allow one to determine not only the mass function of the binary, but also the binary inclination and hence the mass of the binary components.  相似文献   

20.
We estimate the distribution of intrinsic shapes of APM galaxy clusters from the distribution of their apparent shapes. We measure the projected cluster ellipticities using two alternative methods. The first method is based on moments of the discrete galaxy distribution while the second is based on moments of the smoothed galaxy distribution. We study the performance of both methods using Monte Carlo cluster simulations covering the range of APM cluster distances and including a random distribution of background galaxies. We find that the first method suffers from severe systematic biases, whereas the second is more reliable. After excluding clusters dominated by substructure and quantifying the systematic biases in our estimated shape parameters, we recover a corrected distribution of projected ellipticities. We use the non-parametric kernel method to estimate the smooth apparent ellipticity distribution, and numerically invert a set of integral equations to recover the corresponding distribution of intrinsic ellipticities under the assumption that the clusters are either oblate or prolate spheroids. The prolate spheroidal model fits the APM cluster data best.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号