Symbolic analytical developments of the zero pressure cosmological model of the universe

In the present paper, literal analytical solutions in power series forms are developed for the radius of curvature and the expansion velocity of the zero pressure cosmological models of the universe at any time t. Also, we develop literal analytical solutions in power series forms for the inverse problem of the zero pressure cosmological model, that is to find the time $t=\tilde{t}$ (say) at which the radius of curvature of the model $R=\tilde{R}$ (say) is known. The importance of these analytical power series representations is that, they are invariant under many operations because, addition, multiplication, exponent ion, integration, different ion, etc of a power series is also a power series. A fact which provides excellent flexibility in dealing with analytical as well as computational developments of the problems related to zero pressure cosmological models.For computational developments of these solutions, an efficient method using continued fraction theory is provided. By means of the present methods we able to analyze some known zero-pressure cosmological models, of these are Einstein and De Sitter models. In addition we also analyzed some other models by which one can know if the universe keep expanding forever, or will it reach a maximal size and then turn into contraction stage.     

修改引力宇宙学下的大尺度结构非线性演化

主要研究了修改引力论在非线性尺度上的表现。现在已经存在很多修改引力论的模型,如f（R）,DGP,Mond等,但是在工作中并没有使用别人已经提出的修改引力论方法,因为首先不可能去对所有的修改引力论模型都进行研究;其次重点研究的是修改引力对于宇宙中的物质分布在非线性尺度上的影响,而各种修改引力论对于物质聚集的影响是相似的,都可以通过一个参数表达出来。因此笔者提出一种修改引力论的方法,引进了一个简单参数ζ。这个参数ζ定义为G_eff/G/-φ/Ψ,这里G_eff是有效牛顿常数,G是牛顿常数,φ和Ψ则是时空度规的标量扰动方程中的两个势。这个ζ参数在广义相对论的基础上修改了物质粒子受到的加速度,加大或者减小这个参数会相应地改变物质粒子的加速度,进而会加快或者减慢大尺度结构的形成。这种对大尺度结构在线性尺度上的影响可以通过理论计算得到,但是大尺度结构的非线性增长无法通过理论计算得出,现在只有借助已经发展得比较成熟的数值模拟方法来研究这个问题。该文修改了现在已经比较成熟的数值模拟程序GADGET-2,加入了ζ参数。在进行数值模拟之前,首先检验了数值模拟的质量分辨率、模拟尺度,以及初始红移对最终物质功率谱的影响,并且证实修改后的GADGET-2程序是正确的。设定不同ζ参数的值,ζ=0.8,0.9,1.0,1.1,1.2,1.5,这里ζ=1.0对应标准宇宙学,然后对每个ζ值运行了一组数值模拟。为了减少数值效应对最终结果的影响,并且突出显示修改引力论在非线性尺度上的表现,没有简单直接地使用不同ζ值物质功率谱来做比较,而是精巧地设置不同ζ值的输出红移,使不同ζ值在其对应的输出红移处有相同的线性功率谱,并且最终使用了ξ=ζ≠1.0/ζ=1.0这个比值参数来研究修改引力论在非线性尺度上的表现。通过宇宙学尺度上的数值模拟和高精度的小波物质功率谱方法,我们研究了ζ参数对物质在非线性尺度上增长的影响。通过ξ这个参数,在线性功率谱相同的时候,比较了不同ζ参数的功率谱在非线性尺度上相对于标准宇宙学模型的偏离。得到重大的发现:大尺度结构在非线性上的增长并不是想象的那样,较大的ζ参数加大了物质粒子的加速度,进而在线性和非线性尺度上都加快了结构形成,反而是在线性功率谱相同的时候,较大值的ζ给出了较小的非线性功率谱。这个发现也证明HKLM假定的不正确,从而证明基于这个假定的功率谱拟合公式（PD96,Smith2003）不适用于修改引力论,直接将这些功率谱拟合公式推广到修改引力论是不正确的。基于这个ζ参数的修改引力论模型,利用一系列的数值模拟的结果,为修改引力论提出了一个新的拟合公式。这个新的拟合公式不是对现有的这些拟合公式做推广,而是借助于标准宇宙学下的非线性功率谱,通过几个简单的参数就可以得到对应不同ζ值的非线性功率谱,而且其精度保证在5%以内。利用这套数值模拟的结果,我们进而研究ζ修改引力论下的本动速度场的功率谱、暗晕的质量方程,以及暗晕的性质。这些工作为我们从观测上区分修改引力论和广义相对论提供了理论指导。     

