We investigate numerically the configurational statistics of strings. The algorithm models an ensemble of global $U(1)$ cosmic strings, or equivalently vortices in superfluid $^4$He. We use a new method which avoids the specification of boundary conditions on the lattice. We therefore do not have the artificial distinction between short and long string loops or a `second phase' in the string network statistics associated with strings winding around a toroidal lattice. Our lattice is also tetrahedral, which avoids ambiguities associated with the cubic lattices of previous work. We find that the percentage of infinite string is somewhat lower than on cubic lattices, 63\% instead of 80\%. We also investigate the Hagedorn transition, at which infinite strings percolate, controlling the string density by rendering one of the equilibrium states more probable. We measure the percolation threshold, the critical exponent associated with the divergence of a suitably defined susceptibility of the string loops, and that associated with the divergence of the correlation length.

We present analytic solutions to the spatially homogeneous axion field equation, using a model potential which strongly resembles the $(1-\cos \theta)$ standard anharmonic potential. We find that the anharmonicity amplifies early universe axion density fluctuations, but only significantly so for relatively large initial misalignment angles. Enhancements of $\sim$ (2,3,4,13) result for $\theta_{\rm in}\sim (0.85,0.90,0.95,0.99)\times\pi$. The large $\theta$ result agrees with Lyth, and so validates his approximations.

We present analytic solutions to the spatially homogeneous axion field equation, using a model potential which strongly resembles the standard anharmonic $(1-\cos N\theta)$ potential, but contains only a piece-wise second order term. Our exactly soluble model for $\theta(t)$ spans the entire range $[-\pi/N,\pi/N]$. In particular, we are able to confirm (i) Turner's numeric correction factors \cite{Turner} to the adiabatic and harmonic analytic treatments of homogeneous axion oscillations, and (ii) Lyth's estimate \cite{Lyth} valid near the metastable misalignment angle $\pi/N$ at the peak of the potential. We compute the enhancement of axion density fluctuations that occurs when the axion mass becomes significant at $T\sim 1$ GeV. We find that the anharmonicity amplifies density \mbox{f}luctuations, but only significantly for relatively large initial misalignment angles. The enhancement factor is $\sim$ (2,3,4,13) for $\theta_{\rm in}\sim (0.85,0.90,0.95,0.99)\times\pi$.