Cosmic voids are biased tracers of the large-scale structure of the universe. Separate universe simulations (SUS) enable accurate measurements of this biasing relation by implementing the peak-background split (PBS). In this work, we apply the SUS technique to measure the void bias parameters. We confirm that the PBS argument works well for underdense tracers. The response of the void size distribution depends on the void radius. For voids larger (smaller) than the size at the peak of the distribution, the void abundance responds negatively (positively) to a long wavelength mode. The linear bias from the SUS is in good agreement with the cross power spectrum measurement on large scales. Using the SUS, we have detected the quadratic void bias for the first time in simulations. We find that $ b_2 $ is negative when the magnitude of $ b_1 $ is small, and that it becomes positive and increases rapidly when $ |b_1| $ increases. We compare the results from voids identified in the halo density field with those from the dark matter distribution, and find that the results are qualitatively similar, but the biases generally shift to the larger voids sizes.

We study the large-scale bias parameter of cosmic voids with primordial non-Gaussian (PNG) initial conditions of the local type. In this scenario, the dark matter halo bias exhibits a characteristic scale dependence on large scales, which has been recognized as one of the most promising probes of the local PNG. Using a suite of $N$-body simulations with Gaussian and non-Gaussian initial conditions, we find that the void bias features scale-dependent corrections on large scales, similar to its halo counterpart. We find excellent agreement between the numerical measurement of the PNG void bias and the general peak-background split prediction. Contrary to halos, large voids anti-correlate with the dark matter density field, and the large-scale Gaussian void bias ranges from positive to negative values depending on void size and redshift. Thus, the information in the clustering of voids can be complementary to that of the halos. Using the Fisher matrix formalism for multiple tracers, we demonstrate that including the scale-dependent bias information from voids, constraints on the PNG parameter $f_{\rm NL}$ can be tightened by a factor of two compared to the accessible information from halos alone, when the sampling density of tracers reaches $4 \times 10^{-3} \, h^3 \mathrm{Mpc}^{-3} $.

The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $\Omega_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade.

We define and characterise a sample of 1.3 million galaxies extracted from the first year of Dark Energy Survey data, optimised to measure Baryon Acoustic Oscillations in the presence of significant redshift uncertainties. The sample is dominated by luminous red galaxies located at redshifts $z \gtrsim 0.6$. We define the exact selection using color and magnitude cuts that balance the need of high number densities and small photometric redshift uncertainties, using the corresponding forecasted BAO distance error as a figure-of-merit in the process. The typical photo-$z$ uncertainty varies from $2.3\%$ to $3.6\%$ (in units of 1+$z$) from $z=0.6$ to $1$, with number densities from $200$ to $130$ galaxies per deg$^2$ in tomographic bins of width $\Delta z = 0.1$. Next we summarise the validation of the photometric redshift estimation. We characterise and mitigate observational systematics including stellar contamination, and show that the clustering on large scales is robust in front of those contaminants. We show that the clustering signal in the auto-correlations and cross-correlations is generally consistent with theoretical models, which serves as an additional test of the redshift distributions.

We present angular diameter distance measurements obtained by locating the BAO scale in the distribution of galaxies selected from the first year of Dark Energy Survey data. We consider a sample of over 1.3 million galaxies distributed over a footprint of 1318 deg$^2$ with $0.6 < z_{\rm photo} < 1$ and a typical redshift uncertainty of $0.03(1+z)$. This sample was selected, as fully described in a companion paper, using a color/magnitude selection that optimizes trade-offs between number density and redshift uncertainty. We investigate the BAO signal in the projected clustering using three conventions, the angular separation, the co-moving transverse separation, and spherical harmonics. Further, we compare results obtained from template based and machine learning photometric redshift determinations. We use 1800 simulations that approximate our sample in order to produce covariance matrices and allow us to validate our distance scale measurement methodology. We measure the angular diameter distance, $D_A$, at the effective redshift of our sample divided by the true physical scale of the BAO feature, $r_{\rm d}$. We obtain close to a 4 per cent distance measurement of $D_A(z_{\rm eff}=0.81)/r_{\rm d} = 10.75\pm 0.43 $. These results are consistent with the flat $\Lambda$CDM concordance cosmological model supported by numerous other recent experimental results. All data products are publicly available here: https://des.ncsa.illinois.edu/releases/y1a1/bao

We study the methodology and potential theoretical systematics of measuring Baryon Acoustic Oscillations (BAO) using the angular correlation functions in tomographic bins. We calibrate and optimize the pipeline for the Dark Energy Survey Year 1 dataset using 1800 mocks. We compare the BAO fitting results obtained with three estimators: the Maximum Likelihood Estimator (MLE), Profile Likelihood, and Markov Chain Monte Carlo. The fit results from the MLE are the least biased and their derived 1-$\sigma$ error bar are closest to the Gaussian distribution value after removing the extreme mocks with non-detected BAO signal. We show that incorrect assumptions in constructing the template, such as mismatches from the cosmology of the mocks or the underlying photo-$z$ errors, can lead to BAO angular shifts. We find that MLE is the method that best traces this systematic biases, allowing to recover the true angular distance values. In a real survey analysis, it may happen that the final data sample properties are slightly different from those of the mock catalog. We show that the effect on the mock covariance due to the sample differences can be corrected with the help of the Gaussian covariance matrix or more effectively using the eigenmode expansion of the mock covariance. In the eigenmode expansion, the eigenmodes are provided by some proxy covariance matrix. The eigenmode expansion is significantly less susceptible to statistical fluctuations relative to the direct measurements of the covariance matrix because of the number of free parameters is substantially reduced