North and South Hemispheric Solar Activity for Cycles 21 – 23: Asymmetry and Conditional Volatility of Plage Region Areas

We study the dynamic evolution of the time series describing the plage regions areas observed daily at the Observatório Astronómico da Universidade de Coimbra, in each one of the solar hemispheres during solar cycles 21?–?23. The classical ARMA model has proven to be insufficient to describe the time variations seen in the data because of the strong conditional variability. We found that the data are well fitted by ARMA mixed with power-δ TGARCH error models. The power index δ is non-integer; this property has recently been introduced in the literature on time-series analyses and indicates the presence of strong volatility and long memory in the data series. We also detected dynamic asymmetry in the plage region areas observed in the two hemispheres when two different temporal models were obtained to fit them. The finding of a dynamic asymmetry is also supported by the dynamic evolution of the daily difference (north–south) time series, which is significantly different from white noise. This statistical modeling of time series, taking into account new and different characteristics of the solar activity, will be very useful in subsequent forecast developments.     

Formation and evolution of Lakshmi Planum, Venus: Assessment of models using observations from geological mapping

Detailed geological analysis of the Lakshmi Planum region of western Ishtar Terra results in the establishment of the sequence of major events during the formation and evolution of western Ishtar Terra, an important and somewhat unique area on Venus characterized by a raised volcanic plateau surrounded by distinctive folded mountain belts, such as Maxwell Montes. These mapping results and the stratigraphic and structural relationships provide a basis for addressing the complicated problem of Lakshmi Planum formation and for testing the suite of models previously proposed to explain this structure. We review and classify previous models of formation for western Ishtar Terra into “downwelling” models (generally involving convergence and underthrusting) and “upwelling” models (generally involving plume-like upwelling and divergence). The interpreted nature of units and the sequence of events derived from geological mapping are in contrast to the predictions of the divergent models. The major contradictions are as follows: (1) The very likely presence of an ancient (craton-like) tessera massif in the core of Lakshmi, which is inconsistent with the model of formation of Lakshmi due to rise and collapse of a mantle diapir; (2) The absence of rift zones in the interior of Lakshmi that are predicted by the divergent models; (3) The apparent migration of volcanic activity toward the center of Lakshmi, whereas divergent models predict the opposite trend; (4) The abrupt cessation of ridges of the mountain ranges at the edge of Lakshmi Planum and propagation of these ridges over hundreds of kilometers outside Lakshmi; the divergent models predict the opposite progression in the development of major contractional features. In contrast, convergent models of formation and evolution of Lakshmi Planum appear to be more consistent with the observations and explain this structure by collision and underthrusting/subduction of lower-lying plains with the elevated and rigid block of tessera. These models are capable of explaining formation of the major features of western Ishtar (for example, the mountain belts), the sequences of events, and principal volcanic and tectonic trends during the evolution of Lakshmi. To explain the pronounced north-south asymmetry of Lakshmi these models need to consider the likelihood that the major focal points of collision are at the north and north-west margins of the plateau. We note that pure downwelling models, however, face three important difficulties: (1) The possibly unrealistically long time span that appears to be required to produce the major features of Lakshmi; (2) The strong north-south asymmetry of the Planum; the pure downwelling models predict the formation of a more symmetrical structure; and (3) The absence of radial contractional structures (arches and ridges) in the interior of Lakshmi that would represent the predictions of the downwelling models.     

