We determine the distribution of cosmic string loops directly from simulations, rather than determining the loop production function and inferring the loop distribution from that. For a wide range of loop lengths, the results agree well with a power law exponent -2.5 in the radiation era and -2 in the matter era, the universal result for any loop production function that does not diverge at small scales. Our results extend those of Ringeval, Sakellariadou, and Bouchet: we are able to run for 15 times longer in conformal time and simulate a volume 300-2400 times larger. At the times they reached, our simulation is in general agreement with the more negative exponents they found, -2.6 and -2.4. However, our simulations show that this was a transient regime; at later times the exponents decline to the values above. This provides further evidence against models with a rapid divergence of the loop density at small scales, such as ``model 3'' used to analyze LIGO data and predict LISA sensitivity.

A network of cosmic strings would lead to gravitational waves which may be detected by pulsar timing or future interferometers. The details of the gravitational wave signal depend on the distribution of cosmic string loops, which are produced by intercommutations from the scaling network of long strings. We analyze the limits imposed by energy conservation, i.e., by the fact that the total amount of string flowing into loops cannot exceed the amount leaving the long strings. We show that some recent suggestions for the cosmic string loop production rate and distribution are ruled out by these limits. As a result, gravitational waves based on such suggestions, in particular "model 3" used in LIGO data analysis, are not to be expected.

To understand the properties of a possible cosmic string network requires knowledge of the structures on long strings, which control the breaking off of smaller loops. These structures are influenced by backreaction due to gravitational wave emission. To gain an understanding of this process, we calculate the effect of gravitational backreaction on the "helical breather": an infinite cosmic string with a helical standing wave. We do the calculation analytically to first order in the amplitude of the oscillation. Our results for the rate of loss of length agree with the total power emitted from this system as calculated by Sakellariadou. We also find a rotation of the generators of the string that leads to an advancement of the phase of the oscillation in addition to that produced by the loss of length.

We present the results of computational gravitational backreaction on simple models of cosmic string loops. These results give us insight into the general behavior of cusps and kinks on loops, in addition to other features of evolution. Kinks are rounded off via an asymmetric and divergent correction to the string direction. The result is that cusps emerge in the place of kinks but the resulting smooth string section has a small amount of energy. Existing cusps persist, but quickly lose strength as backreaction removes energy from the string surrounding the cusp. Both kinks and cusps have their location in space shifted slightly with each oscillation.

We find the leading-order effect of gravitational back-reaction on cosmic strings for points near kinks and cusps. Near a kink, the effect diverges as the inverse cube root of the distance to the kink, and acts in a direction transverse to the worldsheet. Over time the kink is rounded off, but only regions fairly close to the kink are significantly affected. Near cusps, the effect diverges inverse linearly with the distance to the cusp, and acts against the direction of the cusp motion. This results in a fractional loss of string energy that diverges logarithmically with the distance of closest approach to the cusp.

We analyze the effect of gravitational back reaction on cosmic string loops with kinks, which is an important determinant of the shape, and thus the potential observability, of string loops which may exist in the universe today. Kinks are not rounded off, but may be straightened out. This means that back reaction will only cause loops with kinks to develop cusps after some potentially large fraction of their lifetimes. In some loops, symmetries prevent even this process, so that the loop evaporates in a self-similar fashion and the kinks are unchanged. As an example, we discuss back-reaction on the rectangular Garfinkle-Vachaspati loop.

We do a complete calculation of the stochastic gravitational wave background to be expected from cosmic strings. We start from a population of string loops taken from simulations, smooth these by Lorentzian convolution as a model of gravitational back reaction, calculate the average spectrum of gravitational waves emitted by the string population at any given time, and propagate it through a standard model cosmology to find the stochastic background today. We take into account all known effects, including changes in the number of cosmological relativistic degrees of freedom at early times and the possibility that some energy is in rare bursts that we might never have observed.