Mock catalogues are a crucial tool in the analysis of galaxy surveys data, both for the accurate computation of covariance matrices, and for the optimisation of analysis methodology and validation of data sets. In this paper, we present a set of 1800 galaxy mock catalogues designed to match the Dark Energy Survey Year-1 BAO sample (Crocce et al. 2017) in abundance, observational volume, redshift distribution and uncertainty, and redshift dependent clustering. The simulated samples were built upon HALOGEN (Avila et al. 2015) halo catalogues, based on a $2LPT$ density field with an exponential bias. For each of them, a lightcone is constructed by the superposition of snapshots in the redshift range $0.45<z<1.4$. Uncertainties introduced by so-called photometric redshifts estimators were modelled with a \textit{double-skewed-Gaussian} curve fitted to the data. We also introduce a hybrid HOD-HAM model with two free parameters that are adjusted to achieve a galaxy bias evolution $b(z_{\rm ph})$ that matches the data at the 1-$\sigma$ level in the range $0.6<z_{\rm ph}<1.0$. We further analyse the galaxy mock catalogues and compare their clustering to the data using the angular correlation function $ w(\theta)$, the comoving transverse separation clustering $\xi_{\mu<0.8}(s_{\perp})$ and the angular power spectrum $C_\ell$.

We use data from the first-year (Y1) observations of the DES collaboration to measure the galaxy angular power spectrum (APS), and search for its BAO feature using a template-fitting method. We test our methodology in a sample of 1800 DES Y1-like mock catalogs. The APS is measured with the pseudo-$C_\ell$ method, using pixelized maps constructed from the mock catalogs and the DES mask. The covariance matrix of the $C_\ell$'s in these tests are also obtained from the mock catalogs. We use templates to model the measured spectra and estimate template parameters firstly from the $C_\ell$'s of the mocks using two different methods, a maximum likelihood estimator and a MCMC, finding consistent results with a good reduced $\chi^2$. Robustness tests are performed to estimate the impact of different choices of settings used in our analysis. After these tests on mocks, we apply our method to a galaxy sample constructed from DES Y1 data specifically for LSS studies. This catalog comprises galaxies within an effective area of 1318 deg$^2$ and $0.6<z<1.0$. We fit the observed spectra with our optimized templates, considering models with and without BAO features. We find that the DES Y1 data favors a model with BAO at the $2.6\,\sigma$ C.L. with a best-fit shift parameter of $\alpha=1.023\pm 0.047$. However, the goodness-of-fit is somewhat poor, with $\chi^2/$(dof) = 1.49. We identify a possible cause of this issue and show that using a theoretical covariance matrix obtained from $C_\ell$'s that are better adjusted to data results in an improved value of $\chi^2/$(dof) = 1.36 which is similar to the value obtained with the real-space analysis. Our results correspond to a distance measurement of $D_A(z_{\rm eff} = 0.81)/r_d = 10.65 \pm 0.49$, consistent with the main DES BAO findings. This is a companion paper to the main DES BAO article showing the details of the harmonic space analysis.

Modes with wavelengths larger than the survey window can have significant impact on the covariance within the survey window. The supersample covariance has been recognized as an important source of covariance for the power spectrum on small scales, and it can potentially be important for the bispectrum covariance as well. In this paper, using the response function formalism, we model the supersample covariance contributions to the bispectrum covariance and the cross covariance between the power spectrum and the bispectrum. The supersample covariances due to the long wavelength density and tidal perturbations are investigated, and the tidal contribution is a few orders of magnitude smaller than the density one because in configuration space the bispectrum estimator involves angular averaging and the tidal response function is anisotropic. The impact of the super-survey modes is quantified using numerical measurements with periodic box and subbox setups. For the matter bispectrum, the ratio between the supersample covariance correction and the small scale covariance, which can be computed using a periodic box, is roughly an order of magnitude smaller than that for the matter power spectrum. This is because for the bispectrum, the small scale non-Gaussian covariance is significantly larger than that for the power spectrum. For the cross covariance, the supersample covariance is as important as for the power spectrum covariance. The supersample covariance prediction with the halo model response function is in good agreement with numerical results.

