首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
We analyze the observations of the transient X-ray pulsar 4U 0115+63 with the RXTE and INTEGRAL observatories in a wide X-ray (3–100 keV) energy band during its intense outbursts in 1999 and 2004. The energy of the fundamental harmonic of the cyclotron resonance absorption line near the maximum of the X-ray flux from the source (luminosity range 5 × 1037–2 × 1038 erg s?1) is ~11 keV. When the pulsar luminosity falls below ~5 × 1037 erg s?1, the energy of the fundamental harmonic is displaced sharply toward the high energies, up to ~16 keV. Under the assumption of a dipole magnetic field configuration, this change in cyclotron harmonic energy corresponds to a decrease in the height of the emitting region by ~2 km, while other spectral parameters, in particular, the cutoff energy, remain essentially constant. At a luminosity ~7 × 1037 erg s?1, four almost equidistant cyclotron line harmonics are recorded in the spectrum. This suggests that either the region where the emission originates is compact or the emergent spectrum from different (in height) segments of the accretion column is uniform. We have found significant pulse profile variations with energy, luminosity, and time. In particular, we show that the profile variations from pulse to pulse are not reduced to a simple modulation of the accretion rate specified by external conditions.  相似文献   

2.
We investigate the relation between the optical (g-band) and X-ray (0.5–10 keV) luminosities of accreting nonmagnetic white dwarfs. According to the present-day counts of the populations of star systems in our Galaxy, these systems have the highest space density among the close binary systems with white dwarfs. We show that the dependence of the optical luminosity of accreting white dwarfs on their X-ray luminosity forms a fairly narrow one-parameter curve. The typical half-width of this curve does not exceed 0.2–0.3 dex in optical and X-ray luminosities, which is essentially consistent with the amplitude of the aperiodic flux variability for these objects. At X-ray luminosities L x ~ 1032 erg s?1 or lower, the optical g-band luminosity of the accretion flow is shown to be related to its X-ray luminosity by a factor ~2–3. At even lower X-ray luminosities (L x ? 1030 erg s?1), the contribution from the photosphere of the white dwarf begins to dominate in the optical spectrum of the binary system and its optical brightness does not drop below M g ~ 13–14. Using the latter fact, we show that in current and planned X-ray sky surveys, the family of accreting nonmagnetic white dwarfs can be completely identified to the distance determined by the sensitivity of an optical sky survey in this region. For the Sloan Digital Sky Survey (SDSS) with a limiting sensitivity m g ~ 22.5, this distance is ~400–600 pc.  相似文献   

3.
The peculiarities of non-Hubble bulk motions of galaxies are studied by analyzing a sample of 1271 thin edge-on spirals with distances determined using a multiparametric Tully-Fisher relation that includes the amplitude of the galaxy rotation, the blue and red diameters, surface brightness, and morphological type. In the purely dipole approximation, the bulk motion of galaxies relative to the cosmic microwave background frame can be described by the velocity of 336±96 km s?1 in the direction l=321°, b=?1° within radius R max =10000 km s?1. An analysis of more complex velocity field models shows that the anisotropy of the Hubble expansion described by the quadrupole term is equal to ~5% on scale lengths R max=6000–10000 km s?1. The amplitude within the Local Supercluster (R max=3000 km s?1) is as high as ~20%. The inclusion of the octupole component reduces the dipole amplitude to 134±111 km s?1 on scale lengths of ~8000 km s?1. The most remarkable feature of the galaxy velocity field within R max=8000 km s?1 is the zone of minimum centered on l=80°, b=0° (the constellation of Cygnus) whose amplitude reaches 18% of the mean Hubble velocity.  相似文献   

