首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
The first spectroscopic census of active galactic nuclei (AGNs) associated with late-type galaxies in the Virgo cluster was carried out by observing 213 out of a complete set of 237 galaxies more massive than   M dyn > 108.5 M⊙  . Among them, 77 are classified as AGNs [including 21 transition objects, 47 low-ionization nuclear emission regions (LINERs) and nine Seyferts] and comprise 32 per cent of the late-type galaxies in Virgo. Due to spectroscopic incompleteness, at most 21 AGNs are missed in the survey, so that the fraction would increase up to 41 per cent. Using corollary near-infrared observations that enable us to estimate galaxy dynamical masses, it is found that AGNs are hosted exclusively in massive galaxies, i.e.   M dyn≳ 1010 M⊙  . Their frequency increases steeply with the dynamical mass from zero at   M dyn≈ 109.5 M⊙  to virtually 1 at   M dyn > 1011.5 M⊙  . These frequencies are consistent with those of low-luminosity AGNs found in the general field by the Sloan Digital Sky Survey. Massive galaxies that harbour AGNs commonly show conspicuous r -band star-like nuclear enhancements. Conversely, they often, but not necessarily, contain massive bulges. A few well-known AGNs (e.g. M61, M100, NGC 4535) are found in massive Sc galaxies with little or no bulge. The AGN fraction seems to be only marginally sensitive to galaxy environment. We infer the black hole masses using the known scaling relations of quiescent black holes. No black holes lighter than  ∼106 M⊙  are found active in our sample.  相似文献   

2.
We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 1027 erg s−1 L−1⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 106 M⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L X/SFR  (star formation rate) ratio equals  ∼9 × 1038 (erg s−1)(M⊙ yr−1)−1  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.  相似文献   

3.
Using the spectroscopic sample of the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed approximately 80 per cent of their stars at   z > 1  , while for the less massive ones the value is only approximately 20 per cent. By comparing the BCT with the dark matter merger tree, we find indications that star formation efficiency at   z > 1  had to be approximately a factor of two higher than today (∼10 per cent) in galaxies with present-day stellar mass larger than  2 × 1011 M⊙  , if this early star formation occurred in the main progenitor. Therefore, the λ cold dark matter (LCDM) paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than  1011 M⊙  , efficient star formation seems to have been triggered at   z ∼ 0.2  . We show that there is a characteristic mass  ( M *∼ 1010 M⊙)  for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical models of galaxy formation.  相似文献   

4.
In large spheroidal stellar systems, such as elliptical galaxies, one invariably finds a  106–109 M⊙  supermassive black hole at their centre. In contrast, within dwarf elliptical galaxies one predominantly observes a  105–107 M⊙  nuclear star cluster. To date, few galaxies have been found with both types of nuclei coexisting and even less have had the masses determined for both central components. Here, we identify one dozen galaxies housing nuclear star clusters and supermassive black holes whose masses have been measured. This doubles the known number of such hermaphrodite nuclei – which are expected to be fruitful sources of gravitational radiation. Over the host spheroid (stellar) mass range  108–1011 M⊙  , we find that a galaxy's nucleus-to-spheroid (baryon) mass ratio is not a constant value but decreases from a few per cent to ∼0.3 per cent such that  log[( M BH+ M NC)/ M sph]=−(0.39 ± 0.07) log[ M sph/1010 M⊙]− (2.18 ± 0.07)  . Once dry merging commences and the nuclear star clusters disappear, this ratio is expected to become a constant value.
As a byproduct of our investigation, we have found that the projected flux from resolved nuclear star clusters is well approximated with Sérsic functions having a range of indices from ∼0.5 to ∼3, the latter index describing the Milky Way's nuclear star cluster.  相似文献   

