Literature DB >> 21443262

Optical properties of nanocrystal interfaces in compressed MgO nanopowders.

Keith P McKenna1, David Koller, Andreas Sternig, Nicolas Siedl, Niranjan Govind, Peter V Sushko, Oliver Diwald.   

Abstract

The optical properties and charge trapping phenomena observed on oxide nanocrystal ensembles can be strongly influenced by the presence of nanocrystal interfaces. MgO powders represent a convenient system to study these effects due to the well-defined shape and controllable size distributions of MgO nanocrystals. The spectroscopic properties of nanocrystal interfaces are investigated by monitoring the dependence of absorption characteristics on the concentration of the interfaces in the nanopowders. The presence of interfaces is found to affect the absorption spectra of nanopowders more significantly than changing the size of the constituent nanocrystals and, thus, leading to the variation of the relative abundance of light-absorbing surface structures. We find a strong absorption band in the 4.0-5.5 eV energy range, which was previously attributed to surface features of individual nanocrystals, such as corners and edges. These findings are supported by complementary first-principles calculations. The possibility to directly address such interfaces by tuning the energy of excitation may provide new means for functionalization and chemical activation of nanostructures and can help improve performance and reliability for many nanopowder applications.

Entities:  

Year:  2011        PMID: 21443262      PMCID: PMC3082970          DOI: 10.1021/nn200062d

Source DB:  PubMed          Journal:  ACS Nano        ISSN: 1936-0851            Impact factor:   15.881


Powders of metal oxide nanocrystals find important applications in fields as diverse as catalysis,(1) medicine, cosmetics, and energy generation(2) owing to their unique and often tunable properties. A great deal of research has focused on understanding photoexcited processes in such nanopowders. When discussing how the spectroscopic properties and surface reactivity of nanopowders differ from those of macroscopic materials, it is common to point to their high specific surface area, the size and shape of the constituent nanocrystals, and the high concentration of surface features, such as steps, kinks, and other low coordination sites. Another defining characteristic of nanopowders is the large number of interfaces between nanocrystals they contain. However, the influence of interfaces on the overall optical properties of nanopowders has not so far been studied. As well as being an important issue for the aforementioned applications, understanding the role of interfaces is important for materials characterization, as in order to perform spectroscopic measurements, nanoparticle powders are often pressed into pellets, inevitably generating nanocrystal interfaces. Furthermore, understanding how interfaces affect optical properties may open the way to using optical spectroscopy as a probe of the electronic properties of interfaces between nanocrystals, which are thought to be important for effects such as magnetism in nanopowders,(3) and for understanding the dynamics of electrons and holes under irradiation.[4,5] In this article, we combine experiment and first-principles theoretical calculations to investigate the optical properties of interfaces in metal oxide nanopowders. MgO provides a convenient reference system for this study due to its simple crystal structure, high thermal stability, and the possibility to produce nanocrystals with well-defined size and shape.[6,7] Previous studies, using techniques such as electron paramagnetic resonance, optical spectroscopy, and first-principles calculations, have shown that MgO nanopowders exhibit optical absorption and luminescence bands that are red-shifted compared to the bulk.[6−11] For example, the bulk absorption threshold for MgO is 7.7 eV, but lower energy bands are also observed in nanopowders, which have been attributed to surface features such as corners (4.6 eV) and edges (5.2 eV).[7,9] The approach we have devised to investigate the optical properties of interfaces involves increasing their concentration by mechanical compaction of the nanopowder. We find two substantial changes in the UV diffuse reflectance spectra of MgO nanopowders resulting from compaction: the emergence of new absorption features in the range 4.0−5.5 eV and a depression in intensity at 5.7 eV. With the help of complementary first-principles calculations, these changes are interpreted in terms of the excitation of nanocrystal interfaces. An important conclusion of our study is that photons with low energy, e.g., 4.6 eV, excite electrons not only at MgO surfaces, as previously believed, but also at interfaces inside the powder. This is the first experimental evidence for the emergence of new electronic states that arise from the contact area between nanocrystals.

Results and Discussion

