We consider perturbations of closed Friedmann universes. Perturbation modes of two lowest wavenumbers ($L=0$ and $1$) are generally known to be fictitious, but here we show that both are physical. The issue is more subtle in Einstein static universes where closed background space has a time-like Killing vector with the consequent occurrence of linearization instability. Proper solutions of the linearized equation need to satisfy the Taub constraint on a quadratic combination of first-order variables. We evaluate the Taub constraint in the two available fundamental gauge conditions, and show that in both gauges the $L\geq 1$ modes should accompany the $L=0$ (homogeneous) mode for vanishing sound speed, $c_{s}$. For $c_{s}^{2}>1/5$ (a scalar field supported Einstein static model belongs to this case with $c_s^2 = 1$), the $L\geq 2$ modes are known to be stable. In order to have a stable Einstein static evolutionary stage in the early universe, before inflation and without singularity, although the Taub constraint does not forbid it, we need to find a mechanism to suppress the unstable $L=0$ and $L=1$ modes.

We propose a new method to detect gravitational waves, based on spatial coherence interferometry with stellar light, as opposed to the conventional temporal coherence interferometry with laser sources. The proposed method detects gravitational waves by combining two coherent beams of light from a single distant star measured at two separate detectors, each on-board a satellite. For this experiment, no complex detector system is required, hence any associated acceleration noise can be minimized, allowing for the search of gravitational waves in the lower frequency range of $10^{-7} - 10^{-4}$ Hz, not covered by any experiment at the time of this paper. In this work, we present the detection sensitivity of the proposed stellar interferometer for gravitational waves by taking the detection response and potential noises into account. Furthermore, the proposed experimental setup is capable of searching for primordial black holes and studying the size of the target neutron star, which are also discussed in the paper.

Using the fully nonlinear and exact perturbation formulation with magnetohydrodynamics (MHD) in Minkowski background we derive first-order post-Newtonian (1PN) equations without imposing the slicing (temporal gauge) condition. The 1PN MHD formulation is complementary to our recently presented fully relativistic MHD combined with 0PN gravity available only in the maximal slicing. We present the 1PN MHD equations in two gauge conditions previously used in the literature and provide gauge transformation relations between different gauges. We derive the PN effects on MHD waves in a static homogeneous medium.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

In the circumnuclear region of the nearby radio galaxy NGC~1052. Our $\lambda$3mm VLBI observations with the Korean VLBI Network have spatially resolved the broad HCO$^{+}$ absorption at a velocity range of 1350--1850 km~s$^{-1}$ against a double-sided nuclear jet, and have revealed that the HCO$^{+}$ absorption is concentrated on the receding jet and the nuclear components. The distribution of the HCO$^{+}$ absorbing gas strongly supports the circumnuclear torus surrounding the supermassive black hole. From the estimations of the column density and the volume density of molecular hydrogen, the size of the molecular gas region in the torus is at least 1~pc. The broad spectral profile of HCO$^{+}$ is likely to be a blend of multiple gas clumps with various velocities. The HCO$^{+}$ absorption of NGC~1052 could trace complex kinematics in the vicinity of the active galactic nucleus, such as inflow, outflow, turbulence, and so on.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We present a fully nonlinear and exact perturbation formulation of Einstein's gravity with a general fluid and the ideal magnetohydrodynamics (MHD) without imposing the slicing (temporal gauge) condition. Using this formulation, we derive equations of special relativistic (SR) MHD in the presence of weak gravitation. The equations are consistently derived in the limit of weak gravity and action-at-a-distance in the maximal slicing. We show that in this approximation the relativistic nature of gravity does not affect the SR MHD dynamics, but SR effects manifest themselves in the metric, and thus gravitational lensing. Neglecting these SR effects may lead to an overestimation of lensing masses.

We present hydrodynamic equations with relativistic pressure and velocity in the presence of weak gravity, in a cosmological context. Previously we consistently derived special relativistic hydrodynamic equations with weak gravity in Minkowski background. With the relativistic pressure and velocity one cannot derive the cosmological counterpart by a simple transformation from equations in the Minkowski background. Here we present a proper derivation. We point out the potential importance of relativistic pressure and velocity in gravitational lensing.

We present perturbation initial conditions for two types of scalar field potential often used in the dark energy study: one inverse power-law and the other exponential. The solutions are presented in the presence of the $w =$ constant fluid ($w = 1/3$ for radiation fluid), a minimally coupled scalar field and a sub-dominating zero-pressure fluid (cold dark matter and baryon dust). We consider two gauge conditions, the $w$-fluid comoving gauge and the cold-dark-matter comoving gauge; solutions in the latter gauge are derived by the gauge transformation and this method can be applied to derive solutions in any other gauge condition.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

Late-type galaxies falling into a cluster would evolve being influenced by the interactions with both the cluster and the nearby cluster member galaxies. Most numerical studies, however, tend to focus on the effects of the former with little work done on those of the latter. We thus perform a numerical study on the evolution of a late-type galaxy interacting with neighboring early-type galaxies at high speed, using hydrodynamic simulations. Based on the information obtained from the Coma cluster, we set up the simulations for the case where a Milky Way-like late-type galaxy experiences six consecutive collisions with twice as massive early-type galaxies having hot gas in their halos at the closest approach distances of 15-65 kpc/h at the relative velocities of 1500-1600 km/s. Our simulations show that the evolution of the late-type galaxy can be significantly affected by the accumulated effects of the high-speed multiple collisions with the early-type galaxies, such as on cold gas content and star formation activity of the late-type galaxy, particularly through the hydrodynamic interactions between cold disk and hot gas halos. We find that the late-type galaxy can lose most of its cold gas after the six collisions and have more star formation activity during the collisions. By comparing our simulation results with those of galaxy-cluster interactions, we claim that the role of the galaxy-galaxy interactions on the evolution of late-type galaxies in clusters could be comparable with that of the galaxy-cluster interactions, depending on the dynamical history.

Studies of high-multiplicity, tightly-packed planetary systems suggest that dynamical instabilities are common and affect both the orbits and planet structures, where the compact orbits and typically low densities make physical collisions likely outcomes. Since the structure of many of these planets is such that the mass is dominated by a rocky core, but the volume is dominated by a tenuous gas envelope, the sticky-sphere approximation, used in dynamical integrators, may be a poor model for these collisions. We perform five sets of collision calculations, including detailed hydrodynamics, sampling mass ratios and core mass fractions typical in Kepler Multis. In our primary set of calculations, we use Kepler-36 as a nominal remnant system, as the two planets have a small dynamical separation and an extreme density ratio. We use an N-body code, Mercury 6.2, to integrate initially unstable systems and study the resultant collisions in detail. We use these collisions, focusing on grazing collisions, in combination with realistic planet models created using gas profiles from Modules for Experiments in Stellar Astrophysics and core profiles using equations of state from Seager et al. (2007), to perform hydrodynamic calculations, finding scatterings, mergers, and even a potential planet-planet binary. We dynamically integrate the remnant systems, examine the stability, and estimate the final densities, finding the remnant densities are sensitive to the core masses, and collisions result in generally more stable systems. We provide prescriptions for predicting the outcomes and modeling the changes in mass and orbits following collisions for general use in dynamical integrators.

We present a complete set of exact and fully non-linear equations describing all three types of cosmological perturbations -- scalar, vector and tensor perturbations. We derive the equations in a thoroughly gauge-ready manner, so that any spatial and temporal gauge conditions can be employed. The equations are completely general without any physical restriction except that we assume a flat homogeneous and isotropic universe as a background. We also comment briefly on the application of our formulation to the non-expanding Minkowski background.

We report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets.

High-multiplicity Kepler systems (referred to as Kepler multis) are often tightly packed and may be on the verge of instability. Many systems of this type could have experienced past instabilities, where the compact orbits and often low densities make physical collisions likely outcomes. We use numerical simulations to study the dynamical instabilities and planet-planet interactions in a synthetically generated sample of closely-packed, high-multiplicity systems. We focus specifically on systems resembling Kepler-11, a Kepler multi with six planets, and run a suite of dynamical integrations, sampling the initial orbital parameters around the nominal values reported in Lissauer et al. (2011a), finding that most of the realizations are unstable, resulting in orbit crossings and, eventually, collisions and mergers. We study in detail the dependence of stability on the orbital parameters of the planets and planet-pair characteristics to identify possible precursors to instability, compare the systems that emerge from dynamical instabilities to the observed Kepler sample (after applying observational corrections), and propose possible observable signatures of these instabilities. We examine the characteristics of each planet-planet collision, categorizing collisions by the degree of contact and collision energy, and find that grazing collisions are more common than direct impacts. Since the structure of many planets found in Kepler multis is such that the mass is dominated by a rocky core, but the volume is dominated by a low-density gaseous envelope, the sticky-sphere approximation may not be valid, and we present hydrodynamic calculations of planet-planet collisions clearly deviating from this approximation. Finally, we rerun a subset of our dynamical calculations using instead a modified prescription to handle collisions, finding, in general, higher multiplicity remnant systems.

