We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

We look at the contribution of different galaxy populations to the atomic hydrogen (HI) mass function (HIMF) and the HI density parameter, $\Omega_{\text{HI}}$, in the local Universe. Our analysis is based on a sample of 7857 HI-selected galaxies selected from a volume common to the SDSS and ALFALFA surveys (40$\%$ catalog -- $\alpha.40$). We define different populations of galaxies in the color(u-r)-magnitude($M_{\text{r}}$) plane and compute the HIMF for each of them. Additionally we compute the HIMF for dark galaxies; these are undetected in SDSS and represent $\sim 2\%$ of the total sample. We find that the luminous red population dominates the total HIMF for $\log_{10}(M_{\text{HI}}h^2_{70}/M_{\odot}) \geq 10.4$. The full red population -- luminous and faint -- represents about $\sim 17\%$ of the $\Omega_{\text{HI}}$ budget, while that of the dark population is $\sim 3\%$. The HIMF about the knee, $\log_{10}(M_{\text{HI}}h^2_{70}/M_{\odot}) \in [8,10.4]$, is dominated by the faint and luminous blue populations, the latter dominating at larger masses in this interval. Their total contribution to $\Omega_{\text{HI}}$ is $\sim 55-70\%$, the variation depending on the definition of population. The dominant populations at the low mass end, $\log_{10}(M_{\text{HI}}h^2_{70}/M_{\odot}) \leq 8.0$ are the faint blue and faint bluer populations, the latter's dominance being sensitive to its definition. The full blue (blue--bluer luminous and faint) population represents $\sim 80\%$ of $\Omega_{\text{HI}}$. A bimodal HIMF suggested by our results is however not seen since the amplitude of the HIMF of the luminous red population is small compared to that of the luminous blue population.

