首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 22 毫秒
1.
R. Helled  P. Bodenheimer 《Icarus》2011,211(2):939-947
Giant protoplanets formed by gravitational instability in the outer regions of circumstellar disks go through an early phase of quasi-static contraction during which radii are large (∼1 AU) and internal temperatures are low (<2000 K). The main source of opacity in these objects is dust grains. We investigate two problems involving the effect of opacity on the evolution of isolated, non-accreting planets of 3, 5, and 7 MJ. First, we pick three different overall metallicities for the planet and simply scale the opacity accordingly. We show that higher metallicity results in slower contraction as a result of higher opacity. It is found that the pre-collapse time scale is proportional to the metallicity. In this scenario, survival of giant planets formed by gravitational instability is predicted to be more likely around low-metallicity stars, since they evolve to the point of collapse to small size on shorter time scales. But metal-rich planets, as a result of longer contraction times, have the best opportunity to capture planetesimals and form heavy-element cores. Second, we investigate the effects of opacity reduction as a result of grain growth and settling, for the same three planetary masses and for three different values of overall metallicity. When these processes are included, the pre-collapse time scale is found to be of order 1000 years for the three masses, significantly shorter than the time scale calculated without these effects. In this case the time scale is found to be relatively insensitive to planetary mass and composition. However, the effects of planetary rotation and accretion of gas and dust, which could increase the timescale, are not included in the calculation. The short time scale we find would preclude metal enrichment by planetesimal capture, as well as heavy-element core formation, over a large range of planetary masses and metallicities.  相似文献   

2.
The giant planets of our solar system possess envelopes consisting mainly of hydrogen and helium but are also significantly enriched in heavier elements relatively to our Sun. In order to better constrain how these heavy elements have been delivered, we quantify the amount accreted during the so-called “late heavy bombardment”, at a time when planets were fully formed and planetesimals could not sink deep into the planets. On the basis of the “Nice model”, we obtain accreted masses (in terrestrial units) equal to for Jupiter, and for Saturn. For the two other giant planets, the results are found to depend mostly on whether they switched position during the instability phase. For Uranus, the accreted mass is with an inversion and without an inversion. Neptune accretes in models in which it is initially closer to the Sun than Uranus, and otherwise. With well-mixed envelopes, this corresponds to an increase in the enrichment over the solar value of 0.033±0.001 and 0.074±0.007 for Jupiter and Saturn, respectively. For the two other planets, we find the enrichments to be 2.1±1.4 (w/ inversion) or 1.2±0.7 (w/o inversion) for Uranus, and 2.0±1.2 (w/ inversion) or 2.7±1.6 (w/o inversion) for Neptune. This is clearly insufficient to explain the inferred enrichments of ∼4 for Jupiter, ∼7 for Saturn and ∼45 for Uranus and Neptune.  相似文献   

3.
Ravit Helled  Gerald Schubert 《Icarus》2008,198(1):156-162
Sedimentation rates of silicate grains in gas giant protoplanets formed by disk instability are calculated for protoplanetary masses between 1 MSaturn to 10 MJupiter. Giant protoplanets with masses of 5 MJupiter or larger are found to be too hot for grain sedimentation to form a silicate core. Smaller protoplanets are cold enough to allow grain settling and core formation. Grain sedimentation and core formation occur in the low mass protoplanets because of their slow contraction rate and low internal temperature. It is predicted that massive giant planets will not have cores, while smaller planets will have small rocky cores whose masses depend on the planetary mass, the amount of solids within the body, and the disk environment. The protoplanets are found to be too hot to allow the existence of icy grains, and therefore the cores are predicted not to contain any ices. It is suggested that the atmospheres of low mass giant planets are depleted in refractory elements compared with the atmospheres of more massive planets. These predictions provide a test of the disk instability model of gas giant planet formation. The core masses of Jupiter and Saturn were found to be ∼0.25 M and ∼0.5 M, respectively. The core masses of Jupiter and Saturn can be substantially larger if planetesimal accretion is included. The final core mass will depend on planetesimal size, the time at which planetesimals are formed, and the size distribution of the material added to the protoplanet. Jupiter's core mass can vary from 2 to 12 M. Saturn's core mass is found to be ∼8 M.  相似文献   

4.

Context

Current planet search programs are detecting extrasolar planets at a rate of 60 planets per year. These planets show more diverse properties than was expected.

Aims

