首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 406 毫秒
1.
We present a revised metallicity distribution of dwarfs in the solar neighbourhood. This distribution is centred on solar metallicity. We show that previous metallicity distributions, selected on the basis of spectral type, are biased against stars with solar metallicity or higher. A selection of G-dwarf stars is inherently biased against metal-rich stars and is not representative of the solar neighbourhood metallicity distribution. Using a sample selected on colour, we obtain a distribution where approximately half the stars in the solar neighbourhood have metallicities higher than [Fe/H]=0 . The percentage of mid-metal-poor stars ([Fe/H]<−0.5) is approximately 4 per cent, in agreement with present estimates of the thick disc.
In order to have a metallicity distribution comparable to chemical evolution model predictions, we convert the star fraction to mass fraction, and show that another bias against metal-rich stars affects dwarf metallicity distributions, due to the colour (or spectral type) limits of the samples. Reconsidering the corrections resulting from the increasing thickness of the stellar disc with age, we show that the simple closed-box model with no instantaneous recycling approximation gives a reasonable fit to the observed distribution. Comparisons with the age–metallicity relation and abundance ratios suggest that the simple closed-box model may be a viable model of the chemical evolution of the Galaxy at solar radius.  相似文献   

2.
The analysis of the kinematics of solar neighbourhood stars shows that the low- and high-metallicity tails of the thin disc are populated by objects which orbital properties suggest an origin in the outer and inner Galactic disc, respectively. Signatures of radial migration are identified in various recent samples, and are shown to be responsible for the high-metallicity dispersion in the age–metallicity distribution. Most importantly, it is shown that the population of low-metallicity wanderers of the thin disc (−0.7 < [Fe/H] < −0.3 dex) is also responsible for the apparent hiatus in metallicity with the thick disc (which terminal metallicity is about −0.2 dex). It implies that the thin disc at the solar circle has started to form stars at about this same metallicity. This is also consistent with the fact that 'transition' objects, which have α-element abundance intermediate between that of the thick and thin discs, are found in the range [−0.4, −0.2] dex. Once the metal-poor thin disc stars are recognized for what they are – wanderers from the outer thin disc – the parenthood between the two discs can be identified on stars genuinely formed at the solar circle through an evolutionary sequence in [α/Fe] and [Fe/H]. Another consequence is that stars that can be considered as truly resulting of the chemical evolution at the solar circle have a metallicity restricted to about [−0.2, +0.2] dex, confirming an old idea that most chemical evolution in the Milky Way have preceded the thin disc formation.  相似文献   

3.
The initial mass function (IMF) in the solar neighbourhood is determined on the basis of a recently derived history of the star formation rate (SFR) which shows the presence of a star formation burst about 8 Gyr ago. The observed present-day mass function (PDMF) is considered, and the average past distribution of stars of a given mass is estimated. Two cases are considered, namely (i) constant SFR, and (ii) variable SFR as derived from the new metallicity distribution of G dwarfs. The resulting IMF is compared with previous determinations by Scalo and Kroupa et al., and the variation with stellar mass of the slope of the IMF is compared with reference determinations in the literature.  相似文献   