We present hydrodynamic simulations of gas clouds inflowing from the disk to a few hundred parsec region of the Milky Way. A gravitational potential is generated to include realistic Galactic structures by using thousands of multipole expansions that describe 6.4 million stellar particles of a self-consistent Galaxy simulation. We find that a hybrid multipole expansion model, with two different basis sets and a thick disk correction, accurately reproduces the overall structures of the Milky Way. Through non-axisymmetric Galactic structures of an elongated bar and spiral arms, gas clouds in the disk inflow to the nuclear region and form a central molecular zone (CMZ)-like nuclear ring. We find that the size of the nuclear ring evolves into ~240 pc at T~1500 Myr, regardless of the initial size. For most simulation runs, the rate of gas inflow to the nuclear region is equilibrated to ~0.02 M_sun/yr. The nuclear ring is off-centered, relative to the Galactic center, by the lopsided central mass distribution of the Galaxy model, and thus an asymmetric mass distribution of the nuclear ring arises accordingly. The vertical asymmetry of the the Galaxy model also causes the nuclear ring to be tilted along the Galactic plane. During the first ~100 Myr, the vertical frequency of the gas motion is twice that of the orbital frequency, thus the projected nuclear ring shows a twisted, infinity-like shape.

The homogeneity of matter distribution at large scales, known as the cosmological principle, is a central assumption in the standard cosmological model. The case is testable though, thus no longer needs to be a principle. Here we perform a test for spatial homogeneity using the Sloan Digital Sky Survey Luminous Red Galaxies (LRG) sample by counting galaxies within a specified volume with the radius scale varying up to 300 Mpc/h. We directly confront the large-scale structure data with the definition of spatial homogeneity by comparing the averages and dispersions of galaxy number counts with allowed ranges of the random distribution with homogeneity. The LRG sample shows significantly larger dispersions of number counts than the random catalogues up to 300 Mpc/h scale, and even the average is located far outside the range allowed in the random distribution; the deviations are statistically impossible to be realized in the random distribution. This implies that the cosmological principle does not hold even at such large scales. The same analysis of mock galaxies derived from the N-body simulation, however, suggests that the LRG sample is consistent with the current paradigm of cosmology, thus the simulation is also not homogeneous in that scale. We conclude that the cosmological principle is not in the observed sky and nor is demanded to be there by the standard cosmological world model. This reveals the nature of the cosmological principle adopted in the modern cosmology paradigm, and opens new field of research in theoretical cosmology.

A tensor-type cosmological perturbation, defined as a transverse and traceless spatial fluctuation, is often interpreted as the gravitational waves. While decoupled from the scalar-type perturbations in linear order, the tensor perturbations can be sourced from the scalar-type in the nonlinear order. The tensor perturbations generated by the quadratic combination of linear scalar-type cosmological perturbation are widely studied in the literature, but all previous studies are based on zero-shear gauge without proper justification. Here, we show that, being second order in perturbation, such an induced tensor perturbation is generically gauge dependent. In particular, the gravitational wave power spectrum depends on the hypersurface (temporal gauge) condition taken for the linear scalar perturbation. We further show that, during the matter-dominated era, the induced tensor modes dominate over the linearly evolved primordial gravitational waves amplitude for $k\gtrsim10^{-2}~[h/{\rm Mpc}]$ even for the gauge that gives lowest induced tensor modes with the optimistic choice of primordial gravitational waves ($r=0.1$). The induced tensor modes, therefore, must be modeled correctly specific to the observational strategy for the measurement of primordial gravitational waves from large-scale structure via, for example, parity-odd mode of weak gravitational lensing, or clustering fossils.

The special relativistic hydrodynamics with weak gravity is hitherto unknown in the literature. Whether such an asymmetric combination is possible was unclear. Here, the hydrodynamic equations with Poisson-type gravity considering fully relativistic velocity and pressure under the weak gravity and the action-at-a-distance limit are consistently derived from Einstein's general relativity. Analysis is made in the maximal slicing where the Poisson's equation becomes much simpler than our previous study in the zero-shear gauge. Also presented is the hydrodynamic equations in the first post-Newtonian approximation, now under the {\it general} hypersurface condition. Our formulation includes the anisotropic stress.

We present the first VLBI detection of HCN molecular absorption in the nearby active galactic nucleus NGC 1052. Utilizing the 1 milliarcsecond resolution achieved by the Korean VLBI Network, we have spatially resolved the HCN absorption against a double-sided nuclear jet structure. Two velocity features of HCN absorption are detected significantly at the radial velocity of 1656 and 1719 km/s, redshifted by 149 and 212 km/s with respect to the systemic velocity of the galaxy. The column density of the HCN molecule is estimated to be 10^{15}-10^{16} cm^{-2}, assuming the excitation temperature of 100-230 K. The absorption features show high optical depth localized on the receding jet side, where the free-free absorption occurred due to the circumnuclear torus. The size of the foreground absorbing molecular gas is estimated to be on approximately one-parsec scales, which agrees well with the approximate size of the circumnuclear torus. HCN absorbing gas is likely to be several clumps smaller than 0.1 parsec inside the circumnuclear torus. The redshifted velocities of the HCN absorption features imply that HCN absorbing gas traces ongoing infall motion inside the circumnuclear torus onto the central engine.

We present the leading order non-linear density and velocity power spectra in the complete form; previous studies have omitted the vector- and tensor-type perturbations simultaneously excited by the scalar-type perturbation in nonlinear order. These additional contributions are comparable to the scalar-type purely relativistic perturbations, and thus negligible in the current paradigm of concordance cosmology: i.e., concerning density and velocity perturbations of the pressureless matter in perturbation regime well inside of matter-dominated epoch, we show that pure Einstein's gravity contributions appearing from the third order are entirely negligible (five orders of magnitude smaller than the Newtonian contributions) in all scales. We thus prove that Newtonian perturbation theory is quite reliable in calculating the amplitude of matter fluctuations even in the precision era of cosmology. Therefore, besides the ones imprinted as the initial condition generated in the earlier phase, the other relativistic effect relevant for interpreting observational data must be the projection effect that occurs when mapping galaxies on to the observed coordinate.

We present a new approximation to include fully general relativistic pressure and velocity in Newtonian hydrodynamics. The energy conservation, momentum conservation and two Poisson's equations are consistently derived from Einstein's gravity in the zero-shear gauge assuming weak gravity and action-at-a-distance limit. The equations show proper special relativity limit in the absence of gravity. Our approximation is complementary to the post-Newtonian approximation and the equations are valid in fully nonlinear situations.

We present fully nonlinear and exact cosmological perturbation equations in the presence of multiple components of fluids and minimally coupled scalar fields. We ignore the tensor-type perturbation. The equations are presented without taking the temporal gauge condition in the Friedmann background with general curvature and the cosmological constant. For each fluid component we ignore the anisotropic stress. The multiple component nature, however, introduces the anisotropic stress in the collective fluid quantities. We prove the Newtonian limit of multiple fluids in the zero-shear gauge and the uniform-expansion gauge conditions, present the Newtonian hydrodynamic equations in the presence of general relativistic pressure in the zero-shear gauge, and present the fully nonlinear equations and the third-order perturbation equations of the nonrelativistic pressure fluids in the CDM-comoving gauge.

Newtonian cosmological perturbation equations valid to full nonlinear order are well known in the literature. Assuming the absence of the transverse-tracefree part of the metric, we present the general relativistic counterpart valid to full nonlinear order. The relativistic equations are presented without taking the slicing (temporal gauge) condition. The equations do have the proper Newtonian and first post-Newtonian limits. We also present the relativistic pressure correction terms in the Newtonian hydrodynamic equations.

We investigate aspects of axion as a coherently oscillating massive classical scalar field by analyzing third order perturbations in Einstein's gravity in the axion-comoving gauge. The axion fluid has its characteristic pressure term leading to an axion Jeans scale which is cosmologically negligible for a canonical axion mass. Our classically derived axion pressure term in Einstein's gravity is identical to the one derived in the non-relativistic quantum mechanical context in the literature. We show that except for the axion pressure term, the axion fluid equations are exactly the same as the general relativistic continuity and Euler equations of a zero-pressure fluid up to third order perturbation. The general relativistic density and velocity perturbations of the CDM in the CDM-comoving gauge are exactly the same as the Newtonian perturbations to the second order (in all scales), and the pure general relativistic corrections appearing from the third order are numerically negligible (in all scales as well) in the current paradigm of concordance cosmology. Therefore, here we prove that, in the super-Jeans scale, the classical axion can be handled as the Newtonian CDM fluid up to third order perturbation. We also show that the axion fluid supports the vector-type (rotational) perturbation from the third order. Our analysis includes the cosmological constant.

Massive stars are usually found in binaries, and binaries with periods less than 10 days may have a preference for near equal component masses. In this paper we investigate the evolution of these binaries all the way to contact and the possibility that these systems can be progenitors of double neutron star binaries. The small orbital separations of observed double neutron star binaries suggest that the progenitor systems underwent a common envelope phase at least once during their evolution. Bethe & Brown (1998) proposed that massive binary twins will undergo a common envelope evolution while both components are ascending the red giant branch or asymptotic giant branch simultaneously, also known as double-core evolution. Using models generated from the stellar evolution code Evolve Zero Age Main Sequence, we determine the range of mass ratios resulting in both components simultaneously ascending the RGB or AGB as a function of the difference in birth times, t. We find that, even for a generous t=5 Myr, the minimum mass ratio qmin=0.933 for an 8 Solar Mass primary and increases for larger primaries. We use a hydrodynamics code, StarSmasher, to study specifically the evolution of q=1 common envelope systems as a function of initial component mass, age, and orbital separation. We find the dynamical stability limit, the largest orbital separation where the binary becomes dynamically unstable, as a function of the component mass and age. Finally, we calculate the efficiency of ejecting matter during the inspiral phase to extrapolate the properties of the remnant binary from our numerical results, assuming the common envelope is completely ejected. We find that for the nominal core masses, there is a minimum orbital separation for a given component mass such that the helium cores survive common envelope evolution in a tightly bound binary and are viable progenitors for double neutron stars.