4.
We consider the problem of strange-star (SS) radiation. The bare quark SS surface and electrons on the stellar surface generate an electric field that is strong enough for electron-positron pairs to be produced from a vacuum at a nonzero temperature. The luminosity in pairs is assumed to be within ?1049 erg s?1 from a surface with a characteristic radius of 10 km. We consider the energy transfer from pairs to photons by taking into account the well-studied reactions between e, e +, γ and obtain a change in the photon spectrum with luminosity. Our analysis is restricted to the spherically symmetric case. The magnetic field is disregarded. To solve the problem, we developed a new numerical method of integrating the Boltzmann kinetic equations for pairs and photons. This method is used to calculate the problem up to a luminosity of 1042 erg s?1 This region is difficult to investigate when the optical path for pairs or photons is considerably larger than unity but the two optical depths are not simultaneously much larger than unity (when hydrodynamics with heat conduction is applicable). It turns out that the mean photon energy is approximately equal to $\bar \in _\gamma \approx m_e c^2$ (the annihilation line for pairs) at a modest luminosity, L?1×1037 erg s?1, and decreases to ≈210 keV at L?1038 erg s?1. Hydrodynamic estimates point to an increase in the mean energy $\bar \in _\gamma$ to 1 MeV as the luminosity further increases to L?1049 erg s?1. Our calculations may prove to be useful in interpreting soft gamma repeaters (SGRs) and are of methodological interest.  相似文献   

5.
The apparatus of correlation gamma function (Γ*(r)) is used to analyze volume-limited samples from the DR4 Main Galaxy Sample of the SDSS survey with the aim of determining the characteristic scales of galaxy clustering. Up to 20h ?1 Mpc (H 0 = 65 km s?1 Mpc?1), the distribution of galaxies is described by a power-law density—distance dependence, Γ*(r) ∝ r , with an index γ ≈ 1.0. A change in the state of clustering (a significant deviation from the power law) was found on a scale of (20–25) h ?1 Mpc. The distribution of SDSS galaxies becomes homogeneous (γ ~ 0) from a scale of ~60h ?1 Mpc. The dependence of γ on the luminosity of galaxies in volume-limited samples was obtained. The power-law index γ increases with decreasing absolute magnitude of sample galaxies M abs. At M abs ~ ?21.4, which corresponds to the characteristic value M r * of the SDSS luminosity function, this dependence exhibits a break followed by a more rapid increase in γ.  相似文献   

6.
Here we present the results of panoramic and long-slit observations of eight ULX nebular counterparts performed with the 6m SAO telescope. In two ULX nebulae (ULXNe) we detected for the first time signatures of high excitation ([O III]λ5007 / Hβ > 5). Two of the ULXs were identified with young (T ~ 5–10 Myr) massive star clusters. Four of the eight ULXNe show bright high-excitation lines. This requires existence of luminous (~ 1038 ÷ 1040 erg s?1) UV/EUV sources coinciding with the X-ray sources. The other 4 ULXNe require shock excitation of the gas with shock velocities of 20–100 km s ?1. However, all the studied ULXNe spectra show signatures of shock excitation, but even those ULXNe where the shocks are prevailing show presence of a hard ionizing source with a luminosity of at least ~ 1038 erg s?1. Most likely shock waves, X-ray and EUV ionization act simultaneously in all the ULXNe, but they may be roughly separated in two groups: shock-dominated and photoionization-dominated ULXNe. The ULXs have to produce strong winds and/or jets (~ 1039 erg s?1) for powering their nebulae. Both the wind/jet activity and the existence of a bright UV source are consistent with the suggestion that ULXs are high-mass X-ray binaries with supercritical accretion disks of the SS433 type.  相似文献   

7.
We present the results of measurements of the total X-ray flux from the Andromeda galaxy (M31) in the 3-100 keV band based on data from the RXTE/PCA, INTEGRAL/ISGRI, and SWIFT/BAT space experiments. We show that the total emission from the galaxy has a multicomponent spectrum whose main characteristics are specified by binaries emitting in the optically thick and optically thin regimes. The galaxy’s luminosity at energies 20–100 keV gives about 6% of its total luminosity in the 3–100 keV band. The emissivity of the stellar population in M31 is L 2–20 keV ~ 1.1 × 1029 erg s?1 M ?1 in the 2–20 keV band and L 20–100 keV ~ 8 × 1027 erg s?1 M ?1 in the 20–100 keV band. Since low-mass X-ray binaries at high luminosities pass into a soft state with a small fraction of hard X-ray emission, the detection of individual hard X-ray sources in M31 requires a sensitivity that is tens of times better (up to 10?13 erg s?1 cm?2) than is needed to detect the total hard X-ray emission from the entire galaxy. Allowance for the contribution from the hard spectral component of the galaxy changes the galaxy’s effective Compton temperature approximately by a factor of 2, from ~1.1 to ~2.1 keV.  相似文献   

