首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Small tidal forces in the Earth–Moon system cause detectable changes in the orbit. Tidal energy dissipation causes secular rates in the lunar mean motion n, semimajor axis a, and eccentricity e. Terrestrial dissipation causes most of the tidal change in n and a, but lunar dissipation decreases eccentricity rate. Terrestrial tidal dissipation also slows the rotation of the Earth and increases obliquity. A tidal acceleration model is used for integration of the lunar orbit. Analysis of lunar laser ranging (LLR) data provides two or three terrestrial and two lunar dissipation parameters. Additional parameters come from geophysical knowledge of terrestrial tides. When those parameters are converted to secular rates for orbit elements, one obtains dn/dt = \(-25.97\pm 0.05 ''/\)cent\(^{2}\), da/dt = 38.30 ± 0.08 mm/year, and di/dt = ?0.5 ± 0.1 \(\upmu \)as/year. Solving for two terrestrial time delays and an extra de/dt from unspecified causes gives \(\sim \) \(3\times 10^{-12}\)/year for the latter; solving for three LLR tidal time delays without the extra de/dt gives a larger phase lag of the N2 tide so that total de/dt = \((1.50 \pm 0.10)\times 10^{-11}\)/year. For total dn/dt, there is \(\le \)1 % difference between geophysical models of average tidal dissipation in oceans and solid Earth and LLR results, and most of that difference comes from diurnal tides. The geophysical model predicts that tidal deceleration of Earth rotation is \(-1316 ''\)/cent\(^{2}\) or 87.5 s/cent\(^{2}\) for UT1-AT, a 2.395 ms/cent increase in the length of day, and an obliquity rate of 9 \(\upmu \)as/year. For evolution during past times of slow recession, the eccentricity rate can be negative.  相似文献   

2.
We have studied the spatial and temporal distribution of abundances of chemical elements in large “gradual” solar energetic-particle (SEP) events, and especially the source plasma temperatures, derived from those abundances, using measurements from the Wind and Solar TErrestrial RElations Observatory (STEREO) spacecraft, widely separated in solar longitude. A power-law relationship between abundance enhancements and mass-to-charge ratios [\(A/Q\)] of the ions can be used to determine \(Q\)-values and source plasma temperatures at remote spacecraft with instruments that were not designed for charge-state measurements. We search for possible source variations along the accelerating shock wave, finding one clear case where the accelerating shock wave appears to dispatch ions from \(3.2\pm 0.8~\mbox{MK}\) plasma toward one spacecraft and those from \(1.6\pm 0.2~\mbox{MK}\) plasma toward another, 116° away. The difference persists for three days and then fades away. Three other SEP events show less-extreme variation in source temperatures at different spacecraft, in one case observed over 222° in longitude. This initial study shows how the power-law relation between abundance enhancements and ion \(A/Q\)-values provides a new technique to determine \(Q\) and plasma temperatures in the seed population of SEP ions over a broad region of space using remote spacecraft with instruments that were not originally designed for measurements of ionization states.  相似文献   

3.
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.  相似文献   

