首页 | 官方网站   微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We present an analysis of the geoeffectiveness of corotating interaction regions (CIRs), employing the data recorded from 25 January to 5 May 2005 and throughout 2008. These two intervals in the declining phase of Solar Cycle 23 are characterised by a particularly low number of interplanetary coronal mass ejections (ICMEs). We study in detail how four geomagnetic-activity parameters (the Dst, Ap, and AE indices, as well as the Dst time derivative, \(\mathrm{dDst}/\mathrm{d}t\)) are related to three CIR-related solar wind parameters (flow speed, \(V\), magnetic field, \(B\), and the convective electric field based on the southward Geocentric solar magnetospheric (GSM) magnetic field component, \(\mathit{VB}_{s}\)) on a three-hour time resolution. In addition, we quantify statistical relationships between the mentioned geomagnetic indices. It is found that Dst is correlated best to \(V\), with a correlation coefficient of \(\mathrm{cc}\approx0.6\), whereas there is no correlation between \(\mathrm{dDst}/\mathrm{d}t\) and \(V\). The Ap and AE indices attain peaks about half a day before the maximum of \(V\), with correlation coefficients ranging from \(\mathrm{cc}\approx0.6\) to \(\mathrm{cc}\approx0.7\), depending on the sample used. The best correlations of Ap and AE are found with \(\mathit{VB}_{s}\) with a delay of 3 h, being characterised by \(\mathrm{cc}\gtrsim 0.6\). The Dst derivative \(\mathrm{dDst}/\mathrm{d}t\) is also correlated with \(\mathit{VB}_{s}\), but the correlation is significantly weaker \(\mathrm{cc}\approx 0.4\)?–?0.5, with a delay of 0?–?3 h, depending on the employed sample. Such low values of correlation coefficients indicate that there are other significant effects that influence the relationship between the considered parameters. The correlation of all studied geomagnetic parameters with \(B\) are characterised by considerably lower correlation coefficients, ranging from \(\mathrm{cc}=0.3\) in the case of \(\mathrm{dDst}/\mathrm{d}t\) up to \(\mathrm{cc}=0.56\) in the case of Ap. It is also shown that peak values of geomagnetic indices depend on the duration of the CIR-related structures. The Dst is closely correlated with Ap and AE (\(\mathrm{cc}=0.7\)), Dst being delayed for about 3 h. On the other hand, \(\mathrm{dDst}/\mathrm{d}t\) peaks simultaneously with Ap and AE, with correlation coefficients of 0.48 and 0.56, respectively. The highest correlation (\(\mathrm{cc}=0.81\)) is found for the relationship between Ap and AE.  相似文献   

2.
In this paper, we explore the possibility of accreting primordial black holes as the source of heating for the collapsing gas in the context of the direct collapse black hole scenario for the formation of super-massive black holes (SMBHs) at high redshifts, \(z\sim \) 6–7. One of the essential requirements for the direct collapse model to work is to maintain the temperature of the in-falling gas at \(\approx \)10\(^4\) K. We show that even under the existing abundance limits, the primordial black holes of masses \(\gtrsim \)10\(^{-2}M_\odot \), can heat the collapsing gas to an extent that the \(\mathrm{H}_2\) formation is inhibited. The collapsing gas can maintain its temperature at \(10^4\) K till the gas reaches a critical density \(n_{{c}} \,{\approx }\, 10^3~\hbox {cm}^{-3}\), at which the roto-vibrational states of \(\mathrm{H}_2\) approaches local thermodynamic equilibrium and \(\mathrm{H}_2\) cooling becomes inefficient. In the absence of \(\mathrm{H}_2\) cooling, the temperature of the collapsing gas stays at \(\approx \)10\(^4\) K even as it collapses further. We discuss scenarios of subsequent angular momentum removal and the route to find collapse through either a supermassive star or a supermassive disk.  相似文献   