MgO nanopowders were prepared using a chemical vapor synthesis (CVS) procedure based on metal vapor combustion with oxygen within a flow reactor system (see Methods section for details). We showed previously(6) that by varying CVS parameters one can produce powders consisting of cubic nanocrystals with a well-defined size as small as 3 nm edge length. The, as-formed, powders were cleaned of organic contaminants by heating to 1123 K at a rate of 5 K min−1 and exposing to molecular oxygen at this temperature. Then, the sample temperature was raised to 1173 K at pressures p < 5 × 10−6 mbar and kept at this temperature for 5 h until full dehydroxylation of the sample surface was achieved. This procedure forms low-density powders (less than 1% of the bulk MgO density, ρMgO = 3580 kg m−3), and transmission electron microscopy (TEM) demonstrates they consist of loosely agglomerated cubic nanocrystals.[6,12] To obtain compressed MgO samples, a known mass of the powder is transferred into a 0.05 cm3 cavity and a hydraulic press is used to compact the powder. By varying the pressure applied, powders with a range of densities can be produced. The density of the compressed powder is determined from the known volume of the cavity after pressing. The pressed MgO pellet is then again subjected to a vacuum annealing procedure at 1173 K for 5 h, as well as annealing in 10 mbar of molecular oxygen at 1123 K for 10 min. We now describe how the properties of the MgO nanopowder are modified by compaction. We will focus on the difference between the as-formed powder (<1% ρMgO) and a powder compressed to 50% of the bulk MgO density. We will hereafter refer to these two samples as loose and dense MgO powders, respectively. X-ray diffraction (XRD) confirms that both loose and dense powder samples are crystallized in the cubic rock salt phase (1a). The average nanocrystal size was determined from reflex broadening using the Debye−Scherrer equation. The loose powder contains nanocrystals of average diameter, 7.6 ± 1 nm, which is increased to 11.1 ± 1 nm in the dense powder. For the latter, analysis of the TEM images (1b) shows that the edge length is most commonly distributed between 8 and 12 nm (Supporting Information, Figure S1), consistent with the XRD analysis. To determine the specific surface area (SBET), we measured N2 adsorption isotherms (Supporting Information, Figures. S2a and S2b) and applied the Brunauer, Emmett, and Teller (BET) model.(13) The analysis shows a reduction in specific surface area from 303 m2 g−1 for the loose powder to 108 m2 g−1 for the dense powder. Applying the Barret, Joyner, and Halenda model(14) to the N2 adsorption isotherms also allows one to obtain the pore size distributions. 1c shows that in the case of the loose powder we have a broad pore size distribution, which is consistent with a loosely connected powder with few nanocrystal interfaces. In the case of the dense powder, the pore distribution is peaked at approximately 10 nm, suggesting that a particle network containing numerous interfaces has been formed.
Figure 1

(a) XRD patterns of loose and dense MgO powder with corresponding average nanocrystal size. (b) Electron microscopy images of CVS-MgO after applying uniaxial pressure of 1.7 × 108 Pa and subsequent thermal annealing (T = 1170 K, p < 10−6 mbar). (c) Distribution in pore sizes for the loose and dense powder as determined from analysis of N2-adsorption isotherms. (d) Corrected specific surface area either derived from the average nanocrystal sizes or directly determined from the N2-adsorption isotherms.

(a) XRD patterns of loose and dense MgO powder with corresponding average nanocrystal size. (b) Electron microscopy images of CVS-MgO after applying uniaxial pressure of 1.7 × 108 Pa and subsequent thermal annealing (T = 1170 K, p < 10−6 mbar). (c) Distribution in pore sizes for the loose and dense powder as determined from analysis of N2-adsorption isotherms. (d) Corrected specific surface area either derived from the average nanocrystal sizes or directly determined from the N2-adsorption isotherms. Part of the 195 m2 g−1 reduction in specific surface area that is observed following compaction is a consequence of the increased nanocrystal size. Another important factor is the formation of areas of contact between nanocrystals in the compressed powder, which leads to increased interface area at the expense of surface area. Finally, internal voids that are inaccessible to N2 may be present, which would also contribute to the observed reduction. Given the relatively high porosity of the two powder samples, the contribution from internal voids is expected to be small and is neglected in the following analysis. To help separate the remaining two contributions, we have calculated the specific surface area of nanocrystals using the average edge length determined by XRD (dXRD). For this calculation, we assume isolated cubic nanocrystals and use the following expression: S = 6/dρMgO (1). Direct comparison between the XRD- and BET-determined surface areas is difficult due to the different types of approximation that are involved(15) (see Supporting Information). However, since for the loose powder we expect the assumption of isolated cubic nanocrystals to be valid, we can obtain a normalizing factor by equating SXRD and SBET for the loose powder. In this way, the difference (SXRD − SBET) provides a proportional measure of the interface area (1d, 1). From this analysis, we determine that the area associated with interfaces between nanocrystals is approximately 100 m2 g−1 in the dense powder. In other words, ∼50% of the nanocrystal surfaces in the dense powder are in direct contact with other nanocrystals.
Table 1

