The Huygens-Fresnel principle (named after Dutch physicist Christiaan Huygens and French physicist Augustin-Jean Fresnel) is a method of analysis applied to problems of wave propagation both in the far-field limit and in near-field diffraction and also reflection. It states that every point on a wavefront is itself the source of spherical wavelets, and the secondary wavelets emanating from different points mutually interfere. The sum of these spherical wavelets forms the wavefront.
In 1678, Huygens proposed that every point to which a luminous disturbance reaches becomes a source of a spherical wave; the sum of these secondary waves determines the form of the wave at any subsequent time. He assumed that the secondary waves travelled only in the "forward" direction and it is not explained in the theory why this is the case. He was able to provide a qualitative explanation of linear and spherical wave propagation, and to derive the laws of reflection and refraction using this principle, but could not explain the deviations from rectilinear propagation that occur when light encounters edges, apertures and screens, commonly known as diffraction effects. The resolution of this error was finally explained by David A. B. Miller in 1991. The resolution is that the source is a dipole (not the monopole assumed by Huygens), which cancels in the reflected direction.
In 1818, Fresnel showed that Huygens's principle, together with his own principle of interference could explain both the rectilinear propagation of light and also diffraction effects. To obtain agreement with experimental results, he had to include additional arbitrary assumptions about the phase and amplitude of the secondary waves, and also an obliquity factor. These assumptions have no obvious physical foundation but led to predictions that agreed with many experimental observations, including the Poisson spot.
Poisson was a member of the French Academy, which reviewed Fresnel's work. He used Fresnel's theory to predict that a bright spot ought to appear in the center of the shadow of a small disc, and deduced from this that the theory was incorrect. However, Arago, another member of the committee, performed the experiment and showed that the prediction was correct. (Lisle had observed this fifty years earlier.) This was one of the investigations that led to the victory of the wave theory of light over the then predominant corpuscular theory.
The Huygens-Fresnel principle provides a reasonable basis for understanding and predicting the classical wave propagation of light. However, there are limitations to the principle, namely the same approximations done for deriving the Kirchhoff's diffraction formula and the approximations of Near field due to Fresnel. These can be summarized in the fact that the wavelength of light is much smaller than the dimensions of any optical components encountered.
Kirchhoff's diffraction formula provides a rigorous mathematical foundation for diffraction, based on the wave equation. The arbitrary assumptions made by Fresnel to arrive at the Huygens-Fresnel equation emerge automatically from the mathematics in this derivation.
A simple example of the operation of the principle can be seen when an open doorway connects two rooms and a sound is produced in a remote corner of one of them. A person in the other room will hear the sound as if it originated at the doorway. As far as the second room is concerned, the vibrating air in the doorway is the source of the sound.
Not all experts agree that the Huygens' principle is an accurate microscopic representation of reality. For instance, Melvin Schwartz argued that "Huygens' principle actually does give the right answer but for the wrong reasons".
This can be reflected in the following facts:
The Huygens' principle is essentially compatible with Quantum field theory in the far field approximation, considering effective fields in the center of scattering, considering small perturbations, and in the same sense that quantum optics is compatible with Classical optics, other interpretations are subject of debates and active research.
The Feynman model where every point in an imaginary wave front as large as the room is generating a wavelet, shall also be interpreted in these approximations  and in a probabilistic context, in this context remote points can only contribute minimally to the overall probability amplitude.
Quantum field theory does not include any microscopic model for photon creation and the concept of single photon is also put under scrutiny on a theoretical level.
See Generalized Huygens' principle and the Modern Photon Wave function below.
Consider the case of a point source located at a point P0, vibrating at a frequency f. The disturbance may be described by a complex variable U0 known as the complex amplitude. It produces a spherical wave with wavelength ?, wavenumber k = 2?/?. Within a constant of proportionality, the complex amplitude of the primary wave at the point Q located at a distance r0 from P0 is:
Note that magnitude decreases in inverse proportion to the distance travelled, and the phase changes as k times the distance travelled.
Using Huygens's theory and the principle of superposition of waves, the complex amplitude at a further point P is found by summing the contributions from each point on the sphere of radius r0. In order to get agreement with experimental results, Fresnel found that the individual contributions from the secondary waves on the sphere had to be multiplied by a constant, -i/?, and by an additional inclination factor, K(?). The first assumption means that the secondary waves oscillate at a quarter of a cycle out of phase with respect to the primary wave, and that the magnitude of the secondary waves are in a ratio of 1:? to the primary wave. He also assumed that K(?) had a maximum value when ? = 0, and was equal to zero when ? = ?/2, where ? is the angle between the normal of the primary wave front and the normal of the secondary wave front. The complex amplitude at P, due to the contribute of secondary waves, is then given by:
where S describes the surface of the sphere, and s is the distance between Q and P.
Fresnel used a zone construction method to find approximate values of K for the different zones, which enabled him to make predictions that were in agreement with experimental results. The integral theorem of Kirchhoff includes the basic idea of Huygens-Fresnel principle. Kirchhoff showed that in many cases, the theorem can be approximated to a simpler form that is equivalent to the formation of Fresnel's formulation. 
For an aperture illumination consisting of a single expanding spherical wave, if the radius of the curvature of the wave is sufficiently large, Kirchhoff gave the following expression for K(?):
K has a maximum value at ? = 0 as in the Huygens-Fresnel principle; however, K is not equal to zero at ? = ?/2, but at ? = ?.
Above derivation of K(?) assumed that the diffracting aperture is illuminated by a single spherical wave with a sufficiently large radius of curvature. However, the principle holds for more general illuminations. An arbitrary illumination can be decomposed into a collection of point sources, and the linearity of the wave equation can be invoked to apply the principle to each point source individually. K(?) can be generally expressed as: 
In this case, K satisfies the conditions stated above (maximum value at ? = 0 and zero at ? = ?/2).
Feynman defines the generalized principle in the following way:
"Actually Huygens' principle is not correct in optics. It is replaced by Kirchoff's modification which requires that both the amplitude and its derivative must be known on the adjacent surface. This is a consequence of the fact that the wave equation in optics is second order in the time. The wave equation of quantum mechanics is first order in the time; therefore, Huygens' principle is correct for matter waves, action replacing time."
This clarifies the fact that in this context the generalized principle reflects the linearity of quantum mechanics and the fact that the quantum mechanics equations are first order in time. Finally only in this case the superposition principle fully apply, i.e. the wave function in a point P can be expanded as a superposition of waves on a border surface enclosing P. Wave functions can be interpreted in the usual quantum mechanical sense as probability densities where the formalism of Green's functions and propagators apply. What is note-worthy is that this generalized principle is applicable for "matter waves" and not for light waves any more. The phase factor is now clarified as given by the action and there is no more confusion why the phases of the wavelets are different than the one of the original wave and modified by the additional Fresnel parameters.
As per Greiner  the generalized principle can be expressed for in the form:
Where G is the usual Green function that propagates in time the wave function . This description resembles and generalize the initial Fresnel's formula of the classical model.
Huygens' theory served as a fundamental explanation of the wave nature of light interference and was further developed by Fresnel and Young but did not fully resolve all observations such as the low-intensity double-slit experiment first performed by G. I. Taylor in 1909. It was not until the early and mid-1900s that quantum theory discussions, particularly the early discussions at the 1927 Brussels Solvay Conference, where Louis de Broglie proposed his de Broglie hypothesis that the photon is guided by a wave function. The wave function presents a much different explanation of the observed light and dark bands in a double slit experiment. In this conception, the photon follows a path which is a random choice of one of many possible paths. These possible paths form the pattern: in dark areas, no photons are landing, and in bright areas, many photons are landing. The set of possible photon paths is determined by the surroundings: the photon's originating point (atom), the slit, and the screen. The wave function is a solution to this geometry. The wave function approach was further proven by additional double-slit experiments in Italy and Japan in the 1970s and 1980s with electrons.
Huygens' principle can be seen as a consequence of the homogeneity of space--space is uniform in all locations. Any disturbance created in a sufficiently small region of homogenous space (or in a homogenous medium) propagates from that region in all geodesic directions. The waves produced by this disturbance, in turn, create disturbances in other regions, and so on. The superposition of all the waves results in the observed pattern of wave propagation.
Homogeneity of space is fundamental to quantum field theory (QFT) where the wave function of any object propagates along all available unobstructed paths. When integrated along all possible paths, with a phase factor proportional to the action, the interference of the wave-functions correctly predicts observable phenomena. Every point on the wavefront acts as the source of secondary wavelets that spread out in the light cone with the same speed as the wave. The new wavefront is found by constructing the surface tangent to the secondary wavelets.
In 1900, Jacques Hadamard observed that Huygens' principle was broken when the number of spatial dimensions is even. From this, he developed a set of conjectures that remain an active topic of research. In particular, it has been discovered that Huygens' principle holds on a large class of homogenous spaces derived from the Coxeter group (so, for example, the Weyl groups of simple Lie algebras).