We use $N$-body/smoothed particle hydrodynamics simulations of encounters between an early-type galaxy (ETG) and a late-type galaxy (LTG) to study the effects of hot halo gas on the evolution for a case with the mass ratio of the ETG to LTG of 2:1 and the closest approach distance of $\sim$100 kpc. We find that the dynamics of the cold disk gas in the tidal bridge and the amount of the newly formed stars depend strongly on the existence of a gas halo. In the run of interacting galaxies not having a hot gas halo, the gas and stars accreted into the ETG do not include newly formed stars. However, in the run using the ETG with a gas halo and the LTG without a gas halo, a shock forms along the disk gas tidal bridge and induces star formation near the closest approach. The shock front is parallel to a channel along which the cold gas flows toward the center of the ETG. As a result, the ETG can accrete star-forming cold gas and newly born stars at and near its center. When both galaxies have hot gas halos, a shock is formed between the two gas halos somewhat before the closest approach. The shock hinders the growth of the cold gas bridge to the ETG and also ionizes it. Only some of the disk stars transfer through the stellar bridge. We conclude that the hot halo gas can give significant hydrodynamic effects during distant encounters.

We investigate observational consequences of the early episodically dominating dark energy on the evolution of cosmological structures. For this aim, we introduce the minimally coupled scalar field dark energy model with the Albrecht-Skordis potential which allows a sudden ephemeral domination of dark energy component during the radiation or early matter era. The conventional cosmological parameters in the presence of such an early dark energy are constrained with WMAP and Planck cosmic microwave background radiation data including other external data sets. It is shown that in the presence of such an early dark energy the estimated cosmological parameters can deviate substantially from the currently known $\Lambda\textrm{CDM}$-based parameters, with best-fit values differing by several percents for WMAP and by a percent level for Planck data. For the latter case, only a limited amount of dark energy with episodic nature is allowed since the Planck data strongly favors the $\Lambda\textrm{CDM}$ model. Compared with the conventional dark energy model, the early dark energy dominating near radiation-matter equality or at the early matter era results in the shorter cosmic age or the presence of tensor-type perturbation, respectively. Our analysis demonstrates that the alternative cosmological parameter estimation is allowed based on the same observations even in Einstein's gravity.

We calculate and analyze the distribution of period ratios observed in systems of Kepler exoplanet candidates including studies of both adjacent planet pairs and all planet pairs. These distributions account for both the geometrical bias against detecting more distant planets and the effects of incompleteness due to planets missed by the data reduction pipeline. In addition to some of the known features near first-order mean-motion resonances (MMR), there is a significant excess of planet pairs with period ratios near 2.2. The statistical significance of this feature is assessed using Monte Carlo simulation. We also investigate the distribution of period ratios near first-order MMR and compare different quantities used to measure this distribution. We find that beyond period ratios of ~2.5, the distribution of all period ratios follows a power-law with an exponent -1.26 +/- 0.05. We discuss implications that these results may have on the formation and dynamical evolution of Kepler-like planetary systems---systems of sub-Neptune/super-Earth planets with relatively short orbital periods.

We study the non-linear corrections to the matter and velocity power spectra in the synchronous gauge (SG). We consider the perturbations up to third order in a zero-pressure fluid in flat cosmological background, which is relevant for the non-linear growth of cosmic structure. As a result, we point out that the SG is an inappropriate coordinate choice when handling the non-linear growth of the large-scale structure. Although the equations in the SG happen to coincide with those in the comoving gauge (CG) to linear order, they differ from second order. In particular, the second order hydrodynamic equations in the the SG are apparently in the Lagrangian form, whereas those in the CG are in the Eulerian form. Thus, the non-linear power spectra naively presented in the original SG show strange behavior quite different from the result of the Newtonian theory even on sub-horizon scales. The power spectra in the SG show regularized behaviors only after we introduce convective terms in the second order so that the equations in two gauges coincide to the second order.

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

We present the general relativistic pressure correction terms in Newtonian hydrodynamic equations to the nonlinear order: these are equations (\ref{mass-conservation-Mink})-(\ref{Poisson-eq-Mink}). The derivation is made in the zero-shear gauge based on the fully nonlinear formulation of cosmological perturbation in Einstein's gravity. The correction terms {\it differ} from many of the previously suggested forms in the literature based on hand-waving manners. We confirm our results by comparing with (i) the nonlinear perturbation theory, (ii) the first order post-Newtonian approximation, and (iii) the special relativistic limit, and by checking (iv) the consistency with full Einstein's equation.

We prove that the axion as a coherently oscillating scalar field acts as a cold dark matter (CDM) to the second-order perturbations in all cosmological scales including the super-horizon scale. The proof is made in the axion-comoving gauge. For a canonical mass, the axion pressure term causes deviation from the CDM only on scales smaller than the Solar System size. Beyond such a small scale the equations of the axion fluid are the same as the ones of the CDM based on the CDM-comoving gauge which are exactly identical to the Newtonian equations to the second order. We also show that the axion fluid does not generate the rotational (vector-type) perturbation even to the second order. Thus, in the case of axion fluid, we have the relativistic/Newtonian correspondence to the second order, even considering the rotational perturbation. Our analysis is made in the presence of the cosmological constant, and can be easily extended to the realistic situation including other components of fluids and fields.

We derive the basic equations of the cosmological first-order post-Newtonian approximation from the recently formulated fully nonlinear and exact cosmological perturbation theory in Einstein's gravity. Apparently the latter, being exact, should include the former, and here we use this fact as a new derivation of the former. The complete sets of equations in both approaches are presented without fixing the temporal gauge conditions so that we can use the gauge choice as an advantage. Comparisons between the two approaches are made. Both are potentially important in handling relativistic aspects of nonlinear processes occurring in cosmological structure formation. We consider an ideal fluid and include the cosmological constant.

We prove that in the infinite speed-of-light limit (i.e., non-relativistic and subhorizon limits), the relativistic fully nonlinear cosmological perturbation equations in two gauge conditions, the zero-shear gauge and the uniform-expansion gauge, exactly reproduce the Newtonian hydrodynamic perturbation equations in the cosmological background; as a consequence, in the same two gauge conditions, the Newtonian hydrodynamic equations are exactly recovered in the Minkowsky background.

In 1988 Bardeen has suggested a pragmatic formulation of cosmological perturbation theory which is powerful in practice to employ various fundamental gauge conditions easily depending on the character of the problem. The perturbation equations are presented without fixing the temporal gauge condition and are arranged so that one can easily impose fundamental gauge conditions by simply setting one of the perturbation variables in the equations equal to zero. In this way one can use the gauge degrees of freedom as an advantage in handling problems. Except for the synchronous gauge condition, all the other fundamental gauge conditions completely fix the gauge mode, and consequently, each variable in such a gauge has a unique gauge invariant counterpart, so that we can identify the variable as the gauge-invariant one. Here, we extend Bardeen's linear formulation to fully nonlinear order in perturbations, with the gauge advantage kept intact. Derived equations are exact, and from these we can easily expand to higher order perturbations in a gauge-ready form. We consider scalar- and vector-type perturbations of an ideal fluid in a flat background; we also present the multiple components of ideal fluid case. As applications we present fully nonlinear density and velocity perturbation equations in Einstein's gravity in the zero-pressure medium, vorticity generation from pure scalar-type perturbation, and fluid formulation of a minimally coupled scalar field, all in the comoving gauge. We also present the equation of gravitational waves generated from pure scalar- and vector-type perturbations.

We estimate the power spectra of CMB temperature anisotropy in localized regions on the sky using the WMAP 7-year data. Here, we report that the north hat and the south hat regions at the high Galactic latitude (|b| >= 30 deg) show anomaly in the power spectrum amplitude around the third peak, which is statistically significant up to 3 sigma. We try to figure out the cause of the observed anomaly by analyzing the low Galactic latitude (|b|< 30 deg) regions where the galaxy contamination is expected to be stronger, and regions that are weakly or strongly dominated by the WMAP instrument noise. We also consider the possible effect of unresolved radio point sources. We found another but less statistically significant anomaly in the low Galactic latitude north and south regions whose behavior is opposite to the one at the high latitude. Our analysis shows that the observed north-south anomaly at high latitude becomes weaker on the regions with high number of observations (weak instrument noise), suggesting that the anomaly is significant at sky regions that are dominated by the WMAP instrument noise. We have checked that the observed north-south anomaly has weak dependences on the bin-width used in the power spectrum estimation and the Galactic latitude cut. We have also discussed the possibility that the detected anomaly may hinge on the particular choice of the multipole bin around the third peak. We anticipate that the issue of whether the anomaly is intrinsic one or due to the WMAP instrument noise will be resolved by the forthcoming Planck data.

We construct several Milky Way-like galaxy models containing a gas halo (as well as gaseous and stellar disks, a dark matter halo, and a stellar bulge) following either an isothermal or an NFW density profile with varying mass and initial spin. In addition, galactic winds associated with star formation are tested in some of the simulations. We evolve these isolated galaxy models using the GADGET-3 $N$-body/hydrodynamic simulation code, paying particular attention to the effects of the gas halo on the evolution. We find that the evolution of the models is strongly affected by the adopted gas halo component. The model without a gas halo shows an increasing star formation rate (SFR) at the beginning of the simulation for some hundreds of millions of years and then a continuously decreasing rate to the end of the run at 3 Gyr. On the other hand, the SFRs in the models with a gas halo emerge to be either relatively flat throughout the simulations or increasing over a gigayear and then decreasing to the end. The models with the more centrally concentrated NFW gas halo show overall higher SFRs than those with the isothermal gas halo of the equal mass. The gas accretion from the halo onto the disk also occurs more in the models with the NFW gas halo, however, this is shown to take place mostly in the inner part of the disk and not to contribute significantly to the star formation unless the gas halo has very high density at the central part. The rotation of a gas halo is found to make SFR lower in the model. The SFRs in the runs including galactic winds are found to be lower than the same runs but without winds. We conclude that the effects of a hot gaseous halo on the evolution of galaxies are generally too significant to be simply ignored, and expect that more hydrodynamical processes in galaxies could be understood through numerical simulations employing both gas disk and gas halo components.