8.
The Pulkovo Compilation of Radial Velocities (PCRV) has been made to study the stellar kinematics in the local spiral arm. The PCRV contains weighted mean absolute radial velocities for 35 495 Hipparcos stars of various spectral types and luminosity classes over the entire celestial sphere mainly within 500 pc of the Sun. The median accuracy of the radial velocities obtained is 0.7 km s?1. Results from 203 publications were used in the catalogue. Four of them were used to improve the radial velocities of standard stars from the IAU list. The radial velocities of 155 standard stars turned out to be constant within 0.3 km s?1. These stars were used to analyze 47 768 mean radial velocities for 37 200 stars from 12 major publications (~80% of all the data used). Zero-point discrepancies and systematic dependences on radial velocity, B-V color index, right ascension, and declination were found in radial velocity differences of the form “publication minus IAU list of standards.” These discrepancies and dependences were approximated and taken into account when calculating the weighted mean radial velocities. 1128 stars whose independent radial-velocity determinations were available at least in three of these publications and agreed within 3 km s?1 were chosen as the work list of secondary standards. Radial-velocity differences of the form “publication minus list of secondary standards” were used by analogy to correct the zero points and systematic dependences in the radial velocities from 33 more publications (~ 13% of the data used). In addition, the radial velocities from 154 minor publications (~7% of the data used) pertaining to well-known instruments were used without any corrections.  相似文献   

9.
We have computed a spherically symmetric model for the interaction of matter ejected during the outburst of a classical nova with the stellar wind from its optical component. This model is used to describe the intense X-ray outburst (the peak 3–20 keV flux was ~2 Crab) of the binary system CI Camelopardalis in 1998. According to our model, the stellar wind from the optical component heated by a strong shock wave produced when matter is ejected from the white dwarf as the result of a thermonuclear explosion on its surface is the emission source in the standard X-ray band. Comparison of the calculated and observed time dependences of the mean radiation temperature and luminosity of the binary system during its outburst has yielded very important characteristics of the explosion. We have been able to measure the velocity of the ejected matter immediately after the onset of the explosion for the first time: it follows from our model that the ejected matter had a velocity of ~2700 km s?1 even on 0.1–0.5 day after the outburst onset and it flew with such a velocity for the first 1–1.5 day under an external force, possibly, the radiation pressure from the white dwarf. Subsequently, the matter probably became transparent and began to decelerate. The time dependence of the mean radiation temperature at late expansion phases has allowed us to estimate the mass of the ejected matter, ~10?7–10?6 M . The mass loss rate in the stellar wind required to explain the observed peak luminosity of the binary system during its outburst has been estimated to be \(\dot M\) ~ (1 ? 2) × 10?6 M yr?1.  相似文献   

10.
We carried out the first 21-cm line observations of an extended region around the Wolf-Rayet star WR 102 and the associated nebula G2.4+1.4 with the RATAN-600 radio telescope. An irregular H I shell was identified. Its maximum expansion velocity reaches ~50 km s?1, and its outer diameter (at a distance of 3 kpc) is 56 pc. The mechanical luminosity of the stellar wind required to produce the observed shell is estimated to be ~0.8×1038 erg s?1; the age of the shell is ~3.4×105 yr. We compare the inferred parameters of the H I shell with the structure and kinematics of the ionized nebula and with the dust distribution in the region.  相似文献   

11.
Kormendy’s relation (μe-logR e) is used to investigate the structure of the compact Ursa Major supercluster of galaxies (11h30m+55°, cz=18 000 km s?1); this relation allows the distances of early-type galaxies to be estimated. The relative distances of 13 clusters in the supercluster and their peculiar velocities are determined with a mean statistical accuracy of 6%. In general, the supercluster obeys the Hubble relation between radial velocity and distance. However, there is reason to suggest that the supercluster consists of two subsystems with mean radial velocities of 16 200 and 19 700 km s?1. For a velocity dispersion in the subsystems of ~1100 km s?1, the fact that each of them is gravitationally bound is not ruled out.  相似文献   