The window function for protohalos in Lagrangian space is often assumed to be a tophat in real space. We measure this profile directly and find that it is more extended than a tophat but less extended than a Gaussian; its shape is well-described by rounding the edges of the tophat by convolution with a Gaussian that has a scale length about 5 times smaller. This effective window $W_{\rm eff}$ is particularly simple in Fourier space, and has an analytic form in real space. Together with the excursion set bias parameters, $W_{\rm eff}$ describes the scale-dependence of the Lagrangian halo-matter cross correlation up to $kR_{\rm Lag} \sim 10 $, where $R_{\rm Lag}$ is the Lagrangian size of the protohalo. Moreover, with this $W_{\rm eff}$, all the spectral moments of the power spectrum are finite, allowing a straightforward estimate of the excursion set peak mass function. This estimate requires a prescription of the critical overdensity enclosed within a protohalo if it is to collapse, which we calibrate from simulations. We find that the resulting estimate of halo abundances is only accurate to about 20%, and we discuss why: A tophat in `infall time' towards the protohalo center need not correspond to a tophat in the initial spatial distribution, so models in which infall rather than smoothed overdensity is the relevant variable may be more accurate.

Shot noise is an important ingredient to any measurement or theoretical modeling of discrete tracers of the large scale structure. Recent work has shown that the shot noise in the halo power spectrum becomes increasingly sub-Poissonian at high mass. Interestingly, while the halo model predicts a shot noise power spectrum in qualitative agreement with the data, it leads to an unphysical white noise in the cross halo-matter and matter power spectrum. In this work, we show that absorbing all the halo model sources of shot noise into the halo fluctuation field leads to meaningful predictions for the shot noise contributions to halo clustering statistics and remove the unphysical white noise from the cross halo-matter statistics. Our prescription straightforwardly maps onto the general bias expansion, so that the renormalized shot noise terms can be expressed as combinations of the halo model shot noises. Furthermore, we demonstrate that non-Poissonian contributions are related to volume integrals over correlation functions and their response to long-wavelength density perturbations. This leads to a new class of consistency relations for discrete tracers, which appear to be satisfied by our reformulation of the halo model. We test our theoretical predictions against measurements of halo shot noise bispectra extracted from a large suite of numerical simulations. Our model reproduces qualitatively the observed sub-Poissonian noise, although it underestimates the magnitude of this effect.

The covariance matrix of the matter and halo power spectrum and bispectrum are studied. Using a large suite of simulations, we find that the non-Gaussianity in the covariance is significant already at mildly nonlinear scales. We compute the leading disconnected non-Gaussian correction to the matter bispectrum covariance using perturbation theory, and find that the corrections result in good agreement in the mildly nonlinear regime. The shot noise contribution to the halo power spectrum and bispectrum covariance is computed using the Poisson model, and the model yields decent agreement with simulation results. However, when the shot noise is estimated from the individual realization, which is usually done in reality, we find that the halo covariance is substantially reduced and gets close to the Gaussian covariance. This is because most of the non-Gaussianity in the covariance arises from the fluctuations in the Poisson shot noise. We use the measured non-Gaussian covariance to access the information content of the power spectrum and bispectrum. The signal-to-noise ratio, S/N, of the matter and halo power spectrum levels off in the mildly nonlinear regime, $k \sim 0.1 - 0.2 \,\mathrm{Mpc} h^{-1}$. In the nonlinear regime the S/N of the matter and halo bispectrum increases but much slower than the Gaussian results suggest. We find that both the S/N for power spectrum and bispectrum are overestimated by the Gaussian covariances, but the problem being much more serious for the bispectrum. Because the bispectrum is affected strongly by nonlinearity and shot noise, inclusion of the bispectrum only adds modest amount of S/N compared to that of the power spectrum.

The protohalo patches from which halos form are defined by a number of constraints imposed on the Lagrangian dark matter density field. Each of these constraints contributes to biasing the spatial distribution of the protohalos relative to the matter. We show how measurements of this spatial distribution -- linear combinations of protohalo bias factors -- can be used to make inferences about the physics of halo formation. Our analysis exploits the fact that halo bias factors satisfy consistency relations which encode this physics, and that these relations are the same even for sub-populations in which assembly bias has played a role. We illustrate our methods using a model in which three parameters matter: a density threshold, the local slope and the curvature of the smoothed density field. The latter two are nearly degenerate; our approach naturally allows one to build an accurate effective two-parameter model for which the consistency relations still apply. This, with an accurate description of the smoothing window, allows one to describe the protohalo-matter cross-correlation very well, both in Fourier and configuration space. We then use our determination of the large scale bias parameters together with the consistency relations, to estimate the enclosed density and mean slope on the Lagrangian radius scale of the protohalos. Direct measurements of these quantities, made on smaller scales than those on which the bias parameters are typically measured, are in good agreement.

We study the halo-matter cross bispectrum in the presence of primordial non-Gaussianity of the local type. We restrict ourselves to the squeezed limit, for which the calculation are straightforward, and perform the measurements in the initial conditions of N-body simulations, to mitigate the contamination induced by nonlinear gravitational evolution. Interestingly, the halo-matter cross bispectrum is not trivial even in this simple limit as it is strongly sensitive to the scale-dependence of the quadratic and third-order halo bias. Therefore, it can be used to test biasing prescriptions. We consider three different prescription for halo clustering: excursion set peaks (ESP), local bias and a model in which the halo bias parameters are explicitly derived from a peak-background split. In all cases, the model parameters are fully constrained with statistics other than the cross bispectrum. We measure the cross bispectrum involving one halo fluctuation field and two mass overdensity fields for various halo masses and collapse redshifts. We find that the ESP is in reasonably good agreement with the numerical data, while the other alternatives we consider fail in various cases. This suggests that the scale-dependence of halo bias also is a crucial ingredient to the squeezed limit of the halo bispectrum.

An attempt at understanding downward overshooting in the convective envelopes of post-main-sequence stars has been made on the basis of three-dimensional large-eddy simulations, using artificially modified OPAL opacity and taking into account radiation and ionization in the equation of state. Two types of star, an intermediate-mass star and a massive star, were considered. To avoid a long thermal relaxation time of the intermediate-mass star, we increased the stellar energy flux artificially while trying to maintain a structure close to the one given by a 1D stellar model. A parametric study of the flux factor was performed. For the massive star, no such process was necessary. Numerical results were analysed when the system reached the statistical steady state. It was shown that the penetration distance in pressure scaleheights is of the order of unity. The scaling relations between penetration distance, input flux and vertical velocity fluctuations studied by Singh et al. were checked. The anisotropy of the turbulent convection and the diffusion models of the third-order moments representing the non-local transport were also investigated. These models are dramatically affected by the velocity fields and no universal constant parameters seem to exist. The limitations of the numerical results were also discussed.

We propose a simple model that elucidates the generation of halo velocity bias. The fluid equation approximation is often adopted in modelling the evolution of the halo density field. In this approach, halos are often taken to be point particles even though in reality they are finite-sized objects. In this paper, we generalize the fluid equation approximation to halos to include the finite extent of halos by taking into account the halo profile. We compute the perturbation of the halo density and velocity field to second order and find that the profile correction gives rise to $k^2$ correction terms in Fourier space. These corrections are more important for velocity than for density. In particular, the profile correction generates $k^2$ correction term in the velocity bias and the correction terms do not decay away in the long term limit, but it is not constant. We model the halo profile evolution using the spherical collapse model. We also measure the evolution of proto-halo profile at various redshifts numerically. We find that the spherical collapse model gives a reasonable description of the numerical profile evolution. Static halo profile is often adopted in modelling halos in theories such as the excursion set theory. Our work highlights the importance of including the profile evolution in the calculations.

We study the clustering of voids using $N$-body simulations and simple theoretical models. The excursion-set formalism describes fairly well the abundance of voids identified with the watershed algorithm, although the void formation threshold required is quite different from the spherical collapse value. The void cross bias $b_{\rm c} $ is measured and its large-scale value is found to be consistent with the peak background split results. A simple fitting formula for $b_{\rm c} $ is found. We model the void auto-power spectrum taking into account the void biasing and exclusion effect. A good fit to the simulation data is obtained for voids with radii $\gtrsim$ 30 Mpc/$h$, especially when the void biasing model is extended to 1-loop order. However, the best-fit bias parameters do not agree well with the peak-background split results. Being able to fit the void auto-power spectrum is particularly important not only because it is the direct observable in galaxy surveys, but also our method enables us to treat the bias parameters as nuisance parameters, which are sensitive to the techniques used to identify voids.

We investigate nonlocal Lagrangian bias contributions involving gradients of the linear density field, for which we have predictions from the excursion set peak formalism. We begin by writing down a bias expansion which includes all the bias terms, including the nonlocal ones. Having checked that the model furnishes a reasonable fit to the halo mass function, we develop a 1-point cross-correlation technique to measure bias factors associated with 2-distributed quantities. We validate the method with numerical realizations of peaks of Gaussian random fields before we apply it to N-body simulations. We focus on the lowest (quadratic) order nonlocal contributions. We can reproduce our measurement of \chi_{10} if we allow for an offset between the Lagrangian halo center-of-mass and the peak position. The sign and magnitude of \chi_{10} is consistent with Lagrangian haloes sitting near linear density maxima. The resulting contribution to the halo bias can safely be ignored for M = 10^13 Msun/h, but could become relevant at larger halo masses. For the second nonlocal bias \chi_{01} however, we measure a much larger magnitude than predicted by our model. We speculate that some of this discrepancy might originate from nonlocal Lagrangian contributions induced by nonspherical collapse.

Lagrangian displacement field $\Psi$ is the central object in Lagrangian perturbation theory (LPT). LPT is very successful at high redshifts, but it performs poorly at low redshifts due to severe shell crossing. To understand and quantify the effects of shell crossing, we extract $\Psi$ from N-body simulation and decompose it into scalar and vector parts. We find that at late time the power spectrum of the scalar part agrees with 1-loop results from LPT at large scales, while the power in small scales is much suppressed due to shell crossing. At z=0, the power spectrum of $\Psi$ is 10% lower than the 1-loop results at k = 0.1 h/Mpc. Shell crossing also generates the vector contribution in $\Psi$, although its effect is subdominant in comparison with the power suppression in the scalar part. At z=0, the vector part contributes 10% to the total power spectrum of $\Psi$ at k = 1 h/Mpc, while only 1% is expected from the vector contribution in LPT. We also examine the standard LPT recipes and some of its variants. In one of the variants, we include a power suppression factor in the displacement potential to take into account the power suppression in small scales after shell crossing. However, these simple phenomenological approaches are found to yield limited improvement compared to the standard LPT after the onset of shell crossing.

We explore the scale dependence of halo bias using real space cross-correlation measurements in N-body simulations and in Pinocchio, an algorithm based on Lagrangian Perturbation Theory. Recent work has shown how to interpret such real space measurements in terms of k-dependent bias in Fourier space, and how to remove the k-dependence to reconstruct the k-independent peak-background split halo bias parameters. We compare our reconstruction of the linear bias, which requires no free parameters, with previous estimates from N-body simulations which were obtained directly in Fourier space at large scales, and find very good agreement. Our reconstruction of the quadratic bias is similarly parameter-free, although in this case there are no previous Fourier space measurements to compare with. Our analysis of N-body simulations explicitly tests the predictions of the excursion set peaks (ESP) formalism of Paranjape et al. (2013) for the scale dependence of bias; we find that the ESP predictions accurately describe our measurements. In addition, our measurements in Pinocchio serve as a useful, successful consistency check between Pinocchio and N-body simulations that is not accessible to traditional measurements.

Halos are biased tracers of the dark matter distribution. It is often assumed that the patches from which halos formed are locally biased with respect to the initial fluctuation field, meaning that the halo-patch fluctuation field can be written as a Taylor series in that of the dark matter. If quantities other than the local density influence halo formation, then this Lagrangian bias will generically be nonlocal; the Taylor series must be performed with respect to these other variables as well. We illustrate the effect with Monte-Carlo simulations of a model in which halo formation depends on the local shear (the quadrupole of perturbation theory), and provide an analytic model which provides a good description of our results. Our model, which extends the excursion set approach to walks in more than one dimension, works both when steps in the walk are uncorrelated, as well as when there are correlations between steps. For walks with correlated steps, our model includes two distinct types of nonlocality: one is due to the fact that the initial density profile around a patch which is destined to form a halo must fall sufficiently steeply around it -- this introduces k-dependence to even the linear bias factor, but otherwise only affects the monopole of the clustering signal. The other is due to the surrounding shear field; this affects the quadratic and higher order bias factors, and introduces an angular dependence to the clustering signal. In both cases, our analysis shows that these nonlocal Lagrangian bias terms can be significant, particularly for massive halos; they must be accounted for in analyses of higher order clustering such as the halo bispectrum in Lagrangian or Eulerian space. Although we illustrate these effects using halos, our analysis and conclusions also apply to the other constituents of the cosmic web -- filaments, sheets and voids.

Future NASA X-ray spectroscopy missions will require high throughput, high resolution grating spectrometers. Off-plane reflection gratings are capable of meeting the performance requirements needed to realize the scientific goals of these missions. We have identified a novel grating fabrication method that utilizes common lithographic and microfabrication techniques to produce the high fidelity groove profile necessary to achieve this performance. Application of this process has produced an initial pre-master that exhibits a radial (variable line spacing along the groove dimension), high density (>6000 grooves/mm), laminar profile. This pre-master has been tested for diffraction efficiency at the BESSY II synchrotron light facility and diffracts up to 55% of incident light into usable spectral orders. Furthermore, tests of spectral resolving power show that these gratings are capable of obtaining resolutions well above 1300 ($\lambda/\Delta\lambda$) with limitations due to the test apparatus, not the gratings. Obtaining these results has provided confidence that this fabrication process is capable of producing off-plane reflection gratings for the next generation of X-ray observatories.

We develop a new test of local bias, by constructing a locally biased halo density field from sampling the dark matter-halo distribution. Our test differs from conventional tests in that it preserves the full scatter in the bias relation and it does not rely on perturbation theory. We put forward that bias parameters obtained using a smoothing scale R can only be applied to computing the halo power spectrum at scales k ~ 1/R. Our calculations can automatically include the running of bias parameters and give vanishingly small loop corrections at low-k. Our proposal results in much better agreement of the sampling and perturbation theory results with simulations. In particular, unlike the standard interpretation of local bias in the literature, our treatment of local bias does not generate a constant power in the low-k limit. We search for extra noise in the Poisson corrected halo power spectrum at wavenumbers below its turn-over and find no evidence of significant positive noise (as predicted by the standard interpretation) while we find evidence of negative noise coming from halo exclusion for very massive halos. Using perturbation theory and our non-perturbative sampling technique we also demonstrate that nonlocal bias effects discovered recently in simulations impact the power spectrum only at the few percent level in the weakly nonlinear regime.

We study the scale-dependence of halo bias in generic (non-local) primordial non-Gaussian (PNG) initial conditions of the type motivated by inflation, parametrized by an arbitrary quadratic kernel. We first show how to generate non-local PNG initial conditions with minimal overhead compared to local PNG models for a general class of primordial bispectra that can be written as linear combinations of separable templates. We run cosmological simulations for the local, and non-local equilateral and orthogonal models and present results on the scale-dependence of halo bias. We also derive a general formula for the Fourier-space bias using the peak-background split (PBS) in the context of the excursion set approach to halos and discuss the difference and similarities with the known corresponding result from local bias models. Our PBS bias formula generalizes previous results in the literature to include non-Markovian effects and non-universality of the mass function and are in better agreement with measurements in numerical simulations than previous results for a variety of halo masses, redshifts and halo definitions. We also derive for the first time quadratic bias results for arbitrary non-local PNG, and show that non-linear bias loops give small corrections at large-scales. The resulting well-behaved perturbation theory paves the way to constrain non-local PNG from measurements of the power spectrum and bispectrum in galaxy redshift surveys.

The relationship between galaxy and matter overdensities, bias, is most often assumed to be local. This is however unstable under time evolution, we provide proofs under several sets of assumptions. In the simplest model galaxies are created locally and linearly biased at a single time, and subsequently move with the matter (no velocity bias) conserving their comoving number density (no merging). We show that, after this formation time, the bias becomes unavoidably non-local and non-linear at large scales. We identify the non-local gravitationally induced fields in which the galaxy overdensity can be expanded, showing that they can be constructed out of the invariants of the deformation tensor (Galileons). In addition, we show that this result persists if we include an arbitrary evolution of the comoving number density of tracers. We then include velocity bias, and show that new contributions appear, a dipole field being the signature at second order. We test these predictions by studying the dependence of halo overdensities in cells of fixed matter density: measurements in simulations show that departures from the mean bias relation are strongly correlated with the non-local gravitationally induced fields identified by our formalism. The effects on non-local bias seen in the simulations are most important for the most biased halos, as expected from our predictions. The non-locality seen in the simulations is not fully captured by assuming local bias in Lagrangian space. Accounting for these effects when modeling galaxy bias is essential for correctly describing the dependence on triangle shape of the galaxy bispectrum, and hence constraining cosmological parameters and primordial non-Gaussianity. We show that using our formalism we remove an important systematic in the determination of bias parameters from the galaxy bispectrum, particularly for luminous galaxies. (abridged)

We use N-body simulations to study the nonlinear structure formation in brane-induced gravity, developing a new method that requires alternate use of Fast Fourier Transforms and relaxation. This enables us to compute the nonlinear matter power spectrum and bispectrum, the halo mass function, and the halo bias. From the simulation results, we confirm the expectations based on analytic arguments that the Vainshtein mechanism does operate as anticipated, with the density power spectrum approaching that of standard gravity within a modified background evolution in the nonlinear regime. The transition is very broad and there is no well defined Vainshtein scale, but roughly this corresponds to k_*~ 2 at redshift z=1 and k_*~ 1 at z=0. We checked that while extrinsic curvature fluctuations go nonlinear, and the dynamics of the brane-bending mode C receives important nonlinear corrections, this mode does get suppressed compared to density perturbations, effectively decoupling from the standard gravity sector. At the same time, there is no violation of the weak field limit for metric perturbations associated with C. We find good agreement between our measurements and the predictions for the nonlinear power spectrum presented in paper I, that rely on a renormalization of the linear spectrum due to nonlinearities in the modified gravity sector. A similar prediction for the mass function shows the right trends. Our simulations also confirm the induced change in the bispectrum configuration dependence predicted in paper I.

This paper reports an application of gas-kinetic BGK scheme to the computation of turbulent compressible convection in the stellar interior. After incorporating the Sub-grid Scale (SGS) turbulence model into the BGK scheme, we tested the effects of numerical parameters on the quantitative relationships among the thermodynamic variables, their fluctuations and correlations in a very deep, initially gravity-stratified stellar atmosphere. Comparison indicates that the thermal properties and dynamic properties are dominated by different aspects of numerical models separately. An adjustable Deardorff constant in the SGS model $c_\mu=0.25$ and an amplitude of artificial viscosity in the gas-kinetic BGK scheme $C_2=0$ are appropriate for current study. We also calculated the density-weighted auto- and cross-correlation functions in Xiong's (\cite{xiong77}) turbulent stellar convection theories based on which the gradient type of models of the non-local transport and the anisotropy of the turbulence are preliminarily studied. No universal relations or constant parameters were found for these models.

The constraints on the radion evolution in the Universal Extra Dimension (UED) model from Cosmic Microwave Background (CMB) and Type Ia supernovae (SNe Ia) data are studied. In the UED model, where both the gravity and standard model fields can propagate in the extra dimensions, the evolution of the extra dimensional volume, the radion, induces variation of fundamental constants. We discuss the effects of variation of the relevant constants in the context of UED for CMB power spectrum and SNe Ia data. We then use the three-year WMAP data to constrain the radion evolution at z \sim 1100, and the 2 \sigma constraint on \dot{\rho} / \rho_0 (\rho is a function of the radion, to be defined in the text) is [ -8.8, 6.6] \times 10 ^{-13} yr^-1. The SNe Ia gold sample yields a constraint on \dot{\rho} / \rho_0, for redshift between 0 and 1, to be [-4.7, 14] \times 10^{-13} yr^-1. Furthermore, the constraints from SNe Ia can be interpreted as bounds on the evolution QCD scale parameter, \dot{\Lambda}_{QCD} / \Lambda_{QCD, 0}, [-1.4, 2.8] \times 10^{-11} yr^-1, without reference to the UED model.

We perform Large eddy simulations of turbulent compressible convection in stellar-type convection zones by solving the Navi\'{e}r-Stokes equations in three dimensions. We estimate the extent of penetration into the stable layer above a stellar-type convection zone by varying the rotation rate ({\boldmath$\rm\Omega$}), the inclination of the rotation vector ($\theta$) and the relative stability ($S$) of the upper stable layer. The computational domain is a rectangular box in an f-plane configuration and is divided into two regions of unstable and stable stratification with the stable layer placed above the convectively unstable layer. Several models have been computed and the penetration distance into the stable layer above the convection zone is estimated by determining the position where time averaged kinetic energy flux has the first zero in the upper stable layer. The vertical grid spacing in all the model is non-uniform, and is less in the upper region so that the flows are better resolved in the region of interest. We find that the penetration distance increases as the rotation rate increases for the case when the rotation vector is aligned with the vertical axis. However, with the increase in the stability of the upper stable layer, the upward penetration distance decreases. Since we are not able to afford computations with finer resolution for all the models, we compute a number of models to see the effect of increased resolution on the upward penetration. In addition, we estimate the upper limit on the upward convective penetration from stellar convective cores.

In many versions of brane model, the modulus field of extra dimensions, the radion, could have cosmological evolution, which induces variation of the Higgs vacuum expectation value, $<H>$, resulting in cosmological variation of the electron mass $m_e$. The formation of Cosmic Microwave Background (CMB) anisotropies is thus affected, causing changes both in the peaks positions and amplitudes in the CMB power spectra. Using the three-year Wilkinson Microwave Anisotropies Probe (WMAP) CMB data, with the Hubble parameter $H_0$ fixed to be the Hubble Space Telescope (HST) result 72 km s$^{-1}$ Mpc$^{-1}$, we obtain a constraint on $\rho$, the ratio of the value of $<H>$ at CMB recombination to its present value, to be [0.97, 1.02].

In this work we derive the covariant and gauge invariant perturbation equations in general theories of $f(R)$ gravity in the Palatini formalism to linear order and calculate the cosmic microwave background (CMB) and matter power spectra for an extensively discussed model, $f(R) = R + \alpha(-R)^\beta$, which is a possible candidate for the late-time cosmic accelerating expansion found recently. These spectra are discussed and found to be sensitively dependent on the value of $\beta$. We are thus able to make stringent constraints on $\beta$ from cosmological data on CMB and matter power spectra: The three-year Wilkinson Microwave Anisotropy Probe (WMAP) data alone gives a constraint $|\beta| \lesssim \mathcal{O}(10^{{-3}})$ while the joint WMAP, Supernova Lagacy Survey (SNLS) and Sloan Digital Sky Survey (SDSS) data sets tightens this to $\beta \sim \mathcal{O}(10^{{-6}})$, about an order of magnitude more stringent than the constraint from SDSS data alone, which makes this model practically indistinguishable from the standard $\Lambda\mathrm{CDM}$ paradigm.

We use the Cosmic Microwave Background Anisotropies (CMBA) power spectra to constrain the cosmological variation of gravitational constant G. It is found that the sensitivity of CMBA to the variation of G is enhanced when G is required to converge to its present value. The variations of G from the CMB decoupling epoch z ~ 1000 to the present time are modelled by a step function and a linear function of scale factor $a$ respectively, and the corresponding 95% confidence intervals for G/G_0 are [0.95, 1.05] and [0.89, 1.13], G_0 being the present value. The CMBA constraint is unique in the sense that it entails the range of redshift from z \approx 1000 to 0.

CH Cygni is a symbiotic star consisting of an M giant and an accreting white dwarf, which is known to be a highly variable X-ray source with a complex, two-component, spectra. Here we report on two Suzaku observations of CH Cyg, taken in 2006 January and May, during which the system was seen to be in a soft X-ray bright, hard X-ray faint state. Based on the extraordinary strength of the 6.4 keV fluorescent Fe K-alpha line, we show that the hard X-rays observed with Suzaku are dominated by scattering.

CHICOS (California HIgh school Cosmic ray ObServatory) is presently an array of more than 140 detectors distributed over a large area (~400 km^2) of southern California, and will consist of 180 detectors at 90 locations in the near future. These sites, located at area schools, are equipped with computerized data acquisition and automatic nightly data transfer (via internet) to our Caltech lab. The installed sites make up the largest currently operating ground array for ultra-high energy cosmic ray research in the northern hemisphere. The goal of CHICOS is to provide data related to the flux and distribution of arrival directions for ultra-high energy cosmic rays. We have performed detailed Monte-Carlo calculations to determine the density and arrival-time distribution of charged particles in extensive air showers for the CHICOS array. Calculations were performed for proton primaries with energies 10^18 to 10^21 eV and zenith angles out to 50 degrees. We have developed novel parameterizations for both distributions as functions of distance from the shower axis, primary energy, and incident zenith angle. These parameterizations are used in aperture calculations and reconstruction of shower data, enabling preliminary analysis of ultra-high energy shower data from CHICOS.

Numerical simulations of turbulent stratified convection are used to study models with approximately the same convective flux, but different radiative fluxes. As the radiative flux is decreased, for constant convective flux: the entropy jump at the top of the convection zone becomes steeper, the temperature fluctuations increase and the velocity fluctuations decrease in magnitude, and the distance that low entropy fluid from the surface can penetrate increases. Velocity and temperature fluctuations follow mixing length scaling laws.

Compared to the Sun, the atmospheric structure and convective flow in Procyon A exhibit the following characteristics: (1) the highly superadiabatic transition layer (SAL) is located at much shallower optical depth; it is in a dynamically active region, and its outer region is located part of the time in the optically thin atmosphere; (2) the outer region of the SAL moves from an optically thin region to thick region and back again over a time of 20-30 minutes. This motion, which is driven by the granulation, takes place in a time approximately half the turnover time of the largest granules; The main reason for the radically different radiative-convective behaviour in Procyon A compared to the Sun is the role played by turbulent eddies in determining the overall flow/thermal structure. The turbulent pressure and turbulent kinetic energy can exceed 50 % of the local gas pressure (compared to about 10-20 % in the Sun). The Procyon A simulation thus reveals two distinct timescales - the autocorrelation time of the vertical velocity and the characteristic timescale of the SAL which is tied to granulation. Just below the surface the autocorrelation decay time is about 5 minutes in Procyon A, and the SAL motion timescale is 20-30 mins. When the SAL penetrates the optically thin region there are efficient radiative losses and the peak of the SAL is low. We speculate that these losses damp out the relative amplitudes in luminosity (temperature fluctuations) compared to velocity (Doppler). Although this will not affect the frequencies of the peaks in the power spectrum, it will probably lower the average amplitude of the peaks relative to the noise background.

This paper describes a series of 3D simulations of shallow inefficient convection in the outer layers of the Sun. The computational domain is a closed box containing the convection-radiation transition layer, located at the top of the solar convection zone. The most salient features of the simulations are that: i)The position of the lower boundary can have a major effect on the characteristics of solar surface convection (thermal structure, kinetic energy and turbulent pressure). ii)The width of the box has only a minor effect on the thermal structure, but a more significant effect on the dynamics (rms velocities). iii)Between the surface and a depth of 1 Mm, even though the density and pressure increase by an order of magnitude, the vertical correlation length of vertical velocity is always close to 600 km. iv) In this region the vertical velocity cannot be scaled by the pressure or the density scale height. This casts doubt on the applicability of the mixing length theory, not only in the superadiabatic layer, but also in the adjacent underlying layers. v) The final statistically steady state is not strictly dependent on the initial atmospheric stratification.

We measure the eddy viscosity in the outermost layers of the solar convection zone by comparing the rotation law computed with the Reynolds stress resulting from f-plane simulations of the angular momentum transport in rotating convection with the observed differential rotation pattern. The simulations lead to a negative vertical and a positive horizontal angular momentum transport. The consequence is a subrotation of the outermost layers, as it is indeed indicated both by helioseismology and the observed rotation rates of sunspots. In order to reproduce the observed gradient of the rotation rate a value of about 1.5 x 10^{13} cm/s for the eddy viscosity is necessary. Comparison with the magnetic eddy diffusivity derived from the sunspot decay yields a surprisingly large magnetic Prandtl number of 150 for the supergranulation layer. The negative gradient of the rotation rate also drives a surface meridional flow towards the poles, in agreement with the results from Doppler measurements. The successful reproduction of the abnormally positive horizontal cross correlation (on the northern hemisphere) observed for bipolar groups then provides an independent test for the resulting eddy viscosity.

Far-infrared spectra (43-197um) of 34 nearby galaxies obtained by the Long Wavelength Spectrometer (LWS) aboard the Infrared Space Observatory (ISO) were analyzed to investigate the general properties of interstellar matter in galaxies. The line fluxes of [CII]158um and [NII]122um relative to the total far-infrared flux (FIR) decrease as the far-infrared color becomes bluer, while the ratio of the [OI]63um flux to FIR does not show a systematic trend with the color. The [OIII]88um to FIR ratio shows a large scatter with a weak trend of increase with the color. We estimate the physical conditions of photodissociation regions (PDRs) in the sample galaxies, such as the far-ultraviolet radiation field intensity Go and the gas density n by assuming that all the observed [OI]63um and far-infrared continuum emissions come from PDRs. The present analysis suggests that the decrease in [CII]158um/FIR with the far-infrared color may not be accounted for by the decrease in the photoelectric heating efficiency owing to the increase in positive charges of dust grains because a measure of the efficiency, Go/n, is found to stay constant with the far-infrared color. Instead the decrease can be interpreted in terms of either the increase in the collisional de-excitation of the [CII] transition due to the increase in the gas density or the decrease in the ionized component relative to the far-infrared intensity suggested by the decrease in [NII]122um/FIR. Based on the present analysis, we derive average relations of the far-infrared color with Go and n in galaxies, which can be applied to the investigation of interstellar matter in distant galaxies.

We present results of two simulations of the convection zone, obtained by solving the full hydrodynamic equations in a section of a spherical shell. The first simulation has cylindrical rotation contours (parallel to the rotation axis) and a strong meridional circulation, which traverses the entire depth. The second simulation has isorotation contours about mid-way between cylinders and cones, and a weak meridional circulation, concentrated in the uppermost part of the shell. We show that the solar differential rotation is directly related to a latitudinal entropy gradient, which pervades into the deep layers of the convection zone. We also offer an explanation of the angular velocity shear found at low latitudes near the top. A non-zero correlation between radial and zonal velocity fluctuations produces a significant Reynolds stress in that region. This constitutes a net transport of angular momentum inwards, which causes a slight modification of the overall structure of the differential rotation near the top. In essence, the {\it thermodynamics controls the dynamics through the Taylor-Proudman momentum balance}. The Reynolds stresses only become significant in the surface layers, where they generate a weak meridional circulation and an angular velocity `bump'.