Axion as a coherently oscillating scalar field is known to behave as a cold dark matter in all cosmologically relevant scales. For conventional axion mass with 10^{-5} eV, the axion reveals a characteristic damping behavior in the evolution of density perturbations on scales smaller than the solar system size. The damping scale is inversely proportional to the square-root of the axion mass. We show that the axion mass smaller than 10^{-24} eV induces a significant damping in the baryonic density power spectrum in cosmologically relevant scales, thus deviating from the cold dark matter in the scale smaller than the axion Jeans scale. With such a small mass, however, our basic assumption about the coherently oscillating scalar field is broken in the early universe. This problem is shared by other dark matter models based on the Bose-Einstein condensate and the ultra-light scalar field. We introduce a simple model to avoid this problem by introducing evolving axion mass in the early universe, and present observational effects of present-day low-mass axion on the baryon density power spectrum, the cosmic microwave background radiation (CMB) temperature power spectrum, and the growth rate of baryon density perturbation. In our low-mass axion model we have a characteristic small-scale cutoff in the baryon density power spectrum below the axion Jeans scale. The small-scale deviations from the cold dark matter model in both matter and CMB power spectra clearly differ from the ones expected in the cold dark matter model mixed with the massive neutrinos as a hot dark matter component.

Physical cosmology tries to understand the Universe at large with its origin and evolution. Observational and experimental situations in cosmology do not allow us to proceed purely based on the empirical means. We examine in which sense our cosmological assumptions in fact have shaped our current cosmological worldview with consequent inevitable limits. Cosmology, as other branches of science and knowledge, is a construct of human imagination reflecting the popular belief system of the era. The question at issue deserves further philosophic discussions. In Whitehead's words, "philosophy, in one of its functions, is the critic of cosmologies". (Whitehead 1925)

We compare the cosmological first-order post-Newtonian (1PN) approximation with the relativistic cosmological linear perturbation theory in a zero-pressure medium with the cosmological constant. We compare equations and solutions in several different gauge conditions available in both methods. In the PN method we have perturbation equations for density, velocity and gravitational potential independently of the gauge condition to 1PN order. However, correspondences with these 1PN equations are available only in certain gauge conditions in the perturbation theory. Equations of perturbed velocity and the perturbed gravitational potential in the zero-shear gauge exactly coincide with the Newtonian equations which remain valid even to 1PN order (the same is true for perturbed velocity in the comoving gauge), and equations of perturbed density in the zero-shear gauge and the uniform-expansion gauge coincide to 1PN order. We identify other correspondences available in different gauge conditions of the perturbation theory.

We present the growing mode solutions of cosmological perturbations to the second order in the matter dominated era. We also present several gauge-invariant combinations of perturbation variables to the second order in most general fluid context. Based on the solutions we study the Newtonian correspondence of relativistic perturbations to the second order. In addition to the previously known exact relativistic/Newtonian correspondence of density and velocity perturbations to the second order in the comoving gauge, here we show that in the sub-horizon limit we have the correspondences for density, velocity and potential perturbations in the zero-shear gauge and in the uniform-expansion gauge to the second order. Density perturbation in the uniform curvature gauge also shows the correspondence to the second order in the sub-horizon scale. We also identify the relativistic gravitational potential which shows exact correspondence to the Newtonian one to the second order.

Cold dark matter in the Milky Way halo may have structure defined by flows with low velocity dispersion. The Axion Dark Matter eXperiment high resolution channel is especially sensitive to axions in such low velocity dispersion flows. Results from a combined power spectra analysis of the high resolution channel axion search are presented along with a discussion of the assumptions underlying such an analysis. We exclude KSVZ axion dark matter densities of {\rho} > 0.2 GeV/cm^3 and DFSZ densities of {\rho} > 1.4 GeV/cm^3 over a mass range of m_a = 3.3{\mu}eV to 3.69{\mu}eV for models having velocity dispersions of {\Delta}{\beta} < 3x10^-6.

We present smoothed particle hydrodynamic models of the interactions in the compact galaxy group, Stephan's Quintet. This work is extension of the earlier collisionless N-body simulations of Renaud et al. in which the large-scale stellar morphology of the group was modeled with a series of galaxy-galaxy interactions in the simulations. Including thermohydrodynamic effects in this work, we further investigate the dynamical interaction history and evolution of the intergalactic gas of Stephan's Quintet. The major features of the group, such as the extended tidal features and the group-wide shock, enabled us to constrain the models reasonably well, while trying to reproduce multiple features of the system. We found that reconstructing the two long tails extending from NGC 7319 toward NGC 7320c one after the other in two separate encounters is very difficult and unlikely, because the second encounter usually destroys or distorts the already-generated tidal structure. Our models suggest the two long tails may be formed simultaneously from a single encounter between NGC 7319 and 7320c, resulting in a thinner and denser inner tail than the outer one. The tails then also run parallel to each other as observed. The model results support the ideas that the group-wide shock detected in multi-wavelength observations between NGC 7319 and 7318b and the starburst region north of NGC 7318b are triggered by the high-speed collision between NGC 7318b and the intergalactic gas. Our models show that a gas bridge is formed by the high-speed collision and clouds in the bridge continue to interact for some tens of millions of years after the impact. This produces many small shocks in that region, resulting a much longer cooling time than that of a single impact shock.

We study the non-linear evolution of the curvature perturbations during matter dominated era. We show that regardless of the origin of the primordial perturbation, the Bardeen potential and curvature receive sizable contributions from the classical non-linear evolution effects, and quantify them exactly. On the super-horizon scales we have squeezed peak of the bispectrum with magnitude, in terms of the local non-linear parameters of Bardeen curvature, 1/6 < f_{NL} 19/15, and of Bardeen potential, -1/4 < f_{NL} < 7/5, depending on the configuration of momenta. On the sub-horizon scales the bispectrum show equilateral shape, and can serve as a potential probe of general relativity.

We study characters of recent type Ia supernova (SNIa) data using evolving dark energy models with changing equation of state parameter w. We consider sudden-jump approximation of w for some chosen redshift spans with double transitions, and constrain these models based on Markov Chain Monte Carlo (MCMC) method using the SNIa data (Constitution, Union, Union2) together with baryon acoustic oscillation A parameter and cosmic microwave background shift parameter in a flat background. In the double-transition model the Constitution data shows deviation outside 1 sigma from LCDM model at low (z < 0.2) and middle (0.2 < z < 0.4) redshift bins whereas no such deviations are noticeable in the Union and Union2 data. By analyzing the Union members in the Constitution set, however, we show that the same difference is actually due to different calibration of the same Union sample in the Constitution set, and is not due to new data added in the Constitution set. All detected deviations are within 2 sigma from the LCDM world model. From the LCDM mock data analysis, we quantify biases in the dark energy equation of state parameters induced by insufficient data with inhomogeneous distribution of data points in the redshift space and distance modulus errors. We demonstrate that location of peak in the distribution of arithmetic means (computed from the MCMC chain for each mock data) behaves as an unbiased estimator for the average bias, which is valid even for non-symmetric likelihood distributions.

We show that in the single component situation all perturbation variables in the comoving gauge are conformally invariant to all perturbation orders. Generally we identify a special time slicing, the uniform-conformal transformation slicing, where all perturbations are again conformally invariant to all perturbation orders. We apply this result to the delta N formalism, and show its conformal invariance.

We have discovered long-lived waves in two sets of numerical models of fast (marginally bound or unbound) flyby galaxy collisions, carried out independently with two different codes. In neither simulation set are the spirals the result of a collision-induced bar formation. Although there is variation in the appearance of the waves with time, they do not disappear and reform recurrently, as seen in other cases described in the literature. We also present an analytic theory that can account for the wave structure, not as propagating transients, nor as a fixed pattern propagating through the disc. While these waves propagate through the disc, they are maintained by the coherent oscillations initiated by the impulsive disturbance. Specifically, the analytic theory suggests that they are caustic waves in ensembles of stars pursuing correlated epicyclic orbits after the disturbance. This theory is an extension of that developed by Struck and collaborators for colliding ring galaxies. The models suggest that this type of wave may persist for a couple of Gyr., and galaxy interactions occur on comparable timescales, so waves produced by the mechanism may be well represented in observed spirals. In particular, this mechanism can account for the tightly wound, and presumably long-lived spirals, seen in some nearby early-type galaxies. These spirals are also likely to be common in groups and clusters, where fast encounters between galaxies occur relatively frequently. However, as the spirals become tightly wound, and evolve to modest amplitudes, they may be difficult to resolve unless they are nearby. Nonetheless, the effect may be one of several processes that result from galaxy harassment, and via wave-enhanced star formation contribute to the Butcher-Oemler effect.