4.
The current paper investigates how the empirical, G-dwarf metallicity distribution constrains simple, comoving models of chemical evolution. In doing this, the application of the models to a data sample, performed in a previous paper, is refined and extended. The key idea is that (i) different star formation rates with different mass spectra take place in different phases of evolution, i.e. contraction and equilibrium, and (ii) disk formation begins at a time t = Td and ends at t = Tc, which marks the transition from contraction to equilibrium. In this view, the lowest-metallicity point of the empirical, differential distribution, consistent with a linear fit, is related to the beginning of disk formation, and an apparent discontinuity point to the transition from contraction to equilibrium. In addition, different linear fits hold on the left (early distribution) and on the right (late distribution) of the discontinuity point. Models consistent with the empirical, G-dwarf metallicity distribution are related to linear fits on the early and late side. Homologous solutions during the equilibrium phase are analysed in detail with respect to changes in Tc and Ta, the age of the Galaxy. Then we are left with a single free parameter which is relevant to the chemical evolution, i.e. the mass spectrum exponent during the equilibrium phase. The allowed values for the other parameters, thought as a function of the above mentioned one, are plotted for each case. A Salpeter mass spectrum exponent, p = −2.35, is ruled out by the theoretical, lower stellar mass limit, contrary to a Scalo mass spectrum exponent, p = −2.90, in contrast with previous literature. The reasons for this discrepancy are discussed. Our results are marginally consistent with a same initial mass function during the contraction and equilibrium phase, but in this case the disk mass fraction is of the is same order, or less, than the halo mass fraction. It is also investigated how the empirical age-metallicity relation constrains the duration of the contraction phase, for a reasonable upper limit of Ta. Keeping in mind that the empirical, G-dwarf metallicity distribution has not been corrected for the large cosmic scatter shown by the empirical, age-metallicity relation, we find a duration of disk formation, TcTd = 1.07–1.5 Gyr, by a factor 3–5 less than it is found by use of simple infall models. The reasons of this difference are explained. The idea of a massive, white dwarf halo, which seems to be indicated by microlensing experiments, is ruled out by the empirical, G-dwarf metallicity distribution, in the light of the current model and provided the solar neighbourhood is a typical region of the Galaxy. More refined models involving e.g., the relax of instantaneous recycling would change our results, but the trend is expected to be only slightly altered.  相似文献   

5.
为解释著名的G矮星问题,提出银河系化学演化的三成分模型,即由银晕、厚盘和薄盘所构成的演化模型.相邻演化阶段间隔着一个快速坍缩过程.对不同星族成分的演化过程分别进行模拟,并在总体上得到一个太阳附近区域的G矮星丰度分布函数.检验了三种不同的模型:初始富化模型、比例生成模型和坍缩模型.利用最小二乘拟合得到最佳模型的参数.结果表明,太阳附近区域的化学演化受物质交换的影响较小,至少在银河系演化的晚期,可将太阳附近区域看作封闭系统.同时,单位质量中新合成的重元素比例对三种恒星成分可分别近似为常数,其差别则说明不同星族恒星的初始质量函数存在着显著差异.  相似文献   

6.
The oxygen abundance distribution in solar neighbourhood halo subdwarfs is deduced, using two alternative, known empirical relationships, involving the presence or the absence of [O/Fe] plateau for low [Fe/H] values, from a sample of 372 kinematically selected halo stars, for which the iron abundance distribution has been determined by Ryan & Norris (1991). The data are interpreted by a simple, either homogeneous or inhomogeneous model of chemical evolution, using an updated value of the solar oxygen abundance. The effect of changing the solar oxygen abundance, the power‐law exponent in the initial mass function, and the rate of oxygen nucleosyntesis, keeping the remaining input parameters unchanged, is investigated, and a theorem is stated. In all cases, part of the gas must necessarily be inhibited from forming stars, and no disk contamination has to be advocated for fitting the empirical oxygen abundance distribution in halo subdwarfs of the solar neighbourhood (EGD). Then a theorem is stated, which allows a one‐to‐one correspondence between simple, homogeneous models with and without inhibited gas, related to same independent parameters of chemical evolution, except lower stellar mass limit, real yield, and inhibition parameter. The mutual correlations between the latter parameters are also specified. In addition the starting point, and the point related to the first step, of the theoretical distribution of oxygen abundance (TGD) predicted by simple, inhomogeneous models, is calculated analytically. The mean oxygen abundance of the total and only inhibited gas, respectively, are also determined. Following the idea of a universal, initial mass function (IMF), a power‐law with both an exponent p = 2.9, which is acceptably close to Scalo IMF for mm, and an exponent p = 2.35, i.e. Salpeter IMF, have been considered. In general, both the age‐metallicity relationship and the empirical distribution of oxygen abundance in G dwarfs of the disk solar neighbourhood, are fitted by power‐law IMF exponents in the range 2.35 ≤ p ≤ 2.9. Acceptable models predict about 15% of the total mass in form of long‐lived stars and remnants, at the end of halo evolution, with a mean gas oxygen abundance which is substantially lower than the mean bulge and initial disk oxygen abundance. To avoid this discrepancy, either the existence of a still undetected, baryonic dark halo with about 15% of the total mass, or an equal amount of gas loss during bulge and disk formation, is necessary. The latter alternative implies a lower stellar mass limit close to 0.2 m, which is related to a power‐law IMF exponent close to 2.77. Acceptable models also imply a rapid halo formation, mainly during the first step, Δt = 0.5 Gyr, followed by a period (three steps) where small changes occur. Accordingly, statistical fluctuations are found to produce only minor effects on the evolution.  相似文献   

