首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The evolution of an 0.6 m stellar model during core helium burning is presented. Following the off-center ignition of helium in the core flash, the star remains on the red giant branch for > 106 years, undergoing twelve additional flashes. After leaving the giant branch, the star evolves on the horizontal branch for 8.15×107 years before re turning to the giant branch and undergoing strong helium-shell flashes. The implications for horizontal branch and RR Lyrae stars are discussed.  相似文献   

2.
Cooling of neutron stars is calculated using an exact stellar evolution code. The full general relativistic version of the stellar structure equations are solved, with the best physical input currently available. For neutron stars with a stiff equation of state, we find that the deviation from the isothermality in the interior is significant and that it takes at least a few thousand years to reach the isothermal state. By comparing the most recent theoretical and observational results, we conclude that for Cas A, SN1006, and probably Tycho, standard cooling is inconsistent with the results from the Einstein Observatory, if neutron stars are assumed to be present in these objects. On the other hand, the detection points for RCW103 and the Crab are consistent with these theoretical results.On leave from Department of Physics, Ibaraki University, Japan  相似文献   

3.
Conclusion The EXOSAT observation of PSR 1055-52 supports strongly the idea of a thermally radiating neutron star, either still cooling from its initially high temperature or reheated by unpinning of vortex lines in its superfluid interior. The data do not confirm the Cheng & Helfand conclusions about a synchrotron nebula with a power law spectrum. The discrepancy in the data is, nevertheless, only a 3.5 effect, as no ME data could be used to get a firm upper limit of the flux above 1 keV. With regard to the outstanding theoretical interest in the structure of neutron stars we hope for further observations to clarify this question.  相似文献   

4.
New ultraviolet (1300 A, 3400 A),HST FOC observations have been used to derive the UV color-magnitude diagram (CMD) of R136, with the main scientific goal of studying the upper end of the stellar mass function at ultraviolet wavelengths where the color degeneracy encountered in visual CMDs is less severe. The CMD has been compared to a set of theoretical isochrones, which have been computed using the latest generation of evolutionary models and model atmospheres for early type stars. Wolf-Rayet stars are included. Comparison of theTheoretical andobserved CMD suggests that there are no stars brighter than M130–11. We use the observed main sequence turn-off and the known spectroscopic properties of the stellar population to derive constraints on the most probable age of R136. The presence of WNL stars and the lack of red supergiants suggests a most likely age of 3±1 Myr. A theoretical isochrone of 3±1 Myr is consistent with the observed stellar content of R136 if the most massive stars have initial masses around 50 M.Bases on Observations with the NASA/ESA Hubble Space Telescope, obtained at the STScI, which is operated by AURA, Inc., under NASA contract NAS5-26555.Astrophysics Division, Space Science Department, ESA  相似文献   

5.
A newly formed neutron star in a supernova finds itself in a dense environment, in which the gravitational energy of accreting matter can be lost to neutrinos. For the conditions in SN 1987A, 0.1M may have fallen back onto the central neutron star on a timescale of hours after the explosion, after which the accretion rate is expected to drop sharply. Radiation is trapped in the flow until the mass accretion rate drops to 2×10–4 M yr–1 at which point radiation can begin to escape from the shocked envelope at an Eddington limit luminosity. Between this neutrino limit and the Eddington limit, 3×10–8 M yr–1, there are no steady, spherical solutions for neutron star accretion. SN 1987A should have reached the neutrino limit within a year of the explosion; the current lack of an Eddington luminosity can be attributed to black hole formation or to a clearing of the neutron star envelope. There is no evidence for newly formed neutron stars in supernovae. Radio supernovae, which were initially interpreted as pulsar activity, probably involve circumstellar interaction; SN 1993J shows especially good evidence for outer shock phenomena.  相似文献   

6.
It is commonly accepted that candidates for very high energy -ray sources are neutron stars, binary systems, black holes etc. Close binary systems containing a normal hot star and a neutron star (or a black hole) form an important class of very high energy -ray sources. Such systems are variable in any region of the electromagnetic spectrum and they enable us to study various stages of stellar evolution, accretion processes, mechanisms of particle acceleration, etc. Phenomena connected with this class of very high energy -ray sources are discussed. Particular emphasis has been placed on the TeV energy region.  相似文献   