12.
We obtained U BV R photometric and spectroscopic observations during the outburst of V838 Mon. Before its outburst, the B brightness of the star had been stable ( $\tilde15.^m 85$ ) for 45 years. This was a blue star with the color index $(B - V)_0 = - 0\mathop .\limits^m 03 \pm 0\mathop .\limits^m 1$ and may have been a cataclysmic variable. In the middle of March 2002, the outburst amplitude reached $8\mathop .\limits^m 1$ in B. The star has the counterpart V 1006/7 in M 31 in luminosity at maximum and in spectrum. The unusual spectrum at the premaximum stage originated in the expanding photosphere of a cool K-type giant. The expansion velocity of the photosphere was 150 km s?1; the maximum velocity in the expanding stellar envelope reached 500 km s?1. The absorption components of neutral metal lines were enhanced by a factor of 3 or 4 compared to a normal K-type star. No overabundance of s-process elements was found. One day before the brightness peak, an intense Hα emission line with broad wings, FWZI=3100 km s?1, and numerous lines of ionized metals appeared in V838 Mon, which is characteristic of normal classical novae. We show light, color, and spectral variations of the object.  相似文献   

13.
The ultraviolet spectra of the star RU Lup obtained with the Hubble Space Telescope are analyzed. Emission lines are identified. The presence of absorption components with a nearly zero residual intensity in the Mg II resonance doublet lines is indicative of mass outflow with a velocity V ?300 km s?1. These lines also exhibit a broad (?1400 km s?1 at the base) component originating in the star itself. The profiles of the (optically thin) Si II] and Si III]1892 Å lines for the first time unequivocally prove that these lines originate in an accretion shock wave rather than in the chromosphere, with the gas infall velocity being V 0?400 km s?1. The intensity ratio of the C IV 1550 Å and Si IV 1400 Å resonance doublet components was found to be close to unity, suggesting a high accreted-gas density, logN 0>12.5. Molecular H2 Lyman lines formed in the stellar wind were detected. The H I Lα luminosity of RU Lup was found from their intensities to exceed 10% of L bol. Radiation pressure in the Lα line on atomic hydrogen may play a significant role in the initial acceleration of stellar-wind matter, but the effect of Lα emission on the dynamics of molecular gas is negligible.  相似文献   

14.
We have studied the characteristics of coronal mass ejections (CMEs) associated with Deca-Hectometric (DH) type II radio bursts (1–14 MHz) in the interplanetary medium during the year 1997–2005. The DH CMEs are divided into two parts: (i) DH CMEs (All) and (ii) DH CMEs (Limb). We found that 65% (177/273) of all events have the speed >900 km?s?1 and the remaining 35% (96/273) events have the speed below 900 km?s?1. The average speed of all and limb DH CMEs are 1230 and 1288 km?s?1, respectively, which is nearly three times the average speed of general population of CMEs (473 km?s?1). The average widths of all and limb DH CMEs are 105° and 106°, respectively, which is twice the average width (52°) of the general population of CMEs. We found a better correlation between the speed and width of limb DH CMEs (R=?0.61) than all DH CMEs (R=?0.53). Only 28% (177/637) of fast >900 km?s?1 general population of CMEs are reported with DH type II bursts counterpart. The above results gives that the relation between the CME properties is better for limb events.  相似文献   

15.
Based on the stellar proper motions of the TGAS (Gaia DR1) catalogue, we have analyzed the velocity field of main-sequence stars and red giants from the TGAS catalogue with heliocentric distances up to 1.5 kpc. We have obtained four variants of kinematic parameters corresponding to different methods of calculating the distances from the parallaxes of stars measured with large relative errors. We have established that within the Ogorodnikov–Milne model changing the variant of distances affects significantly only the solar velocity components relative to the chosen centroid of stars, provided that the solution is obtained in narrow ranges of distances (0.1 kpc). The estimates of all the remaining kinematic parameters change little. This allows the Oort coefficients and related Galactic rotation parameters as well as all the remaining Ogorodnikov–Milne model parameters (except for the solar terms) to be reliably estimated irrespective of the parallax measurement accuracy. The main results obtained from main-sequence stars in the range of distances from 0.1 to 1.5 kpc are: A = 16.29 ± 0.06 km s?1 kpc?1, B = ?11.90 ± 0.05 km s?1 kpc?1, C = ?2.99 ± 0.06 km s?1 kpc?1, K = ?4.04 ± 0.16 km s?1 kpc?1, and the Galactic rotation period P = 217.41 ± 0.60 Myr. The analogous results obtained from red giants in the range from 0.2 to 1.6 kpc are: the Oort constants A = 13.32 ± 0.09 km s?1 kpc?1, B = ?12.71 ± 0.06 km s?1 kpc?1, C = ?2.04 ± 0.08 km s?1 kpc?1, K = ?2.72 ± 0.19 km s?1 kpc?1, and the Galactic rotation period P = 236.03 ± 0.98 Myr. The Galactic rotation velocity gradient along the radius vector (the slope of the Galactic rotation curve) is ?4.32 ± 0.08 km s?1 kpc?1 for main-sequence stars and ?0.61 ± 0.11 km s?1 kpc?1 for red giants. This suggests that the Galactic rotation velocity determined from main-sequence stars decreases with increasing distance from the Galactic center faster than it does for red giants.  相似文献   