We present results of three dimensional simulations of the uppermost part of the sun, at 3 stages of its evolution. Each model includes physically realistic radiative-hydrodynamics (the Eddington approximation is used in the optically thin region), varying opacities and a realistic equation of state (full treatment of the ionization of H and He). In each evolution model, we investigate a domain, which starts at the top of the photosphere and ends just inside the convection zone (about 2400 km in the sun model). This includes all of the super-adiabatic layer (SAL). Due to the different positions of the three models in the $log (g) $ vs $log T_{eff}$ plane, the more evolved models have lower density atmospheres. The reduction in density causes the amount of overshoot into the radiation layer, to be greater in the more evolved models.

We have constructed models for the sun at three stages of its evolution: a zero-age main sequence model, the present sun, and a subgiant model. For each model, the turbulent pressure and turbulent kinetic energy were calculated from 3-d radiative hydrodynamical simulations (described in the poster by Robinson et al.), and inserted into the 1-d stellar models. We note that in these simulations, the turbulent pressure is not a free parameter, but can be computed from the resulting velocity field. We show the calculated p-mode frequencies for the model of the present sun, with and without turbulent pressure, and compare them to the observed solar frequencies. When the turbulent pressure is included in the models, the calculated frequencies are brought closer to the observed frequencies in the sun by up to two $\mu Hz$, strictly from structural effects. The effect of including turbulent pressure on p-mode frequencies is also shown for the zero-age main sequence model. Our models also suggest that the importance of turbulent pressure increases as the star evolves into the subgiant region. We discuss the importance of also including realistic turbulence as well as radiation in the non-adiabatic calculation of oscillation frequencies.