Average Nanocrystal Sizes Determined Using the Debye−Scherrer Equation (dXRD), Specific Surface Area Calculated from Average Nanocrystal Sizes (SXRD), Specific Surface Area Determined from N2 Adsorption (SBET), and the Resulting Interface Area

 loose powderdense powder
dXRD (nm)7.6 ± 111.1 ± 1
SXRD (m2 g−1)221151
normalized SXRD (m2 g−1)303207
SBET (m2 g−1)303 ± 30108 ± 11
interface area (SXRDSBET) (m2 g−1) 99
2a compares the UV diffuse reflectance spectra acquired for loose nanopowder samples (with dXRD = 7.6 nm and dXRD = 10.0 nm) and for a dense nanopowder sample (dXRD = 11.1 nm). The relative intensity of the absorption bands in these spectra reflect the relative proportion of the associated structural features responsible for the excitations. For example, for loose powders, previous studies have shown that the 4.6 and 5.2 eV bands are associated with low-coordinated corner and edge ions, respectively.[6,11] As the size of the nanocrystals increases, the concentration of corner ions is decreased, explaining why the absorption band at 4.6 eV is essentially absent for the 10 nm loose powder (bottom trace in 2a). However, slightly larger compressed MgO nanopowders exhibit UV absorption that extends significantly into the energy range below 5.0 eV, indicating substantial changes in the electronic structure induced by the introduction of nanocrystal interfaces. By taking the difference between the two spectra (10 nm loose and 11.1 nm dense), one can separate out the effects of interfaces. 2b shows that as the density increases, a band of excitations in the range 4.0−5.5 eV increases in intensity, while the intensity of a band at 5.7 eV decreases. The differences between the absorption spectra show that the effect of the interfaces in the range 4.0−5.5 eV is much stronger than the effect of the nanocrystal size (see 2 and Figure S4, Supporting Information). Residual surface hydroxyls can be excluded as the origin of the emerging bands, because—if present—such adsorbates would decrease the sample absorption in this energy range rather than increasing it (see Supporting Information). Some of the new bands introduced by interfaces overlap with bands that were previously thought to selectively excite only low-coordinated surface ions.
Figure 2

(a) Experimental UV diffuse reflectance spectra for loose and dense MgO powder samples. Average crystallite sizes (shown) were determined applying the Debye−Scherrer equation to XRD powder patterns. The spectra are displaced vertically for clarity. (b) Difference between the spectra obtained for 10.0 nm loose and 11.1 nm dense MgO powders.

(a) Experimental UV diffuse reflectance spectra for loose and dense MgO powder samples. Average crystallite sizes (shown) were determined applying the Debye−Scherrer equation to XRD powder patterns. The spectra are displaced vertically for clarity. (b) Difference between the spectra obtained for 10.0 nm loose and 11.1 nm dense MgO powders. On the basis of our experimental structural characterization, we adopt the model of powder densification shown in 3a.[16−19] The main effect of compression and annealing is to reduce the size of voids in the powder by forcing the nanocrystals to reorganize, forming a more dense network containing more interfaces in the process. This can involve a number of interrelated processes such as rotation and translation of nanocrystals, mass transfer between the nanocrystals resulting in size change, and growth of grain boundaries between nanocrystals. Following this transformation the constituent nanocrystals remain cubic and crystalline, and there is an increase in their average size, indicating some degree of sintering. It is difficult to extract the structure of interfaces present in the nanopowders directly from the experimental measurements. However, it is expected that the most common type of interface between cubic MgO nanocrystals will involve the commensurate contact between their (001) facets. The reasons for this are as follows: (1) commensurate interfaces are the most stable thermodynamically; (2) they seem to be common in TEM images of nanopowders;(20) and (3) previous TEM studies on larger MgO smoke particles deposited onto a MgO substrate found a clear preference for commensurate interfaces.(21) In this latter study, a smaller fraction of interfaces were also found to be rotated (corresponding to high site coincidence twist grain boundaries).
Figure 3

