This work presents a systematic investigation of the influence of weather conditions on the calibration errors by using Gaussian fitness, least chi-square linear fitness and wavelet transform to analyze the calibration coefficients from observations of the Chinese Solar Broadband Radio Spectrometers (at frequency bands of 1.0-2.0 GHz, 2.6-3.8 GHz, and 5.2-7.6 GHz) during 1997-2007. We found that calibration coefficients are influenced by the local air temperature. Considering the temperature correction, the calibration error will reduce by about $10\%-20\%$ at 2800 MHz. Based on the above investigation and the calibration corrections, we further study the radio emission of the quiet-Sun by using an appropriate hybrid model of the quiet-Sun atmosphere. The results indicate that the numerical flux of the hybrid model is much closer to the observation flux than that of other ones.

We study multifield contributions to the scalar power spectrum in an ensemble of six-field inflationary models obtained in string theory. We identify examples in which inflation occurs by chance, near an approximate inflection point, and we compute the primordial perturbations numerically, both exactly and using an array of truncated models. The scalar mass spectrum and the number of fluctuating fields are accurately described by a simple random matrix model. During the approach to the inflection point, bending trajectories and violations of slow roll are commonplace, and 'many-field' effects, in which three or more fields influence the perturbations, are often important. However, in a large fraction of models consistent with constraints on the tilt the signatures of multifield evolution occur on unobservably large scales. Our scenario is a concrete microphysical realization of quasi-single-field inflation, with scalar masses of order $H$, but the cubic and quartic couplings are typically too small to produce detectable non-Gaussianity. We argue that our results are characteristic of a broader class of models arising from multifield potentials that are natural in the Wilsonian sense.

We study the six-field dynamics of D3-brane inflation for a general scalar potential on the conifold, finding simple, universal behavior. We numerically evolve the equations of motion for an ensemble of more than 7 \times 10^7 realizations, drawing the coefficients in the scalar potential from statistical distributions whose detailed properties have demonstrably small effects on our results. When prolonged inflation occurs, it has a characteristic form: the D3-brane initially moves rapidly in the angular directions, spirals down to an inflection point in the potential, and settles into single-field inflation. The probability of N_{e} e-folds of inflation is a power law, P(N_{e}) \propto N_{e}^{-3}, and we derive the same exponent from a simple analytical model. The success of inflation is relatively insensitive to the initial conditions: we find attractor behavior in the angular directions, and the D3-brane can begin far above the inflection point without overshooting. In favorable regions of the parameter space, models yielding 60 e-folds of expansion arise approximately once in 10^3 trials. Realizations that are effectively single-field and give rise to a primordial spectrum of fluctuations consistent with WMAP, for which at least 120 e-folds are required, arise approximately once in 10^5 trials. The emergence of robust predictions from a six-field potential with hundreds of terms invites an analytic approach to multifield inflation.

The Tree-Particle-Mesh (TPM) N-body algorithm couples the tree algorithm for directly computing forces on particles in an hierarchical grouping scheme with the extremely efficient mesh based PM structured approach. The combined TPM algorithm takes advantage of the fact that gravitational forces are linear functions of the density field. Thus one can use domain decomposition to break down the density field into many separate high density regions containing a significant fraction of the mass but residing in a very small fraction of the total volume. In each of these high density regions the gravitational potential is computed via the tree algorithm supplemented by tidal forces from the external density distribution. For the bulk of the volume, forces are computed via the PM algorithm; timesteps in this PM component are large compared to individually determined timesteps in the tree regions. Since each tree region can be treated independently, the algorithm lends itself to very efficient parallelization using message passing. We have tested the new TPM algorithm (a refinement of that originated by Xu 1995) by comparison with results from Ferrell & Bertschinger's P^3M code and find that, except in small clusters, the TPM results are at least as accurate as those obtained with the well-established P^3M algorithm, while taking significantly less computing time. Production runs of 10^9 particles indicate that the new code has great scientific potential when used with distributed computing resources.

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

We describe an automated method for detecting clusters of galaxies in imaging and redshift galaxy surveys. The Adaptive Matched Filter (AMF) method utilizes galaxy positions, magnitudes, and---when available---photometric or spectroscopic redshifts to find clusters and determine their redshift and richness. The AMF can be applied to most types of galaxy surveys: from two-dimensional (2D) imaging surveys, to multi-band imaging surveys with photometric redshifts of any accuracy (2.5D), to three-dimensional (3D) redshift surveys. The AMF can also be utilized in the selection of clusters in cosmological N-body simulations. The AMF identifies clusters by finding the peaks in a cluster likelihood map generated by convolving a galaxy survey with a filter based on a model of the cluster and field galaxy distributions. In tests on simulated 2D and 2.5D data with a magnitude limit of r' ~ 23.5, clusters are detected with an accuracy of Delta z ~ 0.02 in redshift and ~10% in richness to z < 0.5. Detecting clusters at higher redshifts is possible with deeper surveys. In this paper we present the theory behind the AMF and describe test results on synthetic galaxy catalogs.