Energetics of the aligned pulsar magnetosphere

An energy analysis is performed on two explicit models, due to Jackson, of a pulsar with aligned magnetic and rotational axes. The unknown parameters of these models are determined by calculating the minimum total energy states of the models. It is found that the minimum energy analysis favors states with extended, dynamically active magnetospheres with a high degrees of corotation. By calculating total power input to the magnetosphere via collisions in the stellar crust, and the total power radiated due to azimuthal drift motion, it is determined that the minimum energy states are the only states where a power balance can be achieved. Consideration of a local power balance condition and dissipative flows in the magnetosphere shows that neither model is completely self-consistent, but one is considerably better than the other. Properties of both models and implications for other models are discussed.     

Scale Factor of Isotropic Homogeneous Universe

Four fundamental expressions of isotropic homogeneous universe lead to the scale factor equations. Their solution is achieved by developing the factor as a power series of time, establishing the recursion relations among its coefficients, performing the summation of this power series and obtaining the scale factor as a function of time in a closed form. This one should be used in comparison to observational data with different cosmological models, generally present within this formalism.     

电离层延迟改正模型综述

电离层延迟改正模型通常可以分为广播星历用的预报模型、广域差分用的实时模型、后处理模型3类,不同应用要求需要选择不同的模型。主要比较分析了几种常用的电离层延迟改正模型: 用于广域差分中生成格网模型的三角级数模型、多项式模型、低阶球谐函数模型等都可以获得很好的改正效果,且这3种模型基本上是等价的;电离层延迟谐函数展开模型可以用来分析电离层长时间系列的变化特征;国际电离层参考模型IRI的改正精度一般可以达到60％的效果;而GPS 星历采用的Klobuchar模型的参数设置存在一些不足,对此提出了一些改进措施。     

Power law relating 10.7 cm flux to sunspot number

To investigate the relation between observations of the 10.7 cm flux and the international sunspot number so that a physical unit may be ascribed to historical records, both polynomial and power law models are developed giving the radio flux as a function of sunspot number and vice versa. Bayesian data analysis is used to estimate the model parameters and to discriminate between the models. The effect on the parameter uncertainty and on the relative evidence of normalizing the measure of fit is investigated. The power law giving flux as a function of sunspot number is found to be the most plausible model and may be used to estimate the radio flux from historical sunspot observations.     

On the r-mode spectrum of relativistic stars in the low-frequency approximation

The axial modes for non-barotropic relativistic rotating neutron stars with uniform angular velocity are studied, using the slow-rotation formalism together with the low-frequency approximation, first investigated by Kojima. The time-independent form of the equations leads to a singular eigenvalue problem, which admits a continuous spectrum. We show that for     , it is nevertheless also possible to find discrete mode solutions (the r modes). However, under certain conditions related to the equation of state and the compactness of the stellar model, the eigenfrequency lies inside the continuous band and the associated velocity perturbation is divergent; hence these solutions have to be discarded as being unphysical. We corroborate our results by explicitly integrating the time-dependent equations. For stellar models admitting a physical r-mode solution, it can indeed be excited by arbitrary initial data. For models admitting only an unphysical mode solution, the evolutions do not show any tendency to oscillate with the respective frequency. For higher values of l it seems that in certain cases there are no mode solutions at all.     

基于修正线性组合模型的原子钟钟差预报