We combine new analysis of the stochastic gravitational wave background to be expected from cosmic strings with the latest pulsar timing array (PTA) limits to give an upper bound on the energy scale of the possible cosmic string network, $G\mu < 1.5\times 10^{-11}$ at the 95% confidence level. We also show bounds from LIGO and to be expected from LISA and BBO. Current estimates for the gravitational wave background from supermassive black hole binaries are at the level where a PTA detection is expected. But if PTAs do observe a background soon, it will be difficult in the short term to distinguish black holes from cosmic strings as the source, because the spectral indices from the two sources happen to be quite similar. If PTAs do not observe a background, then the limits on $G\mu$ will improve somewhat, but a string network with $G\mu$ substantially below $10^{-11}$ will produce gravitational waves primarily at frequencies too high for PTA observation, so significant further progress will depend on intermediate-frequency observatories such as LISA, DECIGO and BBO.

We calculate the metric and affine connection due to a piecewise linear cosmic string loop, and the effect of gravitational back reaction for the Garfinkle-Vachaspati loop with four straight segments. As expected, back reaction reduces the size of the loop, in accord with the energy going into gravitational waves. The "square" (maximally symmetric) loop evaporates without changing shape, but for all other loops in this class, the kinks become less sharp and segments between kinks become curved. If the loop is close to the square case, it will evaporate before its kinks are significantly changed, if it is far from square, the opening out of the kinks is much faster than evaporation of the loop.

We analyze the shapes of cosmic string loops found in large-scale simulations of an expanding-universe string network. The simulation does not include gravitational back reaction, but we model that process by smoothing the loop using Lorentzian convolution. We find that loops at formation consist of generally straight segments separated by kinks. We do not see cusps or any cusp-like structure at the scale of the entire loop, although we do see very small regions of string that move with large Lorentz boosts. However, smoothing of the string almost always introduces two cusps on each loop. The smoothing process does not lead to any significant fragmentation of loops that were in non-self-intersecting trajectories before smoothing.

Charge carriers moving at the speed of light along a straight, superconducting cosmic string carry with them a logarithmically divergent slab of electromagnetic field energy. Thus no finite local input can induce a current that travels unimpeded to infinity. Rather, electromagnetic back-reaction must damp this current asymptotically to nothing. We compute this back-reaction and find that the electromagnetic fields and currents decline exactly as rapidly as necessary to prevent a divergence. We briefly discuss the corresponding gravitational situation.

Using recent simulation results, we provide the mass and speed spectrum of cosmic string loops. This is the quantity of primary interest for many phenomenological signatures of cosmic strings, and it can be accurately predicted using recently acquired detailed knowledge of the loop production function. We emphasize that gravitational smoothing of long strings does not play any role in determining the total number of existing loops. We derive a bound on the string tension imposed by recent constraints on the stochastic gravitational wave background from pulsar timing arrays, finding $G\mu \leq 2.8\times 10^{-9}$. We also provide a derivation of the Boltzmann equation for cosmic string loops in the language of differential forms.

Core-collapse supernovae emit on the order of 3x10^53 ergs in high-energy neutrinos over a time of order 10 seconds, and so decrease their mass by about 0.2 solar mass. If the explosion is nearly spherically symmetric, there will be little gravitational wave emission. Nevertheless, the sudden decrease of mass of the progenitor may cause a change in the gravitational time delay of signals from a nearby pulsar. We calculate the change in arrival times as successive pulses pass through the neutrino shell at different times, and find that the effect may be detectable in ideal circumstances.

Given a sufficiently large universe, numerous civilizations almost surely exist. Some of these civilizations will be short-lived and die out relatively early in their development, i.e., before having the chance to spread to other planets. Others will be long-lived, potentially colonizing their galaxy and becoming enormous in size. What fraction of civilizations in the universe are long-lived? The "universal doomsday" argument states that long-lived civilizations must be rare because if they were not, we should find ourselves living in one. Furthermore, because long-lived civilizations are rare, our civilization's prospects for long-term survival are poor. Here, we develop the formalism required for universal doomsday calculations and show that while the argument has some force, our future is not as gloomy as the traditional doomsday argument would suggest, at least when the number of early existential threats is small.

An eternally inflating universe produces an infinite amount of spatial volume, so every possible event happens an infinite number of times, and it is impossible to define probabilities in terms of frequencies. This problem is usually addressed by means of a measure, which regulates the infinities and produces meaningful predictions. I argue that any measure should obey certain general axioms, but then give a simple toy model in which one can prove that no measure obeying the axioms exists. In certain cases of eternal inflation there are measures that obey the axioms, but all such measures appear to be unacceptable for other reasons. Thus the problem of defining sensible probabilities in eternal inflation seems not be solved.