3.
We examine the properties of the viscous dissipative accretion flow around rotating black holes in the presence of mass loss. Considering the thin disc approximation, we self-consistently calculate the inflow-outflow solutions and observe that the mass outflow rates decrease with the increase in viscosity parameter (\(\alpha \)). Further, we carry out the model calculation of quasi-periodic oscillation frequency (\(\nu _{\mathrm{QPO}}\)) that is frequently observed in black hole sources and observe that \(\nu ^\mathrm{max}_{\mathrm{QPO}}\) increases with the increase of black hole spin (\(a_k\)). Then, we employ our model in order to explain the High Frequency Quasi-Periodic Oscillations (HFQPOs) observed in black hole source GRO J1655-40. While doing this, we attempt to constrain the range of \(a_k\) based on observed HFQPOs (\(\sim \)300 Hz and \(\sim \)450 Hz) for the black hole source GRO J1655-40.  相似文献   

4.
We use a formulation of the N-body problem in spaces of constant Gaussian curvature, \({\kappa }\in \mathbb {R}\), as widely used by A. Borisov, F. Diacu and their coworkers. We consider the restricted three-body problem in \(\mathbb {S}^2\) with arbitrary \({\kappa }>0\) (resp. \(\mathbb {H}^2\) with arbitrary \({\kappa }<0\)) in a formulation also valid for the case \({\kappa }=0\). For concreteness when \({\kappa }>0\) we restrict the study to the case of the three bodies at the upper hemisphere, to be denoted as \(\mathbb {S}^2_+\). The main goal is to obtain the totality of relative equilibria as depending on the parameters \({\kappa }\) and the mass ratio \(\mu \). Several general results concerning relative equilibria and its stability properties are proved analytically. The study is completed numerically using continuation from the \({\kappa }=0\) case and from other limit cases. In particular both bifurcations and spectral stability are also studied. The \(\mathbb {H}^2\) case is similar, in some sense, to the planar one, but in the \(\mathbb {S}^2_+\) case many differences have been found. Some surprising phenomena, like the coexistence of many triangular-like solutions for some values \(({\kappa },\mu )\) and many stability changes will be discussed.  相似文献   

5.
In this paper we have studied the anisotropic Kantowski-Sachs, locally rotationally symmetric (LRS) Bianchi type-I and LRS Bianchi type-III geometries filled with dark energy and one dimensional cosmic string in the Saez-Ballester theory of gravitation. To get physically valid solution we take hybrid expansion law of the average scale factor which is a product of power and exponential type of functions that results in time dependent deceleration parameter (\(q\)). The equation of state parameter of dark energy (\(\omega _{\mathit{de}}\)) has been discussed and we have observed that for the three models it crosses the phantom divide line (\(\omega _{\mathit{de}} = -1\)) and shows quintom like behavior. The density of dark energy (\(\rho _{\mathit{de}}\)) is an increasing function of redshift and remains positive throughout the evolution of the universe for the three models. Moreover in Kantowski-Sachs and LRS Bianchi type-I geometries the dark energy density dominates the string tension density (\(\lambda \)) and proper density (\(\rho \)) throughout the evolution of the universe. The physical and geometrical aspects of the statefinder parameters (\(r,s\)), squared speed of sound (\(v_{s}^{2} \)) and \(\omega _{\mathit{de}}\)\(\omega ^{\prime }_{\mathit{de}}\) plane are also discussed.  相似文献   

6.
Profile variations in the \(\hbox {H}\alpha \) and \(\hbox {H}\beta \) lines in the spectra of the star HD14134 are investigated using observations carried out in 2013–2014 and 2016 with the 2-m telescope at the Shamakhy Astrophysical Observatory. The absorption and emission components of the \(\hbox {H}\alpha \) line are found to disappear on some observational days, and two of the spectrograms exhibit inverse P-Cyg profile of \(\hbox {H}\alpha \). It was revealed that when the \(\hbox {H}\alpha \) line disappeared or an inversion of the P-Cyg-type profile is observed in the spectra, the \(\hbox {H}\beta \) line is displaced to the longer wavelengths, but no synchronous variabilities were observed in other spectral lines (CII \( \lambda \) 6578.05 Å, \( \lambda \) 6582.88 Å  and HeI \( \lambda \) 5875.72 Å) formed in deeper layers of the stellar atmosphere. In addition, the profiles of the \(\hbox {H}\alpha \) and \(\hbox {H}\beta \) lines have been analysed, as well as their relations with possible expansion, contraction and mixed conditions of the atmosphere of HD14134. We suggest that the observational evidence for the non-stationary atmosphere of HD14134 can be associated in part with the non-spherical stellar wind.  相似文献   