We study non-linear contributions to the power spectrum of the curvature perturbation on super-horizon scales, produced during slow-roll inflation driven by a canonical single scalar field. We find that on large scales the linear power spectrum completely dominates and leading non-linear corrections remain totally negligible, indicating that we can safely rely on linear perturbation theory to study inflationary power spectrum. We also briefly comment on the infrared and ultraviolet behaviour of the non-linear corrections.

The modified gravity with $f(R)=R^{1+\epsilon}$ ($\epsilon>0$) allows a scaling solution where the density of gravity sector follows the density of the dominant fluid. We present initial conditions of background and perturbation variables during the scaling evolution regime in the modified gravity. As a possible dark energy model we consider a gravity with a form $f(R)=R^{1+\epsilon}+qR^{-n}$ ($-1<n \le 0$) where the second term drives the late-time acceleration. We show that our $f(R)$ gravity parameters are very sensitive to the baryon perturbation growth and baryon density power spectrum, and present observational constraints on the model parameters. Our analysis suggests that only the parameter space extremely close to the $\Lambda\textrm{CDM}$ model is allowed.

Versions of parameterized pseudo-Newtonian gravity theories specially designed for cosmology have been introduced in recent cosmology literature. The modifications demand a zero-pressure fluid in the context of versions of modified Poisson-like equation with two different gravitational potentials. We consider such modifications in the context of relativistic gravity theories where the action is a general algebraic function of the scalar curvature, the scalar field, and the kinetic term of the field. In general it is not possible to isolate the zero-pressure fluid component simultaneously demanding a modification in the Poisson-like equation. Only in the small-scale limit we can realize some special forms of the attempted modifications. We address some loopholes in the possibility of showing non-Einstein gravity nature based on pseudo-Newtonian modifications in the cosmological context. We point out that future observations of gravitational weak lensing together with velocity perturbation can potentially test the validity of Einstein's gravity in cosmology context.

Non-linear nature of Einstein equation introduces genuine relativistic higher order corrections to the usual Newtonian fluid equations describing the evolution of cosmological perturbations. We study the effect of such novel non-linearities on the next-to-leading order matter and velocity power spectra for the case of pressureless, irrotational fluid in a flat Friedmann background. We find that pure general relativistic corrections are negligibly small over all scales. Our result guarantees that, in the current paradigm of standard cosmology, one can safely use Newtonian cosmology even in non-linear regimes.

Using the SEDs of the weak AGNs 35 LINERs presented in a companion paper, we assess whether photoionization by the weak AGN can power the emission-line luminosities measured through the large (few-arcsecond) apertures used in ground-based spectroscopic surveys. Spectra taken through such apertures are used to define LINERs as a class and constrain non-stellar photoionization models for LINERs. Therefore, our energy budget test is a self-consistency check of the idea that the observed emission lines are powered by an AGN. We determine the ionizing luminosities and photon rates by integrating the observed SEDs and by scaling a template SED. Even if all ionizing photons are absorbed by the line-emitting gas, more than half of our LINERs suffer from a deficit of ionizing photons. In 1/3 of LINERs the deficit is severe. If only 10% of the ionizing photons are absorbed by the gas, there is an ionizing photon deficit in 85% of LINERs. We disfavor the possibility that additional electromagnetic power, either obscured or emitted in the unobservable far-UV band, is available from the AGN. We consider other power sources such as mechanical heating by compact jets and photoionization by young or old stars. Photoionization by young stars may be important in a small fraction of cases. Mechanical heating provides enough power in most cases but it is not clear how this power is transferred to the emission-line gas. Photoionization by post-AGB stars is an important power source; it provides more ionizing photons that the AGN in more than half of the LINERs and enough ionizing photons to power the emission lines in 1/3 of the LINERs. It appears likely that the emission-line spectra of LINERs obtained from the ground include the sum of emission from different regions where different power sources dominate.

We present a compilation of spectral energy distributions of 35 weak AGNs in LINERs using recent data from the published literature. We make use of previously published compilations of data, after complementing and extending them with more recent data. The main improvement in the recent data is afforded by high-spatial resolution observations with the Chandra X-Ray Observatory and high-spatial resolution radio observations utilizing a number of facilities. In addition, a considerable number of objects have been observed with the HST in the near-IR through near-UV bands since the earlier compilations were published. The data include upper limits resulting from either non-detections or observations at low spatial resolution that do not isolate the AGN. For the sake of completeness, we also compute and present a number of quantities from the data, such as alpha-ox, bolometric corrections, bolometric luminosities, Eddington ratios, and the average SED. We anticipate that these data will be useful for a number of applications. In a companion paper, we use a subset of these data ourselves to assess the energy budgets of LINERs.

Axions in the micro eV mass range are a plausible cold dark matter candidate and may be detected by their conversion into microwave photons in a resonant cavity immersed in a static magnetic field. The first result from such an axion search using a superconducting first-stage amplifier (SQUID) is reported. The SQUID amplifier, replacing a conventional GaAs field-effect transistor amplifier, successfully reached axion-photon coupling sensitivity in the band set by present axion models and sets the stage for a definitive axion search utilizing near quantum-limited SQUID amplifiers.

We show the importance of properly including the perturbations of the dark energy component in the dynamical dark energy models based on a scalar field and modified gravity theories in order to meet with present and future observational precisions. Based on a simple scaling scalar field dark energy model, we show that observationally distinguishable substantial differences appear by ignoring the dark energy perturbation. By ignoring it the perturbed system of equations becomes inconsistent and deviations in (gauge-invariant) power spectra depend on the gauge choice.

We study a generalized version of Chaplygin gas as unified model of dark matter and dark energy. Using realistic theoretical models and the currently available observational data from the age of the universe, the expansion history based on the type Ia supernovae, the matter power spectrum, the cosmic microwave background radiation anisotropy power spectra, and the perturbation growth factor we put the unified model under observational test. As the model has only two free parameters in the flat Friedmann background [$\Lambda$CDM (cold dark matter) model has only one free parameter] we show that the model is already tightly constrained by currently available observations. The only parameter space extremely close to the $\Lambda$CDM model is allowed in this unified model.

We use smoothed particle hydrodynamics (SPH) models to study the large-scale morphology and dynamical evolution of the intergalactic gas in Stephan's Quintet, and compare to multiwavelength observations. Specifically, we model the formation of the hot X-ray gas, the large-scale shock, and emission line gas as the result of NGC 7318b colliding with the group. We also reproduce the N-body model of Renaud and Appleton for the tidal structures in the group.

We probe the pure Einstein's gravity contributions to the second-order density power spectrum. In the small-scale, we discover that the Einstein's gravity contribution is negligibly small. This guarantees that Newton's gravity is sufficient to handle the baryon acoustic oscillation scale. In the large scale, however, we discover that the Einstein's gravity contribution to the second-order power spectrum dominates the linear-order power spectrum. Thus, pure Einstein gravity contribution appearing in the third-order perturbation leads to an infrared divergence in the power spectrum.

Here we generally prove that the axion as a coherently oscillating scalar field acts as a cold dark matter in nearly all cosmologically relevant scales. The proof is made in the linear perturbation order. Compared with our previous proof based on solutions, here we compare the equations in the axion with the ones in the cold dark matter, thus expanding the valid range of the proof. Deviation from purely pressureless medium appears in very small scale where axion reveals a peculiar equation of state. Our analysis is made in the presence of the cosmological constant, and our conclusions are valid in the presence of other fluid and field components.

Is Newton's gravity sufficient to handle the weakly nonlinear evolution stages of the cosmic large-scale structures? Here we resolve the issue by analytically deriving the density and velocity power spectra to the second order in the context of Einstein's gravity. The recently found pure general relativistic corrections appearing in the third-order perturbation contribute to power spectra to the second order. In this work the complete density and velocity power spectra to the second order are derived. The power transfers among different scales in the density power spectrum are estimated in the context of Einstein's gravity. The relativistic corrections in the density power spectrum are estimated to be smaller than the Newtonian one to the second order, but these could be larger than higher-order nonlinear Newtonian terms.

We present Spitzer infrared, GALEX UV, and SDSS and SARA optical images of the peculiar interacting galaxy pair Arp 285 (NGC 2856/4), and compare with a new numerical model of the interaction. We estimate the ages of clumps of star formation in these galaxies using population synthesis models, carefully considering the uncertainties on these ages. This system contains a striking example of `beads on a string': a series of star formation complexes ~1 kpc apart. These `beads' are found in a tail-like feature that is perpendicular to the disk of NGC 2856, which implies that it was formed from material accreted from the companion NGC 2854. The extreme blueness of the optical/UV colors and redness of the mid-infrared colors implies very young stellar ages (~4 - 20 Myrs) for these star forming regions. Spectral decomposition of these `beads' shows excess emission above the modeled stellar continuum in the 3.6 micron and 4.5 micron bands, indicating either contributions from interstellar matter to these fluxes or a second older stellar population. These clumps have -12.0 < M(B) < -10.6, thus they are less luminous than most dwarf galaxies. Our model suggests that bridge material falling into the potential of the companion overshoots the companion. The gas then piles up at apo-galacticon before falling back onto the companion, and star formation occurs in the pile-up. A luminous (M(B) ~ -13.6) extended (FWHM ~ 1.3 kpc) `bright spot' is visible at the northwestern edge of the NGC 2856 disk, with an intermediate stellar population (400 - 1500 Myrs). Our model suggests that this feature is part of a expanding ripple-like `arc' created by an off-center ring-galaxy-like collision between the two disks.