16.
We have conducted a statistical study 27 coronal mass ejections (CMEs) from January 2007 – June 2008, using the stereoscopic views of STEREO SECCHI A and B combined with SOHO LASCO observations. A flux-rope model, in conjunction with 3D triangulations, has been used to reconstruct the 3D structures and determine the actual speeds of CMEs. The origin and the dynamic evolution of the CMEs are investigated using COR1, COR2 and EUVI images. We have identified four types of solar surface activities associated with CMEs: i) total eruptive prominence (totEP), ii) partially eruptive prominence (PEP), iii) X-ray flare, and iv) X-type magnetic structure (X-line). Among the 27 CMEs, 18.5% (5 of 27) are associated with totEPs, 29.6% (8 of 27) are associated with PEPs, 26% (7 of 27) are flare related, and 26% (7 of 27) are associated with X-line structures, and 43% (3 of 7) are associated with both X-line structures and PEPs. Three (11%) could not be associated with any detectable activity. The mean actual speeds for totEP-CMEs, PEP-CMEs, flare-CMEs, and X-line-CMEs are 404 km?s?1,247 km?s?1,909 km?s?1, and 276 km?s?1, respectively; the average mean values of edge-on and broadside widths for the 27 CMEs are 52 and 85 degrees, respectively. We found that slow CMEs (V≤400 km?s?1) tend to deflect towards and propagate along the streamer belts due to the deflections by the strong polar magnetic fields of corona holes, while some faster CMEs show opposite deflections away from the streamer belts.  相似文献   

17.
Abstract— New model organic microparticles are used to assess the thermal ablation that occurs during aerogel capture at speeds from 1 to 6 km s?1. Commercial polystyrene particles (20 μm diameter) were coated with an ultrathin 20 nm overlayer of an organic conducting polymer, polypyrrole. This overlayer comprises only 0.8% by mass of the projectile but has a very strong Raman signature, hence its survival or destruction is a sensitive measure of the extent of chemical degradation suffered. After aerogel capture, microparticles were located via optical microscopy and their composition was analyzed in situ using Raman microscopy. The ultrathin polypyrrole overlayer survived essentially intact for impacts at ~1 km s?1, but significant surface carbonization was found at 2 km s?1, and major particle mass loss at ≥3 km s?1. Particles impacting at ~6.1 km s?1 (the speed at which cometary dust was collected in the NASA Stardust mission) were reduced to approximately half their original diameter during aerogel capture (i.e., a mass loss of 84%). Thus significant thermal ablation occurs at speeds above a few km s?1. This suggests that during the Stardust mission the thermal history of the terminal dust grains during capture in aerogel may be sufficient to cause significant processing or loss of organic materials. Further, while Raman D and G bands of carbon can be obtained from captured grains, they may well reflect the thermal processing during capture rather than the pre‐impact particle's thermal history.  相似文献   

18.
Microrater frequencies caused by fast (? 3 km s?1) ejecta have been determined using secondary targets in impact experiments. A primary projectile (steel sphere, diam 1.58 mm, mass 1.64 × 10?2 g) was shot in Duran glass with a velocity of 4.1 km s?1 by means of a light gas gun. The angular distribution of the secondary crater number densities shows a primary maximum around 25°, and a secondary maximum at about 60° from the primary target surface. The fraction of mass ejected at velocities of ? 3 km s?1 is only a factor of 7.5 × 10?5 of the primary projectile mass. A conservative calculation shows that the contribution of secondary microcraters (caused by fast ejecta) to primary microcrater densities on lunar rock surfaces (caused by interplanetary particles) is on the statistical average below 1% for any lunar surface orientation. Calculation of the interplanetary dust flux enhancement caused by Moon ejecta turned out to be in good agreement with Lunar Explorer 35in situ measurements.  相似文献   