7.
Thanks to remarkable new tools, such as the Goddard High Resolution Spectrograph (GHRS) on board the HST and the EUVE spectrometer on the interstellar side, and Ulysses particle detectors on the heliospheric side, it is possible now to begin to compare abundances and physical properties of the interstellar matter outside the heliosphere (from absorption features in the stellar spectra), and inside the heliosphere (from in situ or remote detection of the interstellar neutrals or their derivatives, the pick-up ions or the Anomalous Cosmic Rays detected by the two Voyager spacecraft).Ground-based and UV spectra of nearby stars show that the Sun is located between two volumes of gas of different heliocentric velocities V and temperatures T (see also Linsky et al, this issue). One of these clouds has the same velocity (V= 25.6 km s–1 from = 255 and =8) and temperature (6700 K) as the heliospheric helium of interstellar origin probed by Ulysses, and is certainly surrounding our star (and then the Local Interstellar Cloud or LIC). This Identification allows comparisons between interstellar constituents on both sides of the heliospheric interface.Ly-alpha background data (absorption cell and recent HST-GHRS spectra) suggest that the heliospheric neutral H velocity is smaller by 5–6 km s–1 than the local cloud velocity, and therefore that H is decelerated at its entrance into the heliosphere, in agreement with interaction models between the heliosphere and the ISM which include the coupling with the plasma. This is in favor of a non negligible electron density (at least 0.05 cm3). There are other indications of a rather large ionization of the ambient ISM, such as the ionization equilibrium of interstellar magnesium and of sodium. However the resulting range for the plasma density is still broad.The heliospheric neutral hydrogen number density (0.08–0.16 cm–3) is now less precisely determined than the helium density (0.013–0.017 cm–3, see Gloeckler, Witte et al, Mobius, this issue). The comparison between the neutral hydrogen to neutral helium ratios in the ISM (recent EUVE findings) and in the heliosphere, suggests that 15 to 70% of H does not enter the heliosphere. The comparison between the interstellar oxygen relative abundance (with respect to H and He) in the ISM and the heliospheric abundance deduced from pick-up ions is also in favor of some filtration, and thus of a non-negligible ionization.For a significant ISM plasma density, one expects a Hydrogen wall to be present as an intermediate state of the interstellar H around the interface between inside and outside. Since 1993, the two UVS instruments on board Voyager 1 and 2 indeed reveal clearly the existence of an additional Ly-alpha emission, probably due to a combination of light from the compressed H wall, and from a galactic source. On the other hand, the decelerated and heated neutral hydrogen of this H wall has recently been detected in absorption in the spectra of nearby stars (see Linsky, this issue).  相似文献   

8.
Recent studies suggest that when magnetohydrodynamic (MHD) turbulence is excited by stirring a plasma at large scales, the cascade of energy from large to small scales is anisotropic, in the sense that small-scale fluctuations satisfy the inequality k k , where k and k are, respectively, the components of a fluctuations wave vector and to the background magnetic field. Such anisotropic fluctuations are very inefficient at scattering cosmic rays. Results based on the quasilinear approximation for scattering of cosmic rays by anisotropic MHD turbulence are presented and explained. The important role played by molecular-cloud magnetic mirrors in confining and isotropizing cosmic rays when scattering is weak is also discussed.  相似文献   

9.
We present helium and CNO isotopic yields for massive mass-losing stars in the initial mass range 15M M i 50M . We investigate their dependence on assumptions about mass loss rates, internal mixing processes, and metallicity, and specify the contributions from stellar winds and from supernova ejecta.  相似文献   

10.
Since the baryon-to-photon ratio 10 is in some doubt at present, we ignore the constraints on 10 from big bang nucleosynthesis (BBN) and fit the three key cosmological parameters (h, M, 10) to four other observational constraints: Hubble parameter (ho), age of the universe (to), cluster gas (baryon) fraction (fo fGh3/2), and effective shape parameter (o). We consider open and flat CDM models and flat CDM models, testing goodness of fit and drawing confidence regions by the 2 method. CDM models with M = 1 (SCDM models) are accepted only because we allow a large error on ho, permitting h < 0.5. Open CDM models are accepted only for M 0.4. CDM models give similar results. In all of these models, large 10 ( 6) is favored strongly over small 10 ( 2), supporting reports of low deuterium abundances on some QSO lines of sight, and suggesting that observational determinations of primordial 4He may be contaminated by systematic errors. Only if we drop the crucial o constraint are much lower values of M and 10 permitted.  相似文献   