(a) Model for changes in powder structure on compression. (b) Structural features at the interface between nanocrystals that are commensurate or (c) rotated with respect to each other.

To provide atomistic models for interfaces in the MgO powders, we performed theoretical calculations of some possible metastable configurations of two adjoining MgO nanocrystals. These calculations considered nanocrystals containing 1000 ions and described the interionic bonding using the classical shell model.(22) The two cubic MgO nanocrystals were initially oriented and translated with respect to each other before relaxing the structure into various local energy minima. An alternative approach would be to employ molecular dynamics to simulate the sintering of nanocrystals, as has been done recently for TiO2.(23)3b,c shows examples of both a commensurate and a rotated interface that contain various structural features. For example, the corner−terrace feature is formed at the contact point between the corner of one nanocrystal and the (001) terrace of another. In addition to features formed at contact points between nanocrystals, the surface features of individual nanocrystals (i.e., the corner, edge, and terrace) are indicated. While, in general, the structure of nanopowders may be quite complicated, we expect they should be predominantly built from these fundamental structural features. Therefore, these models form the basis for the first-principles calculations that are used to help interpret the experimentally determined absorption spectra. (a) Model for changes in powder structure on compression. (b) Structural features at the interface between nanocrystals that are commensurate or (c) rotated with respect to each other. Excitation spectra have been calculated using time-dependent density functional theory (TD-DFT) implemented within an embedded cluster method. To calculate excitation spectra for the various features identified in 3b,c, quantum clusters were constructed that were sufficiently large to ensure accurate description of the surface and interface electronic states participating in the electronic transitions. The largest system we considered contained 73 atoms treated at an all-electron level (730 electrons in total), spanning a region of ∼1 nm3, which was described using 1723 basis functions (see Supporting Information for more details). For this cluster, we computed 50 excited states using TD-DFT, a calculation that was possible only through the use of the highly parallelized NWChem code.(24) The excitation spectra of the loose powder is simulated by assuming it to be dominated by the excitation of corner, edge, and terrace features. To obtain the total excitation spectra, we sum up the individual contributions weighted in proportion to the number of corner, edge, and terrace sites in the powder. One can either determine these weights from atomistic models or use them as fitting parameters if comparing to an experimental spectrum. 4a (dashed curve) shows the simulated spectra of an isolated 2 nm nanocrystal, where we used the proportion of various surface sites appropriate for cubic nanocrystals.(6) In this spectra, one can identify a very weak feature at 3.7 eV and a more intense peak at 4.7 eV. Both of these features are associated with localized excitations at corner ions. At higher energies, contributions from edges in the range 5.0−5.5 eV partly overlap with terrace excitations peaking near 5.8 eV. This simulated spectrum provides very good agreement with the measured spectrum for the loose powder (2a). To test our hypothesis that the experimentally observed difference spectra are due to the increased number of interfaces in the dense powder, we calculate the excitation spectra of various structural features formed at the interface between MgO nanocrystals, including corner−terrace, edge−terrace, and edge−intersection features (as shown in 3b,c). All interfaces considered contribute excitations in the range 4.4−5.4 eV (4b), consistent with the experimental result. Furthermore, as the density increases, the total area of free terrace is reduced (as shown by N2-sorption measurements); therefore, one should expect a reduction in the intensity associated with terraces. The calculations predict that this depression will occur at 5.8 eV (i.e., the position of the peak in absorption due to terraces, 2a, again in close agreement with experiment (5.7 eV)).
Figure 4

(a) Theoretically simulated spectra for loose powder assuming dominant contributions from corners, edges, and terraces and (b) excitation spectra for various interface features formed between commensurate and rotated (r) nanocrystals.

