Chromospheric activities prior to solar flares provide important clues to solar flare initiation, but are as yet poorly understood. We report a significant and rapid H$\alpha$ line broadening before the solar flare SOL2011-09-29T18:08, that was detected using the unprecedented high-resolution H$\alpha$ imaging spectroscopy with the Fast Imaging Solar Spectrograph (FISS) installed on the 1.6 m New Solar Telescope (NST) at Big Bear Solar Observatory.The strong H$\alpha$ broadening extends as a blue excursion up to -4.5 A and as a red excursion up to 2.0 A, which implies a mixture of velocities in the range of -130 km/s to 38 km/s by applying the cloud model, comparable to the highest chromospheric motions reported before. The H$\alpha$ blueshifted broadening lasts for about 4 minutes, and is temporally and spatially correlated with the start of a rising filament, which is later associated with the main phase of the flare as detected by the Atmosphere Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The potential importance of this H$\alpha$ blueshifted broadening as a preflare chromospheric activity is briefly discussed within the context of the two-step eruption model.

Umbral flashes (UFs) are emissions in the core of chromospheric lines caused by upward propagating waves steepening into shocks. UFs are followed by an expanding blue shifted umbral wave (UW) and red-shifted plasma returning to the initial state. Here we use 5 s cadence images acquired at $\pm$0.04~nm off the Halpha line center by the Visible Imaging Spectrometer (VIS) installed on the Goode Solar Telescope (GST) to detect the origin of UFs and UWs in a sunspot with a uniform umbra free of LBs and clusters of umbral dots. The data showed that UFs do not randomly originate over the umbra. Instead, they appear to be repeatedly triggered at locations with the lowest umbral intensity and the most powerful oscillations of Halpha-0.04 nm intensity. GST magnetic field measurements using Near Infra-Red Imaging Spectropolarimeter (NIRIS) also showed that the dominant location of prevalent UF origin is co-spatial associated with the strongest fields in the umbra. Interface Region Imaging Spectrograph 149.0 nm images showed that no bright UV loops were anchored in the umbra in general and near the UF patches in particular suggesting that UFs and UWs alone can not be responsible for the origin of warm coronal loops. We thus conclude that the existence of locations with prevalent origin of UFs confirms the idea that they may be driven by a sub-surface source located near the axis of a flux rope, while the presence of several UFs trigger centers may indicate the complex structure of a sunspot umbra.

Characterizing the physical properties of the stochastic gravitational waves background (SGWB) is a key step towards identifying the nature of its possible origin. We focus our analysis on SGWB anisotropies. The existence of a non-trivial primordial scalar-tensor-tensor (STT) correlation in the squeezed configuration may be inferred from the effect that a long wavelength scalar mode has on the gravitational wave power spectrum: an anisotropic contribution. Crucially, such contribution is correlated with temperature anisotropies in the cosmic microwave background (CMB). We show that, for inflationary models that generate suitably large STT non-Gaussianities, cross-correlating the CMB with the stochastic background of gravitational waves is a very effective probe of early universe physics. The resulting signal can be a smoking-gun for primordial SGWB anisotropies.

We study the robustness of single-field inflation against inhomogeneities. We derive a simple analytic criterion on the shape of the potential for successful inflation in the presence of inhomogeneities, and demonstrate its validity using full 3+1 dimensional numerical relativity simulations on several classes of popular models of single-field inflation. We find that models with convex potentials are more robust to inhomogeneities than those with concave potentials, and that concave potentials that vary on super-Planckian scales are significantly more robust than those that vary on sub-Planckian scales.

We construct, for the first time, the time-domain gravitational wave strain waveform from the collapse of a strongly gravitating Abelian Higgs cosmic string loop in full general relativity. We show that the strain exhibits a large memory effect during merger, ending with a burst and the characteristic ringdown as a black hole is formed. Furthermore, we investigate the waveform and energy emitted as a function of string width, loop radius and string tension $G\mu$. We find that the mass normalized gravitational wave energy displays a strong dependence on the inverse of the string tension $E_{\mathrm{GW}}/M_0\propto 1/G\mu$, with $E_{\mathrm{GW}}/M_0 \sim {\cal O}(1)\%$ at the percent level, for the regime where $G\mu\gtrsim10^{-3}$. Conversely, we show that the efficiency is only weakly dependent on the initial string width and initial loop radii. Using these results, we argue that gravitational wave production is dominated by kinematical instead of geometrical considerations.

We investigate the physics of black hole formation from the head-on collisions of boosted equal mass Oscillatons (OS) in full numerical relativity, for both the cases where the OS have equal phases or are maximally off-phase (anti-phase). While unboosted OS collisions will form a BH as long as their initial compactness $\mathcal{C}\equiv GM/R$ is above a numerically determined critical value $\mathcal{C}>0.035$, we find that imparting a small initial boost counter-intuitively \emph{prevents} the formation of black holes even if $\mathcal{C}> 0.035$. If the boost is further increased, at very high boosts $\gamma>1/12\mathcal{C}$, BH formation occurs as predicted by the hoop conjecture. These two limits combine to form a "stability band" where collisions result in either the OS "passing through" (equal phase) or "bouncing back" (anti-phase), with a critical point occurring around ${\cal C}\approx 0.07$. We argue that the existence of this stability band can be explained by the competition between the free fall and the interaction timescales of the collision.