11.
We review the possible evolutionary paths from massive stars to explosive endpoints as various types of supernovae associated with Population I and hence with massive stars: Type II-P, Type II-L, Type Ib, Type Ic, and the hybrid events SN 1987K and SN 1993J. We identify SN 1954A as another hybrid event from the evidence for both H and He in its spectrum with velocities nearly the same as SN 1983J. Evidence for ejected56Ni mass of 0.07 M suggests that SN II-P underwent standard iron core collapse, not collapse of an O–Ne–Mg core nor thermonuclear explosion of a C–O core. Most SN II-P presumably arise in single stars or wide binaries of 10–20 M. There may be indirect evidence for duplicity in some cases in the form of strong Ba II lines, such as characterized SN 1987A. SN II-L are recognizably distinct from typical SN II-P and must undergo a significantly different evolution. Despite indications that SN II-L have small envelopes that may be helium enriched, they are also distinct from events like SN 1993J that must have yet again a different evolution. The SN II-L that share a common Luminosity seem to have ejected a small nickel mass and hence may come from stars with O–Ne–Mg cores. The amount of nickel ejected by the exceptionally bright events, SN 1980K and SN 1979C, remains controversial. SN Ib require the complete loss of the H envelope, either to a binary companion or to a wind. The few identified have relatively large ejecta masses. It is not clear what evolutionary processes distinguish SN Ib's evolving in binary systems from hybrid events that retain some H in the envelope. SN Ic events are both H and He deficient. Binary models that can account for transfer of an extended helium envelope from low mass helium cores, 2 to 4 M, imply C–O core masses that are roughly consistent with that deduced from the ejecta mass plus a neutron star, 2 to 3 M. It is possible that the hybrid events are the result of Roche lobe overflow and that the pure events, SN Ib or SN Ic, result from common envelope evolution.  相似文献   

12.
Summary Using values of d, min, and max that Van Riper (1978) has found most promising for a hydrodynamic envelope ejection, we have shown that even a small amount of rotation in the initial core can stop its collapse before nuclear densities are reached. We expected i > 0.02 to produce significant deviations from a spherically symmetric collapse, but have found that i as much as ten times smaller than this will not allow the core to reach densities as high as in the spherical collapse. In no case, however, does the core flatten very much, nor does the value of become very large. Low final 's preclude the formation of an axisymmetric torus. They also indicate that deformation of an iron core into a triaxial configuration or fragmentation of the core during its collapse is an extremely unlikely event. (Note: Classically, must exceed 0.27 before a dynamic instability to non-axisymmetric perturbations is encountered.)The small degree of flattening of the core also suggests that the reduced moment of inertia I of the core will always be relatively small in magnitude and hence that the third time derivative of I, which is proportional to the energy emitted in gravity wave radiation, will not be very significant. Numerically calculated estimates of I- during some of these model evolutions supports this suspicion. If the min and used here are found to be realistic values after the detailed physics of the core collapse is well understood, it is clear that gravitational radiation from a core collapse will be difficult to measure.Finally, we should point out that it is the relatively large values of Ymin (near 4/3) combined with values of d near unity that (a) prevented the core from flattening significantly in these models and (b) prevented the core from reaching high configurations. If realistic values of either one (or both) of these parameters are found to be much smaller in more complete models of the core collapse, then the core will have to become flatter (and denser) before pressure gradients will support it along the rotation axis. All of the conclusions drawn here would be modified accordingly under those circumstances. It should also be noted that in general relativistic models, the critical for spherical collapse is somewhat larger than 4/3 (Van Riper, 1979). Therefore, we predict that when fully general relativistic core collapses are performed including rotation, a given choice of min and i will produce a slightly flatter and slightly denser core than the corresponding model that has been presented here.  相似文献   