Present expansion stage of the universe is believed to be mainly governed by the cosmological constant, collisionless dark matter and baryonic matter. The latter two components are often modeled as zero-pressure fluids. In our previous work we have shown that to the second-order cosmological perturbations, the relativistic equations of the zero-pressure, irrotational, multi-component fluids in a spatially near flat background effectively coincide with the Newtonian equations. As the Newtonian equations only have quadratic order nonlinearity, it is practically interesting to derive the potential third-order perturbation terms in general relativistic treatment which correspond to pure general relativistic corrections. Here, we present pure general relativistic correction terms appearing in the third-order perturbations of the multi-component zero-pressure fluids. We show that, as in a single component situation, the third-order correction terms are quite small (~ 5 x10^{-5} smaller compared with the relativistic/Newtonian second-order terms) due to the weak level anisotropy of the cosmic microwave background radiation. Still, there do exist pure general relativistic correction terms in third-order perturbations which could potentially become important in future development of precision cosmology. We include the cosmological constant in all our analyses.

We present general relativistic correction terms appearing in Newton's gravity to the second-order perturbations of cosmological fluids. In our previous work we have shown that to the second-order perturbations, the density and velocity perturbation equations of general relativistic zero-pressure, irrotational, single-component fluid in a flat background coincide exactly with the ones known in Newton's theory. Here, we present the general relativistic second-order correction terms arising due to (i) pressure, (ii) multi-component, (iii) background curvature, and (iv) rotation. In case of multi-component zero-pressure, irrotational fluids under the flat background, we effectively do not have relativistic correction terms, thus the relativistic result again coincides with the Newtonian ones. In the other three cases we generally have pure general relativistic correction terms. In case of pressure, the relativistic corrections appear even in the level of background and linear perturbation equations. In the presence of background curvature, or rotation, pure relativistic correction terms directly appear in the Newtonian equations of motion of density and velocity perturbations to the second order. In the small-scale limit (far inside the horizon), relativistic equations including the rotation coincide with the ones in Newton's gravity.

We present a detailed numerical study of the evolutions of cosmological linear perturbations through a simple bouncing world model based on two scalar fields. We properly identify the relatively growing and decaying solutions in expanding and collapsing phases. Using a decomposition based on the large-scale limit exact solution of curvature (adiabatic) perturbations with two independent modes, we assign the relatively growing/decaying one in an expanding phase as the C/d-mode. In the collapsing phase, the roles are reversed, and the C/d-mode is relatively decaying/growing. The analytic solution shows that, as long as the large scale and the adiabatic conditions are met, the C- and d-modes preserve their nature throughout the bounce. Here, by using a concrete nonsingular bouncing world model based on two scalar fields, we numerically follow the evolutions of the correctly identified C- and d-modes which preserve their nature through the bounce, thus confirming our previous anticipation based on the analytic solution. Thus, while the large-scale condition is satisfied and the adiabatic condition is met during the bounce, we conclude that it is not possible to obtain the near Harrison-Zel'dovich scale-invariant density spectrum through a bouncing world model as long as the seed fluctuations were generated from quantum fluctuations of the curvature perturbation in the collapsing phase.

Except for the presence of gravitational wave source term, the relativistic perturbation equations of a zero-pressure irrotational fluid in a flat Friedmann world model coincide exactly with the Newtonian ones to the second order in perturbations. Such a relativistic-Newtonian correspondence is available in a special gauge condition (the comoving gauge) in which all the variables are equivalently gauge invariant. In this work we compare our results with the ones in the synchronous gauge which has been used often in the literature. Although the final equations look simpler in the synchronous gauge, the variables have remnant gauge modes. Except for the presence of the gauge mode for the perturbed order variables, however, the equations in the synchronous gauge are gauge invariant and can be exactly identified as the Newtonian hydrodynamic equations in the Lagrangian frame. In this regard, the relativistic equations to the second order in the comoving gauge are the same as the Newtonian hydrodynamic equations in the Eulerian frame. We resolve several issues related to the two gauge conditions often to fully nonlinear orders in perturbations.

Both for the background world model and its linear perturbations Newtonian cosmology coincides with the zero-pressure limits of relativistic cosmology. However, such successes in Newtonian cosmology are not purely based on Newton's gravity, but are rather guided ones by previously known results in Einstein's theory. The action-at-a-distance nature of Newton's gravity requires further verification from Einstein's theory for its use in the large-scale nonlinear regimes. We study the domain of validity of the Newtonian cosmology by investigating weakly nonlinear regimes in relativistic cosmology assuming a zero-pressure and irrotational fluid. We show that, first, if we ignore the coupling with gravitational waves the Newtonian cosmology is exactly valid even to the second order in perturbation. Second, the pure relativistic correction terms start appearing from the third order. Third, the correction terms are independent of the horizon scale and are quite small in the large-scale near the horizon. These conclusions are based on our special (and proper) choice of variables and gauge conditions. In a complementary situation where the system is weakly relativistic but fully nonlinear (thus, far inside the horizon) we can employ the post-Newtonian approximation. We also show that in the large-scale structures the post-Newtonian effects are quite small. As a consequence, now we can rely on the Newtonian gravity in analyzing the evolution of nonlinear large-scale structures even near the horizon volume.

We consider a general relativistic zero-pressure irrotational cosmological medium perturbed to the third order. We assume a flat Friedmann background but include the cosmological constant. We ignore the rotational perturbation which decays in expanding phase. In our previous studies we discovered that, to the second-order perturbation, except for the gravitational wave contributions, the relativistic equations coincide exactly with the previously known Newtonian ones. Since the Newtonian second-order equations are fully nonlinear, any nonvanishing third and higher order terms in the relativistic analyses are supposed to be pure relativistic corrections. In this work we derive such correction terms appearing in the third order. Continuing our success in the second-order perturbations we take the comoving gauge. We discover that the third-order correction terms are of $\phi_v$-order higher than the second-order terms where $\phi_v$ is a gauge-invariant combination related to the three-space curvature perturbation in the comoving gauge; compared with the Newtonian potential we have $\delta \Phi \sim {3 \over 5} \phi_v$ to the linear order. Therefore, the pure general relativistic effects are of $varphi_v$-order higher than the Newtonian ones. The corrections terms are independent of the horizon scale and depend only on the linear order gravitational potential perturbation strength. From the temperature anisotropy of cosmic microwave background we have ${\delta T \over T} \sim {1 \over 3} \delta \Phi \sim {1 \over 5} \phi_v \sim 10^{-5}$. Therefore, our present result reinforces our previous important practical implication that near current era one can use the large-scale Newtonian numerical simulation more reliably even as the simulation scale approaches near the horizon.

It is well known that couplings occur among the scalar-, vector-, and tensor-type perturbations of Friedmann world model in the second perturbational order. Here, we prove that, except for the gravitational wave contribution, the relativistic zero-pressure irrotational fluid perturbed to second order in a flat Friedmann background coincides exactly with the Newtonian result. Since we include the cosmological constant, our results are relevant to currently favoured cosmology. As we prove that the Newtonian hydrodynamic equations are valid in all cosmological scales to the second order, our result has an important practical implication that one can now use the large-scale Newtonian numerical simulation more reliably even as the simulation scale approaches and even goes beyond the horizon. That is, our discovery shows that, in the zero-pressure case, except for the gravitational wave contribution, there are no relativistic correction terms even near and beyond the horizon to the second-order perturbation.

Until now, it has been common to use Newton's gravity to study the non-linear clustering properties of the large-scale structures. Without confirmation from Einstein's theory, however, it has been unclear whether we can rely on the analysis, for example, near the horizon scale. In this work we will provide a confirmation of using Newton's gravity in cosmology based on relativistic analysis of weakly non-linear situations to the third order in perturbations. We will show that, except for the gravitational wave contribution, the relativistic zero-pressure fluid equations perturbed to the second order in a flat Friedmann background coincide exactly with the Newtonian results. We will also present the pure relativistic correction terms appearing in the third order. The third-order correction terms show that these are the linear-order curvature perturbation strength higher than the second-order relativistic/Newtonian terms. Thus, the pure general relativistic corrections in the third order are independent of the horizon scale and are small in the large-scale due to the low-level temperature anisotropy of the cosmic microwave background radiation. Since we include the cosmological constant, our results are relevant to currently favoured cosmology. As we prove that the Newtonian hydrodynamic equations are valid in all cosmological scales to the second order, and that the third-order correction terms are small, our result has a practically important implication that one can now use the large-scale Newtonian numerical simulation more reliably as the simulation scale approaches and even goes beyond the horizon.

The dynamic world model and its linear perturbations were first studied in Einstein's gravity. In the system without pressure the relativistic equations coincide exactly with the later known ones in Newton's gravity. Here we prove that, except for the gravitational wave contribution, even to the second-order perturbations, equations for the relativistic irrotational zero-pressure fluid in a flat Friedmann background coincide exactly with the previously known Newtonian equations. Thus, to the second order, we correctly identify the relativistic density and velocity perturbation variables, and we expand the range of applicability of the Newtonian medium without pressure to all cosmological scales including the super-horizon scale. In the relativistic analyses, however, we do not have a relativistic variable which corresponds to the Newtonian potential to the second order. Mixed usage of different gauge conditions is useful to make such proofs and to examine the result with perspective. We also present the gravitational wave equation to the second order. Since our correspondence includes the cosmological constant, our results are relevant to currently favoured cosmology. Our result has an important practical implication that one can use the large-scale Newtonian numerical simulation more reliably even as the simulation scale approaches near horizon.