原子钟钟差预报在时频工作中起着重要作用.目前常用的预报模型各有优缺点,为了综合各种预报模型特点,可以考虑组合这些模型预报的结果.针对线性组合模型预报钟差时存在的问题,提出了学习权的概念,以有效地利用各种精度信息,建立了修正组合预报模型.为了验证该方法的有效性,选取了4颗GPS卫星IGS(International GN...     

Power-series solutions of the Lane-Emden equation

In this paper the problem of representing the numerical solutions of the Lane-Emden equation analytically by means of a convergent power series has been considered. Our results show that it is possible to represent the numerical solutions of the Lane-Emden equation by means of a power series which can be convergent in the whole interior of a polytropic model.     

The cluster distribution as a test of dark matter models — IV. Topology and geometry

We study the geometry and topology of the large-scale structure traced by galaxy clusters in numerical simulations of a box of side 320 h ⁻¹ Mpc, and compare them with available data on real clusters. The simulations we use are generated by the Zel'dovich approximation, using the same methods as we have used in the first three papers in this series. We consider the following models to see if there are measurable differences in the topology and geometry of the superclustering they produce: (i) the standard cold dark matter model (SCDM); (ii) a CDM model with Ω₀ = 0.2 (OCDM); (iii) a CDM model with a 'tilted' power spectrum having n  = 0.7 (TCDM); (iv) a CDM model with a very low Hubble constant, h  = 0.3 (LOWH); (v) a model with mixed CDM and HDM (CHDM); (vi) a flat low-density CDM model with Ω₀ = 0.2 and a non-zero cosmological Λ term (ΛCDM). We analyse these models using a variety of statistical tests based on the analysis of: (i) the Euler–Poincaré characteristic; (ii) percolation properties; (iii) the minimal spanning tree construction. Taking all these tests together we find that the best-fitting model is ΛCDM and, indeed, the others do not appear to be consistent with the data. Our results demonstrate that despite their biased and extremely sparse sampling of the cosmological density field, it is possible to use clusters to probe subtle statistical diagnostics of models, which go far beyond the low-order correlation functions usually applied to study superclustering.     

Models for jet power in elliptical galaxies: a case for rapidly spinning black holes

The power of jets from black holes is expected to depend on both the spin of the black hole and the structure of the accretion disc in the region of the last stable orbit. We investigate these dependencies using two different physical models for the jet power: the classical Blandford–Znajek (BZ) model and a hybrid model developed by Meier. In the BZ case, the jets are powered by magnetic fields directly threading the spinning black hole while in the hybrid model, the jet energy is extracted from both the accretion disc as well as the black hole via magnetic fields anchored to the accretion flow inside and outside the hole's ergosphere. The hybrid model takes advantage of the strengths of both the Blandford–Payne and BZ mechanisms, while avoiding the more controversial features of the latter. We develop these models more fully to account for general relativistic effects and to focus on advection-dominated accretion flows (ADAFs) for which the jet power is expected to be a significant fraction of the accreted rest mass energy.
We apply the models to elliptical galaxies, in order to see if these models can explain the observed correlation between the Bondi accretion rates and the total jet powers. For typical values of the disc viscosity parameter  α∼ 0.04 –0.3  and mass accretion rates consistent with ADAF model expectations, we find that the observed correlation requires   j ≳ 0.9  ; that is, it implies that the black holes are rapidly spinning. Our results suggest that the central black holes in the cores of clusters of galaxies must be rapidly rotating in order to drive jets powerful enough to heat the intracluster medium and quench cooling flows.     

Some Homogeneous Bianchi type IX Viscous Fluid Cosmological Models with a Varying Λ

Some Bianchi type IX viscous fluid cosmological models are investigated. To get a solution, a supplementary condition between metric potentials is used. The viscosity coefficient of bulk viscous fluid is assumed to be a power function of mass density, whereas the coefficient of shear viscosity is considered as proportional to scale of expansion in the model. The cosmological constant Λ is found to be positive and is a decreasing function of time, which is supported by results from recent supernovae observations. Some physical and geometric properties of the models are also discussed.     

On the structure and composition of Uranus and Neptune

We present a series of models of Uranus and Neptune in which the relative amounts of (1) rock, (2) ices, and (3) hydrogen and helium are allowed to vary. By fitting the density and the gravitational quadruppole moment, the model composition can be determined. Because of the ambiguity in the rotation periods of these planets, several possible models are presented and discussed.     

Modified gravity in a viscous and non-isotropic background

We study the dynamical evolution of an f(R) model of gravity in a viscous and anisotropic background which is given by a Bianchi type-I model of the Universe. We find viable forms of f(R) gravity in which one is exactly the Einsteinian model of gravity with a cosmological constant and other two are power law f(R) models. We show that these two power law models are stable with a suitable choice of parameters. We also examine three potentials which exhibit the potential effect of f(R) models in the context of scalar tensor theory. By solving different aspects of the model and finding the physical quantities in the Jordan frame, we show that the equation of state parameter satisfy the dominant energy condition. At last we show that the two power law f(R) models behave like quintessence model at late times and also the shear coefficient viscosity tends to zero at late times.     

Autoregressive Spectral Estimation for Quasi-Periodic Oscillations

Modern methods of spectral estimation based on parametric time-series models are useful tools in power spectral analysis. We apply the autoregressive (AR) model to study quasi-periodic oscillations (QPOs). An empirical formula to estimate the expectation and standard deviation of the noise AR power densities is derived, which can be used to estimate the statistical significance of an apparent QPO peak in an AR spectrum. An iterative adding-noise algorithm in AR spectral analysis is proposed and applied to studying QPOs in the X-ray binary Cir X-1.     

Warm dark matter versus bumpy power spectra

In this paper we explore the differences between a warm dark matter (WDM) model and a cold dark matter (CDM) model where the power on a certain scale is reduced by introducing a narrow negative feature ('dip'). This dip is placed in a way so as to mimic the loss of power in the WDM model: both models have the same integrated power out to the scale where the power of the dip model rises to the level of the unperturbed CDM spectrum again.
Using N -body simulations we show that some of the large-scale clustering patterns of this new model follow more closely the usual CDM scenario while simultaneously suppressing small-scale structures (within galactic haloes) even more efficiently than WDM. The analysis in the paper shows that the new Dip model appears to be a viable alternative to WDM, but it is based on different physics. Where WDM requires the introduction of a new particle species, the Dip model is based on a non-standard inflationary period. If we are looking for an alternative to the currently challenged standard ΛCDM structure formation scenario, neither the ΛWDM nor the new Dip model can be ruled out based on the analysis presented in this paper. They both make very similar predictions and the degeneracy between them can only be broken with observations yet to come.     

JMAGCM模拟的大气角动量函数Ⅰ.季节性分量(英文)

大气角动量(AAM) 变化和地球自转变化密切相关,它们可应用于全球环流模型的检验.根据日本气象厅全球环流模型的输出数据计算出AAM 函数,将AAM 函数的模拟值与它的观测值以及地球自转变化序列相比较,来检验该模型对大气角动量的模拟.限于篇幅,仅给出AAM 函数的季节性分量模拟的结果.对于半年项,轴向AAM 函数的气压项和风项的模拟值都偏高,赤道向AAM 函数的气压项和风项很小,因而未予考虑;对于周年项,轴向AAM函数气压项的模拟值偏低,而赤道向AAM 函数气压项的模拟值偏高,轴向AAM 函数的风项模拟得最好,而赤道向AAM 函数的风项模拟得较差.总之,高精度的天文观测的地球自转系列可以作为全球环流模型检验的一种重要参考依据.     

Observational tests of x-matter models

We study gravitational lensing statistics, matter power spectra and the angular power spectra of the cosmic microwave background (CMB) radiation in x-matter models. We adopt an equation of state of x-matter which can express a wide range of matter from pressureless dust to the cosmological constant. A new ingredient in this model is the sound speed of the x-component, in addition to the equation of state w ₀ =  p _x0/ρ_x0. Except for the cosmological constant case, the perturbations of x-matter itself are considered. Our primary interest is in the effect of non-zero sound speed on the structure formation and the CMB spectra. It is found that there exist parameter ranges where x-matter models are consistent with all current observations. The x-matter generally leaves imprints in the CMB anisotropy and the matter power spectrum, which should be detectable in future observations.