7.
A three-component chemical evolution model of the Galaxy is presented, which we believe will cast a new light on the G-dwarf problem. The model is based on a scenario of the Galaxy consisting of three major evolutionary phases: halo, thick disk and thin disk, separated by two short interludes of rapid collapse. The evolution of different stellar populations are treated separately, the combination of which yields an overall metallicity distribution function for the solar neighbourhood. We tested three different models using the same set of basic equations: the “prompt initial enrichment” (PIE) model, the “proportional yield” (PPY) model and the “collapse” (CLP) model. Best-fit parameters are derived. The results show that the different populations have remarkably different IMFs, while mass exchange has only minimally affected the chemical evolution in the solar vicinity, so that the solar vicinity can be regarded as a closed system, at least in the late stage of the Galactic evolution.  相似文献   

8.
Models of the chemical evolution of our Galaxy are extended to include radial migration of stars and flow of gas through the disc. The models track the production of both iron and α-elements. A model is chosen that provides an excellent fit to the metallicity distribution of stars in the Geneva–Copenhagen survey (GCS) of the solar neighbourhood and a good fit to the local Hess diagram. The model provides a good fit to the distribution of GCS stars in the age–metallicity plane, although this plane was not used in the fitting process. Although this model's star formation rate is monotonically declining, its disc naturally splits into an α-enhanced thick disc and a normal thin disc. In particular, the model's distribution of stars in the ([O/Fe], [Fe/H]) plane resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's ridge line is entirely due to stellar migration, and there is the characteristic variation of stellar angular momentum along it that has been noted by Haywood in survey data. Radial mixing of stellar populations with high  σ z   from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of  σ z   with age than predicted by theory. The metallicity gradient in the interstellar medium is predicted to be steeper than in earlier models, but appears to be in good agreement with data for both our Galaxy and external galaxies. The models are inconsistent with a cut-off in the star formation rate at low gas surface densities. The absolute magnitude of the disc is given as a function of time in several photometric bands, and radial colour profiles are plotted for representative times.  相似文献   

9.
We show that the explicit assumption of a chemically inhomogeneous interstellar medium allows a better reproduction of the metallicity distribution of G-dwarfs in the solar neighbourhood. The inhomogeneity is considered by assuming that at any time stars are born with a spread in their metallicities, the spread being a Gaussian in the logarithm of the metallicity around the mean metallicity of that epoch. We show that for various simple models of chemical evolution, the fit to the G-dwarf metallicity curve improves considerably once the above assumption is applied. We show that the parameters obtained from the fitting also give acceptable predictions for the age-metallicity relation. We also find that if we use a G-dwarf metallicity function corrected for the scale height inflation of stars, the conventional models of chemical evolution cannot match the shape of the curve, at least under the instantaneous recycling approximation applied to a chemically homogeneous ISM. Under the inhomogeneous ISM approximation, the predicted shapes are found to be better, though not totally satisfactory.  相似文献   