The umbral oscillations of velocity are commonly observed in the chromosphere of a sunspot. Their sources are considered to be either the external p-mode driving or the internal excitation by magnetoconvection. Even though the possibility of the p-mode driving has been often considered, the internal excitation has been rarely investigated. We report the identification of the oscillation patterns that may be closely related to the events of internal excitation from the observations of velocity oscillations in the temperature minimum region of two sunspots. The velocities were determined from the spectral data of the Fe I 5435 A line, a magnetically insensitive line, taken with the Fast Imaging Solar Spectrograph of the 1.6 m Goode Solar Telescope at the Big Bear Solar Observatory. As a result, we identified 4 oscillation patterns of 2.0 * 10^3 km coherent size that were clearly identified for about 7.9 minutes with oscillation amplitude of 9.3 * 10^-2 km/s. The power of the oscillations in these centers was concentrated in the 3 minute band. All the oscillation centers were located above the umbral dots undergoing noticeable morphological and dynamical changes that may be regarded as an observable signature of small-scale magnetoconvection inside the umbrae. Our results support the notion that magnetoconvection associated with umbral dots inside sunspots can drive the 3 minute umbral oscillations.

We present the first fully general relativistic dynamical simulations of Abelian Higgs cosmic strings using 3+1D numerical relativity. Focusing on cosmic string loops, we show that they collapse due to their tension and can either (i) unwind and disperse or (ii) form a black hole, depending on their tension $G\mu$ and initial radius. We show that these results can be predicted using an approximate formula derived using the hoop conjecture, and argue that it is independent of field interactions. We extract the gravitational waveform produced in the black hole formation case and show that it is dominated by the $l=2$ and $m=0$ mode. We also compute the total gravitational wave energy emitted during such a collapse, being $0.5\pm 0.2~ \%$ of the initial total cosmic string loop mass, for a string tension of $G\mu=1.6\times 10^{-2}$ and radius $R=100~M_{pl}^{-1}$. We use our results to put a bound on the production rate of planar cosmic strings loops as $N \lesssim 10^{-2}~\mathrm{Gpc}^{-3}~\mathrm{yr}^{-1}$.

We numerically investigate the gravitational waves generated by the head-on collision of equal-mass, self-gravitating, real scalar field solitons (oscillatons) as a function of their compactness $\mathcal{C}$. We show that there exist three different possible outcomes for such collisions: (1) an excited stable oscillaton for low $\mathcal{C}$, (2) a merger and formation of a black-hole for intermediate $\mathcal{C}$, and (3) a pre-merger collapse of both oscillatons into individual black-holes for large $\mathcal{C}$. For (1), the excited, aspherical oscillaton continues to emit gravitational waves. For (2), the total energy in gravitational waves emitted increases with compactness, and possesses a maximum which is greater than that from the merger of a pair of equivalent mass black-holes. The initial amplitudes of the quasi-normal modes in the post-merger ring-down in this case are larger than that of corresponding mass black-holes -- potentially a key observable to distinguish black-hole mergers with their scalar mimics. For (3), the gravitational wave output is indistinguishable from a similar mass, black-hole--black-hole merger.

We describe the properties of evanescent gravitational waves (EGWs)--wave solutions of Einstein equations which decay exponentially in some direction while propagating in another. Evanescent waves are well known in acoustics and optics and have recently received much attention due to their extraordinary properties such as their transverse spin and spin-momentum locking. We show that EGWs possess similarly remarkable properties, carrying transverse spin angular momenta and driving freely falling test masses along in elliptical trajectories. Hence, test masses on a plane transverse to the direction of propagation exhibit vector and scalar-like deformation, complicating efforts to use detection of such modes as evidence for modified gravity. We demonstrate that EGWs are present and dominant in the vicinity of sub-wavelength sources such as orbiting binaries.

The next generation of instruments designed to measure the polarization of the cosmic microwave background (CMB) will provide a historic opportunity to open the gravitational wave window to the primordial Universe. Through high sensitivity searches for primordial gravitational waves, and tighter limits on the energy released in processes like phase transitions, the CMB polarization data of the next decade has the potential to transform our understanding of the laws of physics underlying the formation of the Universe.

Axions and axion-like particles are compelling candidates for the missing dark matter of the universe. As they undergo gravitational collapse, they can form compact objects such as axion stars or even black holes. In this paper, we study the formation and distribution of such objects. First, we simulate the formation of compact axion stars using numerical relativity with aspherical initial conditions that could represent the final stages of axion dark matter structure formation. We show that the final states of such collapse closely follow the known relationship of initial mass and axion decay constant $f_a$. Second, we demonstrate with a toy model how this information can be used to scan a model density field to predict the number densities and masses of such compact objects. In addition to being detectable by the LIGO/VIRGO gravitational wave interferometer network for axion mass of $10^{-9} < m_a < 10^{-11}$ eV, we show using peak statistics that for $f_a < 0.2M_{pl}$, there exists a "mass gap" between the masses of axion stars and black holes formed from collapse.

Extending our previous work on the robustness of inflation to perturbations in the scalar field, we investigate the effects of perturbations in the transverse traceless part of the extrinsic curvature on the evolution of an inhomogeneous inflaton field. Focusing on small field models, we show that these additional metric inhomogeneities initially reduce the total number of e-folds as the amplitude increases, but that the reduction saturates and even reverses above a certain amplitude. We present an argument that this is due to the presence of a large initial Hubble friction when metric perturbations are large.

We consider the effects of inhomogeneous initial conditions in both the scalar field profile and the extrinsic curvature on different inflationary models. In particular, we compare the robustness of small field inflation to that of large field inflation, using numerical simulations with Einstein gravity in 3+1 dimensions. We find that small field inflation can fail in the presence of subdominant gradient energies, suggesting that it is much less robust to inhomogeneities than large field inflation, which withstands dominant gradient energies. However, we also show that small field inflation can be successful even if some regions of spacetime start out in the region of the potential that does not support inflation. In the large field case, we confirm previous results that inflation is robust if the inflaton occupies the inflationary part of the potential. Furthermore, we show that increasing initial scalar gradients will not form sufficiently massive inflation-ending black holes if the initial hypersurface is approximately flat. Finally, we consider the large field case with a varying extrinsic curvature K, such that some regions are initially collapsing. We find that this may again lead to local black holes, but overall the spacetime remains inflationary if the spacetime is open, which confirms previous theoretical studies.