7.
We investigate the parameters of global solar p-mode oscillations, namely damping width \(\Gamma\), amplitude \(A\), mean squared velocity \(\langle v^{2}\rangle\), energy \(E\), and energy supply rate \(\mathrm{d}E/\mathrm{d}t\), derived from two solar cycles’ worth (1996?–?2018) of Global Oscillation Network Group (GONG) time series for harmonic degrees \(l=0\,\mbox{--}\,150\). We correct for the effect of fill factor, apparent solar radius, and spurious jumps in the mode amplitudes. We find that the amplitude of the activity-related changes of \(\Gamma\) and \(A\) depends on both frequency and harmonic degree of the modes, with the largest variations of \(\Gamma\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le3300~\upmu\mbox{Hz}\) and \(31\le l \le60\) with a minimum-to-maximum variation of \(26.6\pm0.3\%\) and of \(A\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le 3300~\upmu\mbox{Hz}\) and \(61\le l \le100\) with a minimum-to-maximum variation of \(27.4\pm0.4\%\). The level of correlation between the solar radio flux \(F_{10.7}\) and mode parameters also depends on mode frequency and harmonic degree. As a function of mode frequency, the mode amplitudes are found to follow an asymmetric Voigt profile with \(\nu_{\text{max}}=3073.59\pm0.18~\upmu\mbox{Hz}\). From the mode parameters, we calculate physical mode quantities and average them over specific mode frequency ranges. In this way, we find that the mean squared velocities \(\langle v^{2}\rangle\) and energies \(E\) of p modes are anticorrelated with the level of activity, varying by \(14.7\pm0.3\%\) and \(18.4\pm0.3\%\), respectively, and that the mode energy supply rates show no significant correlation with activity. With this study we expand previously published results on the temporal variation of solar p-mode parameters. Our results will be helpful to future studies of the excitation and damping of p modes, i.e., the interplay between convection, magnetic field, and resonant acoustic oscillations.  相似文献   

8.
X-ray stars have been studied since the beginning of X-ray astronomy. Investigating and studying the chromospheric activity from X-ray stellar optical spectra is highly significant in providing insights into stellar magnetic activity. The big data of LAMOST survey provides an opportunity for researching stellar optical spectroscopic properties of X-ray stars. We inferred the physical properties of X-ray stellar sources from the analysis of LAMOST spectra. First, we cross-matched the X-ray stellar catalogue (12254 X-ray stars) from ARXA with LAMOST data release 3 (DR3), and obtained 984 good spectra from 713 X-ray sources. We then visually inspected and assigned spectral type to each spectrum and calculated the equivalent width (EW) of H\(\alpha\) line using the Hammer spectral typing facility. Based on the EW of H\(\alpha\) line, we found 203 spectra of 145 X-ray sources with H\(\alpha\) emission above the continuum. For these spectra we also measured the EWs of H\(\beta\), H\(\gamma\), H\(\delta\) and Ca ii IRT lines of these spectra. After removing novae, planetary nebulae and OB-type stars, we found there are 127 X-ray late-type stars with H\(\alpha\) line emission. By using our spectra and results from the literature, we found 53 X-ray stars showing H\(\alpha\) variability; these objects are Classical T Tauri stars (CTTs), cataclysmic variables (CVs) or chromospheric activity stars. We also found 18 X-ray stars showing obvious emissions in the Ca ii IRT lines. Of the 18 X-ray stars, 16 are CTTs and 2 are CVs. Finally, we discussed the relationships between the EW of H\(\alpha\) line and X-ray flux.  相似文献   