The post-Newtonian (PN) approximation, based on the assumptions of weak gravitational fields and slow motions, provides a way to estimate general relativistic effects in the fully nonlinear evolution stage of the large-scale cosmic structures. We extend Chandrasekhar's first order PN (1PN) hydrodynamics based on the Minkowski background to the Robertson-Walker background. We assume the presence of Friedmann's cosmological spacetime as a background. In the background we include the three-space curvature, the cosmological constant and general pressure. In the Newtonian order and 1PN order we include general pressure, stress, and flux. The Newtonian hydrodynamic equations appear naturally in the 0PN order. The spatial gauge degree of freedom is fixed in a unique manner and the basic equations are arranged without taking the temporal gauge condition. In this way we can conveniently try alternative temporal gauge conditions. We investigate a number of temporal gauge conditions under which all the remaining variables are equivalently gauge-invariant. Our aim is to present the fully nonlinear 1PN equations in a form suitable for implementation in conventional Newtonian hydrodynamic simulations with minimal extensions. The 1PN terms can be considered as relativistic corrections added to the well known Newtonian equations. The proper arrangement of the variables and equations in combination with suitable gauge conditions would allow the possible future 1PN cosmological simulations to become more tractable. Our equations and gauges are arranged for that purpose. We suggest ways of controlling the numerical accuracy. The typical 1PN order terms are about $10^{-6} \sim 10^{-4}$ times smaller than the Newtonian terms.

We present cosmological perturbation theory based on generalized gravity theories including string theory correction terms and a tachyonic complication. The classical evolution as well as the quantum generation processes in these variety of gravity theories are presented in unified forms. These apply both to the scalar- and tensor-type perturbations. Analyses are made based on the curvature variable in two different gauge conditions often used in the literature in Einstein's gravity; these are the curvature variables in the comoving (or uniform-field) gauge and the zero-shear gauge. Applications to generalized slow-roll inflations and its consequent power spectra are derived in unified forms which include wide range of inflationary scenarios based on Einstein's gravity and others.

We consider instability of the Friedmann world model to the second-order in perturbations. We present the perturbed set of equations up to the second-order in the Friedmann background world model with general spatial curvature and the cosmological constant. We consider systems with the completely general imperfect fluids, the minimally coupled scalar fields, the electro-magnetic field, and the generalized gravity theories. We also present the case of null geodesic equations, and the one based on the relativistic Boltzmann equation. In due stage a decomposition is made for the scalar-, vector- and tensor-type perturbations which couple each other to the second-order. Gauge issue is resolved to each order. The basic equations are presented without imposing any gauge condition, thus in a gauge-ready form so that we can use the full advantage of having the gauge freedom in analysing the problems. As an application we show that to the second-order in perturbation the relativistic pressureless ideal fluid of the scalar-type reproduces exactly the known Newtonian result. As another application we rederive the large-scale conserved quantities (of the pure scalar- and tensor-perturbations) to the second order, first shown by Salopek and Bond, now from the exact equations. Several other applications are made as well.

If the linear perturbation theory is valid through the bounce, the surviving fluctuations from the ekpyrotic scenario (cyclic one as well) should have very blue spectra with suppressed amplitude for the scalar-type structure. We derive the same (and consistent) result using the curvature perturbation in the uniform-field (comoving) gauge and in the zero-shear gauge. Previously, Khoury et al. interpreted results from the latter gauge condition incorrectly and claimed the scale-invariant spectrum, thus generating controversy in the literature. We also correct similar errors in the literature based on wrong mode identification and joining condition. No joining condition is needed for the derivation.

We present unified ways of handling the cosmological perturbations in a class of gravity theory covered by a general action in eq. (1). This gravity includes our previous generalized $f(\phi,R)$ gravity and the gravity theory motivated by the tachyonic condensation. We present general prescription to derive the power spectra generated from vacuum quantum fluctuations in the slow-roll inflation era. An application is made to a slow-roll inflation based on the tachyonic condensation with an exponential potential.

We consider evolutions of linear fluctuations as the background Friedmann world model goes from contracting to expanding phases through smooth and non-singular bouncing phases. As long as the gravity dominates over the pressure gradient in the perturbation equation the growing-mode in the expanding phase is characterized by a conserved amplitude, we call it a C-mode. In the spherical geometry with a pressureless medium, we show that there exists a special gauge-invariant combination \Phi which stays constant throughout the evolution from the big-bang to the big-crunch with the same value even after the bounce: it characterizes the coefficient of the C-mode. We show this result by using a bounce model where the pressure gradient term is negligible during the bounce; this requires additional presence of an exotic matter. In such a bounce, even in more general situations of the equation of states before and after the bounce, the C-mode in the expanding phase is affected only by the C-mode in the contracting phase, thus the growing mode in the contracting phase decays away as the world model enters expanding phase. In the case the background curvature has significant role during the bounce, the pressure gradient term becomes important and we cannot trace C-mode in the expanding phase to the one before the bounce. In such situations, perturbations in a fluid bounce model show exponential instability, whereas the ones in a scalar field bounce model show oscillatory behaviors.

We consider the evolution of perturbed cosmological spacetime with multiple fluids and fields in Einstein gravity. Equations are presented in gauge-ready forms, and are presented in various forms using the curvature (\Phi or \phi_\chi) and isocurvature (S_{(ij)} or \delta \phi_{(ij)}) perturbation variables in the general background with K and \Lambda. We clarify the conditions for conserved curvature and isocurvature perturbations in the large-scale limit. Evolutions of curvature perturbations in many different gauge conditions are analysed extensively. In the multi-field system we present a general solution to the linear order in slow-roll parameters.

We point out that a pseudo-Newtonian interpretations of the $(1/3)$ factor in the Sachs-Wolfe effect, which relates the fluctuations in temperature and potential, $(\delta T/ T) = (1/ 3) \delta \Phi$, is not supported by the General Relativistic analysis. Dividing the full gravitational effect into separate parts depends on the choice of time slicing (gauge) and there exist infinitely many different choices. More importantly, interpreting the parts as being due to the gravitational redshift and the time dilation is not justified in the rigorous relativistic perturbation theory. We suggest to regard the $(1/3)$ factor as the General Relativistic result which applies in a restricted situation of adiabatic perturbation in the $K = 0 = \Lambda$ model with the last scattering occuring in the matter dominated era. For an isocurvature initial condition the corresponding result, $(\delta T/ T) = 2 \delta \Phi$, has a different numerical coefficient.

We address the perturbation power spectrum generated in the recently proposed ekpyrotic scenario by Khoury et al. The issue has been raised recently by Lyth who used the conventional method based on a conserved variable in the large-scale limit, and derived different results from Khoury et al. The calculation is straightforward in the uniform-curvature gauge where the generated blue spectrum with suppressed amplitude survives as the final spectrum. Whereas, although the metric fluctuations become unimportant and a scale-invariant spectrum is generated in the zero-shear gauge the mode does not survive the bounce, thus with the same final result. Therefore, an exponential potential leads to a power-law expansion/contraction $a \propto |t|^p$, and the power $p$ dictates the final power spectra of both the scalar and tensor structures. If $p \ll 1$ as one realization of the ekpyrotic scenario suggests, the results are $n_S -1 \simeq 2 \simeq n_T$ and the amplitude of the scalar perturbation is suppressed relative to the one of the gravitational wave by a factor $\sqrt{p}/2$. Both results confirm Lyth's. An observation is made on the constraint on the dynamics of the seed generating stage from the requirement of scale-invariant spectrum.

We present cosmological perturbations of kinetic components based on relativistic Boltzmann equations in the context of generalized gravity theories. Our general theory considers an arbitrary number of scalar fields generally coupled with the gravity, an arbitrary number of mutually interacting hydrodynamic fluids, and components described by the relativistic Boltzmann equations like massive/massless collisionless particles and the photon with the accompanying polarizations. We also include direct interactions among fluids and fields. The background FLRW model includes the general spatial curvature and the cosmological constant. We consider three different types of perturbations, and all the scalar-type perturbation equations are arranged in a gauge-ready form so that one can implement easily the convenient gauge conditions depending on the situation. In the numerical calculation of the Boltzmann equations we have implemented four different gauge conditions in a gauge-ready manner where two of them are new. By comparing solutions solved separately in different gauge conditions we can naturally check the numerical accuracy.

The classical evolution and the quantum generation processes of the scalar- and tensor-type cosmological perturbations in the context of a broad class of generalized gravity theories are presented in unified forms. The exact forms of final spectra of the two types of structures generated during a generalized slow-roll inflation are derived. Results in generalized gravity are characterized by two additional parameters which are the coupling between gravity and field, and the nonminimal coupling in the kinetic part of the field. Our general results include widely studied gravity theories and inflation models as special cases, and show how the well known consistency relation and spectra in ordinary Einstein gravity inflation models are affected by the generalized nature of the gravity theories.

The scalar field with an exponential potential allows a scaling solution where the the density of the field follows the density of the dominating fluid. Such a scaling regime is often used as an important ingredient in many models of quintessence. We analyse evolution of perturbations while the background follows the scaling. As the results, the perturbed scalar field also scales with the perturbed fluid, and the perturbations accompany the adiabatic as well as the isocurvature mode between the fluid and the field.

In this paper we address two points concerning the Sachs-Wolfe effect: (i) the gauge independence of the observable temperature anisotropy, and (ii) a gauge-invariant expression of the effect considering the most general situation of hydrodynamic perturbations. The first result follows because the gauge transformation of the temperature fluctuation at the observation event only contributes to the isotropic temperature change which, in practice, is absorbed into the definition of the background temperature. Thus, we proceed without fixing the gauge condition, and express the Sachs-Wolfe effect using the gauge-invariant variables.