We report a different type of three-minute chromospheric oscillations above a sunspot in association with a small-scale impulsive event in a light bridge. During our observations, we found a transient brightening in the light bridge. The brightening was composed of elementary bursts that may be a manifestation of fast repetitive magnetic reconnections in the light bridge. Interestingly, the oscillations in the nearby sunspot umbra were impulsively excited when the intensity of the brightening reached its peak. The initial period of the oscillations was about 2.3 minutes and then gradually increased to 3.0 minutes with time. In addition, we found that the amplitude of the excited oscillations was twice the amplitude of oscillations before the brightening. Based on our results, we propose that magnetic reconnection occurring in a light bridge can excite of oscillations in the nearby sunspot umbra.

Solar flare ribbons provide an important clue to the magnetic reconnection process and associated magnetic field topology in the solar corona. We detected a large-scale secondary flare ribbon of a circular shape that developed in association with two successive M-class flares and one CME. The ribbon revealed interesting properties such as 1) a quasi-circular shape and enclosing the central active region; 2) the size as large as 500\arcsec\, by 650\arcsec\,, 3) successive brightenings in the clockwise direction at a speed of \kms{160} starting from the nearest position to the flaring sunspots, 4) radial contraction and expansion in the northern and the southern part, respectively at speeds of $\leq$ \kms{10}. Using multi-wavelength data from \textit{SDO}, \textit{RHESSI}, XRT, and Nobeyama, along with magnetic field extrapolations, we found that: 1) the secondary ribbon location is consistent with the field line footpoints of a fan-shaped magnetic structure that connects the flaring region and the ambient decaying field; 2) the second M2.6 flare occurred when the expanding coronal loops driven by the first M2.0 flare encountered the background decayed field. 3) Immediately after the second flare, the secondary ribbon developed along with dimming regions. Based on our findings, we suggest that interaction between the expanding sigmoid field and the overlying fan-shaped field triggered the secondary reconnection that resulted in the field opening and formation of the quasi-circular secondary ribbon. We thus conclude that interaction between the active region and the ambient large-scale fields should be taken into account to fully understand the entire eruption process.

The classical equations of motion for an axion with potential $V(\phi)=m_a^2f_a^2 [1-\cos (\phi/f_a)]$ possess quasi-stable, localized, oscillating solutions, which we refer to as "axion stars". We study, for the first time, collapse of axion stars numerically using the full non-linear Einstein equations of general relativity and the full non-perturbative cosine potential. We map regions on an "axion star stability diagram", parameterized by the initial ADM mass, $M_{\rm ADM}$, and axion decay constant, $f_a$. We identify three regions of the parameter space: i) long-lived oscillating axion star solutions, with a base frequency, $m_a$, modulated by self-interactions, ii) collapse to a BH and iii) complete dispersal due to gravitational cooling and interactions. We locate the boundaries of these three regions and an approximate "triple point" $(M_{\rm TP},f_{\rm TP})\sim (2.4 M_{pl}^2/m_a,0.3 M_{pl})$. For $f_a$ below the triple point BH formation proceeds during winding (in the complex $U(1)$ picture) of the axion field near the dispersal phase. This could prevent astrophysical BH formation from axion stars with $f_a\ll M_{pl}$. For larger $f_a\gtrsim f_{\rm TP}$, BH formation occurs through the stable branch and we estimate the mass ratio of the BH to the stable state at the phase boundary to be $\mathcal{O}(1)$ within numerical uncertainty. We discuss the observational relevance of our findings for axion stars as BH seeds, which are supermassive in the case of ultralight axions. For the QCD axion, the typical BH mass formed from axion star collapse is $M_{\rm BH}\sim 3.4 (f_a/0.6 M_{pl})^{1.2} M_\odot$.

It is well-known that light bridges inside a sunspot produce small-scale plasma ejections and transient brightenings in the chromosphere, but the nature and origin of such phenomena are still unclear. Utilizing the high-spatial and high temporal resolution spectral data taken with the Fast Imaging Solar Spectrograph and the TiO 7057 A broadband filter images installed at the 1.6 meter New Solar Telescope of Big Bear Solar Observatory, we report arcsecond-scale chromospheric plasma ejections (1.7 arc) inside a light bridge. Interestingly, the ejections are found to be a manifestation of upwardly propagating shock waves as evidenced by the sawtooth patterns seen in the temporal-spectral plots of the Ca II 8542 A and H-alpha intensities. We also found a fine-scale photospheric pattern (1 arc) diverging with a speed of about 2 km/s two minutes before the plasma ejections, which seems to be a manifestation of magnetic flux emergence. As a response to the plasma ejections, the corona displayed small-scale transient brightenings. Based on our findings, we suggest that the shock waves can be excited by the local disturbance caused by magnetic reconnection between the emerging flux inside the light bridge and the adjacent umbral magnetic field. The disturbance generates slow-mode waves, which soon develop into shock waves, and manifest themselves as the arcsecond-scale plasma ejections. It also appears that the dissipation of mechanical energy in the shock waves can heat the local corona.

We report three-minute oscillations in the solar chromosphere driven by a strong downflow event in a sunspot. We used the Fast Imaging Solar Spectrograph of the 1.6 m New Solar Telescope and the Interface Region Imaging Spectrograph (IRIS). The strong downflow event is identified in the chromospheric and transition region lines above the sunspot umbra. After the event, oscillations occur at the same region. The amplitude of the Doppler velocity oscillations is 2 km/s, and gradually decreases with time. In addition, the period of the oscillations gradually increases from 2.7 minutes to 3.3 minutes. In the IRIS 1330 slit-jaw images, we identify a transient brightening near the footpoint of the downflow detected in the Ha+0.5A image. The characteristics of the downflowing material are consistent with those of sunspot plumes. Based on our findings, we suggest that the gravitationally stratified atmosphere came to oscillate with three minute period in response to the impulsive downflow event as was theoretically investigated by Chae & Goode(2015).