We present here identification and characterization of the young stellar population associated with an active star-forming site Sh2-242. We used our own new optical imaging and spectroscopic observational data, as well as several archival catalogs, e.g., Pan-STARRS 1, $Gaia$ DR2, IPHAS, WIRCam, 2MASS, and $Spitzer$. Slit spectroscopic results confirm the classification of the main ionizing source BD+26 980 as an early-type star of spectral type B0.5 V. The spectrophotometric distance of the star is estimated as 2.08 $\pm$ 0.24 kpc, which confirms the source as a member of the cluster. An extinction map covering a large area (diameter $\sim$ 50') is generated with $H$ and $K$ photometry toward the region. From the map, three distinct locations of peak extinction complexes ($A_{V}$ $\simeq$ 7$-$17 mag) are identified for the very first time. Using the infrared color excess, a total of 33 Class I and 137 Class II young objects are classified within the region. The IPHAS photometry reveals classification of 36 H$\alpha$ emitting sources, which might be class II objects. Among 36 H$\alpha$ emitting sources, 5 are already identified using infrared excess emission. In total, 201 young objects are classified toward S242 from this study. The membership status of the young sources is further windowed with the inclusion of parallax from the $Gaia$ DR2 catalog. Using the optical and infrared color-magnitude diagrams, the young stellar objects are characterized with an average age of $\sim$ 1 Myr and the masses in the range 0.1$-$3.0 $M_\odot$. The census of the stellar content within the region is discussed using combined photometric and spectroscopic data.

We present optical $I$-band light curves of the stars towards a star-forming region Cygnus OB7 from 17 nights photometric observations. The light curves are generated from a total of 381 image frames with very good photometric precision. From the light curves of 1900 stars and their periodogram analyses, we detect 31 candidate variables including five previously identified. 14 out of 31 objects are periodic and exhibit the rotation rates in the range of 0.15 to 11.60 days. We characterize those candidate variables using optical/infrared colour-colour and colour-magnitude diagrams. From spectral indices of the candidate variables, it turns out that four are probably Classical T-Tauri stars (CTTSs), rest remain unclassified from present data, they are possibly field stars or diskless pre-main sequence stars towards the region. Based on their location on the various colour-magnitude diagrams, the ages of two T Tauri Stars were estimated to be $\sim$ 5 Myr. The light curves indicate at least five of the periodic variables are eclipsing systems. The spatial distribution of young variable candidates on Planck 857 GHz (350 $\mu$m) and 2MASS $K_s$ images suggest that at least two of the CTTSs are part of the active star-forming cloud Lynds 1003.

The {\it Planck} Galactic Cold Clumps (PGCCs) are the possible representations of the initial conditions and the very early stages of star formation. With an objective to understand better the star and star cluster formation, we probe the molecular cloud associated with PGCC G108.37-01.06 (hereafter, PG108.3), which can be traced in a velocity range $-$57 to $-$51 km s$^{-1}$. The IPHAS images reveal H$\alpha$ emission at various locations around PG108.3, and optical spectroscopy of the bright sources in those zones of H$\alpha$ emission disclose two massive ionizing sources with spectral type O8-O9V and B1V. Using the radio continuum, we estimate ionizing gas parameters and find the dynamical ages of \hii regions associated with the massive stars in the range 0.5$-$0.75 Myr. Based on the stellar surface density map constructed from the deep near-infrared CHFT observations, we find two prominent star clusters in PG108.3; of which, the cluster associated with \hii region S148 is moderately massive ($\sim$ 240 M$\sun$). A careful inspection of JCMT $^{13}$CO(3$-$2) molecular data exhibits that the massive cluster is associated with a number of filamentary structures. Several embedded young stellar objects (YSOs) are also identified in the PG108.3 along the length and junction of filaments. We find the evidence of velocity gradient along the length of the filaments. Along with kinematics of the filaments and the distribution of ionized, molecular gas and YSOs, we suggest that the cluster formation is most likely due to the longitudinal collapse of the most massive filament in PG108.3.

We present here the recent results on two distant Galactic H II regions, namely NGC~2282 and Sh2-149, obtained with multiwavelength observations. Our optical spectroscopic analysis of the bright sources have been used to identify the massive members, and to derive the fundamental parameters such as age and distance of these regions. Using IR colour-colour criteria and H$_\alpha$-emission properties, we have identified and classified the candidate young stellar objects (YSOs) in these regions. The $^{12}$CO(1-0) continuum maps along with the {\it K}-band extinction maps, and spatial distribution of YSOs are used to investigate the structure and morphology of the molecular cloud associated with these H II regions. Overall analysis of these regions suggests that the star formation occurs at the locations of the denser gas, and we also find possible evidence of the induced star formation due to the feedback from massive stars to its surrounding molecular medium.

We describe the time-dependent properties of a new spectroscopically confirmed Mira variable, which was discovered in 2013 as MASTER-Net Optical Transient (OT) J212444.87+321738.3 towards the Cygnus constellation. We have performed long-term optical/near-infrared (NIR) photometric and spectroscopic observations to characterize the object. From the optical/NIR light curves, we estimate a variability period of 465 $\pm$ 30 days. The wavelength-dependent amplitudes of the observed light-curves range from $\Delta$I$\sim$4 mag to $\Delta$K$\sim$1.5 mag. The (J-K) color-index varies from 1.78 to 2.62 mag over phases. Interestingly, a phase lag of $\sim$60 days between optical and NIR light curves is also seen, as in other Miras. Our optical/NIR spectra show molecular features of TiO, VO, CO, and strong water bands which are a typical signature of oxygen-rich Mira. We rule out S- or C-type as ZrO bands at 1.03 and 1.06 $\mu$m and $C_2$ band at 1.77 $\mu$m are absent. We estimate the effective temperature of the object from the SED, and distance and luminosity from standard Period-Luminosity relations. The optical/NIR spectra display time-dependent atomic and molecular features (e.g. TiO, NaI, CaI, H$_2$O,CO), as commonly observed in Miras. Such spectroscopic observations are useful for studying pulsation variability in Miras.

We report here CCD $I$-band time-series photometry of a young (2$-$5 Myr) cluster NGC 2282 to identify and understand the variability of pre-main-sequence (PMS) stars. The $I$-band photometry, down to $\sim$ 20.5 mag, enables us to probe the variability towards the lower mass end ($\sim$ 0.1 M$_\odot$) of the PMS stars. From the light curves of 1627 stars, we identified 62 new photometric variable candidates. Their association with the region was established from H$\alpha$ emission and infrared (IR) excess. Among 62 variables, 30 young variables exhibit H$\alpha$ emission, near-IR (NIR)/mid-IR (MIR) excess or both, and they are candidate members of the cluster. Out of 62 variables, 41 are periodic variables with the rotation rate ranging from 0.2 to 7 days. The period distribution exhibits a median period at $\sim$ 1-day as in many young clusters (e.g., NGC~2264, ONC, etc.), but it follows a uni-modal distribution unlike others having bimodality with the slow rotators peaking at $\sim$ 6$-$8 days. To investigate the rotation-disk and variability-disk connection, we derived NIR excess from $\Delta$(I$-$K) and MIR excess from $Spitzer$ [3.6]$-$[4.5] $\mu$m data. No conclusive evidence of slow rotation with the presence of disks around stars and fast rotation for diskless stars is obtained from our periodic variables. A clear increasing trend of the variability amplitude with the IR excess is found for all variables.

This paper deals with the viscous accretion flow of modified Chaplygin gas towards a black hole as the central gravitating object. Modified Chaplygin gas is a particular type of dark energy model which mimics of radiation era to phantom era depending on the different values of its parameters. We compare the drak energy accretion with the flow of adiabatic gas. An accretion disc flow around a black hole is an example of a transonic flow. To make the model, we consider three components of Navier Stokes' equation, the equation of continuity and the modified Chaplygin gas's equation of state. As a transonic flow passes through the sonic point, the velocity gradient being apparently singular there, gives rise to two flow branches : one infalling, the accretion and the other outgoing, the wind. We show that the wind curve is stronger and wind speed reaches to that of light at a finite distance from the black hole when dark energy is considered. Besides, if we increase the viscosity, accretion disc is being shortened. This two process acting together deviates much from adiabatic accretion case. It shows a weakening process for the accretion procedure by the works of viscous system influenced angular momentum transport and the repulsive force of the modified Chaplygin gas together.

The neutron capture cross sections have been studied near the $Z=50$ closed shell for a number of nuclei those take part in heavy element nucleosynthesis, the slow and the rapid neutron capture processes and the proton capture process. An optical model potential is constructed in theoretical approach by folding the density dependent M3y purely real nucleon-nucleon interaction with the target radial matter density in relativistic mean field (RMF) approach. The standard code TALYS1.8 is used for cross-section calculation. We have presented the Maxwellian-averaged capture cross-section (MACS) values and stellar neutron capture reaction rates at astrophysically relevant thermal energies and temperatures.

Proton capture cross sections in the energy range of astrophysical interest for mass region 40-54 have been calculated in the Hauser-Feshbach formalism with reaction code TALYS1.6. The density dependent M3Y effective nucleon-nucleon interaction folded with target radial matter densities from relativistic mean field approach is used to obtain the semi-microscopic optical potential. A definite normalization of potential-well depths has been used over the entire mass region. The $(p,\gamma)$ rates of some reactions, important in the astrophysical scenario, are calculated using the potential in the relevant mass region.

Radiative thermal neutron capture cross sections for nuclei participating in s-process and p-process nucleosynthesis in and around $N=50$ closed neutron shell have been calculated in statistical semi-microscopic Hauser-Feshbach approach for the energy range of astrophysical interest. A folded optical model potential is constructed utilizing the standard DDM3Y real nucleon-nucleon interaction. The folding of the interaction with target radial matter densities, obtained from the relativistic-mean-field approach, is done in coordinate space using the spherical approximation. The standard nuclear reaction code TALYS1.8 is used for cross-section calculation. The cross sections are compared with experimental results and reasonable agreements are found for almost all cases. Maxwellian-averaged cross sections (MACS) for the nuclei are presented at a single thermal energy of 30 keV relevant to s-process. We have also presented the MACS values over a range of energy from 5 to 100 keV for neutron magic nuclei with $(N=50)$.

We present the analysis of the stellar content of NGC~2282, a young cluster in the Monoceros constellation, using deep optical $BVI$ and IPHAS photometry along with infrared (IR) data from UKIDSS and $Spitzer$-IRAC. Based on the stellar surface density analysis using nearest neighborhood method, the radius of the cluster is estimated as $\sim$ 3.15$\arcmin$. From optical spectroscopic analysis of 8 bright sources, we have classified three early B-type members in the cluster, which includes, HD 289120, a previously known B2V type star, a Herbig Ae/Be star (B0.5 Ve) and a B5 V star. From spectrophotometric analyses, the distance to the cluster has been estimated as $\sim$ 1.65 kpc. The $K$-band extinction map is estimated using nearest neighborhood technique, and the mean extinction within the cluster area is found to be A$_V$ $\sim$ 3.9 mag. Using IR colour-colour criteria and H$_\alpha$-emission properties, we have identified a total of 152 candidate young stellar objects (YSOs) in the region, of which, 75 are classified as Class II, 9 are Class I YSOs. Our YSO catalog also includes 50 H$_\alpha$-emission line sources, identified using slitless spectroscopy and IPHAS photometry data. Based on the optical and near-IR colour-magnitude diagram analyses, the cluster age has been estimated to be in the range of 2 $-$ 5 Myr, which is in agreement with the estimated age from disc fraction ($\sim$ 58\%). Masses of these YSOs are found to be $\sim$ 0.1$-$2.0 M$_\odot$. Spatial distribution of the candidate YSOs shows spherical morphology, more or less similar to the surface density map.

We perform a combined perturbation and observational investigation of the scenario of non-minimal derivative coupling between a scalar field and curvature. First we extract the necessary condition that ensures the absence of instabilities, which is fulfilled more sufficiently for smaller coupling values. Then using Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB) observations, we show that, contrary to its significant effects on inflation, the non-minimal derivative coupling term has a negligible effect on the universe acceleration, since it is driven solely by the usual scalar-field potential. Therefore, the scenario can provide a unified picture of early and late time cosmology, with the non-minimal derivative coupling term responsible for inflation, and the usual potential responsible for late-time acceleration. Additionally, the fact that the necessary coupling term does not need to be large, improves the model behavior against instabilities.

We investigate the cosmic evolution of the linear bias in the framework of a flat FLRW spacetime. We consider metric perturbations in the Newtonian gauge, including Hubble scale effects. Making the following assumptions, (i) scale independent current epoch bias $b_0$, (ii) equal accelerations between tracers and matter, (iii) unimportant halo merging effects (which is quite accurate for $z<3$), we analytically derive the scale dependent bias evolution. The identified scale dependence is only due to Hubble scale evolution GR effects, while other scale dependence contributions are ignored. We find that up to galaxy cluster scales the fluctuations of the metric do not introduce a significant scale dependence in the linear bias. Our bias evolution model is then used to derive a connection between the matter growth index $\gamma$ and the observable value of the tracer power spectrum normalization $\sigma_8(z)$. We show how this connection can be used as an observational test of General Relativity on extragalactic scales.

