首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 60 毫秒
1.
N. Hiotelis   《New Astronomy》2002,7(8):531-539
We present density profiles, that are solutions of the spherical Jeans equation, derived under the following two assumptions: (i) the coarse grained phase-density follows a power-law of radius, ρ/σ3r, and (ii) the velocity anisotropy parameter is given by the relation βa(r)=β1+2β2 (r/r*)/[1+(r/r*)2] where β1, β2 are parameters and r* equals twice the virial radius, rvir, of the system. These assumptions are well motivated by the results of N-body simulations. Density profiles have increasing logarithmic slopes γ, defined by γ=−d ln ρ/d ln r. The values of γ at r=10−2.5rvir, a distance where the systems could be resolved by large N-body simulations, lie in the range 1.0–1.6. These inner values of γ increase for increasing β1 and for increasing concentration of the system. On the other hand, slopes at r=rvir lie in the range 2.42–3.82. A model density profile that fits well the results at radial distances between 10−3rvir and rvir and connects kinematic and structural characteristics of spherical systems is described.  相似文献   

2.
We have developed a parallel Particle–Particle, Particle–Mesh (P3M) simulation code for the Cray T3E parallel supercomputer that is well suited to studying the time evolution of systems of particles interacting via gravity and gas forces in cosmological contexts. The parallel code is based upon the public-domain serial Adaptive P3M-SPH (http://coho.astro.uwo.ca/pub/hydra/hydra.html) code of Couchman et al. (1995)[ApJ, 452, 797]. The algorithm resolves gravitational forces into a long-range component computed by discretizing the mass distribution and solving Poisson's equation on a grid using an FFT convolution method, and a short-range component computed by direct force summation for sufficiently close particle pairs. The code consists primarily of a particle–particle computation parallelized by domain decomposition over blocks of neighbour-cells, a more regular mesh calculation distributed in planes along one dimension, and several transformations between the two distributions. The load balancing of the P3M code is static, since this greatly aids the ongoing implementation of parallel adaptive refinements of the particle and mesh systems. Great care was taken throughout to make optimal use of the available memory, so that a version of the current implementation has been used to simulate systems of up to 109 particles with a 10243 mesh for the long-range force computation. These are the largest Cosmological N-body simulations of which we are aware. We discuss these memory optimizations as well as those motivated by computational performance. Performance results are very encouraging, and, even without refinements, the code has been used effectively for simulations in which the particle distribution becomes highly clustered as well as for other non-uniform systems of astrophysical interest.  相似文献   

3.
We compute the big bang nucleosynthesis limit on the number of light neutrino degrees of freedom in a model-independent likelihood analysis based on the abundances of 4He and 7Li. We use the two-dimensional likelihood functions to simultaneously constrain the baryon-to-photon ratio and the number of light neutrinos for a range of 4He abundances Yp = 0.225–0.250, as well as a range in primordial 7Li abundances from (1.6 to 4.1) ×10−10. For (7Li/H)p = 1.6 × 10−10, as can be inferred from the 7Li data from Population II halo stars, the upper limit to Nν based on the current best estimate of the primordial 4He abundance of Yp = 0.238 is Nν < 4.3 and varies from Nν < 3.3 (at 95% C.L.) when Yp = 0.225 to Nν < 5.3 when Yp = 0.250. If 7Li is depleted in these stars the upper limit to Nν is relaxed. Taking (7Li/H)p = 4.1 × 10−10, the limit varies from Nν < 3.9 when Yp = 0.225 to Nν 6 when Yp = 0.250. We also consider the consequences on the upper limit to Nν if recent observations of deuterium in high-redshift quasar absorption-line systems are confirmed.  相似文献   

4.
Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of 13CO, C18O, HCO+ and N2H+ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from 13CO to HCO+, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 103M to 5.1 × 102M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of 13CO to C18O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N2H+ molecules is 3.5 ± 2.5 × 10−10, consistent with the value 1.0 − 5.0 × 10−10 for dark cloud cores and the value 1.2 − 12.8 × 10−10 for massive cores. The abundance of HCO+ molecules is 0.9 ± 0.5 × 10−9, close to the value 1.6 − 2.4 × 10−9 for massive cores. An increase of HCO+ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.  相似文献   

5.
We present sixth- and eighth-order Hermite integrators for astrophysical N-body simulations, which use the derivatives of accelerations up to second-order (snap) and third-order (crackle). These schemes do not require previous values for the corrector, and require only one previous value to construct the predictor. Thus, they are fairly easy to implement. The additional cost of the calculation of the higher-order derivatives is not very high. Even for the eighth-order scheme, the number of floating-point operations for force calculation is only about two times larger than that for traditional fourth-order Hermite scheme. The sixth-order scheme is better than the traditional fourth-order scheme for most cases. When the required accuracy is very high, the eighth-order one is the best. These high-order schemes have several practical advantages. For example, they allow a larger number of particles to be integrated in parallel than the fourth-order scheme does, resulting in higher execution efficiency in both general-purpose parallel computers and GRAPE systems.  相似文献   

6.
《Astroparticle Physics》2005,22(5-6):339-353
Data taken during 1997 with the AMANDA-B10 detector are searched for a diffuse flux of neutrinos of all flavors with energies above 1016 eV. At these energies the Earth is opaque to neutrinos, and thus neutrino induced events are concentrated at the horizon. The background are large muon bundles from down-going atmospheric air shower events. No excess events above the background expectation are observed and a neutrino flux following E−2, with an equal mix of all flavors, is limited to E2Φ(1015 eV < E < 3 × 1018 eV)  0.99 × 10−6 GeV cm−2 s−1 sr−1 at 90% confidence level. This is the most restrictive experimental bound placed by any neutrino detector at these energies. Bounds to specific extraterrestrial neutrino flux predictions are also presented.  相似文献   

7.
Venera 9, 10 measurements of the nightside ionospheric profile and the night airglow were used for investigating ionosphere formation processes. The upper ionospheric layer may be formed by HeI 584 Å radiation; the lower layer by meteorite ionization. Upper limits on the electron energy flux, <4 × 108eV cm−2 s−1, the helium ion flux <107 cm−2 s−1, the nitric oxide mixing ratio, <1.5 × 10−4 and the atomic sulphur mixing ratio, <10−6, are deduced for ionospheric altitudes.  相似文献   

8.
Radio noise observations at frequencies of 0·700 Mc and 2·200 Mc were made at altitudes between 3000 and 11,000 km from a Blue Scout Jr. high-altitude rocket probe on 30 July 1963. A steady background flux of (7·5−3+6) × 10−19 W m−2)(c/s)−1 at 0·700 Mc and (1·8+1.0−0.5 × 10−19 W m−2 (c/s)−1 at 2·200 Mc was observed. Assuming a galactic origin of the observed fluxes at both frequencies, the averaged sky brightnesses are b(0·700 Mc) = (6−3+5) × 10−20 W m−2 (c/s)−1 sr−1b(2·200 Mc) = (1.4+1.0−0.5 × 10−20 W m−2 (c/s)−1 sr−1 The observed brightness at 2·200 Mc is in reasonable agreement with the results of other observers. The apparent brightness at 0·700 Mc is, however, greater than was expected from previous observations. An alternative source of the 0·700 Mc flux in the terrestrial exosphere, as well as characteristics of several noise bursts observed during the flight, is briefly discussed.  相似文献   

9.
A numerical analysis of cyclotron instabilities is carried out by computing the dispersion relation for a three component cold plasma-beam system. Rates of growth and damping for various values of the stream density are calculated from the dispersion relation. The rates of growth and damping increase monotonically as the number density of the proton stream increases. It is found that the frequencies at the rates of maximum growth and the damping decrease slightly to lower frequencies and a sharp peak at these frequencies becomes blunt. The minimum e-folding times of an ion cyclotron wave for (a) σs = 10−4, σi = 10−2 and (b) σs = 10−1, σi = 10−2 are about 3·84 and 0·16 sec respectively in the vicinity of the equatorial plane at 6 Re, where σs and σi are the ratios of the beam density Ns and the helium ion (H6+) density Ni to the total positive ions in the plasma-beam system.  相似文献   

10.
The MSX infrared dark cloud G79.2+0.38 has been observed over a 11′×′ region simultaneously in the J=1-0 rotational transition lines of the 12CO and its isotopic molecules 13CO and 18CO. The dense molecular cores defined by the C18O line are found to be associated with the two high-extinction patches shown in the MSX A-band image. The two dense cores have the column density N (H2) (5 – 12) × 1022 cm−2 and the mean number density n (3 ± 1) × 104 cm−3. Their sizes are 1.7 and 1.2 pc in 13CO(1-0) line, 1.2 and 0.6 pc in C18O(1-0) line, respectively. The masses of these cloud cores are estimated to be in the range from 2 × 102 to 2 × 103 M. The profile of radial mean density of the cloud core can be described by the exponential function ¯n(p) p−0.34±0.02. Compared with the cases of typical optical dark clouds, the abundances of the CO isotopic molecules 13CO and C18O in this MSX infrared dark cloud appear to be depleted by a factor of 4–11, but at present there is no evidence for any obvious variation of the relative abundance ratio X13/18 between 13CO and C18O with the column density.  相似文献   

11.
We used the 13.7 m millimeter-wave telescope of the Qinghai Station of Purple Mountain Observatory at Delingha to carry out simultaneous mapping observations in the 12CO(J = 1 − 0), 13CO(J = 1 − 0) and CO18 (J = 1 − 0) lines towards 11 isolated star-forming molecular cores selected from the sample of Spitzer's c2d program. For all three lines, these sources were all observed to the positions of half the peak intensity of the 13CO(J = 1 − 0) line. The volume density, local thermodynamic equilibrium mass MLTE and virial mass MVIR were derived for each of the molecular cores. The obtained ratio of MVIR to MLTE is 0.85 ± 0.40 for the 13CO(J = 1 − 0) cores and 0.77 ± 0.35 for the CO18 (J = 1 − 0) cores. The density profiles of the cloud cores were also calculated.  相似文献   

12.
The orbit of Intercosmos 13 rocket (1975-22B) has been determined at 103 epochs between 30 April 1975 and 10 April 1980 from almost 7000 observations. One hundred and three values of inclination have been determined and corrections incoporated for the effects due to zonal harmonic, lunisolar and tesseral harmonic perturbations, precession, and solid Earth tides. The modified data have been analysed to yield values of the atmospheric rotation rate, Λ rev day−1, viz. Λ = 0.94 ± 0.10 at an average height of 322 ± 6 km and Λ = 1.27 ± 0.02 at 288 km. Analysis of the inclination near 14th-order resonance has indicated lumped harmonic values 109 1.01.4 = − 76.13 ± 12.47, 109 1,014 = − 29.89 ± 32.64, 109 −1.214 = − 63.11 ± 15.44 109 −1.214 = − 32.52 ± 26.96, for inclination 82.952°.  相似文献   

13.
First generation stars are the oldest stars that were formed in post-big bang, primitive gas, containing no elements heavier than carbon, with ages greater than 14 Gyr and having undergone no evolution so far. Observations over a long time have confirmed that, up to now, no stars with zero metallicity ([Fe/H]) or with [Fe/H]≤ −6 have been found in the Galaxy. To explain this absence, we shall make a theoretical calculation of the probability of detecting first generation stars using Tsuiimoto et al.'s model of chemical evolution of the galactic halo and assuming an initial mass function of the Miller-Scalo form. We use all the observational data on the halo stars to constrain the parameters. Our result is that, if the mass of the cloud that formed the stars is 106–107 M, then the probability of detecting first generation stars is 6.14×10−4–6.14×10−5.  相似文献   

14.
We discuss the possibility of accurately estimating the source number density of ultra-high-energy cosmic rays (UHECRs) using small-scale anisotropy in their arrival distribution. The arrival distribution has information on their source and source distribution. We calculate the propagation of UHE protons in a structured extragalactic magnetic field (EGMF) and simulate their arrival distribution at the Earth using our previously developed method. The source number density that can best reproduce observational results by Akeno Giant Air Shower Array is estimated at about 10−5 Mpc−3 in a simple source model. Despite having large uncertainties of about one order of magnitude, due to small number of observed events in current status, we find that more detection of UHECRs in the Auger era can sufficiently decrease this so that the source number density can be more robustly estimated. Two hundred event observation above 4 × 1019 eV in a hemisphere can discriminate between 10−5 and 10−6 Mpc−3. Number of events to discriminate between 10−4 and 10−5 Mpc−3 is dependent on EGMF strength. We also discuss the same in another source model in this paper.  相似文献   

15.
We calculated the expected neutrino signal in Borexino from a typical Type II supernova at a distance of 10 kpc. A burst of around 110 events would appear in Borexino within a time interval of about 10 s. Most of these events would come from the reaction channel , while about 30 events would be induced by the interaction of the supernova neutrino flux on 12C in the liquid scintillator. Borexino can clearly distinguish between the neutral-current excitations 12C(ν,ν)12C* (15.11 MeV) and the charged-current reactions 12C(νe,e)12N and , via their distinctive event signatures. The ratio of the charged-current to neutral-current neutrino event rates and their time profiles with respect to each other can provide a handle on supernova and non-standard neutrino physics (mass and flavor oscillations).  相似文献   

16.
Recent rocket observations of the N2 V-K (Vegard-Kaplan) system in the aurora have been reinterpreted using an atmospheric model based on mass spectrometer measurements in an aurora of similar intensity at the same time of year. In contrast to the original interpretation, we find that population by cascade from the C3Πu and B3Πg states in the A3Σu+v=0,1 levels, as calculated using recently measured electron excitation cross sections, accurately accounts for the observed relative emission rates (IV-K/12PG0.0). In addition there is no need to change the production rate of A 3 Σ u+ molecules relative to that of C3Πuv=0 as a function of altitude in order to fit the profile of the deactivation probability to the atmospheric model. Quenching of A 3 Σ u+ molecules at high altitudes is dominated by atomic oxygen. The rate constants for the v=0 and v=1 levels are 8 × 10−11 cm3 sec−1 and 1.7 × 10−10 cm3 sec−1 respectively, as determined using the model atmosphere mentioned above. Recent observations with a helium cooled mass spectrometer suggest that conventional mass spectrometer measurements tend to underestimate the atomic oxygen relative concentration. The rate coefficients may therefore be too large by as much as a factor of 3. Below 130 Km we find that it is possible to account for the deactivation in bright auroras by invoking large nitric oxide concentrations, similar to those recently observed mass spectrometrically and using a rate constant of 8 × 10−11 cm3 sec−1 for both the v=1 levels. This rate constant is very nearly the same as that measured in the laboratory (7 × 10−11 cm3 sec−1). Molecular oxygen appears not to play a significant role in deactivating the lower A 3 Σ u+ levels.  相似文献   

17.
The evolution of the cosmic ray primary composition in the energy range 106–107 GeV (i.e. the “knee” region) is studied by means of the e.m. and muon data of the Extensive Air Shower EAS-TOP array (Campo Imperatore, National Gran Sasso Laboratories). The measurement is performed through: (a) the correlated muon number (Nμ) and shower size (Ne) spectra, and (b) the evolution of the average muon numbers and their distributions as a function of the shower size. From analysis (a) the dominance of helium primaries at the knee, and therefore the possibility that the knee itself is due to a break in their energy spectrum (at EkHe=(3.5±0.3)×106 GeV) are deduced. Concerning analysis (b), the measurement accuracies allow the classification in terms of three mass groups: light (p,He), intermediate (CNO), and heavy (Fe). At primary energies E0≈106 GeV the results are consistent with the extrapolations of the data from direct experiments. In the knee region the obtained evolution of the energy spectra leads to: (i) an average steep spectrum of the light mass group (γp,He>3.1), (ii) a spectrum of the intermediate mass group harder than the one of the light component (γCNO2.75, possibly bending at EkCNO≈(6–7)×106 GeV), (iii) a constant slope for the spectrum of the heavy primaries (γFe2.3–2.7) consistent with the direct measurements. In the investigated energy range, the average primary mass increases from lnA=1.6–1.9 at E01.5×106 GeV to lnA=2.8–3.1 at E01.5×107 GeV. The result supports the standard acceleration and propagation models of galactic cosmic rays that predict rigidity dependent cut-offs for the primary spectra of the different nuclei. The uncertainties connected to the hadronic interaction model (QGSJET in CORSIKA) used for the interpretation are discussed.  相似文献   

18.
A plasmon model for the optical-radio knots of M87 is discussed with reference to their observed characteristics and physical conditions. The spectral evolution is investigated under the joint dissipation through expansion, Fermi acceleration and synchrotron radiation. A Fermi acceleration coefficient of 10−9 s−1 was obtained, and it was also shown that a spectral break should occur at 5 × 1014Hz and that this frequency should vary very slowly in time, in agreement with observations.  相似文献   

19.
The diffused gamma halo around our Galaxy recently discovered by EGRET could be produced by annihilations of heavy relic neutrinos N (of fourth generation), whose mass is within a narrow range (MZ/2<mN<MZ). Neutrino annihilation in the halo may lead to either ultrarelativistic electron pairs whose Inverse Compton Scattering on infrared and optical galactic photons could be the source of observed GeV gamma rays, or prompt 100 MeV–1 GeV photons (due to neutral pion secondaries) born by reactions. The consequent gamma flux (10−7–10−6 cm−2 s−1 sr−1) is well comparable to the EGRET observed one, and it is also compatible with the narrow window of neutrino mass 45 GeV <mN<50 GeV, recently required to explain the underground DAMA signals.The presence of heavy neutrinos of fourth generation do not contribute much to solve the dark matter problem of the Universe, but may be easily detectable by outcoming LEP II data.  相似文献   

20.
To study the climatological role of ozone in the Precambrian atmosphere and the consequences of its reduction for the ultraviolet environment of the early biosphere, a coupled one-dimensional radiative-convective and photochemical model has been developed. Oxygen levels between 10−5 and 1 time the present atmospheric level (PAL) are considered. It is shown that when the ice-albedo feedback is taken into account, relatively important temperature decreases are associated with the ozone changes linked to the progressive decrease of the oxygen level from 1 PAL to smaller values.

A similar study is performed for enhanced atmospheric CO2 pressures (Pco2). In these conditions, the ozone column is increased at low O2 concentrations with respect to the Pco2 = 1 PAL case. Consequently, the larger CO2 concentration in the ancient atmosphere could have contributed to strengthen the ultraviolet screening of ozone. The surface temperature response to the ozone decrease, as well as the thermal profiles are also analyzed in these CO2-rich models. A possible evolutionary scenario of atmospheric O2 and CO2 is discussed.

The consequences of these calculations for the ultraviolet environment of the primitive biosphere is discussed with a quantitative model calculating bacterial surviving rates. According to this model, the minimum ozone column being tolerable by unprotected bacteria would fall between 1 × 1018 and 4 × 1018 cm−2, depending on the bacterial species considered and corresponding to an O2 level somewhat lower than 10−2 PAL. For the coccoid blue-green alga Agmenellum quadruplicatum, this minimum ozone column would be of 4.5 × 1018, a value which is only slightly less than the presently observed column in the spring time ozone hole of Antarctica.  相似文献   


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

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