A prerequisite for modern radiotherapy is the ability to accurately determine the absorbed dose (D) that is given to the patient. The subject of this thesis has been to develop and evaluate efficient dose calculation models for high-energy photon beams delivered by linear accelerators. Even though the considered calculation models are general, the work has been focused on quality assurance (QA) tools used to independently verify the dose for individual treatment plans. The purpose of this verification is to guarantee patient safety and to improve the treatment outcome. Furthermore, a vital part of this work has been to explore the prospect of estimating the dose calculation uncertainties associated with individual treatment setups. A discussion on how such uncertainty estimations can facilitate improved clinical QA procedures by providing appropriate action levels has also been included within the scope of this thesis. In order to enable efficient modelling of the physical phenomena that are involved in dose output calculations it is convenient to divide them into two main categories; the first one dealing with the radiation exiting the accelerator’s treatment head and a second one associated with the subsequent energy deposition processes. A multi-source model describing the distribution of energy fluence emitted from the treatment head per delivered monitor unit (MU) is presented and evaluated through comparisons with measurements in multiple photon beams and collimator settings. The calculations show close agreement with the extensive set of experimental data, generally within +/-1% of corresponding measurements. The energy (dose) deposition in the irradiated object has been modelled through a photon pencil kernel solely based on a beam quality index (TPR20,10). This model was evaluated in a similar manner as the multi-source model at three different treatment depths. A separate study was focused on the specific difficulties associated with dose calculations in points located at a distance from the central beam axis. Despite the minimal input data required to characterize individual photon beams, the accuracy proved to be very good when comparing the calculated results with experimental data. The evaluated calculation models were finally used to analyse how well the lateral dose distributions from typical megavoltage photon beams are optimized with respect to the resulting beam flatness characteristics. The results did not reveal any obvious reasons why different manufacturers should provide different lateral dose distributions. Furthermore, the performed lateral optimizations indicate that there is room for improved flatness performance for the investigated linear accelerators.