The theory and application of numerical methods for unstructured meshes have been improved significantly in recent years. Because the grids can be place arbitrarily in space, unstructured meshes can provide much higher spatial resolution than regular meshes. The built-in nature of mesh adaptivity for unstructured meshes gives one way to simulate highly dynamic, hierarchical problems involving both collisionless dark matter and collisional gas dynamics. In this paper, we describe algorithms to construct unstructured meshes from a set of points with periodic boundary conditions through Delaunay triangulation, and algorithms to solve hydrodynamic and N-body problems on an unstructured mesh . A combination of a local transformation algorithm and the traditional Bowyer-Watson algorithm gives an efficient approach to perform Delaunay triangulation. A novel algorithm to solve N-body equations of motion on an unstructured mesh is described. Poisson's equation is solved using the conjugate gradient method. A gas-kinetic scheme based on the BGK model to solve Euler equations is used to evolve the hydrodynamic equations. We apply these algorithms to solve cosmological settings, which involve both dark and baryonic matter. Various cooling and heating processes for primordial baryonic matter are included in the code. The numerical results show that the N-body and hydrodynamic algorithms based on unstructured meshes with mesh refinement are well-suited for hierarchical structure formation problems.

We examine the sensitivity of the spatial morphologies of galaxy clusters to Omega_0 and P(k) using high-resolution N-body simulations with large dynamic range. Variants of the standard CDM model are considered having different spatial curvatures, SCDM (Omega_0=1), OCDM (Omega_0=0.35), LCDM (Omega_0=0.35, lambda_0=0.65), and different normalizations, sigma_8. We also explore critical density models with different spectral indices, n, of the scale-free power spectrum, P(k) propto k^n. Cluster X-ray morphologies are quantified with power ratios (PRs), where we take for the X-ray emissivity j_gas propto rho^2_DM, which we argue is a suitable approximation for analysis of PRs. We find that Omega_0 primarily influences the means of the PR distributions whereas the power spectrum (n and sigma_8) primarily affects their variances: log_10(P_3/P_0) is the cleanest probe of Omega_0 since its mean is very sensitive to Omega_0 but very insensitive to P(k). The PR means easily distinguish the SCDM and OCDM models, while the SCDM and LCDM means show a more modest, but significant, difference (3 sigma). The OCDM and LCDM models are largely indistinguishable in terms of the PRs. Finally, we compared these models to a sample of ROSAT clusters and find that the cluster morphologies strongly disfavor Omega_0=1, CDM while favoring low density, CDM models with or without a cosmological constant.

Weak gravitational lensing by large-scale structure affects the determination of the cosmological deceleration parameter $q_0$. We find that the lensing induced dispersions on truly standard candles are $0.04$ and $0.02$ mag at redshift $z=1$ and $z=0.5$, respectively, in a COBE-normalized cold dark matter universe with $\Omega_0=0.40$, $\Lambda_0=0.6$, $H=65$km/s/Mpc and $\sigma_8=0.79$. It is shown that one would observe $q_0=-0.44^{+0.17}_{-0.05}$ and $q_0=-0.45^{+0.10}_{-0.03}$ (the errorbars are $2\sigma$ limits) with standard candles with zero intrinsic dispersion at redshift $z=1$ and $z=0.5$, respectively, compared to the truth of $q_0=-0.40$ in this case, i.e., a 10\% error in $q_0$ will be made. A standard COBE normalized $\Omega_0=1$ CDM model would produce three times as much variance and a mixed (hot and cold) dark matter model would lead to an intermediate result. One unique signature of this dispersion effect is its non Gaussianity. Although the lensing induced dispersion at lower redshift is still significantly smaller than the currently best observed (total) dispersion of $0.12$ mag in a sample of type Ia supernovae, selected with the multicolor light curve shape method, it becomes significant at higher redshift. We show that there is an optimal redshift, in the range $z\sim 0.5-2.0$ depending on the amplitude of the intrinsic dispersion of the standard candles, at which $q_0$ can be most accurately determined.

We have developed a gravity solver based on combining the well developed Particle-Mesh (PM) method and TREE methods. It is designed for and has been implemented on parallel computer architectures. The new code can deal with tens of millions of particles on current computers, with the calculation done on a parallel supercomputer or a group of workstations. Typically, the spatial resolution is enhanced by more than a factor of 20 over the pure PM code with mass resolution retained at nearly the PM level. This code runs much faster than a pure TREE code with the same number of particles and maintains almost the same resolution in high density regions. Multiple time step integration has also been implemented with the code, with second order time accuracy. The performance of the code has been checked in several kinds of parallel computer configuration, including IBM SP1, SGI Challenge and a group of workstations, with the speedup of the parallel code on a 32 processor IBM SP2 supercomputer nearly linear (efficiency $\approx 80\%$) in the number of processors. The computation/communication ratio is also very high ($\sim 50$), which means the code spends $95\%$ of its CPU time in computation.