9.
The kinetic Alfven waves in the presence of homogeneous magnetic field plasma with multi-ions effect are investigated. The dispersion relation and normalised damping rate are derived for low-\(\beta\) plasma using kinetic theory. The effect of density variation of \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions is observed on frequency and damping rate of the wave. The variation of frequency (\(\omega\)) and normalised damping rate (\(\gamma / \varOmega_{H^{ +}} \)) of the wave are studied with respect to \(k_{ \bot} \rho_{j}\), where \(k_{ \bot} \) is the perpendicular wave number, \(\rho_{j}\) is the ion gyroradius and \(j \) denotes \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions. The variation with \(k_{ \bot} \rho_{j}\) is considered over wide range. The parameters appropriate to cusp region are used for the explanation of results. It is found that with hydrogen and helium ions gyration, the frequency of wave is influenced by the density variation of \(\text{H}^{+}\) and \(\text{He}^{+}\) ions but remains insensitive to the change in density of \(\text{O}^{+}\) ions. For oxygen ion gyration, the frequency of wave varies over a short range only for \(\text{O}^{+}\) ion density variation. The wave shows damping at lower altitude due to variation in density of lighter \(\text{H}^{+}\) and \(\text{He}^{+}\) ions whereas at higher altitude only heavy \(\text{O}^{+}\) ions contribute in wave damping. The damping of wave may be due to landau damping or energy transfer from wave to particles. The present study signifies that the both lighter and heavier ions dominate differently to change the characteristics of kinetic Alfven wave and density variation is also an important parameter to understand wave phenomena in cusp region.  相似文献   

10.
It is first proposed a theoretical scaling law respectively for the coronal magnetic field strength \(B\) and electron power-law index \(\delta\) versus frequency and coronal height in solar microwave burst sources. Based on the non-thermal gyro-synchrotron radiation model (Ramaty in Astrophys. J. 158:753, 1969), \(B\) and \(\delta\) are uniquely solved by the observable optically-thin spectral index and turnover (peak) frequency, the other parameters (plasma density, temperature, view angle, low and high energy cutoffs, etc.) are relatively insensitive to the calculations, thus taken as some typical values. Both of \(B\) and \(\delta\) increase with increasing of radio frequency but with decreasing of coronal height above photosphere, and well satisfy a square or cubic logarithmic fitting.  相似文献   

11.
In this note, by using Smale’s \(\alpha \)-theorem on the convergence of Newton’s method, the \(\alpha \)-sets of convergence of some starters of solving the elliptic Kepler’s equation are derived. For each starter we compute the exact \(\alpha \)-set in the eccentricity-main anomaly \((e,M)\in [0,1)\times [0,\pi ]\), showing that these sets are larger than those derived by Avendaño et al. (Celest Mech Dyn Astron 119:27–44, 2014). Further, new convergence tests based on the Newton–Kantorowitch theorem are given comparing with the derived from Smale’s \(\alpha \)-test.  相似文献   

12.
In this paper, multiwavelength chromatic luminosity at radio (\(\log L _{\mathrm{R}}\)), optical (\(\log L_{\mathrm{O}}\)), X-ray (\(\log L _{\mathrm{X}}\)), and \(\gamma \)-rays (\(\log L_{\gamma }\)) for a sample of 442 Fermi blazars with known redshifts are collected from Fan et al. (2016), to study the correlations between the \(\gamma \)-rays and the low-energy bands using a multiple linear regression analysis. In this way, we can see which band is more important for the \(\gamma \)-ray emissions. Mutual correlation analysis is also used to discuss the correlations between the \(\gamma \)-ray and the low energy bands for the whole sample and subclasses. We come to following conclusions:
  1. 1.
    The multiple linear correlation indicates that the \(\gamma \)-rays are correlated with the radio, optical and the X-ray emissions for the whole sample and the subclasses of flat spectrum radio quasars (FSRQs) and BL Lac objects (BL Lacs), the correlation between the \(\gamma \)-rays and the radio emissions is the strongest one.
     
  2. 2.
    For BL Lacs, the optical emissions are more important than the X-rays for the \(\gamma \)-rays, while the X-ray emissions are more important than optical ones in FSRQs.
     
  3. 3.
    The \(\gamma \)-ray emissions in HBL are from an synchrotron self-Compton, while those in FSRQs may be from external Compton and synchrotron self-Compton as well.
     
  相似文献   