Using the Mid-Infrared Spectrometer on board the Infrared Telescope in Space and the low-resolution grating spectrometer (PHT-S) on board the Infrared Space Observatory, we obtained 820 mid-infrared (5 to 12 $\mu$m) spectra of the diffuse interstellar medium (DIM) in the Galactic center, W51, and Carina Nebula regions. These spectra indicate that the emission is dominated by the unidentified infrared (UIR) emission bands at 6.2, 7.7, 8.6, and 11.2 $\mu$m. The relative band intensities (6.2/7.7 $\mu$m, 8.6/7.7 $\mu$m, and 11.2/7.7 $\mu$m) were derived from these spectra, and no systematic variation in these ratios was found in our observed regions, in spite of the fact that the incident radiation intensity differs by a factor of 1500. Comparing our results with the polycyclic aromatic hydrocarbons (PAHs) model for the UIR band carriers, PAHs in the DIM have no systematic variation in their size distribution, their degree of dehydrogenation is independent of the strength of UV radiation field, and they are mostly ionized. The latter finding is incompatible with past theoretical studies, in which a large fraction of neutral PAHs is predicted in this kind of environment. A plausible resolution of this discrepancy is that the recombination coefficients for electron and large PAH positive ion are by at least an order of magnitude less than those adopted in past theoretical studies. Because of the very low population of neutral state molecules, photoelectric emission from interstellar PAHs is probably not the dominant source of heating of the diffuse interstellar gas. The present results imply constant physical and chemical properties of the carriers of the UIR emission bands in the DIM.