The coronagraph is an instrument enables the investigation of faint features in the vicinity of the Sun, particularly coronal mass ejections (CMEs). So far coronagraphic observations have been mainly used to determine the geometric and kinematic parameters of CMEs. Here, we introduce a new method for the determination of CME temperature using a two filter (4025 A and 3934 A) coronagraph system. The thermal motion of free electrons in CMEs broadens the absorption lines in the optical spectra that are produced by the Thomson scattering of visible light originating in the photosphere, which affects the intensity ratio at two different wavelengths. Thus the CME temperature can be inferred from the intensity ratio measured by the two filter coronagraph system. We demonstrate the method by invoking the graduated cylindrical shell (GCS) model for the 3 dimensional CME density distribution and discuss its significance.

We studied temporal changes of morphological and magnetic properties of a succession of four confined flares followed by an eruptive flare using the high-resolution New Solar Telescope (NST) operating at the Big Bear Solar Observatory (BBSO), Helioseismic and Magnetic Imager (HMI) magnetograms and Atmospheric Image Assembly (AIA) EUV images provided by Solar Dynamics Observatory (SDO). From the NST/Halpha and the SDO/AIA~304 A observations we found that each flare developed a jet structure that evolved in a manner similar to evolution of the blowout jet : 1) an inverted-Y shape jet appeared and drifted away from its initial position; 2) jets formed a curtain-like structure that consisted of many fine threads accompanied with subsequent brightenings near the footpoints of the fine threads; and finally 3) the jet showed a twisted structure visible near the flare maximum. Analysis of the HMI data showed that both the negative magnetic flux and the magnetic helicity have been gradually increasing in the positive polarity region indicating the continuous injection of magnetic twist before and during the series of flares. Based on these results, we suggest that the continuous emergence of twisted magnetic flux played an important role in producing a successive flares and developing a series of blowout jets.

Theories of gravitation without Lorentz invariance are candidates of low-energy descriptions of quantum gravity. In this review we will describe the phenomenological consequences of the candidates associated to the existence of a preferred time direction

Two of the most commonly used tools to constrain the primordial non-Gaussianity are the bispectrum and the Minkowski functionals of CMB temperature anisotropies. These two measures of non-Gaussianity in principle provide distinct (though correlated) information, but in the past constraints from them have only been loosely compared, and not statistically combined. In this work we evaluate, for the first time, the covariance matrix between the local non-Gaussianity coefficient fnl estimated through the bispectrum and Minkowski functionals. We find that the estimators are positively correlated, with corerlation coefficient r ~ 0.3. Using the WMAP7 data to combine the two measures and accounting for the point-source systematics, we find the combined constraint fnl=37+/-28, which has a ~20% smaller error than either of the individual constraints.

The Planck mission has conclusively detected lensing of the Cosmic Microwave Background (CMB) radiation from foreground sources to an overall significance of greater than $25\sigma$. The high precision of this measurement motivates the development of a more complete formulation of the calculation of this effect. While most effects on the CMB anisotropies are widely studied through direct solutions of the Boltzmann equation, the non-linear effect of CMB lensing is formulated through the solutions of the geodesic equation. In this paper, we present a new formalism to the calculation of the lensing effect by \emph{directly solving the Boltzmann equation}, as we did in the calculation of the CMB anisotropies at recombination. In particular, we developed a diagrammatic approach to efficiently keep track of all the interaction terms and calculate all possible non-trivial correlations to arbitrary high orders. Using this formalism, we explicitly articulate the approximations required to recover the usual remapping approach used in current studies of the weak lensing. In addition, we point out additional unexplored corrections that are manifest in our formalism to which experiments may be sensitive. As an example, we calculate the correction to the CMB temperature power spectrum for the \emph{lens-lens} coupling effects which are neglected in standard calculations. We find that the correction is $\lesssim 0.1\%$ of the CMB temperature power spectrum for $\ell$ up to 3000 and thus is comparable to the cosmic variance.

We construct a perturbative framework for understanding the collision of solitons (more precisely, solitary waves) in relativistic scalar field theories. Our perturbative framework is based on the suppression of the space-time interaction area proportional to $1/(\gamma v)$, where $v$ is the relative velocity of an incoming solitary wave and $\gamma = 1/\sqrt{1-v^2} \gg 1$. We calculate the leading order results for collisions of (1+1) dimensional kinks in periodic potentials, and provide explicit, closed form expressions for the phase shift and the velocity change after the collisions. We find excellent agreement between our results and detailed numerical simulations. Crucially, our perturbation series is controlled by a kinematic parameter, and hence not restricted to small deviations around integrable cases such as the Sine-Gordon model.

We derive a closed-form expression for the phase shift experienced by 1+1 dimensional kinks colliding at ultra-relativistic velocities (gamma v >> 1), valid for arbitrary periodic potentials. Our closed-form expression is the leading order result of a more general scattering theory of solitary waves described in a companion paper [Amin, Lim and Yang, arXiv:1308.0606]. This theory relies on a small kinematic parameter 1/(gamma v) << 1 rather than a small parameter in the Lagrangian. Our analytic results can be directly extracted from the Lagrangian without solving the equation of motion. Based on our closed-form expression, we prove that kink-kink and kink-antikink collisions have identical phase shifts at leading order.

We prove that the linear instability in a non-degenerate higher derivative theory, the Ostrogradski instability, can only be removed by the addition of constraints if the original theory's phase space is reduced.