10.
The mean metallicity of the Milky Way thin disc in the solar neighbourhood is still a matter of debate, and we recently proposed an upward revision. Our star sample was drawn from a set of solar neighbourhood dwarfs with photometric metallicities. In a very recent study, it has been suggested that our metallicity calibration, based on Geneva photometry, is biased. We show here that the effect detected is not a consequence of our adopted metallicity scale, and we confirm that our findings are robust. On the contrary, the application to Strömgren photometry of the Schuster & Nissen metallicity scale is problematic. Systematic discrepancies of  ∼0.1–0.3 dex  affect the photometric metallicity determination of metal-rich stars, on the colour interval  0.22 < b − y < 0.59  , i.e. including F and G stars. For F stars, it is shown that this is a consequence of a mismatch between the standard sequence   m 1( b − y )  of the Hyades used by Schuster & Nissen to calibrate their metallicity scale, and the system of Olsen. It means that although the calibration of Schuster & Nissen and Olsen's photometry are intrinsically correct, they are mutually incompatible for metal-rich F-type stars. For G stars, the discrepancy is most probably the continuation of the same problem, albeit worsened by the lack of spectroscopic calibrating stars. A corrected calibration is proposed that renders the calibration of Schuster & Nissen applicable to the catalogues of Olsen. We also give a simpler calibration referenced to the Hyades sequence, valid over the same colour and metallicity ranges.  相似文献   

11.
The current Swift sample of gamma-ray bursts (GRBs) with measured redshifts allows us to test the assumption that GRBs trace star formation in the Universe. Some authors have claimed that the rate of GRBs increases with cosmic redshift faster than the star formation rate, whose cause is not yet known. In this paper, I investigate the possibility of interpreting the observed discrepancy between the GRB rate history and the star formation rate history using cosmic metallicity evolution. I am motivated by the observation that cosmic metallicity evolves with redshift and GRBs tend to occur in low-metallicity galaxies. First, I derive a star formation history up to redshift   z = 7.4  from an updated sample of star formation rate densities. This is obtained by adding the new ultraviolet measurements of Bouwens et al. and the new ultraviolet and infrared measurements of Reddy et al. to the existing sample compiled by Hopkins & Beacom. Then, adopting a simple model for the relation between GRB production and the cosmic metallicity history as proposed by Langer & Norman, I show that the observed redshift distribution of the Swift GRBs can be reproduced with good accuracy. Although the results are limited by the small size of the GRB sample and the poorly understood selection biases in detection and localization of GRBs and in redshift determination, they suggest that GRBs trace both star formation and metallicity evolution. If the star formation history can be accurately measured with other approaches, which is presumably achievable in the near future, it will be possible to determine the cosmic metallicity evolution using the study of the redshift distribution of GRBs.  相似文献   

12.
We show that there is a relationship between the age excess, defined as the difference between the stellar isochrone and chromospheric ages, and the metallicity as measured by the index [Fe/H] for late-type dwarfs. The chromospheric age tends to be lower than the isochrone age for metal-poor stars, and the opposite occurs for metal-rich objects. We suggest that this could be an effect of neglecting the metallicity dependence of the calibrated chromospheric emission–age relation. We propose a correction to account for this dependence. We also investigate the metallicity distributions of these stars, and show that there are distinct trends according to the chromospheric activity level. Inactive stars have a metallicity distribution which resembles the metallicity distribution of solar neighbourhood stars, while active stars appear to be concentrated in an activity strip on the log  R 'HK × [Fe/H] diagram. We provide some explanations for these trends, and show that the chromospheric emission–age relation probably has different slopes on the two sides of the Vaughan–Preston gap.  相似文献   

13.
The evolution of the Solar neighbourhood is followed using a unique, magnitude complete and kinematically unbiased sample of 14,000 F, G, and K dwarfs. Metallicity, age, space motion and galactic orbits have been determined for all stars. The result is a detailed view of the complex evolution of the local Milky Way, which must be matched by any model for the chemical and dynamical evolution of the Galactic disk. E.g., such models must explain the shape and large scatter in the age-metallity relation as well as the overall metallicity distribution; the evolution of stellar kinematics with age; the distribution of stars in velocity space; and the contributions from the thick and thin disks to all these relations. This revised version was published online in August 2006 with corrections to the Cover Date.  相似文献   

14.
The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for  <10−4 Z  . With initial mass  <0.1 M  , if fragments are ejected in embryo from the host clouds by multibody interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.  相似文献   