4.
In this study, we present CCD UBV photometry of poorly studied open star clusters, Dolidze 36, NGC 6728, NGC 6800, NGC 7209, and Platais 1, located in the first and second Galactic quadrants. Observations were obtained with T100, the 1-m telescope of the TÜB?TAK National Observatory. Using photometric data, we determined several astrophysical parameters such as reddening, distance, metallicity and ages and from them, initial mass functions, integrated magnitudes and colours. We took into account the proper motions of the observed stars to calculate the membership probabilities. The colour excesses and metallicities were determined independently using two-colour diagrams. After obtaining the colour excesses of the clusters Dolidze 36, NGC 6728, NGC 6800, NGC 7209, and Platais 1 as \(0.19\pm0.06\), \(0.15\pm0.05\), \(0.32\pm0.05\), \(0.12\pm 0.04\), and \(0.43\pm0.06\) mag, respectively, the metallicities are found to be \(0.00\pm0.09\), \(0.02\pm0.11\), \(0.03\pm0.07\), \(0.01\pm0.08\), and \(0.01\pm0.08\) dex, respectively. Furthermore, using these parameters, distance moduli and age of the clusters were also calculated from colour-magnitude diagrams simultaneously using PARSEC theoretical models. The distances to the clusters Dolidze 36, NGC 6728, NGC 6800, NGC 7209, and Platais 1 are \(1050\pm90\), \(1610\pm190\), \(1210\pm150\), \(1060\pm90\), and \(1710\pm250\) pc, respectively, while corresponding ages are \(400\pm100\), \(750\pm150\), \(400\pm100\), \(600\pm100\), and \(175\pm50\) Myr, respectively. Our results are compatible with those found in previous studies. The mass function of each cluster is derived. The slopes of the mass functions of the open clusters range from 1.31 to 1.58, which are in agreement with Salpeter’s initial mass function. We also found integrated absolute magnitudes varying from ?4.08 to ?3.40 for the clusters.  相似文献   

5.
To investigate the \(M_\bullet -\sigma \) relation, we consider realistic elliptical galaxy profiles that are taken to follow a single power-law density profile given by \(\rho (r) = \rho _{0}(r/ r_{0})^{-\gamma }\) or the Nuker intensity profile. We calculate the density using Abel’s formula in the latter case by employing the derived stellar potential; in both cases. We derive the distribution function f(E) of the stars in the presence of the supermassive black hole (SMBH) at the center and hence compute the line-of-sight (LoS) velocity dispersion as a function of radius. For the typical range of values for masses of SMBH, we obtain \(M_{\bullet } \propto \sigma ^{p}\) for different profiles. An analytical relation \(p = (2\gamma + 6)/(2 + \gamma )\) is found which is in reasonable agreement with observations (for \(\gamma = 0.75{-}1.4\), \(p = 3.6{-}5.3\)). Assuming that a proportionality relation holds between the black hole mass and bulge mass, \(M_{\bullet } =f M_\mathrm{b}\), and applying this to several galaxies, we find the individual best fit values of p as a function of f; also by minimizing \(\chi ^{2}\), we find the best fit global p and f. For Nuker profiles, we find that \(p = 3.81 \pm 0.004\) and \(f = (1.23 \pm 0.09)\times 10^{-3}\) which are consistent with the observed ranges.  相似文献   

6.
We examine the average magnetic field magnitude (\(| \boldsymbol{B} | \equiv B\)) within magnetic clouds (MCs) observed by the Wind spacecraft from 1995 to July 2015 to understand the difference between this \(B\) and the ideal \(B\)-profiles expected from using the static, constant-\(\alpha\), force-free, cylindrically symmetric model for MCs of Lepping, Jones, and Burlaga (J. Geophys. Res. 95, 11957, 1990, denoted here as the LJB model). We classify all MCs according to an assigned quality, \(Q_{0}\) (\(= 1, 2, 3\), for excellent, good, and poor). There are a total of 209 MCs and 124 when only \(Q_{0} = 1\), 2 cases are considered. The average normalized field with respect to the closest approach (\(\mathit{CA}\)) is stressed, where we separate cases into four \(\mathit{CA}\) sets centered at 12.5 %, 37.5 %, 62.5 %, and 87.5 % of the average radius; the averaging is done on a percentage-duration basis to treat all cases the same. Normalized \(B\) means that before averaging, the \(B\) for each MC at each point is divided by the LJB model-estimated \(B\) for the MC axis, \(B_{0}\). The actual averages for the 209 and 124 MC sets are compared to the LJB model, after an adjustment for MC expansion (e.g. Lepping et al. in Ann. Geophys. 26, 1919, 2008). This provides four separate difference-relationships, each fitted with a quadratic (Quad) curve of very small \(\sigma\). Interpreting these Quad formulae should provide a comprehensive view of the variation in normalized \(B\) throughout the average MC, where we expect external front and rear compression to be part of its explanation. These formulae are also being considered for modifying the LJB model. This modification will be used in a scheme for forecasting the timing and magnitude of magnetic storms caused by MCs. Extensive testing of the Quad formulae shows that the formulae are quite useful in correcting individual MC \(B\)-profiles, especially for the first \({\approx\,}1/3\) of these MCs. However, the use of this type of \(B\) correction constitutes a (slight) violation of the force-free assumption used in the original LJB MC model.  相似文献   