13.
In this paper, we compile the very-high-energy and high-energy spectral indices of 43 BL Lac objects from the literature. Based on a simple math model, \(\Delta \Gamma_{obs}= \alpha {{{z}}}+\beta \), we present evidence for the origin of an observed spectral break that is denoted by the difference between the observed very-high-energy and high-energy spectral indices, \(\Delta \Gamma_{obs}\). We find by linear regression analysis that \(\alpha \ne 0\) and \(\beta \ne 0\). These results suggest that the extragalactic background light attenuation and the intrinsic curvature dominate on the GeV-TeV \(\gamma \)-ray energy spectral break of BL Lac objects. We argue that the extragalactic background light attenuation is an exclusive explanation for the redshift evolution of the observed spectral break.  相似文献   

14.
We have studied the variability of S5 0716+714 at radio 15 GHz and \(\gamma\)-ray band using three different methods. A possible periodicity of \(P_{15~\text{GHz}}=266.0\pm11.5\) and \(P_{\gamma}=344.0 \pm16.4\) days are obtained for radio 15 GHz and \(\gamma\)-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio 15 GHz and \(\gamma \)-ray may be due to that the emission of radio 15 GHz is produced via the synchrotron process, while the \(\gamma\)-ray is produced by both the SSC and EC processes.  相似文献   

15.
The physical parameters of the solar wind observed in-situ near 1 AU have been studied for several decades, and relationships between them, such as the positive correlation between the solar wind plasma temperature, \(T\), and velocity, \(V\), and the negative correlation between density, \(N\), and velocity, \(V\), are well known. However, the magnetic field intensity, \(B\), does not appear to be well correlated with any individual plasma parameter. In this article, we discuss previously under-reported correlations between \(B\) and the combined plasma parameters \(\sqrt{N V^{2}} \) as well as between \(B\) and \(\sqrt{NT}\). These two correlations are strong during periods of corotating interaction regions and high-speed streams, and moderate during intervals of slow solar wind. The results indicate that the magnetic pressure in the solar wind is well correlated both with the plasma dynamic pressure and the thermal pressure.  相似文献   

16.
In this paper we address an \(n+1\)-body gravitational problem governed by the Newton’s laws, where n primary bodies orbit on a plane \(\varPi \) and an additional massless particle moves on the perpendicular line to \(\varPi \) passing through the center of mass of the primary bodies. We find a condition for the described configuration to be possible. In the case when the primaries are in a rigid motion, we classify all the motions of the massless particle. We study the situation when the massless particle has a periodic motion with the same minimal period as the primary bodies. We show that this fact is related to the existence of a certain pyramidal central configuration.  相似文献   

17.
This work is a continuation of our previous articles (Yermolaev et al. in J. Geophys. Res.120, 7094, 2015 and Yermolaev et al. in Solar Phys.292, 193, 2017), which describe the average temporal profiles of interplanetary plasma and field parameters in large-scale solar-wind (SW) streams: corotating interaction regions (CIRs), interplanetary coronal mass ejections (ICMEs, including both magnetic clouds (MCs) and ejecta), and sheaths as well as interplanetary shocks (ISs). Changes in the longitude angle, \(\varphi\), in CIRs from ?2 to \(2^{\circ}\) agree with earlier results (e.g. Gosling and Pizzo, 1999). We have also analyzed the average temporal profiles of the bulk velocity angles in sheaths and ICMEs. We have found that the angle \(\varphi\) in ICMEs changes from 2 to \(-2^{\circ}\), while in sheaths it changes from ?2 to \(2^{\circ}\) (similar to the change in CIRs), i.e. the angle in CIRs and sheaths deflects in the opposite sense to ICMEs. When averaging the latitude angle \(\vartheta\) on all the intervals of the chosen SW types, the angle \(\vartheta\) is almost constant at \({\sim}\,1^{\circ}\). We made for the first time a selection of SW events with increasing and decreasing \(\vartheta\) and found that the average \(\vartheta\) temporal profiles in the selected events have the same “integral-like” shape as for \(\varphi\). The difference in \(\varphi\) and \(\vartheta\) average profiles is explained by the fact that most events have increasing profiles for the angle in the ecliptic plane as a result of solar rotation, while for the angle in the meridional plane, the numbers of events with increasing and decreasing profiles are equal.  相似文献   