We try to get an overview of possible gas giant (proto-) planets for a full range of orbital periods and stellar masses. This allows the prediction of the full range of possible planetary properties which might be discovered in the near future.

Methods

We calculate the purely hydrostatic structure of the envelopes of proto-planets that are embedded in protoplanetary disks for all conceivable locations: combinations of different planetesimal accretion rates, host star masses, and orbital separations. At each location all hydrostatic equilibrium solutions to the planetary structure equations are determined by variation of core mass and pressure over many orders of magnitude. For each location we analyze the distribution of planetary masses.

Results

We get a wide spectrum of core-envelope structures. However, practically all calculated proto-planets are in the planetary mass range. Furthermore, the planet masses show a characteristic bimodal, sometimes trimodal, distribution. For the first time, we identify three physical processes that are responsible for the three characteristic planet masses: self-gravity in the Hill sphere, compact objects, and a region of very low adiabatic pressure gradient in the hydrogen equation of state. Using these processes, we can explain the dependence of the characteristic masses on the planet’s location: orbital period, host star mass, and planetesimal accretion rate (luminosity). The characteristic mass caused by the self-gravity effect at close proximity to the host star is typically one Neptune mass, thus producing the so-called hot Neptunes.

Conclusions

Our results suggest that hot Jupiters with orbital period less than 64 days (the exact location of the boundary depends on stellar type and accretion rate) have quite distinct properties which we expect to be reflected in a different mass distribution of these planets when compared to the “normal” planetary population. We use our theoretical survey to produce an upper mass limit for embedded planets: the maximum embedded equilibrium mass (MEEM). This naturally explains the lack of high mass planets between 3 and 64 days orbital period.  相似文献   

5.
6.
Junko Kominami  Shigeru Ida 《Icarus》2004,167(2):231-243
We have performed N-body simulations on final accretion stage of terrestrial planets, including the eccentricity and inclination damping effect due to tidal interaction with a gas disk. We investigated the dependence on a depletion time scale of the disk, and the effect of secular perturbations by Jupiter and Saturn. In the final stage, terrestrial planets are formed through coagulation of protoplanets of about the size of Mars. They would collide and grow in a decaying gas disk. Kominami and Ida [Icarus 157 (2002) 43-56] showed that it is plausible that Earth-sized, low-eccentricity planets are formed in a mostly depleted gas disk. In this paper, we investigate the formation of planets in a decaying gas disk with various depletion time scales, assuming disk surface density of gas component decays exponentially with time scale of τgas. Fifteen protoplanets with are initially distributed in the terrestrial planet regions. We found that Earth-sized planets with low eccentricities are formed, independent of initial gas surface density, when the condition (τcross+τgrowth)/2?τgas?τcross is satisfied, where τcross is the time scale for initial protoplanets to start orbit crossing in a gas-free case and τgrowth is the time scale for Earth-sized planets to accrete during the orbit crossing stage. In the cases satisfying the above condition, the final masses and eccentricities of the largest planets are consistent with those of Earth and Venus. However, four or five protoplanets with the initial mass remain. In the final stage of terrestrial planetary formation, it is likely that Jupiter and Saturn have already been formed. When Jupiter and Saturn are included, their secular perturbations pump up eccentricities of protoplanets and tend to reduce the number of final planets in the terrestrial planet regions. However, we found that the reduction is not significant. The perturbations also shorten τcross. If the eccentricities of Jupiter and Saturn are comparable to or larger than present values (∼0.05), τcross become too short to satisfy the above condition. As a result, eccentricities of the planets cannot be damped to the observed value of Earth and Venus. Hence, for the formation of terrestrial planets, it is preferable that the secular perturbations from Jupiter and Saturn do not have significant effect upon the evolution. Such situation may be reproduced by Jupiter and Saturn not being fully grown, or their eccentricities being smaller than the present values during the terrestrial planets' formation. However, in such cases, we need some other mechanism to eliminate the problem that numerous Mars-sized planets remain uncollided.  相似文献   

7.
系外类地行星是目前搜寻地外生命的主要目标.随着观测仪器的发展,现在已经能探测到低于10个地球质量的系外行星.该文简要回顾了系外类地行星的形成与演化,介绍了当前研究它们内部结构的模型和方法,以及由此得出的类地行星质量-半径关系.同时,对应不同的行星初始物质成分,讨论了各种可能的大气结构.最后介绍了未来的空间任务在相关方面的工作.  相似文献   