7.
We analyze the families of central configurations of the spatial 5-body problem with four masses equal to 1 when the fifth mass m varies from 0 to \(+\infty \). In particular we continue numerically, taking m as a parameter, the central configurations (which all are symmetric) of the restricted spatial (\(4+1\))-body problem with four equal masses and \(m=0\) to the spatial 5-body problem with equal masses (i.e. \(m=1\)), and viceversa we continue the symmetric central configurations of the spatial 5-body problem with five equal masses to the restricted (\(4+1\))-body problem with four equal masses. Additionally we continue numerically the symmetric central configurations of the spatial 5-body problem with four equal masses starting with \(m=1\) and ending in \(m=+\infty \), improving the results of Alvarez-Ramírez et al. (Discrete Contin Dyn Syst Ser S 1: 505–518, 2008). We find four bifurcation values of m where the number of central configuration changes. We note that the central configurations of all continued families varying m from 0 to \(+\infty \) are symmetric.  相似文献   

8.
We study the solar-cycle variation of subsurface flows from the surface to a depth of 16 Mm. We have used ring-diagram analysis to analyze Dopplergrams obtained with the Michelson Doppler Imager (MDI) Dynamics Program, the Global Oscillation Network Group (GONG), and the Helioseismic and Magnetic Imager (HMI) instrument. We combined the zonal and meridional flows from the three data sources and scaled the flows derived from MDI and GONG to match those from HMI observations. In this way, we derived their temporal variation in a consistent manner for Solar Cycles 23 and 24. We have corrected the measured flows for systematic effects that vary with disk positions. Using time-depth slices of the corrected subsurface flows, we derived the amplitudes and times of the extrema of the fast and slow zonal and meridional flows during Cycles 23 and 24 at every depth and latitude. We find an average difference between maximum and minimum amplitudes of \(8.6 \pm0.4~\mbox{m}\,\mbox{s}^{-1}\) for the zonal flows and \(7.9 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the meridional flows associated with Cycle 24 averaged over a depth range from 2 to 12 Mm. The corresponding values derived from GONG data alone are \(10.5 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the zonal and \(10.8 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the meridional flow. For Cycle 24, the flow patterns are precursors of the magnetic activity. The timing difference between the occurrence of the flow pattern and the magnetic one increases almost linearly with increasing latitude. For example, the fast zonal and meridional flow appear \(2.1 \pm 0.6\) years and \(2.5\pm 0.6\) years, respectively, before the magnetic pattern at \(30^{\circ}\) latitude in the northern hemisphere, while in the southern hemisphere, the differences are \(3.2 \pm 1.2\) years and \(2.6 \pm 0.6\) years. The flow patterns of Cycle 25 are present and have reached \(30^{\circ}\) latitude. The amplitude differences of Cycle 25 are about 22% smaller than those of Cycle 24, but are comparable to those of Cycle 23. Moreover, polynomial fits of meridional flows suggest that equatorward meridional flows (counter-cells) might exist at about \(80^{\circ}\) latitude except during the declining phase of the solar cycle.  相似文献   

9.
A statistical study of the chromospheric ribbon evolution in H\(\alpha\) two-ribbon flares was performed. The data set consists of 50 confined (62%) and eruptive (38%) flares that occurred from June 2000 to June 2015. The flares were selected homogeneously over the H\(\alpha\) and Geostationary Operational Environmental Satellite (GOES) classes, with an emphasis on including powerful confined flares and weak eruptive flares. H\(\alpha\) filtergrams from the Kanzelhöhe Observatory in combination with Michelson Doppler Imager (MDI) and Helioseismic and Magnetic Imager (HMI) magnetograms were used to derive the ribbon separation, the ribbon-separation velocity, the magnetic-field strength, and the reconnection electric field. We find that eruptive flares reveal statistically larger ribbon separation and higher ribbon-separation velocities than confined flares. In addition, the ribbon separation of eruptive flares correlates with the GOES SXR flux, whereas no clear dependence was found for confined flares. The maximum ribbon-separation velocity is not correlated with the GOES flux, but eruptive flares reveal on average a higher ribbon-separation velocity (by ≈?10 km?s?1). The local reconnection electric field of confined (\(cc=0.50 \pm0.02\)) and eruptive (\(cc=0.77 \pm0.03\)) flares correlates with the GOES flux, indicating that more powerful flares involve stronger reconnection electric fields. In addition, eruptive flares with higher electric-field strengths tend to be accompanied by faster coronal mass ejections.  相似文献   

10.
In this work we consider the Kepler problem with linear drag, and prove the existence of a continuous vector-valued first integral, obtained taking the limit as \(t\rightarrow +\infty \) of the Runge–Lenz vector. The norm of this first integral can be interpreted as an asymptotic eccentricity \(e_{\infty }\) with \(0\le e_{\infty } \le 1\). The orbits satisfying \(e_{\infty } <1\) approach the singularity by an elliptic spiral and the corresponding solutions \(x(t)=r(t)e^{i\theta (t)}\) have a norm r(t) that goes to zero like a negative exponential and an argument \(\theta (t)\) that goes to infinity like a positive exponential. In particular, the difference between consecutive times of passage through the pericenter, say \(T_{n+1} -T_n\), goes to zero as \(\frac{1}{n}\).  相似文献   

11.
We analyzed the BVR photometry of comet C/2014 S2 obtained between March and June 2016, in observatories installed in Europe and the United States. Using the Lomb–Scargle periodogram, we found that the most probable periodicity deduced from the V-band magnitudes is 2.70 days, suggesting that it is the period of rotation of the nucleus of this comet is \(2.70 \pm 0.07\) days or \(68 \pm 2\) h, with a peak-to-peak light curve amplitude of \(0.4 \pm 0.1\) magnitudes. We verify that the absolute magnitude \(H_0\) and the activity index n differ from each other when they are calculated from the visual or CCD magnitudes. Considering the absolute magnitude \(H_{v0}=\) 6.0, obtained from visual magnitudes, we estimate that the lower limit of nuclear radius is 1.3 km. Analyzing the variation of magnitude R with the photometric aperture, we suggest that the coma of this object was in steady-state within the time limits of our observational interval. The coma had a mean color index B–V\(=0.79\pm 0.22\), which is typical of active comets. Additionally, we have shown that the use of a variable photometric aperture, linked to geocentric distance, is probably unnecessary for the comet PANSTARRS .  相似文献   

12.
On 27 June 2012, an eruptive solar prominence was observed in the extreme ultraviolet (EUV) and radio wavebands. At the Aalto University Metsähovi Radio Observatory (MRO) it was observed at 37 GHz. It was the first time that the MRO followed a radio prominence with dense sampling in the millimetre wavelengths. This prompted us to study the connection of the 37 GHz event with other wavelength domains. At 37 GHz, the prominence was tracked to a height of around \(1.6~\mathrm{R}_{\odot}\), at which the loop structure collapsed. The average velocity of the radio prominence was \(55 \pm 6~\mbox{km}\,\mbox{s}^{-1}\). The brightness temperature of the prominence varied between \(800 \pm 100\) K and \(3200 \pm 100\) K. We compared our data with the Solar Dynamic Observatory (SDO)/Atmospheric Imaging Assembly (AIA) instrument’s 304 Å EUV data, and found that the prominence behaves very similarly in both wavelengths. The EUV data also reveal flaring activity nearby the prominence. We present a scenario in which this flare works as a trigger that causes the prominence to move from a stable stage to an acceleration stage.  相似文献   

13.
Data of geomagnetic indices (aa, Kp, Ap, and Dst) recorded near 1 AU over the period 1967–2016, have been studied based on the asymmetry between the interplanetary magnetic field (IMF) directions above and below of the heliospheric current sheet (HCS). Our results led to the following conclusions: (i) Throughout the considered period, 31 random years (62%) showed apparent asymmetries between Toward (\(\mathbf{T}\)) and Away (\(\mathbf{A}\)) polarity days and 19 years (38%) exhibited nearly a symmetrical behavior. The days of \(\mathbf{A}\) polarity predominated over the \(\mathbf{T}\) polarity days by 4.3% during the positive magnetic polarity epoch (1991–1999). While the days of \(\mathbf{T}\) polarity exceeded the days of \(\mathbf{A}\) polarity by 5.8% during the negative magnetic polarity epoch (2001–2012). (ii) Considerable yearly North–South (N–S) asymmetries of geomagnetic indices observed throughout the considered period. (iii) The largest toward dominant peaks for \(aa\) and \(Ap\) indices occurred in 1995 near to minimum of solar activity. Moreover, the most substantial away dominant peaks for \(aa\) and \(Ap\) indices occurred in 2003 (during the descending phase of the solar cycle 23) and in 1991 (near the maximum of solar activity cycle) respectively. (iv) The N–S asymmetry of \(Kp\) index indicated a most significant away dominant peak occurred in 2003. (v) Four of the away dominant peaks of Dst index occurred at the maxima of solar activity in the years 1980, 1990, 2000, and 2013. The largest toward dominant peak occurred in 1991 (at the reversal of IMF polarity). (vi) The geomagnetic indices (aa, Ap, and \(Kp\)) all have northern dominance during positive magnetic polarity epoch (1971–1979), while the asymmetries shifts to the southern solar hemisphere during negative magnetic polarity epoch (2001–2012).  相似文献   

14.
The outer gravitational potential V of the level ellipsoid of revolution T is uniquely determined by two quantities: the eccentricity \(\varepsilon \) of the ellipsoid and Clairaut parameter q, proportional to the angular velocity of rotation squared and inversely proportional to the mean density of the ellipsoid. Quantities \(\varepsilon \) and q are independent, though they lie in a rather strict two-dimensional domain. It follows that Stokes coefficients \(I_n\) of Laplace series representing the outer potential of T are uniquely determined by \(\varepsilon \) and q. In this paper, we have found explicit expressions for Stokes coefficients via \(\varepsilon \) and q, as well as their asymptotics when \(n\rightarrow \infty \). If T does not coincide with a Maclaurin ellipsoid, then \(|I_n|\sim B\varepsilon ^n/n\) with a certain constant B. Let us compare this asymptotics with one of \(I_n\) for ellipsoids constrained by the only condition of increasing (even nonstrict) of oblateness from the centre to the periphery: \(|I_n|\sim \bar{B}\varepsilon ^n/(n^2)\). Hence, level ellipsoids with ellipsoidal equidensites do not exist. The only exception represents Maclaurin ellipsoids. It should be recalled that we confine ourselves by ellipsoids of revolution.  相似文献   

15.
Based on energetic particle observations made at \({\approx}\,1\) AU, we present a catalogue of 46 wide-longitude (\({>}\,45^{\circ}\)) solar energetic particle (SEP) events detected at multiple locations during 2009?–?2016. The particle kinetic energies of interest were chosen as \({>}\,55\) MeV for protons and 0.18?–?0.31 MeV for electrons. We make use of proton data from the Solar and Heliospheric Observatory/Energetic and Relativistic Nuclei and Electron Experiment (SOHO/ERNE) and the Solar Terrestrial Relations Observatory/High Energy Telescopes (STEREO/HET), together with electron data from the Advanced Composition Explorer/Electron, Proton, and Alpha Monitor (ACE/EPAM) and the STEREO/Solar Electron and Proton Telescopes (SEPT). We consider soft X-ray data from the Geostationary Operational Environmental Satellites (GOES) and coronal mass ejection (CME) observations made with the SOHO/Large Angle and Spectrometric Coronagraph (LASCO) and STEREO/Coronagraphs 1 and 2 (COR1, COR2) to establish the probable associations between SEP events and the related solar phenomena. Event onset times and peak intensities are determined; velocity dispersion analysis (VDA) and time-shifting analysis (TSA) are performed for protons; TSA is performed for electrons. In our event sample, there is a tendency for the highest peak intensities to occur when the observer is magnetically connected to solar regions west of the flare. Our estimates for the mean event width, derived as the standard deviation of a Gaussian curve modelling the SEP intensities (protons \({\approx}\,44^{\circ}\), electrons \({\approx}\,50^{\circ}\)), largely agree with previous results for lower-energy SEPs. SEP release times with respect to event flares, as well as the event rise times, show no simple dependence on the observer’s connection angle, suggesting that the source region extent and dominant particle acceleration and transport mechanisms are important in defining these characteristics of an event. There is no marked difference between the speed distributions of the CMEs related to wide events and the CMEs related to all near-Earth SEP events of similar energy range from the same time period.  相似文献   

16.
We will show that the period T of a closed orbit of the planar circular restricted three body problem (viewed on rotating coordinates) depends on the region it encloses. Roughly speaking, we show that, \(2 T=k\pi +\int _\Omega g\) where k is an integer, \(\Omega \) is the region enclosed by the periodic orbit and \(g:{\mathbb {R}}^2\rightarrow {\mathbb {R}}\) is a function that only depends on the constant C known as the Jacobian constant; it does not depend on \(\Omega \). This theorem has a Keplerian flavor in the sense that it relates the period with the space “swept” by the orbit. As an application we prove that there is a neighborhood around \(L_4\) such that every periodic solution contained in this neighborhood must move clockwise. The same result holds true for \(L_5\).  相似文献   

17.
In a two-component jet model, the emissions are the sum of the core and extended emissions: \(S^{\mathrm{ob}}=S_{\mathrm{core}}^{\mathrm{ob}}+S_{\mathrm{ext}}^{\mathrm{ob}}\), with the core emissions, \(S_{\mathrm{core}}^{\mathrm{ob}}= f S_{\mathrm{ext}}^{\mathrm{ob}}\delta ^{q}\) being a function of the Doppler factor \(\delta \), the extended emission \(S_{\mathrm{ext}}^{\mathrm{ob}}\), the jet type dependent factor q, and the ratio of the core to the extended emissions in the comoving frame, f. The f is an unobservable but important parameter. Following our previous work, we collect 65 blazars with available Doppler factor \(\delta \), superluminal velocity \(\beta _{\mathrm{app}}\), and core-dominance parameter, R, and calculated the ratio, f, and performed statistical analyses. We found that the ratio, f, in BL Lacs is on average larger than that in FSRQs. We suggest that the difference of the ratio f between FSRQs and BL Lacs is one of the possible reasons that cause the difference of other observed properties between them. We also find some significant correlations between \(\log f\) and other parameters, including intrinsic (de-beamed) peak frequency, \(\log \nu _{\mathrm{p}}^{\mathrm{in}}\), intrinsic polarization, \(\log P^{\mathrm{in}}\), and core-dominance parameter, \(\log R\), for the whole sample. In addition, we show that the ratio, f, can be estimated by R.  相似文献   

18.
We aim to probe the dynamic structure of the extended Solar neighborhood by calculating the radial metallicity gradients from orbit properties, which are obtained for axisymmetric and non-axisymmetric potential models, of red clump (RC) stars selected from the RAdial Velocity Experiment’s Fourth Data Release. Distances are obtained by assuming a single absolute magnitude value in near-infrared, i.e. \(M_{Ks}=-1.54\pm0.04\) mag, for each RC star. Stellar orbit parameters are calculated by using the potential functions: (i) for the MWPotential2014 potential, (ii) for the same potential with perturbation functions of the Galactic bar and transient spiral arms. The stellar age is calculated with a method based on Bayesian statistics. The radial metallicity gradients are evaluated based on the maximum vertical distance (\(z_{max}\)) from the Galactic plane and the planar eccentricity (\(e_{p}\)) of RC stars for both of the potential models. The largest radial metallicity gradient in the \(0< z_{max} \leq0.5\) kpc distance interval is \(-0.065\pm0.005~\mbox{dex}\,\mbox{kpc}^{-1}\) for a subsample with \(e_{p}\leq0.1\), while the lowest value is \(-0.014\pm0.006~\mbox{dex}\,\mbox{kpc}^{-1}\) for the subsample with \(e_{p}\leq0.5\). We find that at \(z_{max}>1\) kpc, the radial metallicity gradients have zero or positive values and they do not depend on \(e_{p}\) subsamples. There is a large radial metallicity gradient for thin disc, but no radial gradient found for thick disc. Moreover, the largest radial metallicity gradients are obtained where the outer Lindblad resonance region is effective. We claim that this apparent change in radial metallicity gradients in the thin disc is a result of orbital perturbation originating from the existing resonance regions.  相似文献   

19.
We carried out high-precision photometric observations of three eclipsing ultrashort-period contact binaries (USPCBs). Theoretical models were fitted to the light curves by means of the Wilson-Devinney code. The solutions suggest that the three targets have evolved to a contact phase. The photometric results are as follows: (a) 1SWASP?J030749.87?365201.7, \(q=0.439\pm0.003\), \(f=0.0\pm3.6\%\); (b) 1SWASP?J213252.93?441822.6, \(q=0.560\pm0.003\), \(f=14.2\pm1.9\%\); (c) 1SWASP?J200059.78+054408.9, \(q=0.436\pm0.008\), \(f=58.4\pm1.8\%\). The light curves show O’Connell effects, which can be modeled by the assumed cool spots. The cool spots models are strongly supported by the night-to-night variations in the \(I\)-band light curves of 1SWASP?J030749.87?365201.7. For a comparative study, we collected the whole set of 28 well-studied USPCBs with \(P < 0.24\) day. Thus, we found that most of them (17 of 28) are in shallow contact (i.e. fill-out factors \(f<20\%\)). Only four USPCBs have deep fill-out factors (i.e. \(f>50\%\)). Generally, contact binaries with deep fill-out factors are going to merge, but it is believed that USPCBs have just evolved to a contact phase. Hence, the deep USPCB 1SWASP?J200059.78+054408.9 seems to be a contradiction, making it very interesting. Particularly, 1SWASP?J030749.87?365201.7 is a zero contact binary in thermal equilibrium, implying that it should be a turn-off sample as predicted by the thermal relaxation oscillation (TRO) theory.  相似文献   

20.
New optimal, contractivity-preserving (CP), explicit, d-derivative, k-step Hermite–Obrechkoff series methods of order p up to \(p=20\), denoted by CP HO(d, k, p), with nonnegative coefficients are constructed. These methods are used to solve nonstiff first-order initial value problems \(y'=f(t,y)\), \(y(t_0)=y_0\). The upper bound \(p_u\) of order p of HO(d, k, p) can reach, approximately, as high as 2.4 times the number of derivatives d. The stability regions of HO(d, k, p) have generally a good shape and grow with decreasing \(p-d\). We, first, note that three selected CP HO methods: 4-derivative 7-step HO of order 13, denoted by HO(4, 7, 13), 5-derivative 6-step HO of order 13, denoted by HO(5, 6, 13), and 9-derivative 2-step HO of order 13, denoted by CMDAHO(13) compare favorably with Adams–Cowell of order 13, denoted by AC(13), in solving standard N-body problems over an interval of 1000 periods on the basis of the relative error of energy as a function of the CPU time. Next, the three HO methods compare positively with AC(13) in solving standard N-body problems on the basis of the growth of relative positional error and relative energy error over 10, 000 periods of integration. Finally, these three methods compare also well with P-stable methods of Cash and Franco et al. on some quasi periodic, second-order linear and nonlinear problems. The coefficients of selected HO methods are listed in the appendix.  相似文献   

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

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