Using a new parallel computing technique, we have run the largest cosmic string simulations ever performed. Our results confirm the existence of a long transient period where a non-scaling distribution of small loops is produced at lengths depending on the initial correlation scale. As time passes, this initial population gives way to the true scaling regime, where loops of size approximately equal to one-twentieth the horizon distance become a significant component. We observe similar behavior in matter and radiation eras, as well as in flat space. In the matter era, the scaling population of large loops becomes the dominant component; we expect this to eventually happen in the other eras as well.

We develop a "semi-parallel" simulation technique suggested by Pretorius and Lehner, in which the simulation spacetime volume is divided into a large number of small 4-volumes which have only initial and final surfaces. Thus there is no two-way communication between processors, and the 4-volumes can be simulated independently without the use of MPI. This technique allows us to simulate much larger volumes than we otherwise could, because we are not limited by total memory size. No processor time is lost waiting for other processors. We compare a cosmic string simulation we developed using the semi-parallel technique with our previous MPI-based code for several test cases and find a factor of 2.6 improvement in the total amount of processor time required to accomplish the same job for strings evolving in the matter-dominated era.

Superconducting cosmic strings naturally emit highly boosted charge carriers from cusps. This occurs when a cosmic string or a loop moves through a magnetic field and develops an electric current. The charge carriers and the products of their decay, including protons, photons and neutrinos, are emitted as a narrow jets with opening angle $\theta \sim 1/\gamma_c$, where $\gamma_c$ is the Lorentz factor of the cusp. The excitation of electric currents in strings occurs mostly in clusters of galaxies, which are characterized by magnetic fields $B \sim 10^{-6}$ G and a filling factor $f_B \sim 10^{-3}$. Two string parameters determine the emission of the particles: the symmetry breaking scale $\eta$, which for successful applications should be of order $10^9$--$10^{12}$ GeV, and the dimensionless parameter $i_c$, which determines the maximum induced current as $J_{max} =i_c e \eta$ and the energy of emitted charge carriers as $\epsilon_x \sim i_c \gamma_c \eta$, where $e$ is the electric charge of a particle. For the parameters $\eta $ and $B$ mentioned above, the Lorentz factor reaches $\gamma_c \sim 10^{12}$ and the maximum particle energy can be as high as $\gamma_c\eta \sim 10^{22}$ GeV. The diffuse fluxes of UHE neutrinos are close to the cascade upper limit, and can be detected by future neutrino observatories. The signatures of this model are: very high energies of neutrinos, in excess of $10^{20}$ eV, correlation of neutrinos with clusters of galaxies, simultaneous appearance of several neutrino-produced showers in the field of view of very large detectors, such as JEM-EUSO, and 10 TeV gamma radiation from the Virgo cluster. The flux of UHE protons from cusps may account for a large fraction of the observed events at the highest energies.

A sequence of two symmetry breaking transitions in the early universe may produce monopoles whose flux is confined into two strings each, which thus assemble into "necklaces" with monopoles as beads. Such "cosmic necklaces" have been proposed as a source of ultra-high-energy cosmic rays. We analyze the evolution of these systems and show that essentially all monopoles annihilate or leave the string at early times, after which cosmic necklaces evolve in a similar way to a network of ordinary cosmic strings. We investigate several modifications to the basic picture, but in nearly all cases we find that too few monopoles remain on the necklaces to produce any observable cosmic rays. There may be a small window for superconducting condensates to prevent annihilations, but only if both the string and the condensate scale are very high.

We study the formation of cosmic strings by confining a stochastic magnetic field into flux tubes in a numerical simulation. We use overdamped evolution in a potential that is minimized when the flux through each face in the simulation lattice is a multiple of the fundamental flux quantum. When the typical number of flux quanta through a correlation-length-sized region is initially about 1, we find a string network similar to that generated by the Kibble-Zurek mechanism. With larger initial flux, the loop distribution and the Brownian shape of the infinite strings remain unchanged, but the fraction of length in infinite strings is increased. A 2D slice of the network exhibits bundles of strings pointing in the same direction, as in earlier 2D simulations. We find, however, that strings belonging to the same bundle do not stay together in 3D for much longer than the correlation length. As the initial flux per correlation length is decreased, there is a point at which infinite strings disappear, as in the Hagedorn transition.