On the basis of the combined experimental and theoretical investigation, we find that a wide range of interfaces in MgO nanopowders, even those that are highly coordinated, can contribute low-energy excitations in the 4.0−5.5 eV range. The theoretical calculations explain the origin of this effect. At interfaces, there is a significant perturbation in the structure, and as a consequence, the Madelung potential on interface ions can be significantly reduced. For example, the Madelung potential near the corner−terrace feature is up to 0.5 V lower than at the free edge. It is this electrostatic perturbation that is responsible for the presence of interfacial electronic states and gives rise to excitation energies significantly lower than those of isolated particles. For example, one of the dominant excitations for the edge−intersection feature at 4.7 eV involves transition of an electron from oxygen p-states on the edge of one cube to the edge of the adjacent cube. This novel type of interparticle excitation may have interesting consequences, for example, leading to the development of nonequilibrium electric fields within nanopowders under irradiation. It has been proposed that such electric fields may weaken chemical bonds in small molecules and, therefore, facilitate chemical reactions between them.[25,26] (a) Theoretically simulated spectra for loose powder assuming dominant contributions from corners, edges, and terraces and (b) excitation spectra for various interface features formed between commensurate and rotated (r) nanocrystals. In summary, we have shown, by monitoring the optical properties of a nanopowder as the density of the powder is increased, how one can determine the excitation/absorption characteristics of interfaces. While we have considered MgO as a model system, the results suggest that understanding exciton generation at interfaces, and the associated electron−hole dynamics, may be important for other materials and applications. By combining theory and experiment one can identify regions of the spectra that are associated with particular types of interface. The possibility to directly address such interfaces by tuning the energy of excitation is key to understanding exciton generation at interfaces and associated charge separation. Moreover, it may provide new means for functionalization of nanostructures and chemical activation and can help improve performance and reliability for many nanopowder applications, such as photocatalysts and solar cells. For example, there is considerable interest in porous nanocrystalline ceramics for catalytic applications,[27,28] batteries, solar cells,(29) and gas sensing, where interfaces may play an important role. The key result of this study is that the effect of interfaces in powders of ionic nanocrystals on the optical properties may be as important as those specific for low-coordinated adsorption sites.

Methods

For the production of MgO nanocrystals we use chemical vapor synthesis, which allows for the controlled evaporation and subsequent oxidation of alkaline earth metals under reduced pressures.[30,31] Further details can be found in the Supporting Information. After production, the MgO nanocrystal powder is transferred into quartz glass cells, which allow one to carry out thermal activation of the nanocrystal powders in defined gas atmospheres on one hand and spectroscopy measurements on the other. To obtain well-defined cubic MgO nanocrystals, the as-obtained MgO powder is cleaned of organic contaminants by heating to 1123 K at a rate of 5 K min−1 and exposing to molecular oxygen at this temperature. Then, the sample temperature was raised to 1173 K at pressures p < 5 × 10−6 mbar and kept at this temperature for 5 h until full dehydroxylation of the sample surface was achieved. X-ray diffraction patterns were collected on a PANalytical X'Pert PRO diffractometer using Cu Kα radiation. The nanocrystal sizes were derived on the basis of the Debye−Scherrer equation with the assumption of cubically shaped nanocrystals. We used pseudo-Voigt functions to determine the full width half-maximum of the three main reflexes to calculate an average nanocrystal size (dXRD). To obtain the transmission electron microscopy images, small amounts of the metal oxide powders were cast on a holey carbon grid for investigation with a TECNAI F20 analytical transmission electron microscope equipped with a field emission gun and an S-twin objective lens. Nitrogen adsorption isotherms were obtained at 77 K using an adsorption porosimeter (Micromeritics ASAP 2020). Samples were outgassed for 6 h in the degas unit of the adsorption apparatus at 473 K under vacuum prior to analysis. The BET surface area (SBET) was evaluated using adsorption data in a relative pressure range p/p0 of 0.05−0.2. UV diffuse reflectance spectra for the loose and dense powders were acquired in the presence of 10 mbar of O2 using a PerkinElmer Lambda 15 spectrophotometer, equipped with an integrating sphere, and then converted to absorption spectra by the Kubelka−Munk transform procedure. There are differences in the absolute light penetration depth due to the different densities of the samples; therefore, reported spectra are normalized with respect to their maximum reflectance value. Excitation spectra have been calculated using time-dependent density functional theory implemented within an embedded cluster method. This technique consistently combines quantum mechanical and classical methods to enable large systems to be studied and to ensure accurate description of long-range polarization effects. Classical ions are described using empirical pair potentials,(22) and quantum mechanical ions are treated using the B3LYP hybrid density functional(32) and the Gaussian 6-311G* basis set. This method is implemented in the Guess code,(10) which was interfaced with the NWChem code[24,33] for the quantum mechanical part of the calculation.
  12 in total

