We present the results of a detailed investigation of the low-temperature
properties of the vortex system in strongly anisotropic layered superconductors
with a random array of columnar pinning centers. Our method involves numerical
minimization of a free energy functional in terms of the time-averaged local
vortex density. It yields the detailed vortex density distribution for all
local free-energy minima, and therefore allows the computation of any desired
correlation function of the time-averaged local vortex density. Results for the
phase diagram in the temperature vs. pin concentration plane at constant
magnetic induction are presented. We confirm that for very low pin
concentrations, the low-temperature phase is a Bragg glass, which melts into an
interstitial liquid phase via two first-order steps, separated by a Bose glass
phase. At higher concentrations, however, the low-temperature phase is a Bose
glass, and the melting transition becomes continuous. The transition is then
characterized by the onset of percolation of liquid-like regions across the
sample. Inhomogeneous local melting of the Bose glass is found to occur. There
is also a depinning crossover between the interstitial liquid and a completely
unpinned liquid at higher temperatures. At sufficiently large pin
concentrations, the depinning line merges with the Bose glass to interstitial
liquid transition. Many of the features we find have been observed
experimentally and in simulations. We discuss the implications of our results
for future experimental and theoretical work.Comment: 15 pages including Figure