Systems that exhibit pattern formation are typically driven and dissipative. In the early universe, parametric resonance can drive explosive particle production called preheating. The fields that are populated then decay quantum mechanically if their particles are unstable. Thus, during preheating, a driven-dissipative system exists. We have shown previously that pattern formation can occur in two dimensions in a self-coupled inflaton system undergoing parametric resonance. In this paper, we provide evidence of pattern formation for more realistic initial conditions in both two and three dimensions. In the one-field case, we have the novel interpretation that these patterns can be thought of as a network of domain walls. We also show that the patterns are spatio-temporal, leading to a distinctive, but probably low-amplitude peak in the gravitational wave spectrum. In the context of a two-field model, we discuss putting power from resonance into patterns on cosmological scales, in particular to explain the observed excess power at 100 h^{-1}Mpc, but why this seems unlikely in the absence of a period of post-preheating inflation. Finally we note our model is similar to that of the decay of DCCs and therefore pattern formation may also occur at RHIC and LHC.

We study the interaction between monopoles and embedded domain walls in a O(3) linear sigma model. We discover that there is an attractive force between the monopole and the wall. We provide evidence that after the monopole and domain wall collide, the monopole unwinds on the wall, and that the winding number spreads out on the surface. These results support the suggestion by Dvali et al. (hep-ph/9710301) that the monopole problem can be solved by monopole-domain wall interactions.

The formation of regular patterns is a well-known phenomenon in condensed matter physics. Systems that exhibit pattern formation are typically driven and dissipative with pattern formation occurring in the weakly non-linear regime and sometimes even in more strongly non-linear regions of parameter space. In the early universe, parametric resonance can drive explosive particle production called preheating. The fields that are populated then decay quantum mechanically if their particles are unstable. Thus, during preheating, a driven-dissipative system exists. In this paper, we show that a self-coupled inflaton oscillating in its potential at the end of inflation can exhibit pattern formation.

We numerically investigate the decay, via parametric resonance, of the inflaton with an m^2 phi^2 potential into a scalar matter field with a symmetry breaking potential. We consider the case where symmetry breaking takes place during inflation. We show that when expansion is not taken into account symmetry restoration and non-thermal defect production during reheating is possible. However in an expanding universe the fields do not spend sufficient time in the instability bands; thus symmetry restoration and subsequent domain wall production do not occur.

We present results of a cosmological hydrodynamical study of gravitational accretion in cosmic string wakes and filaments. Cosmic string wakes are formed by fast moving (v ~ c) strings. A conical deficit angle in the string spacetime induces a velocity perturbation in the background matter and a two-dimensional wake accretes in the path of the string. Filaments are formed by slow moving strings with a large amount of small scale structure. The major gravitational perturbation from slow moving strings is due to the Newtonian field induced by the effective mass of the wiggles. In cosmic string wakes, cool streams of baryons collide and are trapped at the center of the wake causing an enhancement of baryons versus dark matter by a factor of 2.4. In filaments, a high pressure is induced at the filament core and baryonic matter is expelled leading to a baryon deficit in the center of the filament.

We study the accretion of cold and hot dark matter onto a cosmic string wake using a high-resolution numerical simulation. We verify previous analytical calculations predicting the radius of bound matter around wakes and inflow velocities of the dark matter, as well as assumptions about the self-similarity of the solution. In cold wakes, we show that self-similarity is approached quickly after a `binding' transition. Hot wakes approach self-similarity rapidly once the free streaming `pressure' falls below a critical value and accretion begins. We also analytically calculate the size of the overdensity in wakes with cold dark matter and compare the results to simulations. We remark that the results derived may be used in testing gravitational numerical codes in the non-linear regime.

The ordering of scalar fields after a phase transition in which a group $G$ of global symmetries is spontaneously broken to a subgroup $H$ provides a possible explanation for the origin of structure in the universe, as well as leading to observable effects in condensed matter systems. The field dynamics can depend in principle on the geometry and topology of the vacuum manifold G/H; for example, texture configurations which collapse and unwind will exist if the third homotopy group $\pi_3(G/H)$ is nontrivial. We numerically simulate the evolution of texture-like configurations in a number of different models, in order to determine the extent to which the geometry and topology of the vacuum manifold influences the field evolution. We find that the dynamics is affected by whether or not the theory supports strings or monopoles [characterized by $\pi_1(G/H)$ and $\pi_2(G/H)$, respectively]. In some of the theories studied, configurations with initially spherically symmetric energy densities are unstable to nonspherical collapse; these theories are also found to nucleate defects during the collapse. Models that do not support monopoles or strings behave similarly to each other, regardless of the specific vacuum manifold.

We present a method for integrating the cosmological hydrodynamical equations including a collisionless dark matter component. For modelling the baryonic matter component, we use the Piecewise Parabolic Method (PPM) which is a high-accuracy shock capturing technique. The dark matter component is modeled using gravitationally interacting particles whose evolution is determined using standard particle-in-cell techniques. We discuss details of the inclusion of gravity and expansion in the PPM code and give results of a number of tests of the code. This code has been developed for a massively parallel, SIMD supercomputer: the MasPar MP-2 parallel processor. We present details of the techniques we have used to implement the code for this architecture and discuss performance of the code on the MP-2. The code processes $5.0 \times 10^4$ grid zones per second and requires 53 seconds of machine time for a single timestep in a $128^3$ simulation.

The clustering of baryons and cold dark matter induced by a single moving string is analyzed numerically making use of a new three-dimensional Eulerian cosmological hydro code$^{1)}$ which is based on the PPM method to track the baryons and the PIC method to evolve the dark matter particles. A long straight string moving with a speed comparable to $c$ induces a planar overdensity (a``wake"). Since the initial perturbation is a velocity kick towards the plane behind the string and there is no initial Newtonian gravitational line source, the baryons are trapped in the center of the wake, leading to an enhanced baryon to dark matter ratio. The cold coherent flow leads to very low post--shock temperatures of the baryonic fluid. In contrast, long strings with a lot of small-scale structure (which can be described by adding a Newtonian gravitational line source) move slowly and form filamentary objects. The large central pressure due to the gravitational potential causes the baryons to be expelled from the central regions and leads to a relative deficit in the baryon to dark matter ratio. In this case, the velocity of the baryons is larger, leading to high post-shock temperatures.

This study presents a simplified approach to studying the dynamics of global texture collapse. We derive equations of motion for a spherically symmetric field configuration using a two parameter ansatz. Then we analyse the effective potential for the resulting theory to understand possible trajectories of the field configuration in the parameter space of the ansatz. Numerical results are given for critical winding and collapse time in spatially flat non-expanding, and flat expanding universes. In addition, the open non-expanding and open-expanding cases are studied.

A theory of gravitation is constructed in which all homogeneous and isotropic solutions are nonsingular, and in which all curvature invariants are bounded. All solutions for which curvature invariants approach their limiting values approach de Sitter space. The action for this theory is obtained by a higher derivative modification of Einstein's theory. We expect that our model can easily be generalized to solve the singularity problem also for anisotropic cosmologies.

We give an upper estimate for the number of gamma ray bursts from ordinary (non-superconducting) cosmic strings expected to be observed at terrestrial detectors. Assuming that cusp annihilation is the mechanism responsible for the bursts we consider strings arising at a GUT phase transition and compare our estimate with the recent BATSE results. Further we give a lower limit for the effective area of future detectors designed to detect the cosmic string induced flux of gamma ray bursts.