5.
The cluster 3C 129 is classified as a rich cluster. An analysis of the properties of the cluster 3C 129 from ROSAT PSPC and HRI, Einstein IPC, and EXOSAT ME observations is presented. The mean temperature from a joint fit of the ROSAT PSPC and EXOSAT ME data is 5.5(±0.2) keV. The luminosity is 0.6×1044 erg s−1 in 0.2–2.4 keV and 2.7×1044 erg s−1 in 0.2–10 keV. We find a cooling flow with a rate of ∼84 M⊙ yr−1. The central gas density is 6×10−3 cm−3, and the ICM mass is 3.6×1013 M⊙. The total cluster mass is ∼5×1014 M⊙. The X-ray morphology shows an east–west elongation, which is evidence for a recent merger event. The radio source 3C 129.1 is located near the X-ray centre. Another cluster member galaxy (the radio galaxy 3C 129) is a prototype of head-tailed radio galaxies, and is located in the west part of the cluster. The tail points along the gradient of intracluster gas pressure. There are no significant point X-ray sources associated with the AGNs of the two radio galaxies.  相似文献   

6.
We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a nearly constant dry merger rate of  ∼6 × 10−5 Mpc−3 Gyr−1  at   z ≤ 1  and a steep decline at larger z . Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early- and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of   M *, c ∼ 6.3 × 1010 M⊙  by mergers only. As a consequence of this transition, we find that around   M *, c   , the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggests that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z . We however find that only 10–20 per cent of galaxies more massive than   M *, c   experience dry major mergers within their last Gyr at any given redshift   z ≤ 1  .  相似文献   

7.
We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to  7 × 109 M⊙  . Most friends-of-friends haloes in the simulation have a distinct central galaxy, typically 2–5 times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive haloes form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5–1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Typically, satellites in our simulation are 0.1–0.2 mag bluer than in models that assume no gas accretion on to satellites after a halo merger. Since   z = 1  , 27 per cent of central galaxies (above  3 × 1010 M⊙  ) and 22 per cent of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain 'central' objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.  相似文献   

8.
We present the results of Australia Telescope Compact Array (ATCA) H  i line and 20-cm radio continuum observations of the galaxy quartet NGC 6845. The H  i emission extends over all four galaxies but can only be associated clearly with the two spiral galaxies, NGC 6845A and B, which show signs of strong tidal interaction. We derive a total H  i mass of at least  1.8 × 1010 M⊙  , most of which is associated with NGC 6845A, the largest galaxy of the group. We investigate the tidal interaction between NGC 6845A and B by studying the kinematics of distinct H  i components and their relation to the known H  ii regions. No H  i emission is detected from the two lenticular galaxies, NGC 6845C and D. A previously uncatalogued dwarf galaxy, ATCA  J2001−4659  , was detected 4.4 arcmin NE from NGC 6845B and has an H  i mass of  ∼5 × 108 M⊙  . No H  i bridge is visible between the group and its newly detected companion. Extended 20-cm radio continuum emission is detected in NGC 6845A and B as well as in the tidal bridge between the two galaxies. We derive star formation rates of  15–40 M⊙ yr−1  .  相似文献   

9.
Using simple stellar population synthesis, we model the bulge stellar contribution in the optical spectrum of a narrow-line Seyfert 1 galaxy, RE J1034+396. We find that its bulge stellar velocity dispersion is  67.7 ± 8 km s−1  . The supermassive black hole (SMBH) mass is about  (1–4) × 106 M⊙  if it follows the well-known   M BH–σ*  relation found in quiescent galaxies. We also derive the SMBH mass from the Hβ second moment, which is consistent with that from its bulge stellar velocity dispersion. The SMBH mass of (1–4)  × 106 M⊙  implies that the X-ray quasi-periodic oscillation (QPO) of RE J1034+396 can be scaled to a high-frequency QPO at 27–108 Hz found in Galactic black hole binaries with a  10-M⊙  black hole. With the mass distribution in different age stellar populations, we find that the mean specific star formation rate (SSFR) over the past 0.1 Gyr is  0.0163 ± 0.0011  Gyr−1, the stellar mass in the logarithm is  10.155 ± 0.06  in units of solar mass and the current star formation rate is  0.23 ± 0.016 M⊙ yr−1  . For RE J1034+396, there is no relation between the Eddington ratio and the SSFR as suggested by Chen et al., despite a larger scatter in their relation. We also suggest that about 7.0 per cent of the total Hα luminosity and 50 per cent of the total [O  ii ] luminosity come from the star formation process.  相似文献   