8.
By considering the physical and orbital characteristics of G type stars and their exoplanets, we examine the association between stellar mass and its metallicity that follows a power law. Similar relationship is also obtained in case of single and multiplanetary stellar systems suggesting that, \(\hbox {Sun}^{\prime }\)s present mass is about 1% higher than the estimated value for its metallicity. Further, for all the stellar systems with exoplanets, association between the planetary mass and the stellar metallicity is investigated, that suggests planetary mass is independent of stellar metallicity. Interestingly, in case of multiplanetary systems, planetary mass is linearly dependent on the stellar absolute metallicity, that suggests, metal rich stars produce massive (\(\ge \)1 Jupiter mass) planets compared to metal poor stars. This study also suggests that there is a solar system planetary missing mass of \({\sim }\)0.8 Jupiter mass. It is argued that probably 80% of missing mass is accreted onto the Sun and about 20% of missing mass might have been blown off to the outer solar system (beyond the present Kuiper belt) during early history of solar system formation. We find that, in case of single planetary systems, planetary mass is independent of stellar metallicity with an implication of their non-origin in the host star’s protoplanetary disk and probably are captured from the space. Final investigation of dependency of the orbital distances of planets on the host stars metallicity reveals that inward migration of planets is dominant in case of single planetary systems supporting the result that most of the planets in single planetary systems are captured from the space.  相似文献   