We study the production of loops in the cosmic string network in the expanding background by means of a numerical simulation exact in the flat-spacetime limit and first-order in the expansion rate. We find an initial regime characterized by production of small loops at the scale of the initial correlation length, but later we see the emergence of a scaling regime of loop production. This qualitatively agrees with earlier expectations derived from the results of flat-spacetime simulations. In the final scaling regime we find that the characteristic length of loops scales as $\sim 0.1 t$ in both radiation and matter eras.

We study the spectrum of loops as a part of a complete network of cosmic strings in flat spacetime. After a long transient regime, characterized by production of small loops at the scale of the initial conditions, it appears that a true scaling regime takes over. In this final regime the characteristic length of loops scales as $0.1 t$, in contrast to earlier simulations which found tiny loops. We expect the expanding-universe behavior to be qualitatively similar. The large loop sizes have important cosmological implications. In particular, the nucleosynthesis bound becomes $G\mu \lesssim 10^{-7}$, much tighter than before.

Loops formed from a cosmic string network at early times would act as seeds for early formation of halos, which would form galaxies and lead to early reionization. With reasonable guesses about astrophysical and string parameters, the cosmic string scale $G\mu$ must be no more than about $3\times 10^{-8}$ to avoid conflict with the reionization redshift found by WMAP. The bound is much stronger for superstring models with a small string reconnection probability. For values near the bound, cosmic string loops may explain the discrepancy between the WMAP value and theoretical expectations.

We present general solutions to the equations of motion for a superconducting relativistic chiral string that satisfy the unit magnitude constraint in terms of products of rotations. From this result we show how to construct a general family of odd harmonic superconducting chiral loops. We further generalise the product of rotations to an arbitrary number of dimensions.

We present a method for the calculation of the gravitational back reaction cutoff on the smallest scales of cosmic string networks taking into account that not all modes on strings interact with all other modes. This results in a small scale structure cutoff that is sensitive to the initial spectrum of perturbations present on strings. From a simple model, we compute the cutoffs in radiation- and matter-dominated universes.

We calculate the gravitational radiation emitted by an infinite cosmic string with two oppositely moving wave-trains, in the small amplitude approximation. After comparing our result to the previously studied cases we extend the results to a new regime where the wavelengths of the opposing wave-trains are very different. We show that in this case the amount of power radiated vanishes exponentially. This means that small excitations moving in only one direction may be very long lived, and so the size of the smallest scales in a string network might be much smaller than what one would expect from gravitational back reaction. This result allows for a potential host of interesting cosmological possibilities involving ultra-high energy cosmic rays, gamma ray bursts and gravitational wave bursts.

It is difficult for conventional sources to accelerate cosmic ray particles to the highest energies that have been observed. Topological defects such as monopoles and strings overcome this difficulty, because their natural energy scale is at or above the observed energies. Monopoles connected by strings are a particularly attractive source, because they would cluster in the galactic halo and thus explain the absence of the GZK cutoff. Heavy monopoles connected by light strings could last for the age of the universe as required. Further observations might support this model by detection of the anisotropy due to the halo, or might refute such models if strong clustering of arrival directions or correlations with known astrophysical objects are confirmed. All top-down models must contend with recent claims that the percentage of photons among the cosmic rays is smaller than such models predict.

We rederive, using an elementary formalism, the general solution to the equations of motion for a superconducting string with a chiral (null) neutral current, earlier obtained by Carter and Peter. We apply this solution to show that the motion of such string loops is strictly periodic and analyze cusp-like behavior and vorton solutions of arbitrary shape. We argue that this solution can be used to approximately describe the dynamics of superconducting strings with small non-chiral currents. We use this description to estimate the electromagnetic radiation power from such strings.

In about half of all near-cusp events in superconducting cosmic strings with chiral currents (and probably with general currents as well), the string intersects itself near the cusp. Intercommutation causes the conversion of the string near the cusp into a vorton (in the chiral case) with very high Lorentz boost. We demonstrate how to analyze the cusp shape in a Lorentz frame that makes the motion simple, by the use of a 5-dimensional procedure, and analyze the resulting production of vortons.