We study cosmological perturbations by solving the governing Boltzmann and Einstein Field equations up to second order, and calculate the corresponding CMB bispectrum during recombination. We include all the second-order Liouville and collision terms, truncating the multipole hierarchy at $l=10$, consistently including all $m\neq 0$ terms when calculating the bispectrum in the flat-sky limit. At this stage, we focus on contributions at recombination and we neglect 2nd-order vector and tensor terms, lensing effects, and late time non-linear ISW. We find that the signal-to-noise for the bispectrum is 0.69 for $l_{max}=2000$, yielding an overall signal $F_{\rm NL} = 3.19$ (normalised relative to the local model). We find that the effective $f_{NL}$'s of the equilateral and local type are 5.11 and 0.88 respectively. This recombination signal will have to be taken into account in a quantitative analysis of the Planck data.

The formation and the temporal evolution of a bipolar moving magnetic feature (MMF) was studied with high spatial and temporal resolution. The photometric properties were observed with the New Solar Telescope at Big Bear Solar Observatory using a broadband TiO filter (705.7 nm), while the magnetic field was analyzed using the spectropolarimetric data obtained by Hinode. For the first time, we observed a bipolar MMF simultaneously in intensity images and magnetic field data, and studied the details of its structure. The vector magnetic field and the Doppler velocity of the MMF were also studied. A bipolar MMF having its positive polarity closer to the negative penumbra formed being accompanied by a bright, filamentary structure in the TiO data connecting the MMF and a dark penumbral filament. A fast downflow (<2km/s) was detected at the positive polarity. The vector magnetic field obtained from the full Stokes inversion revealed that a bipolar MMF has a U-shaped magnetic field configuration. Our observations provide a clear intensity counterpart of the observed MMF in the photosphere, and strong evidence of the connection between the MMF and the penumbral filament as a serpentine field.

We investigate the effect of modifed gravity with screening mechanisms, such as the chameleon or symmetron models, upon the structure of main sequence stars. We find that unscreened stars can be significantly more luminous and ephemeral than their screened doppelgangers. By embedding these stars into dwarf galaxies, which can be unscreened for values of the parameters not yet ruled out observationally, we show that the cumulative effect of their increased luminosity can enhance the total galactic luminosity. We estimate this enhancement and find that it can be considerable given model parameters that are still under experimental scrutiny. By looking for systematic offsets between screened dwarf galaxies in clusters and unscreened galaxies in voids, these effects could form the basis of an independent observational test that can potentially lower the current experimental bounds on the model independent parameters of these theories by and order of magnitude or more.

We provide analytic solutions for the power spectrum and bispectrum of curvature fluctuations produced by a step feature in the inflaton potential, valid in the limit that the step is short and sharp. In this limit, the bispectrum is strongly scale dependent and its effective non-linearity attains a large oscillatory amplitude. The perturbations to the curvature power spectrum, on the other hand, remain a small component on top of the usual spectrum of fluctuations generated by slow roll. We utilize our analytic solutions to assess the observability of the predicted non-Gaussian signatures and show that, if present, only very sharp steps on scales larger than ~ 2 Gpc are likely to be able to be detected by Planck. Such features are not only consistent with WMAP7 data, but can also improve its likelihood by 2 Delta ln L ~ 12 for two extra parameters, the step location and height. If this improvement were due to a slow roll violating step as considered here, a bispectrum or corresponding polarization power spectrum detection would provide definitive checks as to its primordial origin.

In this paper, we show that Minkowski Functionals (MFs) of weak gravitational lensing (WL) convergence maps contain significant non-Gaussian, cosmology-dependent information. To do this, we use a large suite of cosmological ray-tracing N-body simulations to create mock WL convergence maps, and study the cosmological information content of MFs derived from these maps. Our suite consists of 80 independent 512^3 N-body runs, covering seven different cosmologies, varying three cosmological parameters Omega_m, w, and sigma_8 one at a time, around a fiducial LambdaCDM model. In each cosmology, we use ray-tracing to create a thousand pseudo-independent 12 deg^2 convergence maps, and use these in a Monte Carlo procedure to estimate the joint confidence contours on the above three parameters. We include redshift tomography at three different source redshifts z_s=1, 1.5, 2, explore five different smoothing scales theta_G=1, 2, 3, 5, 10 arcmin, and explicitly compare and combine the MFs with the WL power spectrum. We find that the MFs capture a substantial amount of information from non-Gaussian features of convergence maps, i.e. beyond the power spectrum. The MFs are particularly well suited to break degeneracies and to constrain the dark energy equation of state parameter w (by a factor of ~ three better than from the power spectrum alone). The non-Gaussian information derives partly from the one-point function of the convergence (through V_0, the "area" MF), and partly through non-linear spatial information (through combining different smoothing scales for V_0, and through V_1 and V_2, the boundary length and genus MFs, respectively). In contrast to the power spectrum, the best constraints from the MFs are obtained only when multiple smoothing scales are combined.

Recent studies have shown that the number counts of convergence peaks N(kappa) in weak lensing (WL) maps, expected from large forthcoming surveys, can be a useful probe of cosmology. We follow up on this finding, and use a suite of WL convergence maps, obtained from ray-tracing N-body simulations, to study (i) the physical origin of WL peaks with different heights, and (ii) whether the peaks contain information beyond the convergence power spectrum P_ell. In agreement with earlier work, we find that high peaks (with amplitudes >~ 3.5 sigma, where sigma is the r.m.s. of the convergence kappa) are typically dominated by a single massive halo. In contrast, medium-height peaks (~0.5-1.5 sigma) cannot be attributed to a single collapsed dark matter halo, and are instead created by the projection of multiple (typically, 4-8) halos along the line of sight, and by random galaxy shape noise. Nevertheless, these peaks dominate the sensitivity to the cosmological parameters w, sigma_8, and Omega_m. We find that the peak height distribution and its dependence on cosmology differ significantly from predictions in a Gaussian random field. We directly compute the marginalized errors on w, sigma_8, and Omega_m from the N(kappa) + P_ell combination, including redshift tomography with source galaxies at z_s=1 and z_s=2. We find that the N(kappa) + P_ell combination has approximately twice the cosmological sensitivity compared to P_ell alone. These results demonstrate that N(kappa) contains non-Gaussian information complementary to the power spectrum.