10.
The evolution of number density, size and intrinsic colour is determined for a volume-limited sample of visually classified early-type galaxies selected from the Hubble Space Telescope /Advanced Camera for Surveys images of the Great Observatories Origins Deep Survey (GOODS) North and South fields (version 2). The sample comprises 457 galaxies over 320 arcmin2 with stellar masses above  3 × 1010 M⊙  in the redshift range  0.4 < z < 1.2  . Our data allow a simultaneous study of number density, intrinsic colour distribution and size. We find that the most massive systems  (≳3 × 1011 M⊙)  do not show any appreciable change in comoving number density or size in our data. Furthermore, when including the results from 2dF galaxy redshift survey, we find that the number density of massive early-type galaxies is consistent with no evolution between   z = 1.2  and 0, i.e. over an epoch spanning more than half of the current age of the Universe. We find large discrepancies between the predictions of semi-analytic models. Massive galaxies show very homogeneous intrinsic colour distributions, with nearly flat radial colour gradients, but with a significant negative correlation between stellar mass and colour gradient, such that red cores appear predominantly in massive galaxies. The distribution of half-light radii – when compared to   z ∼ 0  and   z > 1  samples – is compatible with the predictions of semi-analytic models relating size evolution to the amount of dissipation during major mergers.  相似文献   

11.
We present a new determination of the local volume-averaged star formation rate from the 1.4-GHz luminosity function of star forming galaxies. Our sample, taken from the   B ≤12  Revised Shapley–Ames catalogue (231 normal spiral galaxies over an effective area of 7.1 sr) has ≃100 per cent complete radio detections and is insensitive to dust obscuration and cirrus contamination. After removal of known active galaxies, the best-fitting Schechter function has a faint-end slope of  −1.27±0.07  in agreement with the local H α luminosity function, characteristic luminosity   L ∗=(2.6±0.7)×1022 W Hz−1  and density   φ ∗=(4.8±1.1)×10−4 Mpc−3.  The inferred local radio luminosity density of  (1.73±0.37±0.03)×1019 W Hz−1 Mpc−3  (Poisson noise, large-scale structure fluctuations) implies a volume-averaged star formation rate ∼2 times larger than the Gallego et al. H α estimate, i.e.   ρ 1.4 GHz=(2.10±0.45±0.04)×10−2 M⊙ yr−1 Mpc−3  for a Salpeter initial mass function from  0.1–125 M⊙  and Hubble constant of 50 km s−1 Mpc−1. We demonstrate that the Balmer decrement is a highly unreliable extinction estimator, and argue that optical–ultraviolet (UV) star formation rates (SFRs) are easily underestimated, particularly at high redshift.  相似文献   

12.
39 galaxies are now known, from follow-up of faint IRAS sources and from submillimetre observations of high-redshift AGN, with far-infrared luminosities >1013 L⊙. 13 of these, which have been found in 60- or 850-μm surveys, form an important unbiased subsample. 12 have been found by comparison of 60-μm surveys with quasar or radio galaxy catalogues, or from infrared surveys with colour selection biased towards AGN, while a further 14 have been found through submillimetre observations of known high-redshift AGN. In this paper I argue, on the basis of detailed modelling of the spectral energy distributions of hyperluminous galaxies with accurate radiative transfer models, and from evidence of high gas mass in several cases, that the bulk of the emission from these galaxies at rest frame wavelengths ≥50 μm is caused by star formation. Even after correction for the effects of lensing, hyperluminous galaxies with emission peaking at rest frame wavelengths ≥50 μm are therefore undergoing star formation at rates >103 M⊙ yr−1 and are strong candidates for being primeval galaxies, in the process of a major episode of star formation.  相似文献   

13.
We present discovery images, together with follow-up imaging and spectroscopy, of two large-separation gravitational lenses found by our survey for wide arcs [the CAmbridge Sloan Survey Of Wide ARcs in the skY (CASSOWARY)]. The survey exploits the multicolour photometry of the Sloan Digital Sky Survey to find multiple blue components around red galaxies. CASSOWARY 2 (or 'the Cheshire Cat') is composed of two massive early-type galaxies at   z = 0.426  and 0.432, respectively, lensing two background sources, the first a star-forming galaxy at   z = 0.97  and the second a high -redshift galaxy  ( z > 1.4)  . There are at least three images of the former source and probably four or more of the latter, arranged in two giant arcs. The mass enclosed within the larger arc of radius ∼11 arcsec is  ∼33 × 1012 M⊙  . CASSOWARY 3 comprises an arc of three bright images of a   z = 0.725  source, lensed by a foreground elliptical at   z = 0.274  . The radius of the arc is ∼4 arcsec and the enclosed mass is  ∼2.5 × 1012 M⊙  . Together with earlier discoveries like the Cosmic Horseshoe and the 8 o'clock Arc, these new systems, with separations intermediate between the arcsecond-separation lenses of typical strong galaxy lensing and arcminute-separation cluster lenses, probe the very high end of the galaxy mass function.  相似文献   