13.
Present status of the theories for presupernova evolution and triggering mechanisms of supernova explosions are summarized and discussed from the standpoint of the theory of stellar structure and evolution. It is not intended to collect every detail of numerical results thus far obtained, but to extract physically clear-cut understanding from complexities of the numerical stellar models. For this purpose the evolution of stellar cores is discussed in a generalized fashion. The following types of the supernova explosions are discussed. The carbon deflagration supernova of intermediate mass star which results in the total disruption of the star. Massive star evolves into a supernova triggered by photo-dissociation of iron nuclei which results in a formation of a neutron star or a black hole depending on its mass. These two are typical types of the sueprnovae. Between them there remains a range of mass for which collapse of the stellar core is triggered by electron captures, which has been recently shown to leave a neutron star despite oxygen deflagration competing with the electron captures. Also discussed are combustion and detonation of helium or carbon which take place in accreting white dwarfs, and the collapse which is triggered by electron-pair creation in very massive stars.Appendix: Notations A mass number of atomic nucleus - B v(a, b) incomplete beta function - c p specific heat at constant pressure - c p sound velocity - c(sub) center of the star - E e mean energy of an electron captured by nucleus - E n nuclear energy release from unit mass of the nuclear fuel specified by n - E thr threshold energy (9.3) - E thr,0 energy difference between the ground states of daughter nucleus and parent nucleus (9.1) - E energy of gamma ray emitted from daughter nucleus (9.1) - E v mean energy of a neutrino emitted by electron capture (9.1) - f flatness parameter (2.17) - g local gravitational acceleration (2.16) - H atomic mass unit - H p scale height of pressure (2.22) - H (sub) hydrogen-burning shell - k Boltzmann constant - l mixing length of convection - L cr(M r ) local Eddington's critical luminosity (4.3) - L n integrated nuclear energy generation rate by nuclear fuel specified by n - L v neutrino luminosity - L v, cr(M r ) local Eddington's critical neutrino luminosity (11.2) - M (current) mass of a star - m M core mass contained interior to the carbon-burning shell - M Ch Chandrasekhar's limiting mass (9.6) - M H core mass contained interior to the hydrogen-burning shell - M He core mass contained interior to the helium-burning shell - M ms mass of a star at its zero-age min-sequence - M O core mass contained interior to the oxygen-burning shell - M r mass contained interior to a shell at r - M Si core mass contained interior to the silicon-burning shell - M WD mass of white dwarf (7.1) - M 0 normalization factor to the non-dimensional mass (3.3) - M 1 core mass (3.6) - N polytropic index between pressure and density (2.3) - n polytropic index between pressure and temperature (10.1) - N A Avogadro number - N ad adiabatic polytropic index - N e number of electrons in unit mass of matter - NSE nuclear statistical equilibrium - P pressure - ph (sub) photosphere - Q e mass fraction of the envelope exterior of the shell e (2.14) - R stellar radius - r radial distance of a shell - r 0 normalization factor to the non-dimensional radius (3.2) - s specific entropy - S i specific entropy of ions - T temperature - U homology invariant defined by (2.1) - u gas specific internal energy of gas - u rad energy of the radiation field per volume in which unit mass of gas is contained (6.4) - V homology invariant defined by (2.2) - def velocity of deflagration front (6.10) - X concentration by weight of hydrogen - Y concentration by weight of helium - Y e mole number of electrons in one gram of matter (9.7) - Y v mole number of neutrinos in one gram of matter - Z concentration by weight of the elements other than hydrogen and helium - z shock strength (6.6) - 1 (sub) usually denotes the core edge (2.13) - ratio of the mixing length to the scale height of pressure (l/H p ) - ratio of gas pressure to the total pressure - ratio of the specific heats - gD locus of singularity in U-V plane (2.5) - M(H p ) mass contained within unit scale height of pressure (4.4) - ec energy rate by electron captures (9.5) - n nuclear energy generation rate by the nuclear fuel specified by n - v neutrino loss rate - L v (D) neutrino loss rate excluding the neutrinos from the electron captures (9.4) - non-dimensional density (3.1) - P/, not the non-dimensional temperature (2.7) - W Weinberg's angle (5.8) - opacity - v neutrino opacity (11.2) - describes the effect of electron degeneracy in equation of state (2.19) - ec rate of electron capture - mean molecular weight - e mean molecular weight of electrons - e chemical potential of an electron excluding the rest mass (8.1) - i mean molecular weight of ions - non-dimensional radius (3.1) - non-dimensional pressure (3.1) - matter density - cr GR critical density above which the general relativistic instability sets in - cr critical density for reimplosion of the core by beta processes (Section 5) - ign density at the ignition - nse density above which the deflagrated matter results in NSE composition - e non-dimensional entropy of electron-per one electron in units of k(9.2) - ff timescale of free fall (6.2) - h (H p ) timescale of heat transport over unit scale height of pressure (4.4) - n nuclear timescale for a change in temperature (6.1) - non-dimensional mass (3.1) - e chemical potential of an electron in units of kT (8.1)  相似文献   