15.
We use nearby K dwarf stars to measure the helium-to-metal enrichment ratio  Δ Y /Δ Z   , a diagnostic of the chemical history of the solar neighbourhood. Our sample of K dwarfs has homogeneously determined effective temperatures, bolometric luminosities and metallicities, allowing us to fit each star to the appropriate stellar isochrone and determine its helium content indirectly. We use a newly computed set of Padova isochrones which cover a wide range of helium and metal content.
Our theoretical isochrones have been checked against a congruous set of main-sequence binaries with accurately measured masses, to discuss and validate their range of applicability. We find that the stellar masses deduced from the isochrones are usually in excellent agreement with empirical measurements. Good agreement is also found with empirical mass-luminosity relations.
Despite fitting the masses of the stars very well, we find that anomalously low helium content (lower than primordial helium) is required to fit the luminosities and temperatures of the metal-poor K dwarfs, while more conventional values of the helium content are derived for the stars around solar metallicity.
We have investigated the effect of diffusion in stellar models and the assumption of local thermodynamic equilibrium (LTE) in deriving metallicities. Neither of these is able to resolve the low-helium problem alone and only marginally if the cumulated effects are included, unless we assume a mixing-length which is strongly decreasing with metallicity. Further work in stellar models is urgently needed.
The helium-to-metal enrichment ratio is found to be  Δ Y /Δ Z = 2.1 ± 0.9  around and above solar metallicity, consistent with previous studies, whereas open problems still remain at the lowest metallicities. Finally, we determine the helium content for a set of planetary host stars.  相似文献   

16.
Based on a new metallicity distribution for G dwarfs within 25 pc of the sun this paper discusses in detail four galactic chemical evolution models in an inhomogeneous interstellar medium. It is shown that both the so-called simple model and the collapse model deviate greatly from the observations, while the two models labelled PIE (Prompt Initial Enrichment) and PPY (Proportional Yield) not only give much better fit, but also alleviate the so-called G-dwarf problem. The indications are that the outer halo has little influence on the chemical evolution near the Sun.  相似文献   

17.
Infall models for the evolution of the local galactic disk were studied and confronted with a large number of observational constraints from the solar vicinity, inclusive of the white dwarf luminosity function. The models are characterized as follows: 1. The key-functions (SFR, IMF, gas infall rate) are not prescribed by simple laws, but are directly derived from observational constraints. 2. A scatter in the metallicity at fixed age is considered which partly reflects inhomogeous chemical evolution. 3. Special attention is drawn to the internal consistency of the models. 4. In addition to infall of low-metallicity gas, metal-enriched outflows are allowed. The “best” model is characterized by a disk age of ≈︁ 12 Gyr, a SFR which is decreasing over the first half and is nearly constant over the second half of the disk evolution, and by a similar temporal run of the gas infall rate. Moderate metal-enriched outflow can not be excluded.  相似文献   

18.
A sample containing 1 026 stars of spectral types F, G, and K, mainly dwarfs, from the solar neighbourhood with available space velocities and metallicities is treated. The treatment comprises a statistical analysis of the metallicity and velocity data and calculation of galactocentric orbits. Sample stars identified as members of the galactic halo are detached from the rest of the sample based on the values of their metallicities, velocity components and galactocentric orbits. In identifying halo stars a new, kinematical, criterion is proposed. Except one, these halo stars are the metal‐poorest ones in the sample. Besides, they have very high velocities with respect to LSR. On the other hand, the separation between the thin disc and thick one is done statistically based on LSR space velocities, membership probability (Schwarzschild distribution with assumed parameters) and galactocentric orbits. In the metallicity these two groups are not much different. For each of the three subsamples the mean motion and velocity ellipsoid are calculated. The elements of the velocity ellipsoids agree well with the values found in the literature, especially for the thin disc. The fractions of the subsystems found for the present sample are: thin disc 93%, thick disc 6%, halo 1%. The sample stars established to be members of the thin disc are examined for existence of star streams. Traces of both, known and unknown, star streams are not found (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)  相似文献   