1.  Local Structure Evolution in Particle Network Formation Studied by Brownian Dynamics Simulation.

Authors: 
Journal:  J Colloid Interface Sci       Date:  2000-11-15       Impact factor: 8.128

2.  Dye-sensitized solar cells.

Authors:  Anders Hagfeldt; Gerrit Boschloo; Licheng Sun; Lars Kloo; Henrik Pettersson
Journal:  Chem Rev       Date:  2010-09-10       Impact factor: 60.622

3.  Optical surface properties and morphology of MgO and CaO nanocrystals.

Authors:  Slavica Stankic; Johannes Bernardi; Oliver Diwald; Erich Knözinger
Journal:  J Phys Chem B       Date:  2006-07-20       Impact factor: 2.991

4.  The mechanism of dihydrogen activation by frustrated Lewis pairs revisited.

Authors:  Stefan Grimme; Holger Kruse; Lars Goerigk; Gerhard Erker
Journal:  Angew Chem Int Ed Engl       Date:  2010-02-15       Impact factor: 15.336

5.  Electron-trapping polycrystalline materials with negative electron affinity.

Authors:  Keith P McKenna; Alexander L Shluger
Journal:  Nat Mater       Date:  2008-10-12       Impact factor: 43.841

6.  The strange magnetism of oxides and carbons.

Authors:  Marshall Stoneham
Journal:  J Phys Condens Matter       Date:  2010-02-03       Impact factor: 2.333

7.  Understanding the surface potential of water.

Authors:  Shawn M Kathmann; I-Feng William Kuo; Christopher J Mundy; Gregory K Schenter
Journal:  J Phys Chem B       Date:  2011-03-30       Impact factor: 2.991

8.  BaO clusters on MgO nanocubes: a quantitative analysis of optical-powder properties.

Authors:  Andreas Sternig; Markus Müller; Mark McCallum; Johannes Bernardi; Oliver Diwald
Journal:  Small       Date:  2010-02-22       Impact factor: 13.281

9.  Engineering of efficient panchromatic sensitizers for nanocrystalline TiO(2)-based solar cells.

Authors:  M K Nazeeruddin; P Péchy; T Renouard; S M Zakeeruddin; R Humphry-Baker; P Comte; P Liska; L Cevey; E Costa; V Shklover; L Spiccia; G B Deacon; C A Bignozzi; M Grätzel
Journal:  J Am Chem Soc       Date:  2001-02-28       Impact factor: 15.419

10.  Molecular dynamics simulation of the aggregation of titanium dioxide nanocrystals: preferential alignment.

Authors:  Mozhgan Alimohammadi; Kristen A Fichthorn
Journal:  Nano Lett       Date:  2009-12       Impact factor: 11.189

View more
  4 in total

1.  Surface Decoration of MgO Nanocubes with Sulfur Oxides: Experiment and Theory.

Authors:  Andreas Sternig; Oliver Diwald; Silvia Gross; Peter V Sushko
Journal:  J Phys Chem C Nanomater Interfaces       Date:  2013-03-21       Impact factor: 4.126

2.  Electronic and chemical properties of a surface-terminated screw dislocation in MgO.

Authors:  Keith P McKenna
Journal:  J Am Chem Soc       Date:  2013-12-05       Impact factor: 15.419

3.  A Label-Free Photoluminescence Genosensor Using Nanostructured Magnesium Oxide for Cholera Detection.

Authors:  Manoj Kumar Patel; Md Azahar Ali; Sadagopan Krishnan; Ved Varun Agrawal; AbdulAziz A Al Kheraif; H Fouad; Z A Ansari; S G Ansari; Bansi D Malhotra
Journal:  Sci Rep       Date:  2015-11-27       Impact factor: 4.379

4.  Thin water films and particle morphology evolution in nanocrystalline MgO.

Authors:  Daniel Thomele; Amir R Gheisi; Matthias Niedermaier; Michael S Elsässer; Johannes Bernardi; Henrik Grönbeck; Oliver Diwald
Journal:  J Am Ceram Soc       Date:  2018-05-30       Impact factor: 3.784

  4 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.