In a class of non-singular cosmologies derived from higher-order corrections to the low-energy bosonic string action, we derive evolution equations for the most general cosmological scalar, vector and tensor perturbations. In the large scale limit, the evolutions of both scalar and tensor perturbations are characterised by conserved quantities, the usual curvature perturbation in the uniform-field gauge and the tensor-type perturbed metric. The vector perturbation is not affected, being described by the conservation of the angular momentum of the fluid component in the absence of any additional dissipative process. For the scalar- and tensor-type perturbations, we show how, given a background evolution during kinetic driven inflation of the dilaton field, we can obtain the final power spectra generated from the vacuum quantum fluctuations of the metric and the dilaton field during the inflation.

We consider the large-scale cosmic structure generation by an inflation based on a pure f(R)-type gravity theory. Comparison with recent CMBR observations gives the following results: (1) The near Zel'dovich spectral conditions uniquely choose R^2-type gravity. (2) The R^2 gravity predicts specific nearly scale-invariant Zel'dovich spectra for both the scalar- and tensor-type perturbations. Thus, (3) the considered model survives current observational data. The COBE-DMR quadrupole data (4) give constraints on the coupling constant and the energy scale during inflation, and (5) require the gravitational wave contribution to be suppressed. (6) Therefore, future observations of (2) and (5) can provide strong tests of the inflation scenario based on R^2 gravity. Parallel analyses made in the conformally transformed Einstein frame give the observationally identical results.

We present the power spectra of the scalar- and tensor-type structures generated in an inflation model based on the nonminimally coupled scalar field with a self coupling. By comparing the contributions of these structures to the anisotropy of the cosmic microwave background radiation with the four year COBE DMR data we derive strong constraints on model parameters and the inflation model.

%auto-ignore This paper has been withdrawn by the author(s) since somebody do not agree with my idea.

This paper has been withdrawn by the author(s) since somebody do not agree with my idea.

%auto-ignore This paper has been withdrawn by the author(s) because somebody do not agree with my idea.

We consider the evolution of perturbed cosmological spacetime with multiple scalar fields in Einstein gravity. A complete set of scalar-type perturbation equations is presented in a gauge-ready form, and we derived the closed set of second-order differential equations in several useful forms. Conserved behaviors of the perturbed three-space curvature in the comoving gauge, $\phi_v$, under several conditions are clarified. Under the slow-roll conditions, the adiabatic and isocurvature modes decouple from each other, and in the large-scale limit we have (i) the adiabatic mode is generally conserved, (ii) for a couple of special potential the isocurvature modes decouple from each other and are described by conserved quantities, (iii) in the two field system, the isocurvature mode is described by a conserved quantity for the general potential.

A generic form of low-energy effective action of superstring theories with one-loop quantum correction is well known. Based on this action we derive the complete perturbation equations and general analytic solutions in the cosmological spacetime. Using the solutions we identify conserved quantities characterizing the perturbations: the amplitude of gravitational wave and the perturbed three-space curvature in the uniform-field gauge both in the large-scale limit, and the angular-momentum of rotational perturbation are conserved independently of changing gravity sector. Implications for calculating perturbation spectra generated in the inflation era based on the string action are presented.

The evolutions of linear structures in a spatially homogeneous and isotropic world model are characterized by some conserved quantities: the amplitude of gravitational wave is conserved in the super-horizon scale, the perturbed three-space curvature in the comoving gauge is conserved in the super-sound-horizon scale, and the angular momentum of rotational perturbation is generally conserved.

We derive power spectra of the scalar- and tensor-type structures generated in an inflation model based on a massive non-minimally coupled scalar field with the strong coupling assumption. We make analyses in both the original-frame and the conformally transformed Einstein-frame. We derive contributions of both structures to the anisotropy of the cosmic microwave background radiation, and compare the contributions with the four-year COBE-DMR data. Previous study showed that sufficient amount of inflation requires a small coupling parameter. In such a case the spectra become near Zeldovich spectra, and the gravitational wave contribution becomes negligible compared with the scalar-type contribution which is testable in future CMBR experiments.

We introduce a scalar-type perturbation variable $\Phi$ which is conserved in the large-scale limit considering general sign of three-space curvature ($K$), the cosmological constant ($\Lambda$), and time varying equation of state. In a pressureless medium $\Phi$ is {\it exactly conserved} in all scales.

Relativistic cosmological perturbation analyses can be made based on several different fundamental gauge conditions. In the pressureless limit the variables in certain gauge conditions show the correct Newtonian behaviors. Considering the general curvature ($K$) and the cosmological constant ($\Lambda$) in the background medium, the perturbed density in the comoving gauge, and the perturbed velocity and the perturbed potential in the zero-shear gauge show the same behavior as the Newtonian ones in general scales. In the first part, we elaborate these Newtonian correspondences. In the second part, using the identified gauge-invariant variables with correct Newtonian correspondences, we present the relativistic results with general pressures in the background and perturbation. We present the general super-sound-horizon scale solutions of the above mentioned variables valid for general $K$, $\Lambda$, and generally evolving equation of state. We show that, for vanishing $K$, the super-sound-horizon scale evolution is characterised by a conserved variable which is the perturbed three-space curvature in the comoving gauge. We also present equations for the multi-component hydrodynamic situation and for the rotation and gravitational wave.

We examine the effects of the axion coupling to $R\tilde{R}$ on the evolution of cosmic structures. It is shown that the evolutions of the scalar- and vector-type perturbations are not affected by this axion coupling. However the axion coupling causes an asymmetric evolution of the two polarization states of the tensor-type perturbation, which may lead to a sizable polarization asymmetry in the cosmological gravitational wave if inflation involves a period in which the axion coupling is important. The polarization asymmetry produced during inflation are conserved over the subsequent evolution as long as the scales remain in the large-scale limit, and thus this may lead to an observable trace in the cosmic microwave background radiation.

We present a set of equations describing the evolution of the scalar-type cosmological perturbation in a gravity with general quadratic order curvature coupling terms. Equations are presented in a gauge ready form, thus are ready to implement various temporal gauge conditions depending on the problems. The Ricci-curvature square term leads to a fourth-order differential equation for describing the spacetime fluctuations in a spatially homogeneous and isotropic cosmological background.

We derive the power spectra of the scalar- and tensor-type structures generated in a chaotic-type inflation model based on non-minimally coupled scalar field with a self interaction. By comparing contributions of both-structures to the anisotropy of the cosmic microwave background radiation with the four-year COBE-DMR data on quadrupole anisotropy we derive constraints on the ratio of self-coupling and non-minimal coupling constants, and the expansion rate in the inflation era, Eq. (12). The requirement of successful amount of inflation further constrains the relative amount of tensor-type contribution, Eq. (16).

We calculate the density and gravitational wave spectrums generated in a version of string cosmology termed pre-big bang scenario. The large scale structures are originated from quantum fluctuations of the metric and dilaton field during a pole-like inflation stage driven by a potential-less dilaton field realized in the low-energy effective action of string theory. The generated classical density field and the gravitational wave in the second horizon crossing epoch show tilted spectrums with n ~ 4 and n_T ~ 3, respectively. These differ from the observed spectrum of the large angular scale anisotropy of the cosmic microwave background radiation which supports scale invariant ones with n ~ 1 and n_T ~ 0. This suggests that the pre-big bang stage is not suitable for generating the present day observable large scale structures, and suggests the importance of investigating the quantum generation processes during stringy era with higher order quantum correction terms.

Relativistic cosmological perturbation analyses can be made based on several different fundamental gauge conditions. In the pressureless limit the variables in certain gauge conditions show the correct Newtonian behaviors. We consider the general curvature and the cosmological constant in the background medium. The perturbed density in the comoving gauge, and the perturbed velocity and the perturbed potential in the zero-shear gauge show the same behavior as the Newtonian ones in a general scale. Far inside horizon, except for the uniform-density gauge, density perturbations in all the fundamental gauge conditions show the correct Newtonian behavior. In this paper we elaborate these Newtonian correspondences. We also present the relativistic results considering general pressures in the background and perturbation.

A coherently oscillating scalar field, an axion as an example, is known to behave as a cold dark matter. The arguments were usually made in the Newtonian context. Ratra proved the case in relativistic context using the synchronous gauge. In this paper we present another proof based on a more suitable gauge choice, the uniform-curvature gauge, which fits the problem. By a proper time averaging the perturbed oscillating scalar field behaves as a cold dark matter on the relevant scales including the superhorizon scale.

In a class of generalized Einstein's gravity theories we derive the equations and general asymptotic solutions describing the evolution of the perturbed universe in unified forms. Our gravity theory considers general couplings between the scalar field and the scalar curvature in the Lagrangian, thus includes broad classes of generalized gravity theories resulting from recent attempts for the unification. We analyze both the scalar-type mode and the gravitational wave in analogous ways. For both modes the large scale evolutions are characterized by the same conserved quantities which are valid in the Einstein's gravity. This unified and simple treatment is possible due to our proper choice of the gauges, or equivalently gauge invariant combinations.

Fundamental formulae representing the density perturbation spectrums generated in inflationary scenarios are rigorously proved. Quantum fluctuations as initial conditions for structure generation in some inflationary era can be calculated in general scale. Thus, we could avoid the matching of the large and the small scale results at the horizon crossing epoch. This is possible because the perturbed scalar field equation in a special gauge has some remarkable nice properties. In the uniform-curvature gauge the contributions from the metric fluctuations effectively disappear in certain expansion stages, and also we can derive a simple general solution in the large scale limit.