14.
The purpose of this work is to study the various -ray-production mechanisms in solar flares and to calculate the flux, the spectrum, and the decay curves of radiation. Using the continuity equation and taking into account the energy losses for solar-flare-accelerated particles, we obtain the time-dependent particle distribution and thus the time behavior of the resulting rays. The important processes for producing rays in solar flares are found to be nonthermal electron bremsstrahlung, decay of neutral mesons, positron annihilation, neutron capture, and decay of excited nuclei. The results are applied to several known solar flares. For a large flare such as the class 3+ on February 23, 1956, continuous rays with energies up to 100 MeV from electron bremsstrahlung and neutral meson decays are observable at the orbit of the Earth by existing -ray detectors. Line rays from positron annihilation (0.51 MeV), neutron capture (2.23 MeV), and deexcitation of excited nuclei O16 (6.14 and 7.12 MeV) and C12 (4.43 MeV) are particularly strong and well above the continuous -ray background due to electron bremsstrahlung. These lines can be detected at the Earth.NASA-NRC Resident Research Associate.  相似文献   

15.
In this paper we discuss theoretical expressions, determining the difference of Doppler shifts of various coherent radiowave frequencies emitted by a radiator moving in the ionosphere or interplanetary medium. The rotating Doppler effect (Faraday effect) caused by the Doppler shifts ±H of the ordinary and extraordinary waves is also considered. In a three-dimensional inhomogeneous ionosphere, stationary in time (N/t = 0), is determined in the general case, by an equation with three variables. The equation for proper depends only on the local value of the electron concentration N c around the radiator and on integral values, determining, by means of additional calculations, the angle of refraction or its components, the horizontal gradients of electron concentration N/x and N/y, and in some cases, the integral electron concentration 0 zcN dz. We describe the analysis of the measurements, made with the satellites Cosmos I, II and partially XI, assuming that N/t = N/y = 0, with a two variables equation. The expected errors are considered. The results coincide well for different points (Moscow, The Crimea, Sverdlovsk) and thus agree with the measurements of H and with height-frequency ionospheric characteristics. The curve giving electron concentration versus height N (z) in the outer ionosphere (above the maximum of F2), shows a new maximum higher than the main maximum of the ionosphere N MF2 at 120–140 km. At this maximum the value of N (z) is (0.9–0.95) N MF2. The new data on the large-scale horizontal inhomogeneities of the ionosphere, exceed the previous ones by about a factor 10. By means of the irregular variations of the spectrum W() of the inhomogenous formation is determined. Three unknown constant maxima with values 16 to 18 km, 28 to 32 km and 100 to 120 km are found. The spectrum W () mainly characterizes the local properties of the ionosphere along the orbit of the satellite.  相似文献   