18.
We use the data for the \(\text{H}\beta\) emission-line, far-ultraviolet (FUV) and mid-infrared 22 μm continuum luminosities to estimate star formation rates \(\langle \mbox{SFR} \rangle \) averaged over the galaxy lifetime for a sample of about 14000 bursting compact star-forming galaxies (CSFGs) selected from the Data Release 12 (DR12) of the Sloan Digital Sky Survey (SDSS). The average coefficient linking \(\langle \mbox{SFR} \rangle \) and the star formation rate \(\mbox{SFR}_{0}\) derived from the \(\text{H}\beta\) luminosity at zero starburst age is found to be 0.04. We compare \(\langle \mbox{SFR} \rangle \mbox{s}\) with some commonly used SFRs which are derived adopting a continuous star formation during a period of \({\sim}\,100~\mbox{Myr}\), and find that the latter ones are 2–3 times higher. It is shown that the relations between SFRs derived using a geometric mean of two star-formation indicators in the UV and IR ranges and reduced to zero starburst age have considerably lower dispersion compared to those with single star-formation indicators. We suggest that our relations for \(\langle \mbox{SFR} \rangle \) determination are more appropriate for CSFGs because they take into account a proper temporal evolution of their luminosities. On the other hand, we show that commonly used SFR relations can be applied for approximate estimation within a factor of \({\sim}\,2\) of the \(\langle \mbox{SFR} \rangle \) averaged over the lifetime of the bursting compact galaxy.  相似文献   

19.
In extremely dense neutrino environments like in supernova core, the neutrino-neutrino refraction may give rise to self-induced flavor conversion. These neutrino flavor oscillations are well understood from the idea of the exponentially growing modes of the interacting oscillators in the flavor space. Until recently, the growth rates of these modes were found to be of the order of the vacuum oscillation frequency \(\Delta m^2/2E\) [\(\mathcal {O}(1~\mathrm{km}^{-1})\)] and were considered slow growing. However, in the last couple of years it was found that if the system was allowed to have different zenith-angle distributions for the emitted \(\nu _e\) and \(\bar{\nu }_e\) beams then the fastest growing modes of the interacting oscillators grew at the order of \(\mu =\sqrt{2} G_\mathrm{F}n_{\nu }\), a typical \(\nu \)\(\nu \) interaction energy [\(\mathcal {O}(10^5~\mathrm{km}^{-1})\)]. Thus the growth rates are very large in comparison to the so-called ‘slow oscillations’ and can result in neutrino flavor conversion on a much faster scale. In fact, the point that the growth rates are no longer dependent on the vacuum oscillation frequency \(\Delta m^2/2E\), makes these ‘fast flavor conversions’ independent of \(\Delta m^2\) (thus mass) and energy. This is a surprising result as neutrino flavor conversions are considered to be the ultimate proof of massive neutrinos. However, the importance of this effect in the realistic astrophysical scenarios still remains to be understood.  相似文献   

20.
The aim of this paper is to study new holographic dark energy (HDE) model in modified \(f(R,T)\) gravity theory within the framework of a flat Friedmann-Robertson-Walker model with bulk viscous matter content. It is thought that the negative pressure caused by the bulk viscosity can play the role of dark energy component, and drive the accelerating expansion of the universe. This is the motive of this paper to observe such phenomena with bulk viscosity. In the specific model \(f(R,T)=R+\lambda T\), where \(R\) is the Ricci scalar, \(T\) the trace of the energy-momentum tensor and \(\lambda \) is a constant, we find the solution for non-viscous and viscous new HDE models. We analyze new HDE model with constant bulk viscosity, \(\zeta =\zeta _{0}= \text{const.}\) to explain the present accelerated expansion of the universe. We classify all possible scenarios (deceleration, acceleration and their transition) with possible positive and negative ranges of \(\lambda \) over the constraint on \(\zeta _{0}\) to analyze the evolution of the universe. We obtain the solutions of scale factor and deceleration parameter, and discuss the evolution of the universe. We observe the future finite-time singularities of type I and III at a finite time under certain constraints on \(\lambda \). We also investigate the statefinder and \(\mathit{Om}\) diagnostics of the viscous new HDE model to discriminate with other existing dark energy models. In late time the viscous new HDE model approaches to \(\varLambda \mathit{CDM}\) model. We also discuss the thermodynamics and entropy of the model and find that it satisfies the second law of thermodynamics.  相似文献   

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

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

京公网安备 11010802026262号