19.
This paper has two parts: one about observational constraints related to the empirical differential oxygen abundance distribution (EDOD), and the other about inhomogeneous models of chemical evolution, in particular the theoretical differential oxygen abundance distribution (TDOD). In the first part, the EDOD is deduced from subsamples related to two different samples involving (i) N=532 solar neighbourhood (SN) stars within the range, −1.5<[Fe/H]<0.5, for which the oxygen abundance has been determined both in presence and in absence of the local thermodynamical equilibrium (LTE) approximation (Ramirez et al. in Astron. Astrophys. 465:271, 2007); and (ii) N=64 SN thick disk, SN thin disk, and bulge K-giant stars within the range, −1.7<[Fe/H]<0.5, for which the oxygen abundance has been determined (Melendez et al. in Astron. Astrophys. 484:L21, 2008). A comparison is made with previous results implying use of [O/H]–[Fe/H] empirical relations (Caimmi in Astron. Nachr. 322:241, 2001b; New Astron. 12:289, 2007) related to (iii) 372 SN halo subdwarfs (Ryan and Norris in Astron. J. 101:1865, 1991); and (iv) 268 K-giant bulge stars (Sadler et al. in Astron. J. 112:171, 1996). The EDOD of the SN thick + thin disk is determined by weighting the mass, for assumed SN thick to thin disk mass ratio within the range, 0.1–0.9. In the second part, inhomogeneous models of chemical evolution for the SN thick disk, the SN thin disk, the SN thick + thin disk, the SN halo, and the bulge, are computed assuming the instantaneous recycling approximation. The EDOD data are fitted, to an acceptable extent, by their TDOD counterparts with the exception of the thin or thick + thin disk, where two additional restrictions are needed: (i) still undetected, low-oxygen abundance thin disk stars exist, and (ii) a single oxygen overabundant star is removed from a thin disk subsample. In any case, the (assumed power-law) stellar initial mass function (IMF) is universal but gas can be inhibited from, or enhanced in, forming stars at different rates with respect to a selected reference case. Models involving a strictly universal IMF (i.e. gas neither inhibited from, nor enhanced in, forming stars with respect to a selected reference case) can also reproduce the data to an acceptable extent. Our main conclusions are (1) different models are necessary to fit the (incomplete) halo sample, which is consistent with the idea of two distinct halo components: an inner, flattened halo in slow prograde rotation, and an outer, spherical halo in net retrograde rotation (Carollo et al. in Nature 450:1020, 2007); (2) the oxygen enrichment within the inner SN halo, the SN thick disk, and the bulge, was similar and coeval within the same metallicity range, as inferred from observations (Prochaska et al. in Astron. J. 120:2513, 2000); (3) the fit to thin disk data implies an oxygen abundance range similar to its thick disk counterpart, with the extension of conclusion (2) to the thin disk, and the evolution of the thick + thin disk as a whole (Haywood in Mon. Not. R. Astron. Soc. 388:1175, 2008) cannot be excluded; (4) leaving outside the outer halo, a fit to the data related to different environments is provided by models with a strictly universal IMF but different probabilities of a region being active, which implies different global efficiencies of the star formation rate; (5) a special case of stellar migration across the disk can be described by models with enhanced star formation, where a fraction of currently observed SN stars were born in situ and a comparable fraction is due to the net effect of stellar migration, according to recent results based on high-resolution N-body + smooth particle hydrodynamics simulations (Roškar et al. in Astrophys. J. Lett. 684:L79, 2008).  相似文献   

20.
The formation of a disk galaxy within a slowly growing dark halo is simulated with a new chemo-dynamical model. The model describes the evolution of the stellar populations, the multi-phase ISM and all important interaction. I find, that the galaxy forms radially from inside-out and vertically from top-to-bottom. The derived stellar age distributions show that the inner halo is the oldest component, followed by the outer halo, the triaxial bulge, the halo-disk transition region and the disk. Despite the still idealized model, the final galaxy resembles present-day disk galaxies in many aspects. In particular, the stellar metallicity distribution in the halo of the model resembles the one of M31. The bulge in the model shows, at least two stellar subpopulations, an early collapse population and a population that formed later out of accreted disk mass. In the stellar metallicity distribution of the disk, I find a pronounced ‘G-dwarf problem’ which is the result of a pre-enrichment of the disk ISM with metal-rich gas from the bulge. This revised version was published online in September 2006 with corrections to the Cover Date.  相似文献   

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

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