A compact analytical solution obtained in the paraxial approximation is used to investigate focused and unfocused vortex beams radiated by a source with a Gaussian amplitude distribution. Comparisons with solutions of the Helmholtz equation are conducted to determine bounds on the parameter space in which the paraxial approximation is accurate. A linear relation is obtained for the dependence of the vortex ring radius on the topological charge, characterized by its orbital number, in the far field of an unfocused beam and in the focal plane of a focused beam. For a focused beam, it is shown that as the orbital number increases, the vortex ring not only increases in radius but also moves out of the focal plane in the direction of the source. For certain parameters, it is demonstrated that with increasing orbital number, the maximum amplitude in a focused beam becomes localized along a spheroidal surface enclosing a shadow zone in the prefocal region. This field structure is described analytically by ray theory developed in the present work, showing that the spheroidal surface in the prefocal region coincides with a simple expression for the coordinates of the caustic surface formed in a focused vortex beam.