16.
There is a warm tenuous partially ionized cloud (T104 K,n(HI)0.1 cm–3,n(Hii 0.22–0.44 cm–3) surrounding the solar system which regulates the environment of the solar system, determines the structure of the heliopause region, and feeds neutral interstellar gas into the inner solar system. The velocity (V–20 km s–1 froml335°,b0° in the local standard of rest) and enhanced Caii and Feii abundances of this cloud suggest an origin as evaporated gas from cloud surfaces in the Scorpius-Centaurus Association. Although the soft X-ray emission attributed to the Local Bubble is enigmatic, optical and ultraviolet data are consistent with bubble formation caused by star formation epochs in the Scorpius-Centaurus Association as regulated by the nearby spiral arm configuration. The cloud surrounding the solar system (the local fluff) appears to be the leading region of an expanding interstellar structure (the squall line) which contains a magnetic field causing polarization of the light of nearby stars, and also absorption features in nearby upwind stars. The velocity vectors of the solar system and local fluff are perpendicular in the local standard of rest. Combining this information with the low column densities seen towards Sirius in the anti-apex direction, and the assumption that the cloud velocity vector is parallel to the surface normal, suggests that the Sun entered the local fluff within the historical past (less than 10 000 years ago) and is skimming the surface of the cloud. Comparison of magnesium absorption lines towards Sirius and anomalous cosmic-ray data suggest the local fluff is in ionization equilibrium.Reason has moons, but moons not hers, Lie mirror'd on her sea, Confounding her astronomers, But, O! delighting me.Ralph Hodgson  相似文献   

17.
Summary The most striking aspect of the Cep and 53 Per stars is their complexity. Whereas in Cepheid-type variables, a dominant mechanism excites a dominant mode (or two at most) of a dominant kind of pulsation, in Cep stars, a number of mechanisms and processes are at work. Still there is hope that the mystery will soon be unraveled by the careful application of a combination of observational and theoretical techniques. These same techniques will provide a better understanding of B stars in general: their interior and atmosphere, mass loss and coronal heating.  相似文献   

18.
Høg  E.  Pagel  B.E.J.  Portinari  L.  Thejll  P.A.  Macdonald  J.  Girardi  L. 《Space Science Reviews》1998,84(1-2):115-126
The primordial helium abundance YP is important for cosmology and the ratio Y/Z of the changes relative to primordial abundances constrains models of stellar evolution. While the most accurate estimates of YP come from emission lines in extragalactic H II regions, they involve an extrapolation to zero metallicity which itself is closely tied up with the slope Y/Z. Recently certain systematic effects have come to light in this exercise which make it useful to have an independent estimate of Y/Z from fine structure in the main sequence of nearby stars. We derive such an estimate from Hipparcos data for stars with Z Z and find values between 2 and 3, which are consistent with stellar models, but still have a large uncertainty.  相似文献   

19.
Methods are discussed for establishing the optical identification of X ray sources in the medium and deep X-ray surveys of the Einstein Observatory. Of the 63 X-ray sources with a statistical significance of 5 in the medium survey (Maccacaro et al. 1981), optical identification work is summarized for 51, of which identifications have been made with 30 active galactic nuclei. The optical properties of some of these X-ray selected objects are briefly discussed.The Einstein deep survey of Pavo (Griffiths et al. 1981) is used to illustrate the problems and methods used for securing optical identifications for X-ray sources in the deep survey fields. Identifications have been made with 4 QSOs at the bright end of the optical candidate distribution (together with 3 G stars) and it is shown that a further 7 fainter objects are also likely to be QSOs.  相似文献   

20.
Bursts of massive star formation can temporarily dominate the luminosity of galaxies spanning a wide range of morphological types. This review is concerned primarily with such events in the central 1 kpc region of spiral galaxies which result from bar driven inflows of gas triggered by interactions or mergers. Most of the stellar radiant luminosity of such bursts is absorbed by dust and re-emitted in the far-infrared and is accompanied by radio and X-ray emission from supernova remnants which can also act collectively to drive galaxy scale outflows. Both evolutionary stellar models and estimates of the gas depletion times are consistent with typical burst durations of 107–8 yr. Spatially-resolved studies of nearby starburst galaxies reveal that the activity is distributed over many individual star forming complexes within rings and other structures organized by interactions between bars and the disc over a range of scales. More distant and extreme examples associated with mergers of massive spirals have luminosities > 1013 L and molecular gas masses > 1010 M implying star formation rates > 1000 M yr–1 which can only be sustained for 107 yr. In the most luminous merging systems, however, the relative importance of starburst and AGN activity and their possible evolutionary connection is still a hotly debated issue. Also controversial are suggestions that starbursts in addition to a black hole are required to account for the properties of AGNs or that starbursts alone may be sufficient under certain conditions. In a wider context, starbursts must clearly have played an important role in galaxy formation and evolution at earlier times. Recent detections of high redshift galaxies show that star formation was underway at z 4 but do not support a continuing increase of the strong evolution in the co-moving star formation density seen out to z l. Primeval starburst pre-cursors of spheroidal systems also remain elusive. The most distant candidates are radio galaxies and quasars at z = 4–5 and a possible population of objects responsible for an isotropic sub-mm wave background tentatively claimed to have been detected by the COBE satellite.  相似文献   

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

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