19.
Two CCD spectra of the star BM Ori were obtained with the echelle spectrograph of the 6-m telescope. In one of the spectra, a large proportion of lines are distorted by emission. The emission component is blueshifted by 50 km s?1, suggesting hot-gas outflow from the atmosphere. The equivalent-width ratio of measured lines in the spectra outside and during eclipse is consistent with the assumption that ~2/3 of the primary star’s area is obscured during eclipse, as follows from light curves. Measured line equivalent widths were used to estimate atmospheric parameters of the secondary star, T eff=7300 K, log g=5.2, and microturbulence ξt=6 km s?1, and to determine its chemical composition. The C, Na, Al, Si, S, Ca, Fe, Ni, and Zn abundances are solar, within the error limits. Li, Sc, Ti, V, Cr, Mn, Co, and Y are overabundant, while Mg, Cu, and Ba are underabundant. In general, the secondary is similar in chemical composition to the star V 1016 Ori. Based on the secondary’s mass determined by solving the radial-velocity curve and on log g estimated spectroscopically from iron ionization equilibrium, we calculated its photospheric radius, R 2 = 0.5R . However, the spectroscopic log g=5.2 disagrees with log g=3.5 calculated from the luminosity and effective temperature and with log g=3.0 calculated from light and radial-velocity curves. If the secondary’s photospheric radius is indeed small; this argues for the hypothesis that the eclipsing body is a dust envelope. The radial velocities measured from the two spectra are systematically higher than those calculated from the radial-velocity curve by +34 and +24 km s?1. It is likely that the secondary’s atmosphere occasionally shrinks.  相似文献   

20.
We consider two samples of OB stars with different distance scales that we have studied previously. The first and second samples consist of massive spectroscopic binaries with photometric distances and distances determined from interstellar calcium lines, respectively. The OB stars are located at heliocentric distances up to 7 kpc. We have identified them with the Gaia DR1 catalogue. Using the proper motions taken from the Gaia DR1 catalogue is shown to reduce the random errors in the Galactic rotation parameters compared to the previously known results. By analyzing the proper motions and parallaxes of 208 OB stars from the Gaia DR1 catalogue with a relative parallax error of less than 200%, we have found the following kinematic parameters: (U, V) = (8.67, 6.63)± (0.88, 0.98) km s?1, Ω0 = 27.35 ± 0.77 km s?1 kpc?1, Ω′0 = ?4.13 ± 0.13 km s?1 kpc?2, and Ω″0 = 0.672 ± 0.070 km s?1 kpc?3, the Oort constants are A = ?16.53 ± 0.52 km s?1 kpc?1 and B = 10.82 ± 0.93 km s?1 kpc?1, and the linear circular rotation velocity of the local standard of rest around the Galactic rotation axis is V 0 = 219 ± 8 km s?1 for the adopted R 0 = 8.0 ± 0.2 kpc. Based on the same stars, we have derived the rotation parameters only from their line-of-sight velocities. By comparing the estimated values of Ω′0, we have found the distance scale factor for the Gaia DR1 catalogue to be close to unity: 0.96. Based on 238 OB stars of the combined sample with photometric distances for the stars of the first sample and distances in the calcium distance scale for the stars of the second sample, line-of-sight velocities, and proper motions from the Gaia DR1 catalogue, we have found the following kinematic parameters: (U, V, W) = (8.19, 9.28, 8.79)± (0.74, 0.92, 0.74) km s?1, Ω0 = 31.53 ± 0.54 km s?1 kpc?1, Ω′0 = ?4.44 ± 0.12 km s?1 kpc?2, and Ω″0 = 0.706 ± 0.100 km s?1 kpc?3; here, A = ?17.77 ± 0.46 km s?1 kpc?1, B = 13.76 ± 0.71 km s?1 kpc?1, and V 0 = 252 ± 8 km s?1.  相似文献   

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

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