We report the detection of a broad 22 microns emission feature in the Carina nebula H II region by the Infrared Space Observatory (ISO) Short Wavelength Spectrometer. The feature shape is similar to that of the 22 microns emission feature of newly synthesized dust observed in the Cassiopeia A supernova remnant. This finding suggests that both of the features are arising from the same carrier, and that supernovae are probably the dominant production source of this new interstellar grain. A similar broad emission dust feature is also found in the spectra of two starburst galaxies from the ISO archival data. This new dust grain could be an abundant component of interstellar grains and can be used to trace the supernova rate or star formation rate in external galaxies. The existence of the broad 22 microns emission feature complicates the dust model for starburst galaxies and must be taken into account correctly in the derivation of dust color temperature. Mg protosilicate has been suggested as the carrier of the 22 microns emission dust feature observed in Cassiopeia A. The present results provide useful information in studies on chemical composition and emission mechanism of the carrier.

Using the Mid-Infrared Spectrometer (MIRS) on board the Infrared Telescope in Space (IRTS) we obtained the 4.5 to 11.7 micron spectra of the stellar populations and diffuse interstellar medium in the Galactic bulge (l = 8.7 degree, b = 2.9, 4.0, 4.7, and 5.7 degree). Below galactic latitudes of 4.0 degree, the mid-infrared background spectra in the bulge are similar to the spectra of M and K giants. The UIR emission bands (6.2, 7.7, 8.6, and 11.3 micron) are also detected in these regions and likely arise from the diffuse interstellar medium in the disk. Above galactic latitudes of 4.0 degree, the mid-infrared background spectra are similar to the spectra of those oxygen-rich evolved stars with high mass-loss rates detected by IRAS. One likely interpretation is that this background emission arises predominantly from those stars with very low luminosities that have not been detected by IRAS. The age for such low-luminosity evolved stars could be 15 Gyr, and the existence of a large number of evolved stars with high mass-loss rates in the bulge has a significant impact on our understanding of the stellar content in the Galactic bulge.