Cusps in superconducting cosmic strings produce strongly beamed electromagnetic radiation. To calculate the energy emitted requires taking into account the effect of the charge carriers on the string motion, which has previously been done only heuristically. Here, we use the known exact solution to the equations of motion for the case where the current is chiral to update previous calculations for the total energy, spectrum and angular distribution in that case. We analyze the dependence of the radiated energy on the cusp parameters, and discuss which types of cusp dominate the total radiation emitted from an ensemble.

We perform numerical simulations of cosmic necklaces (systems of monopoles connected to two strings each) and investigate the conditions under which monopoles annihilate. When the total monopole energy is large compared to the string energy, we find that the string motion is no longer periodic, and thus the strings will be chopped up by self intersection. When the total monopole energy is much smaller than the string energy, the string motion is periodic, but that of the monopoles is not, and thus the monopoles travel along the string and annihilate with each other.

We discuss the large-scale structure of the universe in inflationary cosmology and the implications that it may have for the long-term future of civilizations. Although each civilization is doomed to perish, it may be possible to transmit its accumulated knowledge to future civilizations. We consider several scenarios of this sort. If the cosmological constant is positive, it eventually dominates the universe and bubbles of inflationary phase begin to nucleate at a constant rate. Thermalized regions inside these inflating bubbles will give rise to new galaxies and civilizations. It is possible in principle to send a message to one of them. It might even be possible to send a device whose purpose is to recreate an approximation of the original civilization in the new region. However, the message or device will almost certainly be intercepted by black holes, which nucleate at a much higher rate than inflating bubbles. Formation of new inflating regions can also be triggered by gravitational collapse, but again the probability is low, and the number of attempts required for a positive outcome is enormous. The probability can be higher if the energy scale of inflation is closer to the Planck scale, but a high energy scale produces a tight bound on the amount of information that can be transmitted. One can try to avoid quantum tunneling altogether, but this requires a violation of quantum inequalities which constrain the magnitude of negative energy densities. However, the limits of validity of quantum inequalities are not clear, and future research may show that the required violation is in fact possible. Therein lies the hope for the future of civilizations.

We have simulated large-amplitude standing waves on an Abelian-Higgs cosmic string in classical lattice field theory. The radiation rate falls exponentially with wavelength, as one would expect from the field profile around a gauge string. Our results agree with those of Shellard and Moore, but not those of Vincent, Antunes, and Hindmarsh. The radiation rate falls too rapidly to sustain a scaling solution via direct radiation of particles from string length. There is thus reason to doubt claims of strong constraints on cosmic string theories from cosmic ray observations.

The electromagnetic decay and final annihilation of magnetic monopole-antimonopole pairs formed in the early universe has been proposed as a possible mechanism to produce the highest energy cosmic rays. We show that for a monopole abundance saturating the Parker limit, the density of magnetic monopolonium formed is many orders of magnitude less than that required to explain the observed cosmic ray flux. We then propose a different scenario in which the monopoles and antimonopoles are connected by strings formed at a low energy phase transition (~ 100 GeV). The bound states decay by gravitational radiation, with lifetimes comparable with the age of the universe. This mechanism avoids the problems of the standard monopolonium scenario, since the binding of monopoles and antimonopoles is perfectly efficient.

We have performed a lattice field theory simulation of cusps in Abelian-Higgs cosmic strings. The results are in accord with the theory that the portion of the strings which overlaps near the cusp is released as radiation. The radius of the string cores which must touch to produce the evaporation is approximately $r = 1$ in natural units. In general, the modifications to the string shape due to the cusp may produce many cusps later in the evolution of a string loop, but these later cusps will be much smaller in magnitude and more closely resemble kinks.

We classify the possible shapes of cosmic string cusps and how they transform under Lorentz boosts. A generic cusp can be brought into a form in which the motion of the cusp tip lies in the plane of the cusp. The cusp whose motion is perpendicular to this plane, considered by some authors, is a special case and not the generic situation. We redo the calculation of the energy in the region where the string overlaps itself near a cusp, which is the maximum energy that can be released in radiation. We take into account the motion of a generic cusp and the resulting Lorentz contraction of the string core. The result is that the energy scales as $\sqrt {rL}$ instead of the usual value of $r^{1/3} L^{2/3}$, where $r$ is the string radius and $L$ and is the typical length scale of the string. Since $r << L $ for cosmological strings, the radiation is strongly suppressed and could not be observed.