We studied flux emergence events of sub-granular scale in a solar active region. New Solar Telescope (NST) of Big Bear Solar Observatory made it possible to clearly observe the photospheric signature of flux emergence with very high spatial (0".11 at 7057{\AA}) and temporal (15 s) resolution. From TiO observations with the pixel scale of 0".0375, we found several elongated granule-like features (GLFs) stretching from the penumbral filaments of a sunspot at a relatively high speed of over 4 km s-1. After a slender arched darkening appeared at a tip of a penumbral filament, a bright point (BP) developed and quickly moved away from the filament forming and stretching a GLF. The size of a GLF was approximately 0.5" wide and 3" long. The moving BP encountered nearby structures after several minutes of stretching, and a well-defined elongated shape of a GLF faded away. Magnetograms from SDO/HMI and NST/IRIM revealed that those GLFs are photospheric indicators of small-scale flux emergence, and their disappearance is related to magnetic cancellation. From two well-observed events, we describe detailed development of the sub-structures of GLFs, and different cancellation processes that each of the two GLFs underwent.

In this paper, we investigate the utility of Minkowski Functionals as a probe of cold/hot disk-like structures in the CMB. In order to construct an accurate estimator, we resolve a long-standing issue with the use of Minkowski Functionals as probes of the CMB sky -- namely that of systematic differences ("residuals") when numerical and analytical MF are compared. We show that such residuals are in fact by-products of binning, and not caused by pixelation or masking as originally thought. We then derive a map-independent estimator that encodes the effects of binning, applicable to beyond our present work. Using this residual-free estimator, we show that small disk-like effects (as claimed by Vielva et al.) can be detected only when a large sample of such maps are averaged over. In other words, our estimator is noise-dominated for small disk sizes at WMAP resolution. To confirm our suspicion, we apply our estimator to the WMAP7 data to obtain a null result.

A magnetic channel - a series of polarity reversals separating elongated flux threads with opposite polarities - may be a manifestation of a highly non-potential magnetic configuration in active regions. To understand its formation we have carried out a detailed analysis of the magnetic channel in AR 10930 using data taken by the Solar Optical Telescope/Hinode. As a result, we found upflows (-0.5 to -1.0 km/s) and downflows (+1.5 to +2.0 km/s) inside and at both tips of the thread respectively, and a pair of strong vertical currents of opposite polarity along the channel. Moreover, our analysis of the nonlinear force-free fields constructed from the photospheric magnetic field indicates that the current density in the lower corona may have gradually increased as a result of the continuous emergence of the highly sheared flux along the channel. With these results, we suggest that the magnetic channel originates from the emergence of the twisted flux tube that has formed below the surface before the emergence.

It has recently been shown in high resolution numerical simulations that relativistic collisions of bubbles in the context of a multi-vacua potential may lead to the creation of bubbles in a new vacuum. In this paper, we show that scalar fields with only potential interactions behave like free fields during high-speed collisions; the kick received by them in a collision can be deduced simply by a linear superposition of the bubble wall profiles. This process is equivalent to the scattering of solitons in 1+1 dimensions. We deduce an expression for the field excursion (shortly after a collision), which is related simply to the field difference between the parent and bubble vacua, i.e. contrary to expectations, the excursion cannot be made arbitrarily large by raising the collision energy. There is however a minimum energy threshold for this excursion to be realized. We verify these predictions using a number of 3+1 and 1+1 numerical simulations. A rich phenomenology follows from these collision induced excursions - they provide a new mechanism for scanning the landscape, they might end/begin inflation, and they might constitute our very own big bang, leaving behind a potentially observable anisotropy.

We introduce a novel class of field theories where energy always flows along timelike geodesics, mimicking in that respect dust, yet which possess non-zero pressure. This theory comprises two scalar fields, one of which is a Lagrange multiplier enforcing a constraint between the other's field value and derivative. We show that this system possesses no wave-like modes but retains a single dynamical degree of freedom. Thus, the sound speed is always identically zero on all backgrounds. In particular, cosmological perturbations reproduce the standard behaviour for hydrodynamics with vanishing sound speed. Using all these properties we propose a model unifying Dark Matter and Dark Energy in a single degree of freedom. In a certain limit this model exactly reproduces the evolution history of Lambda-CDM, while deviations away from the standard expansion history produce a potentially measurable difference in the evolution of structure.

We consider the 3-pt function (i.e. the bispectrum or non-Gaussianity) for stochastic backgrounds of gravitational waves. We estimate the amplitude of this signal for the primordial inflationary background, gravitational waves generated during preheating, and for gravitational waves produced by self-ordering scalar fields following a global phase transition. To assess detectability, we describe how to extract the 3-pt signal from an idealized interferometric experiment and compute the signal to noise ratio as a function of integration time. The 3-pt signal for the stochastic gravitational wave background generated by inflation is unsurprisingly tiny. For gravitational radiation generated by purely causal, classical mechanisms we find that, no matter how non-linear the process is, the 3-pt correlations produced vanish in direct detection experiments. On the other hand, we show that in scenarios where the B-mode of the CMB is sourced by gravitational waves generated by a global phase transition, a strong 3-pt signal among the polarization modes could also be produced. This may provide another method of distinguishing inflationary B-modes. To carry out this computation, we have developed a diagrammatic approach to the calculation of stochastic gravitational waves sourced by scalar fluids, which has applications beyond the present scenario.

Given a scalar field with metastable minima, bubbles nucleate quantum mechanically. When bubbles collide, energy stored in the bubble walls is converted into kinetic energy of the field. This kinetic energy can facilitate the classical nucleation of new bubbles in minima that lie below those of the "parent" bubbles. This process is efficient and classical, and changes the dynamics and statistics of bubble formation in models with multiple vacua, relative to that derived from quantum tunneling.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We explore the consequences of the existence of a very large number of light scalar degrees of freedom in the early universe. We distinguish between participator and spectator fields. The former have a small mass, and can contribute to the inflationary dynamics; the latter are either strictly massless or have a negligible VEV. In N-flation and generic assisted inflation scenarios, inflation is a co-operative phenomenon driven by N participator fields, none of which could drive inflation on their own. We review upper bounds on N, as a function of the inflationary Hubble scale H. We then consider stochastic and eternal inflation in models with N participator fields showing that individual fields may evolve stochastically while the whole ensemble behaves deterministically, and that a wide range of eternal inflationary scenarios are possible in this regime. We then compute one-loop quantum corrections to the inflationary power spectrum. These are largest with N spectator fields and a single participator field, and the resulting bound on N is always weaker than those obtained in other ways. We find that loop corrections to the N-flation power spectrum do not scale with N, and thus place no upper bound on the number of participator fields. This result also implies that, at least to leading order, the theory behaves like a composite single scalar field. In order to perform this calculation, we address a number of issues associated with loop calculations in the Schwinger-Keldysh "in-in" formalism.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Filaments which form either between or around active regions (ARs) are called intermediate filaments. In spite of various theoretical studies, the origin of the chirality of filaments is still uncovered. We investigated how intermediate filaments are related to their associated ARs, especially from the point of view of magnetic helicity and the orientation of polarity inversion lines (PILs). The chirality of filaments has been determined based on the orientations of barbs observed in BBSO full-disk Halpha images taken during the rising phase of solar cycle 23. The sign of magnetic helicity of ARs has been determined using S/inverse-S shaped sigmoids from Yohkoh SXT images. As a result, we have found a good correlation between the chirality of filaments and the magnetic helicity sign of ARs. Among 45 filaments, 42 filaments have shown the same sign as helicity sign of nearby ARs. It has been also confirmed that the role of both the orientation and the relative direction of PILs to ARs in determining the chirality of filaments is not significant, against a theoretical prediction. These results suggest that the chirality of intermediate filaments may originate from magnetic helicity of their associated ARs.

We consider the impact of thermal inflation -- a short, secondary period of inflation that can arise in supersymmetric scenarios -- on the stochastic gravitational wave background. We show that while the primordial inflationary gravitational wave background is essentially unchanged at CMB scales, it is massively diluted at solar system scales and would be unobservable by a BBO style experiment. Conversely, bubble collisions at the end of thermal inflation can generate a new stochastic background. We calculate the likely properties of the bubbles created during this phase transition, and show that the expected amplitude and frequency of this signal would fall within the BBO range.

Inflation driven by a single, minimally coupled, slowly rolling field generically yields a negligible primordial non-Gaussianity. We discuss two distinct mechanisms by which a non-trivial potential can generate large non-Gaussianities. Firstly, if the inflaton traverses a feature in the potential, or if the inflationary phase is short enough so that initial transient contributions to the background dynamics have not been erased, modes near horizon-crossing can acquire significant non-Gaussianities. Secondly, potentials with small-scale structure may induce significant non-Gaussianities while the relevant modes are deep inside the horizon. The first case includes the "step" potential we previously analyzed while the second "resonance" case is novel. We derive analytic approximations for the 3-point terms generated by both mechanisms written as products of functions of the three individual momenta, permitting the use of efficient analysis algorithms. Finally, we present a significantly improved approach to regularizing and numerically evaluating the integrals that contribute to the 3-point function.

Parametric resonance or preheating is a plausible mechanism for bringing about the transition between the inflationary phase and a hot, radiation dominated universe. This epoch results in the rapid production of heavy particles far from thermal equilibrium and could source a significant stochastic background of gravitational radiation. Here, we present a numerical algorithm for computing the contemporary power spectrum of gravity waves generated in this post-inflationary phase transition for a large class of scalar-field driven inflationary models. We explicitly calculate this spectrum for both quartic and quadratic models of chaotic inflation, and low-scale hybrid models. In particular, we consider hybrid models with an ``inverted'' potential. These models have a very short and intense period of resonance which is qualitatively different from previous examples studied in this context, but we find that they lead to a similar spectrum of gravitational radiation.

We consider gravitational wave production due to parametric resonance at the end of inflation, or ``preheating''. This leads to large inhomogeneities which source a stochastic background of gravitational waves at scales inside the comoving Hubble horizon at the end of inflation. We confirm that the present amplitude of these gravitational waves need not depend on the inflationary energy scale. We analyze an explicit model where the inflationary energy scale is ~10^9 GeV, yielding a signal close to the sensitivity of Advanced LIGO and BBO. This signal highlights the possibility of a new observational ``window'' into inflationary physics, and provides significant motivation for searches for stochastic backgrounds of gravitational waves in the Hz to GHz range, with an amplitude on the order of \Omega_{gw}(k)h^2 ~ 10^-11. Finally, the strategy used in our numerical computations is applicable to the gravitational waves generated by many inhomogeneous processes in the early universe.

We compute the 3-point correlation function for a general model of inflation driven by a single, minimally coupled scalar field. Our approach is based on the numerical evaluation of both the perturbation equations and the integrals which contribute to the 3-point function. Consequently, we can analyze models where the potential has a "feature", in the vicinity of which the slow roll parameters may take on large, transient values. This introduces both scale and shape dependent non-Gaussianities into the primordial perturbations. As an example of our methodology, we examine the ``step'' potentials which have been invoked to improve the fit to the glitch in the $<TT>$ $C_l$ for $l \sim 30$, present in both the one and three year WMAP data sets. We show that for the typical parameter values, the non-Gaussianities associated with the step are far larger than those in standard slow roll inflation, and may even be within reach of a next generation CMB experiment such as Planck. More generally, we use this example to explain that while adding features to potential can improve the fit to the 2-point function, these are generically associated with a greatly enhanced signal at the 3-point level. Moreover, this 3-point signal will have a very nontrivial shape and scale dependence, which is correlated with the form of the 2-point function, and may thus lead to a consistency check on the models of inflation with non-smooth potentials.

We consider the long standing puzzle of how to obtain meaningful probabilities in eternal inflation. We demonstrate a new algorithm to compute the probability distribution of pocket universe types, given a multivacua inflationary potential. The computed probability distribution is finite and manifestly gauge-independent. We argue that in some scenarios this technique can be applied to disfavor some eternally inflating potentials.

In many models of inflation, the period of accelerated expansion ends with preheating, a highly non-thermal phase of evolution during which the inflaton pumps energy into a specific set of momentum modes of field(s) to which it is coupled. This necessarily induces large, transient density inhomogeneities which can source a significant spectrum of gravitational waves. In this paper, we consider the generic properties of gravitational waves produced during preheating, perform detailed calculations of the spectrum for several specific inflationary models, and identify problems that require further study. In particular, we argue that if these gravitational waves exist they will necessarily fall within the frequency range that is feasible for direct detection experiments -- from laboratory through to solar system scales. We extract the gravitational wave spectrum from numerical simulations of preheating after $\lambda \phi^4$ and $m_{\phi}^2 \phi^2$ inflation, and find that they lead to a gravitational wave amplitude of around $\Omega_{gw}h^2\sim 10^{-10}$. This is considerably higher than the amplitude of the primordial gravitational waves produced during inflation. However, the typical wavelength of these gravitational waves is considerably shorter than LIGO scales, although in extreme cases they may be visible at scales accessible to the proposed BBO mission. We survey possible experimental approaches to detecting any gravitational wave background generated during preheating.

We study gravitationally bound static and spherically symmetric configurations of k-essence fields. In particular, we investigate whether these configurations can reproduce the properties of dark matter haloes. The classes of Lagrangians we consider lead to non-isotropic fluids with barotropic and polytropic equations of state. The latter include microscopic realizations of the often-considered Chaplygin gases, which we find can cluster into dark matter halo-like objects with flat rotation curves, while exhibiting a dark energy-like negative pressure on cosmological scales. We complement our studies with a series of formal general results about the stability and initial value formulation of non-canonical scalar field theories, and we also discuss a new class of de Sitter solutions with spacelike field gradients.

We consider the gravitational effects of a single, fixed-norm, Lorentz-violating timelike vector field. In a cosmological background, such a vector field acts to rescale the effective value of Newton's constant. The energy density of this vector field precisely tracks the energy density of the rest of the universe, but with the opposite sign, so that the universe experiences a slower rate of expansion for a given matter content. This vector field similarly rescales Newton's constant in the Newtonian limit, although by a different factor. We put constraints on the parameters of the theory using the predictions of primordial nucleosynthesis, demonstrating that the norm of the vector field should be less than the Planck scale by an order of magnitude or more.

We investigate the perturbation theory of a fixed-norm, timelike Lorentz-violating vector field. After consistently quantizing the vector field to put constraints on its parameters, we compute the primordial spectra of perturbations generated by inflation in the presence of this vector field. We find that its perturbations are sourced by the perturbations of the inflaton; without the inflaton perturbation the vector field perturbations decay away leaving no primordial spectra of perturbations. Since the inflaton perturbation does not have a spin-1 component, the vector field generically does not generate any spin-1 ``vector-type'' perturbations. Nevertheless, it will modify the amplitude of both the spin-0 ``scalar-type'' and spin-2 ``tensor-type'' perturbation spectra, leading to violations of the inflationary consistency relationship.

We investigate the ability to constrain oscillatory features in the primordial power spectrum using current and future cosmic microwave background observations. In particular, we study the observability of an oscillation arising from imprints of physics at the cut-off energy scale. We perform a likelihood analysis on the WMAP data set, and find that the current data set constrains the amplitude of the oscillations to be less than 0.77 at 2-sigma, consistent with a power spectrum without oscillations. In addition, we investigate the fundamental limitations in the measurement of oscillation parameters by studying the constraints from a cosmic variance limited experiment. We find that such an experiment is capable of constraining the amplitude of such oscillations to be below 0.005, implying that reasonable models with cut-off energy scales Lambda>200 H_infl are unobservable through the microwave background.

We propose a new alternative mechanism to seed a scale invariant spectrum of primordial density perturbations that does not rely on inflation. In our scenario, a perfect fluid dominates the early stages of an expanding, non-inflating universe. Because the speed of sound of the fluid decays, perturbations are left frozen behind the sound horizon, with a spectral index that depends on the fluid equation of state. We explore here a toy model that realizes this idea. Although the model can explain an adiabatic, Gaussian, scale invariant spectrum of primordial perturbations, it turns out that in its simplest form it cannot account for the observed amplitude of the primordial density perturbations.

In the presence of a short-distance cutoff, the choice of a vacuum state in an inflating, non-de Sitter universe is unavoidably ambiguous. The ambiguity is related to the time at which initial conditions for the mode functions are specified and to the way the expansion of the universe affects those initial conditions. In this paper we study the imprint of these uncertainties on the predictions of inflation. We parametrize the most general set of possible vacuum initial conditions by two phenomenological variables. We find that the generated power spectrum receives oscillatory corrections whose amplitude is proportional to the Hubble parameter over the cutoff scale. In order to further constrain the phenomenological parameters that characterize the vacuum definition, we study gravitational particle production during different cosmological epochs.