9.
We present results from 44 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary “embryos” (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 44 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 “habitable” planets between 0.9 and 1.1 AU. These planets range from dry worlds to “water worlds” with 100+oceans of water (1 ocean=1.5×1024 g), and vary in mass between 0.23M and 3.85M. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call “super embryos,” 1-2M protoplanets which can serve as the accretion seeds for 2+M planets with large water contents.  相似文献   

10.
Orbital resonances may have played an important role in determining the locations where the planetesimal swarm eventually accreted into full-size planets. Several pairs of planets do indeed have commensurable orbital periods at present, but the case for control of planet formation by resonances is weakened by the fact that many pairs are not commensurable and that those which are do not necessarily exist at the strongest resonances. However, the mass loss and redistribution that occurred in the early solar system evolution can substantially alter the positions of planets and planetary embryos within the swarm. A cascaded resonance structure is hypothesized where planetesimal growth was accelerated at 2:1 interior and 1:2 exterior resonances with an early-formed Jupiter producing runaway growth of planetary embryos. These embryos produce their own resonances which, in turn, lead to additional embryos in a process that successively propagates inward and outward to generate a resonant configuration of embryos. In this manner, the early presence of Jupiter imposed a harmonic structure on the accumulating planetesimal swarm. For an accretion disk with surface density obeying a power law of index ?1.2 the positions of the planetary embryos can be moved into a reasonably good agreement with most of the present planetary positions that is as good as that given by the Titius-Bode law.  相似文献   

11.
As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct N-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about ∼0.1 Earth mass within 107 years. This brings another difficulty in formation of gas giant planets, which requires cores with ∼10 Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos’ collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the Solar System.  相似文献   

12.
New numerical simulations of the formation and evolution of Jupiter are presented. The formation model assumes that first a solid core of several M accretes from the planetesimals in the protoplanetary disk, and then the core captures a massive gaseous envelope from the protoplanetary disk. Earlier studies of the core accretion-gas capture model [Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y., 1996. Icarus 124, 62-85] demonstrated that it was possible for Jupiter to accrete with a solid core of 10-30 M in a total formation time comparable to the observed lifetime of protoplanetary disks. Recent interior models of Jupiter and Saturn that agree with all observational constraints suggest that Jupiter's core mass is 0-11 M and Saturn's is 9-22 M [Saumon, G., Guillot, T., 2004. Astrophys. J. 609, 1170-1180]. We have computed simulations of the growth of Jupiter using various values for the opacity produced by grains in the protoplanet's atmosphere and for the initial planetesimal surface density, σinit,Z, in the protoplanetary disk. We also explore the implications of halting the solid accretion at selected core mass values during the protoplanet's growth. Halting planetesimal accretion at low core mass simulates the presence of a competing embryo, and decreasing the atmospheric opacity due to grains emulates the settling and coagulation of grains within the protoplanet's atmosphere. We examine the effects of adjusting these parameters to determine whether or not gas runaway can occur for small mass cores on a reasonable timescale. We compute four series of simulations with the latest version of our code, which contains updated equation of state and opacity tables as well as other improvements. Each series consists of a run without a cutoff in planetesimal accretion, plus up to three runs with a cutoff at a particular core mass. The first series of runs is computed with an atmospheric opacity due to grains (hereafter referred to as ‘grain opacity’) that is 2% of the interstellar value and . Cutoff runs are computed for core masses of 10, 5, and 3 M. The second series of Jupiter models is computed with the grain opacity at the full interstellar value and . Cutoff runs are computed for core masses of 10 and 5 M. The third series of runs is computed with the grain opacity at 2% of the interstellar value and . One cutoff run is computed with a core mass of 5 M. The final series consists of one run, without a cutoff, which is computed with a temperature dependent grain opacity (i.e., 2% of the interstellar value for ramping up to the full interstellar value for ) and . Our results demonstrate that reducing grain opacities results in formation times less than half of those for models computed with full interstellar grain opacity values. The reduction of opacity due to grains in the upper portion of the envelope with has the largest effect on the lowering of the formation time. If the accretion of planetesimals is not cut off prior to the accretion of gas, then decreasing the surface density of planetesimals lowers the final core mass of the protoplanet, but increases the formation timescale considerably. Finally, a core mass cutoff results in a reduction of the time needed for a protoplanet to evolve to the stage of runaway gas accretion, provided the cutoff mass is sufficiently large. The overall results indicate that, with reasonable parameters, it is possible that Jupiter formed at 5 AU via the core accretion process in 1 Myr with a core of 10 M or in 5 Myr with a core of 5 M.  相似文献   

13.
Most stars reside in binary/multiple star systems; however, previous models of planet formation have studied growth of bodies orbiting an isolated single star. Disk material has been observed around both components of some young close binary star systems. Additionally, it has been shown that if planets form at the right places within such disks, they can remain dynamically stable for very long times. Herein, we numerically simulate the late stages of terrestrial planet growth in circumbinary disks around ‘close’ binary star systems with stellar separations 0.05 AU?aB?0.4 AU and binary eccentricities 0?eB?0.8. In each simulation, the sum of the masses of the two stars is 1 M, and giant planets are included. The initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet formation within our Solar System by Chambers [Chambers, J.E., 2001. Icarus 152, 205-224], and around each individual component of the α Centauri AB binary star system by Quintana et al. [Quintana, E.V., Lissauer, J.J., Chambers, J.E., Duncan, M.J., 2002. Astrophys. J. 576, 982-996]. Multiple simulations are performed for each binary star system under study, and our results are statistically compared to a set of planet formation simulations in the Sun-Jupiter-Saturn system that begin with essentially the same initial disk of protoplanets. The planetary systems formed around binaries with apastron distances QB≡aB(1+eB)?0.2 AU are very similar to those around single stars, whereas those with larger maximum separations tend to be sparcer, with fewer planets, especially interior to 1 AU. We also provide formulae that can be used to scale results of planetary accretion simulations to various systems with different total stellar mass, disk sizes, and planetesimal masses and densities.  相似文献   

14.
We compute the growth of isolated gaseous giant planets for several values of the density of the protoplanetary disk, several distances from the central star and two values for the (fixed) radii of accreted planetesimals. Calculations were performed in the frame of the core instability mechanism and the solids accretion rate adopted is that corresponding to the oligarchic growth regime. We find that for massive disks and/or for protoplanets far from the star and/or for large planetesimals, the planetary growth occurs smoothly. However, notably, there are some cases for which we find an envelope instability in which the planet exchanges gas with the surrounding protoplanetary nebula. The timescale of this instability shows that it is associated with the process of planetesimals accretion. The presence of this instability makes it more difficult the formation of gaseous giant planets.  相似文献   

15.
The accumulation of giant planets involves processes typical for terrestrial planet formation as well as gasdynamic processes that were previously known only in stars. The condensible element cores of the gas-giants grow by solid body accretion while envelope formation is governed by stellar-like equilibria and the dynamic departures thereof. Two hypotheses for forming Uranus/Neptune-type planets — at sufficiently large heliocentric distances while allowing accretion of massive gaseous envelopes, i.e. Jupiter-type planets at intermediate distances — have been worked out in detailed numerical calculations: (1) Hydrostatic gas-accretion models with time-dependent solid body accretion-rates show a slow-down of core-accretion at the appropriate masses of Uranus and Neptune. As a consequence, gas-accretion also stagnates and a window is opened for removing the solar nebula during a time of roughly constant envelope mass. (2) Gasdynamic calculations of envelope accretion for constant planetesimal accretion-rates show a dynamic transition to new envelope equilibria at the so called critical mass. For a wide range of solar nebula conditions the new envelopes have respective masses similar to those of Uranus and Neptune and are more tightly bound to the cores. The transitions occur under lower density conditions typical for the outer parts of the solar nebula, whereas for higher densities, i.e. closer to the Sun, gasdynamic envelope accretion sets in and is able to proceed to Jupiter-masses.  相似文献   

16.
17.
According to current observational data, planets of many exoplanetary systems have resonant motion. The formation of resonance configurations is studied within a unified model of planetary migration. Planets in the observed systems 24 Sex, HD 37124, HD 73526, HD 82943, HD 128311, HD 160691, Kepler 9, NN Ser, which are moving in the 2: 1 resonance, could have been captured into this resonance due to both the Type I and II migration with a wide range of parameters. The migration conditions are defined for the formation of HD 45364 and HD 200964 that are in the 3: 2 and 4: 3 first-order resonances, correspondingly. The results obtained for HD 200964 show that planets can be captured in the first-order resonances, when the outer-to-inner orbital period ratios for the planets are less than 3: 2, only if Type I migration rates are large, and the mass of at least one planet is substantially less than the modern masses of the observed giant planets. The formation of the HD 102272, HD 108874, HD 181433 and HD 202206 systems with planets in high-order resonances is considered. The capture into these resonances can be realized with very slow Type II migration. Possible bounds for migration parameters are considered. In particular, it has been found that the capture of HD 108874 into the 4: 1 resonance is possible only if the angle between the plane of planetary orbits and the plane of sky is appreciably less than 90°, i.e., the planetary masses are a few times larger than the minimum values. The capture of HD 202206 into the 5: 1 resonance is possible at low migration rates; however, another mechanism is required to explain the high observed eccentricity of the inner planet (for example, strong gravitational interaction between the planets). Resonant configurations can be disrupted due to the interaction between planets and remaining fragments of the planetesimal disk as, for example, may occur in the three-planet system 47 UMa. The specific orbital features observed for this system are explained.  相似文献   

18.
The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.  相似文献   

19.
When protoplanets growing by accretion of planetesimals have atmospheres, small planetesimals approaching the protoplanets lose their energy by gas drag from the atmospheres, which leads them to be captured within the Hill sphere of the protoplanets. As a result, growth rates of the protoplanets are enhanced. In order to study the effect of an atmosphere on planetary growth rates, we performed numerical integration of orbits of planetesimals for a wide range of orbital elements and obtained the effective accretion rates of planetesimals onto planets that have atmospheres. Numerical results are obtained as a function of planetesimals’ eccentricity, inclination, planet’s radius, and non-dimensional gas-drag parameters which can be expressed by several physical quantities such as the radius of planetesimals and the mass of the protoplanet. Assuming that the radial distribution of the gas density near the surface can be approximated by a power-law, we performed analytic calculation for the loss of planetesimals’ kinetic energy due to gas drag, and confirmed agreement with numerical results. We confirmed that the above approximation of the power-law density distribution is reasonable for accretion rate of protoplanets with 1-10 Earth masses, unless the size of planetesimals is too small. We also calculated the accretion rates of planetesimals averaged over a Rayleigh distribution of eccentricities and inclinations, and derived a semi-analytical formula of accretion rates, which reproduces the numerical results very well. Using the obtained expression of the accretion rate, we examined the growth of protoplanets in nebular gas. We found that the effect of atmospheric gas drag can enhance the growth rate significantly, depending on the size of planetesimals.  相似文献   

20.
Planets orbiting a planetesimal circumstellar disc can migrate inward from their initial positions because of dynamical friction between planets and planetesimals. The migration rate depends on the disc mass and on its time evolution. Planets that are embedded in long-lived planetesimal discs, having total mass of 10−4– 0.01 M , can migrate inward a large distance and can survive only if the inner disc is truncated or as a result of tidal interaction with the star. In this case the semimajor axis, a , of the planetary orbit is less than 0.1 au. Orbits with larger a are obtained for smaller values of the disc mass or for a rapid evolution (depletion) of the disc. This model may explain not only several of the orbital features of the giant planets that have been discovered in recent years orbiting nearby stars, but also the metallicity enhancement found in several stars associated with short-period planets.  相似文献   

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

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