An analytic sun pillar model is developed which indicates that large hexagonal ice columns can cause sun pillars. The model shows that the Stuchtey method for explaining sun pillars is not incorrect, only incomplete. The model uses an expression for the intensity of the sun pillar that considers both the optical mapping of light for a single crystal orientation and for the probability of having that particular orientation. The model, which is extended to a sun of finite size, clearly shows that the portion of the pillar that is seen above the horizon is brightest when the sun is below the horizon.
Alistair B. Fraser and Gary J. Thompson, "Analytic sun pillar model," J. Opt. Soc. Am. 70, 1145-1148 (1980)