14.
A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 μm has been carried out to map the regions around several 850-μm-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared (IR) luminosities and hence star formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350-μm flux densities for 24 850-μm sources, seven of which are detected at ≥2.5σ within a 10 arcsec search radius of the 850-μm positions. When results from the literature are included the total number of 350-μm flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-IR (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich  ( M dust≃ 9 × 108 M⊙)  , luminous  ( L FIR≃ 2 × 1012 L⊙)  star-forming galaxies with intrinsic dust temperatures of ≃35 K and star formation rates of  ≃400 M⊙ yr−1  . We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350-μm source counts, N 350(>25 mJy) ∼ 200–500 deg−2.  相似文献   

15.
A comparison between published field galaxy stellar mass functions (GSMFs) shows that the cosmic stellar mass density is in the range 4–8 per cent of the baryon density (assuming  Ωb= 0.045  ). There remain significant sources of uncertainty for the dust correction and underlying stellar mass-to-light ratio even assuming a reasonable universal stellar initial mass function. We determine the   z < 0.05  GSMF using the New York University Value-Added Galaxy Catalog sample of 49 968 galaxies derived from the Sloan Digital Sky Survey and various estimates of stellar mass. The GSMF shows clear evidence for a low-mass upturn and is fitted with a double Schechter function that has  α2≃−1.6  . At masses below  ∼108.5 M⊙  , the GSMF may be significantly incomplete because of missing low-surface-brightness galaxies. One interpretation of the stellar mass–metallicity relation is that it is primarily caused by a lower fraction of available baryons converted to stars in low-mass galaxies. Using this principle, we determine a simple relationship between baryonic mass and stellar mass and present an 'implied baryonic mass function'. This function has a faint-end slope,  α2≃−1.9  . Thus, we find evidence that the slope of the low-mass end of the galaxy mass function could plausibly be as steep as the halo mass function. We illustrate the relationship between halo baryonic mass function → galaxy baryonic mass function → GSMF. This demonstrates the requirement for peak galaxy formation efficiency at baryonic masses  ∼1011 M⊙  corresponding to a minimum in feedback effects. The baryonic-infall efficiency may have levelled off at lower masses.  相似文献   

16.
We present a ROSAT and ASCA study of the Einstein source X-9 and its relation to a shock-heated shell-like optical nebula in a tidal arm of the M81 group of interacting galaxies. Our ASCA observation of the source shows a flat and featureless X-ray spectrum well described by a multicolour disc blackbody model. The source most likely represents an optically thick accretion disc around an intermediate-mass black hole  ( M ∼102 M⊙)  in its high/soft state, similar to other variable ultraluminous X-ray sources observed in nearby disc galaxies. Using constraints derived from both the innermost stable orbit around a black hole and the Eddington luminosity, we find that the black hole is fast-rotating and that its mass is between ∼80 M⊙–1.5×102 M⊙. The inferred bolometric luminosity of the accretion disc is ∼(1.1×1040 erg s−1)/(cos  i ). Furthermore, we find that the optical nebula is very energetic and may contain large amounts of hot gas, accounting for a soft X-ray component as indicated by archival ROSAT PSPC data. The nebula is apparently associated with X-9; the latter may be powering the former and/or they could be formed in the same event (e.g. a hypernova). Such a connection, if confirmed, could have strong implications for understanding both the birth of intermediate-mass black holes and the formation of energetic interstellar structures.  相似文献   

17.
We find that at redshifts   z ≳ 10, HD  line cooling allows strongly shocked primordial gas to cool to the temperature of the cosmic microwave background (CMB). This temperature is the minimum value attainable via radiative cooling. Provided that the abundance of HD, normalized to the total number density, exceeds a critical level of  ∼10−8  , the CMB temperature floor is reached in a time which is short in comparison to the Hubble time. We estimate the characteristic masses of stars formed out of shocked primordial gas in the wake of the first supernovae, and resulting from the virialization of dark matter haloes during hierarchical structure formation to be  ∼10 M⊙  . In addition, we show that cooling by HD enables the primordial gas in relic H  ii regions to cool to temperatures considerably lower than those reached via H2 cooling alone. We confirm that HD cooling is unimportant in cases where the primordial gas does not go through an ionized phase, as in the formation process of the very first stars in   z ≳ 20  minihaloes of mass  ∼106 M⊙  .  相似文献   

18.
We explore the implications of a possible cosmic-ray (CR) background generated during the first supernova explosions that end the brief lives of massive Population III stars. We show that such a CR background could have significantly influenced the cooling and collapse of primordial gas clouds in minihaloes around redshifts of   z ∼ 15–20  , provided the CR flux was sufficient to yield an ionization rate greater than about 10−19 s−1 near the centre of the minihalo. The presence of CRs with energies  ≲107  eV would indirectly enhance the molecular cooling in these regions, and we estimate that the resulting lower temperatures in these minihaloes would yield a characteristic stellar mass as low as  ∼10 M⊙  . CRs have a less-pronounced effect on the cooling and collapse of primordial gas clouds inside more massive dark matter haloes with virial masses  ≳108 M⊙  at the later stages of cosmological structure formation around   z ∼ 10–15  . In these clouds, even without CR flux the molecular abundance is already sufficient to allow cooling to the floor set by the temperature of the cosmic microwave background.  相似文献   

19.
It has been recently shown that the dynamical V -band mass-to-light ratios of compact stellar systems with masses from 106 to  108 M⊙  are not consistent with the predictions from simple stellar population models. Top-heavy stellar initial mass functions (IMFs) in these so-called ultra-compact dwarf galaxies (UCDs) offer an attractive explanation for this finding, the stellar remnants and retained stellar envelopes providing the unseen mass. We therefore construct a model which quantifies by how much the IMFs of UCDs would have to deviate in the intermediate- and high-mass range from the canonical IMF in order to account for the enhanced   M / LV   ratio of the UCDs. The deduced high-mass IMF in the UCDs depends on the age of the UCDs and the number of faint products of stellar evolution retained by them. Assuming that the IMF in the UCDs is a three-part power law equal to the canonical IMF in the low-mass range and taking 20 per cent as a plausible choice for the fraction of the remnants of high-mass stars retained by UCDs, the model suggests the exponent of the high-mass IMF to be ≈1.6 if the UCDs are  13 Gyr  old (i.e. almost as old as the Universe) or ≈1.0 if the UCDs are  7 Gyr  old, in contrast to 2.3 for the Salpeter–Massey IMF. If the IMF was as top heavy as suggested here, the stability of the UCDs might have been threatened by heavy mass loss induced by the radiation and evolution of massive stars. The central densities of UCDs must have been in the range  106 to 107 M⊙ pc−3  when they formed with star formation rates of  10 to 100 M⊙ yr−1  .  相似文献   

20.
We use the very large Millennium Simulation of the concordance Λ cold dark matter cosmogony to calibrate the bias and error distribution of Timing Argument estimators of the masses of the Local Group and of the Milky Way. From a large number of isolated spiral–spiral pairs similar to the Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to true mass to be a factor of 1.8, while the 5 and 95 per cent points of the distribution of this ratio are separated by a factor of 5.7. Here, we define true mass as the sum of the 'virial' masses, M 200, of the two dominant galaxies. For present best values of the distance and approach velocity of Andromeda, this leads to a median likelihood estimate of the true mass of the Local Group of  5.27 × 1012 M⊙  or  log  M LG/M⊙= 12.72  , with an interquartile range of [12.58, 12.83] and a 5–95 per cent range of [12.26, 13.01]. Thus, a 95 per cent lower confidence limit on the true mass of the Local Group is  1.81 × 1012 M⊙  . A timing estimate of the Milky Way's mass based on the large recession velocity observed for the distant satellite Leo I works equally well, although with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of  2.43 × 1012 M⊙  with a 95 per cent lower confidence limit of  0.80 × 1012 M⊙  .  相似文献   

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

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