We investigate cosmological scenarios of generalized Chaplygin gas in a universe governed by Horava-Lifshitz gravity. We consider both the detailed and non-detailed balance versions of the gravitational background, and we include the baryonic matter and radiation sectors. We use observational data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis (BBN), to constrain the parameters of the model, and we provide the corresponding likelihood contours. We deduce that the present scenario is compatible with observations. Additionally, examining the evolution of the total equation-of-state parameter, we find in a unified way the succession of the radiation, matter, and dark energy epochs, consistently with the thermal history of the universe.

We show that the f(T) gravitational paradigm, in which gravity is described by an arbitrary function of the torsion scalar, can provide a mechanism for realizing bouncing cosmologies, thereby avoiding the Big Bang singularity. After constructing the simplest version of an f(T) matter bounce, we investigate the scalar and tensor modes of cosmological perturbations. Our results show that metric perturbations in the scalar sector lead to a background-dependent sound speed, which is a distinguishable feature from Einstein gravity. Additionally, we obtain a scale-invariant primordial power spectrum, which is consistent with cosmological observations, but suffers from the problem of a large tensor-to-scalar ratio. However, this can be avoided by introducing extra fields, such as a matter bounce curvaton.

We examine the evolution of quintessence models with potentials satisfying (V'/V)^2<<1 and V"/V<<1, in the case where the initial field velocity is nonzero. We derive an analytic approximation for the evolution of the equation of state parameter, w, for the quintessence field. We show that such models are characterized by an initial rapid freezing phase, in which the equation of state parameter w decreases with time, followed by slow thawing evolution, for which w increases with time. These models resemble constant-V models at early times but diverge at late times. Our analytic approximation gives results in excellent agreement with exact numerical evolution.

We investigate f(T) cosmology in both the background, as well as in the perturbation level, and we present the general formalism for reconstructing the equivalent one-parameter family of f(T) models for any given dynamical dark energy scenario. Despite the completely indistinguishable background behavior, the perturbations break this degeneracy and the growth histories of all these models differ from one another. As an application we reconstruct the f(T) equivalent for quintessence, and we show that the deviation of the matter overdensity evolution is strong for small scales and weak for large scales, while it is negligible for large redshifts.

We investigate the cosmological perturbations in f(T) gravity. Examining the pure gravitational perturbations in the scalar sector using a diagonal vierbien, we extract the corresponding dispersion relation, which provides a constraint on the f(T) ansatzes that lead to a theory free of instabilities. Additionally, upon inclusion of the matter perturbations, we derive the fully perturbed equations of motion, and we study the growth of matter overdensities. We show that f(T) gravity with f(T) constant coincides with General Relativity, both at the background as well as at the first-order perturbation level. Applying our formalism to the power-law model we find that on large subhorizon scales (O(100 Mpc) or larger), the evolution of matter overdensity will differ from LCDM cosmology. Finally, examining the linear perturbations of the vector and tensor sectors, we find that (for the standard choice of vierbein) f(T) gravity is free of massive gravitons.

We propose a scale dependent analytic approximation to the exact linear growth of density perturbations in Scalar-Tensor (ST) cosmologies. In particular, we show that on large subhorizon scales, in the Newtonian gauge, the usual scale independent subhorizon growth equation does not describe the growth of perturbations accurately, as a result of scale-dependent relativistic corrections to the Poisson equation. A comparison with exact linear numerical analysis indicates that our approximation is a significant improvement over the standard subhorizon scale independent result on large subhorizon scales. A comparison with the corresponding results in the Synchronous gauge demonstrates the validity and consistency of our analysis.

Models that lead to a cosmological stiff fluid component, with a density $\rho_S$ that scales as $a^{-6}$, where $a$ is the scale factor, have been proposed recently in a variety of contexts. We calculate numerically the effect of such a stiff fluid on the primordial element abundances. Because the stiff fluid energy density decreases with the scale factor more rapidly than radiation, it produces a relatively larger change in the primordial helium-4 abundance than in the other element abundances, relative to the changes produced by an additional radiation component. We show that the helium-4 abundance varies linearly with the density of the stiff fluid at a fixed fiducial temperature. Taking $\rho_{S10}$ and $\rho_{R10}$ to be the stiff fluid energy density and the standard density in relativistic particles, respectively, at $T = 10$ MeV, we find that the change in the primordial helium abundance is well-fit by $\Delta Y_p = 0.00024(\rho_{S10}/\rho_{R10})$. The changes in the helium-4 abundance produced by additional radiation or by a stiff fluid are identical when these two components have equal density at a "pivot temperature", $T_*$, where we find $T_* = 0.55$ MeV. Current estimates of the primordial $^4$He abundance give the constraint on a stiff fluid energy density of $\rho_{S10}/\rho_{R10} < 30$.

A $\Lambda$CDM model with dark matter that decays into inert relativistic energy on a timescale longer than the Hubble time will produce an expansion history that can be misinterpreted as stable dark matter with time-varying dark energy. We calculate the corresponding spurious equation of state parameter, $\widetilde w_\phi$, as a function of redshift, and show that the evolution of $\widetilde w_\phi$ depends strongly on the assumed value of the dark matter density, erroneously taken to scale as $a^{-3}$. Depending on the latter, one can obtain models that mimic quintessence ($\widetilde w_\phi > -1$), phantom models ($\widetilde w_\phi < -1$) or models in which the equation of state parameter crosses the phantom divide, evolving from $\widetilde w_\phi > -1$ at high redshift to $\widetilde w_\phi < -1$ at low redshift. All of these models generically converge toward $\widetilde w_\phi \approx -1$ at the present. The degeneracy between the $\Lambda$CDM model with decaying dark matter and the corresponding spurious quintessence model is broken by the growth of density perturbations.

We propose a class of simple dark energy models which predict a late-time dark radiation component and a distinctive time-dependent equation of state $w(z)$ for redshift $z < 3$. The dark energy field can be coupled strongly enough to Standard Model particles to be detected in colliders, and the model requires only modest additional particle content and little or no fine-tuning other than a new energy scale of order milli-electron volts.

We use observational data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis (BBN), to constrain the running parameter $\lambda$ of Horava-Lifshitz gravity, which determines the flow between the Ultra-Violet and the Infra-Red. We consider both the detailed and non-detailed balance versions of the gravitational sector, and we include the matter and radiation sectors. Allowing for variation of all the parameters of the theory, we construct the likelihood contours and we conclude that in $1\sigma$ confidence $\lambda$ is restricted to $|\lambda-1|\lesssim0.02$, while its best fit value is $|\lambda_{b.f}-1|\approx0.002$. Although this observational analysis restricts the running parameter $\lambda$ very close to its IR value 1, it does not enlighten the discussion about the theory's possible conceptual and theoretical problems.

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of co-moving side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, which evolves from a sub-Mpc scale at the beginning of reionization to R > 10 Mpc towards the end. We present a detailed statistical description of our results, and compare them with a numerical scheme based on the analytic model by Furlanetto, Zaldarriaga, and Hernquist. We find that the analytic calculation reproduces the size distribution of HII regions and the 21 cm power spectrum of the radiative transfer simulation remarkably well. The ionization field from the simulation, however, has more small scale structure than the analytic calculation, owing to Poisson scatter in the simulated abundance of galaxies on small scales. We propose and validate a simple scheme to incorporate this scatter into our calculations. Our results suggest that analytic calculations are sufficiently accurate to aid in predicting and interpreting the results of future 21 cm surveys. In particular, our fast numerical scheme is useful for forecasting constraints from future 21 cm surveys, and in constructing mock surveys to test data analysis procedures.

The third US Naval Observatory (USNO) CCD Astrograph Catalog, UCAC3 was released at the IAU General Assembly on 2009 August 10. It is the first all-sky release in this series and contains just over 100 million objects, about 95 million of them with proper motions, covering about R = 8 to 16 magnitudes. Current epoch positions are obtained from the observations with the 20 cm aperture USNO Astrograph's "red lens", equipped with a 4k by 4k CCD. Proper motions are derived by combining these observations with over 140 ground- and space-based catalogs, including Hipparcos/Tycho and the AC2000.2, as well as unpublished measures of over 5000 plates from other astrographs. For most of the faint stars in the Southern Hemisphere the Yale/San Juan first epoch plates from the SPM program (YSJ1) form the basis for proper motions. These data are supplemented by all-sky Schmidt plate survey astrometry and photometry obtained from the SuperCOSMOS project, as well as 2MASS near-IR photometry. Major differences of UCAC3 data as compared to UCAC2 include a completely new raw data reduction with improved control over systematic errors in positions, significantly improved photometry, slightly deeper limiting magnitude, coverage of the north pole region, greater completeness by inclusion of double stars and weak detections. This of course leads to a catalog which is not as "clean" as UCAC2 and problem areas are outlined for the user in this paper. The positional accuracy of stars in UCAC3 is about 15 to 100 mas per coordinate, depending on magnitude, while the errors in proper motions range from 1 to 10 mas/yr depending on magnitude and observing history, with a significant improvement over UCAC2 achieved due to the re-reduced SPM data and inclusion of more astrograph plate data unavailable at the time of UCAC2.

We reexamine the evolution of thermal relic particle abundances for the case where the interaction rate depends on the particle velocities. For the case of Sommerfeld enhancement, we show that the standard analytic approximation, modified in a straightforward way, provides an estimate of the relic particle abundance that is accurate to within 10% (in comparison to less than 1% error for the non-Sommerfeld-enhanced case). We examine the effect of kinetic decoupling on relic particle abundances when the interaction rate depends on the velocity. For the case of pure p-wave annihilation, the effect of kinetic decoupling is an increase in the relic abundance, but the effect is negligible when the kinetic decoupling temperature is much less than the chemical decoupling temperature. For the case of Sommerfeld-enhanced s-wave annihilations, after kinetic decoupling occurs, annihilations continue to change the particle abundance down to arbitrarily low temperatures, until either matter domination begins or the Sommerfeld effect cuts off. We derive analytic approximations to give the final relic particle abundances for both of these cases.

Since the violent relaxation in hierarchical merging is incomplete, elliptical galaxies retain a wealth of information about their formation pathways in their present-day orbital structure. A variety of evidence indicates that gas-rich major mergers play an important role in the formation of elliptical galaxies. We simulate 1:1 disk mergers at seven different initial gas fractions ranging from 0 to 40%, using the TreeSPH code Gadget-2. We classify the stellar orbits in each remnant and construct radial profiles of the orbital content, intrinsic shape, and orientation. The dissipationless remnants are typically prolate-triaxial, dominated by box orbits within r_c ~ 1.5Reff, and by tube orbits in their outer parts. As the gas fraction increases, the box orbits within r_c are increasingly replaced by a population of short axis tubes (z-tubes) with near zero net rotation, and the remnants become progressively more oblate and round. The long axis tube (x-tube) orbits are highly streaming and relatively insensitive to the gas fraction, implying that their angular momentum is retained from the dynamically cold initial conditions. Outside r_c, the orbital structure is essentially unchanged by the gas. The 15-20% gas remnants often display disk-like kinematically distinct cores (KDCs). These remnants show an interesting resemblance, in both their velocity maps and intrinsic orbital structure, to the KDC galaxy NGC4365 (van den Bosch et al. 2008). At 30-40% gas, the remnants are rapidly rotating, with sharp embedded disks on ~ 1Reff scales. We predict a characteristic, physically intuitive orbital structure for 1:1 disk merger remnants, with a distinct transition between 1 and 3Reff that will be readily observable with combined data from the 2D kinematics surveys SAURON and SMEAGOL.

We use observational data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis (BBN), to constrain the cosmological scenarios governed by Horava-Lifshitz gravity. We consider both the detailed and non-detailed balance versions of the gravitational sector, and we include the matter and radiation sectors. We conclude that the detailed-balance scenario cannot be ruled out from the observational point of view, however the corresponding likelihood contours impose tight constraints on the involved parameters. The scenario beyond detailed balance is compatible with observational data, and we present the corresponding stringent constraints and contour-plots of the parameters. Although this analysis indicates that Horava-Lifshitz cosmology can be compatible with observations, it does not enlighten the discussion about its possible conceptual and theoretical problems.

The amplitude of fluctuations in the Ly-a forest on small spatial scales is sensitive to the temperature of the IGM and its spatial fluctuations. The temperature of the IGM and its spatial variations contain important information about hydrogen and helium reionization. We present a new measurement of the small-scale structure in the Ly-a forest from 40 high resolution, high signal-to-noise, VLT spectra at z=2.2-4.2. We convolve each Ly-a forest spectrum with a suitably chosen wavelet filter, which allows us to extract the amount of small-scale structure in the forest as a function of position across each spectrum. We compare these measurements with high resolution hydrodynamic simulations of the Ly-a forest which track more than 2 billion particles. This comparison suggests that the IGM temperature close to the cosmic mean density (T_0) peaks near z=3.4, at which point it is greater than 20,000 K at 2-sigma confidence. The temperature at lower redshift is consistent with the fall-off expected from adiabatic cooling ($T_0 \propto (1+z)^2$), after the peak temperature is reached near z=3.4. At z=4.2 our results favor a temperature of T_0 = 15-20,000 K. However, owing mostly to uncertainties in the mean transmitted flux at this redshift, a cooler IGM model with T_0 = 10,000 K is only disfavored at the 2-sigma level here, although such cool IGM models are strongly discrepant with the z ~ 3-3.4 measurement. We do not detect large spatial fluctuations in the IGM temperature at any redshift covered by our data set. The simplest interpretation of our measurements is that HeII reionization completes sometime near z ~ 3.4, although statistical uncertainties are still large [Abridged].

We examine a dark energy model where a scalar unparticle degree of freedom plays the role of quintessence. In particular, we study a model where the unparticle degree of freedom has a standard kinetic term and a simple mass potential, the evolution is slowly rolling and the field value is of the order of the unparticle energy scale ($\lambda_u$). We study how the evolution of $w$ depends on the parameters $B$ (a function of the unparticle scaling dimension $d_u$), the initial value of the field $\phi_i$ (or equivalently, $\lambda_u$) and the present matter density $\Omega_{m0}$. We use observational data from Type Ia supernovae, BAO and CMB to constrain the model parameters and find that these models are not ruled out by the observational data. From a theoretical point of view, an unparticle dark energy model is very attractive, since unparticles (being bound states of fundamental fermions) are protected from radiative corrections. Further, coupling of unparticles to the standard model fields can be arbitrarily suppressed by raising the fundamental energy scale $M_F$, making the unparticle dark energy model free of most of the problems that plague conventional scalar field quintessence models.

We derive slow-roll conditions for thawing k-essence with a separable Lagrangian $p(X,\phi)=F(X)V(\phi)$. We examine the evolution of the equation of state parameter, $w$, as a function of the scale factor $a$, for the case where $w$ is close to -1. We find two distinct cases, corresponding to $X \approx 0$ and $F_X \approx 0$, respectively. For the case where $X\approx0$ the evolution of $\phi$ and hence $w$ is described by only two parameters, and $w(a)$ is model-independent and coincides with similar behavior seen in thawing quintessence models. This result also extends to non-separable Lagrangians where $X\approx0$. For the case $F_X \approx 0$, an expression is derived for $w(a)$, but this expression depends on the potential $V(\phi)$, so there is no model-independent limiting behavior. For the $X \approx 0$ case, we derive observational constraints on the two parameters of the model, $w_0$ (the present-day value of $w$), and the $K$, which parametrizes the curvature of the potential. We find that the observations sharply constrain $w_0$ to be close to -1, but provide very poor constraints on $K$.

The properties of elliptical galaxies are broadly consistent with simulated remnants of gas-rich mergers between spirals, motivating more detailed studies of the imprint of this formation mechanism on the remnant distribution function. Gas has a strong impact on the non-Gaussian shapes of the line-of-sight velocity distributions (LOSVDs) of the merger remnant, owing to the embedded disk that forms out of the gas that retains its angular momentum during the merger, and the strong central mass concentration from the gas that falls to the center. The deviations from Gaussianity are parametrized by the Gauss-Hermite moments h_3 and h_4, which are related to the skewness and kurtosis of the LOSVDs. We quantify the dependence of the (h_3,h_4)-v/sigma relations on the initial gas fraction of the progenitor disks in 1:1 mergers, using Gadget-2 simulations including star formation, radiative cooling, and feedback from supernovae and AGN. For gas fractions f_gas < ~15% the overall correlation between h_3 and v/sigma is weak, consisting of a flat negatively correlated component arising from edge-on viewing angles plus a steep positively correlated part from face-on projections. The spread in v/sigma values decreases toward high positive h_4, and there is a trend toward lower h_4 as the gas fraction increases from 0 to 15%. For f_gas > ~20% the (h_3,4)- v/sigma distributions look quite different - there is a tight negative h_3- v/sigma correlation, and a wide spread in v/sigma values at all h_4, in better agreement with observations. Re-mergers of the high-f_gas remnants (dry mergers) produce slowly rotating systems with nearly Gaussian LOSVDs. We explain all of these trends in terms of the underlying orbit structure of the remnants, as molded by their dissipative formation histories.

We demonstrate the scale dependence of the growth function of cosmological perturbations in dark energy models based on General Relativity. This scale dependence is more prominent on cosmological scales of $100h^{-1}Mpc$ or larger. We derive a new scale dependent parametrization which generalizes the well known Newtonian approximation result $f_0(a)\equiv \frac{d\ln \delta_0}{d\ln a}=\om(a)^\gamma$ ($\gamma ={6/11}$ for \lcdm) which is a good approximation on scales less than $50h^{-1}Mpc$. Our generalized parametrization is of the form $f(a)=\frac{f_0(a)}{1+\xi(a,k)}$ where $\xi(a,k)=\frac{3 H_{0}^{2} \omm}{a k^2}$. We demonstrate that this parametrization fits the exact result of a full general relativistic evaluation of the growth function up to horizon scales for both \lcdm and dynamical dark energy. In contrast, the scale independent parametrization does not provide a good fit on scales beyond 5% of the horizon scale ($k\simeq 0.01 h^{-1}Mpc$).

We examine dark energy models in which a phantom field $\phi$ is rolling near a local minimum of its potential $V(\phi)$.We require that $(1/V)(dV/d\phi) \ll 1$, but $(1/V)(d^2 V/d\phi^2)$ can be large. Using techniques developed in the context of hilltop quintessence, we derive a general expression for $w$ as a function of the scale factor, and as in the hilltop case, we find that the dynamics of the field depend on the value of $(1/V)(d^2 V/d\phi^2)$ near the mimimum. Our general result gives a value for $w$ that is within 1% of the true (numerically-derived) value for all of the particular cases examined. Our expression for $w(a)$ reduces to the previously-derived phantom slow-roll result of Sen and Scherrer in the limit where the potential is flat, $(1/V)(dV/d\phi) \ll 1$.

We propose a simple dark energy model with the following properties: the model predicts a late-time dark radiation component that is not ruled out by current observational data, but which produces a distinctive time-dependent equation of state w(z) for z < 3. The dark energy field can be coupled strongly enough to Standard Model particles to be detected in colliders, and the model requires only modest additional particle content and little or no fine-tuning other than a new energy scale of order milli-electron volts.

We examine dark energy models in which a quintessence or a phantom field, $\phi$, rolls near the vicinity of a local minimum or maximum, respectively, of its potential $V(\phi)$. Under the approximation that $(1/V)(dV/d\phi) \ll 1$, [although $(1/V)(d^2 V/d\phi^2)$ can be large], we derive a general expression for the equation of state parameter $w$ as a function of the scale factor for these models. The dynamics of the field depends on the value of $(1/V)(d^2 V/d\phi^2)$ near the extremum, which describes the potential curvature. For quintessence models, when $(1/V)(d^2 V/d\phi^2)<3/4$ at the potential minimum, the equation of state parameter $w(a)$ evolves monotonically, while for $(1/V)(d^2 V/d\phi^2)>3/4$, $w(a)$ has oscillatory behavior. For phantom fields, the dividing line between these two types of behavior is at $(1/V)(d^2 V/d\phi^2) = -3/4$. Our analytical expressions agree within 1% with the exact (numerically-derived) behavior, for all of the particular cases examined, for both quintessence and phantom fields. We present observational constraints on these models.

We examine the accuracy of the growth equation $\ddot{\delta} + 2H\dot{\delta} - 4\pi G\rho\delta = 0$, which is ubiquitous in the cosmological literature, in the context of the Newtonian gauge. By comparing the growth predicted by this equation to a numerical solution of the linearized Einstein equations in the $\Lambda$CDM scenario, we show that while this equation is a reliable approximation on small scales ($k\gtrsim $h Mpc$^{-1}$), it can be disastrously inaccurate ($\sim 10^4% $) on larger scales in this gauge. We propose a modified version of the growth equation for the Newtonian gauge, which while preserving the simplicity of the original equation, provides considerably more accurate results. We examine the implications of the failure of the growth equation on a few recent studies, aimed at discriminating general relativity from modified gravity, which use this equation as a starting point. We show that while the results of these studies are valid on small scales, they are not reliable on large scales or high redshifts, if one works in the Newtonian gauge. Finally, we discuss the growth equation in the synchronous gauge and show that the corrections to the Poisson equation are exactly equivalent to the difference between the overdensities in the synchronous and Newtonian gauges.

Observations of the intergalactic medium (IGM) suggest that quasars reionize HeII in the IGM at z ~ 3. We have run a set of 190 and 430 comoving Mpc simulations of HeII being reionized by quasars to develop an understanding of the nature of HeII reionization and its potential impact on observables. We find that HeII reionization heats regions in the IGM by as much as 25,000 K above the temperature that is expected otherwise, with the volume-averaged temperature increasing by ~ 12,000 K and with large temperature fluctuations on ~ 50 Mpc scales. Much of the heating occurs far from QSOs by hard photons. We find a temperature-density equation of state of gamma -1 ~ 0.3 during HeII reionization, but with a wide dispersion in this relation having sigma ~ 10^4 K. HeII reionization by the observed population of quasars cannot produce an inverted relation (gamma - 1 < 0). Our simulations are consistent with the observed evolution in the mean transmission of the HeII Ly-alpha forest. We argue that the heat input due to HeII reionization is unable to cause the observed depression at z = 3.2 in the HI Ly-alpha forest opacity as has been suggested. We investigate how uncertainties in the properties of QSOs and of HeII Lyman-limit systems influence our predictions.

We examine hilltop quintessence models, in which the scalar field is rolling near a local maximum in the potential, and w is close to -1. We first derive a general equation for the evolution of the scalar field in the limit where w is close to -1. We solve this equation for the case of hilltop quintessence to derive w as a function of the scale factor; these solutions depend on the curvature of the potential near its maximum. Our general result is in excellent agreement (delta w < 0.5%) with all of the particular cases examined. It works particularly well (delta w < 0.1%) for the pseudo-Nambu-Goldstone Boson potential. Our expression for w(a) reduces to the previously-derived slow-roll result of Sen and Scherrer in the limit where the curvature goes to zero. Except for this limiting case, w(a) is poorly fit by linear evolution in a.

The effect of quintessence perturbations on the ISW effect is studied for a mixed dynamical scalar field dark energy (DDE) and pressureless perfect fluid dark matter. A new and general methodology is developed to track the growth of the perturbations, which uses only the equation of state (EoS) parameter $w_{\rm DDE} (z) \equiv p_{\rm DDE}/\rho_{\rm DDE}$ of the scalar field DDE, and the initial values of the the relative entropy perturbation (between the matter and DDE) and the intrinsic entropy perturbation of the scalar field DDE as inputs. We also derive a relation between the rest frame sound speed $\hat{c}_{s,{\rm DDE}}^2$ of an arbitrary DDE component and its EoS $w_{\rm DDE} (z)$. We show that the ISW signal differs from that expected in a $\Lambda$CDM cosmology by as much as +20% to -80% for parameterizations of $w_{\rm DDE}$ consistent with SNIa data, and about $\pm$ 20% for parameterizations of $w_{\rm DDE}$ consistent with SNIa+CMB+BAO data, at 95% confidence. Our results indicate that, at least in principle, the ISW effect can be used to phenomenologically distinguish a cosmological constant from DDE.

We study the origin and properties of 'extra' or 'excess' central light in the surface brightness profiles of cusp or power-law ellipticals. Dissipational mergers give rise to two-component profiles: an outer profile established by violent relaxation acting on stars present in the progenitors prior to the final merger, and an inner stellar population comprising the extra light, formed in a compact starburst. Combining a large set of hydrodynamical simulations with data that span a broad range of profiles and masses, we show that this picture is borne out -- cusp ellipticals are indeed 'extra light' ellipticals -- and examine how the properties of this component scale with global galaxy properties. We show how to robustly separate the 'extra' light, and demonstrate that observed cusps are reliable tracers of the degree of dissipation in the spheroid-forming merger. We show that the typical degree of dissipation is a strong function of stellar mass, tracing observed disk gas fractions at each mass. We demonstrate a correlation between extra light content and effective radius at fixed mass: systems with more dissipation are more compact. The outer shape of the light profile does not depend on mass, with a mean outer Sersic index ~2.5. We explore how this relates to shapes, kinematics, and stellar population gradients. Simulations with the gas content needed to match observed profiles also reproduce observed age, metallicity, and color gradients, and we show how these can be used as tracers of the degree of dissipation in spheroid formation.

Oscillating scalar fields, with an oscillation frequency much greater than the expansion rate, have been proposed as models for dark energy. We examine these models, with particular emphasis on the evolution of the ratio of the oscillation frequency to the expansion rate. We show that this ratio always increases with time if the dark energy density declines less rapidly than the background energy density. This allows us to classify oscillating dark energy models in terms of the epoch at which the oscillation frequency exceeds the expansion rate, which is effectively the time at which rapid oscillations begin. There are three basic types of behavior: early oscillation models, in which oscillations begin during the matter-dominated era, late oscillation models, in which oscillations begin after scalar-field domination, and non-oscillating models. We examine a representative set of models (those with power-law potentials) and determine the parameter range giving acceptable agreement with the supernova observations. We show that a subset of all three classes of models can be consistent with the observational data.

It is possible that the properties of HII regions during reionization depend sensitively on many poorly constrained quantities (the nature of the ionizing sources, the clumpiness of the gas in the IGM, the degree to which photo-ionizing feedback suppresses the abundance of low mass galaxies, etc.), making it extremely difficult to interpret upcoming observations of this epoch. We demonstrate that the actual situation is more encouraging, using a suite of radiative transfer simulations, post-processed on outputs from a 1024^3, 94 Mpc N-body simulation. Analytic prescriptions are used to incorporate small-scale structures that affect reionization, yet remain unresolved in the N-body simulation. We show that the morphology of the HII regions is most dependent on the global ionization fraction x_i. This is not to say that the bubble morphology is completely independent of all parameters besides x_i. The next most important dependence is that of the nature of the ionizing sources. The rarer the sources, the larger and more spherical the HII regions become. The typical bubble size can vary by as much as a factor of 4 at fixed x_i between different possible source prescriptions. The final relevant factor is the abundance of minihalos or of Lyman-limit systems. These systems suppress the largest bubbles from growing, and the magnitude of this suppression depends on the thermal history of the gas as well as the rate at which minihalos are photo-evaporated. We find that neither source suppression owing to photo-heating nor gas clumping significantly affect the large-scale structure of the HII regions. We discuss how observations of the 21cm line with MWA and LOFAR can constrain properties of reionization, and we study the effect patchy reionization has on the statistics of Lyman-alpha emitting galaxies. [abridged]

We study the origin and properties of 'extra' or 'excess' central light in the surface brightness profiles of gas-rich merger remnants. Combining a large set of hydrodynamical simulations with data on observed mergers (spanning a broad range of profiles at various masses and degrees of relaxation), we show how to robustly separate the physically meaningful extra light -- stellar populations formed in a compact central starburst during a gas-rich merger -- from the outer profile established by violent relaxation acting on stars already present in the progenitors prior to the final merger. This separation is sensitive to the profile treatment, and we demonstrate that certain fitting procedures can yield physically misleading results. We show that our method reliably recovers the younger starburst population, and examine how the properties of this component scale with mass, gas content, and other aspects of the progenitors. We consider the time evolution of profiles in different bands, and estimate biases introduced by observational studies at different times and wavelengths. We show that extra light is ubiquitous in observed and simulated gas-rich merger remnants, with sufficient mass (~3-30% of the stellar mass) to explain the discrepancy in the maximum phase-space densities of ellipticals and their progenitor spirals. The nature of this central component provides powerful new constraints on the formation histories of observed systems.

We show that near-infrared observations of the red side of the Ly-alpha line from a single gamma ray burst (GRB) afterglow cannot be used to constrain the global neutral fraction of the intergalactic medium (IGM), x_H, at the GRB's redshift to better than ~0.3. Some GRB sight-lines will encounter more neutral hydrogen than others at fixed x_H owing to the patchiness of reionisation. GRBs during the epoch of reionization will often bear no discernible signature of a neutral IGM in their afterglow spectra. We discuss the constraints on x_H from the z = 6.3 burst, GRB050904, and quantify the probability of detecting a neutral IGM using future spectroscopic observations of high-redshift, near-infrared GRB afterglows. Assuming an observation with signal-to-noise similar to the Subaru FOCAS spectrum of GRB050904 and that the column density distribution of damped Ly-alpha absorbers is the same as measured at lower redshifts, a GRB from an epoch when x_H = 0.5 can be used to detect a partly neutral IGM at 98% confidence level 10% of the time (and, for an observation with three times the sensitivity, 30% of the time).

We investigate the clustering properties of a dynamical dark energy component. In a cosmic mix of a pressureless fluid and a light scalar field, we follow the linear evolution of spherical matter perturbations. We find that the scalar field tends to form underdensities in response to the gravitationally collapsing matter. We thoroughly investigate these voids for a variety of initial conditions, explain the physics behind their formation and consider possible observational implications. Detection of dark energy voids will clearly rule out the cosmological constant as the main source of the present acceleration.

We show that observations of high-redshift Ly-alpha emitters (LAEs) have the potential to provide definitive evidence for reionization in the near future. Using 200 Mpc radiative transfer simulations, we calculate the effect that patchy reionization has on the line profile, on the luminosity function, and, most interestingly, on the clustering of emitters for several realistic models of reionization. Reionization increases the measured clustering of emitters, and we show that this enhancement would be essentially impossible to attribute to anything other than reionization. Our results motivate looking for the signature of reionization in existing LAE data. We find that for stellar reionization scenarios the angular correlation function of the 58 LAEs in the Subaru Deep Field z = 6.6 photometric sample is more consistent with a fully ionized universe (mean volume ionized fraction x_i = 1) than a universe with x_i < 0.5 at >2-sigma confidence level. Measurements in the next year on Subaru will increase their z = 6.6 LAE sample by a factor of five and tighten these limits. If the clustering signature of reionization is detected in a LAE survey, a comparison with a Lyman-break or a H-alpha survey in the same field would confirm the reionization hypothesis. We discuss the optimal LAE survey specifications for detecting reionization, with reference to upcoming programs.

Lyman-alpha forest absorption spectra towards quasars at z ~ 6 show regions of enhanced transmission close to their source. Several authors have argued that the apparently small sizes of these regions indicate that quasar ionization fronts at z >~ 6 expand into a largely or partly neutral intergalactic medium (IGM). Assuming that the typical region in the IGM is reionized by z <= 6, as is suggested by Ly-a forest observations, we argue that at {\em least} 50% of the volume of the IGM was reionized before the highest redshift quasars turned on. Further, even if the IGM is as much as 50% neutral at quasar turn-on, the quasars are likely born into large galaxy-generated HII regions. The HII regions during reionization are themselves clustered, and using radiative transfer simulations, we find that long skewers through the IGM towards quasar progenitor halos pass entirely through ionized bubbles, even when the IGM is half neutral. These effects have been neglected in most previous analyses of quasar proximity zones, which assumed a spatially {\em uniform} neutral fraction. We model the subsequent ionization from a quasar, and construct mock Ly-a forest spectra. Our mock absorption spectra are more sensitive to the level of small-scale structure in the IGM than to the volume-averaged neutral fraction, and suggest that existing proximity-zone size measurements are compatible with a fully ionized IGM. However, we mention several improvements in our modeling that are necessary to make more definitive conclusions.

We consider the contribution of 3rd and 4th order terms to the power spectrum of 21 cm brightness temperature fluctuations during the epoch of reionization, which arise because the 21 cm brightness temperature involves a product of the hydrogenic neutral fraction and the gas density. The 3rd order terms vanish for Gaussian random fields, and have been previously neglected or ignored. We measure these higher order terms from radiative transfer simulations and estimate them using cosmological perturbation theory. In our simulated models, the higher order terms are significant: neglecting them leads to a >~100% error in 21 cm power spectrum predictions on scales of k >~ 1 Mpc^{-1} when the neutral fraction is <x_H>~0.5. The higher order terms have a simple physical interpretation. On small scales they are produced by gravitational mode coupling. Small scale structure grows more readily in large-scale overdense regions, but the same regions tend to be ionized and hence do not contribute to the 21 cm signal. This acts to suppress the influence of non-linear density fluctuations and the small-scale amplitude of the 21 cm power spectrum. These results modify earlier intuition that the 21 cm power spectrum simply traces the density power spectrum on scales smaller than that of a typical bubble, and imply that small scale measurements contain more information about the nature of the ionizing sources than previously believed. On large scales, higher order moments are not directly related to gravity. They are non-zero because over-dense regions tend to ionize first and are important in magnitude at late times owing to the large fluctuations in the neutral fraction. (Abridged)

We use numerical simulations to study the kinematic structure of remnants formed from mergers of equal-mass disk galaxies. In particular, we show that remnants of dissipational mergers, which include the radiative cooling of gas, star formation, feedback from supernovae, and the growth of supermassive black holes, are smaller, rounder, have, on average, a larger central velocity dispersion, and show significant rotation compared to remnants of dissipationless mergers. The increased rotation speed of dissipational remnants owes its origin to star formation that occurs in the central regions during the galaxy merger. We have further quantified the anisotropy, three-dimensional shape, minor axis rotation, and isophotal shape of each merger remnant, finding that dissipational remnants are more isotropic, closer to oblate, have the majority of their rotation along their major axis, and are more disky than dissipationless remnants. Individual remnants display a wide variety of kinematic properties. A large fraction of the dissipational remnants are oblate isotropic rotators. Many dissipational, and all of the dissipationless, are slowly rotating and anisotropic. The remnants of gas-rich major mergers can well-reproduce the observed distribution of projected ellipticities, rotation parameter (V/\sigma)*, kinematic misalignments, Psi, and isophotal shapes. The dissipationless remnants are a poor match to this data. Our results support the merger hypothesis for the origin of low-luminosity elliptical galaxies provided that the progenitor disks are sufficiently gas-rich, however our remnants are a poor match to the bright ellipticals that are slowly rotating and uniformly boxy.

Computing the perturbation spectrum in the recently proposed Island Cosmology remains an open problem. In this paper we present a classical computation of the perturbations generated in this scenario by assuming that the NEC-violating field behaves as a classical phantom field. Using an exactly-solvable potential, we show that the model generates a scale-invariant spectrum of scalar perturbations, as well as a scale-invariant spectrum of gravitational waves. The scalar perturbations can have sufficient amplitude to seed cosmological structure, while the gravitational waves have a vastly diminished amplitude.

We propose an alternate cosmological model in which our observable universe is an island in a cosmological constant sea. Initially the universe is filled with cosmological constant of the currently observed value but is otherwise empty. In this eternal or semi-eternal de Sitter spacetime, we show that local quantum fluctuations (upheavals) can violate the null energy condition and create islands of matter. The perturbation spectra of quantum fields other than that responsible for the upheaval, are shown to be scale invariant. With further cosmic evolution the island disappears and the local universe returns to its initial cosmological constant dominated state.

Neutrino processes in semi-transparent supernova matter, opaque to semi-transparent protoneutron star matter, and catalyzed neutron stars are discussed.

Explicit expressions for the differential and total rates and emissivities of neutrino pairs from the photo-neutrino process $e^\pm + \gamma \to e^\pm + \nu + \bar\nu$ in hot and dense matter are derived. Full information about the emitted neutrinos is retained by evaluating the squared matrix elements for this process which was hitherto bypassed through the use of Lenard's identity in obtaining the total neutrino emissivities. Accurate numerical results are presented for widely varying conditions of temperature and density. Analytical results helpful in understanding the qualitative behaviors of the rates and emissivities in limiting situations are derived. The corresponding production and absorption kernels in the source term of the Boltzmann equation for neutrino transport are developed. The appropriate Legendre coefficients of these kernels, in forms suitable for multigroup flux-limited diffusion schemes are also provided.

The differential rates and emissivities of neutrino pairs from an equilibrium plasma are calculated for the wide range of density and temperature encountered in astrophysical systems. New analytical expressions are derived for the differential emissivities which yield total emissivities in full agreement with those previously calculated. The photon and plasmon pair production and absorption kernels in the source term of the Boltzmann equation for neutrino transport are provided. The appropriate Legendre coefficients of these kernels, in forms suitable for multi-group flux-limited diffusion schemes are also computed.

In scenarios with large extra dimensions and TeV scale gravity ultrahigh energy neutrinos produce black holes in their interactions with the nucleons. We show that ICECUBE and OWL may observe large number of black hole events and provide valuable information about the fundamental Planck scale and the number of extra dimensions. OWL is especially well suited to observe black hole events produced by neutrinos from the interactions of cosmic rays with the 3 K background radiation. Depending on the parameters of the scenario of large extra dimensions and on the flux model, as many as 28 events per year are expected for a Planck scale of 3 TeV.

We evaluate the upward shower and muon event rates for two characteristic four neutrino mixing models for extragalactic neutrinos, as well as for the atmospheric neutrinos, with energy thresholds of 1 TeV, 10 TeV and 100 TeV. We show that by comparing the shower to muon event rates, one can distinguish between oscillation and no-oscillation models. By measuring shower and muon event rates for energy thresholds of 10 TeV and 100 TeV, and by considering their ratio, it is possible to use extragalactic neutrino sources to determine the type of four-flavor mixing pattern. We find that one to ten years of data taking with kilometer-size detector has a very good chance of providing valuable information about the physics beyond the Standard Model.

The photonuclear contribution to charged lepton energy loss has been re-evaluated taking into account HERA results on real and virtual photon interactions with nucleons. With large $Q^2$ processes incorporated, the average muon range in rock for muon energies of $10^9$ GeV is reduced by only 5% as compared with the standard treatment. We have calculated the tau energy loss for energies up to $10^9$ GeV taking into consideration the decay of the tau. A Monte Carlo evaluation of tau survival probability and range show that at energies below $10^7-10^8$ GeV, depending on the material, only tau decays are important. At higher energies the tau energy losses are significant, reducing the survival probability of the tau. We show that the average range for tau is shorter than its decay length and reduce to 17 km in water for an incident tau energy of $10^9$ GeV, as compared with its decay length of 49 km at that energy. In iron, the average tau range is 4.7 km for the same incident energy.

The appearance of high energy tau neutrinos due to $\nu_\mu \to \nu_\tau$ oscillations of extragalactic neutrinos can be observed by measuring the neutrino induced upward hadronic and electromagnetic showers and upward muons. We evaluate quantitatively the tau neutrino regeneration in the Earth for a variety of extragalactic neutrino fluxes. Charged-current interactions of the upward tau neutrinos below and in the detector, and the subsequent tau decay create muons or hadronic and electromagnetic showers. The background for these events are muon neutrino and electron neutrino charged-current and neutral-current interactions, where in addition to extragalactic neutrinos, we consider atmospheric neutrinos. We find significant signal to background ratios for the hadronic/electromagnetic showers with energies above 10 TeV to 100 TeV initiated by the extragalactic neutrinos. We show that the tau neutrinos from point sources also have the potential for discovery above a 1 TeV threshold. A kilometer-size neutrino telescope has a very good chance of detecting the appearance of tau neutrinos when both muon and hadronic/electromagnetic showers are detected.

The frequency shift of spectral lines is most often explained by the Doppler Effect in terms of relative motion, whereas the Doppler broadening of a particular line mainly depends on the absolute temperature. The Wolf effect on the other hand deals with the correlation induced spectral change and explains both the broadening and shift of the spectral lines. In this framework a relation between the width of the spectral line is related to the redshift z for the line and hence with the distance. For smaller values of z a relation similar to the Tully-Fisher relation can be obtained and for larger values of z a more general relation can be constructed. The derivation of this kind of relation based on dynamic multiple scattering theory may play a significant role in explaining the overall spectra of quasi stellar objects. We emphasize that this mechanism is not applicable for nearby galaxies, $z \leq 1$.

We study the persistence of warps in galactic discs in the presence of massive halos. A disc is approximated by a set of massive rings, while a halo is represented by a conventional n-body simulation. We confirm the conclusion of Nelson & Tremaine (1995) that a halo responds strongly to an embedded precessing disc. This response invalidates the approximations made by in the derivation of classical `modified tilt' modes. We show that the response of the halo causes the line of nodes of a disc that starts from a modified tilt mode to wind up within a few dynamical times. We explain this finding in terms of the probable spectrum of true normal modes of a combined disc-halo system.

The shape of galaxies depends on their orbital populations. These populations change through capture into and escape from resonance. Capture problems fall into distinct cases depending upon the shape of the potential well. To visualise the effective potential well for orbital capture, a diagrammatic approach to the resonant perturbation theory of Born is presented. These diagrams we call equiaction sections. To illustrate their use, we present examples drawn from both galactic and Solar System dynamics. The probability of capture for generic shapes of the potential well is calculated. A number of predictions are made. First, there are barred galaxies that possess two outer rings of gas and stars (type R1R2). We show how to relate changes in the pattern speed and amplitude of the bar to the strength of the two rings. Secondly, under certain conditions, small disturbances can lead to dramatic changes in orbital shape. This can be exploited as a mechanism to pump counter-rotating stars and gas into the nuclei of disk galaxies. Tidal resonant forcing of highly inclined orbits around a central mass causes a substantial increase in the likelihood of collision. Thirdly, the angular momentum of a potential well is changed by the passage of stars across or capture into the well. This can lead to the creation of holes, notches and high velocity tails in the stellar distribution function, whose form we explicitly calculate.

We investigate the formation of clusters of galaxies in an expanding universe using a new code that regrids at a region of high density. In particular we investigate two models for the initial conditions, both with the standard CDM power spectrum - one has $\Omega = 1$ and the other $\Omega = 0.2$. Both universes have $H_0 = 100 \kms\mpc^{-1}$ and $\Lambda = 0$. The level of substructure in the final cluster can be used as a discriminant of the cosmic density. We discuss various statistics which can be measured observationally from clusters of galaxies, that can be used to discriminate between the two models. We find that most of the statistics that use the clusters' member galaxies may not be the best measures of substructure. Statistics that rely more on X-ray maps and other observables depending more directly on the mass distribution could be better discriminants of $\Omega$.

We investigate two possible effects of the tidal field induced by a spherical cluster on its elliptical galaxy members: the modification of the ellipticity of a spherical galaxy and the isophotal alignment in the cluster radial direction of a misaligned prolate galaxy. Numerical N-body simulations have been performed for radial and circular galactic orbits. The properties of the stars' zero--velocity surfaces in the perturbed galaxies are explored briefly, and the adiabaticity of the galaxy to the external field is discussed. For a choice of parameters characteristic of rich clusters we find that the induced ellipticity on a spherical galaxy is below or close to the detectability level. But we find that the tidal torque can result in significant isophotal alignment of the galaxies' major axis with the cluster radial direction if the galaxy is outside the cluster core radius. The time required for the alignment is very short compared with the Hubble time. A significant increase in the ellipticity of the outer isophotes of the prolate model is also found, but with no observable isophotal twisting. Our main prediction is an alignment segregation of the elliptical galaxy population according to whether their orbits lie mostly outside or inside the cluster core radius. These results also suggest that galactic alignment in rich clusters is not incompatible with a bottom-up galaxy formation scenario.